Dissecting genome-wide association signals for loss-of-function phenotypes in sorghum flavonoid pigmentation traits

Identification and analysis of novel recessive alleles for Tan1 and Tan2 in sorghum

Background

The identification and analysis of allelic variation are important bases for crop diversity research, trait domestication and molecular marker development. Grain tannin content is a very important quality trait in sorghum. Higher tannin levels in sorghum grains are usually required when breeding varieties resistant to bird damage or those used for brewing liquor. Non-tannin-producing or low-tannin-producing sorghum accessions are commonly used for food and forage. Tan1 and Tan2, two important cloned genes, regulate tannin biosynthesis in sorghum, and mutations in one or two genes will result in low or no tannin content in sorghum grains. Even if sorghum accessions contain dominant Tan1 and Tan2, the tannin contents are distributed from low to high, and there must be other new alleles of the known regulatory genes or new unknown genes contributing to tannin production.

Methods

The two parents 8R306 and 8R191 did not have any known recessive alleles for Tan1 and Tan2, and it was speculated that they probably both had dominant Tan1 and Tan2 genotypes. However, the phenotypes of two parents were different; 8R306 had tannins and 8R191 had non-tannins in the grains, so these two parents were constructed as a RIL population. Bulked segregant analysis (BSA) was used to determine other new alleles of Tan1 and Tan2 or new Tannin locus. Tan1 and Tan2 full-length sequences and tannin contents were detected in wild sorghum resources, landraces and cultivars.

Results

We identified two novel recessive tan1-d and tan1-e alleles and four recessive Tan2 alleles, named as tan2-d, tan2-e, tan2-f, and tan2-g. These recessive alleles led to loss of function of Tan1 and Tan2, and low or no tannin content in sorghum grains. The loss-of-function alleles of tan1-e and tan2-e were only found in Chinese landraces, and other alleles were found in landraces and cultivars grown all around the world. tan1-a and tan1-b were detected in foreign landraces, Chinese cultivars and foreign cultivars, but not in Chinese landraces.

Conclusion

These results implied that Tan1 and Tan2 recessive alleles had different geographically distribution in the worldwide, but not all recessive alleles had been used in breeding. The discovery of these new alleles provided new germplasm resources for breeding sorghum cultivars for food and feed, and for developing molecular markers for low-tannin or non-tannin cultivar-assisted breeding in sorghum.

Major impacts of widespread structural variation on sorghum

Genetic diversity is critical to crop breeding and improvement, and dissection of the genomic variation underlying agronomic traits can both assist breeding and give insight into basic biological mechanisms. Although recent genome analyses in plants reveal many structural variants (SVs), most current studies of crop genetic variation are dominated by single-nucleotide polymorphisms (SNPs). The extent of the impact of SVs on global trait variation, as well as their utility in genome-wide selection, is not yet understood. In this study, we built an SV data set based on whole-genome resequencing of diverse sorghum lines (n = 363), validated the correlation of photoperiod sensitivity and variety type, and identified SV hotspots underlying the divergent evolution of cellulosic and sweet sorghum. In addition, we showed the complementary contribution of SVs for heritability of traits related to sorghum adaptation. Importantly, inclusion of SV polymorphisms in association studies revealed genotype-phenotype associations not observed with SNPs alone. Three-way genome-wide association studies (GWAS) based on whole-genome SNP, SV, and integrated SNP + SV data sets showed substantial associations between SVs and sorghum traits. The addition of SVs to GWAS substantially increased heritability estimates for some traits, indicating their important contribution to functional allelic variation at the genome level. Our discovery of the widespread impacts of SVs on heritable gene expression variation could render a plausible mechanism for their disproportionate impact on phenotypic variation. This study expands our knowledge of SVs and emphasizes the extensive impacts of SVs on sorghum.

Novel allelic variations in Tannin1 and Tannin2 contribute to tannin absence in sorghum

Sorghum is an important food crop commonly used for brewing, feed, and bioenergy. Certain genotypes of sorghum contain high concentrations of condensed tannins in seeds, which are beneficial, such as protecting grains from herbivore bird pests, but also impair grain quality and digestibility. Previously, we identified Tannin1 and Tannin2, each with three recessive causal alleles, regulate tannin absence in sorghum. In this study, via characterizing 421 sorghum accessions, we further identified three novel recessive alleles from these two genes. The tan1-d allele contains a 12-bp deletion at position 659 nt and the tan1-e allele contains a 10-bp deletion at position 771 nt in Tannin1. The tan2-d allele contains a C-to-T transition, which results in a premature stop codon before the bHLH domain in Tannin2, and was predominantly selected in China. We further developed KASP assays targeting these identified recessive alleles to efficiently genotype large populations. These studies provide new insights in sorghum domestication and convenient tools for breeding programs.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-024-01463-y.

Deleterious mutations predicted in the sorghum (Sorghum bicolor) Maturity (Ma) and Dwarf (Dw) genes from whole-genome resequencing

In sorghum [Sorghum bicolor (L.) Moench] the Maturity (Ma1, Ma2, Ma3, Ma4, Ma5, Ma6) and Dwarf (Dw1, Dw2, Dw3, Dw4) loci, encode genes controlling flowering time and plant height, respectively, which are critical for designing sorghum ideotypes for a maturity timeframe and a harvest method. Publicly available whole-genome resequencing data from 860 sorghum accessions was analyzed in silico to identify genomic variants at 8 of these loci (Ma1, Ma2, Ma3, Ma5, Ma6, Dw1, Dw2, Dw3) to identify novel loss of function alleles and previously characterized ones in sorghum germplasm. From ~ 33 million SNPs and ~ 4.4 million InDels, 1445 gene variants were identified within these 8 genes then evaluated for predicted effect on the corresponding encoded proteins, which included newly identified mutations (4 nonsense, 15 frameshift, 28 missense). Likewise, most accessions analyzed contained predicted loss of function alleles (425 ma1, 22 ma2, 40 ma3, 74 ma5, 414 ma6, 289 dw1, 268 dw2 and 45 dw3) at multiple loci, but 146 and 463 accessions had no predicted ma or dw mutant alleles, respectively. The ma and dw alleles within these sorghum accessions represent a valuable source for manipulating flowering time and plant height to develop the full range of sorghum types: grain, sweet and forage/biomass.

Recent advancements in the breeding of sorghum crop: current status and future strategies for marker-assisted breeding

