Genomic selection performs as effectively as phenotypic selection for increasing seed yield in soybean.

Increasing the rate of genetic gain for seed yield remains the primary breeding objective in both public and private soybean [Glycine max (L.) Merr.] breeding programs. Genomic selection (GS) has the potential to accelerate the rate of genetic gain for soybean seed yield. Limited studies to date have validated GS accuracy and directly compared GS with phenotypic selection (PS), and none have been reported in soybean. This study conducted the first empirical validation of GS for increasing seed yield using over 1,500 lines and over 7 yr (2010-2016) of replicated experiments in the University of Nebraska-Lincoln soybean breeding program. The study was designed to capture the varying genetic relatedness of the training population to three validation sets: two large biparental populations (TBP-1 and TBP-2) and a large validation set comprised of 457 preselected advanced lines derived from 45 biparental populations (TMP). We found that prediction accuracy (.54) realized in our validation experiments was comparable with what we obtained from a series of cross-validation experiments (.64). Both GS and PS were more effective for increasing the population mean performance compared with random selection (RS). We found a selection advantage of GS over PS, where higher genetic gain and identification of top-performing lines was maximized at 10-20% selected proportion. Genomic selection led to small increases in genetic similarity when compared with PS and RS presumably because of a significant shift on allelic frequencies toward the extremes, suggesting that it could erode genetic diversity more quickly. Overall, we found that GS can perform as effectively as PS but that measures should be considered to protect against loss of genetic variance when using GS.

Multi-trait genomic prediction improves selection accuracy for enhancing seed mineral concentrations in pea.

Multi-trait genomic selection (MT-GS) has the potential to improve predictive ability by maximizing the use of information across related genotypes and genetically correlated traits. In this study, we extended the use of sparse phenotyping method into the MT-GS framework by split testing of entries to maximize borrowing of information across genotypes and predict missing phenotypes for targeted traits without additional phenotyping expenditure. Using 300 advanced breeding lines from North Dakota State University (NDSU) pulse breeding program and ∼200 USDA accessions that were evaluated for 10 nutritional traits, our results show that the proposed sparse phenotyping aided MT-GS can further improve predictive ability by >12% across traits compared with univariate (UNI) genomic selection. The proposed strategy departed from the previous reports that weak genetic correlation is a limitation to the advantage of MT-GS over UNI genomic selection, which was evident in the partially balanced phenotyping-enabled MT-GS. Our results point to heritability and genetic correlation between traits as possible metrics to optimize and further improve the estimation of model parameters, and ultimately, prediction performance. Overall, our study offers a new approach to optimize the prediction performance using the MT-GS and further highlight strategy to maximize the efficiency of GS in a plant breeding program. The sparse-testing-aided MT-GS proposed in this study can be further extended to multi-environment, multi-trait GS to improve prediction performance and further reduce the cost of phenotyping and time-consuming data collection process.

We extended the use of sparse phenotyping into the multi-trait genomic selection (MT-GS) framework by split testing of entries. The sparse-phenotyping-aided MT-GS can increase predictive ability by >12% across traits. Heritability and genetic correlation are possible metrics to optimize and further improve prediction performance of MT-GS. The sparse-testing-aided MT-GS can be further extended to multi-environment, multi-trait GS framework.

Genetic mapping and genomic prediction of sclerotinia stem rot resistance to rapeseed/canola (Brassica napus L.) at seedling stage.

