Genetic and Environmental Patterns Underlying Phenotypic Plasticity in Flowering Time and Plant Height in Sorghum.

Phenotypic plasticity is the property of a genotype to produce different phenotypes under different environmental conditions. Understanding genetic and environmental factors behind phenotypic plasticity helps answer some longstanding biology questions and improve phenotype prediction. In this study, we investigated the phenotypic plasticity of flowering time and plant height with a set of diverse sorghum lines evaluated across 14 natural field environments. An environmental index was identified to quantitatively connect the environments. Reaction norms were then obtained with the identified indices for genetic dissection of phenotypic plasticity and performance prediction. Genome-wide association studies (GWAS) detected different sets of loci for reaction-norm parameters (intercept and slope), including 10 new genomic regions in addition to known maturity (Ma1) and dwarfing genes (Dw1, Dw2, Dw3, Dw4 and qHT7.1). Cross-validations under multiple scenarios showed promising results in predicting diverse germplasm in dynamic environments. Additional experiments conducted at four new environments, including one from a site outside of the geographical region of the initial environments, further validated the predictions. Our findings indicate that identifying the environmental index enriches our understanding of gene-environmental interplay underlying phenotypic plasticity, and that genomic prediction with the environmental dimension facilitates prediction-guided breeding for future environments.

Sorghum Association Panel whole-genome sequencing establishes cornerstone resource for dissecting genomic diversity.

Association mapping panels represent foundational resources for understanding the genetic basis of phenotypic diversity and serve to advance plant breeding by exploring genetic variation across diverse accessions. We report the whole-genome sequencing (WGS) of 400 sorghum (Sorghum bicolor (L.) Moench) accessions from the Sorghum Association Panel (SAP) at an average coverage of 38× (25-72×), enabling the development of a high-density genomic marker set of 43 983 694 variants including single-nucleotide polymorphisms (approximately 38 million), insertions/deletions (indels) (approximately 5 million), and copy number variants (CNVs) (approximately 170 000). We observe slightly more deletions among indels and a much higher prevalence of deletions among CNVs compared to insertions. This new marker set enabled the identification of several novel putative genomic associations for plant height and tannin content, which were not identified when using previous lower-density marker sets. WGS identified and scored variants in 5-kb bins where available genotyping-by-sequencing (GBS) data captured no variants, with half of all bins in the genome falling into this category. The predictive ability of genomic best unbiased linear predictor (GBLUP) models was increased by an average of 30% by using WGS markers rather than GBS markers. We identified 18 selection peaks across subpopulations that formed due to evolutionary divergence during domestication, and we found six Fst peaks resulting from comparisons between converted lines and breeding lines within the SAP that were distinct from the peaks associated with historic selection. This population has served and continues to serve as a significant public resource for sorghum research and demonstrates the value of improving upon existing genomic resources.

Variation in Root Exudate Composition Influences Soil Microbiome Membership and Function.

Root exudation is one of the primary processes that mediate interactions between plant roots, microorganisms, and the soil matrix, yet the mechanisms by which exudation alters microbial metabolism in soils have been challenging to unravel. Here, utilizing distinct sorghum genotypes, we characterized the chemical heterogeneity between root exudates and the effects of that variability on soil microbial membership and metabolism. Distinct exudate chemical profiles were quantified and used to formulate synthetic root exudate treatments: a high-organic-acid treatment (HOT) and a high-sugar treatment (HST). To parse the response of the soil microbiome to different exudate regimens, laboratory soil reactors were amended with these root exudate treatments as well as a nonexudate control. Amplicon sequencing of the 16S rRNA gene illustrated distinct microbial diversity patterns and membership in response to HST, HOT, or control amendments. Exometabolite changes reflected these microbial community changes, and we observed enrichment of organic and amino acids, as well as possible phytohormones in the HST relative to the HOT and control. Linking the metabolic capacity of metagenome-assembled genomes in the HST to the exometabolite patterns, we identified microorganisms that could produce these phytohormones. Our findings emphasize the tractability of high-resolution multiomics tools to investigate soil microbiomes, opening the possibility of manipulating native microbial communities to improve specific soil microbial functions and enhance crop production. IMPORTANCE Decrypting the chemical interactions between plant roots and the soil microbiome is a gateway for future manipulation and management of the rhizosphere, a soil compartment critical to promoting plant fitness and yields. Our experimental results demonstrate how soil microbial community and genomic diversity is influenced by root exudates of differing chemical compositions and how changes in this microbiome result in altered production of plant-relevant metabolites. Together, these findings demonstrate the tractability of high-resolution multiomics tools to investigate soil microbiomes and provide new information on plant-soil environments useful for the development of efficient and precise microbiota management strategies in agricultural systems.

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 C4 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 new reference genome for Sorghum bicolor reveals high levels of sequence similarity between sweet and grain genotypes: implications for the genetics of sugar metabolism.