Sorghum is emerging as a model crop for functional genetics and genomics of tropical grasses with abundant uses, including food, feed, and fuel, among others. It is currently the fifth most significant primary cereal crop. Crops are subjected to various biotic and abiotic stresses, which negatively impact on agricultural production. Developing high-yielding, disease-resistant, and climate-resilient cultivars can be achieved through marker-assisted breeding. Such selection has considerably reduced the time to market new crop varieties adapted to challenging conditions. In the recent years, extensive knowledge was gained about genetic markers. We are providing an overview of current advances in sorghum breeding initiatives, with a special focus on early breeders who may not be familiar with DNA markers. Advancements in molecular plant breeding, genetics, genomics selection, and genome editing have contributed to a thorough understanding of DNA markers, provided various proofs of the genetic variety accessible in crop plants, and have substantially enhanced plant breeding technologies. Marker-assisted selection has accelerated and precised the plant breeding process, empowering plant breeders all around the world.

Genetic analysis of seed traits in Sorghum bicolor that affect the human gut microbiome

Prebiotic fibers, polyphenols and other molecular components of food crops significantly affect the composition and function of the human gut microbiome and human health. The abundance of these, frequently uncharacterized, microbiome-active components vary within individual crop species. Here, we employ high throughput in vitro fermentations of pre-digested grain using a human microbiome to identify segregating genetic loci in a food crop, sorghum, that alter the composition and function of human gut microbes. Evaluating grain produced by 294 sorghum recombinant inbreds identifies 10 loci in the sorghum genome associated with variation in the abundance of microbial taxa and/or microbial metabolites. Two loci co-localize with sorghum genes regulating the biosynthesis of condensed tannins. We validate that condensed tannins stimulate the growth of microbes associated with these two loci. Our work illustrates the potential for genetic analysis to systematically discover and characterize molecular components of food crops that influence the human gut microbiome.

Diet affects the human gut microbiome, but studies linking crop genetics to seed traits that influence the human gut microbiome are lacking. Here, the authors develop an in vitro microbiome screening method and reveal the association between sorghum genes regulating condensed tannin biosynthesis and human gut microbiome.

Sorghum genetic, genomic, and breeding resources

Sorghum research has entered an exciting and fruitful era due to the genetic, genomic, and breeding resources that are now available to researchers and plant breeders. As the world faces the challenges of a rising population and a changing global climate, new agricultural solutions will need to be developed to address the food and fiber needs of the future. To that end, sorghum will be an invaluable crop species as it is a stress-resistant C₄ plant that is well adapted for semi-arid and arid regions. Sorghum has already remained as a staple food crop in many parts of Africa and Asia and is critically important for animal feed and niche culinary applications in other regions, such as the United States. In addition, sorghum has begun to be developed into a promising feedstock for forage and bioenergy production. Due to this increasing demand for sorghum and its potential to address these needs, the continuous development of powerful community resources is required. These resources include vast collections of sorghum germplasm, high-quality reference genome sequences, sorghum association panels for genome-wide association studies of traits involved in food and bioenergy production, mutant populations for rapid discovery of causative genes for phenotypes relevant to sorghum improvement, gene expression atlas, and online databases that integrate all resources and provide the sorghum community with tools that can be used in breeding and genomic studies. Used in tandem, these valuable resources will ensure that the rate, quality, and collaborative potential of ongoing sorghum improvement efforts is able to rival that of other major crops.

SorGSD: updating and expanding the sorghum genome science database with new contents and tools

BACKGROUND

As the fifth major cereal crop originated from Africa, sorghum (Sorghum bicolor) has become a key C4 model organism for energy plant research. With the development of high-throughput detection technologies for various omics data, much multi-dimensional and multi-omics information has been accumulated for sorghum. Integrating this information may accelerate genetic research and improve molecular breeding for sorghum agronomic traits.

RESULTS

We updated the Sorghum Genome SNP Database (SorGSD) by adding new data, new features and renamed it to Sorghum Genome Science Database (SorGSD). In comparison with the original version SorGSD, which contains SNPs from 48 sorghum accessions mapped to the reference genome BTx623 (v2.1), the new version was expanded to 289 sorghum lines with both single nucleotide polymorphisms (SNPs) and small insertions/deletions (INDELs), which were aligned to the newly assembled and annotated sorghum genome BTx623 (v3.1). Moreover, phenotypic data and panicle pictures of critical accessions were provided in the new version. We implemented new tools including ID Conversion, Homologue Search and Genome Browser for analysis and updated the general information related to sorghum research, such as online sorghum resources and literature references. In addition, we deployed a new database infrastructure and redesigned a new user interface as one of the Genome Variation Map databases. The new version SorGSD is freely accessible online at http://ngdc.cncb.ac.cn/sorgsd/ .

CONCLUSIONS

SorGSD is a comprehensive integration with large-scale genomic variation, phenotypic information and incorporates online data analysis tools for data mining, genome navigation and analysis. We hope that SorGSD could provide a valuable resource for sorghum researchers to find variations they are interested in and generate customized high-throughput datasets for further analysis.

Genetic characterization of a Sorghum bicolor multiparent mapping population emphasizing carbon-partitioning dynamics

Sorghum bicolor, a photosynthetically efficient C₄ grass, represents an important source of grain, forage, fermentable sugars, and cellulosic fibers that can be utilized in myriad applications ranging from bioenergy to bioindustrial feedstocks. Sorghum’s efficient fixation of carbon per unit time per unit area per unit input has led to its classification as a preferred biomass crop highlighted by its designation as an advanced biofuel by the U.S. Department of Energy. Due to its extensive genetic diversity and worldwide colonization, sorghum has considerable diversity for a range of phenotypes influencing productivity, composition, and sink/source dynamics. To dissect the genetic basis of these key traits, we present a sorghum carbon-partitioning nested association mapping (NAM) population generated by crossing 11 diverse founder lines with Grassl as the single recurrent female. By exploiting existing variation among cellulosic, forage, sweet, and grain sorghum carbon partitioning regimes, the sorghum carbon-partitioning NAM population will allow the identification of important biomass-associated traits, elucidate the genetic architecture underlying carbon partitioning and improve our understanding of the genetic determinants affecting unique phenotypes within Poaceae. We contrast this NAM population with an existing grain population generated using Tx430 as the recurrent female. Genotypic data are assessed for quality by examining variant density, nucleotide diversity, linkage decay, and are validated using pericarp and testa phenotypes to map known genes affecting these phenotypes. We release the 11-family NAM population along with corresponding genomic data for use in genetic, genomic, and agronomic studies with a focus on carbon-partitioning regimes.

Meta-Analysis Identifies Pleiotropic Loci Controlling Phenotypic Trade-offs in Sorghum