KEY MESSAGE: GWAS detected ninety-eight significant SNPs associated with Sclerotinia sclerotiorum resistance. Six statistical models resulted in medium to high predictive ability, depending on trait, indicating potential of genomic prediction for disease resistance breeding. The lack of complete host resistance and a complex resistance inheritance nature between rapeseed/canola and Sclerotinia sclerotiorum often limits the development of functional molecular markers that enable breeding for sclerotinia stem rot (SSR) resistance. However, genomics-assisted selection has the potential to accelerate the breeding for SSR resistance. Therefore, genome-wide association (GWA) mapping and genomic prediction (GP) were performed using a diverse panel of 337 rapeseed/canola genotypes. Three-week-old seedlings were screened using the petiole inoculation technique (PIT). Days to wilt (DW) up to 2 weeks and lesion phenotypes (LP) at 3, 4, and 7 days post-inoculation (dpi) were recorded. A strong correlation (r = - 0.90) between DW and LP_4dpi implied that a single time point scoring at four days could be used as a proxy trait. GWA analyses using single-locus (SL) and multi-locus (ML) models identified a total of 41, and 208 significantly associated SNPs, respectively. Out of these, ninety-eight SNPs were identified by a combination of the SL model and any of the ML models, at least two ML models, or two traits. These SNPs explained 1.25-12.22% of the phenotypic variance and considered as significant, could be associated with SSR resistance. Eighty-three candidate genes with a function in disease resistance were associated with the significant SNPs. Six GP models resulted in moderate to high (0.42-0.67) predictive ability depending on SSR resistance traits. The resistant genotypes and significant SNPs will serve as valuable resources for future SSR resistance breeding. Our results also highlight the potential of genomic selection to improve rapeseed/canola breeding for SSR resistance.

New genomic regions associated with white mold resistance in dry bean using a MAGIC population.

Dry bean (Phaseolus vulgaris L.) production in many regions is threatened by white mold (WM) [Sclerotinia sclerotiorum (Lib.) de Bary]. Seed yield losses can be up to 100% under conditions favorable for the pathogen. The low heritability, polygenic inheritance, and cumbersome screening protocols make it difficult to breed for improved genetic resistance. Some progress in understanding genetic resistance and germplasm improvement has been accomplished, but cultivars with high levels of resistance are yet to be released. A WM multiparent advanced generation inter-cross (MAGIC) population (n = 1060) was developed to facilitate mapping and breeding efforts. A seedling straw test screening method provided a quick assay to phenotype the population for response to WM isolate 1980. Nineteen MAGIC lines were identified with improved resistance. For genome-wide association studies (GWAS), the data was transformed into three phenotypic distributions-quantitative, polynomial, and binomial-and coupled with ∼52,000 single-nucleotide polymorphisms (SNPs). The three phenotypic distributions identified 30 significant genomic intervals [-log10 (P value) ≥ 3.0]. However, across distributions, four new genomic regions as well as two regions previously reported were found to be associated with resistance. Cumulative R2 values were 57% for binomial distribution using 13 genomic intervals, 41% for polynomial using eight intervals, and 40% for quantitative using 11 intervals. New resistant germplasm as well as new genomic regions associated with resistance are now available for further investigation.

Phenotyping stomatal closure by thermal imaging for GWAS and TWAS of water use efficiency-related genes.

Stomata allow CO2 uptake by leaves for photosynthetic assimilation at the cost of water vapor loss to the atmosphere. The opening and closing of stomata in response to fluctuations in light intensity regulate CO2 and water fluxes and are essential for maintaining water-use efficiency (WUE). However, a little is known about the genetic basis for natural variation in stomatal movement, especially in C4 crops. This is partly because the stomatal response to a change in light intensity is difficult to measure at the scale required for association studies. Here, we used high-throughput thermal imaging to bypass the phenotyping bottleneck and assess 10 traits describing stomatal conductance (gs) before, during and after a stepwise decrease in light intensity for a diversity panel of 659 sorghum (Sorghum bicolor) accessions. Results from thermal imaging significantly correlated with photosynthetic gas exchange measurements. gs traits varied substantially across the population and were moderately heritable (h2 up to 0.72). An integrated genome-wide and transcriptome-wide association study identified candidate genes putatively driving variation in stomatal conductance traits. Of the 239 unique candidate genes identified with the greatest confidence, 77 were putative orthologs of Arabidopsis (Arabidopsis thaliana) genes related to functions implicated in WUE, including stomatal opening/closing (24 genes), stomatal/epidermal cell development (35 genes), leaf/vasculature development (12 genes), or chlorophyll metabolism/photosynthesis (8 genes). These findings demonstrate an approach to finding genotype-to-phenotype relationships for a challenging trait as well as candidate genes for further investigation of the genetic basis of WUE in a model C4 grass for bioenergy, food, and forage production.