BACKGROUND: The process of crop domestication often consists of two stages: initial domestication, where the wild species is first cultivated by humans, followed by diversification, when the domesticated species are subsequently adapted to more environments and specialized uses. Selective pressure to increase sugar accumulation in certain varieties of the cereal crop Sorghum bicolor is an excellent example of the latter; this has resulted in pronounced phenotypic divergence between sweet and grain-type sorghums, but the genetic mechanisms underlying these differences remain poorly understood. RESULTS: Here we present a new reference genome based on an archetypal sweet sorghum line and compare it to the current grain sorghum reference, revealing a high rate of nonsynonymous and potential loss of function mutations, but few changes in gene content or overall genome structure. We also use comparative transcriptomics to highlight changes in gene expression correlated with high stalk sugar content and show that changes in the activity and possibly localization of transporters, along with the timing of sugar metabolism play a critical role in the sweet phenotype. CONCLUSIONS: The high level of genomic similarity between sweet and grain sorghum reflects their historical relatedness, rather than their current phenotypic differences, but we find key changes in signaling molecules and transcriptional regulators that represent new candidates for understanding and improving sugar metabolism in this important crop.

Genetic architecture of kernel composition in global sorghum germplasm.

BACKGROUND: Sorghum [Sorghum bicolor (L.) Moench] is an important cereal crop for dryland areas in the United States and for small-holder farmers in Africa. Natural variation of sorghum grain composition (protein, fat, and starch) between accessions can be used for crop improvement, but the genetic controls are still unresolved. The goals of this study were to quantify natural variation of sorghum grain composition and to identify single-nucleotide polymorphisms (SNPs) associated with variation in grain composition concentrations. RESULTS: In this study, we quantified protein, fat, and starch in a global sorghum diversity panel using near-infrared spectroscopy (NIRS). Protein content ranged from 8.1 to 18.8%, fat content ranged from 1.0 to 4.3%, and starch content ranged from 61.7 to 71.1%. Durra and bicolor-durra sorghum from Ethiopia and India had the highest protein and fat and the lowest starch content, while kafir sorghum from USA, India, and South Africa had the lowest protein and the highest starch content. Genome-wide association studies (GWAS) identified quantitative trait loci (QTL) for sorghum protein, fat, and starch. Previously published RNAseq data was used to identify candidate genes within a GWAS QTL region. A putative alpha-amylase 3 gene, which has previously been shown to be associated with grain composition traits, was identified as a strong candidate for protein and fat variation. CONCLUSIONS: We identified promising sources of genetic material for manipulation of grain composition traits, and several loci and candidate genes that may control sorghum grain composition. This survey of grain composition in sorghum germplasm and identification of protein, fat, and starch QTL contributes to our understanding of the genetic basis of natural variation in sorghum grain nutritional traits.

Integration of Experiments across Diverse Environments Identifies the Genetic Determinants of Variation in Sorghum bicolor Seed Element Composition.