Community association populations are composed of phenotypically and genetically diverse accessions. Once these populations are genotyped, the resulting marker data can be reused by different groups investigating the genetic basis of different traits. Because the same genotypes are observed and scored for a wide range of traits in different environments, these populations represent a unique resource to investigate pleiotropy. Here we assembled a set of 234 separate trait datasets for the Sorghum Association Panel, a group of 406 sorghum genotypes widely employed by the sorghum genetics community. Comparison of genome wide association studies conducted with two independently generated marker sets for this population demonstrate that existing genetic marker sets do not saturate the genome and likely capture only 35-43% of potentially detectable loci controlling variation for traits scored in this population. While limited evidence for pleiotropy was apparent in cross-GWAS comparisons, a multivariate adaptive shrinkage approach recovered both known pleiotropic effects of existing loci and new pleiotropic effects, particularly significant impacts of known dwarfing genes on root architecture. In addition, we identified new loci with pleiotropic effects consistent with known trade-offs in sorghum development. These results demonstrate the potential for mining existing trait datasets from widely used community association populations to enable new discoveries from existing trait datasets as new, denser genetic marker datasets are generated for existing community association populations.

Genome-wide association analysis reveals seed protein loci as determinants of variations in grain mold resistance in sorghum

GWAS analysis revealed variations at loci harboring seed storage, late embryogenesis abundant protein, and a tannin biosynthesis gene associated with sorghum grain mold resistance. Grain mold is the most important disease of sorghum [Sorghum bicolor (L.) Moench]. It starts at the early stages of grain development due to concurrent infection by multiple fungal species. The genetic architecture of resistance to grain mold is poorly understood. Using a diverse set of 635 Ethiopian sorghum accessions, we conducted a multi-stage disease rating for resistance to grain mold under natural infestation in the field. Through genome-wide association analyses with 173,666 SNPs and multiple models, two novel loci were identified that were consistently associated with grain mold resistance across environments. Sequence variation at new loci containing sorghum KAFIRIN gene encoding a seed storage protein affecting seed texture and LATE EMBRYOGENESIS ABUNDANT 3 (LEA3) gene encoding a protein that accumulates in seeds, previously implicated in stress tolerance, were significantly associated with grain mold resistance. The KAFIRIN and LEA3 loci were also significant factors in grain mold resistance in accessions with non-pigmented grains. Moreover, we consistently detected the known SNP (S4_62316425) in TAN1 gene, a regulator of tannin accumulation in sorghum grain to be significantly associated with grain mold resistance. Identification of loci associated with new mechanisms of resistance provides fresh insight into genetic control of the trait, while the highly resistant accessions can serve as sources of resistance genes for breeding. Overall, our association data suggest the critical role of loci harboring seed protein genes and implicate grain chemical and physical properties in sorghum grain mold resistance.

Genome-wide variations analysis of sorghum cultivar Hongyingzi for brewing Moutai liquor

Background

Hongyingzi is a sorghum (Sorghum bicolor L. Moench) cultivar for brewing Moutai liquor. For an overall understanding of the whole genome of Hongyingzi, we performed whole-genome resequencing technology to reveal its comprehensive variations.

Results

Compared with the BTx623 reference genome, we uncovered 1,885,774 single nucleotide polymorphisms (SNPs), 309,381 small fragments insertions and deletions (Indels), 31,966 structural variations (SVs), and 217,273 copy number variations (CNVs). These alterations conferred 29,614 gene variations. It was also predicted that 35 gene variations were related to the multidrug and toxic efflux (MATE) transporter, chalcone synthase (CHS), ATPase isoform 10 (AHA10) transporter, dihydroflavonol-4-reductase (DFR), the laccase 15 (LAC15), flavonol 3′-hydroxylase (F3′H), flavanone 3-hydroxylase (F3H), O-methyltransferase (OMT), flavonoid 3′5′ hydroxylase (F3′5′H), UDP-glucose:sterol-glucosyltransferase (SGT), flavonol synthase (FLS), and chalcone isomerase (CHI) involved in the tannin synthesis.

Conclusions

These results would provide theoretical supports for the molecular markers developments and gene function studies related to the tannin synthesis, and the genetic improvement of liquor-making sorghum based on the genome editing technology.

Genomics-assisted breeding in minor and pseudo-cereals

Minor and pseudo-cereals, which can grow with lower input and often produce specific nutrients compared to major cereal crops, are attracting worldwide attention. Since these crops generally have a large genetic diversity in a breeding population, rapid genetic improvement can be possible by the application of genomics-assisted breeding methods. In this review, we discuss studies related to biparental quantitative trait locus (QTL) mapping, genome-wide association study, and genomic selection for minor and pseudo-cereals. Especially, we focus on the current progress in a pseudo-cereal, buckwheat. Prospects for the practical utilization of genomics-assisted breeding in minor and pseudo-cereals are discussed including the issues to overcome especially for these crops.

Genetic Architecture of Chilling Tolerance in Sorghum Dissected with a Nested Association Mapping Population

Dissecting the genetic architecture of stress tolerance in crops is critical to understand and improve adaptation. In temperate climates, early planting of chilling-tolerant varieties could provide longer growing seasons and drought escape, but chilling tolerance (<15°) is generally lacking in tropical-origin crops. Here we developed a nested association mapping (NAM) population to dissect the genetic architecture of early-season chilling tolerance in the tropical-origin cereal sorghum (Sorghum bicolor [L.] Moench). The NAM resource, developed from reference line BTx623 and three chilling-tolerant Chinese lines, is comprised of 771 recombinant inbred lines genotyped by sequencing at 43,320 single nucleotide polymorphisms. We phenotyped the NAM population for emergence, seedling vigor, and agronomic traits (>75,000 data points from ∼16,000 plots) in multi-environment field trials in Kansas under natural chilling stress (sown 30-45 days early) and normal growing conditions. Joint linkage mapping with early-planted field phenotypes revealed an oligogenic architecture, with 5-10 chilling tolerance loci explaining 20-41% of variation. Surprisingly, several of the major chilling tolerance loci co-localize precisely with the classical grain tannin (Tan1 and Tan2) and dwarfing genes (Dw1 and Dw3) that were under strong directional selection in the US during the 20^th century. These findings suggest that chilling sensitivity was inadvertently selected due to coinheritance with desired nontannin and dwarfing alleles. The characterization of genetic architecture with NAM reveals why past chilling tolerance breeding was stymied and provides a path for genomics-enabled breeding of chilling tolerance.

Allelochemicals targeted to balance competing selections in African agroecosystems

