Imputation of 3 million SNPs in the Arabidopsis regional mapping population

Trait genetic architecture and population structure determine model selection for genomic prediction in natural Arabidopsis thaliana populations

Genomic prediction applies to any agro- or ecologically relevant traits, with distinct ontologies and genetic architectures. Selecting the most appropriate model for the distribution of genetic effects and their associated allele frequencies in the training population is crucial. Linear regression models are often preferred for genomic prediction. However, linear models may not suit all genetic architectures and training populations. Machine learning approaches have been proposed to improve genomic prediction owing to their capacity to capture complex biology including epistasis. However, the applicability of different genomic prediction models, including non-linear, non-parametric approaches, has not been rigorously assessed across a wide variety of plant traits in natural outbreeding populations. This study evaluates genomic prediction sensitivity to trait ontology and the impact of population structure on model selection and prediction accuracy. Examining 36 quantitative traits in 1,000+ natural genotypes of the model plant Arabidopsis thaliana, we assessed the performance of penalized regression, random forest, and multilayer perceptron at producing genomic predictions. Regression models were generally the most accurate, except for biochemical traits where random forest performed best. We link this result to the genetic architecture of each trait-notably that biochemical traits have simpler genetic architecture than macroscopic traits. Moreover, complex macroscopic traits, particularly those related to flowering time and yield, were strongly correlated to population structure, while molecular traits were better predicted by fewer, independent markers. This study highlights the relevance of machine learning approaches for simple molecular traits and underscores the need to consider ancestral population history when designing training samples.

Population-aware permutation-based significance thresholds for genome-wide association studies

Motivation

Permutation-based significance thresholds have been shown to be a robust alternative to classical Bonferroni significance thresholds in genome-wide association studies (GWAS) for skewed phenotype distributions. The recently published method permGWAS introduced a batch-wise approach to efficiently compute permutation-based GWAS. However, running multiple univariate tests in parallel leads to many repetitive computations and increased computational resources. More importantly, traditional permutation methods that permute only the phenotype break the underlying population structure.

Results

We propose permGWAS2, an improved method that does not break the population structure during permutations and uses an elegant block matrix decomposition to optimize computations, thereby reducing redundancies. We show on synthetic data that this improved approach yields a lower false discovery rate for skewed phenotype distributions compared to the previous version and the commonly used Bonferroni correction. In addition, we re-analyze a dataset covering phenotypic variation in 86 traits in a population of 615 wild sunflowers (Helianthus annuus L.). This led to the identification of dozens of novel associations with putatively adaptive traits, and removed several likely false-positive associations with limited biological support.

Availability and implementation

permGWAS2 is open-source and publicly available on GitHub for download: https://github.com/grimmlab/permGWAS.

Bivariate GWA mapping reveals associations between aliphatic glucosinolates and plant responses to thrips and heat stress

Although plants harbor a huge phytochemical diversity, only a fraction of plant metabolites is functionally characterized. In this work, we aimed to identify the genetic basis of metabolite functions during harsh environmental conditions in Arabidopsis thaliana. With machine learning algorithms we predicted stress-specific metabolomes for 23 (a)biotic stress phenotypes of 300 natural Arabidopsis accessions. The prediction models identified several aliphatic glucosinolates (GLSs) and their breakdown products to be implicated in responses to heat stress in siliques and herbivory by Western flower thrips, Frankliniella occidentalis. Bivariate GWA mapping of the metabolome predictions and their respective (a)biotic stress phenotype revealed genetic associations with MAM, AOP, and GS-OH, all three involved in aliphatic GSL biosynthesis. We, therefore, investigated thrips herbivory on AOP, MAM, and GS-OH loss-of-function and/or overexpression lines. Arabidopsis accessions with a combination of MAM2 and AOP3, leading to 3-hydroxypropyl dominance, suffered less from thrips feeding damage. The requirement of MAM2 for this effect could, however, not be confirmed with an introgression line of ecotypes Cvi and Ler, most likely due to other, unknown susceptibility factors in the Ler background. However, AOP2 and GS-OH, adding alkenyl or hydroxy-butenyl groups, respectively, did not have major effects on thrips feeding. Overall, this study illustrates the complex implications of aliphatic GSL diversity in plant responses to heat stress and a cell-content-feeding herbivore.