Seedling establishment and seed nutritional quality require the sequestration of sufficient element nutrients. The identification of genes and alleles that modify element content in the grains of cereals, including sorghum (Sorghum bicolor), is fundamental to developing breeding and selection methods aimed at increasing bioavailable element content and improving crop growth. We have developed a high-throughput work flow for the simultaneous measurement of multiple elements in sorghum seeds. We measured seed element levels in the genotyped Sorghum Association Panel, representing all major cultivated sorghum races from diverse geographic and climatic regions, and mapped alleles contributing to seed element variation across three environments by genome-wide association. We observed significant phenotypic and genetic correlation between several elements across multiple years and diverse environments. The power of combining high-precision measurements with genome-wide association was demonstrated by implementing rank transformation and a multilocus mixed model to map alleles controlling 20 element traits, identifying 255 loci affecting the sorghum seed ionome. Sequence similarity to genes characterized in previous studies identified likely causative genes for the accumulation of zinc, manganese, nickel, calcium, and cadmium in sorghum seeds. In addition to strong candidates for these five elements, we provide a list of candidate loci for several other elements. Our approach enabled the identification of single-nucleotide polymorphisms in strong linkage disequilibrium with causative polymorphisms that can be evaluated in targeted selection strategies for plant breeding and improvement.

Genetic dissection of sorghum grain quality traits using diverse and segregating populations.

KEY MESSAGE: Coordinated association and linkage mapping identified 25 grain quality QTLs in multiple environments, and fine mapping of the Wx locus supports the use of high-density genetic markers in linkage mapping. There is a wide range of end-use products made from cereal grains, and these products often demand different grain characteristics. Fortunately, cereal crop species including sorghum [Sorghum bicolor (L.) Moench] contain high phenotypic variation for traits influencing grain quality. Identifying genetic variants underlying this phenotypic variation allows plant breeders to develop genotypes with grain attributes optimized for their intended usage. Multiple sorghum mapping populations were rigorously phenotyped across two environments (SC Coastal Plain and Central TX) in 2 years for five major grain quality traits: amylose, starch, crude protein, crude fat, and gross energy. Coordinated association and linkage mapping revealed several robust QTLs that make prime targets to improve grain quality for food, feed, and fuel products. Although the amylose QTL interval spanned many megabases, the marker with greatest significance was located just 12 kb from waxy (Wx), the primary gene regulating amylose production in cereal grains. This suggests higher resolution mapping in recombinant inbred line (RIL) populations can be obtained when genotyped at a high marker density. The major QTL for crude fat content, identified in both a RIL population and grain sorghum diversity panel, encompassed the DGAT1 locus, a critical gene involved in maize lipid biosynthesis. Another QTL on chromosome 1 was consistently mapped in both RIL populations for multiple grain quality traits including starch, crude protein, and gross energy. Collectively, these genetic regions offer excellent opportunities to manipulate grain composition and set up future studies for gene validation.

The Association of Neighborhood Gene-Environment Susceptibility with Cortisol and Blood Pressure in African-American Adults.

BACKGROUND: African-American adults are disproportionately affected by stress-related chronic conditions like high blood pressure (BP), and both environmental stress and genetic risk may play a role in its development. PURPOSE: This study tested whether the dual risk of low neighborhood socioeconomic status (SES) and glucocorticoid genetic sensitivity interacted to predict waking cortisol and BP. METHODS: Cross-sectional waking cortisol and BP were collected from 208 African-American adults who were participating in a follow-up visit as part of the Positive Action for Today's Health trial. Three single-nucleotide polymorphisms were genotyped, salivary cortisol samples were collected, and neighborhood SES was calculated using 2010 Census data. RESULTS: The sample was mostly female (65 %), with weight classified as overweight or obese (M BMI = 32.74, SD = 8.88) and a mean age of 55.64 (SD = 15.21). The gene-by-neighborhood SES interaction predicted cortisol (B = 0.235, p = .001, r (2) = .036), but not BP. For adults with high genetic sensitivity, waking cortisol was lower with lower SES but higher with higher SES (B = 0.87). Lower neighborhood SES was also related to higher systolic BP (B = -0.794, p = .028). CONCLUSIONS: Findings demonstrated an interaction whereby African-American adults with high genetic sensitivity had high levels of waking cortisol with higher neighborhood SES, and low levels with lower neighborhood SES. This moderation effect is consistent with a differential susceptibility gene-environment pattern, rather than a dual-risk pattern. These findings contribute to a growing body of evidence that demonstrates the importance of investigating complex gene-environment relations in order to better understand stress-related health disparities.

Population genomic and genome-wide association studies of agroclimatic traits in sorghum.