Among major cereals domesticated as staple food, only sorghum has a high proportion of cultivars with condensed tannins in grain, which can trigger bitter taste perception in animals by binding to type 2 taste receptors (TAS2Rs). Here, we report the completion of uncovering of a pair of duplicate recessive genes (Tannin1 and Tannin2) underlying tannin presence. Three loss-of-function alleles from each gene were identified in non-tannin sorghum desired as palatable food. Condensed tannins effectively prevented sparrows from consuming sorghum grain. Parallel geographic distributions between tannin sorghum and Quelea quelea supported the role of tannins in fighting against this major herbivore threat. Association between geographic distributions of human TAS2R variants and tannin sorghum across Africa suggested that different causes had probably driven this bidirectional selection according to varied local herbivore threats and human taste sensitivity. Our investigation uncovered coevolution among humans, plants and environments linked by allelochemicals.

Genome-Wide Association Mapping of Anthracnose (Colletotrichum sublineolum) Resistance in NPGS Ethiopian Sorghum Germplasm

The National Plant Germplasm System (NPGS) Ethiopian sorghum [Sorghum bicolor (L.) Moench] collection of the United States is an important genetic resource for sorghum improvement. Anthracnose (Colletotrichum sublineolum) is one of the most harmful fungal diseases in humid sorghum production regions. Although multiple resistance sources have been identified in temperate-adapted germplasm in the Sorghum Association Panel (SAP), these resistance loci explain a limited portion of the total variation, and sources of resistance from tropical germplasm are not available for breeding programs at temperate regions. Using a core set of 335 previously genotyped NPGS Ethiopian accessions, we identified 169 accessions resistant to anthracnose. To identify resistance loci, we merged the genotypic and anthracnose response data for both NPGS Ethiopian germplasm and the SAP and performed genome-wide association scans using 219,037 single nucleotide polymorphisms and 617 accessions. The integrated data set enabled the detection of a locus on chromosome 9 present in the SAP at a low frequency. The locus explains a limited portion of the observed phenotypic variation (r² = 0.31), suggesting the presence of other resistance loci. The locus in chromosome 9 was constituted by three R genes clustered within a 47-kb region. The presence of multiple sources of resistance in NPGS Ethiopian germplasm and SAP requires the inclusion of other resistance response evaluation that could revealed others low frequency resistance alleles in the panel.

Genome-Wide Association Mapping of Grain Mold Resistance in the US Sorghum Association Panel

Sorghum [ (L.) Moench] production in warm and humid regions is limited by grain mold disease, which can be caused by a complex of >40 pathogenic and opportunistic fungi. The identification of resistant plants within temperate-adapted germplasm is imperative for the development of better-adapted varieties. The performance of 331 accessions from the previously genotyped sorghum association panel (SAP) was evaluated in four tropical environments. Only 18 accessions showed low seed deterioration and high emergence rates. The resistant accessions showed high variation in seed tannin contents and panicle shape, indicating that grain mold resistance is not associated with a single phenotypic trait. Seed mycoflora analysis recovered pathogenic fungi , , and in both resistant and susceptible accessions. By genome-wide association scans using 268,289 single nucleotide polymorphisms (SNPs), we identified two loci associated with low seed deterioration and another associated with emergence rate. Candidate genes within these loci included one gene () and two genes ( and ) with domains associated with systemic acquired resistance, suggesting that resistance involved pathogen recognition and downstream signaling cascades. This study provides insight into the genetic control of grain mold resistance as well as valuable accessions for breeding programs in temperate environments.

GWAS analysis of sorghum association panel lines identifies SNPs associated with disease response to Texas isolates of Colletotrichum sublineola

SNPs identify prospective genes related to response to Colletotrichum sublineola (anthracnose) in the sorghum association panel lines. Sorghum association panel (SAP) lines were scored over several years for response to Colletotrichum sublineola, the causal agent of the disease anthracnose. Known resistant and susceptible lines were included each year to verify successful inoculation. Over 79,000 single-nucleotide polymorphic (SNP) loci from a publicly available genotype by sequencing dataset available for the SAP lines were used with TASSEL association mapping software to identify chromosomal locations associated with differences in disease response. When the top-scoring SNPs were mapped to the published sorghum genome, in each case, the nearest annotated gene has precedence for a role in host defense.

An Integrated Genotyping-by-Sequencing Polymorphism Map for Over 10,000 Sorghum Genotypes

Mining crop genomic variation can facilitate the genetic research of complex traits and molecular breeding. In sorghum [ L. (Moench)], several large-scale single nucleotide polymorphism (SNP) datasets have been generated using genotyping-by-sequencing of KI reduced representation libraries. However, data reuse has been impeded by differences in reference genome coordinates among datasets. To facilitate reuse of these data, we constructed and characterized an integrated 459,304-SNP dataset for 10,323 sorghum genotypes on the version 3.1 reference genome. The SNP distribution showed high enrichment in subtelomeric chromosome arms and in genic regions (48% of SNPs) and was highly correlated ( = 0.82) to the distribution of KI restriction sites. The genetic structure reflected population differences by botanical race, as well as familial structure among recombinant inbred lines (RILs). Faster linkage disequilibrium decay was observed in the diversity panel than in the RILs, as expected, given the greater opportunity for recombination in diverse populations. To validate the quality and utility of the integrated SNP dataset, we used genome-wide association studies (GWAS) of genebank phenotype data, precisely mapping several known genes (e.g and ) and identifying novel associations for other traits. We further validated the dataset with GWAS of new and published plant height and flowering time data in a nested association mapping population, precisely mapping known genes and identifying epistatic interactions underlying both traits. These findings validate this integrated SNP dataset as a useful genomics resource for sorghum genetics and breeding.

Genetic and genomic resources of sorghum to connect genotype with phenotype in contrasting environments

With the recent development of genomic resources and high-throughput phenotyping platforms, the 21st century is primed for major breakthroughs in the discovery, understanding and utilization of plant genetic variation. Significant advances in agriculture remain at the forefront to increase crop production and quality to satisfy the global food demand in a changing climate all while reducing the environmental impacts of the world's food production. Sorghum, a resilient C₄ grain and grass important for food and energy production, is being extensively dissected genetically and phenomically to help connect the relationship between genetic and phenotypic variation. Unlike genetically modified crops such as corn or soybean, sorghum improvement has relied heavily on public research; thus, many of the genetic resources serve a dual purpose for both academic and commercial pursuits. Genetic and genomic resources not only provide the foundation to identify and understand the genes underlying variation, but also serve as novel sources of genetic and phenotypic diversity in plant breeding programs. To better disseminate the collective information of this community, we discuss: (i) the genomic resources of sorghum that are at the disposal of the research community; (ii) the suite of sorghum traits as potential targets for increasing productivity in contrasting environments; and (iii) the prospective approaches and technologies that will help to dissect the genotype-phenotype relationship as well as those that will apply foundational knowledge for sorghum improvement.

