In comparison to mitochondrial genomes from other members of the Ericaceae, this assembly is similar to that of Vaccinium macrocarpon at 459 678 bp but considerably smaller than that of the saprophytic Hypopitys monotropa or Rhododendron simsii .Our assembly is the first sequenced genome in the genus Arctostaphylos, representing the initial step toward understanding genetics and adaptation in this highly diversified genus. In addition to enhancing ongoing phylogenetic and conservation research, this assembly will enable investigations of the genetic basis of adaptations including drought tolerance and fire resilience. These traits, which are ecological hallmarks of manzanitas, are of growing importance in the context of the increased drought and fire frequency and intensity that are occurring in California as a result of climate change. This assembly can also serve as a reference for studying the diversification and population genetics of Arctostaphylos, which may shed light on aspects of diversification in other complex groups of the CFP. Arctostaphylos is the third genus with a genome assembly in the heath family, Ericaceae, following release of assembled genome sequences of Rhododendron and Vaccinium . These two genera are of significant economic importance: many species, hybrids, and cultivars of Rhododendron, including rhododendron and azalea, are important landscape and ornamental plants, and the fruits of many Vaccinium species, container raspberries which include cranberry, blueberry, and huckleberry, are consumed by humans and other animals. Species in the Ericaceae are also notable for their ability to tolerate acidic and nutrient-poor soils that often characterize boreal forests and bogs, allowing them to thrive in habitats that are inaccessible to most species.

Their tolerance for these conditions is due in part to the formation of mutualistic associations between the roots of the plants and soil fungi of a type unique to the heath family known as ericoid mycorrhizae. Ericoid mycorrhizae are distinct from common mycorrhizal associations found in most angiosperms, and are far less well understood . Complete genome sequences from three genera in this family will provide a strong foundation for investigating the basis of this unique mutualism and its ability to promote survival in inhospitable soils. The size of the A. glauca assembly is 547Mb, which is similar to the two Rhododendron genomes and half that of V corymbosum, which is tetraploid . The tetraploid nature of V. corymbosum also explains the vastly greater number of duplicated genes in its assembly compared to the two diploid assemblies. The scaffold N50 of the A. glauca assembly is longer than R. williamsianum, and close to R. simsii and V. corymbosum, suggesting that the contiguity of Arctostaphylos is comparable to the other taxa . Analysis using RepeatModeler indicated that 57.71% of the A. glauca genome is composed of different categories of repetitive elements . In contrast, analysis using RepeatModeler identified only 26%, 47.5%, and 44.3% of the genome comprising repeat elements in R. williamsianum, R. simsii, and V. corymbosum respectively. The BUSCO completeness assessment of the A. glauca assembly is higher than R. williamsianum and close to the V. corymbosum , indicating that our final assembly is high quality . Overall, the A. glauca, R. simsii and V. corymbosum genomes are of comparably high contiguity and completeness. The lower contiguity and completeness of the R. williamsianum genome may be due to the lack of HiFi or other longread data in the assembly. This explanation is consistent with other studies demonstrating improved assembly with the inclusion of longer reads . Although the big berry manzanita is a common and widespread species, nearly half of the 60+ manzanita species are rare or threatened.

Many are now represented by only one or two populations, and are thus vulnerable to complete eradication by the increasingly common and intense wildfires experienced across California each year. Our manzanita genome sequence will help fulfill the overall goal of the CCGP, serving as a key resource to assess genetic diversity in these threatened endemics and move forward with coordinated conservation programs.The increasing number of diversified small-scale farms and outdoor-raised livestock in the United States , reflects growing consumer interest and demand for sustainably-produced or organic local foods, including animal products such as meat and eggs. However, there is a lack of research evaluating the unique agricultural management practices of DSSF and if these types of farming operations involve risk factors that affect the transmission of food borne pathogens in the food supply. California is the top agricultural production state in the US with annual sales over $50 billion from 69,900 total farms. California is also the sole producer in the US of 17 crops including figs, artichokes and almonds. Additionally, California leads organic sales, accounting for 40% of all organic crop production and 18.16% of United States Department of Agriculture certified organic farms. Many organic farms are small-scale and diversified and these type of DSSF sell food directly to consumers through marketing channels, such as farmers markets or Community Supported Agriculture programs. To suit the unique characteristics of California’s diverse, and within interior valley regions, year-round growing environment, we adjusted the USDA-ERS definition of “small-scale farm” to encompass operations that gross less than $500,000 annually and market directly to consumers through farmers markets, farm stands, CSAs, etc. DSSF are defined as those operations that grow a combination of livestock and specialty crops or raise multiple livestock species with the intent of selling sustainably-raised animal products directly to consumers. Some diversified farms integrate livestock and crop production by using their animals to graze crop residues or cover crops before planting with fresh market crops. Grazing enhances soil fertility, recycles farm nutrients and animals provide another source of revenue through fiber or food products.

Many consumers perceive small-scale farms or outdoor-raised livestock as more natural and safer than food grown on large-scale conventional farms or meat animals raised in confinement systems, but animals naturally harbor food borne pathogens that can cause severe human illness, like Salmonella spp., Campylobacter spp. and Shiga toxin-producing Escherichia coli . For instance, a study by Patterson et al calculated a 4.17% prevalence of STEC in sheep raised on a diversified organic farm in California. Animals are intermittent shedders of enteric pathogens and shedding may increase under certain conditions, such as during periods of stress or due to husbandry practices . Many food borne pathogens can exist in the soil for extended periods of time and can be transmitted to humans through direct contact with feces or animals, or indirect contact with a contaminated environment or through ingestion of produce, meat or water. STEC is consistently one of the major pathogens involved in food borne outbreaks in the US. Vegetables and fruit consumed raw, including spinach, tomatoes and melons, are especially considered high-risk foods. In the summer of 2011, a small family farm u-pick berry operation was the center of an E. coli O157:H7 outbreak when strawberries were contaminated by wild deer feces. Six of the fifteen cases were hospitalized and two people died. More recently, several nationwide outbreaks of E. coli O157:H7 have occurred through consumption of romaine lettuce, including one outbreak traced back to California farms and linked to cattle grazing upstream from the lettuce fields. STEC outbreaks associated with DSSF might be under reported, due to their smaller volume of sales compared to large farms. A study by Harvey et al identified six STEC outbreaks connected with organic agricultural operations between 1992-2014. All these food borne outbreaks underscore the need to conduct prevalence studies of food borne pathogens on DSSF. Swine are especially a livestock species of concern because they are reservoirs for zoonotic diseases like swine influenza and brucellosis and food borne pathogens like STEC and Campylobacter spp. Although most swine production in the US occurs inside buildings with high levels of bio-security, draining pots the US is currently experiencing a resurgence of outdoor-based swine operations, due to consumer demand for sustainably-raised animal products. Still primarily considered a niche production system in the US, outdoor-raised pig operations  are numerous and broadly distributed within California. The USDA National Animal Health Monitoring System swine report defined “small-enterprise operations” as those raising fewer than 100 pigs. Approximately 68.8 – 78.9% of the small operations included in this nationwide survey raised domestic pigs with some level of outside access. A challenge in raising domestic pigs outdoors is the increasing risk of their directly or indirectly interacting with wildlife like feral pigs, and a subsequent potential increase for pathogen sharing, especially as feral pig abundance and distribution grows throughout the US. Feral pigs are considered an invasive species as they only need water, food and shrub cover to survive, can double their population in four months and are difficult to eradicate., Moreover, if an area contains favorable habitat for feral pigs, then their population numbers can be maintained or increase over time. They also have the widest geographic distribution and one of the broadest habitat ranges of any large mammal except humans.

The wide distribution of feral swine is in part due to their ability to adapt to many ecological habitats and their opportunistic omnivore diet. California has one of the largest and widest distributions of feral pigs. Feral pigs in the US are a mix of introduced Eurasian wild boars, which are native to Asia and Europe, and domestically-raised pigs turned feral. As early as 2005, Corn et. al described the disease implications of expanding feral pig populations in the US, because they can serve as a vehicle for pathogen transmission to domestic pigs, and they could play a significant role in the transmission and maintenance of transboundary animal diseases , that may be introduced or re-introduced to North America , classical swine fever). Other studies have also reported the disease risks from the expanding distribution of feral pigs in the US, including their future role as spreaders of TBD like ASF, which was recently found in the Dominican Republic, a mere 700 miles from Florida, US. Additionally, eradicated diseases in indoor-pig herds have been documented in feral swine populations in California, for example, feral pig samples collected by the California Department of Fish and Wildlife from 1978- 2013 showed feral pigs testing positive for Brucella suis, Leptospira spp. and swine influenza virus . Contact between feral pigs and outdoor-raised pig herds increases the risk for the transmission of these diseases in domestic swine. Many studies have reported that feral pigs maintain and transmit zoonotic and food borne pathogens; however, only a small subset of these studies focused on the risk of pathogen sharing between feral pigs and outdoor raised pigs. Transmission of pathogens between feral pigs and outdoor raised pigs has been documented in the US. For example, a 2016 human brucellosis case on a NewYork State farm began with a feral pig infecting domestic pigs reared outdoors. Swine sold from this index farm led to Brucella suis positive domestic swine in nine other herds in multiple states. Additionally, one swine brucellosis case each in Texas, Iowa and Georgia in 2005 also involved domestic swine being exposed to feral pigs through inadequate bio-security or wildlife controls. Feral pigs are known to forage on farmland and some California farmers and ranchers regularly experience feral pig intrusions in their crop fields and/or contact between outdoor raised pigs and feral pigs. According to studies conducted in California and Texas, contact has been documented between feral pigs and outdoor-raised pigs. A 2012 spatial study by Wyckoff et al reported that feral pigs are attracted to agricultural habitats as food sources, which may facilitate pathogen transmission to livestock raised outdoors and humans or contaminate crops. The authors assessed habitat and movement of feral swine within 10 miles of outdoor domestic pig operations in Texas and calculated that at least 50% of these facilities were surrounded by suitable feral pig habitat. Another Wyckoff et al study assessed the disease transmission risk of feral pigs near domestic pigs facilities in Texas. This 2009 study used GPS collars to quantify contact between feral and domestic pigs and detected evidence of direct contact, as well as antibodies for the same diseases in both swine groups. They concluded that feral swine are an increasing risk for the reintroduction of eradicated diseases as well as emerging TBD, especially for operations that allowed domestic swine outdoor access, as male feral pigs are attracted to female pens. International studies have also assessed the risk of disease transmission at the wild boar-domestic pig interface.

This is consistent with the higher NCED3, NCED5 and BAM1 transcript abundance in RNO berries . Thus, there are complex responses of ABA metabolism and signaling. It would appear that there may be two different ABA pathways affecting ABA concentrations and signaling: one involved with embryo development and one involved with the water status in the skins. Auxin is also involved with ABA signaling during the late stages of embryo development in the seeds. Auxin signaling responses are complex. ABF5 is an auxin receptor that degrades Aux/IAA proteins, which are repressors of ARF transcriptional activity. Thus, a rise in auxin concentration releases Aux/IAA repression of ARF transcription factors, activating auxin signaling. In the berry skins, there was a diversity of transcriptional responses of Aux/IAA and ARF genes in the two locations, some with increased transcript abundance and others with decreased transcript abundance. As with ABA signaling, there may be multiple auxin signaling pathways operating simultaneously. One pathway appears to involve seed dormancy. ARF2 had a higher transcript abundance in BOD berries. ARF2 promotes dormancy through the ABA signaling pathway. This is consistent with the hypothesis that BOD berries reach maturity at a lower sugar level than RNO berries.Grapevines have very dynamic gene expression responses to pathogens. The top 150 DEGs for BOD berries were highly enriched with biotic stress genes. The BOD vineyard site had a higher rainfall and higher relative humidity than RNO and these conditions are likely to be more suitable for fungi to grow. We detected a much higher transcript abundance of powdery mildew-responsive genes in BOD berries and this may be connected to a higher transcript abundance of ethylene and phenylpropanoid genes as part of a defense response.

The transcript abundance profiles of some of these genes are remarkably similar. Increased ethylene signaling in grapevines has been associated with powdery mildew infection and phenylpropanoid metabolism and appears to provide plant protection against the fungus. Genes involved with phenylpropanoid metabolism, blueberry pot size especially PAL and STS genes, appear to be quite sensitive to multiple stresses in the environment . In Arabidopsis there are four PAL genes. These PAL genes appear to be involved with flavonoid biosynthesis and pathogen resistance in Arabidopsis. Ten different PAL1 and two PAL2 orthologs had higher transcript abundance in BOD berry skins; many STS genes also had a higher transcript abundance in BOD berry skins . Stilbenes are phytoalexins and provide pathogen resistance in grapes and STS genes are strongly induced by pathogens. Thus, the higher transcript abundance of powdery mildew genes may be associated with the higher transcript abundance of genes in the ethylene and phenylpropanoid pathways.The transcript abundance of a number of iron homeostasis genes were significantly different in the two locations and there was a difference in soil available iron concentrations in the two locations. However, iron uptake and transport in plants is complicated depending on multiple factors, such as pH, soil redox state, organic matter composition, solubility in the phloem, etc. Thus, it is impossible to predict iron concentrations in the berry without direct measurements. The roles of these genes in iron homeostasis and plant physiological functions are diverse. Iron supply can affect anthocyanin concentrations and the transcript abundance of genes in the phenylpropanoid pathway in Cabernet Sauvignon berry skins. One of the DEGs, SIA1, is located in the chloroplast in Arabidopsis and appears to function in plastoglobule formation and iron homeostasis signaling in concert with ATH13. Another DEG, YSL3, is involved in iron transport. It acts in the SA signaling pathway and appears to be involved in defense responses to pathogens.

