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

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.

Epistasis and pleiotropy-induced variation for plant breeding

Epistasis refers to nonallelic interaction between genes that cause bias in estimates of genetic parameters for a phenotype with interactions of two or more genes affecting the same trait. Partitioning of epistatic effects allows true estimation of the genetic parameters affecting phenotypes. Multigenic variation plays a central role in the evolution of complex characteristics, among which pleiotropy, where a single gene affects several phenotypic characters, has a large influence. While pleiotropic interactions provide functional specificity, they increase the challenge of gene discovery and functional analysis. Overcoming pleiotropy-based phenotypic trade-offs offers potential for assisting breeding for complex traits. Modelling higher order nonallelic epistatic interaction, pleiotropy and non-pleiotropy-induced variation, and genotype × environment interaction in genomic selection may provide new paths to increase the productivity and stress tolerance for next generation of crop cultivars. Advances in statistical models, software and algorithm developments, and genomic research have facilitated dissecting the nature and extent of pleiotropy and epistasis. We overview emerging approaches to exploit positive (and avoid negative) epistatic and pleiotropic interactions in a plant breeding context, including developing avenues of artificial intelligence, novel exploitation of large-scale genomics and phenomics data, and involvement of genes with minor effects to analyse epistatic interactions and pleiotropic quantitative trait loci, including missing heritability.

Genomic co-localization of variation affecting agronomic and human gut microbiome traits in a meta-analysis of diverse sorghum

Substantial functional metabolic diversity exists within species of cultivated grain crops that directly or indirectly provide more than half of all calories consumed by humans around the globe. While such diversity is the molecular currency used for improving agronomic traits, diversity is poorly characterized for its effects on human nutrition and utilization by gut microbes. Moreover, we know little about agronomic traits' potential tradeoffs and pleiotropic effects on human nutritional traits. Here, we applied a quantitative genetics approach using a meta-analysis and parallel genome-wide association studies of Sorghum bicolor traits describing changes in the composition and function of human gut microbe communities, and any of 200 sorghum seed and agronomic traits across a diverse sorghum population to identify associated genetic variants. A total of 15 multiple-effect loci (MEL) were initially found where different alleles in the sorghum genome produced changes in seed that affected the abundance of multiple bacterial taxa across 2 human microbiomes in automated in vitro fermentations. Next, parallel genome-wide studies conducted for seed, biochemical, and agronomic traits in the same population identified significant associations within the boundaries of 13/15 MEL for microbiome traits. In several instances, the colocalization of variation affecting gut microbiome and agronomic traits provided hypotheses for causal mechanisms through which variation could affect both agronomic traits and human gut microbes. This work demonstrates that genetic factors affecting agronomic traits in sorghum seed can also drive significant effects on human gut microbes, particularly bacterial taxa considered beneficial. Understanding these pleiotropic relationships will inform future strategies for crop improvement toward yield, sustainability, and human health.

Nonphotochemical quenching kinetics GWAS in sorghum identifies genes that may play conserved roles in maize and Arabidopsis thaliana photoprotection

Photosynthetic organisms must cope with rapid fluctuations in light intensity. Nonphotochemical quenching (NPQ) enables the dissipation of excess light energy as heat under high light conditions, whereas its relaxation under low light maximizes photosynthetic productivity. We quantified variation in NPQ kinetics across a large sorghum (Sorghum bicolor) association panel in four environments, uncovering significant genetic control for NPQ. A genome-wide association study (GWAS) confidently identified three unique regions in the sorghum genome associated with NPQ and suggestive associations in an additional 61 regions. We detected strong signals from the sorghum ortholog of Arabidopsis thaliana Suppressor Of Variegation 3 (SVR3) involved in plastid-nucleus signaling. By integrating GWAS results for NPQ across maize (Zea mays) and sorghum-association panels, we identified a second gene, Non-yellowing 1 (NYE1), originally studied by Gregor Mendel in pea (Pisum sativum) and involved in the degradation of photosynthetic pigments in light-harvesting complexes. Analysis of nye1 insertion alleles in A. thaliana confirmed the effect of this gene on NPQ kinetics in eudicots. We extended our comparative genomics GWAS framework across the entire maize and sorghum genomes, identifying four additional loci involved in NPQ kinetics. These results provide a baseline for increasing the accuracy and speed of candidate gene identification for GWAS in species with high linkage disequilibrium.

Multi-trait association mapping for phosphorous efficiency reveals flexible root architectures in sorghum

BACKGROUND