A Large-Scale Genome-Wide Association Analyses of Ethiopian Sorghum Landrace Collection Reveal Loci Associated With Important Traits

The eastern Africa region, Ethiopia and its surroundings, is considered as the center of origin and diversity for sorghum, and has contributed to global sorghum genetic improvement. The germplasm from this region harbors enormous genetic variation for various traits but little is known regarding the genetic architecture of most traits. Here, 1425 Ethiopian landrace accessions were phenotyped under field conditions for presence or absence of awns, panicle compactness and shape, panicle exsertion, pericarp color, glume cover, plant height and smut resistance under diverse environmental conditions in Ethiopia. In addition, F1 hybrids obtained from a subset of 1341 accessions crossed to an A1 cytoplasmic male sterile line, ATx623, were scored for fertility/sterility reactions. Subsequently, genotyping-by-sequencing generated a total of 879,407 SNPs from which 72,190 robust SNP markers were selected after stringent quality control (QC). Pairwise distance-based hierarchical clustering identified 11 distinct groups. Of the genotypes assigned to either one of the 11 sub-populations, 65% had high ancestry membership coefficient with the likelihood of more than 0.60 and the remaining 35% represented highly admixed accessions. A genome-wide association study (GWAS) identified loci and SNPs associated with aforementioned traits. GWAS based on compressed mixed linear model (CMLM) identified SNPs with significant association (FDR ≤ 0.05) to the different traits studied. The percentage of total phenotypic variation explained with significant SNPs across traits ranged from 2 to 43%. Candidate genes showing significant association with different traits were identified. The sorghum bHLH transcription factor, ABORTED MICROSPORES was identified as a strong candidate gene conditioning male fertility. Notably, sorghum CLAVATA1 receptor like kinase, known for regulation of plant growth, and the ETHYLENE RESPONSIVE TRANSCRIPTION FACTOR gene RAP2-7, known to suppress transition to flowering, were significantly associated with plant height. In addition, the YELLOW SEED1 like MYB transcription factor and TANNIN1 showed strong association with pericarp color validating previous observations. Overall, the genetic architecture of natural variation representing the complex Ethiopian sorghum germplasm was established. The study contributes to the characterization of genes and alleles controlling agronomic traits, and will serve as a source of markers for molecular breeding.

Genome-Wide Association Mapping of Anthracnose (Colletotrichum sublineolum) Resistance in the U.S. Sorghum Association Panel

The productivity and profitability of sorghum [ (L.) Moench] is reduced by susceptibility to fungal diseases, such as anthracnose ( P. Henn.). A limited number of resistant accessions are present in the temperate-adapted germplasm; other exotic sources of resistance are not currently available for breeding programs. Among 335 accessions available to breeders from a previously genotyped sorghum association panel (SAP), we found that 75 were resistant to anthracnose. A phylogenetic analysis of these accessions showed high genetic diversity and multiple resistance sources. Genome-wide association scans (GWAS) were conducted using 268,289 single-nucleotide polymorphisms to identify loci associated with anthracnose resistance. Using logistic regressions for binary measures of resistance responses, we identified three loci within a region on chromosome 5 that have been previously associated with three sources of anthracnose resistance. A GWAS limited to Caudatum germplasm identified an association with a region on chromosome 1 and with the same previous region on chromosome 5. Candidate genes within these loci were related to R-gene families, signaling cascades, and transcriptional reprogramming, suggesting that the resistance response is controlled by multiple defense mechanisms. The strategic integration of exotic resistant germplasm into the SAP is needed to identify additional rare resistance alleles via GWAS.

Population genomics of sorghum (Sorghum bicolor) across diverse agroclimatic zones of Niger

Improving adaptation of staple crops in developing countries is important to ensure food security. In the West African country of Niger, the staple crop sorghum (Sorghum bicolor) is cultivated across diverse agroclimatic zones, but the genetic basis of local adaptation has not been described. The objectives of this study were to characterize the genomic diversity of sorghum from Niger and to identify genomic regions conferring local adaptation to agroclimatic zones and farmer preferences. We analyzed 516 Nigerien accessions for which local variety name, botanical race, and geographic origin were known. We discovered 144 299 single nucleotide polymorphisms (SNPs) using genotyping-by-sequencing (GBS). We performed discriminant analysis of principal components (DAPC), which identified six genetic groups, and performed a genome scan for loci with high discriminant loadings. The highest discriminant coefficients were on chromosome 9, near the putative ortholog of maize flowering time adaptation gene Vgt1. Next, we characterized differentiation among local varieties and used a genome scan of pairwise F_ST values to identify SNPs associated with specific local varieties. Comparison of varieties named for light- versus dark-grain identified differentiation near Tannin1, the major gene responsible for grain tannins. These findings could facilitate genomics-assisted breeding of locally adapted and farmer-preferred sorghum varieties for Niger.

Genomic Prediction and Association Mapping of Curd-Related Traits in Gene Bank Accessions of Cauliflower

Genetic resources are an important source of genetic variation for plant breeding. Genome-wide association studies (GWAS) and genomic prediction greatly facilitate the analysis and utilization of useful genetic diversity for improving complex phenotypic traits in crop plants. We explored the potential of GWAS and genomic prediction for improving curd-related traits in cauliflower (Brassica oleracea var. botrytis) by combining 174 randomly selected cauliflower gene bank accessions from two different gene banks. The collection was genotyped with genotyping-by-sequencing (GBS) and phenotyped for six curd-related traits at two locations and three growing seasons. A GWAS analysis based on 120,693 single-nucleotide polymorphisms identified a total of 24 significant associations for curd-related traits. The potential for genomic prediction was assessed with a genomic best linear unbiased prediction model and BayesB. Prediction abilities ranged from 0.10 to 0.66 for different traits and did not differ between prediction methods. Imputation of missing genotypes only slightly improved prediction ability. Our results demonstrate that GWAS and genomic prediction in combination with GBS and phenotyping of highly heritable traits can be used to identify useful quantitative trait loci and genotypes among genetically diverse gene bank material for subsequent utilization as genetic resources in cauliflower breeding.

A High-Throughput, Field-Based Phenotyping Technology for Tall Biomass Crops