The benefits of permutation-based genome-wide association studies

Linear mixed models (LMMs) are a commonly used method for genome-wide association studies (GWAS) that aim to detect associations between genetic markers and phenotypic measurements in a population of individuals while accounting for population structure and cryptic relatedness. In a standard GWAS, hundreds of thousands to millions of statistical tests are performed, requiring control for multiple hypothesis testing. Typically, static corrections that penalize the number of tests performed are used to control for the family-wise error rate, which is the probability of making at least one false positive. However, it has been shown that in practice this threshold is too conservative for normally distributed phenotypes and not stringent enough for non-normally distributed phenotypes. Therefore, permutation-based LMM approaches have recently been proposed to provide a more realistic threshold that takes phenotypic distributions into account. In this work, we discuss the advantages of permutation-based GWAS approaches, including new simulations and results from a re-analysis of all publicly available Arabidopsis phenotypes from the AraPheno database.

Natural genetic variation underlying the negative effect of elevated CO2 on ionome composition in Arabidopsis thaliana

The elevation of atmospheric CO2 leads to a decline in plant mineral content, which might pose a significant threat to food security in coming decades. Although few genes have been identified for the negative effect of elevated CO2 on plant mineral composition, several studies suggest the existence of genetic factors. Here, we performed a large-scale study to explore genetic diversity of plant ionome responses to elevated CO2, using six hundred Arabidopsis thaliana accessions, representing geographical distributions ranging from worldwide to regional and local environments. We show that growth under elevated CO2 leads to a global decrease of ionome content, whatever the geographic distribution of the population. We observed a high range of genetic diversity, ranging from the most negative effect to resilience or even to a benefit in response to elevated CO2. Using genome-wide association mapping, we identified a large set of genes associated with this response, and we demonstrated that the function of one of these genes is involved in the negative effect of elevated CO2 on plant mineral composition. This resource will contribute to understand the mechanisms underlying the effect of elevated CO2 on plant mineral nutrition, and could help towards the development of crops adapted to a high-CO2 world.

Natural variation of warm temperature-induced raffinose accumulation identifies TREHALOSE-6-PHOSPHATE SYNTHASE 1 as a modulator of thermotolerance

High-temperature stress limits plant growth and reproduction. Exposure to high temperature, however, also elicits a physiological response, which protects plants from the damage evoked by heat. This response involves a partial reconfiguration of the metabolome including the accumulation of the trisaccharide raffinose. In this study, we explored the intraspecific variation of warm temperature-induced raffinose accumulation as a metabolic marker for temperature responsiveness with the aim to identify genes that contribute to thermotolerance. By combining raffinose measurements in 250 Arabidopsis thaliana accessions following a mild heat treatment with genome-wide association studies, we identified five genomic regions that were associated with the observed trait variation. Subsequent functional analyses confirmed a causal relationship between TREHALOSE-6-PHOSPHATE SYNTHASE 1 (TPS1) and warm temperature-dependent raffinose synthesis. Moreover, complementation of the tps1-1 null mutant with functionally distinct TPS1 isoforms differentially affected carbohydrate metabolism under more severe heat stress. While higher TPS1 activity was associated with reduced endogenous sucrose levels and thermotolerance, disruption of trehalose 6-phosphate signalling resulted in higher accumulation of transitory starch and sucrose and was associated with enhanced heat resistance. Taken together, our findings suggest a role of trehalose 6-phosphate in thermotolerance, most likely through its regulatory function in carbon partitioning and sucrose homoeostasis.

Genome-wide association study unveils ascorbate regulation by PAS/LOV PROTEIN during high light acclimation

