Coordinated Regulation of Biosynthetic and Regulatory Genes Coincides with Anthocyanin Accumulation in Developing Eggplant Fruit
Click on author name to view affiliation information
Click on author name to view affiliation information
Violet to black pigmentation of eggplant (Solanum melongena L.) fruit is caused by anthocyanin accumulation. Model systems demonstrate the role of regulatory genes in the control of anthocyanin biosynthesis. Anthocyanin structural gene transcription requires the expression of at least one member of each of three transcription factor families: MYB, MYC, and WD. To determine the molecular genetic basis for anthocyanin pigmentation in eggplant fruit, we used real-time polymerase chain reaction (PCR) to evaluate the expression of anthocyanin biosynthetic (Chs, Dfr, Ans) and regulatory (Myc, MybB, MybC, Wd) genes in S. melongena genotypes that produce fruit with dark violet (‘Classic’) or white (‘Ghostbuster’) coloration, respectively. Transcript levels and anthocyanin content were evaluated in fruit at various stages of development ranging from small post-anthesis fruit to full-sized marketable fruit. Anthocyanin content increased 9-fold in developing violet-colored ‘Classic’ fruit, whereas low but detectable concentrations were found in white ‘Ghostbuster’ fruit. Chs, Dfr, and Ans as well as MybC and Myc transcript levels were significantly higher in ‘Classic’ in comparison with ‘Ghostbuster’ fruit at comparable stages of fruit development with greatest differences observed for Ans transcript levels. MybC and Myc transcript levels increased in developing ‘Classic’ fruit coincident with increasing anthocyanin content. MybB and Wd transcript levels were not coordinated with changes in biosynthetic transcript levels or anthocyanin concentration.
Eggplant fruit color is influenced by the presence or absence of anthocyanins and chlorophylls and the concentration as well as distribution of these compounds. Anthocyanins are located in fruit epidermal cells and chlorophylls are found in subepidermal cell layers. Cis- and trans-isomers of the acylated anthocyanin delphinidin-3-p-coumaroylrutinoside-5-glucoside are the predominant anthocyanins commonly reported in eggplant cultivars (Ichiyanagi et al., 2005). Genotypes containing the unacylated delphinidin-3-rutinoside as the major anthocyanin have also been characterized (Sadilova et al., 2006; Tanchev et al., 1970; Wu and Prior, 2005). Dependent on pigment concentrations, the presence of both anthocyanins and chlorophylls results in violet to black fruit color. In the absence of chlorophyll pigments in these developing fruit, varying violet hues occur and are dependent on anthocyanin concentration. Conversely, in the absence of anthocyanin pigments, or if their concentrations are very low, varying fruit chlorophyll concentration results in fruit color ranging from light green to dark green. White fruit color occurs in the absence or very low concentrations of both pigment classes. Concentration of both pigment classes is genotype-dependent and subject to environmental influence. Whereas Western breeders have focused on development of cultivars with dark purple to black- and green-colored fruit, the full range of possible fruit colors are valued in diverse world markets (Daunay et al., 2004).
Anthocyanins are secondary metabolites derived from the phenylpropanoid pathway with important physiological roles in plant tissues that also influence crop quality and nutritive value. In flowers, anthocyanins may function as attractants for pollinators and aid in the dispersal of seeds and fruit (Harborne, 1994; Schemske and Bradshaw, 1999). They also function in plant tissues as antioxidants and ultraviolet and visible light protectants (Close and Beadle, 2003; Gross, 1991). This class of compounds has generated considerable interest for their presumptive health-promoting effects in mitigating or preventing chronic diseases. Anthocyanins are considered functional dietary components with antiinflammatory and anticarcinogenic activity and have roles in cardiovascular disease prevention, obesity control, and diabetes alleviation properties (He and Giusti, 2010; Pojer et al., 2013).
Anthocyanin biosynthesis has been studied extensively. The enzymes in the anthocyanin biosynthetic pathway are functionally conserved across diverse plant species and the genes that encode these enzymes share high sequence similarity across species. Six enzymes are commonly involved in anthocyanin biosynthesis (Griesbach, 2005). Chalcone synthase (CHS) is the first and key regulatory enzyme of flavonoid biosynthesis. Chalcone is isomerized by chalcone isomerase to naringenin and flavanone 3-hydroxylase converts naringen into dihydrokaempferol. The dihydroflavonols are subsequently converted to a colorless leucoanthocyanidin by dihydroflavonol 4-reductase (DFR), the first committed enzyme of anthocyanin biosynthesis. Leucoanthocyanidins are converted to colored anthocyanidins by anthocyanidin synthase (ANS). UDP-glucose-flavonoid 3-O-glucosyltransferase finally produces the anthocyanin-3-glucoside.
Tissue- or developmental-specific anthocyanin biosynthesis is controlled by a MYB-MYC-WD40 regulatory protein complex (Koes et al., 2005). A common set of proteins, comprised of 1) MYB proteins characterized by R2 and R3 imperfect repeats of the conserved DNA-binding motif; 2) MYC proteins defined by a conserved basic helix-loop-helix domain consisting of a presumed DNA binding basic region and a protein–protein dimerization motif; and 3) WD40 repeat proteins with a characteristic 40 residue core containing a glycine-histidine dipeptide at the N-terminus and a tryptophan-aspartate dipeptide at the C-terminus.
Gene families with multiple copies of anthocyanin biosynthetic and regulatory genes have been described within various species (Dooner et al., 1991; Du et al., 2012; Griesbach, 2005; Mol et al., 1996). Evolutionary studies indicate that these multiple copies likely arose from gene duplication followed by sequence divergence (Purugganan and Wessler, 1994). Resulting variation in function of these genes, particularly for regulatory genes, is responsible for differences in pigmentation and patterns of tissue- or developmental-specific pigmentation within and between species. We previously characterized tissue-specific regulation of the anthocyanin biosynthetic pathway in pepper (Capsicum annuum L.), demonstrating differential expression of biosynthetic and regulatory genes coincident with anthocyanin pigmentation in accessions with divergent pigmentation (Stommel et al., 2009). Here, we extend those findings on anthocyanin gene regulation to eggplant, a related Solanaceous species.
The S. melongena hybrid cultivar Classic bears fruit with characteristic purple pigmentation restricted to the fruit peel and the hybrid ‘Ghostbuster’ bears uniform white-colored fruit. ‘Classic’ and ‘Ghostbuster’ were produced from seed and grown at the Beltsville Agricultural Research Center, Beltsville, MD, during the summer (June to September) season under standard greenhouse conditions without supplementary lighting. Six 7-week-old plants of each cultivar were transplanted to 5.7-L pots containing a 65% to75% sphagnum peatmoss plus perlite soil-free mix (Pro-Mix HP Mycorrhizae; Premier Tech Horticulture, Quakertown, PA). Pots were distributed in a glasshouse (lat. 39°1′38.96″ N, long. 76°55′35.67″ W) with climate control (heating started at temperatures below 20 °C and cooling at temperatures above 24 °C) using a completely randomized design. Plants were trained on bamboo canes and spaced 0.8 m between rows and 0.7 m apart within the row. Fertilization was applied on a constant feed basis with drip irrigation throughout the growing cycle. Final concentrations of 144 mg·L−1 nitrogen, 42 mg·L−1 phosphorus, 221 mg·L−1 potassium, 106 mg·L−1 calcium, 50 mg·L−1 magnesium, and 34 mg·L−1 sulfur in the dilute fertilizer solution were supplemented with micronutrients applied in final concentrations of 2.8 mg·L−1 iron, 1.3 mg·L−1 boron, 0.18 mg·L−1 copper, 0.62 mg·L−1 manganese, 0.18 mg·L−1 zinc, 0.09 mg·L−1 molybdenum, and 8.8 mg·L−1 chlorine.
Fruit peel was sliced from developing fruit, frozen in liquid nitrogen, and stored at –80 °C. Replicate samples were collected for each of four fruit developmental stages ranging from small fruit enclosed in the calyx (stage w) to glossy full size marketable fruit (stage z) (mean fruit weight ± se of stages w, x, y, and z for ‘Classic’ and ‘Ghostbuster’, respectively: w = 2.7 ± 0.6 g, 1.9 ± 0.3 g, stage x = 7.3 ± 0.9 g, 7.7 ± 0.7 g, stage y = 50.0 ± 13.5 g, 37.7 ± 0.2 g, stage z = 645.8 ± 27.7 g, 350.8 ± 19.1 g). Each developmental stage was represented by four replicates obtained from individual fruit of each cultivar. Total RNA for individual samples was isolated from ≈1000 mg frozen tissue using the Rneasy Plant Mini Kit with the optional DNase digest (Qiagen, Valencia, CA). RNA quantity was recorded using a spectrophotometer (ND1000; NanoDrop Technologies, Wilmington, DE) and integrity and quantification confirmed using a 2100 Bioanalyzer (Agilent Technologies, Santa Clara, CA).
Real-time PCR was used to compare flavonoid gene expression (Chs, Dfr, Ans, Myb, Myc, Wd) between anthocyanin-pigmented and non-pigmented tissues collected from contrasting genotypes. Degenerate primer sets from Petunia ×hybrida Vilm. (Griesbach and Beck, 2005) were used to generate PCR products. PCR products amplified using degenerate primers were cloned and sequenced to identify gene products with homology to published flavonoid gene sequences. The sequences were used to design S. melongena gene-specific primers (Table 1). Multiple Myb products with sequence homology to flavonoid related Myb clones from potato (Solanum tuberosum L.) and tomato (Solanum lycopersicum L.) were identified and designated MybB and MybC, respectively. The sequences of amplified products were used to develop S. melongena gene-specific primers (Table 1).