Finding Needles in a Haystack: Using Geo-References to Enhance the Selection and Utilization of Landraces in Breeding for Climate-Resilient Cultivars of Upland Cotton (Gossypium hirsutum L.).

The genetic uniformity of cultivated cotton as a consequence of domestication and modern breeding makes it extremely vulnerable to abiotic challenges brought about by major climate shifts. To sustain productivity amidst worsening agro-environments, future breeding objectives need to seriously consider introducing new genetic variation from diverse resources into the current germplasm base of cotton. Landraces are genetically heterogeneous, population complexes that have been primarily selected for their adaptability to specific localized or regional environments. This makes them an invaluable genetic resource of novel allelic diversity that can be exploited to enhance the resilience of crops to marginal environments. The utilization of cotton landraces in breeding programs are constrained by the phenology of the plant and the lack of phenotypic information that can facilitate efficient selection of potential donor parents for breeding. In this review, the genetic value of cotton landraces and the major challenges in their utilization in breeding are discussed. Two strategies namely Focused Identification of Germplasm Strategy and Environmental Association Analysis that have been developed to effectively screen large germplasm collections for accessions with adaptive traits using geo-reference-based, mathematical modelling are highlighted. The potential applications of both approaches in mining available cotton landrace collections are also presented.

Comparative evolutionary genetics of deleterious load in sorghum and maize.

Sorghum and maize share a close evolutionary history that can be explored through comparative genomics1,2. To perform a large-scale comparison of the genomic variation between these two species, we analysed ~13 million variants identified from whole-genome resequencing of 499 sorghum lines together with 25 million variants previously identified in 1,218 maize lines. Deleterious mutations in both species were prevalent in pericentromeric regions, enriched in non-syntenic genes and present at low allele frequencies. A comparison of deleterious burden between sorghum and maize revealed that sorghum, in contrast to maize, departed from the domestication-cost hypothesis that predicts a higher deleterious burden among domesticates compared with wild lines. Additionally, sorghum and maize population genetic summary statistics were used to predict a gene deleterious index with an accuracy greater than 0.5. This research represents a key step towards understanding the evolutionary dynamics of deleterious variants in sorghum and provides a comparative genomics framework to start prioritizing these variants for removal through genome editing and breeding.

Harnessing Genetic Diversity in the USDA Pea Germplasm Collection Through Genomic Prediction.

Phenotypic evaluation and efficient utilization of germplasm collections can be time-intensive, laborious, and expensive. However, with the plummeting costs of next-generation sequencing and the addition of genomic selection to the plant breeder's toolbox, we now can more efficiently tap the genetic diversity within large germplasm collections. In this study, we applied and evaluated genomic prediction's potential to a set of 482 pea (Pisum sativum L.) accessions-genotyped with 30,600 single nucleotide polymorphic (SNP) markers and phenotyped for seed yield and yield-related components-for enhancing selection of accessions from the USDA Pea Germplasm Collection. Genomic prediction models and several factors affecting predictive ability were evaluated in a series of cross-validation schemes across complex traits. Different genomic prediction models gave similar results, with predictive ability across traits ranging from 0.23 to 0.60, with no model working best across all traits. Increasing the training population size improved the predictive ability of most traits, including seed yield. Predictive abilities increased and reached a plateau with increasing number of markers presumably due to extensive linkage disequilibrium in the pea genome. Accounting for population structure effects did not significantly boost predictive ability, but we observed a slight improvement in seed yield. By applying the best genomic prediction model (e.g., RR-BLUP), we then examined the distribution of genotyped but nonphenotyped accessions and the reliability of genomic estimated breeding values (GEBV). The distribution of GEBV suggested that none of the nonphenotyped accessions were expected to perform outside the range of the phenotyped accessions. Desirable breeding values with higher reliability can be used to identify and screen favorable germplasm accessions. Expanding the training set and incorporating additional orthogonal information (e.g., transcriptomics, metabolomics, physiological traits, etc.) into the genomic prediction framework can enhance prediction accuracy.