Recent advances in omics technologies have not been accompanied by equally efficient, cost-effective, and accurate phenotyping methods required to dissect the genetic architecture of complex traits. Even though high-throughput phenotyping platforms have been developed for controlled environments, field-based aerial and ground technologies have only been designed and deployed for short-stature crops. Therefore, we developed and tested Phenobot 1.0, an auto-steered and self-propelled field-based high-throughput phenotyping platform for tall dense canopy crops, such as sorghum (Sorghum bicolor). Phenobot 1.0 was equipped with laterally positioned and vertically stacked stereo RGB cameras. Images collected from 307 diverse sorghum lines were reconstructed in 3D for feature extraction. User interfaces were developed, and multiple algorithms were evaluated for their accuracy in estimating plant height and stem diameter. Tested feature extraction methods included the following: (1) User-interactive Individual Plant Height Extraction (UsIn-PHe) based on dense stereo three-dimensional reconstruction; (2) Automatic Hedge-based Plant Height Extraction (Auto-PHe) based on dense stereo 3D reconstruction; (3) User-interactive Dense Stereo Matching Stem Diameter Extraction; and (4) User-interactive Image Patch Stereo Matching Stem Diameter Extraction (IPaS-Di). Comparative genome-wide association analysis and ground-truth validation demonstrated that both UsIn-PHe and Auto-PHe were accurate methods to estimate plant height, while Auto-PHe had the additional advantage of being a completely automated process. For stem diameter, IPaS-Di generated the most accurate estimates of this biomass-related architectural trait. In summary, our technology was proven robust to obtain ground-based high-throughput plant architecture parameters of sorghum, a tall and densely planted crop species.

Increased Power To Dissect Adaptive Traits in Global Sorghum Diversity Using a Nested Association Mapping Population

Adaptation of domesticated species to diverse agroclimatic regions has led to abundant trait diversity. However, the resulting population structure and genetic heterogeneity confounds association mapping of adaptive traits. To address this challenge in sorghum [Sorghum bicolor (L.) Moench]-a widely adapted cereal crop-we developed a nested association mapping (NAM) population using 10 diverse global lines crossed with an elite reference line RTx430. We characterized the population of 2214 recombinant inbred lines at 90,000 SNPs using genotyping-by-sequencing. The population captures ∼70% of known global SNP variation in sorghum, and 57,411 recombination events. Notably, recombination events were four- to fivefold enriched in coding sequences and 5' untranslated regions of genes. To test the power of the NAM population for trait dissection, we conducted joint linkage mapping for two major adaptive traits, flowering time and plant height. We precisely mapped several known genes for these two traits, and identified several additional QTL. Considering all SNPs simultaneously, genetic variation accounted for 65% of flowering time variance and 75% of plant height variance. Further, we directly compared NAM to genome-wide association mapping (using panels of the same size) and found that flowering time and plant height QTL were more consistently identified with the NAM population. Finally, for simulated QTL under strong selection in diversity panels, the power of QTL detection was up to three times greater for NAM vs. association mapping with a diverse panel. These findings validate the NAM resource for trait mapping in sorghum, and demonstrate the value of NAM for dissection of adaptive traits.

Genomic characterization of a core set of the USDA-NPGS Ethiopian sorghum germplasm collection: implications for germplasm conservation, evaluation, and utilization in crop improvement

BACKGROUND

The USDA Agriculture Research Service National Plant Germplasm System (NPGS) preserves the largest sorghum germplasm collection in the world, which includes 7,217 accessions from the center of diversity in Ethiopia. The characterization of this exotic germplasm at a genome-wide scale will improve conservation efforts and its utilization in research and breeding programs. Therefore, we phenotyped a representative core set of 374 Ethiopian accessions at two locations for agronomic traits and characterized the genomes.

RESULTS

Using genotyping-by-sequencing, we identified 148,476 single-nucleotide polymorphism (SNP) markers distributed across the entire genome. Over half of the alleles were rare (frequency < 0.05). The genetic profile of each accession was unique (i.e., no duplicates), and the average genetic distance among accessions was 0.70. Based on population structure and cluster analyses, we separated the collection into 11 populations with pairwise F _ST values ranging from 0.11 to 0.47. In total, 198 accessions (53%) were assigned to one of these populations with an ancestry membership coefficient of larger than 0.60; these covered 90% of the total genomic variation. We characterized these populations based on agronomic and seed compositional traits. We performed a cluster analysis with the sorghum association panel based on 26,026 SNPs and determined that nine of the Ethiopian populations expanded the genetic diversity in the panel. Genome-wide association analysis demonstrated that these low-coverage data and the observed population structure could be employed for the genomic dissection of important phenotypes in this core set of Ethiopian sorghum germplasm.

CONCLUSIONS

The NPGS Ethiopian sorghum germplasm is a genetically and phenotypically diverse collection comprising 11 populations with high levels of admixture. Genetic associations with agronomic traits can be used to improve the screening of exotic germplasm for selection of specific populations. We detected many rare alleles, suggesting that this germplasm contains potentially useful undiscovered alleles, but their discovery and characterization will require extensive effort. The genotypic data available for these accessions provide a valuable resource for sorghum breeders and geneticists to effectively improve crops.

Genome-wide association analysis of seedling traits in diverse Sorghum germplasm under thermal stress

BACKGROUND

Climate variability due to fluctuation in temperature is a worldwide concern that imperils crop production. The need to understand how the germplasm variation in major crops can be utilized to aid in discovering and developing breeding lines that can withstand and adapt to temperature fluctuations is more necessary than ever. Here, we analyzed the genetic variation associated with responses to thermal stresses in a sorghum association panel (SAP) representing major races and working groups to identify single nucleotide polymorphisms (SNPs) that are associated with resilience to temperature stress in a major cereal crop.

RESULTS

The SAP exhibited extensive variation for seedling traits under cold and heat stress. Genome-wide analyses identified 30 SNPs that were strongly associated with traits measured at seedling stage under cold stress and tagged genes that act as regulators of anthocyanin expression and soluble carbohydrate metabolism. Meanwhile, 12 SNPs were significantly associated with seedling traits under heat stress and these SNPs tagged genes that function in sugar metabolism, and ion transport pathways. Evaluation of co-expression networks for genes near the significantly associated SNPs indicated complex gene interactions for cold and heat stresses in sorghum. We focused and validated the expression of four genes in the network of Sb06g025040, a basic-helix-loop-helix (bHLH) transcription factor that was proposed to be involved in purple color pigmentation of leaf, and observed that genes in this network were upregulated during cold stress in a moderately tolerant line as compared to the more sensitive line.

CONCLUSION

This study facilitated the tagging of genome regions associated with variation in seedling traits of sorghum under cold and heat stress. These findings show the potential of genotype information for development of temperature resilient sorghum cultivars and further characterization of genes and their networks responsible for adaptation to thermal stresses. Knowledge on the gene networks from this research can be extended to the other cereal crops to better understand the genetic basis of resilience to temperature fluctuations during plant developmental stages.