For real time PCR, cDNA synthesis was performed with 1 μg RNA using the iScriptTM cDNA Synthesis Kit (Bio-Rad, Hercules, CA). Real-time PCR amplification of the cDNA was carried out as described by Lightbourn et al. (2007) using the iCycler 1QTM Multicolor Real-Time Detection System (Bio-Rad). Results were normalized to the expression of β-tubulin (Tub) where degenerate primers were designed based on gene sequences reported in GenBank (Table 1; Lightbourn et al., 2007). The efficiency of reaction kinetics and R values were tabulated and found to fall within the range of 95% to 105% (efficiency) and 0.98 to 1.0 (R value). Threshold values [Ct (iQTM 5 Version 2.0; Bio-Rad)] were manually checked for sd > 0.5 within triplicate reactions. sd of Tub values across fruit development stages and cultivars were ≤ 1.0. Relative gene expression was normalized to the reference gene as described (Lightbourn et al., 2007).
Analytical high-performance liquid chromatography (HPLC) (Griesbach et al., 1991) was used to quantify anthocyanin pigments in fruit peel tissue. Sample collection for anthocyanins was equivalent to that described for gene expression studies except duplicate samples were collected for each replicate from each of the four fruit developmental stages. Anthocyanins were extracted in acidified methanol (1% HCl), evaporated to dryness, washed with acetone, and re-solubilized in acidified methanol. HPLC characterization of anthocyanidins was performed on a 7.8 × 300-mm column of 5-μm Bondapak C18 (Waters Corp., Milford, MA) using a 30-min linear gradient of 0% to 10% (v/v) acetonitrile in aqueous 1.5% (v/v) phosphoric acid and 15% (v/v) acetic acid followed by a 10-min linear increase to 20% (v/v) acetonitrile and held at 20% (v/v) acetonitrile for an additional 20 min. Flow rate was 1.0 mL·min−1 and detection was by absorption at 540 nm. Anthocyanins were identified as previously described (Lightbourn et al., 2008). Anthocyanin content was determined by measuring sample peak areas obtained from replicate samples of each cultivar using Maxima software (Waters Corp.). The relative amount of pigment was calculated as the mean of the area of absorption per gram fresh weight.
The S. melongena cultivars Classic and Ghostbuster produce characteristically dark violet and white-colored fruit, respectively (Fig. 1). Delphinidin-3-p-coumaroylrutinoside-5-glucoside was identified as the major anthocyanin present in fruit peels of ‘Classic’ and ‘Ghostbuster’. Anthocyanin concentration increased 9-fold during early fruit development (stages w, x, y) in ‘Classic’ and subsequently declined 5-fold in full-sized market quality fruit [stage z (Table 2)]. Very low but detectable anthocyanin concentrations were found in ‘Ghostbuster’. Anthocyanin concentration in fruit peels of this white-pigmented genotype did not change significantly during fruit development and were ≈145- to 1000-fold lower relative to those found in ‘Classic’ at comparable stages of fruit development.