Varying light conditions elicit metabolic responses as part of acclimation with changes in ascorbate levels being an important component. Here, we adopted a genome-wide association-based approach to characterize the response in ascorbate levels on high light (HL) acclimation in a panel of 315 Arabidopsis (Arabidopsis thaliana) accessions. These studies revealed statistically significant SNPs for total and reduced ascorbate under HL conditions at a locus in chromosome 2. Ascorbate levels under HL and the region upstream and within PAS/LOV PROTEIN (PLP) were strongly associated. Intriguingly, subcellular localization analyses revealed that the PLPA and PLPB splice variants co-localized with VITAMIN C DEFECTIVE2 (VTC2) and VTC5 in both the cytosol and nucleus. Yeast 2-hybrid and bimolecular fluorescence complementation analyses revealed that PLPA and PLPB interact with VTC2 and that blue light diminishes this interaction. Furthermore, PLPB knockout mutants were characterized by 1.5- to 1.7-fold elevations in their ascorbate levels, whereas knockout mutants of the cry2 cryptochromes displayed 1.2- to 1.3-fold elevations compared to WT. Our results collectively indicate that PLP plays a critical role in the elevation of ascorbate levels, which is a signature response of HL acclimation. The results strongly suggest that this is achieved via the release of the inhibitory effect of PLP on VTC2 upon blue light illumination, as the VTC2-PLPB interaction is stronger under darkness. The conditional importance of the cryptochrome receptors under different environmental conditions suggests a complex hierarchy underpinning the environmental control of ascorbate levels. However, the data we present here clearly demonstrate that PLP dominates during HL acclimation.

Somatic mutations inferred from RNA-seq data highlight the contribution of replication timing to mutation rate variation in a model plant

Variation in the rates and characteristics of germline and somatic mutations across the genome of an organism is informative about DNA damage and repair processes and can also shed light on aspects of organism physiology and evolution. We adapted a recently developed method for inferring somatic mutations from bulk RNA-seq data and applied it to a large collection of Arabidopsis thaliana accessions. The wide range of genomic data types available for A. thaliana enabled us to investigate the relationships of multiple genomic features with the variation in the somatic mutation rate across the genome of this model plant. We observed that late replicated regions showed evidence of an elevated rate of somatic mutation compared to genomic regions that are replicated early. We identified transcriptional strand asymmetries, consistent with the effects of transcription-coupled damage and/or repair. We also observed a negative relationship between the inferred somatic mutation count and the H3K36me3 histone mark which is well documented in the literature of human systems. In addition, we were able to support previous reports of an inverse relationship between inferred somatic mutation count and guanine-cytosine content as well as a positive relationship between inferred somatic mutation count and DNA methylation for both cytosine and noncytosine mutations.

Identification of gene function based on models capturing natural variability of Arabidopsis thaliana lipid metabolism

Lipids play fundamental roles in regulating agronomically important traits. Advances in plant lipid metabolism have until recently largely been based on reductionist approaches, although modulation of its components can have system-wide effects. However, existing models of plant lipid metabolism provide lumped representations, hindering detailed study of component modulation. Here, we present the Plant Lipid Module (PLM) which provides a mechanistic description of lipid metabolism in the Arabidopsis thaliana rosette. We demonstrate that the PLM can be readily integrated in models of A. thaliana Col-0 metabolism, yielding accurate predictions (83%) of single lethal knock-outs and 75% concordance between measured transcript and predicted flux changes under extended darkness. Genome-wide associations with fluxes obtained by integrating the PLM in diel condition- and accession-specific models identify up to 65 candidate genes modulating A. thaliana lipid metabolism. Using mutant lines, we validate up to 40% of the candidates, paving the way for identification of metabolic gene function based on models capturing natural variability in metabolism.

Forecasting trait responses in novel environments to aid seed provenancing under climate change

Revegetation projects face the major challenge of sourcing optimal plant material. This is often done with limited information about plant performance and increasingly requires factoring resilience to climate change. Functional traits can be used as quantitative indices of plant performance and guide seed provenancing, but trait values expected under novel conditions are often unknown. To support climate-resilient provenancing efforts, we develop a trait prediction model that integrates the effect of genetic variation with fine-scale temperature variation. We train our model on multiple field plantings of Arabidopsis thaliana and predict two relevant fitness traits-days-to-bolting and fecundity-across the species' European range. Prediction accuracy was high for days-to-bolting and moderate for fecundity, with the majority of trait variation explained by temperature differences between plantings. Projection under future climate predicted a decline in fecundity, although this response was heterogeneous across the range. In response, we identified novel genotypes that could be introduced to genetically offset the fitness decay. Our study highlights the value of predictive models to aid seed provenancing and improve the success of revegetation projects.