Genome-Wide Association Studies of Grain Yield Components in Diverse Sorghum Germplasm

Grain yield and its primary determinants, grain number and weight, are important traits in cereal crops that have been well studied; however, the genetic basis of and interactions between these traits remain poorly understood. Characterization of grain yield per primary panicle (YPP), grain number per primary panicle (GNP), and 1000-grain weight (TGW) in sorghum [ (L.) Moench], a hardy C cereal with a genome size of ∼730 Mb, was implemented in a diversity panel containing 390 accessions. These accessions were genotyped to obtain 268,830 single-nucleotide polymorphisms (SNPs). Genome-wide association studies (GWAS) were performed to identify loci associated with each grain yield component and understand the genetic interactions between these traits. Genome-wide association studies identified associations across the genome with YPP, GNP, and TGW that were located within previously mapped sorghum QTL for panicle weight, grain yield, and seed size, respectively. There were no significant associations between GNP and TGW that were within 100 kb, much greater than the average linkage disequilibrium (LD) in sorghum. The identification of nonoverlapping loci for grain number and weight suggests these traits may be manipulated independently to increase the grain yield of sorghum. Following GWAS, genomic regions surrounding each associated SNP were mined for candidate genes. Previously published expression data indicated several TGW candidate genes, including an ethylene receptor homolog, were primarily expressed within developing seed tissues to support GWAS. Furthermore, maize ( L.) homologs of identified TGW candidates were differentially expressed within the seed between small- and large-kernel lines from a segregating maize population.

A Genomic Resource for the Development, Improvement, and Exploitation of Sorghum for Bioenergy

With high productivity and stress tolerance, numerous grass genera of the Andropogoneae have emerged as candidates for bioenergy production. To optimize these candidates, research examining the genetic architecture of yield, carbon partitioning, and composition is required to advance breeding objectives. Significant progress has been made developing genetic and genomic resources for Andropogoneae, and advances in comparative and computational genomics have enabled research examining the genetic basis of photosynthesis, carbon partitioning, composition, and sink strength. To provide a pivotal resource aimed at developing a comparative understanding of key bioenergy traits in the Andropogoneae, we have established and characterized an association panel of 390 racially, geographically, and phenotypically diverse Sorghum bicolor accessions with 232,303 genetic markers. Sorghum bicolor was selected because of its genomic simplicity, phenotypic diversity, significant genomic tools, and its agricultural productivity and resilience. We have demonstrated the value of sorghum as a functional model for candidate gene discovery for bioenergy Andropogoneae by performing genome-wide association analysis for two contrasting phenotypes representing key components of structural and non-structural carbohydrates. We identified potential genes, including a cellulase enzyme and a vacuolar transporter, associated with increased non-structural carbohydrates that could lead to bioenergy sorghum improvement. Although our analysis identified genes with potentially clear functions, other candidates did not have assigned functions, suggesting novel molecular mechanisms for carbon partitioning traits. These results, combined with our characterization of phenotypic and genetic diversity and the public accessibility of each accession and genomic data, demonstrate the value of this resource and provide a foundation for future improvement of sorghum and related grasses for bioenergy production.

Exploiting Phenylpropanoid Derivatives to Enhance the Nutraceutical Values of Cereals and Legumes

Phenylpropanoids are a diverse chemical class with immense health benefits that are biosynthesized from the aromatic amino acid L-phenylalanine. This article reviews the progress for accessing variation in phenylpropanoids in germplasm collections, the genetic and molecular basis of phenylpropanoid biosynthesis, and the development of cultivars dense in seed-phenylpropanoids. Progress is also reviewed on high-throughput assays, factors that influence phenylpropanoids, the site of phenylpropanoids accumulation in seed, Genotype × Environment interactions, and on consumer attitudes for the acceptance of staple foods rich in phenylpropanoids. A paradigm shift was noted in barley, maize, rice, sorghum, soybean, and wheat, wherein cultivars rich in phenylpropanoids are grown in Europe and North and Central America. Studies have highlighted some biological constraints that need to be addressed for development of high-yielding cultivars that are rich in phenylpropanoids. Genomics-assisted breeding is expected to facilitate rapid introgression into improved genetic backgrounds by minimizing linkage drag. More research is needed to systematically characterize germplasm pools for assessing variation to support crop genetic enhancement, and assess consumer attitudes to foods rich in phenylpropanoids.

A Genome-Wide Association Study on the Seedless Phenotype in Banana (Musa spp.) Reveals the Potential of a Selected Panel to Detect Candidate Genes in a Vegetatively Propagated Crop

Banana (Musa sp.) is a vegetatively propagated, low fertility, potentially hybrid and polyploid crop. These qualities make the breeding and targeted genetic improvement of this crop a difficult and long process. The Genome-Wide Association Study (GWAS) approach is becoming widely used in crop plants and has proven efficient to detecting candidate genes for traits of interest, especially in cereals. GWAS has not been applied yet to a vegetatively propagated crop. However, successful GWAS in banana would considerably help unravel the genomic basis of traits of interest and therefore speed up this crop improvement. We present here a dedicated panel of 105 accessions of banana, freely available upon request, and their corresponding GBS data. A set of 5,544 highly reliable markers revealed high levels of admixture in most accessions, except for a subset of 33 individuals from Papua. A GWAS on the seedless phenotype was then successfully applied to the panel. By applying the Mixed Linear Model corrected for both kinship and structure as implemented in TASSEL, we detected 13 candidate genomic regions in which we found a number of genes potentially linked with the seedless phenotype (i.e. parthenocarpy combined with female sterility). An additional GWAS performed on the unstructured Papuan subset composed of 33 accessions confirmed six of these regions as candidate. Out of both sets of analyses, one strong candidate gene for female sterility, a putative orthologous gene to Histidine Kinase CKI1, was identified. The results presented here confirmed the feasibility and potential of GWAS when applied to small sets of banana accessions, at least for traits underpinned by a few loci. As phenotyping in banana is extremely space and time-consuming, this latest finding is of particular importance in the context of banana improvement.

The Sorghum Gene for Leaf Color Changes upon Wounding (P) Encodes a Flavanone 4-Reductase in the 3-Deoxyanthocyanidin Biosynthesis Pathway