It also functions in iron transport into seeds. FER1 is one of a family of ferritins found in Arabidopsis . VIT1 and NRAMP3 are vacuolar iron transporters and are also involved in iron storage in seeds. Other DEGs are also responsive to iron supply. IREG3 appears to be involved in iron transport in plastids; its transcript abundance increases with increasing iron concentrations. ABCI8 is an iron-stimulated ATPase located in the chloroplast that functions in iron homeostasis. It is unclear what specific roles these iron homeostasis genes are playing in grape berry skins, but they appear to be involved in iron storage in seeds and protection against oxidative stress responses. One possible explanation for the transcript abundance profiles in the BOD and RNO berry skins is that ferritins are known to bind iron and are thought to reduce the free iron concentrations in the chloroplast, thus, reducing ROS production that is caused by the Fenton reaction. As chloroplasts senesce during berry ripening, iron concentrations mayrise as a result of the catabolism of iron-containing proteins in the thylakoid membranes; thus, berry skins may need higher concentrations of ferritins to keep free iron concentrations low. This might explain the increase in ferritin transcript abundance with increasing sugar levels. Most soils contain 2 to 5% iron including available and unavailable iron; soils with 15 and 25 μg g− 1 of available iron are considered moderate for grapevines, but soils with higher concentrations are not considered toxic. Therefore, for both soils in this study, iron concentrations can be considered to be very high but not toxic. The higher available iron concentrations in the BOD vineyard may be associated with the wetter conditions and the lower soil pH.Other researchers using Omics approaches have identified environmental factors that influence grape berry transcript abundance and metabolites. One study investigated the differences in transcript abundance in berries of Corvina in 11 different vineyards within the same region over 3 years. They determined that approximately 18% of the berry transcript abundance was affected by the environment.

Climate had an overwhelming effect but viticultural practices were also significant. Phenylpropanoid metabolism was very sensitive to the environment and PAL transcript abundance was associated with STS transcript abundance. In another study of a white grape cultivar, Garganega, berries were analyzed by transcriptomic and metabolomic approaches. Berries were selected from vineyards at different altitudes and soil types. Again, phenylpropanoid metabolism was strongly influenced by the environment. Carotenoid and terpenoid metabolism were influenced as well. Two studies investigated the grape berry transcriptomes during the ripening phase in two different regions of China, a dry region in Western China and a wet region in Eastern China. These two locations mirror some of the differences in our conditions in our study, namely moisture, light and elevation, plant raspberry in container although the dry China western region has higher night temperatures and more rainfall than the very dry RNO location. In the Cabernet Sauvignon study, they compared the berry transcriptomes from the two regions at three different stages: pea size, veraison and maturity. The TSS at maturity was slightly below 20°Brix. Similar to our study, the response to stimulus, phenylpropanoid and diterpenoid metabolism GO categories were highly enriched in mature berries between the two locations. Differences in the transcript abundance of NCED and PR proteins were also noted. Like in our study, the authors associated the transcript abundance of these proteins to the dry and wet locations, respectively. In the second study comparing these two regions in China , the effects of the environment on the metabolome and transcriptome of Muscat Blanc à Petits Grains berries were investigated over two seasons; specifically, terpenoid metabolism was targeted. Like in our study, the transcripts in terpenoid were in higher abundance in the wetter location. The transcript abundances were correlated with terpenoid concentrations and a coexpression network was constructed. A specific set of candidate regulatory genes were identified including some terpene synthases , glycosyl transferases and 1-hydroxy-2-methyl-2-butenyl 4-diphosphate reductase . We examined the transcript abundance of some of these candidate genes in our own data but did not find significant differences between our two locations. The contrasting results between our study and Wen et al. could be for a variety of reasons such as different cultivar responses, berry versus skin samples, or different environmental conditions that affect terpenoid production. Terpenoid metabolism is influenced by the microclimate and is involved in plant defense responses to pathogens and insects. Light exposure to Sauvignon Blanc grapes was manipulated by removing adjacent leaves without any detectable differences in berry temperatures. Increased light exposure increased specific carotenoid and terpene concentrations in the berry. Some effect of temperature was associated with carotenoid and terpenoid production, but to a lesser extent than light. Higher concentrations of rotundone, a sesquiterpene, have been associated with cooler temperatures. Water deficit can also alter carotenoid and terpenoid metabolism in grapes.

Terpenes can act as signals for insect attacks and attract insect predators. Thus, terpenoid metabolism is highly sensitive to the environment and influenced by many factors. In contrast to these studies, excess light and heat can affect transcript abundance and damage berry quality. In addition to a higher rate of malate catabolism, anthocyanin concentrations and some of the transcript abundances associated with them are decreased as well.BOD berries reached maturity at a lower °Brix level than RNO berries; the cause is likely to be the warmer days and cooler nights in RNO. Higher day temperature may increase photosynthesis and sugar transport and cooler night temperatures may reduce fruit respiration. °Brix or TSS approximates the % sugar in a berry and is a reliable marker of berry maturity in any given location ; however, TSS is an unreliable marker of berry maturity when comparing grapes from very different climates. The differences in TSS between BOD and RNO are consistent with other studies on the temperature effects on berry development. Indirect studies have associated gradual warming over the last century to accelerated phenology and increased sugar concentrations in the grape berries . Increasing temperature can accelerate metabolism, including sugar biosynthesis and transport, but the increase in metabolism is not uniform. For example, the increase in anthocyanin concentration during the ripening phase is not affected as much as the increase in sugar concentration. These responses vary with the cultivar, complicating this kind of analysis even further. Direct studies of temperature effects on Cabernet Sauvignon berry composition also are consistent with our data. In one study, the composition of Cabernet Sauvignon berries was altered substantially for vines grown in phytotrons at 20 or 30 °C temperatures. Cooler temperatures promoted anthocyanin development and malate concentrations and higher temperatures promoted TSS and proline concentrations. In a second study, vines were grown at 20 or 30 °C day temperatures with night temperatures 5 °C cooler than the day. In this study, higher temperatures increased berry volume and veraison started earlier by about 3 to 4 weeks. The authors concluded that warmer temperatures hastened berry development. In a third study, Cabernet Sauvignon berry composition was affected in a similar manner by soil temperatures that differed by 13 °C. TSS concentrations are also affected by light and the vine water status. Light is generally not a factor because there is usually a large enough leaf area and sufficient light levels to saturate this source to sink relationship. Sun-exposed Cabernet Sauvignon berries in the vineyard had higher TSS than shaded berries. This sunlight effect was attributed largely to an increase in berry temperature rather than an increase in the fluence rate per se. A higher grapevine water status results in larger berry size and lower sugar concentrations and water deficit is known to increase sugar concentrations in Cabernet Sauvignon. However, temperature is thought to have the largest effect on sugar concentrations. Other transcriptomic data in the present study indicated that BOD berries were more mature at a lower sugar level than RNO berries. These included the transcript abundance profiles of genes involved in autophagy, auxin and ABA signaling, iron homeostasis and seed development. Many of these DEGs had an accelerated rate of change in BOD berries. While these transcripts are in the skins, they may be influenced by signals coming from the seed. In addition, there was a higher transcript abundance for most genes involved with the circadian clock in BOD berries. PHYB can regulate the circadian clock and PHYB activity is very sensitive to night temperatures ; PHYB reversion is accelerated to the inactive form at warmer temperatures . The inactivity of phytochrome promotes the expression of RVE1, which promotes auxin concentrations and seed dormancy. Thus, all things considered, it is likely that temperature and/or the temperature differentials between day and night significantly contributed to the differences in the rate of berry development and sugar accumulation in the two locations.

Grapevine berry ripening can be divided into three major stages. In stage 1, berry size increases sigmoidally. Stage 2 is known as a lag phase where there is no increase in berry size. Stage 3 is considered the ripening stage. Veraison is at the beginning of the ripening stage and is characterized by the initiation of color development, softening of the berry and rapid accumulation of the hexoses, glucose and fructose. Berry growth is sigmoidal in Stage 3 and the berries double in size. Many of the flavor compounds and volatile aromas are derived from the skin and synthesized at the end of this stage. Many grape flavor compounds are produced as glycosylated, cysteinylated and glutathionylated precursors and phenolics and many of the precursors of the flavor compounds are converted to various flavors by yeast during the fermentation process of wine. Nevertheless, there are distinct fruit flavors and aromas that are produced and can be tasted in the fruit, many of which are derived from terpenoids, fatty acids and amino acids. Terpenes are important compounds for distinguishing important cultivar fruit characteristics. There are 69 putatively functional, 20 partial and 63 partial pseudogenes in the terpene synthase family that have been identified in the Pinot Noir reference genome. Terpene synthases are multi-functional enzymes using multiple substrates and producing multiple products. More than half of the putatively functional terpene synthases in the Pinot Noir reference genome have been functionally annotated experimentally and distinct differences have been found in some of these enzymes amongst three grape varieties: Pinot Noir, Cabernet Sauvignon and Gewürztraminer. Other aromatic compounds also contribute significant cultivar characteristics. C13-norisoprenoids are flavor compounds derived from carotenoids by the action of the carotenoid cleavage dioxygenase enzymes.

Cabernet Sauvignon, Sauvignon Blanc and Cabernet Franc are characterized by specific volatile thiols and methoxypyrazines. Enzymes involved in the production of these aromas have been recently characterized. Phenolic compounds play a central role in the physical mouthfeel properties of red wine; recent work relates quality with tannin levels. While the grape genotype has a tremendous impact on tannin content, blackberries in containers the environment also plays a very large role in grape composition. The pathway for phenolic biosynthesis is well known, but the mechanisms of environmental influence are poorly understood. Ultimately, there is an interaction between molecular genetics and the environment. Flavor is influenced by climate, topography and viticultural practices. For example, water deficit alters gene expression of enzymes involved in aroma biosynthesis in grapes, which is genotype dependent, and may lead to increased levels of compounds, such as terpenes and hexyl acetate, that contribute to fruity volatile aromas. The grapevine berry can be subdivided into the skin, pulp and seeds. The skin includes the outer epidermis and inner hypodermis . A thick waxy cuticle covers the epidermis. The hypodermal cells contain chloroplasts, which lose their chlorophyll at veraison and become modified plastids; they are the sites of terpenoid biosynthesis and carotenoid catabolism. Anthocyanins and tannins accumulate in the vacuoles of hypodermal cells. Pulp cells are the main contributors to the sugar and organic acid content of the berries. Pulp cells also have a much higher set of transcripts involved in carbohydrate metabolism, but a lower set of transcripts involved in lipid, amino acid, vitamin, nitrogen and sulfur metabolism than in the skins. Hormones can influence berry development and ripening. Concentrations of auxin, cytokinins and gibberellins tend to increase in early fruit development of the first stage. At veraison, these hormone concentrations have declined concomitant with a peak in abscisic acid concentration just before veraison.