Analysis of dog breed diversity using a composite selection index

During breed development, domestic dogs have undergone genetic bottlenecks and sustained selective pressures, as a result distinctive genomic diversity occurs to varying degrees within and between breed groups. This diversity can be identified using standard methods or combinations of these methods. This study explored the application of a combined selection index, composite selection signals (CSS), derived from multiple methods to an existing genotype dataset from three breed groups developed in distinct regions of Asia: Qinghai-Tibet plateau dogs (adapted to living at altitude), Xi dogs (with superior running ability) and Mountain hounds (used for hunting ability). The CSS analysis confirmed top ranked genomic regions on CFA10 and CFA21 in Qinghai-Tibet plateau dogs, CFA1 in Xi dogs and CFA5 in Mountain hounds. CSS analysis identified additional significant genomic regions in each group, defined by a total of 1,397, 1,475 and 1,675 significant SNPs in the Qinghai-Tibetan Plateau dogs, Xi dogs and Mountain hounds, respectively. Chitinase 3 Like 1 (CHI3L1) and Leucine Rich Repeat Containing G Protein-Coupled Receptor 6 (LGR6) genes were located in the top ranked region on CFA7 (0.02-1 Mb) in the Qinghai-Tibetan Plateau dogs. Both genes have been associated with hypoxia responses or altitude adaptation in humans. For the Xi dogs, the top ranked region on CFA25 contained the Transient Receptor Potential Cation Channel Subfamily C Member 4 (TRPC4) gene. This calcium channel is important for optimal muscle performance during exercise. The outstanding signals in the Mountain dogs were on CFA5 with 213 significant SNPs that spanned genes involved in cardiac development, sight and generation of biochemical energy. These findings support the use of the combined index approach for identifying novel regions of genome diversity in dogs. As with other methods, the results do not prove causal links between these regions and phenotypes, but they may assist in focusing future studies that seek to identify functional pathways that contribute to breed diversity.

A genome-wide association study identifies novel players in Na and Fe homeostasis in Arabidopsis thaliana under alkaline-salinity stress

In nature, multiple stress factors occur simultaneously. The screening of natural diversity panels and subsequent Genome-Wide Association Studies (GWAS) is a powerful approach to identify genetic components of various stress responses. Here, the nutritional status variation of a set of 270 natural accessions of Arabidopsis thaliana grown on a natural saline-carbonated soil is evaluated. We report significant natural variation on leaf Na (LNa) and Fe (LFe) concentrations in the studied accessions. Allelic variation in the NINJA and YUC8 genes is associated with LNa diversity, and variation in the ALA3 is associated with LFe diversity. The allelic variation detected in these three genes leads to changes in their mRNA expression and correlates with plant differential growth performance when plants are exposed to alkaline salinity treatment under hydroponic conditions. We propose that YUC8 and NINJA expression patters regulate auxin and jasmonic signaling pathways affecting plant tolerance to alkaline salinity. Finally, we describe an impairment in growth and leaf Fe acquisition associated with differences in root expression of ALA3, encoding a phospholipid translocase active in plasma membrane and the trans Golgi network which directly interacts with proteins essential for the trafficking of PIN auxin transporters, reinforcing the role of phytohormonal processes in regulating ion homeostasis under alkaline salinity.

Predicting Gene Regulatory Interactions Using Natural Genetic Variation

Genome-wide association studies (GWAS) are a powerful tool to elucidate the genotype-phenotype map. Although GWAS are usually used to assess simple univariate associations between genetic markers and traits of interest, it is also possible to infer the underlying genetic architecture and to predict gene regulatory interactions. In this chapter, we describe the latest methods and tools to perform GWAS by calculating permutation-based significance thresholds. For this purpose, we first provide guidelines on univariate GWAS analyses that are extended in the second part of this chapter to more complex models that enable the inference of gene regulatory networks and how these networks vary.

Imputation of low-density marker chip data in plant breeding: Evaluation of methods based on sugar beet