Transcriptome-Wide Association Supplements Genome-Wide Association in Zea mays.

Modern improvement of complex traits in agricultural species relies on successful associations of heritable molecular variation with observable phenotypes. Historically, this pursuit has primarily been based on easily measurable genetic markers. The recent advent of new technologies allows assaying and quantifying biological intermediates (hereafter endophenotypes) which are now readily measurable at a large scale across diverse individuals. The usefulness of endophenotypes for delineating the regulatory landscape of the genome and genetic dissection of complex trait variation remains underexplored in plants. The work presented here illustrated the utility of a large-scale (299-genotype and seven-tissue) gene expression resource to dissect traits across multiple levels of biological organization. Using single-tissue- and multi-tissue-based transcriptome-wide association studies (TWAS), we revealed that about half of the functional variation acts through altered transcript abundance for maize kernel traits, including 30 grain carotenoid abundance traits, 20 grain tocochromanol abundance traits, and 22 field-measured agronomic traits. Comparing the efficacy of TWAS with genome-wide association studies (GWAS) and an ensemble approach that combines both GWAS and TWAS, we demonstrated that results of TWAS in combination with GWAS increase the power to detect known genes and aid in prioritizing likely causal genes. Using a variance partitioning approach in the largely independent maize Nested Association Mapping (NAM) population, we also showed that the most strongly associated genes identified by combining GWAS and TWAS explain more heritable variance for a majority of traits than the heritability captured by the random genes and the genes identified by GWAS or TWAS alone. This not only improves the ability to link genes to phenotypes, but also highlights the phenotypic consequences of regulatory variation in plants.

Metabolome-Scale Genome-Wide Association Studies Reveal Chemical Diversity and Genetic Control of Maize Specialized Metabolites.

Cultivated maize (Zea mays) has retained much of the genetic diversity of its wild ancestors. Here, we performed nontargeted liquid chromatography-mass spectrometry metabolomics to analyze the metabolomes of the 282 maize inbred lines in the Goodman Diversity Panel. This analysis identified a bimodal distribution of foliar metabolites. Although 15% of the detected mass features were present in >90% of the inbred lines, the majority were found in <50% of the samples. Whereas leaf bases and tips were differentiated by flavonoid abundance, maize varieties (stiff-stalk, nonstiff-stalk, tropical, sweet maize, and popcorn) showed differential accumulation of benzoxazinoid metabolites. Genome-wide association studies (GWAS), performed for 3,991 mass features from the leaf tips and leaf bases, showed that 90% have multiple significantly associated loci scattered across the genome. Several quantitative trait locus hotspots in the maize genome regulate the abundance of multiple, often structurally related mass features. The utility of maize metabolite GWAS was demonstrated by confirming known benzoxazinoid biosynthesis genes, as well as by mapping isomeric variation in the accumulation of phenylpropanoid hydroxycitric acid esters to a single linkage block in a citrate synthase-like gene. Similar to gene expression databases, this metabolomic GWAS data set constitutes an important public resource for linking maize metabolites with biosynthetic and regulatory genes.

Deleterious Mutation Burden and Its Association with Complex Traits in Sorghum (Sorghum bicolor).

Sorghum (Sorghum bicolor L.) is a major food cereal for millions of people worldwide. The sorghum genome, like other species, accumulates deleterious mutations, likely impacting its fitness. The lack of recombination, drift, and the coupling with favorable loci impede the removal of deleterious mutations from the genome by selection. To study how deleterious variants impact phenotypes, we identified putative deleterious mutations among ∼5.5 M segregating variants of 229 diverse biomass sorghum lines. We provide the whole-genome estimate of the deleterious burden in sorghum, showing that ∼33% of nonsynonymous substitutions are putatively deleterious. The pattern of mutation burden varies appreciably among racial groups. Across racial groups, the mutation burden correlated negatively with biomass, plant height, specific leaf area (SLA), and tissue starch content (TSC), suggesting that deleterious burden decreases trait fitness. Putatively deleterious variants explain roughly one-half of the genetic variance. However, there is only moderate improvement in total heritable variance explained for biomass (7.6%) and plant height (average of 3.1% across all stages). There is no advantage in total heritable variance for SLA and TSC. The contribution of putatively deleterious variants to phenotypic diversity therefore appears to be dependent on the genetic architecture of traits. Overall, these results suggest that incorporating putatively deleterious variants into genomic models slightly improves prediction accuracy because of extensive linkage. Knowledge of deleterious variants could be leveraged for sorghum breeding through either genome editing and/or conventional breeding that focuses on the selection of progeny with fewer deleterious alleles.