Auxin prolongs the Stage 2 lag phase and inhibits anthocyanin biosynthesis and color development in Stage 3. Grapevine, a non-climacteric fruit, is not very sensitive to ethylene; however, ethylene appears to be necessary for normal fruit ripening. Ethylene concentration is highest at anthesis, but declines to low levels upon fruit set; ethylene concentrations rise slightly thereafter and peak just before veraison then decline to low levels by maturity. Ethylene also plays a role in the ripening of another non-climacteric fruit, strawberry. ABA also appears to be important in grape berry ripening during veraison when ABA concentrations increase resulting in increased expression of anthocyanin biosynthetic genes and anthocyanin accumulation in the skin. ABA induces ABF2, a transcription factor that affects berry ripening by stimulating berry softening and phenylpropanoid accumulation. In addition, ABA affects sugar accumulation in ripening berries by stimulating acid invertase activity and the induction of sugar transporters. It is not clear whether ABA directly affects flavor volatiles , but there could be indirect effects due to competition for common precursors in the carotenoid pathway. Many grape berry ripening studies have focused on targeted sampling over a broad range of berry development stages, but generally with an emphasis around veraison, when berry ripening is considered to begin. In thisstudy, a narrower focus is taken on the late ripening stages where many berry flavors are known to develop in the skin. We show that that the abundance of transcripts involved in ethylene signaling is increased along with those associated with terpenoid and fatty acid metabolism, particularly in the skin.Cabernet Sauvignon clusters were harvested in 2008 from a commercial vineyard in Paso Robles, California at various times after veraison with a focus on targeting °Brix levels near maturity. Dates and metabolic details that establish the developmental state of the berries at each harvest are presented in Additional file 1. Berries advanced by harvest date with the typical developmental changes for Cabernet Sauvignon: decreases in titratable acidity and 2- isobutyl-3-methoxypyrazine concentrations and increases in sugar and color . Transcriptomic analysis focused on four harvest dates having average cluster °Brix levels of 22.6, 23.2, 25.0 and 36.7. Wines made in an earlier study from grapes harvested at comparable levels of sugars or total soluble solids to those in the present study showed clear sensory differences. Six biological replicates, comprising two clusters each, were separated into skins and pulp in preparation for RNA extraction and transcriptomic analysis using the NimbleGen Grape Whole-Genome Microarray. Thus, a 4 × 2 factorial experimental design was established. After standard microarray processing and data normalization, two-way ANOVA indicated that the transcript abundance of 16,280 transcripts statistically significantly changed across the °Brix levels below the adjusted p-value of 0.05 , the transcript abundance of 10,581 transcripts changed significantly across Tissue types, and the abundance of 2053 transcripts changed significantly with respect to the °Brix x Tissue interaction term p-value column: adjBrix, adjTissue or adjTissue*Brix. A note of caution must be added here. There are high similarities amongst members in certain Vitis gene families , making it very likely that cross-hybridization can occur with probes on the microarray with high similarity to other genes. We estimate approximately 13,000 genes have the potential for cross-hybridization, with at least one probe of a set of four unique probes for that gene on the microarray potentially cross-hybridizing with probes for another gene on the microarray. Genes with the potential for cross hybridization have been identified and are highlighted in light red in Additional file 2. The rationale to include them is that although individual genes can not be uniquely separated, the probe sets can identify a gene and its highly similar gene family members, thus, providing some useful information about the biological responses of the plant. An additional approach was taken, removing cross-hybridizing probes before quantitative data analysis . Many of the significant genes were unaffected by this processing, blackberry container but 3600 genes were completely removed from the analysis. Thus, it was felt that valuable information was lost using such a stringent approach.

The less stringent approach allowing for analysis of genes with potential cross hybridization was used here in the rest of the analyses. To assess the main processes affected by these treatments, the gene ontologies of significantly affected transcripts were analyzed for statistical significance using BinGO. Based on transcripts that had significant changes in abundance with °Brix level, 230 biological processes were significantly over represented in this group . The three top over represented processes were response to abiotic stress, bio-synthetic process, and response to chemical stimulus, a rather generic set of categories. Tissue differences were more revealing at the stage when flavors peak; 4865 transcripts that were significantly higher in skins compared to pulp at 23.2 °Brix were tested for over represented GO functional categories . Some of the top GO categories included photosynthesis, isoprenoid biosynthesis, and pigment biosynthesis . Some of the transcripts with the largest differences between skin and pulp at 23.2 °Brix are β-ketoacyl-CoA synthase , taxane 10-β-hydroxylase , wax synthase, a lipase, an ABC transporter, and phenylalanine ammonia-lyase . The abundance of 5716 transcripts was significantly higher in pulp than skin at 23.2 °Brix . Some of the top GO categories over represented were a variety of transport processes and small GTPase mediated signal transduction . Some of the transcripts with the largest differences in abundance with pulp greater than skin at 23.2 °Brix were polygalacturonase , flavonol synthase, stachyose synthase, an amino acid transporter, a potassium channel , and HRE2 . The transcript abundance of 2053 genes had significantly differential expression across °Brix levels and tissues . The top GOcategories over represented in this set involved photosynthesis and phenylpropanoid metabolism, both associated with the berry skin . Other flavorcentric categories of the 57 categories over represented include aromatic compound biosynthesis, fatty acid metabolism and alcohol catabolism. This transcript set was further analyzed by dividing into 10 clusters using k-means clustering . The over represented GO categories were determined for each cluster . Eight of the 10 clusters had distinct over represented GO categories; two clusters did not have any over represented GO categories, meaning that the genes in these two clusters were assigned to GO categories of expected proportions when compared to the entire NimbleGen array. Clusters 1, 8, 9 and 10 had a large number of over represented categories. Many GO categories within a cluster are subsets of others in that cluster and were grouped together. For example, cluster 4 had four over represented GO categories, oxygen transport, gas transport, heat acclimation and response to heat. The four categories could be grouped into two, as two are subsets of the others; this is how they were listed in Table 1.As we were interested in compounds associated with berry flavors as they develop or change in the late stages of berry ripening, we took a more targeted approach for analysis with this in mind. Berries at 24° Brix are known to be near-optimal for flavor, thus we took a simple approach to look for genes that were peaking around this stage. We found some significant and large increases in transcript abundance between the 22.6 and 23.2 °Brix levels. A group of VviERF6 transcription factor paralogs represented 6dance from 22.6 to 23.2 °Brix in the skin, but not in the pulp . These VviERF6 TFs were also found in Cluster 8 . This is very interesting since many flavor compounds are derived from the skin and ERF TFs are known to be responsive to ethylene, a known fruit-ripening hormone. These VviERF TFs were named ERF105 in the annotation by Grimplet et al., however they are more orthologous with AtERF6 as determined by a more comprehensive phylogenetic method using many plant species at Gramene . Annotation details of the V1 gene models of the VviAP2/ERF super family can be found in Additional file 8 including updated Vvi symbols according to its closest Arabidopsis ortholog as instructed by the Grapevine Gene Nomenclature System developed by the International Grape Genome Program Supernomenclature committee. This renaming of the AP2/ERF super family should facilitate comparative analyses and functions with other species, particularly Arabidopsis. To properly annotate the AP2/ERF super family of Vitis vinifera according to the IGGP Supernomenclature committee instructions, a phylogenetic tree was generated for the AP2/ERF super family of Arabidopsis thaliana and Vitis vinifera using the TAIR 10 and V1 gene models, respectively . The labeled family classifications were derived from the Arabidopsis naming scheme by Nakano et al.. There are 130 members in the VitisAP2/ERF super-family in the Pinot Noir reference genome. However, the six paralogs of ERF6 discussed above belong to a Vitis vinifera clade in subfamily IX and are distinctly different or separate from any Arabidopsis subfamily IX ERF TFs .

To date, quantum anomalous Hall effects have been observed at band fillings ν = 1 and ν = 3 in various graphene heterostructures, where ν = An corresponds to the number of electrons per unit cell area A with n the carrier density. Although orbital magnetism is generally expected theoretically in twisted bilayer graphene, no direct experimental probes of magnetism have been reported because of the relative scarcity of magnetic samples, their small size, and the low expected magnetization density they are predicted to have.We perform spatially resolved magnetometry to image the submicron magnetic structure of the same sample presented in the previous chapter and in Ref. , which consists of a twisted graphene bilayer aligned to one of the hexagonal boron nitride encapsulating layers. Figure 5.1B shows a schematic representation of our experimental setup. We use a superconducting quantum interference device fabricated on the tip of a quartz tube from cryogenically deposited indium with a magnetic field sensitivity of approximately 15 nT/Hz1/2 at select outof-plane magnetic fields of less than 50 mT. The SQUID is mounted to a quartz tuning fork and rastered in a 2D plane parallel to, and at a fixed height above, the tBLG heterostructure. Here n is the nominal density inferred from a parallel plate capacitor model, with the capacitance determined from the lowfield Hall density . Images are acquiredin the same background magnetic field B = 22 mT but on opposite branches of the hysteresis loop shown in Fig. 5.1A. As demonstrated in Fig. 8.16, and explained in more detail in the appendix, raspberry container the measured BT F contains contributions from magnetic signals as well as other effects arising from electric fields or thermal gradients. To isolate the magnetic structure that gives rise to the observed hysteretic transport, we subtract data from Figs. 5.1, C and D, from each other.

The result is shown in Fig. 5.1E, which depicts the gradient magnetometry signal associated with the fully polarized orbital ferromagnet. To reconstruct the static out-of-plane magnetic field, Bz, we integrate BT F along ˆa from the lower and left boundaries of the image . We infer the total magnetization density m from the Bz data using standard Fourier domain techniques as shown in Fig. 5.1G. The algorithm used for this process is covered in depth in the associatd publication. Figures 5.1, H and I, show a comparison of BT F and m plotted along the contours indicated in Figs. 5.1, E and G. The shaded regions in Fig. 5.1I denote confidence intervals, whose absolute error bounds form the dominant systematic uncertainty in |ˆa|. Our measurements are taken close to ν = 3, equivalent to a single hole per unit cell relative to the nonmagnetic state at ν = 4 that corresponds to full filling of the lowest energy bands. We find that the magnetization density is considerably larger than 1 µB per unit cell area A ≈ 130/nm2 , where we have taken g = 2 as appropriate for graphene, in which spin orbit coupling is negligible. Without any assumptions about the nature of the broken symmetries, this state has a maximum spin magnetization of 1 µB per moir´e unit cell. Our data reject this hypothesis, finding instead a maximum magnetization density of m in the range 2 − 4 µB per moir´e unit cell corresponding to an orbital magnetization of 1.8-3.6× 10−4µB/carbon atom. We conclude that the magnetic moment associated with the Chern magnet phase in tBLG is dominated by its orbital component.In an intrinsic orbital magnet in which all moments arise from conduction electrons, the magnetization depends strongly on the density. Additional density dependence arises from the fact that contributions to the orbital magnetization from both wave-packet angular momentum and Berry curvature need not be uniformly distributed within the Brillouin zone.

Transport observations of aquantum anomalous Hall effect measure only the total Berry curvature of a completely filled band. At partial band filling, however, extrinsic contributions from scattering complicate the relationship between transport and band properties. In contrast, measuring m provides direct information about the density-dependent occupation of the Bloch states in momentum space. Although crystalline defects on the atomic scale are unlikely in tBLG thanks to the high quality of the constituent graphene and hBN layers, the thermodynamic instability of magic angle twisted bilayer graphene makes it highly susceptible to inhomogeneity at scales larger than the moir´e period, as shown in prior spatially resolved studies. For example, the twist angle between the layers as well as their registry to the underlying hBN substrate may all vary spatially, providing potential pinning sites. Moir´e disorder may thus be analogous to crystalline disorder in conventional ferromagnets, which gives rise to Barkhausen noise as it was originally described. A subtler issue raised by our data is the density dependence of magnetic pinning; as shown in Fig. 5.3, Bc does not simply track 1/m across the entire density range, in particular failing to collapse with the rise in m in the Chern magnet gap. This suggests nontrivial dependence of either the pinning potential or the magnetocrystalline anisotropy energy on the realized many body state. Understanding the pinning dynamics is critical for stabilizing magnetism in tBLG and the growing class of related orbital magnets, which includes both moir´e systems as well as more traditional crystalline systems such as rhombohedral graphite. In order to understand the microscopic mechanism behind magnetic grain boundaries in the Chern magnet phase in tBLG/hBN, we used nanoSQUID magnetometry to map the local moir´e superlattice unit cell area, and thus the local twist angle, in this device, using techniques discussed in the literature.

This technique involves applying a large magnetic field to the tBLG/hBN device and then using the chiral edge state magnetization of the Landau levels produced by the gap between the moir´e band and the dispersive bands to extract the electron density at which full filling of the moir´e superlattice band occurs . The strength of this Landau level’s magnetization can be mapped in real space , and the density at which maximum magnetization occurs can be processed into a local twist angle as a function of position . It was noted in that the moir´e superlattice twist angle distribution in tBLG is characterized by slow long length scale variations interspersed with thin wrinkles, across which the local twist angle changes rapidly. These are also present in the sample imaged here . The magnetic grain boundaries we extracted by observing the domain dynamics of the Chern magnet appear to correspond to a subset of these moir´e superlattice wrinkles. It may thus be the case that these wrinkles serve a function in moir´e superlattice magnetism analogous to that of crystalline grain boundaries in more traditional transition metal magnets, pinning magnetic domain walls to structural disorder and producing Barkhausen noise in measurements of macroscopic properties. In tBLG, a set of moir´e subbands is created through rotational misalignment of a pair of identical graphene monolayers. In twisted monolayer-bilayer graphene a set of moir´e subbands is created through rotational misalignment of a graphene monolayer and a graphene bilayer. These systems both support Chern magnets. Both systems are also members of a class of moir´e superlattices known as homobilayers; in these systems, large plastic pots for plants the 2D crystals forming the moir´e superlattice share the same lattice constant, and the moir´e superlattice appears as a result of rotational misalignment, as illustrated in Fig. 5.17A. Homobilayers have many desirable properties; the most important one is that the twist angle can easily be used as a variational parameter for minimizing the bandwidth of the moir´e subbands, producing the so-called ‘magic angle’ tBLG and tMBG systems. Homobilayers do, however, have some undesirable properties. Although local variations in electron density are negligible in these devices, the local filling factor of the moir´e superlattice varies with the moir´e unit cell area, and thus with the relative twist angle. The tBLG moir´e superlattice is shown for two different twist angles in 5.17B-C across the magic angle regime; it is clear that the unit cell area couples strongly to twist angle in this regime, illustrating the sensitivity of these devices to twist angle disorder. The relative twist angle of the two crystals in moir´e superlattice devices is never uniform. Imaging studies have clearly shown that local twist angle variations provide the dominant source of disorder in tBLG . It is hard to exaggerate the significance of this problem to the study of moir´e superlattices. Orbital magnetism at B = 0 has only been realized in a handful of tBLG devices, and quantization of the anomalous Hall resistance has only been demonstrated in a single tBLG device, in spite of years of sustained effort by several research groups. A mixture of careful device design limiting the active area of devices and the use of local probes has allowed researchers to make many important discoveries while sidestepping the twist angle disorder issue- indeed, some exotic phases are known in tBLG only from a single device, or even from individual scanning probe experiments- but if the field is ever to realize sophisticated devices incorporating these exotic electronic ground states the problem needs to be addressed.