Low-density genotyping followed by imputation reduces genotyping costs while still providing high-density marker information. An increased marker density has the potential to improve the outcome of all applications that are based on genomic data. This study investigates techniques for 1k to 20k genomic marker imputation for plant breeding programs with sugar beet (Beta vulgaris L. ssp. vulgaris) as an example crop, where these are realistic marker numbers for modern breeding applications. The generally accepted 'gold standard' for imputation, Beagle 5.1, was compared with the recently developed software AlphaPlantImpute2 which is designed specifically for plant breeding. For Beagle 5.1 and AlphaPlantImpute2, the imputation strategy as well as the imputation parameters were optimized in this study. We found that the imputation accuracy of Beagle could be tremendously improved (0.22 to 0.67) by tuning parameters, mainly by lowering the values for the parameter for the effective population size and increasing the number of iterations performed. Separating the phasing and imputation steps also improved accuracies when optimized parameters were used (0.67 to 0.82). We also found that the imputation accuracy of Beagle decreased when more low-density lines were included for imputation. AlphaPlantImpute2 produced very high accuracies without optimization (0.89) and was generally less responsive to optimization. Overall, AlphaPlantImpute2 performed relatively better for imputation whereas Beagle was better for phasing. Combining both tools yielded the highest accuracies.

Combined Approach of GWAS and Phylogenetic Analyses to Identify New Candidate Genes That Participate in Arabidopsis thaliana Primary Root Development Using Cellular Measurements and Primary Root Length

Genome-wide association studies (GWAS) have allowed the identification of different loci associated with primary root (PR) growth, and Arabidopsis is an excellent model for these studies. The PR length is controlled by cell proliferation, elongation, and differentiation; however, the specific contribution of proliferation and differentiation in the control of PR growth is still poorly studied. To this end, we analyzed 124 accessions and used a GWAS approach to identify potential causal genomic regions related to four traits: PR length, growth rate, cell proliferation and cell differentiation. Twenty-three genes and five statistically significant SNPs were identified. The SNP with the highest score mapped to the fifth exon of NAC048 and this change makes a missense variant in only 33.3% of the accessions with a large PR, compared with the accessions with a short PR length. Moreover, we detected five more SNPs in this gene and in NAC3 that allow us to discover closely related accessions according to the phylogenetic tree analysis. We also found that the association between genetic variants among the 18 genes with the highest scores in our GWAS and the phenotypic classes into which we divided our accessions are not straightforward and likely follow historical patterns.

A comparison of classical and machine learning-based phenotype prediction methods on simulated data and three plant species

Genomic selection is an integral tool for breeders to accurately select plants directly from genotype data leading to faster and more resource-efficient breeding programs. Several prediction methods have been established in the last few years. These range from classical linear mixed models to complex non-linear machine learning approaches, such as Support Vector Regression, and modern deep learning-based architectures. Many of these methods have been extensively evaluated on different crop species with varying outcomes. In this work, our aim is to systematically compare 12 different phenotype prediction models, including basic genomic selection methods to more advanced deep learning-based techniques. More importantly, we assess the performance of these models on simulated phenotype data as well as on real-world data from Arabidopsis thaliana and two breeding datasets from soy and corn. The synthetic phenotypic data allow us to analyze all prediction models and especially the selected markers under controlled and predefined settings. We show that Bayes B and linear regression models with sparsity constraints perform best under different simulation settings with respect to explained variance. Further, we can confirm results from other studies that there is no superiority of more complex neural network-based architectures for phenotype prediction compared to well-established methods. However, on real-world data, for which several prediction models yield comparable results with slight advantages for Elastic Net, this picture is less clear, suggesting that there is a lot of room for future research.

Efficient permutation-based genome-wide association studies for normal and skewed phenotypic distributions

MOTIVATION

Genome-wide association studies (GWAS) are an integral tool for studying the architecture of complex genotype and phenotype relationships. Linear mixed models (LMMs) are commonly used to detect associations between genetic markers and a trait of interest, while at the same time allowing to account for population structure and cryptic relatedness. Assumptions of LMMs include a normal distribution of the residuals and that the genetic markers are independent and identically distributed-both assumptions are often violated in real data. Permutation-based methods can help to overcome some of these limitations and provide more realistic thresholds for the discovery of true associations. Still, in practice, they are rarely implemented due to the high computational complexity.

RESULTS