Accelerating crop improvement in sorghum, a staple food for people in semiarid regions across the developing world, is key to ensuring global food security in the context of climate change. To facilitate gene discovery and molecular breeding in sorghum, we have characterized ~265,000 single nucleotide polymorphisms (SNPs) in 971 worldwide accessions that have adapted to diverse agroclimatic conditions. Using this genome-wide SNP map, we have characterized population structure with respect to geographic origin and morphological type and identified patterns of ancient crop diffusion to diverse agroclimatic regions across Africa and Asia. To better understand the genomic patterns of diversification in sorghum, we quantified variation in nucleotide diversity, linkage disequilibrium, and recombination rates across the genome. Analyzing nucleotide diversity in landraces, we find evidence of selective sweeps around starch metabolism genes, whereas in landrace-derived introgression lines, we find introgressions around known height and maturity loci. To identify additional loci underlying variation in major agroclimatic traits, we performed genome-wide association studies (GWAS) on plant height components and inflorescence architecture. GWAS maps several classical loci for plant height, candidate genes for inflorescence architecture. Finally, we trace the independent spread of multiple haplotypes carrying alleles for short stature or long inflorescence branches. This genome-wide map of SNP variation in sorghum provides a basis for crop improvement through marker-assisted breeding and genomic selection.

The Sorghum bicolor genome and the diversification of grasses.

Sorghum, an African grass related to sugar cane and maize, is grown for food, feed, fibre and fuel. We present an initial analysis of the approximately 730-megabase Sorghum bicolor (L.) Moench genome, placing approximately 98% of genes in their chromosomal context using whole-genome shotgun sequence validated by genetic, physical and syntenic information. Genetic recombination is largely confined to about one-third of the sorghum genome with gene order and density similar to those of rice. Retrotransposon accumulation in recombinationally recalcitrant heterochromatin explains the approximately 75% larger genome size of sorghum compared with rice. Although gene and repetitive DNA distributions have been preserved since palaeopolyploidization approximately 70 million years ago, most duplicated gene sets lost one member before the sorghum-rice divergence. Concerted evolution makes one duplicated chromosomal segment appear to be only a few million years old. About 24% of genes are grass-specific and 7% are sorghum-specific. Recent gene and microRNA duplications may contribute to sorghum's drought tolerance.

A unified mixed-model method for association mapping that accounts for multiple levels of relatedness.

As population structure can result in spurious associations, it has constrained the use of association studies in human and plant genetics. Association mapping, however, holds great promise if true signals of functional association can be separated from the vast number of false signals generated by population structure. We have developed a unified mixed-model approach to account for multiple levels of relatedness simultaneously as detected by random genetic markers. We applied this new approach to two samples: a family-based sample of 14 human families, for quantitative gene expression dissection, and a sample of 277 diverse maize inbred lines with complex familial relationships and population structure, for quantitative trait dissection. Our method demonstrates improved control of both type I and type II error rates over other methods. As this new method crosses the boundary between family-based and structured association samples, it provides a powerful complement to currently available methods for association mapping.

A Sorghum bicolor expression atlas reveals dynamic genotype-specific expression profiles for vegetative tissues of grain, sweet and bioenergy sorghums.

BACKGROUND: Effective improvement in sorghum crop development necessitates a genomics-based approach to identify functional genes and QTLs. Sequenced in 2009, a comprehensive annotation of the sorghum genome and the development of functional genomics resources is key to enable the discovery and deployment of regulatory and metabolic genes and gene networks for crop improvement. RESULTS: This study utilizes the first commercially available whole-transcriptome sorghum microarray (Sorgh-WTa520972F) to identify tissue and genotype-specific expression patterns for all identified Sorghum bicolor exons and UTRs. The genechip contains 1,026,373 probes covering 149,182 exons (27,577 genes) across the Sorghum bicolor nuclear, chloroplast, and mitochondrial genomes. Specific probesets were also included for putative non-coding RNAs that may play a role in gene regulation (e.g., microRNAs), and confirmed functional small RNAs in related species (maize and sugarcane) were also included in our array design. We generated expression data for 78 samples with a combination of four different tissue types (shoot, root, leaf and stem), two dissected stem tissues (pith and rind) and six diverse genotypes, which included 6 public sorghum lines (R159, Atlas, Fremont, PI152611, AR2400 and PI455230) representing grain, sweet, forage, and high biomass ideotypes. CONCLUSIONS: Here we present a summary of the microarray dataset, including analysis of tissue-specific gene expression profiles and associated expression profiles of relevant metabolic pathways. With an aim to enable identification and functional characterization of genes in sorghum, this expression atlas presents a new and valuable resource to the research community.