There is another way to make a moir´e superlattice. Two different 2D crystals with different lattice constants will form a moir´e superlattice without a relative twist angle; these systems are known as heterobilayers . These systems do not have ‘magic angles’ in the same sense that tBLG and tMBG do, and as a result there is no meaningful sense in which they are flat band systems, but interactions are so strong that they form interaction-driven phases at commensurate filling of the moir´e superlattice anyway. Indeed, many of the interaction-driven insulators these systems support survive to temperatures well above 100 K. The most important way in which heterobilayers differ from homobilayers, however, is in their insensitivity to twist angle disorder. In the small angle regime, the moir´e unit cell area of a heterobilayer is almost completely independent of twist angle, as illustrated in 5.17E-F. A new intrinsic Chern magnet was discovered in one of these systems, a heterobilayer moir´e superlattice formed through alignment of MoTe2 and WSe2 monolayers. The researchers who discovered this phase measured a well-quantized QAH effect in electronic transport in several devices, demonstrating much better repeatability than was observed in tBLG. The unit cell area as a function of twist angle is plotted for three moir´e superlattices that support Chern insulators in 5.17G, with the magic angle regime highlighted for the homobilayers, demonstrating greatly diminished sensitivity of unit cell area to local twist angle in the heterobilayer AB-MoTe2/WSe2. MoTe2/WSe2 does have its own sources of disorder, but it is now clear that the insensitivity of this system to twist angle disorder has solved the replication issue for Chern magnets in moir´e superlattices. Dozens of MoTe2/WSe2 devices showing well-quantized QAH effects have now been fabricated, and these devices are all considerably larger and more uniform than the singular tBLG device that was shown to support a QAH effect, and was discussed in the previous chapters. The existence of reliable, high-yield fabrication processes for repeatably realizing uniform intrinsic Chern magnets is an important development, and this has opened the door to a wide variety of devices and measurements that would not have been feasible in tBLG/hBN.The basic physics of this electronic phase differs markedly from the systems we have so far discussed, and we will start our discussion of MoTe2/WSe2 by comparing and contrasting it with graphene moir´e superlattices. In tBLG/hBN and its cousins, valley and spin degeneracy and the absence of significant spin-orbit coupling combine to make the moir´e subbands fourfold degenerate. When inversion symmetry is broken the resulting valley subbands can have finite Chern numbers, so that when the system forms a valley ferromagnet a Chern magnet naturally appears. Spin order may be present but is not necessary to realize the Chern magnet; it need not have any meaningful relationship with the valley order, since spin-orbit coupling is absent. MoTe2/WSe2 has strong spin-orbit coupling, and as a result, the spin order is locked to the valley degree of freedom. This manifests most obviously as a reduction of the degeneracy of the moir´e subbands; these are twofold degenerate in MoTe2/WSe2 and all other TMD-based moir´e superlattices. The closest imaginable analog of the tBLG/hBN Chern magnet in this system is one in which interactions favor the formation of a valley-polarized ferromagnet, at which point the finite Chern number of the valley subbands would produce a Chern magnet.

What they do provide is a convenient way for us to produce two dimensional monocrystalline devices with exceptionally low disorder for which electron density and band structure can be conveniently accessed as independent variables. That is valuable for furthering our understanding of condensed matter phenomena, independent of whether the fabrication procedures for making these material systems can ever be scaled up enough to be viable for use in technologies. The properties of crystals differ from the properties of atoms floating in free space because the atomic orbitals of the atoms in a crystal are close enough to those of adjacent atoms for electrons to hop between atoms. The resulting hybridization of atomic orbitals produces quantum states delocalized over the entire crystal with the capacity to carry momentum. This situation is shown in schematic form in Fig. 1.9A. For quantum states delocalized over the entire crystal, position ceases to be a useful basis. Instead, under these conditions we label electronic wave functions by their momenta, kx and ky. The atomic orbitals that prior to hybridization had discrete energy spectra now have energy spectra given by discrete functions of momentum, f. We call these functions electronic bands. Electrons loaded into the electronic bands of a two dimensional crystal will occupy the quantum states with the lowest available energies, so we can specify a maximum energy at which we expect to find electrons for any given electron density. We call that energy EF , the Fermi level. We can raise the Fermi level by adding additional electrons to the crystal, as shown in Fig. 1.9B. We have already discussed how two dimensional crystals naturally allow for manipulation of the electron density, square plant pots and thus the Fermi level. We have also already discussed how the application of an out-of-plane electric field to a two dimensional crystal will change the structure of the atomic orbitals supported by that crystal.

It naturally follows that atomic orbitals so modified will produce different electronic bands, as shown in Fig. 1.9C. It is relatively straightforward to compute how electronic bands will respond to the application of a displacement field. We will be using the momentum and energy basis for the rest of this document; this basis is known as momentum space. The simplest experiment we can perform to probe the electronic properties of a two dimensional crystal in this geometry is an electronic transport experiment, in which a voltage is applied to a region of the crystal with another region grounded, so that electrical current flows through the crystal. We can check if the crystal supports any electrical transport at all, and if it does we can measure the electrical resistance of the crystal this way, in close analogy to how this is done for three dimensional crystals. Crystals will only accept and thus conduct electrons if there are available quantum states at the Fermi level; we call these crystals metals , and they can be identified in band structure diagrams by the intersection of the Fermi level with an electronic band . Crystals without empty quantum states at the Fermi level will not accept and conduct electrons , and they can be identified in band structure diagrams with crystals for which the Fermi level does not intersect with an electronic band . There exists a variety of other experiments we can perform on two dimensional crystals in order to understand their properties. Two dimensional crystals can support electronic transport in the in-plane direction if they are metals, as shown in Fig. 1.10. Capacitors can also support electronic transport in the out-of-plane direction, as long as that electronic transport occurs at finite frequency. The same structure that we use to modify the electron density and ambient out-of-plane electric field of a two dimensional crystal can also be used as a capacitive AC conductor, as illustrated in Fig. 1.11A.

The conductance will depend only on the frequency at which an AC voltage is applied and the geometry of the parallel plate capacitor. However, if a two dimensional crystal is added in series, the capacitance of the top gate to the bottom gate may be substantially modified. If the two dimensional crystal is an insulator, electric fields will penetrate it and the capacitance between the two gates will not change. However, if the two dimensional crystal is a metal it will accept electrons and cancel the applied electric field, dramatically reducing the capacitance between the top and bottom gates and neutralizing the AC current through the capacitor. This technique can be used to measure the electronic properties of specifically the bulk of a two dimensional crystal; it is the property that was both calculated and measured in Fig. 1.2C and D. These two techniques are the bread and butter of the experimental study of two dimensional crystals, because they require only the ability to create stacks of two dimensional crystals and access to tools common to the study of all other microelectronic systems. We will discuss a considerable amount of electronic transport and capacitance data as well. However, the primary focus of this thesis will be on systems for which the nanoSQUID microscope can provide important information that is inaccessible to these techniques, and so we will discuss a few such systems next. Consider the following procedure: we obtain a pair of identical two dimensional atomic crystals. We slightly rotate one relative to the other, and then place the rotated crystal on top of the other . The resulting pattern brings the top layer atoms in alignment with the bottom layer atoms periodically, but with a lattice constant that is different from and in practice often much larger than the lattice constant of the original two atomic lattices. We call the resulting lattice a ‘moir´e superlattice.’ The idea to do this with two dimensional materials is relatively new, but the notion of a moir´e pattern is much older, and it applies to many situations outside of condensed matter physics. Pairs of incommensurate lattices will always produce moir´e patterns, and there are many situations in daily life in which we are exposed to pairs of incommensurate lattices, like when we look out a window through two slightly misaligned screens, or try to take pictures of televisions or computer screens with our camera phones. Of course these ‘crystals’ differ pretty significantly from the vast majority of crystals with which we have practical experience, so we’ll have to tread carefully while working to understand their properties. To start with, if we attempt to proceed as we normally would- by assigning atomicorbitals to all of the atoms in the unit cell, computing overlap integrals, and then diagonalizing the resulting matrix to extract the hybridized eigenstates of the system- we would immediately run into problems, because the unit cell has far too many atoms for this calculation to be feasible.

There exist clever approximations that allow us to sidestep this issue, and these have been developed into very powerful tools over the past few years, but they are mostly beyond the scope of this document. I’d like to instead focus on conclusions we can draw about these systems using much simpler arguments. The physical arguments justifying the existence of electronic bands apply wherever and whenever an electron is exposed to an electric potential that is periodic, and thus has a set of discrete translation symmetries. For this reason, even though the moir´e superlattice is not an atomic crystal, we can always expect it to support electronic band structure for the same reason that we canal ways expect atomic crystals to support band structure. Two crystals with identical crystal symmetries will always produce moir´e superlattices with the same crystal symmetry, so we don’t need to worry about putting two triangular lattices together and ending up with something else.Another property we can immediately notice is that the electron density required to fill a moir´e superlattice band is not very large. This can be made clear by simply comparing the original atomic lattice to a moir´e superlattice in real space . Full depletion of a band in an atomic crystal requires removing an electron for every unit cell , plastic pots for planting and full filling of the band occurs when we have added an electron for every unit cell. We have already discussed how this is not possible for the vast majority of materials using only electrostatic gating, because the resulting charge densities are immense. Full depletion of the moir´e band, on the other hand, requires removing one electron per moir´e unit cell, and the moir´e unit cell contains many atoms . So the difference in charge density between full filling and full depletion of an electronic band in a moir´e superlattice is actually not so great , and indeed this is easily achievable with available technology. Before we go on, I want to make a few of the limitations of this argument clear. There are two things this argument does not necessarily imply: the moir´e bands we produce might not be near the Fermi level of the system at charge neutrality, and the bandwidth of the moir´e superlattice need not be small. In the first case, we won’t be apply to modify the electron density enough to reach the moir´e band, and in the latter, we won’t be able to fill the moir´e band’s highest energy levels using our electrostatic gate. We know of examples of real systems with moir´e superlattice bands that fail each of those criteria. But if these moir´e superlattice bands are near charge neutrality, and if their bandwidths are small, then we should be able to easily fill and deplete them with an electrostic gate.Finally, moir´e superlattice bands inherit any electronic degeneracies- like, for example, electron spin- that came with the original lattice. We haven’t discussed electronic degeneracies yet, and we will shortly. So if a moir´e superlattice satisfies all of these criteria, then it will provide a set of electronic bands that can be completely filled or depleted with an electronic gate. I’m sure this seems to the reader like a pretty niche system, and that’s more or less because it is. There aren’t too many material systems that need their atomic bonds aligned with a mechanical goniometer, and it’s hard to imagine ever integrating such a procedure into an industrial fabrication line. However, it’s tough to adequately express how hard it would be to replicate the properties of a moir´e superlattice band in an atomic crystal. I made an attempt to do so in the introduction to this thesis; suffice to say the control we have over the properties of these systems is more or less unprecedented within experimental condensed matter physics, and this means that we can perform experiments on electronic phases in these systems that would be difficult or impossible in atomic crystals. A variety of scanning probe microscopy techniques have been developed for examining condensed matter systems. It’s easy to justify why magnetic imaging might be interesting in gate-tuned two dimensional crystals, but magnetic properties of materials form only a small subset of the properties in which we are interested. Scanning tunneling microscopy is capable of probing the atomic-scale topography of a crystal as well as its local density of states, and a variety of scanning probe electrometry techniques exist as well, mostly based on single electron transistors. It’s worth pointing out that if you’re interested specifically in performing a scanning probe microscopy experiment on a dual-gated device, then these techniques both struggle, because the top gate both blocks tunnel current and screens out the electric fields to which a single electron transistor would be sensitive. Magnetic fields have an important advantage over electric fields: most materials have very low magnetic susceptibility, and thus magnetic fields pass unmodified through the vast majority of materials . This means that magnetic imaging is more than just one of many interesting things one can do with a dual-gated device; in these systems, magnetic imaging is a member of a very short list of usable scanning probe microscopy techniques. The simplest way in which we can use our nanoSQUID magnetometry microscope is as a DC magnetometer, probing the static magnetic field at a particular position in space .