Dissecting the Genetic Basis of Local Adaptation in Soybean.

Soybean (Glycine max) is the most widely grown oilseed in the world and is an important source of protein for both humans and livestock. Soybean is widely adapted to both temperate and tropical regions, but a changing climate demands a better understanding of adaptation to specific environmental conditions. Here, we explore genetic variation in a collection of 3,012 georeferenced, locally adapted landraces from a broad geographical range to help elucidate the genetic basis of local adaptation. We used geographic origin, environmental data and dense genome-wide SNP data to perform an environmental association analysis and discover loci displaying steep gradients in allele frequency across geographical distance and between landrace and modern cultivars. Our combined application of methods in environmental association mapping and detection of selection targets provide a better understanding of how geography and selection may have shaped genetic variation among soybean landraces. Moreover, we identified several important candidate genes related to drought and heat stress, and revealed important genomic regions possibly involved in the geographic divergence of soybean.

Genome-wide Association Mapping of Qualitatively Inherited Traits in a Germplasm Collection.

Genome-wide association (GWA) has been used as a tool for dissecting the genetic architecture of quantitatively inherited traits. We demonstrate here that GWA can also be highly useful for detecting many major genes governing categorically defined phenotype variants that exist for qualitatively inherited traits in a germplasm collection. Genome-wide association mapping was applied to categorical phenotypic data available for 10 descriptive traits in a collection of ∼13,000 soybean [ (L.) Merr.] accessions that had been genotyped with a 50,000 single nucleotide polymorphism (SNP) chip. A GWA on a panel of accessions of this magnitude can offer substantial statistical power and mapping resolution, and we found that GWA mapping resulted in the identification of strong SNP signals for 24 classical genes as well as several heretofore unknown genes controlling the phenotypic variants in those traits. Because some of these genes had been cloned, we were able to show that the narrow GWA mapping SNP signal regions that we detected for the phenotypic variants had chromosomal bp spans that, with just one exception, overlapped the bp region of the cloned genes, despite local variation in SNP number and nonuniform SNP distribution in the chip set.

Allelic variants of OsHKT1;1 underlie the divergence between indica and japonica subspecies of rice (Oryza sativa) for root sodium content.

Salinity is a major factor limiting crop productivity. Rice (Oryza sativa), a staple crop for the majority of the world, is highly sensitive to salinity stress. To discover novel sources of genetic variation for salt tolerance-related traits in rice, we screened 390 diverse accessions under 14 days of moderate (9 dS·m-1) salinity. In this study, shoot growth responses to moderate levels of salinity were independent of tissue Na+ content. A significant difference in root Na+ content was observed between the major subpopulations of rice, with indica accessions displaying higher root Na+ and japonica accessions exhibiting lower root Na+ content. The genetic basis of the observed variation in phenotypes was elucidated through genome-wide association (GWA). The strongest associations were identified for root Na+:K+ ratio and root Na+ content in a region spanning ~575 Kb on chromosome 4, named Root Na+ Content 4 (RNC4). Two Na+ transporters, HKT1;1 and HKT1;4 were identified as candidates for RNC4. Reduced expression of both HKT1;1 and HKT1;4 through RNA interference indicated that HKT1;1 regulates shoot and root Na+ content, and is likely the causal gene underlying RNC4. Three non-synonymous mutations within HKT1;1 were present at higher frequency in the indica subpopulation. When expressed in Xenopus oocytes the indica-predominant isoform exhibited higher inward (negative) currents and a less negative voltage threshold of inward rectifying current activation compared to the japonica-predominant isoform. The introduction of a 4.5kb fragment containing the HKT1;1 promoter and CDS from an indica variety into a japonica background, resulted in a phenotype similar to the indica subpopulation, with higher root Na+ and Na+:K+. This study provides evidence that HKT1;1 regulates root Na+ content, and underlies the divergence in root Na+ content between the two major subspecies in rice.