Identification of pleiotropic loci mediating structural and non-structural carbohydrate accumulation within the sorghum bioenergy association panel using high-throughput markers.

Molecular characterization of diverse germplasm can contribute to breeding programs by increasing genetic gain for sorghum [Sorghum bicolor (L.) Moench] improvement. Identifying novel marker-trait associations and candidate genes enriches the existing genomic resources and can improve bioenergy-related traits using genomic-assisted breeding. In the current scenario, identifying the genetic loci underlying biomass and carbon partitioning is vital for ongoing efforts to maximize each carbon sink's yield for bioenergy production. Here, we have processed a high-density genomic marker (22 466 550) data based on whole-genome sequencing (WGS) using a set of 365 accessions from the bioenergy association panel (BAP), which includes ~19.7 million (19 744 726) single nucleotide polymorphism (SNPs) and 2.7 million (~2 721 824) insertion deletions (indels). A set of high-quality filtered SNP (~5.48 million) derived markers facilitated the assessment of population structure, genetic diversity, and genome-wide association studies (GWAS) for various traits related to biomass and its composition using the BAP. The phenotypic traits for GWAS included seed color (SC), plant height (PH), days to harvest (DTH), fresh weight (FW), dry weight (DW), brix content % (BRX), neutral detergent fiber (NDF), acid detergent fiber (ADF), non-fibrous carbohydrate (NFC), and lignin content. Several novel loci and candidate genes were identified for bioenergy-related traits, and some well-characterized genes for plant height (Dw1 and Dw2) and the YELLOW SEED1 locus (Y1) were validated. We further performed a multi-variate adaptive shrinkage analysis to identify pleiotropic QTL, which resulted in several shared marker-trait associations among bioenergy and compositional traits. Significant marker-trait associations with pleiotropic effects can be used to develop molecular markers for trait improvement using a marker-assisted breeding approach. Significant nucleotide diversity and heterozygosity were observed between photoperiod-sensitive and insensitive individuals of the panel. This diverse bioenergy panel with genomic resources will provide an excellent opportunity for further genetic studies, including selecting parental lines for superior hybrid development to improve biomass-related traits in sorghum.

Functional genomic effects of indels using Bayesian genome-phenome wide association studies in sorghum.

High-throughput genomic and phenomic data have enhanced the ability to detect genotype-to-phenotype associations that can resolve broad pleiotropic effects of mutations on plant phenotypes. As the scale of genotyping and phenotyping has advanced, rigorous methodologies have been developed to accommodate larger datasets and maintain statistical precision. However, determining the functional effects of associated genes/loci is expensive and limited due to the complexity associated with cloning and subsequent characterization. Here, we utilized phenomic imputation of a multi-year, multi-environment dataset using PHENIX which imputes missing data using kinship and correlated traits, and we screened insertions and deletions (InDels) from the recently whole-genome sequenced Sorghum Association Panel for putative loss-of-function effects. Candidate loci from genome-wide association results were screened for potential loss of function using a Bayesian Genome-Phenome Wide Association Study (BGPWAS) model across both functionally characterized and uncharacterized loci. Our approach is designed to facilitate in silico validation of associations beyond traditional candidate gene and literature-search approaches and to facilitate the identification of putative variants for functional analysis and reduce the incidence of false-positive candidates in current functional validation methods. Using this Bayesian GPWAS model, we identified associations for previously characterized genes with known loss-of-function alleles, specific genes falling within known quantitative trait loci, and genes without any previous genome-wide associations while additionally detecting putative pleiotropic effects. In particular, we were able to identify the major tannin haplotypes at the Tan1 locus and effects of InDels on the protein folding. Depending on the haplotype present, heterodimer formation with Tan2 was significantly affected. We also identified major effect InDels in Dw2 and Ma1, where proteins were truncated due to frameshift mutations that resulted in early stop codons. These truncated proteins also lost most of their functional domains, suggesting that these indels likely result in loss of function. Here, we show that the Bayesian GPWAS model is able to identify loss-of-function alleles that can have significant effects upon protein structure and folding as well as multimer formation. Our approach to characterize loss-of-function mutations and their functional repercussions will facilitate precision genomics and breeding by identifying key targets for gene editing and trait integration.