Upon wounding or pathogen invasion, leaves of sorghum [Sorghum bicolor (L.) Moench] plants with the P gene turn purple, whereas leaves with the recessive allele turn brown or tan. This purple phenotype is determined by the production of two 3-deoxyanthocyanidins, apigeninidin and luteolinidin, which are not produced by the tan-phenotype plants. Using map-based cloning in progeny from a cross between purple Nakei-MS3B (PP) and tan Greenleaf (pp) cultivars, we isolated this gene, which was located in a 27-kb genomic region around the 58.1 Mb position on chromosome 6. Four candidate genes identified in this region were similar to the maize leucoanthocyanidin reductase gene. None of them was expressed before wounding, and only the Sb06g029550 gene was induced in both cultivars after wounding. The Sb06g029550 protein was detected in Nakei-MS3B, but only slightly in Greenleaf, in which it may be unstable because of a Cys252Tyr substitution. A recombinant Sb06g029550 protein had a specific flavanone 4-reductase activity, and converted flavanones (naringenin or eriodictyol) to flavan-4-ols (apiforol or luteoforol) in vitro Our data indicate that the Sb06g029550 gene is involved in the 3-deoxyanthocyanidin synthesis pathway.

The Sorghum Gene for Leaf Color Changes upon Wounding (P) Encodes a Flavanone 4-Reductase in the 3-Deoxyanthocyanidin Biosynthesis Pathway

Upon wounding or pathogen invasion, leaves of sorghum [Sorghum bicolor (L.) Moench] plants with the P gene turn purple, whereas leaves with the recessive allele turn brown or tan. This purple phenotype is determined by the production of two 3-deoxyanthocyanidins, apigeninidin and luteolinidin, which are not produced by the tan-phenotype plants. Using map-based cloning in progeny from a cross between purple Nakei-MS3B (PP) and tan Greenleaf (pp) cultivars, we isolated this gene, which was located in a 27-kb genomic region around the 58.1 Mb position on chromosome 6. Four candidate genes identified in this region were similar to the maize leucoanthocyanidin reductase gene. None of them was expressed before wounding, and only the Sb06g029550 gene was induced in both cultivars after wounding. The Sb06g029550 protein was detected in Nakei-MS3B, but only slightly in Greenleaf, in which it may be unstable because of a Cys252Tyr substitution. A recombinant Sb06g029550 protein had a specific flavanone 4-reductase activity, and converted flavanones (naringenin or eriodictyol) to flavan-4-ols (apiforol or luteoforol) in vitro. Our data indicate that the Sb06g029550 gene is involved in the 3-deoxyanthocyanidin synthesis pathway.

Transcriptome profiling and validation of gene based single nucleotide polymorphisms (SNPs) in sorghum genotypes with contrasting responses to cold stress

BACKGROUND

Sorghum is a versatile cereal crop, with excellent heat and drought tolerance. However, it is susceptible to early-season cold stress (12-15 °C) which limits stand-establishment and seedling growth. To gain further insights on the molecular mechanism of cold tolerance in sorghum we performed transcriptome profiling between known cold sensitive and tolerant sorghum lines using RNA sequencing technology under control and cold stress treatments.

RESULTS

Here we report on the identification of differentially expressed genes (DEGs) between contrasting sorghum genotypes, HongkeZi (cold tolerant) and BTx623 (cold sensitive) under cool and control temperatures using RNAseq approach to elucidate the molecular basis of sorghum response to cold stress. Furthermore, we validated bi-allelic variants in the form of single nucleotide polymorphism (SNPs) between the cold susceptible and tolerant lines of sorghum. An analysis of transcriptome profile showed that in response to cold, a total of 1910 DEGs were detected under cold and control temperatures in both genotypes. We identified a subset of genes under cold stress for downstream analysis, including transcription factors that exhibit differential abundance between the sensitive and tolerant genotypes. We identified transcription factors including Dehydration-responsive element-binding factors, C-repeat binding factors, and Ethylene responsive transcription factors as significantly upregulated during cold stress in cold tolerant HKZ. Additionally, specific genes such as plant cytochromes, glutathione s-transferases, and heat shock proteins were found differentially regulated under cold stress between cold tolerant and susceptible genotype of sorghum. A total of 41,603 SNP were identified between the cold sensitive and tolerant genotypes with minimum read of four. Approximately 89 % of the 114 SNP sites selected for evaluation were validated using endpoint genotyping technology.

CONCLUSION

A new strategy which involved an integrated analysis of differential gene expression and identification of bi-allelic single nucleotide polymorphism (SNP) was conducted to determine and analyze differentially expressed genes and variation involved in cold stress response of sorghum. The results gathered provide an insight into the complex mechanisms associated with cold response in sorghum, which involve an array of transcription factors and genes which were previously related to abiotic stress response. This study also offers resource for gene based SNP that can be applied towards targeted genomic studies of cold tolerance in sorghum and other cereal crops.

Genome-environment associations in sorghum landraces predict adaptive traits

Improving environmental adaptation in crops is essential for food security under global change, but phenotyping adaptive traits remains a major bottleneck. If associations between single-nucleotide polymorphism (SNP) alleles and environment of origin in crop landraces reflect adaptation, then these could be used to predict phenotypic variation for adaptive traits. We tested this proposition in the global food crop Sorghum bicolor, characterizing 1943 georeferenced landraces at 404,627 SNPs and quantifying allelic associations with bioclimatic and soil gradients. Environment explained a substantial portion of SNP variation, independent of geographical distance, and genic SNPs were enriched for environmental associations. Further, environment-associated SNPs predicted genotype-by-environment interactions under experimental drought stress and aluminum toxicity. Our results suggest that genomic signatures of environmental adaptation may be useful for crop improvement, enhancing germplasm identification and marker-assisted selection. Together, genome-environment associations and phenotypic analyses may reveal the basis of environmental adaptation.

Genome-wide association study of grain polyphenol concentrations in global sorghum [Sorghum bicolor (L.) Moench] germplasm

Identifying natural variation of health-promoting compounds in staple crops and characterizing its genetic basis can help improve human nutrition through crop biofortification. Some varieties of sorghum, a staple cereal crop grown worldwide, have high concentrations of proanthocyanidins and 3-deoxyanthocyanidins, polyphenols with antioxidant and anti-inflammatory properties. We quantified total phenols, proanthocyanidins, and 3-deoxyanthocyanidins in a global sorghum diversity panel (n = 381) using near-infrared spectroscopy (NIRS), and characterized the patterns of variation with respect to geographic origin and botanical race. A genome-wide association study (GWAS) with 404,628 SNP markers identified novel quantitative trait loci for sorghum polyphenols, some of which colocalized with homologues of flavonoid pathway genes from other plants, including an orthologue of maize (Zea mays) Pr1 and a homologue of Arabidopsis (Arabidopsis thaliana) TT16. This survey of grain polyphenol variation in sorghum germplasm and catalog of flavonoid pathway loci may be useful to guide future enhancement of cereal polyphenols.