On tropical regions, phosphorus (P) fixation onto aluminum and iron oxides in soil clays restricts P diffusion from the soil to the root surface, limiting crop yields. While increased root surface area favors P uptake under low-P availability, the relationship between the three-dimensional arrangement of the root system and P efficiency remains elusive. Here, we simultaneously assessed allelic effects of loci associated with a variety of root and P efficiency traits, in addition to grain yield under low-P availability, using multi-trait genome-wide association. We also set out to establish the relationship between root architectural traits assessed in hydroponics and in a low-P soil. Our goal was to better understand the influence of root morphology and architecture in sorghum performance under low-P availability.

RESULT

In general, the same alleles of associated SNPs increased root and P efficiency traits including grain yield in a low-P soil. We found that sorghum P efficiency relies on pleiotropic loci affecting root traits, which enhance grain yield under low-P availability. Root systems with enhanced surface area stemming from lateral root proliferation mostly up to 40 cm soil depth are important for sorghum adaptation to low-P soils, indicating that differences in root morphology leading to enhanced P uptake occur exactly in the soil layer where P is found at the highest concentration.

CONCLUSION

Integrated QTLs detected in different mapping populations now provide a comprehensive molecular genetic framework for P efficiency studies in sorghum. This indicated extensive conservation of P efficiency QTL across populations and emphasized the terminal portion of chromosome 3 as an important region for P efficiency in sorghum. Increases in root surface area via enhancement of lateral root development is a relevant trait for sorghum low-P soil adaptation, impacting the overall architecture of the sorghum root system. In turn, particularly concerning the critical trait for water and nutrient uptake, root surface area, root system development in deeper soil layers does not occur at the expense of shallow rooting, which may be a key reason leading to the distinctive sorghum adaptation to tropical soils with multiple abiotic stresses including low P availability and drought.

Strategies to identify and dissect trade-offs in plants

Trade-offs between traits arise and reflect constraints imposed by the environment and physicochemical laws. Trade-off situations are expected to be highly relevant for sessile plants, which have to respond to changes in the environment to ensure survival. Despite increasing interest in determining the genetic and molecular basis of plant trade-offs, there are still gaps and differences with respect to how trade-offs are defined, how they are measured, and how their genetic architecture is dissected. The first step to fill these gaps is to establish what is meant by trade-offs. In this review we provide a classification of the existing definitions of trade-offs according to: (1) the measures used for their quantification, (2) the dependence of trade-offs on environment, and (3) experimental designed used (i.e. a single individual across different environments or a population of individuals in single or multiple environments). We then compare the approaches for quantification of trade-offs based on phenotypic, between-individual, and genetic correlations, and stress the need for developing further quantification indices particularly for trade-offs between multiple traits. Lastly, we highlight the genetic mechanisms underpinning trade-offs and experimental designs that facilitate their discovery in plants, with focus on usage of natural variability. This review also offers a perspective for future research aimed at identification of plant trade-offs, dissection of their genetic architecture, and development of strategies to overcome trade-offs, with applications in crop breeding.

The genomic basis of nitrogen utilization efficiency and trait plasticity to improve nutrient stress tolerance in cultivated sunflower

Maintaining crop productivity is challenging as population growth, climate change, and increasing fertilizer costs necessitate expanding crop production to poorer lands whilst reducing inputs. Enhancing crops' nutrient use efficiency is thus an important goal, but requires a better understanding of related traits and their genetic basis. We investigated variation in low nutrient stress tolerance in a diverse panel of cultivated sunflower genotypes grown under high and low nutrient conditions, assessing relative growth rate (RGR) as performance. We assessed variation in traits related to nitrogen utilization efficiency (NUtE), mass allocation, and leaf elemental content. Across genotypes, nutrient limitation generally reduced RGR. Moreover, there was a negative correlation between vigor (RGR in control) and decline in RGR in response to stress. Given this trade-off, we focused on nutrient stress tolerance independent of vigor. This tolerance metric correlated with the change in NUtE, plasticity for a suite of morphological traits, and leaf element content. Genome-wide associations revealed regions associated with variation and plasticity in multiple traits, including two regions with seemingly additive effects on NUtE change. Our results demonstrate potential avenues for improving sunflower nutrient stress tolerance independent of vigor, and highlight specific traits and genomic regions that could play a role in enhancing tolerance.

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.

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.

A common resequencing-based genetic marker data set for global maize diversity