Loss of integrity of the cuticle and epithelial cells of the integumentary and tracheal systems in T. molitor larvae due to densovirus infection may not only cause dehydration via excessive water loss, but also secondary infection by bacteria and fungi as a result of the loss of their primary defense barrier against infection. One limitation encountered in this study was the interpretation of a dark discoloration presented by densovirus-infected T. molitor larvae. Further investigation into the mechanism of this change is suggested as there remains uncertainty to whether this is in fact hemocyte granule prophenoloxidase-mediated melanization or some other possible explanation, such as the collection of debris over an unshed and desiccated cuticle. InI is the most important microscopical feature of the densovirus infection in T. molitor. InI observed on bright-field microscopy represents the DRAC observed on electron microscopy. The infection with TmDNV induces homogeneous basophilic- and eosinophilic-texturized InIs. Feulgen staining confirmed the DNA-dominant nature of the basophilic InIs, which correlates with DRACs accumulating maturing and mature virions on TEM as is shown in Figs. 6d and 7a. On the other hand, a great number of eosinophilic-texturized InIs presented with a fainter staining after the Feulgen reaction, which correlates with DRACs shown in TEM photomicrographs in Fig. 7c–f. Eosinophilic-texturized inclusions represented DRACs that ultrastructurally displayed an evolving mixture of heterogeneous viral matrix, maturing and mature virions, and ribosomes. Viral matrix observed in DARC was composed of electron-dense fibrillar small-aggregating-to-long-anastomosing streams, which evolved into a large compact mass that resembles the shape of a “rotary dial” . Similar structures were described in early studies of G. mellonella infected with densovirus. Later, H-1 parvovirus studies indicated that H-1 parvovirus–associated replication bodies, square plastic planter which are not uncommonly observed in DRAC in T. molitor , are composed of nucleolar fibrils and serve as the location of viral DNA replication. Differentiating maturing and mature virions from ribosomes in the infected nucleus on conventional plastic-embedded preparations posed a challenge.

However, the large number of virus particles recovered on direct TEM supports the observation that large quantities of round vesicles/ particles of about 18 to 24 nm in diameter found within the DRAC were in fact maturing and mature densovirus virions rather than ribosomes. The DRAC displayed a very pleomorphic morphology, which varied from accumulation of maturing and mature virions to dense virus matrix that mimicked nucleolar structures and paracrystalline arrays. This heterogeneous appearance of the DRAC is an important feature that characterizes and facilitates TmDNV diagnosis at the ultrastructural level. Initial TEM studies on the G. mellonella densovirus demonstrated that the nucleoplasm is completely replaced by virions during the first 20 hours of infection. As in other autonomous parvoviruses that depend on the S phase of the infected cells, TmDNV replication and assembly is a unique and dynamic process that demands a deeper understanding. Nonetheless, crystallization of the DRAC forming “paracrystalline arrays” was interpreted as the final stage of nuclear infection with densovirus in T. molitor. Numerous infected cells presented cytoplasmic membrane–bound vesicles containing virus particles or paracrystalline arrays. Intracytoplasmic paracrystalline arrays have been considered a distinct feature of Densovirinae infections. The origin of the structures is not well understood, but it should be considered as a process of infected nuclear breakdown, heterophagia, or autophagia. Ultrastructural changes of the nuclei and cytoplasm in TmDNV–infected epithelial cells suggest apoptosis as the mechanism of cell death, although studies on viruses of the Parvoviridae family suggest that the mechanisms of Parvoviridae-induced cell death are unique and complex, warranting further investigation. Fruit, nut, and berry crops are commonly grouped into one of three categories: temperate, subtropical, and tropical. Temperate zone crops include almond, apple, apricot, peach, grape, blueberry, and strawberry . Avocado, citrus, and guava are considered to be subtropical, while banana, cashew, and pineapple are tropical . Generally, temperate and subtropical crops can be grown in San Mateo and San Francisco Counties , but tropical crops are rarely successful.

This publication focuses on temperate and subtropical crops. Temperate zone crops generally require a period of cold temperature during the winter months for successful flower and fruit development. This cold temperature period is measured in “chill hours” . Some crops require many chill hours, while others require few. This is called the crop’s “chill requirement.” When selecting temperate zone crops, it is important to choose only those crops that have a chill requirement that will be met at your location. Subtropical crops, such as citrus, loquat, and guava, require little or no chilling. Native to warm-climate regions, these crops can be injured by cold temperatures during winter and spring months, and they require heat during the growing season for fruit maturation and flavor.Selecting climate zones and meeting chill requirements are not the only factors necessary for good fruit production. Pollination, sunlight, heat accumulation, and wind are all important considerations. For fruit development to occur, flowers have to be pollinated; that is, pollen must move from the male organs to the female organs . Pollen transfer can be facilitated by bees, beetles, flies, butterflies, moths, birds, bats, wind, and water. The pollen can come from flowers on the same tree , or from flowers of other trees of the same species . In some cases, a crop may require another tree of a specific variety for pollination. For information on pollination requirements, see the “Crop and Variety Selection Table” at the end of this publication. For further details, see “Notes on Crops and Varieties” below. Note that even in self fertile crops, cross-pollination can increase fruit set. Also, poor pollination can occur as a result of insufficient pollinator activity, such as during cool and wet weather in the spring when bees may be less active . For most crops, plentiful sunlight and warm temperatures during the growing season are needed. Factors such as fog, wind, elevation, aspect , and distance from the ocean affect sunlight level and temperature.In recent years, a variety of novel two-dimensional van der Waals magnets have been discovered, founding the active field of 2D magnetism.

Among these prospective compounds, binary chromium tellurides Cr1–δTe are attractive owing to their rich magnetic properties, as well as inherent chemical and structural compatibility when forming heterostructures with other topological systems, such as tetradymite-type topological insulators or chalcogenide-based Dirac/ Weyl semimetals. Furthermore, the broken time-reversal symmetry and spin-orbit coupling offer unique opportunities for the interplay between spin configurations and reciprocal-space topology. In this regard, ferromagnetic Cr2Te3 with strong perpendicular magnetic anisotropy is an intriguing platform to host non-trivial topological physics, particularly for the high-quality thin films grown by molecular beam epitaxy. The intrinsic AHE is topological in nature and a hallmark of itinerant ferromagnets, which has also been observed in more exotic systems even without a net magnetization, such as spin liquids, antiferromagnets, and Weyl semimetals. When SOC coexists with long-range magnetic order, square plastic plant pot the Berry curvature can be significantly influenced near avoided band crossings, rendering the system an incredibly rich playground combining topology and magnetism. Here, we report the unique magnetotransport signatures of highquality quasi-2D Cr2Te3 MBE-grown thin films governed by non-trivial band topologies. Via synergetic structural, magnetic, and transport measurements, together with first-principles simulations, we have uncovered novel Berry-curvature-induced magnetism featuring an extraordinary sign reversal of the AHE as we modulate the temperature and the strain for the thin films containing 3–24 unit cells on Al2O3 or SrTiO3 substrates. Moreover, a hump-shaped Hall feature emerges, most likely due to the presence of multiple magnetic layers/domains under different levels of interfacial strain. This work identifies pristine ferromagnetic Cr2Te3 thin films as a fascinating platform for further engineering topological effects, given their nontrivial Berry curvature physics.The crystalline structure of Cr2Te3 thin films is described first, followed by the development of strain at the substrate/film interface by the epitaxy. Bulk Cr2Te3 crystalizes in a three-dimensional lattice with space group P31c ðD2 3d,No:163Þ, as shown in Fig. 1a–c, where each unit cell contains four vertically stacked hexagonal layers of Cr. There are three symmetrically unique sites for Cr, labeled Cr1, Cr2, and Cr3, respectively: The Cr1 atoms are sparsely arranged in a weakly antiferromagnetic sublattice, while the Cr2/Cr3 atoms form ferromagnetic layers similar to those in CrTe2. Since the Cr1 sites are often only partially filled , Cr2Te3 behaves essentially as a quasi-2D magnet. This quasi-2D nature of Cr2Te3 allows for high-quality, layer-by-layer epitaxial growth of c-oriented films on a variety of substrates. The hexagonal c axis is the easy magnetic axis, leading to PMA for the films. The sixfold in-plane symmetry is seen in the honeycombs visualized by atomic resolution scanning tunneling microscopy and scanning transmission electron microscopy high-angle annular dark-field imaging, as well as in the reflection high-energy electron diffraction and X-ray diffraction patterns.

The sharp substrate/film interface is confirmed by the cross-sectional HAADF and the corresponding integrated differential phase contrast images. The intrinsic random distribution of Cr atoms on the Cr1 sites is resolved in the enlarged view of the atoms in Fig. 1k–m, shown overlaid with red circles, while the overall chemical composition of the thin film is uniform within the resolution of energy dispersive X-ray spectroscopy . Figure 1f illustrates the basic sample architecture, where the strain in the Cr2Te3 thin films is governed by the interface with the substrate.Upon reducing the thickness t, films grown on Al2O3 can develop an IP compressive strain up to −0.15%, as determined by XRD and summarized in Fig. 1d. A higher strain level can be sustained using SrTiO3 substrates. Such control of strain is well suited for exploring interface-sensitive properties in Cr2Te3 thin films.The magnetic properties of Cr2Te3 thin films with selected thicknesses were assessed using vibrating sample magnetometry . Figure 2a shows the temperature dependence of the magnetization M for a t = 24 u.c. film on Al2O3 substrate with an out-of-plane applied magnetic field μ0H = 0.1 T. Under the field-cool condition, M rises below the Curie temperature TC ~ 180 K, reaching M ~2.50 μB per Cr at 2 K in the 0.1 T field. The zero-FC scan, on the other hand, deviates from the FC curve below the blocking temperature Tb, signaling the freezing out of domains in a random direction in the absence of an aligning field. As illustrated in Fig. 2b, Cr2Te3 favors PMA with coercive field μ0Hc = 0.76 T and saturation magnetization Ms ~2.83 μB per Cr at 2 K for t = 24 u.c., whereas the IP measurements have weaker ferromagnetic hysteresis loops. The low-T zero-field kink feature in the OOP M becomes more prominent with reduced thickness . The multistep hysteresis attests to the presence of varied layer-dependent magnetic anisotropies, despite the overall chemical and phase homogeneity of the films34. This is consistent with the interfacial strain-driven magnetic profiles revealed by the depth-sensitive polarized neutron reflectometry as described below. The PNR experiments, responsive to the IP magnetization, were carried out at chosen T and H on samples with t = 24 and 6 u.c. to uncover the impact of interfacial strain and the details of the stepwise hysteresis loops due to the interplay between anisotropy and the Zeeman energies in an applied external magnetic field. The PNR spin asymmetry ratio SA = /, measured as a function of the wave vector transfer Q = 4πsin/λ with R+ and R− being the reflectivity for the neutron spin parallel or antiparallel to the external field, evidently confirms the magnetization . By simultaneously refining PNR and X-ray reflectivity data, the depth profiles of nuclear and magnetic scattering length densities at μ0H = 1 T, 0.8 T and 0.05 T for t = 24 u.c. were obtained and are shown in Fig. 2c. The uniform MSLD profile at the IP saturation field μ0H = 1 T attests to the high quality of the magnetic Cr2Te3 film with well-defined interfaces of 0.5 nm roughness. Remarkably, at reduced IP field μ0H = 0.8 T and 0.05 T, M develops a non-uniform depth-dependent profile with two distinct regions, possessing a lower IP magnetization value close to the substrate. Given that the NSLD depth profile of the Cr2Te3 layer is uniform and no changes are detected in the structure and chemical composition of the film, we attribute the reduced IP magnetization approaching the substrate to a canting of the magnetization vector towards the OOP direction .

Those differences could be explained by the diploid phasing of the Cabernet Sauvignon genome assembly and that multiple ISNT transcripts might correspond to a single gene locus. Nonetheless, similar relative amounts of Biological Process GO terms were found among the differentially expressed genes, confirming that the transcriptome obtained using Iso-Seq captured the transcriptional reprogramming underlying the main physiological and biochemical changes during grape berry development. In addition, gene expression analysis revealed that some private isoforms are significantly modulated during berry development, indicating that in addition to identifying the private gene space, the ISNT reference makes it possible to observe its expression. In conclusion, this study demonstrates that Iso-Seq data can be used to create and refine a comprehensive reference transcriptome that represents most genes expressed in a tissue undergoing extensive transcriptional reprogramming during development. In grapes, this approach can aid developing transcriptome references and is particularly valuable given diverse cultivars with private transcripts and accessions that are genetically distant from available genome references, like the non-vinifera Vitis species used as rootstocks or for breeding. The pipeline described here can be useful in efforts to reconstruct the gene space in plant species with large and complex genomes still unresolved.While San Joaquin Valley vineyards are currently fertilized with boron through the soil and foliage , some growers have expressed interest in applying boron via drip irrigation or “fertigation.” Fertigation is a relatively simple, cost-effective and efficient way to apply nutrients. However, irrigation water with more than 1 part per million boron can lead to vine toxicity, so the safety of boron fertigation is also a concern. Our research evaluates the safety and efficacy of boron fertigation in grapevines using drip irrigation. Boron is unique among the micronutrients due to the narrow range between deficiency and toxicity in soil and plant tissues. For grapevines, 25 liter plant pot this range is 0.15 ppm to 1 ppm in saturated soil extracts, and 30 ppm to 80 ppm in leaf tissue.