Genomic prediction of hybrid performance for agronomic traits in sorghum.

Hybrid breeding in sorghum [Sorghum bicolor (L.) Moench] utilizes the cytoplasmic-nuclear male sterility (CMS) system for seed production and subsequently harnesses heterosis. Since the cost of developing and evaluating inbred and hybrid lines in the CMS system is costly and time-consuming, genomic prediction of parental lines and hybrids is based on genetic data genotype. We generated 602 hybrids by crossing two female (A) lines with 301 diverse and elite male (R) lines from the sorghum association panel and collected phenotypic data for agronomic traits over two years. We genotyped the inbred parents using whole genome resequencing and used 2,687,342 high quality (minor allele frequency > 2%) single nucleotide polymorphisms for genomic prediction. For grain yield, the experimental hybrids exhibited an average mid-parent heterosis of 40%. Genomic best linear unbiased prediction (GBLUP) for hybrid performance yielded an average prediction accuracy of 0.76-0.93 under the prediction scenario where both parental lines in validation sets were included in the training sets (T2). However, when only female tester was shared between training and validation sets (T1F), prediction accuracies declined by 12-90%, with plant height showing the greatest decline. Mean accuracies for predicting the general combining ability of male parents ranged from 0.33 to 0.62 for all traits. Our results showed hybrid performance for agronomic traits can be predicted with high accuracy, and optimizing genomic relationship is essential for optimal training population design for genomic selection in sorghum breeding.

Development and characterization of a sorghum multi-parent advanced generation intercross (MAGIC) population for capturing diversity among seed parent gene pool.

Multiparent advanced eneration inter-cross (MAGIC) populations improve the precision of quantitative trait loci (QTL) mapping over biparental populations by incorporating increased diversity and opportunities to reduce linkage disequilibrium among variants. Here, we describe the development of a MAGIC B-Line (MBL) population from an inter-cross among 4 diverse founders of grain sorghum [Sorghum bicolor (L.) Moench] across different races (kafir, guinea, durra, and caudatum). These founders were selected based on genetic uniqueness and several distinct qualitative features including panicle architecture, plant color, seed color, endosperm texture, and awns. A whole set of MBL (708 F6) recombinant inbred lines along with their founders were genotyped using Diversity Arrays Technology (DArTseq) and 5,683 single-nucleotide polymorphisms (SNPs) were generated. A genetic linkage map was constructed using a set of polymorphic, quality-filtered markers (2,728 SNPs) for QTL interval-mapping. For population validation, 3 traits (seed color, plant color, and awns) were used for QTL mapping and genome-wide association study (GWAS). QTL mapping and GWAS identified 4 major genomic regions located across 3 chromosomes (Chr1, Chr3, and Chr6) that correspond to known genetic loci for the targeted traits. Founders of this population consist of the fertility maintainer (A/B line) gene pool and derived MBL lines could serve as female/seed parents in the cytoplasmic male sterility breeding system. The MBL population will serve as a unique genetic and genomic resource to better characterize the genetics of complex traits and potentially identify superior alleles for crop improvement efforts to enrich the seed parent gene pool.

Discovering useful genetic variation in the seed parent gene pool for sorghum improvement.