Citation: Journal of the American Society for Horticultural Science J. Amer. Soc. Hort. Sci. 140, 2; 10.21273/JASHS.140.2.129

The presence of measureable anthocyanin concentration in epidermal tissues of violet as well as white-colored fruit of these two genotypes indicated that functional anthocyanin biosynthetic pathways were present. The expression of regulatory genes (Myc, MybB, MybC, Wd) that comprise the MYB-MYC-WD transcription factor complex and key biosynthetic genes (Chs, Dfr, Ans) were evaluated to characterize the genetic mechanism underlying differences observed in fruit pigmentation of these cultivars.
With the exception of Chs transcript level in stage z fruit of ‘Classic’ and ‘Ghostbuster’, significant differences between ‘Classic’ and ‘Ghostbuster’ Chs, Dfr, and Ans transcript occurred at all stages of fruit development. Chs and Dfr transcript levels were ≈750- and 5700-fold greater, respectively, in ‘Classic’ vs. ‘Ghostbuster’ in stage y fruit, coincident with peak measures of anthocyanin concentration (Table 3). The greatest fold difference in structural gene expression was for Ans, where transcript levels were 353,000-fold greater in stage y fruit of ‘Classic’ in comparison with ‘Ghostbuster’. In violet-pigmented epidermal tissue of ‘Classic’ fruit, transcript levels for all three genes increased between stage w through stage y coincident with increases in anthocyanin concentration. ‘Classic’ Chs and Dfr transcript levels increased 4-fold, and Ans transcript increased 6-fold in small developing fruit between stage w and stage y. Biosynthetic transcript levels for ‘Classic’ decreased 4.5- to 14-fold between stage y and stage z, coincident with declines in anthocyanin concentration of full-sized fruit. Similar to anthocyanin concentration in developing ‘Ghostbuster’ fruit, transcript levels for all three biosynthetic genes were relatively unchanged throughout fruit development.