We propose permGWAS, an efficient LMM reformulation based on 4D tensors that can provide permutation-based significance thresholds. We show that our method outperforms current state-of-the-art LMMs with respect to runtime and that permutation-based thresholds have lower false discovery rates for skewed phenotypes compared to the commonly used Bonferroni threshold. Furthermore, using permGWAS we re-analyzed more than 500 Arabidopsis thaliana phenotypes with 100 permutations each in less than 8 days on a single GPU. Our re-analyses suggest that applying a permutation-based threshold can improve and refine the interpretation of GWAS results.

AVAILABILITY AND IMPLEMENTATION

permGWAS is open-source and publicly available on GitHub for download: https://github.com/grimmlab/permGWAS.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Genome-wide association mapping within a local Arabidopsis thaliana population more fully reveals the genetic architecture for defensive metabolite diversity

A paradoxical finding from genome-wide association studies (GWAS) in plants is that variation in metabolite profiles typically maps to a small number of loci, despite the complexity of underlying biosynthetic pathways. This discrepancy may partially arise from limitations presented by geographically diverse mapping panels. Properties of metabolic pathways that impede GWAS by diluting the additive effect of a causal variant, such as allelic and genetic heterogeneity and epistasis, would be expected to increase in severity with the geographical range of the mapping panel. We hypothesized that a population from a single locality would reveal an expanded set of associated loci. We tested this in a French Arabidopsis thaliana population (less than 1 km transect) by profiling and conducting GWAS for glucosinolates, a suite of defensive metabolites that have been studied in depth through functional and genetic mapping approaches. For two distinct classes of glucosinolates, we discovered more associations at biosynthetic loci than the previous GWAS with continental-scale mapping panels. Candidate genes underlying novel associations were supported by concordance between their observed effects in the TOU-A population and previous functional genetic and biochemical characterization. Local populations complement geographically diverse mapping panels to reveal a more complete genetic architecture for metabolic traits. This article is part of the theme issue 'Genetic basis of adaptation and speciation: from loci to causative mutations'.

The megabase-scale crossover landscape is largely independent of sequence divergence

Meiotic recombination frequency varies along chromosomes and strongly correlates with sequence divergence. However, the causal relationship between recombination landscapes and polymorphisms is unclear. Here, we characterize the genome-wide recombination landscape in the quasi-absence of polymorphisms, using Arabidopsis thaliana homozygous inbred lines in which a few hundred genetic markers were introduced through mutagenesis. We find that megabase-scale recombination landscapes in inbred lines are strikingly similar to the recombination landscapes in hybrids, with the notable exception of heterozygous large rearrangements where recombination is prevented locally. In addition, the megabase-scale recombination landscape can be largely explained by chromatin features. Our results show that polymorphisms are not a major determinant of the shape of the megabase-scale recombination landscape but rather favour alternative models in which recombination and chromatin shape sequence divergence across the genome.

The frequency of recombination varies along chromosomes and highly correlates with sequence divergence. Here, the authors show that polymorphisms are not a major determinant of the megabase-scale recombination landscape in Arabidopsis, which is rather determined by chromatin accessibility and DNA methylation.

Adaptive significance of flowering time variation across natural seasonal environments in Arabidopsis thaliana

The relevance of flowering time variation and plasticity to climate adaptation requires a comprehensive empirical assessment. We investigated natural selection and the genetic architecture of flowering time in Arabidopsis through field experiments in Europe across multiple sites and seasons. We estimated selection for flowering time, plasticity and canalization. Loci associated with flowering time, plasticity and canalization by genome-wide association studies were tested for a geographic signature of climate adaptation. Selection favored early flowering and increased canalization, except at the northernmost site, but was rarely detected for plasticity. Genome-wide association studies revealed significant associations with flowering traits and supported a substantial polygenic inheritance. Alleles associated with late flowering, including functional FRIGIDA variants, were more common in regions experiencing high annual temperature variation. Flowering time plasticity to fall vs spring and summer environments was associated with GIGANTEA SUPPRESSOR 5, which promotes early flowering under decreasing day length and temperature. The finding that late flowering genotypes and alleles are associated with climate is evidence for past adaptation. Real-time phenotypic selection analysis, however, reveals pervasive contemporary selection for rapid flowering in agricultural settings across most of the species range. The response to this selection may involve genetic shifts in environmental cuing compared to the ancestral state.