Multi-parent populations contain valuable genetic material for dissecting complex, quantitative traits and provide a unique opportunity to capture multi-allelic variation compared to the biparental populations. A multi-parent advanced generation inter-cross (MAGIC) B-line (MBL) population composed of 708 F6 recombinant inbred lines (RILs), was recently developed from four diverse founders. These selected founders strategically represented the four most prevalent botanical races (kafir, guinea, durra, and caudatum) to capture a significant source of genetic variation to study the quantitative traits in grain sorghum [Sorghum bicolor (L.) Moench]. MBL was phenotyped at two field locations for seven yield-influencing traits: panicle type (PT), days to anthesis (DTA), plant height (PH), grain yield (GY), 1000-grain weight (TGW), tiller number per meter (TN) and yield per panicle (YPP). High phenotypic variation was observed for all the quantitative traits, with broad-sense heritabilities ranging from 0.34 (TN) to 0.84 (PH). The entire population was genotyped using Diversity Arrays Technology (DArTseq), and 8,800 single nucleotide polymorphisms (SNPs) were generated. A set of polymorphic, quality-filtered markers (3,751 SNPs) and phenotypic data were used for genome-wide association studies (GWAS). We identified 52 marker-trait associations (MTAs) for the seven traits using BLUPs generated from replicated plots in two locations. We also identified desirable allelic combinations based on the plant height loci (Dw1, Dw2, and Dw3), which influences yield related traits. Additionally, two novel MTAs were identified each on Chr1 and Chr7 for yield traits independent of dwarfing genes. We further performed a multi-variate adaptive shrinkage analysis and 15 MTAs with pleiotropic effect were identified. The five best performing MBL progenies were selected carrying desirable allelic combinations. Since the MBL population was designed to capture significant diversity for maintainer line (B-line) accessions, these progenies can serve as valuable resources to develop superior sorghum hybrids after validation of their general combining abilities via crossing with elite pollinators. Further, newly identified desirable allelic combinations can be used to enrich the maintainer germplasm lines through marker-assisted backcross breeding.

Dissecting the Genetic Architecture of Carbon Partitioning in Sorghum Using Multiscale Phenotypes.

Carbon partitioning in plants may be viewed as a dynamic process composed of the many interactions between sources and sinks. The accumulation and distribution of fixed carbon is not dictated simply by the sink strength and number but is dependent upon the source, pathways, and interactions of the system. As such, the study of carbon partitioning through perturbations to the system or through focus on individual traits may fail to produce actionable developments or a comprehensive understanding of the mechanisms underlying this complex process. Using the recently published sorghum carbon-partitioning panel, we collected both macroscale phenotypic characteristics such as plant height, above-ground biomass, and dry weight along with microscale compositional traits to deconvolute the carbon-partitioning pathways in this multipurpose crop. Multivariate analyses of traits resulted in the identification of numerous loci associated with several distinct carbon-partitioning traits, which putatively regulate sugar content, manganese homeostasis, and nitrate transportation. Using a multivariate adaptive shrinkage approach, we identified several loci associated with multiple traits suggesting that pleiotropic and/or interactive effects may positively influence multiple carbon-partitioning traits, or these overlaps may represent molecular switches mediating basal carbon allocating or partitioning networks. Conversely, we also identify a carbon tradeoff where reduced lignin content is associated with increased sugar content. The results presented here support previous studies demonstrating the convoluted nature of carbon partitioning in sorghum and emphasize the importance of taking a holistic approach to the study of carbon partitioning by utilizing multiscale phenotypes.

Genomic patterns of structural variation among diverse genotypes of Sorghum bicolor and a potential role for deletions in local adaptation.

Genomic structural mutations, especially deletions, are an important source of variation in many species and can play key roles in phenotypic diversification and evolution. Previous work in many plant species has identified multiple instances of structural variations (SVs) occurring in or near genes related to stress response and disease resistance, suggesting a possible role for SVs in local adaptation. Sorghum [Sorghum bicolor (L.) Moench] is one of the most widely grown cereal crops in the world. It has been adapted to an array of different climates as well as bred for multiple purposes, resulting in a striking phenotypic diversity. In this study, we identified genome-wide SVs in the Biomass Association Panel, a collection of 347 diverse sorghum genotypes collected from multiple countries and continents. Using Illumina-based, short-read whole-genome resequencing data from every genotype, we found a total of 24,648 SVs, including 22,359 deletions. The global site frequency spectrum of deletions and other types of SVs fit a model of neutral evolution, suggesting that the majority of these mutations were not under any types of selection. Clustering results based on single nucleotide polymorphisms separated the genotypes into eight clusters which largely corresponded with geographic origins, with many of the large deletions we uncovered being unique to a single cluster. Even though most deletions appeared to be neutral, a handful of cluster-specific deletions were found in genes related to biotic and abiotic stress responses, supporting the possibility that at least some of these deletions contribute to local adaptation in sorghum.