Maize (Zea mays ssp. mays) populations exhibit vast ranges of genetic and phenotypic diversity. As sequencing costs have declined, an increasing number of projects have sought to measure genetic differences between and within maize populations using whole-genome resequencing strategies, identifying millions of segregating single-nucleotide polymorphisms (SNPs) and insertions/deletions (InDels). Unlike older genotyping strategies like microarrays and genotyping by sequencing, resequencing should, in principle, frequently identify and score common genetic variants. However, in practice, different projects frequently employ different analytical pipelines, often employ different reference genome assemblies and consistently filter for minor allele frequency within the study population. This constrains the potential to reuse and remix data on genetic diversity generated from different projects to address new biological questions in new ways. Here, we employ resequencing data from 1276 previously published maize samples and 239 newly resequenced maize samples to generate a single unified marker set of approximately 366 million segregating variants and approximately 46 million high-confidence variants scored across crop wild relatives, landraces as well as tropical and temperate lines from different breeding eras. We demonstrate that the new variant set provides increased power to identify known causal flowering-time genes using previously published trait data sets, as well as the potential to track changes in the frequency of functionally distinct alleles across the global distribution of modern maize.

Gene-by-environment interactions in plants: Molecular mechanisms, environmental drivers, and adaptive plasticity

Plants demonstrate a broad range of responses to environmental shifts. One of the most remarkable responses is plasticity, which is the ability of a single plant genotype to produce different phenotypes in response to environmental stimuli. As with all traits, the ability of plasticity to evolve depends on the presence of underlying genetic diversity within a population. A common approach for evaluating the role of genetic variation in driving differences in plasticity has been to study genotype-by-environment interactions (G × E). G × E occurs when genotypes produce different phenotypic trait values in response to different environments. In this review, we highlight progress and promising methods for identifying the key environmental and genetic drivers of G × E. Specifically, methodological advances in using algorithmic and multivariate approaches to understand key environmental drivers combined with new genomic innovations can greatly increase our understanding about molecular responses to environmental stimuli. These developing approaches can be applied to proliferating common garden networks that capture broad natural environmental gradients to unravel the underlying mechanisms of G × E. An increased understanding of G × E can be used to enhance the resilience and productivity of agronomic systems.

Association mapping across a multitude of traits collected in diverse environments in maize

Classical genetic studies have identified many cases of pleiotropy where mutations in individual genes alter many different phenotypes. Quantitative genetic studies of natural genetic variants frequently examine one or a few traits, limiting their potential to identify pleiotropic effects of natural genetic variants. Widely adopted community association panels have been employed by plant genetics communities to study the genetic basis of naturally occurring phenotypic variation in a wide range of traits. High-density genetic marker data-18M markers-from 2 partially overlapping maize association panels comprising 1,014 unique genotypes grown in field trials across at least 7 US states and scored for 162 distinct trait data sets enabled the identification of of 2,154 suggestive marker-trait associations and 697 confident associations in the maize genome using a resampling-based genome-wide association strategy. The precision of individual marker-trait associations was estimated to be 3 genes based on a reference set of genes with known phenotypes. Examples were observed of both genetic loci associated with variation in diverse traits (e.g., above-ground and below-ground traits), as well as individual loci associated with the same or similar traits across diverse environments. Many significant signals are located near genes whose functions were previously entirely unknown or estimated purely via functional data on homologs. This study demonstrates the potential of mining community association panel data using new higher-density genetic marker sets combined with resampling-based genome-wide association tests to develop testable hypotheses about gene functions, identify potential pleiotropic effects of natural genetic variants, and study genotype-by-environment interaction.

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 F_st 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.

Comparing Deep Learning Approaches for Understanding Genotype × Phenotype Interactions in Biomass Sorghum

We explore the use of deep convolutional neural networks (CNNs) trained on overhead imagery of biomass sorghum to ascertain the relationship between single nucleotide polymorphisms (SNPs), or groups of related SNPs, and the phenotypes they control. We consider both CNNs trained explicitly on the classification task of predicting whether an image shows a plant with a reference or alternate version of various SNPs as well as CNNs trained to create data-driven features based on learning features so that images from the same plot are more similar than images from different plots, and then using the features this network learns for genetic marker classification. We characterize how efficient both approaches are at predicting the presence or absence of a genetic markers, and visualize what parts of the images are most important for those predictions. We find that the data-driven approaches give somewhat higher prediction performance, but have visualizations that are harder to interpret; and we give suggestions of potential future machine learning research and discuss the possibilities of using this approach to uncover unknown genotype × phenotype relationships.

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.

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.