Development of a versatile resource for post-genomic research through consolidating and characterizing 1500 diverse wild and cultivated soybean genomes

BACKGROUND

With advances in next-generation sequencing technologies, an unprecedented amount of soybean accessions has been sequenced by many individual studies and made available as raw sequencing reads for post-genomic research.

RESULTS

To develop a consolidated and user-friendly genomic resource for post-genomic research, we consolidated the raw resequencing data of 1465 soybean genomes available in the public and 91 highly diverse wild soybean genomes newly sequenced. These altogether provided a collection of 1556 sequenced genomes of 1501 diverse accessions (1.5 K). The collection comprises of wild, landraces and elite cultivars of soybean that were grown in East Asia or major soybean cultivating areas around the world. Our extensive sequence analysis discovered 32 million single nucleotide polymorphisms (32mSNPs) and revealed a SNP density of 30 SNPs/kb and 12 non-synonymous SNPs/gene reflecting a high structural and functional genomic diversity of the new collection. Each SNP was annotated with 30 categories of structural and/or functional information. We further identified paired accessions between the 1.5 K and 20,087 (20 K) accessions in US collection as genomic "equivalent" accessions sharing the highest genomic identity for minimizing the barriers in soybean germplasm exchange between countries. We also exemplified the utility of 32mSNPs in enhancing post-genomics research through in-silico genotyping, high-resolution GWAS, discovering and/or characterizing genes and alleles/mutations, identifying germplasms containing beneficial alleles that are potentially experiencing artificial selection.

CONCLUSION

The comprehensive analysis of publicly available large-scale genome sequencing data of diverse cultivated accessions and the newly in-house sequenced wild accessions greatly increased the soybean genome-wide variation resolution. This could facilitate a variety of genetic and molecular-level analyses in soybean. The 32mSNPs and 1.5 K accessions with their comprehensive annotation have been made available at the SoyBase and Ag Data Commons. The dataset could further serve as a versatile and expandable core resource for exploring the exponentially increasing genome sequencing data for a variety of post-genomic research.

Genome-wide association of the metabolic shifts underpinning dark-induced senescence in Arabidopsis

Dark-induced senescence provokes profound metabolic shifts to recycle nutrients and to guarantee plant survival. To date, research on these processes has largely focused on characterizing mutants deficient in individual pathways. Here, we adopted a time-resolved genome-wide association-based approach to characterize dark-induced senescence by evaluating the photochemical efficiency and content of primary and lipid metabolites at the beginning, or after 3 or 6 days in darkness. We discovered six patterns of metabolic shifts and identified 215 associations with 81 candidate genes being involved in this process. Among these associations, we validated the roles of four genes associated with glycine, galactinol, threonine, and ornithine levels. We also demonstrated the function of threonine and galactinol catabolism during dark-induced senescence. Intriguingly, we determined that the association between tyrosine contents and TYROSINE AMINOTRANSFERASE 1 influences enzyme activity of the encoded protein and transcriptional activity of the gene under normal and dark conditions, respectively. Moreover, the single-nucleotide polymorphisms affecting the expression of THREONINE ALDOLASE 1 and the amino acid transporter gene AVT1B, respectively, only underlie the variation in threonine and glycine levels in the dark. Taken together, these results allow us to present a very detailed model of the metabolic aspects of dark-induced senescence, as well as the process itself.

Genome wide association study reveals plant loci controlling heritability of the rhizosphere microbiome

Host genetics has recently been shown to be a driver of plant microbiome composition. However, identifying the underlying genetic loci controlling microbial selection remains challenging. Genome-wide association studies (GWAS) represent a potentially powerful, unbiased method to identify microbes sensitive to the host genotype and to connect them with the genetic loci that influence their colonization. Here, we conducted a population-level microbiome analysis of the rhizospheres of 200 sorghum genotypes. Using 16S rRNA amplicon sequencing, we identify rhizosphere-associated bacteria exhibiting heritable associations with plant genotype, and identify significant overlap between these lineages and heritable taxa recently identified in maize. Furthermore, we demonstrate that GWAS can identify host loci that correlate with the abundance of specific subsets of the rhizosphere microbiome. Finally, we demonstrate that these results can be used to predict rhizosphere microbiome structure for an independent panel of sorghum genotypes based solely on knowledge of host genotypic information.