Myb products with sequence homology to flavonoid-related Myb clones from potato and tomato were identified in developing eggplant fruit and designated MybB and MybC, respectively. MybB transcript levels were up to 10-fold greater in developing fruit of ‘Classic’ vs. ‘Ghostbuster’ but exhibited little or no change throughout fruit development in either cultivar and were present at low levels relative to those for MybC (Table 4). Significant differences between ‘Classic’ and ‘Ghostbuster’ transcript levels were found at all stages of fruit development. MybC and Myc transcript levels were up to 14.5- and 288-fold greater in ‘Classic’ vs. ‘Ghostbuster’ fruit at comparable stages of fruit development. MybC transcript level in ‘Classic’ fruit increased 2.6-fold between stage w and stage y, and Myc transcript increased over 7.7-fold in developing fruit coincident with increasing biosynthetic transcript level and anthocyanin concentration. Transcript for both regulatory genes subsequently declined in market-sized fruit (stage z) coincident with declines in biosynthetic transcript level and anthocyanin concentration. Little or no change in MybC or Myc was observed in ‘Ghostbuster’. Wd transcript increased ≈2-fold in ‘Classic’ between stage w and stage y as biosynthetic transcript level and anthocyanin concentration increased. However, transcript levels in ‘Classic’ vs. ‘Ghostbuster’ were not significantly different at comparable stages of fruit development. Wd transcript level was relatively unchanged during fruit development in ‘Ghostbuster’ and increased in market-sized fruit of both cultivars relative to smaller stage y fruit.