The goal of boron fertilization of grapevines is to keep tissue levels within this narrow range, since both deficiency and toxicity can have serious negative effects on vine growth and production. Fertilization amounts must be precise to avoid toxicity while providing adequate boron to satisfy grapevine requirements . On the east side of the San Joaquin Valley, boron deficiency of grapevines occurs on soils formed from igneous rocks of the Sierra Nevada. This parent material is low in total boron, which is crystallized in borosilicate minerals that are highly resistant to weathering. Boron deficiency is often associated with sandy soils and vineyard areas with excessive leaching, such as in low spots or near leaky irrigation valves. Vine symptoms of boron deficiency are more widespread and pronounced following high rainfall years, when greater amounts of soluble boron are leached from the root zone. In addition, snow melt water has very low levels of boron, and vineyards irrigated primarily with this water have a greater risk of deficiency. Boron is required for the germination and growth of pollen during flowering, and vines that are deficient in this micronutrient will have clusters that set numerous shot berries, small berries with a distinctive pumpkin shape. When boron deficiency is severe, vines produce almost no crop. Foliar symptoms appear in the spring: shoots have shortened, swollen internodes and their tips sometimes die, and leaves have irregular, yellowish mottling between the veins. Grapevines are also sensitive to too much boron. Toxicity is common on the west side of the San Joaquin Valley, where most soils are derived from marine sedimentary and meta sedimentary parent material that is rich in easily weathered boron minerals. Symptoms of boron toxicity include leaves that are cupped downward in the spring and that develop brown spots adjacent to the leaf margin in middle or late summer, intensifying and leading to necrosis as boron accumulates. Yields are reduced, the result of diminished vine vigor and canopy development. When foliar boron sprays are applied in excess in the spring, juvenile leaves become cupped within 2 weeks; however, vines quickly recover and yields are usually unaffected. Toxicity also occurs when boron fertilizer is applied in excess, regardless of the soil type, and this can lead to yield loss. Over-fertilization is the sole reason for boron toxicity on the east side of the San Joaquin Valley, so it is critical to establish how much boron fertilizer can be applied safely and effectively.

Our research investigated the uptake of boron by grapevines when fertilizer was applied with a drip-irrigation system.Research was conducted from 1998 to 2001 in a mature ‘Thompson Seedless’ raisin vineyard near Woodlake in Tulare County. The vineyard was planted in Cajon sandy loam on a recent alluvial fan associated with the Kaweah River. This soil is derived from granitic parent materials, and the surface soil is highly micaceous with a slight to moderate amount of lime. The underlying soil has a coarse, sandy texture. At the onset of this study, the vineyard’s boron status was in the questionable range for deficiency. The vine’s leaf petioles and blades contained about 30 ppm boron. While the foliage had no symptoms of boron deficiency, in the past the grower had observed sticking caps and pumpkin-shaped shot berries, which are indicative of boron deficiency. During the course of the research, the vineyard was drip-irrigated from April through October. The vineyard canopy covered 60% of the land surface during summer months and about 20 inches of water was applied during the season. Boron treatments consisted of applying fertilizer in varying amounts 3 weeks prior to bloom on May 18, 1998, and then again 3 weeks prior to bloom the following year, on May 3, 1999. Growers who fertigate grapevines with a drip system generally inject material into the irrigation water over a 45-to- 60-minute period at the beginning of an irrigation set. We simulated fertigation by applying Solubor, a soluble boron product , to a shovel-sized hole beneath drippers during the first hour of the irrigation set. By doing this, precise amounts of boron could be applied to each plot and plot size could be reduced. This technique has been used successfully in previous research with other nutrients . The experiment was designed as a randomized complete block with five treatments, five blocks and five vine plots . To evaluate the rate of boron uptake and accumulation in tissue with consecutive years of fertilization, grape tissue samples were collected in 1998 and 1999 at bloom and then again about 6 weeks later during veraison. Veraison is the stage of development where berries begin to soften and/or color. To evaluate carryover, leaf tissue samples were also collected at this Tulare County site at both bloom and veraison in 2001, 2 years after the fertilization was discontinued. In each case, 100 petioles and 50 blades were sampled per plot from the center three vines. Petioles and blades were taken opposite inflorescences during bloom, and recently matured leaves were sampled at veraison. Samples were oven-dried, ground in a Wiley mill and sent to the UC Davis DANR Analytical Laboratory for analysis of total boron. Statistical analysis was by ANOVA using least significance difference to separate treatment means. A second experiment was conducted in 1998 in Fresno County near Selma, in a mature Thompson Seedless raisin vineyard planted on Pollaski sandy loam and drip-irrigated.

The soil was formed in place from the weathering of softly to moderately consolidated granitic sediments. The particle size distribution of the surface soil is 63% sand, 25% silt and 12% clay. At the onset of the experiment, boron tissue levels were in the adequate range, 40 ppm. In both experiments, drip irrigations during the season were based on a schedule using historical evapotranspiration and developed for raisin vineyards in the San Joaquin Valley . The irrigation source was high-quality pump water with a boron concentration less than 0.1 ppm. The experimental design and methods were identical in both vineyards, except that the 1/16-poundper-acre boron treatment was omitted in the second vineyard. The Fresno County trial was discontinued after tissue samples were taken at bloom and veraison in 1998.At both the Tulare and Fresno county sites, black plastic plant pots boron uptake was rapid when fertilizer was applied in the spring. In both vineyards, applying boron at 2/3 or 1 pound per acre increased the boron concentration in blades by bloom, 3 weeks after application. Boron increased further in blades by veraison . In the Tulare County vineyard, boron in bloom tissue increased from a questionable deficiency range to adequate; at the Fresno location, boron in bloom tissue increased from 40 ppm to 54 ppm, a dramatic increase considering boron fertilizer was applied just 3 weeks prior. This indicates that boron uptake is rapid. None of the fertigation treatments resulted in either symptoms of boron toxicity or deficiency. Applying boron at 1/3 pound per acre or less did not significantly increase boron tissue levels by bloom or veraison at either site the first year. Fertigation over consecutive years was evaluated at the Tulare County location. Boron in grapevine tissue continued to increase with consecutive years of application. At the higher fertilizer rate , boron levels in blades increased from 35 ppm in control vines to about 60 ppm. We speculate that continuing with annual applications of 1 pound boron per acre would result in toxicity within 4 to 5 years. The 1/3-pound-per-acre rate significantly elevated boron in blades by veraison of the second year to adequate levels . There were no visual signs of toxicity in any of the fertilizer treatments, even when boron was applied at 2 pounds per acre in a single application. Boron levels in tissue remained unchanged 2 years after fertilization was discontinued at the Tulare County location . This indicates substantial treatment longevity with fertigation of a drip-irrigated vineyard. Rainfall during this experiment was below normal, which helped minimize leaching. Also, well-managed drip irrigation minimizes leaching. Under drip irrigation, salts tend to accumulate near the soil surface and 2 to 3 feet away from the drip line, with minimal water and salt movement below the root zone when irrigations are accurately scheduled . Boron concentrated more in the blades than in the petioles in response to fertilization. At the onset of the Tulare County experiment, boron concentrations in petioles and blades were similar at 31 ppm and 34 ppm, respectively. Fertilizing with 1 pound boron per acre for 2 consecutive years resulted in a 25% increase of boron in petioles but a 76% increase in blades . All fertilizer treatments increased boron in blades more than in petioles, indicating that blades should be sampled when monitoring the vines’ boron status following fertilization. Potential boron toxicity values at the time of sampling during the bloom period are 80 ppm for petioles and 120 ppm for blades, and in mid- to late summer are 100 ppm for petioles and 300 ppm for blades.Annual boron fertigation at 1/3 pound per acre elevated grapevine tissue levels from questionable to the adequate range within 2 years . In addition, tissue boron levels remained unchanged 2 years after fertilization was discontinued. This is probably because leaching was reduced by two factors: below-normal rainfall and accurately scheduled drip irrigations. After fertilization, boron was concentrated more in blades than in petioles, indicating that blades are the best choice for monitoring toxicity. Blade samples should be monitored on a routine basis and fertilizer amounts should be adjusted accordingly to avoid boron toxicity or deficiency. The results of this research can be applied to other drip-irrigated vineyards in the San Joaquin Valley under similar conditions: rapidly drained soils, high quality irrigation water, and low boron content in soil, water and vine tissue. In other regions of the state where winter rainfall is much higher, there would presumably be more leaching of boron fertilizer during winter months and less carryover time after fertilization is discontinued. In contrast, less leaching and greater carryover of boron would be expected in areas of less rainfall or on soils with finer texture and higher water-holding capacity. These variables underscore the importance of monitoring boron in tissue when developing a long-term fertilization program.The study of crop domestication has long been used as a proxy for studying evolutionary processes, such as the genetic effects of bottlenecks and the detection of selection to identify agronomically important loci .

To further explore links between management and soil fertility, we used the results from the PCA to formalize a gradient in management across all farms, and then used this gradient as the basis for comparison between Field A and Field B across all indicators for soil fertility. Using the ggplot and tidy verse packages , we displayed the difference in values between Field A and Field B for each indicator for soil fertility sampled at each farm using bar plots. We also included error bars to show the range of uncertainty in these indicators for soil fertility. Lastly, we further compared Field A and Field B for each farm using radar plots. To generate the radar plots, we first scaled each soil indicator from 0 to 1. Using Jenks natural breaks optimization, we then grouped each farm based on low, medium, and high N-based fertilizer application, as this soil management metric was the strongest coefficient loading from the first principal component . Using the fmsb package in R , we used an averaging approach for each level of N-based fertilizer application to create three radar plots that each compared Field A and Field B across the eight indicators for soil fertility. Farmers provided an overview of their farm operation, including farm size , the total number of crops each farm planted per growing season at the whole farm level, the types of crops planted in their field during the initial field visit , the type and amount of nitrogen-based fertilizer they applied on farm, and key aspects of soil health in their own words . Farm sizes ranged from 15 to 800 acres, with about one third of farms in the 15 – 50-acre range, drainage pot another third in the 100 – 450-acre range, and roughly a final third in the 500 – 800 acre-range. Farmers grew primarily summer crops, including tomato, a variety of cucurbits, strawberry, herbs, nightshades, root vegetables, and sunflower/safflower for oil.

Farmers reported applying a range of external N-based organic fertilizers, including fish emulsion, Wiserg , pelleted chicken manure, and seabird guano, at varying rates . On the low end, farmers applied <1 kg-N/acre, and on the high end, farmers applied 90 – 180 kg-N/acre per season. About a third of farmers applied 2 – 25 kg-N/acre of N-based fertilizer. Farmer responses for describing key aspects of soil health were relatively similar and overlapped considerably in content and language . Specifically, farmers usually emphasized the importance of maintaining soil life and/or soil biology, promoting diversity, limiting soil compaction and minimizing disturbance to soil, and maintaining good soil structure and moisture. Several farmers also touched on the importance of using crops as indicators for monitoring soil health and the importance of limiting pests and disease. Discussion of the importance of promoting soil life, soil biology, and microbial and fungal activity had the highest count among farmers with ten mentions across the 13 farmers interviewed. Next to this topic, minimizing tillage and soil disturbance was the second most discussed with six of 13 farmers highlighting this key aspect of soil health. In addition to discussing soil health more broadly, farmers also provided in-depth responses to a series of questions related to soil fertility—such as key nutrients of interest on their farm, details about their fertility program, and the usefulness of soil tests in their farm operation— summarized in Table 2. When asked to elaborate on the extent to which they considered key nutrients, a handful of farmers readily listed several nutrients, including nitrogen, phosphorous, potassium , and other general macronutrients as well as one micronutrient .

Among these farmers that responded with a list of key nutrients, some talked about having their nutrients “lined up” as part of their fertility program. This approach involved keeping nutrients “in balance,” such as for example, monitoring pH to ensure magnesium levels did not impact calcium availability to plants. These farmers also emphasized that though nitrogen represented a key nutrient and was important to consider in their farm operation, tracking soil nitrogen levels was less important than other aspects of soil management, such as promoting soil biological processes, maintaining adequate soil moisture and aeration, or planting cover crops regularly. As one farmer put it, “if you add nutrients to the soil, and the biology is not right, the plants will not be able to absorb it.” Or, as another farmer emphasized, “It’s not about adding more [nitrogen]… I try to cover crop more too.” A third farmer emphasized, that “I don’t use any fertilizers because I honestly don’t believe in adding retroactively to fix a plant from the top down.” This same farmer relied on planting a cover crop once per year in each field, and discing that cover crop into the ground to ensure his crops were provided with adequate nitrogen for the following two seasons. While most farmers readily listed key nutrients, several farmers shifted conversation away from focusing on nutrients. These farmers generally found that this interview question missed the mark with regards to soil fertility. One farmer responded, “I’m not really a nutrient guy.” This same farmer added that he considered [soil fertility] a soil biology issue as much as a chemistry issue.” The general sentiment among these farmers emphasized that soil fertility was not about measuring and “lining up” nutrients, but about taking a more holistic approach. This approach focused on facilitating conditions in the soil and on-farm that promoted a soil-plant-microbe environment ideal for crop health and vigor. For example, the same farmer quoted above mentioned the importance of establishing and maintaining crop root systems, emphasizing that “if the root systems of a crop are not well established, that’s not something I can overcome just by dumping more nitrogen on the plants.” Another farmer similarly emphasized that they simply created the conditions for plants to “thrive,” and “have pretty much just stepped back and let our system do what it does; specifically, we feed our chickens whey-soaked wheat berries and then we rotate our chickens on the field prior to planting.