A Population Structure and Genome-Wide Association Analysis on the USDA Soybean Germplasm Collection.

Population structure analyses and genome-wide association studies (GWAS) conducted on crop germplasm collections provide valuable information on the frequency and distribution of alleles governing economically important traits. The value of these analyses is substantially enhanced when the accession numbers can be increased from ∼1,000 to ∼10,000 or more. In this research, we conducted the first comprehensive analysis of population structure on the collection of 14,000 soybean accessions [Glycine max (L.) Merr. and G. soja Siebold & Zucc.] using a 50K-SNP chip. Accessions originating from Japan were relatively homogenous and distinct from the Korean accessions. As a whole, both Japanese and Korean accessions diverged from the Chinese accessions. The ancestry of founders of the American accessions derived mostly from two Chinese subpopulations, which reflects the composition of the American accessions as a whole. A 12,000 accession GWAS conducted on seed protein and oil is the largest reported to date in plants and identified single nucleotide polymorphisms (SNPs) with strong signals on chromosomes 20 and 15. A chromosome 20 region previously reported to be important for protein and oil content was further narrowed and now contains only three plausible candidate genes. The haplotype effects show a strong negative relationship between oil and protein at this locus, indicating negative pleiotropic effects or multiple closely linked loci in repulsion phase linkage. The vast majority of accessions carry the haplotype allele conferring lower protein and higher oil. Our results provide a fuller understanding of the distribution of genetic variation contained within the USDA soybean collection and how it relates to phenotypic variation for economically important traits.

Multi-parent advanced generation inter-cross (MAGIC) populations in rice: progress and potential for genetics research and breeding.

BACKGROUND: This article describes the development of Multi-parent Advanced Generation Inter-Cross populations (MAGIC) in rice and discusses potential applications for mapping quantitative trait loci (QTLs) and for rice varietal development. We have developed 4 multi-parent populations: indica MAGIC (8 indica parents); MAGIC plus (8 indica parents with two additional rounds of 8-way F1 inter-crossing); japonica MAGIC (8 japonica parents); and Global MAGIC (16 parents - 8 indica and 8 japonica). The parents used in creating these populations are improved varieties with desirable traits for biotic and abiotic stress tolerance, yield, and grain quality. The purpose is to fine map QTLs for multiple traits and to directly and indirectly use the highly recombined lines in breeding programs. These MAGIC populations provide a useful germplasm resource with diverse allelic combinations to be exploited by the rice community. RESULTS: The indica MAGIC population is the most advanced of the MAGIC populations developed thus far and comprises 1328 lines produced by single seed descent (SSD). At the S4 stage of SSD a subset (200 lines) of this population was genotyped using a genotyping-by-sequencing (GBS) approach and was phenotyped for multiple traits, including: blast and bacterial blight resistance, salinity and submergence tolerance, and grain quality. Genome-wide association mapping identified several known major genes and QTLs including Sub1 associated with submergence tolerance and Xa4 and xa5 associated with resistance to bacterial blight. Moreover, the genome-wide association study (GWAS) results also identified potentially novel loci associated with essential traits for rice improvement. CONCLUSION: The MAGIC populations serve a dual purpose: permanent mapping populations for precise QTL mapping and for direct and indirect use in variety development. Unlike a set of naturally diverse germplasm, this population is tailor-made for breeders with a combination of useful traits derived from multiple elite breeding lines. The MAGIC populations also present opportunities for studying the interactions of genome introgressions and chromosomal recombination.