Despite the importance of eggplant worldwide as a food crop and considerable interest in eggplant fruit secondary metabolites for their nutritive value, little research has been conducted to understand the molecular mechanisms that control anthocyanin accumulation in this crop. Early genetic studies describing the inheritance of eggplant fruit pigmentation (Tigchelaar et al., 1968) and more recent anthocyanin pigmentation-linked quantitative trait loci studies (Barchi et al., 2012; Cericola et al., 2014) support the action of multiple genetic factors that influence anthocyanin pigmentation. Using a genome-wide association mapping approach, Cericola et al. (2014) identified 56 single nucleotide polymorphism (SNP) markers associated with anthocyanin pigmentation. Ten of these SNP loci, scattered over five chromosomes, were associated with fruit color. Synteny between eggplant and other Solanaceous species enabled speculation on possible orthologs associated with fruit color. Eggplant fruit color-associated markers identified on chromosome 10 were syntenic with chromosomal regions in tomato and potato harboring anthocyanin-related genes, namely, UDP-glucose:flavonoid-3-O-glucosyltransferase (3GT) and two MYB transcription factors that regulate the genes encoding CHS and DFR and influence tissue specific anthocyanin expression (Cericola et al., 2014).
We report differential expression of Chs, Dfr, and Ans biosynthetic genes in eggplant fruit with violet vs. white fruit color as a result of high vs. very low anthocyanin concentration, respectively. All three biosynthetic genes were highly expressed in violet fruit relative to white fruit. Upregulation of both MybC and Myc regulatory genes, coincident with increased expression of anthocyanin biosynthetic genes and anthocyanin concentration in developing fruit post-anthesis through market maturity stages of development, suggest that differential expression of MybC and Myc are important determinants of anthocyanin accumulation in these genotypes. Whereas MybC with sequence homology to a tomato Myb clone was differentially expressed in eggplant fruit, MybB with sequence homology to a potato Myb clone did not exhibit differential expression and was expressed at very low levels relative to MybC. Zhang et al. (2014) recently reported upregulation of anthocyanin biosynthetic genes and Myb in violet anthocyanin-pigmented eggplant fruit 12 to 15 d post-anthesis. The investigators did not observe differential expression of Myc. Anthocyanin pigmentation of small developing fruit enclosed in the calyx suggests that ‘Classic’ carries the dominant Puc (Pigment under calyx) allele (Tigchelaar et al., 1968) and that anthocyanin synthesis in this cultivar is not light-dependent.
Our focus on Chs, Dfr, and Ans biosynthetic genes does not preclude differential regulation of other anthocyanin biosynthetic genes as determinants of observed fruit pigmentation. However, differential MybC and Myc expression in violet- vs. white-colored fruit suggests that regulation of anthocyanin biosynthesis and resulting fruit color in these genotypes is associated with expression of these transcription factors. In well-studied plant systems including Arabidopsis thaliana (L.) Heynh., Zea mays L., P. ×hybrida, and Vitis vinifera L. (Du et al., 2012; Dubos et al., 2010; Feller et al., 2011; Fournier-Level et al., 2010; Griesbach, 2005; Matus et al., 2008), multiple MYB encoding genes have been reported. In A. thaliana, 13 or more MYBs are related to the regulation of flavonoid metabolism and control specific anthocyanin biosynthetic enzymes and/or tissue specificity and may be redundant with other members of the MYB gene family (Dubos et al., 2010). Similarly, in grape, numerous anthocyanin MYBs were identified. However, only two MYBs located on chromosome 2 differentiated presence vs. absence of anthocyanin-pigmented fruit color (Fournier-Level et al., 2010; Matus et al., 2008). More recently, a candidate gene association study of grape accessions with a diverse range of fruit color identified fruit anthocyanin-associated transcription factors that included a MYC encoded on chromosome 2 and two MYBs, one identified on chromosome 9 and a second of unknown location (Cardoso et al., 2012). Similarly, discrete regulation of fruit anthocyanin content by one or several anthocyanin MYBs and MYCs has been reported for other Solanaceous species including pepper (Stommel et al., 2009) and tomato (Gonzali et al., 2009) and roseaceous fruits including apple (Malus domestica Borkh.) and pear (Pyrus communis L.) (Jaakola, 2013). Our results demonstrating differential regulation of MYB and MYC coincident with fruit anthocyanin pigmentation are congruent with these results. Functional studies are needed to confirm the biological roles of these eggplant transcription factors.
The role of the MYC-MYB-WD40 regulatory protein complex is well accepted for its role in anthocyanin biosynthesis (Koes et al., 2005). Numerous investigations have demonstrated that all three transcription factors are required for anthocyanin production. Myb and Myc provide one of the best examples of combinatorial plant gene regulation. Together with the WD40 protein, MYB and MYC provide extensive regulatory function in plants and result in new phenotypic traits by virtue of variation in the MYC or, more commonly, the MYB component (Ramsay and Glover, 2005). Z. mays MybC1 and MycLc are two of the alleles that are well studied. They have been expressed in a number of plant species (Bovy et al., 2002; Bradley et al., 1998; Goldsbrough et al., 1996; Lloyd et al., 1992; Quattrocchio et al., 1993) in which they either enhanced the amount of anthocyanin produced or activated biosynthesis in unpigmented tissues. We reported differential regulation of both Myb and Myc in pepper genotypes with dichotomous anthocyanin pigmentation (Stommel et al., 2009). Ma et al. (2008) reported that Myb as well as Myc expression were required to elicit high levels of anthocyanin accumulation in Phalaenopsis amabilis (L.) Blume flowers. In P. ×hybrida, the combination MycAn1/MybAn2 induces anthocyanin pigmentation in the flower; MycAn1/MybAn4 induces anthocyanin pigmentation in the anthers; whereas the MycAn1/MybPh4 combination induces vacuolar acidification (Quattrocchio et al., 2006). Neither regulatory gene alone was sufficient to influence anthocyanin content. Differences in the threshold requirements for these regulatory genes to activate anthocyanin production also differ in different plants. In mammalian cells, both the MYC protein and mRNA have relatively short half-lives (Sharova et al., 2009). As a result, its temporal level and functionality in transcription factor complexes is dependent on transcription rate or promoter strength.
Numerous studies have documented the role of plant MYB proteins in conferring tissue- and cell-specific anthocyanin regulation and their influence on regulatory control of multiple aspects of epidermal cell differentiation (Koes et al., 2005; Quattrocchio et al., 2006). Relative to MYB proteins, plant MYC proteins have been studied less extensively. This class of transcription factors was first identified in Z. mays where they determine the timing, amount, and distribution of anthocyanin pigmentation (Ludwig et al., 1989). The N-terminal region of plant Myc bHLH transcription factors has an ≈190 amino acid N-terminus that defines the MYB-interacting region. After this interaction domain is an activation domain that is responsible for transactivation of biosynthetic gene promoters in the anthocyanin pathway. Pattanaik et al. (2006) demonstrated that a lysine to methionine substitution in the interaction domain of MYC-RP of Perilla frutescens (L.) Britt. resulted in a 50-fold increase in transactivation activity. Ectopic expression of MYC-RP in transgenic tobacco (Nicotiana tabacum L.) resulted in increased anthocyanin accumulation relative to expression of the wild-type gene. Pattanaik et al. (2006) also documented MYC variants from Antirrhinum majus L. that showed increased transactivation of the A. thaliana Chs promoter. Changes in the interaction domain likely influence anthocyanin biosynthesis through 1) increased transactivation strength of respective MYC variants; 2) increased transactivation strength and improved interaction of MYC variants with MYB proteins; or 3) more stable bHLH-MYB-WD complexes resulting in enhanced anthocyanin production (Pattanaik et al., 2008).
In the current study, Wd expression was constitutive in both violet anthocyanin-pigmented eggplant fruit and white fruit with very low anthocyanin concentration. Wd transcript levels were comparable in both genotypes and were not differentially regulated coincident with anthocyanin accumulation. Similarly, differential Wd expression was not observed in P. amabilis (Ma et al., 2009) or pepper (Stommel et al., 2009). Despite lack of differential regulation, the essential role of WD40 proteins in the transcription factor complex is highlighted in recent studies (Aguilar-Barragan and Ochoa-Alejo, 2014) where silencing of Wd resulted in significant reduction of anthocyanin biosynthetic genes. All of these observations are consistent with the proposed model for the MYC-MYB-WD40 regulatory protein complex where WD40 serves as a platform for interaction with MYCs that in turn interact with MYBs (Koes et al., 2005). MYB and MYC appear to be sufficient in providing specificity in the regulatory complex, which binds to structural gene promoters to initiate gene expression.
Available studies document diversity in anthocyanin-related MYB, MYC, and WD40 transcription factors. DNA as well as protein binding sites of MYB, MYC, and WD40 proteins impose differential constraints on variability that can be tolerated in these proteins so that structure and functionality can be maintained (Ramsay and Glover, 2005; Streisfeld et al., 2011). Nonetheless, sufficient variation exists among one or multiple components of the transcription factor complex to produce varied anthocyanin composition and concentration and tissue and temporal specificity observed for anthocyanin biosynthesis. Our results on coordinated regulation of anthocyanin biosynthetic and regulatory genes in S. melongena further illustrate diversity within and among these transcription factors in Solanaceous species and more distant genera.

Contributor Notes
Mention of trade names or commercial products is solely for the purpose of providing specific information and does not imply recommendation or endorsement by the U.S. Department of Agriculture.
Corresponding author. E-mail: john.stommel@ars.usda.gov.