And we cover crop.” A third farmer also maintained that their base fertility program—a combination of planting a cover crop two seasons per year, an ICLS chicken rotation program, minimal liquid N-based fertilizer addition, and occasionally compost application—all worked together to “synergize with biology in the soil.” This synergy in the soil created by management practices—rather than focusing on nutrient levels—guided this farmer’s approach to building and assessing soil fertility on-farm. Another farmer called this approach “place-based” farming. This particular farmer elaborated on this concept, saying “I think the best style of farming is one where you come up with a routine [meaning like a fertility program] that uses resources you have: cover crops, waste materials beneficial to crops, animals” in order to build organic matter, which “seems to buffer some of the problems” that this farmer encountered on their farm. Similar to other farmers, this farmer asserted that adding more nitrogen-based fertilizer did not lead to better soil fertility or increase yields, drainage planter pot in their direct experience. Regardless of whether farmers listed key nutrients, a majority of farmers voiced that nitrogen was not a big concern for them on their farm. This sentiment was shared among most farmers in part because they felt the amount of nitrogen additions from fertilizers they added were insignificant compared to nitrogen additions by conventional farms. Farmers also emphasized that the amount of nitrogen they were adding was not enough to cause environmental harm; relatedly, a few farmers noted the absurdity and added economic burden of the recent nitrogen management plan requirements—specifically among organic farms with very low N-based fertilizer application. The majority of farmers also expressed that their use of cover crops and the small amount of N-based fertilizer additions as part of their soil fertility program ensured on-farm nitrogen demands were met for their crops. Across all farmers interviewed, cover cropping served as the baseline and heart of each fertility program, and was considered more effective than additional N-based fertilizers at maintaining and building soil fertility. Farmers used a range of cover crop species and often applied a mix of cover crops, including vetches and other legumes like red clover and cowpea , grains and cereals like oats . Farmers cited several reasons for the effectiveness of cover cropping, such as increased organic matter content, more established root systems, greater microbial activity, better aeration and crumble in their soils, greater number of earthworms and arthropods, improved drainage in their soils, and more bio-available N. Whereas farmers agreed that “more is not better” with regards to N-based fertilizers, farmers did agree that allocating more fields for planting cover crops over the course of the year was beneficial in terms of soil fertility. However, as one farmer pointed out, while cover crops provide the best basis for an effective soil fertility program, this approach is not always economically viable or physically possible. Several farmers expressed concern because they often must allocate more fields to cover crops than cash crops in any given season, which means that their farm operation requires more land to be able to produce the same amount of vegetables than if they had all their fields in cash crops.

Farmers also shared that in some circumstances, such as in early spring, they are not able to realize the full potential of a winter cover crop if they are forced to mow the cover crop early to plant cash crops and ensure the harvest timeline of a high-value summer vegetable crop. The cover crop approach to soil fertility takes “persistence,” as one farmer emphasized; another farmer similarly pointed out that the benefits of cover cropping “are not always realized in the crop year. You’re in it [organic agriculture] for the long haul, there is no quick fix.” Indeed, farmers who choose to regularly plant cover crops to build soil fertility, rather than just add N-based fertilizers, reported that they came up against issues of land tenure and access to land, market pressures, and long-term economic sustainability. To build on conversations about soil fertility, farmers also provided responses to interview questions that asked them to elaborate on the usefulness of available soil tests to gauge soil fertility more broadly—and then more specifically, the usefulness of soil tests in informing their soil fertility program and/or management approaches on-farm. Overall, only three of 13 farmers reported regularly using and relying on soil tests to inform their soil fertility program or aspects of their farm operation. These farmers offered very short responses and did not elaborate. For example, one farmer shared that they “test twice a year in general,” and that they “rely on the results of the soil tests to tweak [their] fertility program.” Another farmer said briefly, “We use soil tests… we utilize them to decide what to do to try to improve the soil.” A third farmer admitted that though he “used to do a soil test every year, literally used to spend hundreds of dollars per year on soil tests,” he found that the results of soil tests did not change year-to-year and were, as he put it, very “stable.” This particular farmer no longer regularly uses or relies on soil testing for their farm operation. The remaining ten farmers confirmed that they had previously submitted a soil test, usually once and most often to a local commercial lab in the region. These farmers expressed a range of sentiments when asked about the usefulness of soil tests, including disappointment, distrust, or both, particularly in the capacity of soil tests to inform soil fertility on their farm. Some farmers said directly, “I just don’t trust soil tests,” or “frankly, I don’t believe a lot in soil testing because it’s too standardized,” while other farmers initially stated they had used “limited” or “infrequent” soil tests, and then later admitted that they did not use or rely on soil tests on their farm operation. These farmers tended to focus on the limitations of soil tests that they encountered for their particular farm application. Limitations of soil tests discussed by farmers varied.

Across all farmers interviewed, including both first- and second-generation farmers, farmers stressed the steep learning curves associated with learning to farm alternatively and/or organically. While these farmers represent a case study for building a successful, organic farm within one generations, the results of this study beg the question: What advancements in farm management and soil management could be possible with multiple generations of farmer knowledge transfer on the same land? Rather than re-learning the ins and outs of farming every generation or two, as new farmers arrive on new land, farmers could have the opportunity to build on existing knowledge from a direct line of farmers before them, and in this way, potentially contribute to breakthroughs in alternative farming. In this sense, moving forward agriculture in the US has a lot to learn from agroecological farming approaches with a deep multi-generational history . To this end, in most interviews—particularly among older farmers—there was a deep concern over the future of their farm operation beyond their lifetime. Many farmers lamented that no one is slated to take over their farm operation and that all the knowledge they had accumulated would not pass on. There exists a need to fill this gap in knowledge transfer between shifting generations of farmers in order to safeguard farmer knowledge and promote adaptations in alternative agriculture into the future.Most studies often speak to the scalability of approach or generalizability of the information presented. While aspects of this study are generalizable particularly to similar farming systems in California such as the Central Coast region, the farmer knowledge presented in this study is not generalizable and not scalable to other regions in the US. To access farmer knowledge, pots with drainage holes relationship building with individual farmers leading up to interviews as well as the in-depth interviews themselves require considerable time and energy.

While surveys often provide a way to overcome time and budget constraints to learn about farmer knowledge, this study shows that to achieve specificity and depth in analysis of farmer knowledge requires an interactive approach that includes—at a minimum—relationship building, multiple field visits, and in-depth, multi-hour interviews. Accessing farmer knowledge necessitates locally interactive research; this knowledge may not be immediately generalizable or scalable without further locally interactive assessment in other farming regions. Local knowledge among farmers in US alternative agriculture has often been dismissed or overlooked by the scientific community, policymakers, and agricultural industry experts alike; however, this study makes the case for inclusion of farmer knowledge in these arenas. In-depth interviews established that farmers provide an important role in translating theoretical aspects of agricultural knowledge into practice. It is for this reason that farmer knowledge must be understood in the context of working farms and the local landscapes they inhabit. As one of the first systematic assessments of farmer knowledge of soil management in the US, this research contributes key insights to design future studies on farmer knowledge and farmer knowledge of soil. Specifically, this study suggests that research embedded in local farming communities provides one of the most direct ways to learn about the substance of farmer knowledge; working with the local UCCE advisor in combination with community referrals provided avenues to build rapport and relationships with individual farmers—relationships that were essential to effective research of farmer knowledge. Farmer knowledge of soil management for maintaining healthy soils and productive, resilient agriculture represents an integral knowledge base in need of further scientific research. This study provides a place-based case study as a starting point for documenting this extensive body of knowledge among farmers.

It is our hope that this research will inspire future studies on farmer knowledge in other contexts so that research in alternative agriculture can widen its frame to encompass a more complete understanding of farming systems and management motivations—from theory to practice. A fundamental challenge in agriculture is to limit the environmental impacts of nitrogen losses while still supplying adequate nitrogen to crops and achieving a farm’s expected yields . To balance among such environmental, ecological, and agronomic demands, it is essential to establish actual availability of nitrogen to crops . A holistic, functional understanding of plant N availability is particularly imperative in organic agriculture, as in this farming context, synthetic fertilizers are not applied and instead, production of inorganic N—the dominant form of N available to crops—depends on internal soil processes . In organic agricultural systems, farmers may seasonally apply cover crops or integrate livestock as alternative sources of nitrogen to crops—in addition to or in place of using organic fertilizers. In applying these alternative sources of nitrogen to soil, organic farmers rely on the activity of soil microbes to transform organic N into inorganic forms of N that are more readily available for crop uptake . Currently, the predominant way crop available N is measured in organic agricultural systems tends to examine pools of inorganic N in soil . Inorganic N, or more specifically ammonium and nitrate , represents the predominant forms of N taken up by crop species in ecosystems where N is relatively available, such as in non-organic agricultural systems that apply inorganic fertilizers . However, in organic systems, crop available N is largely controlled by complex soil processes not adequately captured by simply measuring pools of ammonium and nitrate. First, because nitrogen made available to crops is controlled by soil microbes—wherein crops only have access to inorganic forms of N after microbial N transformations occur to first meet microbial N demand—pinpointing the flow of N moving through inorganic N pools as a result of these microbial N transformations is necessary to accurately measure actual N availability to crops . Second, extensive recycling of N among components of the plant-soil-microbe system complicates relying solely on measurements of inorganic N pools, which do not reflect these dynamics .

As an example, one previous study in organic vegetable systems showed examples where inorganic N pool sizes in the soil were measured to be low, yet there existed high production and consumption rates of inorganic N . This outcome highlighted that if the turnover of inorganic N is high—for instance, high rates of soil ammonium production exist in the soil with simultaneously high rates of immobilization by soil microbes and high rates of uptake by plants—measured pools of inorganic N may still be low . This study also showed that conversely, there may also exist situations when inorganic N pools are low and rates of ammonium and nitrate production are also low, in which case N availability would limit crop production. In organic systems especially, higher carbon availability as a result of organic management can increase these microbially mediated gross N flows, thereby increasing N cycling and turnover of inorganic N . Thus, we hypothesize that measuring total production of ammonium from organic N, or gross N mineralization, and subsequent total production of nitrate from ammonium, or gross N nitrification, may provide a more complete characterization of crop available N in the context of organic systems . Though the application of such diverse management practices on organic farms is known to affect rates of N cycling in soil , drainage pot measuring N flow rates as a proxy for crop available N is currently uncommon on working organic farms. The current historical emphasis on measuring inorganic pools of N in organic agriculture was originally imported from non-organic farming, wherein the Sprengel-Liebig Law of the Minimum was a widely accepted agronomic principle . In practice, this Law of the Minimum placed particular importance on using artificial fertilizers to overcome so-called “limiting” nutrients—namely, inorganic forms of N. Because inorganic N is relatively straightforward to measure, focus on quantifying pools of inorganic N has since become common practice among agronomists and agricultural researchers . However, the continued acceptance of the Law of the Minimum in organic agriculture underscores the gap in a functional understanding of organic agricultural systems, in particular the role of soil microbes in mediating N cycling. To understand crop available N more holistically, there is a need to measure actual flow rates of soil N—in addition to—static pools of inorganic N . Soil indicators that adequately capture N availability to crops are therefore necessary to move beyond the legacy of the Law of the Minimum in organic agriculture. Unpacking the soil processes that mediate flows of N may ultimately provide a more accurate characterization of soil N cycling and in turn, N availability to crops. Unfortunately, gross N mineralization and nitrification rates are very difficult to measure in practice, particularly on working organic farms . While net N flows are easier to measure in comparison to gross N flows and can provide a useful measure of N cycling dynamics as a complement to measurements of inorganic N pools, net N flows still pose serious limitations— namely that net rates cannot detect plant-soil-microbe interactions and therefore are not adequate as metrics for determining crop available N . In particular, relying on net N flows as a measure of N availability does not account for the ability of plants to compete for inorganic N, and assumes plants take up inorganic N only after microbial N demands are satisfied .