Improving Genomic Prediction Using High-Dimensional Secondary Phenotypes

In the past decades, genomic prediction has had a large impact on plant breeding. Given the current advances of high-throughput phenotyping and sequencing technologies, it is increasingly common to observe a large number of traits, in addition to the target trait of interest. This raises the important question whether these additional or “secondary” traits can be used to improve genomic prediction for the target trait. With only a small number of secondary traits, this is known to be the case, given sufficiently high heritabilities and genetic correlations. Here we focus on the more challenging situation with a large number of secondary traits, which is increasingly common since the arrival of high-throughput phenotyping. In this case, secondary traits are usually incorporated through additional relatedness matrices. This approach is however infeasible when secondary traits are not measured on the test set, and cannot distinguish between genetic and non-genetic correlations. An alternative direction is to extend the classical selection indices using penalized regression. So far, penalized selection indices have not been applied in a genomic prediction setting, and require plot-level data in order to reliably estimate genetic correlations. Here we aim to overcome these limitations, using two novel approaches. Our first approach relies on a dimension reduction of the secondary traits, using either penalized regression or random forests (LS-BLUP/RF-BLUP). We then compute the bivariate GBLUP with the dimension reduction as secondary trait. For simulated data (with available plot-level data), we also use bivariate GBLUP with the penalized selection index as secondary trait (SI-BLUP). In our second approach (GM-BLUP), we follow existing multi-kernel methods but replace secondary traits by their genomic predictions, with the advantage that genomic prediction is also possible when secondary traits are only measured on the training set. For most of our simulated data, SI-BLUP was most accurate, often closely followed by RF-BLUP or LS-BLUP. In real datasets, involving metabolites in Arabidopsis and transcriptomics in maize, no method could substantially improve over univariate prediction when secondary traits were only available on the training set. LS-BLUP and RF-BLUP were most accurate when secondary traits were available also for the test set.

Plant-ImputeDB: an integrated multiple plant reference panel database for genotype imputation

Genotype imputation is a process that estimates missing genotypes in terms of the haplotypes and genotypes in a reference panel. It can effectively increase the density of single nucleotide polymorphisms (SNPs), boost the power to identify genetic association and promote the combination of genetic studies. However, there has been a lack of high-quality reference panels for most plants, which greatly hinders the application of genotype imputation. Here, we developed Plant-ImputeDB (http://gong_lab.hzau.edu.cn/Plant_imputeDB/), a comprehensive database with reference panels of 12 plant species for online genotype imputation, SNP and block search and free download. By integrating genotype data and whole-genome resequencing data of plants from various studies and databases, the current Plant-ImputeDB provides high-quality reference panels of 12 plant species, including ∼69.9 million SNPs from 34 244 samples. It also provides an easy-to-use online tool with the option of two popular tools specifically designed for genotype imputation. In addition, Plant-ImputeDB accepts submissions of different types of genomic variations, and provides free and open access to all publicly available data in support of related research worldwide. In general, Plant-ImputeDB may serve as an important resource for plant genotype imputation and greatly facilitate the research on plant genetic research.

The patterns of deleterious mutations during the domestication of soybean

Globally, soybean is a major protein and oil crop. Enhancing our understanding of the soybean domestication and improvement process helps boost genomics-assisted breeding efforts. Here we present a genome-wide variation map of 10.6 million single-nucleotide polymorphisms and 1.4 million indels for 781 soybean individuals which includes 418 domesticated (Glycine max), 345 wild (Glycine soja), and 18 natural hybrid (G. max/G. soja) accessions. We describe the enhanced detection of 183 domestication-selective sweeps and the patterns of putative deleterious mutations during domestication and improvement. This predominantly selfing species shows 7.1% reduction of overall deleterious mutations in domesticated soybean relative to wild soybean and a further 1.4% reduction from landrace to improved accessions. The detected domestication-selective sweeps also show reduced levels of deleterious alleles. Importantly, genotype imputation with this resource increases the mapping resolution of genome-wide association studies for seed protein and oil traits in a soybean diversity panel.