It is also possible that measuring soil organic matter pools could help indicate N availability because SOM supports microbial abundance and activity, and because SOM is also the source of substrates for N mineralization . Several studies have proposed measuring soil organic matter levels to complement measuring inorganic N pools, understand soil N cycling, and infer N availability . Assessing the total quantity of organic carbon and nitrogen within soil organic matter represents one established method for measuring levels of soil organic matter, and is more readily measurable than gross N rates. Additional indicators for quantifying “labile” pools of organic matter, such as POXC and soil protein, have also become more widely studied in recent years, and applied on organic farms as well . When used in combination with more established soil indicators that measure organic C and N pools , this suite of indicators may potentially provide added insight to understanding crop available N . Importantly, applied together these four indicators for soil organic matter levels may also more readily and accurately serve as a proxy for soil quality—generally defined as a soil’s ability to perform essential ecological functions key to sustaining a farm operation . Despite the availability of these soil indicators, very few studies have systematically examined the way in which SOM levels on working farms compare to N cycling processes, and specifically how SOM levels compare to microbially mediated gross N rates. Further, it is still unclear to what degree the interactions between soil edaphic characteristics and soil management influence N cycling and N availability to crops . For instance, soil texture may play a mediating role in N cycling, where soils high in clay content may limit substrate availability as well as access to oxygen, which in turn, may restrict the efficiency of N cycling . In this sense, it is important to understand the role that soil edaphic characteristics play in order to identify the underlying baseline limits imposed by the soil itself. Equally important to consider is the role of soil management in mediating N cycling. Compared to controlled experiments, soil management regimes on working farms can be more complex and nonlinear in nature due to multiple interacting practices applied over the span of several years, and even multiple decades. To date, a handful of studies conducted on working farms have examined tradeoffs among different management systems , though few such studies examine the cumulative effects of multiple management practices across a gradient of working organic farms. However, understanding the cumulative effects of management practices is key to link soil management to N cycling on working farms . Likewise, it is important to examine the ways in which local soil edaphic characteristics may limit farmers’ ability to improve soil quality through management practices. Though underutilized in this context, the development of farm typologies presents a useful approach to quantitatively integrate the heterogeneity in management on working organic farms . Broadly, typologies allow for the categorization of different types of organic agriculture and provide a way to synthesize the complexity of agricultural systems . Previous studies that make use of farm typologies found that differences in total soil N across farms are largely defined by levels of soil organic matter.

Blue elderberries were determined to be ripe when the berries in a cyme were deep purple, with or without the white bloom, and had no green berries present. Ripe elderberries were harvested by hand from all four quadrants of the elderberry shrub, totaling approximately 3 kg of elderberries. The berries were placed in clear plastic bags, stored on ice, and transported to the laboratory. A subsample was separated for moisture analysis, while the rest was de-stemmed and stored at -20 °C until analyzed. HPLC grade methanol , acetonitrile , phosphoric acid, ethanol, hydrochloric acid, and sodium hydroxide were purchased from Fisher Scientific . Ascorbic acid was purchased from Acros Organics . Formic acid, gallic acid, sucrose, chlorogenic acid, rutin, and catechin, and Folin-Ciocalteu reagent were purchased from Sigma Aldrich . Cyanidin-3-glucoside was purchased from Extrasynthese . Ultra-pure water was obtained from a Milli-Q water system . From each shrub, about 250 g of frozen berries were thawed in a glass container overnight at 4 °C. The following day, the thawed berries were mashed for 4 min by hand using a plastic pestle, then homogenized , and centrifuged at 3,000 rpm for 7 min. The supernatant was strained, collected, and weighed. A 15 mL aliquot of sample extract was stored in a 15 mL plastic tube at -20 °C for total monomeric anthocyanin analysis. The remaining supernatant was used to determine soluble solids, pH, and titratable acidity . One pooled sample was analyzed from each shrub and analyzed in duplicate analytical repetitions. A refractometer was used to determine soluble solids. It was calibrated with standard solutions of 5°, 10°, vertical gardening in greenhouse and 15° Brix made with sucrose and water. A 150 µL aliquot of elderberry juice was placed on the prism and read using the automatic setting. The pH was determined using a SevenMulti pH meter .

It was calibrated before each use using buffers at pH 4.0, 7.0, and 10.0. To determine TA, 10 mL of elderberry juice was diluted to 100 ml with nanopure water and mixed. This dilute juice was titrated to pH 8.2 using 0.1 N NaOH. The volume of 0.1 N NaOH used to achieve the desired pH was used to calculate the mg citric acid per 100 g fresh weight . For each of these analyses, duplicates were run on each juice. Elderberries were extracted by combining 5 g frozen berries with 25 mL MeOH:formic acid in a conical tube. The contents were homogenized, placed in a shaker without water at speed 7.5 for 20 min, then centrifuged at 3,000 rpm for 7 min. The supernatant was transferred to a 15 mL plastic tube and stored at -80 °C for no more than two weeks prior to analysis. Duplicate extracts were made from each shrub. TPC was determined using the Folin-Ciocalteu method. First, elderberry phenolic extract was diluted 1:4 with water. Each extract was analyzed in duplicate and averaged. In 10 mL glass tubes, 6 mL water was combined with 100 µL sample and 500 µL Folin-Ciocalteu reagent. After mixing and incubating for 8 min at room temperature, 1.5 mL 20% aqueous sodium carbonate was added. The tubes were mixed, covered with foil to avoid light exposure, placed in a water bath at 40 °C for 40 min, then cooled at room temperature for 15 min. The samples were read by a UV visible spectrophotometer at 765 nm and quantified using an external standard curve prepared with gallic acid . TPC is expressed as mg gallic acid equivalents per 100 g FW. Five grams of frozen berries were mixed with 25 mL of in a conical tube, which was then homogenized for 1 min at 7,000 rpm . The mixture was stored at 4 °C overnight, then in the morning, centrifuged at 4,000 rpm for 7 min. The supernatant was used directly for analysis. Three pooled samples were made for each hedgerow, each consisting of even amounts of berries from three distinct shrubs. Eachpooled sample was extracted once to give 3 biological replicates, and each extract was run in duplicate . The concentration of phenolic compounds in blue elderberry followed the method by Giardello et al. with some modifications.

Briefly, samples were analyzed via reversed-phase liquid chromatography on an Agilent 1200 with a diode array detector and fluorescence detector . The column used was a PLRP-S 100A 3 µm 150 x 4.6 mm at 35 °C, and the injection volume was 10.0 µl. Mobile phase A was water with 1.5 % phosphoric acid, while mobile phase B was 80%/20% acetonitrile/ mobile phase A. The gradient used was 0 min 6% B, 73 to 83 min 31% B, 90 to 105 min 6% B. The DAD was used to monitor hydroxybenzoic acids at 280 nm, hydroxycinnamic acids at 320 nm, flavonols at 360 nm, and anthocyanins at 520 nm. The FLD was used to monitor flavan-3-ols, with excitation at 230 nm and emission at 321 nm. External calibration curves were prepared using chlorogenic acid for phenolic acids, rutin for flavonols, and cyanidin-3-glucoside for anthocyanins , at the following concentrations: 200, 150, 100, 75, 50, 25, 10, 5, and 2.5 mg/L. Catechin was used to quantify flavan-3-ols and standards were run at 150, 100, 75, 50, 25, 10, 5, 2.5, and 1 mg/L. Compounds were identified based on retention time and spectral comparisons with standards. Information about the linear equations and lower limits of detection and quantitation can be found in Table S1 in the supplementary material. The LLOD was calculated as 3.3 times the standard deviation of the y-intercept of the curve divided the slope, while the LLOQ was calculated as 10 times those values.Several peaks appeared in the HPLC chromatograms that could not be identified using the above parameters. Chromatographic eluents of these peaks were collected individually and dried under vacuum. These extracts were reconstituted with mobile phase A, and 5 µL were injected into the HPLC- QTOF-MS/MS for accurate mass analysis . A Poroshell 120 EC-C18 column was used at 35 °C. Mobile phase A was 1% formic acid in distilled water, and mobile phase B was 1% formic acid in acetonitrile. The gradient used was 0 min 3% B, 30 min 50% B, 31-32 min 95% B, 33-38 min 3% B. The mass spectrometer was used in negative mode, and the mass range for MS was 100 to 1000 m/z while the range for MS/MS was 20-700 m/z. Collision energies at 10, 20, and 40 V were applied. The drying gas was set to a flow of 12 L/min at 250 °C, while the sheath gas was set to 11 L/min at 350 °C. The nebulizer was set to 40 psig, the capillary voltage was 3500 V, the nozzle was set to 500 V, and the fragmentor was set to 100 V. Data was analyzed using Agilent MassHunter Workstation Qualitative Analysis 10.0 .

Tentative identification was achieved by comparing the mass to charge ratio of the precursor and fragment ions to online libraries of compounds as well as using formula generation for the peaks in the spectra. The composition of blue elderberries is presented for the first time, which is key to understanding how this subspecies of Sambucus nigra compares to commercialized elderberry subspecies, S. nigra ssp. nigra and S. nigra ssp. canadensis. These data help to establish the blue elderberry grown in hedgerows in California as a viable source of berries and bioactive compounds. Data for the compositional assays is presented for the 2018 and 2019 harvest years as the average of all shrubs sampled in Table 2. The average moisture for the blue elderberries was 79.5 ± 1.5% in 2018 and 79.5 ± 1.6% in 2019, which is very similar to the levels found in wild elderberries in Spain 95. The average soluble solids found in blue elderberry ranged from 11.94 ± 2.08 to 14.95 ± 1.02 g per 100 g FW in 2018 and from 12.64 ± 1.86 to 17.09 ± 1.60 g per 100 g FW in 2019. These values are slightly higher than the soluble solids found in S. nigra ssp. cerulea grown in Slovenia29 and American elderberries grown in Ohio52. Compared to European and American elderberries evaluated in other studies, blue elderberries have similar levels of soluble solids 8,18,29,49,50,95. In the present study, the overall average content of soluble solids was significantly different between years, greenhouse vertical farming as blue elderberries harvested in 2019 had significantly higher average soluble solids than the elderberries harvested in 2018 . The pH in the blue elderberry ranged from 3.44 to 3.86 in 2018 and from 3.46 to 3.79 in 2019, with no significant difference found between harvest years. These values are slightly lower than the values found in European elderberry, which ranged from 3.9 ± 0.06 to 4.1 ± 0.04 with an average pH of 3.9 ± 0.2, and American elderberry, which ranged from 3.9 ± 0.04 to 4.5 ± 0.03 with an average pH of 4.2 ± 0.2 49 Another evaluation of pH in American elderberries had a range of 4.5 ± 0.08 to 4.9 ±0.12, higher than those found in the blue elderberry.52 The higher sugar and lower pH levels in blue elderberry could potentially impact taste and performance in food and beverages as compared with the European and American species. The average titratable acidity in blue elderberries ranged from 0.45 ± 0.08 to 0.77 ± 0.03 g citric acid per 100 g FW in 2018 and from 0.54 ± 0.06 to 0.77 ± 0.11 g citric acid per 100 g FW in 2019 with no significant difference found between harvest years. These values are lower than the total acids found by Mikulic-Petkovsek et al. 29 in S. nigra ssp. cerulea , but they are similar to the levels found in European elderberry 8,18,49,50 . Anthocyanins are a class of phenolics that contribute red, purple, and blue hues to fruits and vegetables, act as attractants for pollinators, and are potent antioxidants. European and American elderberries are well-known for containing high levels of anthocyanins 8,18,49. The anthocyanin content of elderberries strongly correlates to the antioxidant potential of the fruit, which may confer health-promoting properties 50,89, which is one reason why elderberries are used in supplements and value-added products. Elderberry is also used as a source of natural food colorants due to the levels of anthocyanins35. Understanding the levels of anthocyanins in the blue elderberry grown in hedgerows is critical towards establishing this native fruit as an additional and more sustainable elderberry. The average TMA measured in blue elderberry ranged from 34.2 ± 9.7 to 113.4 ± 18.2 mg CGE per 100 g FW in 2018 and from 43.1 ± 11.5 to 121.5 ± 11.5 mg CGE per 100 g FW in 2019 . TMA was variable between hedgerows in both years of harvest, with relative standard deviation values between 16% and 30%, yet there was not a significant difference in the overall average TMA between 2018 and 2019 . Furthermore, most hedgerows were not significantly different from the other hedgerows harvested that year despite significant differences in TMA values found between farms in both years . Regarding the age of the elderberry shrub, hedgerows 2 and 14 had two of the three highest concentrations of TMA in 2019 . This suggests that blue elderberries can be harvested from plants as young as two years without a significant loss of TMA concentrations. TMA values for the blue elderberries are lower than those found in other elderberry subspecies. In European elderberries, TMA levels range from 170 ± 12 to 343 ± 11 with an average of 239 ± 94 mg CGE per 100 g FW 49. A study of American elderberry grown in Ohio showed a range from 354 ± 59 to 595 ± 26 mg CGE per 100 g FW.52 In the present study, bare root prerooted cuttings of American elderberries were planted, along with blue elderberries, on Farm 1 in 2018, and three shrubs were harvested in 2019. These American elderberries had an average TMA value of 263 ± 5.4 mg CGE per 100 g FW, which is more similar to what has been observed in other studies on this subspecies. This suggests it is a subspecies difference contributing to the lower anthocyanin concentration in the blue elderberry and not the difference in growing conditions.