From association to prediction: statistical methods for the dissection and selection of complex traits in plants

Genomic prediction of cereal crop architectural traits using models informed by gene regulatory circuitries from maize

Plant architecture is a major determinant of planting density, which enhances productivity potential for crops per unit area. Genomic prediction is well positioned to expedite genetic gain of plant architectural traits since they are typically highly heritable. Additionally, the adaptation of genomic prediction models to query predictive abilities of markers tagging certain genomic regions could shed light on the genetic architecture of these traits. Here, we leveraged transcriptional networks from a prior study that contextually described developmental progression during tassel and leaf organogenesis in maize (Zea mays) to inform genomic prediction models for architectural traits. Since these developmental processes underlie tassel branching and leaf angle, 2 important agronomic architectural traits, we tested whether genes prioritized from these networks quantitatively contribute to the genetic architecture of these traits. We used genomic prediction models to evaluate the ability of markers in the vicinity of prioritized network genes to predict breeding values of tassel branching and leaf angle traits for 2 diversity panels in maize and diversity panels from sorghum (Sorghum bicolor) and rice (Oryza sativa). Predictive abilities of markers near these prioritized network genes were similar to those using whole-genome marker sets. Notably, markers near highly connected transcription factors from core network motifs in maize yielded predictive abilities that were significantly greater than expected by chance in not only maize but also closely related sorghum. We expect that these highly connected regulators are key drivers of architectural variation that are conserved across closely related cereal crop species.

Genome-wide association study and genotypic variation for the major tocopherol content in rice grain

Rice tocopherols, vitamin E compounds with antioxidant activity, play essential roles in human health. Even though the key genes involved in vitamin E biosynthetic pathways have been identified in plants, the genetic architecture of vitamin E content in rice grain remains unclear. A genome-wide association study (GWAS) on 179 genotypically diverse rice accessions with 34,323 SNP markers was conducted to detect QTLs that define total and α- tocopherol contents in rice grains. Total and α-tocopherol contents had a strong positive correlation and varied greatly across the accessions, ranging from 0.230-31.76 and 0.011-30.83 (μg/g), respectively. A total of 13 QTLs were identified, which were spread across five of the rice chromosomes. Among the 13 QTLs, 11 were considered major with phenotypic variation explained (PVE) greater than 10%. Twelve transcription factor (TF) genes, one microprotein (miP), and a transposon were found to be associated with the QTLs with putative roles in controlling tocopherol contents. Moreover, intracellular transport proteins, ABC transporters, nonaspanins, and SNARE, were identified as associated genes on chromosomes 1 and 8. In the vicinity of seven QTLs, protein kinases were identified as key signaling factors. Haplotype analysis revealed the QTLs qAlph1.1, qTot1.1, qAlph2.1, qAlph6.1, qTot6.1, and qTot8.3 to have significant haplogroups. Quantitative RT-PCR validated the expression direction and magnitude of WRKY39 (Os02g0265200), PIP5Ks (Os08g0450800), and MADS59 (Os06g0347700) in defining the major tocopherol contents. This study provides insights for ongoing biofortification efforts to breed and/or engineer vitamin E and antioxidant levels in rice and other cereals.

Mapping of flumioxazin tolerance in a snap bean diversity panel leads to the discovery of a master genomic region controlling multiple stress resistance genes

Introduction

Effective weed management tools are crucial for maintaining the profitable production of snap bean (Phaseolus vulgaris L.). Preemergence herbicides help the crop to gain a size advantage over the weeds, but the few preemergence herbicides registered in snap bean have poor waterhemp (Amaranthus tuberculatus) control, a major pest in snap bean production. Waterhemp and other difficult-to-control weeds can be managed by flumioxazin, an herbicide that inhibits protoporphyrinogen oxidase (PPO). However, there is limited knowledge about crop tolerance to this herbicide. We aimed to quantify the degree of snap bean tolerance to flumioxazin and explore the underlying mechanisms.

Methods

We investigated the genetic basis of herbicide tolerance using genome-wide association mapping approach utilizing field-collected data from a snap bean diversity panel, combined with gene expression data of cultivars with contrasting response. The response to a preemergence application of flumioxazin was measured by assessing plant population density and shoot biomass variables.

Results

Snap bean tolerance to flumioxazin is associated with a single genomic location in chromosome 02. Tolerance is influenced by several factors, including those that are indirectly affected by seed size/weight and those that directly impact the herbicide's metabolism and protect the cell from reactive oxygen species-induced damage. Transcriptional profiling and co-expression network analysis identified biological pathways likely involved in flumioxazin tolerance, including oxidoreductase processes and programmed cell death. Transcriptional regulation of genes involved in those processes is possibly orchestrated by a transcription factor located in the region identified in the GWAS analysis. Several entries belonging to the Romano class, including Bush Romano 350, Roma II, and Romano Purpiat presented high levels of tolerance in this study. The alleles identified in the diversity panel that condition snap bean tolerance to flumioxazin shed light on a novel mechanism of herbicide tolerance and can be used in crop improvement.

Genome-wide association study identifies novel loci and candidate genes for rust resistance in wheat (Triticum aestivum L.)

BACKGROUND

Wheat rusts are important biotic stresses, development of rust resistant cultivars through molecular approaches is both economical and sustainable. Extensive phenotyping of large mapping populations under diverse production conditions and high-density genotyping would be the ideal strategy to identify major genomic regions for rust resistance in wheat. The genome-wide association study (GWAS) population of 280 genotypes was genotyped using a 35 K Axiom single nucleotide polymorphism (SNP) array and phenotyped at eight, 10, and, 10 environments, respectively for stem/black rust (SR), stripe/yellow rust (YR), and leaf/brown rust (LR).

RESULTS

Forty-one Bonferroni corrected marker-trait associations (MTAs) were identified, including 17 for SR and 24 for YR. Ten stable MTAs and their best combinations were also identified. For YR, AX-94990952 on 1A + AX-95203560 on 4A + AX-94723806 on 3D + AX-95172478 on 1A showed the best combination with an average co-efficient of infection (ACI) score of 1.36. Similarly, for SR, AX-94883961 on 7B + AX-94843704 on 1B and AX-94883961 on 7B + AX-94580041 on 3D + AX-94843704 on 1B showed the best combination with an ACI score of around 9.0. The genotype PBW827 have the best MTA combinations for both YR and SR resistance. In silico study identifies key prospective candidate genes that are located within MTA regions. Further, the expression analysis revealed that 18 transcripts were upregulated to the tune of more than 1.5 folds including 19.36 folds (TraesCS3D02G519600) and 7.23 folds (TraesCS2D02G038900) under stress conditions compared to the control conditions. Furthermore, highly expressed genes in silico under stress conditions were analyzed to find out the potential links to the rust phenotype, and all four genes were found to be associated with the rust phenotype.

CONCLUSION

The identified novel MTAs, particularly stable and highly expressed MTAs are valuable for further validation and subsequent application in wheat rust resistance breeding. The genotypes with favorable MTA combinations can be used as prospective donors to develop elite cultivars with YR and SR resistance.

Major impacts of widespread structural variation on sorghum

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

Transcription factor encoding gene OsC1 regulates leaf sheath color through anthocyanidin metabolism in Oryza rufipogon and Oryza sativa

Carbohydrates, proteins, lipids, minerals and vitamins are nutrient substances commonly seen in rice grains, but anthocyanidin, with benefit for plant growth and animal health, exists mainly in the common wild rice but hardly in the cultivated rice. To screen the rice germplasm with high intensity of anthocyanidins and identify the variations, we used metabolomics technique and detected significant different accumulation of anthocyanidins in common wild rice (Oryza rufipogon, with purple leaf sheath) and cultivated rice (Oryza sativa, with green leaf sheath). In this study, we identified and characterized a well-known MYB transcription factor, OsC1, through phenotypic (leaf sheath color) and metabolic (metabolite profiling) genome-wide association studies (pGWAS and mGWAS) in 160 common wild rice (O. rufipogon) and 151 cultivated (O. sativa) rice varieties. Transgenic experiments demonstrated that biosynthesis and accumulation of cyanidin-3-Galc, cyanidin 3-O-rutinoside and cyanidin O-syringic acid, as well as purple pigmentation in leaf sheath were regulated by OsC1. A total of 25 sequence variations of OsC1 constructed 16 functional haplotypes (higher accumulation of the three anthocyanidin types within purple leaf sheath) and 9 non-functional haplotypes (less accumulation of anthocyanidins within green leaf sheath). Three haplotypes of OsC1 were newly identified in our germplasm, which have potential values in functional genomics and molecular breeding of rice. Gene-to-metabolite analysis by mGWAS and pGWAS provides a useful and efficient tool for functional gene identification and omics-based crop genetic improvement.

Multi-model genome-wide association studies for appearance quality in rice

Improving the quality of the appearance of rice is critical to meet market acceptance. Mining putative quality-related genes has been geared towards the development of effective breeding approaches for rice. In the present study, two SL-GWAS (CMLM and MLM) and three ML-GWAS (FASTmrEMMA, mrMLM, and FASTmrMLM) genome-wide association studies were conducted in a subset of 3K-RGP consisting of 198 rice accessions with 553,831 SNP markers. A total of 594 SNP markers were identified using the mixed linear model method for grain quality traits. Additionally, 70 quantitative trait nucleotides (QTNs) detected by the ML-GWAS models were strongly associated with grain aroma (AR), head rice recovery (HRR, %), and percentage of grains with chalkiness (PGC, %). Finally, 39 QTNs were identified using single- and multi-locus GWAS methods. Among the 39 reliable QTNs, 20 novel QTNs were identified for the above-mentioned three quality-related traits. Based on annotation and previous studies, four functional candidate genes (LOC_Os01g66110, LOC_Os01g66140, LOC_Os07g44910, and LOC_Os02g14120) were found to influence AR, HRR (%), and PGC (%), which could be utilized in rice breeding to improve grain quality traits.

Genome-wide association study on resistance of cultivated soybean to Fusarium oxysporum root rot in Northeast China

BACKGROUND

Fusarium oxysporum is a prevalent fungal pathogen that diminishes soybean yield through seedling disease and root rot. Preventing Fusarium oxysporum root rot (FORR) damage entails on the identification of resistance genes and developing resistant cultivars. Therefore, conducting fine mapping and marker development for FORR resistance genes is of great significance for fostering the cultivation of resistant varieties. In this study, 350 soybean germplasm accessions, mainly from Northeast China, underwent genotyping using the SoySNP50K Illumina BeadChip, which includes 52,041 single nucleotide polymorphisms (SNPs). Their resistance to FORR was assessed in a greenhouse. Genome-wide association studies utilizing the general linear model, mixed linear model, compressed mixed linear model, and settlement of MLM under progressively exclusive relationship models were conducted to identify marker-trait associations while effectively controlling for population structure.

RESULTS

The results demonstrated that these models effectively managed population structure. Eight SNP loci significantly associated with FORR resistance in soybean were detected, primarily located on Chromosome 6. Notably, there was a strong linkage disequilibrium between the large-effect SNPs ss715595462 and ss715595463, contributing substantially to phenotypic variation. Within the genetic interval encompassing these loci, 28 genes were present, with one gene Glyma.06G088400 encoding a protein kinase family protein containing a leucine-rich repeat domain identified as a potential candidate gene in the reference genome of Williams82. Additionally, quantitative real-time reverse transcription polymerase chain reaction analysis evaluated the gene expression levels between highly resistant and susceptible accessions, focusing on primary root tissues collected at different time points after F. oxysporum inoculation. Among the examined genes, only this gene emerged as the strongest candidate associated with FORR resistance.

CONCLUSIONS

The identification of this candidate gene Glyma.06G088400 improves our understanding of soybean resistance to FORR and the markers strongly linked to resistance can be beneficial for molecular marker-assisted selection in breeding resistant soybean accessions against F. oxysporum.

Genetic architecture of head rice and rice chalky grain percentages using genome-wide association studies

High head rice and low chalky grain percentages are key grain quality traits selected in developing rice cultivars. The objectives of this research were to characterize the phenotypic variation of head rice and chalky grain percentages in a diverse collection of rice accessions, identify single nucleotide polymorphism (SNP) markers associated with each of these traits using genome-wide association studies (GWAS), and identify putative candidate genes linked to the SNPs identified by GWAS. Diverse rice varieties, landraces, and breeding lines were grown at the Texas A&M AgriLife Research Center in Beaumont. Head rice percentages (HRP) and chalky grain percentages (CGP) of 195 and 199 non-waxy accessions were estimated in 2018 and 2019, respectively. Phenotypic data were analyzed along with 854,832 SNPs using three statistical models: mixed linear model (MLM), multi-locus mixed model (MLMM), and fixed and random model circulating probability unification (FarmCPU). Significant variations in HRP and CGP were observed between rice accessions. Two significant marker-trait associations (MTAs) were detected on chromosomes 1 and 2, respectively, based on best linear unbiased prediction (BLUP) values in 2018, while in 2019, one SNP was significantly associated with HRP in each of chromosomes 6, 8, 9, and 11, and two in chromosome 7. CGP was significantly associated with five SNPs located in chromosomes 2, 4, 6, and 8 in the 2018 study and ten SNPs in chromosomes 1, 2, 3, 4, 7, 8, 11, and 12 in the 2019 study. The SNPs are located within or linked to putative candidate genes involved in HRP and CGP. This study reports five and ten novel MTAs for HRP and CGP, respectively, while three and five MTAs co-located with previously reported quantitative trait loci for HRP and CGP, respectively. The validation of candidate genes for their roles in determining HRP and CGP is necessary to design functional molecular markers that can be used to effectively develop rice cultivars with desirable grain quality.

Regularized multi-trait multi-locus linear mixed models for genome-wide association studies and genomic selection in crops

BACKGROUND

We consider two key problems in genomics involving multiple traits: multi-trait genome wide association studies (GWAS), where the goal is to detect genetic variants associated with the traits; and multi-trait genomic selection (GS), where the emphasis is on accurately predicting trait values. Multi-trait linear mixed models build on the linear mixed model to jointly model multiple traits. Existing estimation methods, however, are limited to the joint analysis of a small number of genotypes; in fact, most approaches consider one SNP at a time. Estimating multi-dimensional genetic and environment effects also results in considerable computational burden. Efficient approaches that incorporate regularization into multi-trait linear models (no random effects) have been recently proposed to identify genomic loci associated with multiple traits (Yu et al. in Multitask learning using task clustering with applications to predictive modeling and GWAS of plant varieties. arXiv:1710.01788 , 2017; Yu et al in Front Big Data 2:27, 2019), but these ignore population structure and familial relatedness (Yu et al in Nat Genet 38:203-208, 2006).

RESULTS

This work addresses this gap by proposing a novel class of regularized multi-trait linear mixed models along with scalable approaches for estimation in the presence of high-dimensional genotypes and a large number of traits. We evaluate the effectiveness of the proposed methods using datasets in maize and sorghum diversity panels, and demonstrate benefits in both achieving high prediction accuracy in GS and in identifying relevant marker-trait associations.

CONCLUSIONS

The proposed regularized multivariate linear mixed models are relevant for both GWAS and GS. We hope that they will facilitate agronomy-related research in plant biology and crop breeding endeavors.

Genome-wide association studies using multi-models and multi-SNP datasets provide new insights into pasmo resistance in flax

Introduction

Flax (Linum usitatissimum L.) is an economically important crop due to its oil and fiber. However, it is prone to various diseases, including pasmo caused by the fungus Septoria linicola.

Methods

In this study, we conducted field evaluations of 445 flax accessions over a five-year period (2012-2016) to assess their resistance to pasmo A total of 246,035 single nucleotide polymorphisms (SNPs) were used for genetic analysis. Four statistical models, including the single-locus model GEMMA and the multi-locus models FarmCPU, mrMLM, and 3VmrMLM, were assessed to identify quantitative trait nucleotides (QTNs) associated with pasmo resistance.

Results

We identified 372 significant QTNs or 132 tag QTNs associated with pasmo resistance from five pasmo resistance datasets (PAS2012-PAS2016 and the 5-year average, namely PASmean) and three genotypic datasets (the all SNPs/ALL, the gene-based SNPs/GB and the RGA-based SNPs/RGAB). The tag QTNs had R2 values of 0.66-16.98% from the ALL SNP dataset, 0.68-20.54%from the GB SNP dataset, and 0.52-22.42% from the RGAB SNP dataset. Of these tag QTNs, 93 were novel. Additionally, 37 resistance gene analogs (RGAs)co-localizing with 39 tag QTNs were considered as potential candidates for controlling pasmo resistance in flax and 50 QTN-by-environment interactions(QEIs) were identified to account for genes by environmental interactions. Nine RGAs were predicted as candidate genes for ten QEIs.

Discussion

Our results suggest that pasmo resistance in flax is polygenic and potentially influenced by environmental factors. The identified QTNs provide potential targets for improving pasmo resistance in flax breeding programs. This study sheds light on the genetic basis of pasmo resistance and highlights the importance of considering both genetic and environmental factors in breeding programs for flax.

Genome-wide association study of fiber yield-related traits uncovers the novel genomic regions and candidate genes in Indian upland cotton (Gossypium hirsutum L.)

Upland cotton (Gossypium hirsutum L.) is a major fiber crop that is cultivated worldwide and has significant economic importance. India harbors the largest area for cotton cultivation, but its fiber yield is still compromised and ranks 22nd in terms of productivity. Genetic improvement of cotton fiber yield traits is one of the major goals of cotton breeding, but the understanding of the genetic architecture underlying cotton fiber yield traits remains limited and unclear. To better decipher the genetic variation associated with fiber yield traits, we conducted a comprehensive genome-wide association mapping study using 117 Indian cotton germplasm for six yield-related traits. To accomplish this, we generated 2,41,086 high-quality single nucleotide polymorphism (SNP) markers using genotyping-by-sequencing (GBS) methods. Population structure, PCA, kinship, and phylogenetic analyses divided the germplasm into two sub-populations, showing weak relatedness among the germplasms. Through association analysis, 205 SNPs and 134 QTLs were identified to be significantly associated with the six fiber yield traits. In total, 39 novel QTLs were identified in the current study, whereas 95 QTLs overlapped with existing public domain data in a comparative analysis. Eight QTLs, qGhBN_SCY_D6-1, qGhBN_SCY_D6-2, qGhBN_SCY_D6-3, qGhSI_LI_A5, qGhLI_SI_A13, qGhLI_SI_D9, qGhBW_SCY_A10, and qGhLP_BN_A8 were identified. Gene annotation of these fiber yield QTLs revealed 2,509 unique genes. These genes were predominantly enriched for different biological processes, such as plant cell wall synthesis, nutrient metabolism, and vegetative growth development in the gene ontology (GO) enrichment study. Furthermore, gene expression analysis using RNAseq data from 12 diverse cotton tissues identified 40 candidate genes (23 stable and 17 novel genes) to be transcriptionally active in different stages of fiber, ovule, and seed development. These findings have revealed a rich tapestry of genetic elements, including SNPs, QTLs, and candidate genes, and may have a high potential for improving fiber yield in future breeding programs for Indian cotton.

Effective use of legacy data in a genome-wide association studies improves the credibility of quantitative trait loci detection in rice

Genome-wide association studies (GWASs) are used to detect quantitative trait loci (QTL) using genomic and phenotypic data as inputs. While genomic data are obtained with high throughput and low cost, obtaining phenotypic data requires a large amount of effort and time. In past breeding programs, researchers and breeders have conducted a large number of phenotypic surveys and accumulated results as legacy data. In this study, we conducted a GWAS using phenotypic data of temperate japonica rice (Oryza sativa) varieties from a public database. The GWAS using the legacy data detected several known agriculturally important genes, indicating reliability of the legacy data for GWAS. By comparing the GWAS using legacy data (L-GWAS) and a GWAS using phenotypic data that we measured (M-GWAS), we detected reliable QTL for agronomically important traits. These results suggest that an L-GWAS is a strong alternative to replicate tests to confirm the reproducibility of QTL detected by an M-GWAS. In addition, because legacy data have often been accumulated for many traits, it is possible to evaluate the pleiotropic effect of the QTL identified for the specific trait that we focused on with respect to various other traits. This study demonstrates the effectiveness of using legacy data for GWASs and proposes the use of legacy data to accelerate genomic breeding.

Genome-wide association study, haplotype analysis, and genomic prediction reveal the genetic basis of yield-related traits in soybean (Glycine max L.)

Identifying the genetic components underlying yield-related traits in soybean is crucial for improving its production and productivity. Here, 211 soybean genotypes were evaluated across six environments for four yield-related traits, including seed yield per plant (SYP), number of pods per plant number of seeds per plant and 100-seed weight (HSW). Genome-wide association study (GWAS) and genomic prediction (GP) analyses were performed using 12,617 single nucleotide polymorphism markers from NJAU 355K SoySNP Array. A total of 57 SNPs were significantly associated with four traits across six environments and a combined environment using five Genome-wide association study models. Out of these, six significant SNPs were consistently identified in more than three environments using multiple GWAS models. The genomic regions (±670 kb) flanking these six consistent SNPs were considered stable QTL regions. Gene annotation and in silico expression analysis revealed 15 putative genes underlying the stable QTLs that might regulate soybean yield. Haplotype analysis using six significant SNPs revealed various allelic combinations regulating diverse phenotypes for the studied traits. Furthermore, the GP analysis revealed that accurate breeding values for the studied soybean traits is attainable at an earlier generation. Our study paved the way for increasing soybean yield performance within a short breeding cycle.

A multi-trait multi-locus stepwise approach for conducting GWAS on correlated traits

The ability to accurately quantify the simultaneous effect of multiple genomic loci on multiple traits is now possible due to current and emerging high-throughput genotyping and phenotyping technologies. To date, most efforts to quantify these genotype-to-phenotype relationships have focused on either multi-trait models that test a single marker at a time or multi-locus models that quantify associations with a single trait. Therefore, the purpose of this study was to compare the performance of a multi-trait, multi-locus stepwise (MSTEP) model selection procedure we developed to (a) a commonly used multi-trait single-locus model and (b) a univariate multi-locus model. We used real marker data in maize (Zea mays L.) and soybean (Glycine max L.) to simulate multiple traits controlled by various combinations of pleiotropic and nonpleiotropic quantitative trait nucleotides (QTNs). In general, we found that both multi-trait models outperformed the univariate multi-locus model, especially when analyzing a trait of low heritability. For traits controlled by either a combination of pleiotropic and nonpleiotropic QTNs or a large number of QTNs (i.e., 50), our MSTEP model often outperformed at least one of the two alternative models. When applied to the analysis of two tocochromanol-related traits in maize grain, MSTEP identified the same peak-associated marker that has been reported in a previous study. We therefore conclude that MSTEP is a useful addition to the suite of statistical models that are commonly used to gain insight into the genetic architecture of agronomically important traits.

Average semivariance directly yields accurate estimates of the genomic variance in complex trait analyses

Many important traits in plants, animals, and microbes are polygenic and challenging to improve through traditional marker-assisted selection. Genomic prediction addresses this by incorporating all genetic data in a mixed model framework. The primary method for predicting breeding values is genomic best linear unbiased prediction, which uses the realized genomic relationship or kinship matrix (K) to connect genotype to phenotype. Genomic relationship matrices share information among entries to estimate the observed entries’ genetic values and predict unobserved entries’ genetic values. One of the main parameters of such models is genomic variance (σg2), or the variance of a trait associated with a genome-wide sample of DNA polymorphisms, and genomic heritability (hg2); however, the seminal papers introducing different forms of K often do not discuss their effects on the model estimated variance components despite their importance in genetic research and breeding. Here, we discuss the effect of several standard methods for calculating the genomic relationship matrix on estimates of σg2 and hg2. With current approaches, we found that the genomic variance tends to be either overestimated or underestimated depending on the scaling and centering applied to the marker matrix (Z), the value of the average diagonal element of K, and the assortment of alleles and heterozygosity (H) in the observed population. Using the average semivariance, we propose a new matrix, KASV, that directly yields accurate estimates of σg2 and hg2 in the observed population and produces best linear unbiased predictors equivalent to routine methods in plants and animals.

Phenotypic variation and genome-wide association studies of main culm panicle node number, maximum node production rate, and degree-days to heading in rice

Background

Grain yield is a complex trait that results from interaction between underlying phenotypic traits and climatic, edaphic, and biotic variables. In rice, main culm panicle node number (MCPNN; the node number on which the panicle is borne) and maximum node production rate (MNPR; the number of leaves that emerge per degree-day > 10°C) are primary phenotypic plant traits that have significant positive direct effects on yield-related traits. Degree-days to heading (DDTH), which has a significant positive effect on grain yield, is influenced by the interaction between MCPNN and MNPR. The objective of this research is to assess the phenotypic variation of MCPNN, MNPR, and DDTH in a panel of diverse rice accessions, determine regions in the rice genome associated with these traits using genome-wide association studies (GWAS), and identify putative candidate genes that control these traits.

Results

Considerable variation was observed for the three traits in a 220-genotype diverse rice population. MCPNN ranged from 8.1 to 20.9 nodes in 2018 and from 9.9 to 21.0 nodes in 2019. MNPR ranged from 0.0097 to 0.0214 nodes/degree day > 10°C in 2018 and from 0.0108 to 0.0193 nodes/degree-day > 10°C in 2019. DDTH ranged from 713 to 2,345 degree-days > 10°C in 2018 and from 778 to 2,404 degree-days > 10°C in 2019.

Thirteen significant (P < 2.91 x 10^-7) trait-single nucleotide polymorphism (SNP) associations were identified using the multilocus mixed linear model for GWAS. Significant associations between MCPNN and three SNPs in chromosome 2 (S02_12032235, S02_11971745, and S02_12030176) were detected with both the 2018 and best linear unbiased prediction (BLUP) datasets. Nine SNPs in chromosome 6 (S06_1970442, S06_2310856, S06_2550351, S06_1968653, S06_2296852, S06_1968680, S06_1968681, S06_1970597, and S06_1970602) were significantly associated with MNPR in the 2019 dataset. One SNP in chromosome 11 (S11_29358169) was significantly associated with the DDTH in the BLUP dataset.

Conclusions

This study identifies SNP markers that are putatively associated with MCPNN, MNPR, and DDTH. Some of these SNPs were located within or near gene models, which identify possible candidate genes involved in these traits. Validation of the putative candidate genes through expression and gene editing analyses are necessary to confirm their roles in regulating MCPNN, MNPR, and DDTH. Identifying the underlying genetic basis for primary phenotypic traits MCPNN and MNPR could lead to the development of fast and efficient approaches for their estimation, such as marker-assisted selection and gene editing, which is essential in increasing breeding efficiency and enhancing grain yield in rice. On the other hand, DDTH is a resultant variable that is highly affected by nitrogen and water management, plant density, and several other factors.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08629-y.

Genome Wide Association Study Identifies Candidate Genes Related to the Earlywood Tracheid Properties in Picea crassifolia Kom

Picea crassifolia Kom. is one of the timber and ecological conifers in China and its wood tracheid traits directly affect wood formation and adaptability under harsh environment. Molecular studies on P. crassifolia remain inadequate because relatively few genes have been associated with these traits. To identify markers and candidate genes that can potentially be used for genetic improvement of wood tracheid traits, we examined 106 clones of P. crassifolia, and investigated phenotypic data for 14 wood tracheid traits before specific-locus amplified fragment sequencing (SLAF-seq) was employed to perform a genome wide association study (GWAS). Subsequently, the results were used to screen single nucleotide polymorphism (SNP) loci and candidate genes that exhibited a significant correlation with the studied traits. We developed 4,058,883 SLAF-tags and 12,275,765 SNP loci, and our analyses identified a total of 96 SNP loci that showed significant correlations with three earlywood tracheid traits using a mixed linear model (MLM). Next, candidate genes were screened in the 100 kb zone (50 kb upstream, 50 kb downstream) of each of the SNP loci, whereby 67 candidate genes were obtained in earlywood tracheid traits, including 34 genes of known function and 33 genes of unknown function. We provide the most significant SNP for each trait-locus combination and candidate genes occurring within the GWAS hits. These resources provide a foundation for the development of markers that could be used in wood traits improvement and candidate genes for the development of earlywood tracheid in P. crassifolia.

Genome-wide association study of cassava starch paste properties

An understanding of cassava starch paste properties (CSPP) can contribute to the selection of clones with differentiated starches. This study aimed to identify genomic regions associated with CSPP using different genome-wide association study (GWAS) methods (MLM, MLMM, and Farm-CPU). The GWAS was performed using 23,078 single-nucleotide polymorphisms (SNPs). The rapid viscoanalyzer (RVA) parameters were pasting temperature (PastTemp), peak viscosity (PeakVisc), hot-paste viscosity (Hot-PVisc), cool-paste viscosity (Cold-PVisc), final viscosity (FinalVis), breakdown (BreDow), and setback (Setback). Broad phenotypic and molecular diversity was identified based on the genomic kinship matrix. The broad-sense heritability estimates (h²) ranged from moderate to high magnitudes (0.66 to 0.76). The linkage disequilibrium (LD) declined to between 0.3 and 2.0 Mb (r² <0.1) for most chromosomes, except chromosome 17, which exhibited an extensive LD. Thirteen SNPs were found to be significantly associated with CSPP, on chromosomes 3, 8, 17, and 18. Only the BreDow trait had no associated SNPs. The regional marker-trait associations on chromosome 18 indicate a LD block between 2907312 and 3567816 bp and that SNP S18_3081635 was associated with SetBack, FinalVis, and Cold-PVisc (all three GWAS methods) and with Hot-PVisc (MLM), indicating that this SNP can track these four traits simultaneously. The variance explained by the SNPs ranged from 0.13 to 0.18 for SetBack, FinalVis, and Cold-PVisc and from 0.06 to 0.09 for PeakVisc and Hot-PVisc. The results indicated additive effects of the genetic control of Cold-PVisc, FinalVis, Hot-PVisc, and SetBack, especially on the large LD block on chromosome 18. One transcript encoding the glycosyl hydrolase family 35 enzymes on chromosome 17 and one encoding the mannose-p-dolichol utilization defect 1 protein on chromosome 18 were the most likely candidate genes for the regulation of CSPP. These results underline the potential for the assisted selection of high-value starches to improve cassava root quality through breeding programs.

Interpretation of Manhattan Plots and Other Outputs of Genome-Wide Association Studies

With increasing marker density, estimation of recombination rate between a marker and a causal mutation using linkage analysis becomes less important. Instead, linkage disequilibrium (LD) becomes the major indicator for gene mapping through genome-wide association studies (GWAS). In addition to the linkage between the marker and the causal mutation, many other factors may contribute to the LD, including population structure and cryptic relationships among individuals. As statistical methods and software evolve to improve statistical power and computing speed in GWAS, the corresponding outputs must also evolve to facilitate the interpretation of input data, the analytical process, and final association results. In this chapter, our descriptions focus on (1) considerations in creating a Manhattan plot displaying the strength of LD and locations of markers across a genome; (2) criteria for genome-wide significance threshold and the different appearance of Manhattan plots in single-locus and multiple-locus models; (3) exploration of population structure and kinship among individuals; (4) quantile-quantile (QQ) plot; (5) LD decay across the genome and LD between the associated markers and their neighbors; (6) exploration of individual and marker information on Manhattan and QQ plots via interactive visualization using HTML. The ultimate objective of this chapter is to help users to connect input data to GWAS outputs to balance power and false positives, and connect GWAS outputs to the selection of candidate genes using LD extent.

Preparation and Curation of Omics Data for Genome-Wide Association Studies

With the development of large-scale molecular phenotyping platforms, genome-wide association studies have greatly developed, being no longer limited to the analysis of classical agronomic traits, such as yield or flowering time, but also embracing the dissection of the genetic basis of molecular traits. Data generated by omics platforms, however, pose some technical and statistical challenges to the classical methodology and assumptions of an association study. Although genotyping data are subject to strict filtering procedures, and several advanced statistical approaches are now available to adjust for population structure, less attention has been instead devoted to the preparation of omics data prior to GWAS. In the present chapter, we briefly present the methods to acquire profiling data from transcripts, proteins, and small molecules, and discuss the tools and possibilities to clean, normalize, and remove the unwanted variation from large datasets of molecular phenotypic traits prior to their use in GWAS.

Recent Advances in the Analysis of Cold Tolerance in Maize

Maize (Zea mays L.) is an annual grass that originated in tropical and subtropical regions of the New World. Maize is highly sensitive to cold stress during seed gemination and the seedling phase, which can lead to reductions in plant vigor and grain production. There are large differences in the morphological and physiological changes caused by cold stress among maize varieties. In general, cold tolerant varieties have a stronger ability to maintain such changes in traits related to seed germination, root phenotypes, and shoot photosynthesis. These morphological and physiological characteristics have been widely used to evaluate the cold tolerance of maize varieties in genetic analyses. In recent years, considerable progress has been made in elucidating the mechanisms of maize in response to cold tolerance. Several QTL, GWAS, and transcriptomic analyses have been conducted on various maize genotypes and populations that show large variations in cold tolerance, resulting in the discovery of hundreds of candidate cold regulation genes. Nevertheless, only a few candidate genes have been functionally characterized. In the present review, we summarize recent progress in molecular, physiological, genetic, and genomic analyses of cold tolerance in maize. We address the advantages of joint analyses that combine multiple genetic and genomic approaches to improve the accuracy of identifying cold regulated genes that can be further used in molecular breeding. We also discuss the involvement of long-distance signaling in plant cold tolerance. These novel insights will provide a better mechanistic understanding of cold tolerance in maize.

GWAS for genetics of complex quantitative traits: Genome to pangenome and SNPs to SVs and k-mers

The development of improved methods for genome-wide association studies (GWAS) for genetics of quantitative traits has been an active area of research during the last 25 years. This activity initially started with the use of mixed linear model (MLM), which was variously modified. During the last decade, however, with the availability of high throughput next generation sequencing (NGS) technology, development and use of pangenomes and novel markers including structural variations (SVs) and k-mers for GWAS has taken over as a new thrust area of research. Pangenomes and SVs are now available in humans, livestock, and a number of plant species, so that these resources along with k-mers are being used in GWAS for exploring additional genetic variation that was hitherto not available for analysis. These developments have resulted in significant improvement in GWAS methodology for detection of marker-trait associations (MTAs) that are relevant to human healthcare and crop improvement.

The genetic architecture of leaf stable carbon isotope composition in Zea mays and the effect of transpiration efficiency on leaf elemental accumulation

With increased demand on freshwater resources for agriculture, it is imperative that more water-use efficient crops are developed. Leaf stable carbon isotope composition, δ13C, is a proxy for transpiration efficiency and a possible tool for breeders, but the underlying mechanisms effecting δ13C in C4 plants are not known. It has been suggested that differences in specific leaf area (SLA), which potentially reflects variation in internal CO2 diffusion, can impact leaf δ13C. Furthermore, although it is known that water movement is important for elemental uptake, it is not clear how manipulation of transpiration for increased water-use efficiency may impact nutrient accumulation. Here, we characterize the genetic architecture of leaf δ13C and test its relationship to SLA and the ionome in five populations of maize. Five significant QTL for leaf δ13C were identified, including novel QTL as well as some that were identified previously in maize kernels. One of the QTL regions contains an Erecta-like gene, the ortholog of which has been shown to regulate transpiration efficiency and leaf δ13C in Arabidopsis. QTL for δ13C were located in the same general chromosome region, but slightly shifted, when comparing data from two different years. Our data does not support a relationship between δ13C and SLA, and of the 19 elements analyzed, only a weak correlation between molybdenum and δ13C was detected. Together these data add to the genetic understanding of leaf δ13C in maize and suggest that improvements to plant water use may be possible without significantly influencing elemental homeostasis.

Genome-wide association study for deoxynivalenol production and aggressiveness in wheat and rye head blight by resequencing 92 isolates of Fusarium culmorum

Background

Fusarium culmorum is an important pathogen causing head blight of cereals in Europe. This disease is of worldwide importance leading to reduced yield, grain quality, and contamination by mycotoxins. These mycotoxins are harmful for livestock and humans; therefore, many countries have strict regulatory limits for raw materials and processed food. Extensive genetic diversity is described among field populations of F. culmorum isolates for aggressiveness and production of the trichothecene mycotoxin deoxynivalenol (DON). However, the causes for this quantitative variation are not clear, yet. We analyzed 92 isolates sampled from different field populations in Germany, Russia, and Syria together with an international collection for aggressiveness and DON production in replicated field experiments at two locations in two years with two hosts, wheat and rye. The 30x coverage whole-genome resequencing of all isolates resulted in the identification of 130,389 high quality single nucleotide polymorphisms (SNPs) that were used for the first genome-wide association study in this phytopathogenic fungus.

Results

In wheat, 20 and 27 SNPs were detected for aggressiveness and DON content, respectively, of which 10 overlapped. Additionally, two different SNPs were significantly associated with aggressiveness in rye that were among those SNPs being associated with DON production in wheat. Most of the SNPs explained only a small proportion of genotypic variance (p_G), however, four SNPs were associated with major quantitative trait loci (QTLs) with p_G ranging from 12 to 48%. The QTL with the highest p_G was involved in DON production and associated with a SNP most probably located within the Tri4 gene.

Conclusions

The diversity of 92 isolates of F. culmorum were captured using a heuristic approach. Key phenotypic traits, SNPs, and candidate genes underlying aggressiveness and DON production were identified. Clearly, many QTLs are responsible for aggressiveness and DON content in wheat, both traits following a quantitative inheritance. Several SNPs involved in DON metabolism, among them the Tri4 gene of the trichothecene pathway, were inferred as important source of variation in fungal aggressiveness. Using this information underlying the phenotypic variation will be of paramount importance in evaluating strategies for successful resistance breeding.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07931-5.

High-resolution genome-wide association study pinpoints metal transporter and chelator genes involved in the genetic control of element levels in maize grain

Despite its importance to plant function and human health, the genetics underpinning element levels in maize grain remain largely unknown. Through a genome-wide association study in the maize Ames panel of nearly 2,000 inbred lines that was imputed with ∼7.7 million SNP markers, we investigated the genetic basis of natural variation for the concentration of 11 elements in grain. Novel associations were detected for the metal transporter genes rte2 (rotten ear2) and irt1 (iron-regulated transporter1) with boron and nickel, respectively. We also further resolved loci that were previously found to be associated with one or more of five elements (copper, iron, manganese, molybdenum, and/or zinc), with two metal chelator and five metal transporter candidate causal genes identified. The nas5 (nicotianamine synthase5) gene involved in the synthesis of nicotianamine, a metal chelator, was found associated with both zinc and iron and suggests a common genetic basis controlling the accumulation of these two metals in the grain. Furthermore, moderate predictive abilities were obtained for the 11 elemental grain phenotypes with two whole-genome prediction models: Bayesian Ridge Regression (0.33–0.51) and BayesB (0.33–0.53). Of the two models, BayesB, with its greater emphasis on large-effect loci, showed ∼4–10% higher predictive abilities for nickel, molybdenum, and copper. Altogether, our findings contribute to an improved genotype-phenotype map for grain element accumulation in maize.

The importance of dominance and genotype-by-environment interactions on grain yield variation in a large-scale public cooperative maize experiment

High-dimensional and high-throughput genomic, field performance, and environmental data are becoming increasingly available to crop breeding programs, and their integration can facilitate genomic prediction within and across environments and provide insights into the genetic architecture of complex traits and the nature of genotype-by-environment interactions. To partition trait variation into additive and dominance (main effect) genetic and corresponding genetic-by-environment variances, and to identify specific environmental factors that influence genotype-by-environment interactions, we curated and analyzed genotypic and phenotypic data on 1918 maize (Zea mays L.) hybrids and environmental data from 65 testing environments. For grain yield, dominance variance was similar in magnitude to additive variance, and genetic-by-environment variances were more important than genetic main effect variances. Models involving both additive and dominance relationships best fit the data and modeling unique genetic covariances among all environments provided the best characterization of the genotype-by-environment interaction patterns. Similarity of relative hybrid performance among environments was modeled as a function of underlying weather variables, permitting identification of weather covariates driving correlations of genetic effects across environments. The resulting models can be used for genomic prediction of mean hybrid performance across populations of environments tested or for environment-specific predictions. These results can also guide efforts to incorporate high-throughput environmental data into genomic prediction models and predict values in new environments characterized with the same environmental characteristics.

An integrated framework reinstating the environmental dimension for GWAS and genomic selection in crops

Identifying mechanisms and pathways involved in gene-environment interplay and phenotypic plasticity is a long-standing challenge. It is highly desirable to establish an integrated framework with an environmental dimension for complex trait dissection and prediction. A critical step is to identify an environmental index that is both biologically relevant and estimable for new environments. With extensive field-observed complex traits, environmental profiles, and genome-wide single nucleotide polymorphisms for three major crops (maize, wheat, and oat), we demonstrated that identifying such an environmental index (i.e., a combination of environmental parameter and growth window) enables genome-wide association studies and genomic selection of complex traits to be conducted with an explicit environmental dimension. Interestingly, genes identified for two reaction-norm parameters (i.e., intercept and slope) derived from flowering time values along the environmental index were less colocalized for a diverse maize panel than for wheat and oat breeding panels, agreeing with the different diversity levels and genetic constitutions of the panels. In addition, we showcased the usefulness of this framework for systematically forecasting the performance of diverse germplasm panels in new environments. This general framework and the companion CERIS-JGRA analytical package should facilitate biologically informed dissection of complex traits, enhanced performance prediction in breeding for future climates, and coordinated efforts to enrich our understanding of mechanisms underlying phenotypic variation.

Diversifying maize genomic selection models

Genomic selection (GS) is one of the most powerful tools available for maize breeding. Its use of genome-wide marker data to estimate breeding values translates to increased genetic gains with fewer breeding cycles. In this review, we cover the history of GS and highlight particular milestones during its adaptation to maize breeding. We discuss how GS can be applied to developing superior maize inbreds and hybrids. Additionally, we characterize refinements in GS models that could enable the encapsulation of non-additive genetic effects, genotype by environment interactions, and multiple levels of the biological hierarchy, all of which could ultimately result in more accurate predictions of breeding values. Finally, we suggest the stages in a maize breeding program where it would be beneficial to apply GS. Given the current sophistication of high-throughput phenotypic, genotypic, and other -omic level data currently available to the maize community, now is the time to explore the implications of their incorporation into GS models and thus ensure that genetic gains are being achieved as quickly and efficiently as possible.

Validation of Functional Polymorphisms Affecting Maize Plant Height by Unoccupied Aerial Systems (UAS) Discovers Novel Temporal Phenotypes

Plant height (PHT) in maize (Zea mays L.) has been scrutinized genetically and phenotypically due to relationship with other agronomically valuable traits (e.g. yield). Heritable variation of PHT is determined by many discovered quantitative trait loci (QTLs); however, phenotypic effects of such loci often lack validation across environments and genetic backgrounds, especially in the hybrid state grown by farmers rather than the inbred state more often used by geneticists. A previous genome wide association study using a topcrossed hybrid diversity panel identified two novel quantitative trait variants (QTVs) controlling both PHT and grain yield. Here, heterogeneous inbred families demonstrated that these two loci, characterized by two single nucleotide polymorphisms (SNPs), cause phenotypic variation in inbred lines, but that size of these effects were variable across four different genetic backgrounds, ranging from 1 to 10 cm. Weekly unoccupied aerial system flights demonstrated the two SNPs had larger effects, varying from 10 to 25 cm, in early growth while effects decreased towards the end of the season. These results show that allelic effect sizes of economically valuable loci are both dynamic in temporal growth and dynamic across genetic backgrounds, resulting in informative phenotypic variability overlooked following traditional phenotyping methods. Public genotyping data shows recent favorable allele selection in elite temperate germplasm with little change across tropical backgrounds. As these loci remain rarer in tropical germplasm, with effects most visible early in growth, they are useful for breeding and selection to expand the genetic basis of maize.

Status and prospects of genome-wide association studies in plants

Genome-wide association studies (GWAS) have developed into a powerful and ubiquitous tool for the investigation of complex traits. In large part, this was fueled by advances in genomic technology, enabling us to examine genome-wide genetic variants across diverse genetic materials. The development of the mixed model framework for GWAS dramatically reduced the number of false positives compared with naïve methods. Building on this foundation, many methods have since been developed to increase computational speed or improve statistical power in GWAS. These methods have allowed the detection of genomic variants associated with either traditional agronomic phenotypes or biochemical and molecular phenotypes. In turn, these associations enable applications in gene cloning and in accelerated crop breeding through marker assisted selection or genetic engineering. Current topics of investigation include rare-variant analysis, synthetic associations, optimizing the choice of GWAS model, and utilizing GWAS results to advance knowledge of biological processes. Ongoing research in these areas will facilitate further advances in GWAS methods and their applications.

Integrated genomic selection for rapid improvement of crops

An increase in the rate of crop improvement is essential for achieving sustained food production and other needs of ever-increasing population. Genomic selection (GS) is a potential breeding tool that has been successfully employed in animal breeding and is being incorporated into plant breeding. GS promises accelerated breeding cycles through a rapid selection of superior genotypes. Numerous empirical and simulation studies on GS and realized impacts on improvement in the crop yields are recently being reported. For a holistic understanding of the technology, we briefly discuss the concept of genetic gain, GS methodology, its current status, advantages of GS over other breeding methods, prediction models, and the factors controlling prediction accuracy in GS. Also, integration of speed breeding and other novel technologies viz. high throughput genotyping and phenotyping technologies for enhancing the efficiency and pace of GS, followed by its prospective applications in varietal development programs is reviewed.

A large-scale genomic association analysis identifies the candidate causal genes conferring stripe rust resistance under multiple field environments

The incorporation of resistance genes into wheat commercial varieties is the ideal strategy to combat stripe or yellow rust (YR). In a search for novel resistance genes, we performed a large-scale genomic association analysis with high-density 660K single nucleotide polymorphism (SNP) arrays to determine the genetic components of YR resistance in 411 spring wheat lines. Following quality control, 371 972 SNPs were screened, covering over 50% of the high-confidence annotated gene space. Nineteen stable genomic regions harbouring 292 significant SNPs were associated with adult-plant YR resistance across nine environments. Of these, 14 SNPs were localized in the proximity of known loci widely used in breeding. Obvious candidate SNP variants were identified in certain confidence intervals, such as the cloned gene Yr18 and the major locus on chromosome 2BL, despite a large extent of linkage disequilibrium. The number of causal SNP variants was refined using an independent validation panel and consideration of the estimated functional importance of each nucleotide polymorphism. Interestingly, four natural polymorphisms causing amino acid changes in the gene TraesCS2B01G513100 that encodes a serine/threonine protein kinase (STPK) were significantly involved in YR responses. Gene expression and mutation analysis confirmed that STPK played an important role in YR resistance. PCR markers were developed to identify the favourable TraesCS2B01G513100 haplotype for marker-assisted breeding. These results demonstrate that high-resolution SNP-based GWAS enables the rapid identification of putative resistance genes and can be used to improve the efficiency of marker-assisted selection in wheat disease resistance breeding.

Comparison of Single-Trait and Multi-Trait Genome-Wide Association Models and Inclusion of Correlated Traits in the Dissection of the Genetic Architecture of a Complex Trait in a Breeding Program

Unknown genetic architecture makes it difficult to characterize the genetic basis of traits and associated molecular markers because of the complexity of small effect quantitative trait loci (QTLs), environmental effects, and difficulty in phenotyping. Seedling emergence of wheat (Triticum aestivum L.) from deep planting, has a poorly understood genetic architecture, is a vital factor affecting stand establishment and grain yield, and is historically correlated with coleoptile length. This study aimed to dissect the genetic architecture of seedling emergence while accounting for correlated traits using one multi-trait genome-wide association study (MT-GWAS) model and three single-trait GWAS (ST-GWAS) models. The ST-GWAS models included one single-locus model [mixed-linear model (MLM)] and two multi-locus models [fixed and random model circulating probability unification (FarmCPU) and Bayesian information and linkage-disequilibrium iteratively nested keyway (BLINK)]. We conducted GWAS using two populations. The first population consisted of 473 varieties from a diverse association mapping panel phenotyped from 2015 to 2019. The second population consisted of 279 breeding lines phenotyped in 2015 in Lind, WA, with 40,368 markers. We also compared the inclusion of coleoptile length and markers associated with reduced height as covariates in our ST-GWAS models. ST-GWAS found 107 significant markers across 19 chromosomes, while MT-GWAS found 82 significant markers across 14 chromosomes. The FarmCPU and BLINK models, including covariates, were able to identify many small effect markers while identifying large effect markers on chromosome 5A. By using multi-locus model breeding, programs can uncover the complex nature of traits to help identify candidate genes and the underlying architecture of a trait, such as seedling emergence.

Genome-Wide Association Studies in Arabidopsis thaliana: Statistical Analysis and Network-Based Augmentation of Signals

Genome-wide association studies (GWAS) have proven effective at identifying genetic variants and genes that are associated with phenotypes in humans, animals, and plants. Since most phenotypes of plant species are complex traits regulated by many genes and their functional interactions, GWAS are increasing in popularity for genetic dissections of plant phenotypes. For the reference plant, Arabidopsis thaliana, detailed information on genetic variations became available with the completion of the 1001 Genomes Project, enabling highly resolved association mapping between chromosomal loci and complex traits. Improvements have been made in the statistical analysis methods for testing the significance of genotype-to-phenotype associations, thereby substantially reducing the confounding effects of population structures. Furthermore, there have been large efforts toward post-GWAS augmentation of signals via integration with other types of information to overcome the limited statistical power of GWAS. This chapter describes the stepwise procedure of GWAS in Arabidopsis, focusing on data analysis processes including preprocessing of genotype and phenotype data, statistical analysis to identify phenotype-associated chromosomal loci, identification of phenotype-associated genes based on the phenotype-associated loci, and finally network-based augmentation of GWAS signals to identify additional candidate genes for the phenotype.

Association Mapping and Development of Marker-Assisted Selection Tools for the Resistance to White Pine Blister Rust in the Alberta Limber Pine Populations

Since its introduction to North America in the early 1900s, white pine blister rust (WPBR) caused by the fungal pathogen Cronartium ribicola has resulted in substantial economic losses and ecological damage to native North American five-needle pine species. The high susceptibility and mortality of these species, including limber pine (Pinus flexilis), creates an urgent need for the development and deployment of resistant germplasm to support recovery of impacted populations. Extensive screening for genetic resistance to WPBR has been underway for decades in some species but has only started recently in limber pine using seed families collected from wild parental trees in the USA and Canada. This study was conducted to characterize Alberta limber pine seed families for WPBR resistance and to develop reliable molecular tools for marker-assisted selection (MAS). Open-pollinated seed families were evaluated for host reaction following controlled infection using C. ribicola basidiospores. Phenotypic segregation for presence/absence of stem symptoms was observed in four seed families. The segregation ratios of these families were consistent with expression of major gene resistance (MGR) controlled by a dominant R locus. Based on linkage disequilibrium (LD)-based association mapping used to detect single nucleotide polymorphism (SNP) markers associated with MGR against C. ribicola, MGR in these seed families appears to be controlled by Cr4 or other R genes in very close proximity to Cr4. These associated SNPs were located in genes involved in multiple molecular mechanisms potentially underlying limber pine MGR to C. ribicola, including NBS-LRR genes for recognition of C. ribicola effectors, signaling components, and a large set of defense-responsive genes with potential functions in plant effector-triggered immunity (ETI). Interactions of associated loci were identified for MGR selection in trees with complex genetic backgrounds. SNPs with tight Cr4-linkage were further converted to TaqMan assays to confirm their effectiveness as MAS tools. This work demonstrates the successful translation and deployment of molecular genetic knowledge into specific MAS tools that can be easily applied in a selection or breeding program to efficiently screen MGR against WPBR in Alberta limber pine populations.

Introgression Breeding in Barley: Perspectives and Case Studies

Changing production scenarios resulting from unstable climatic conditions are challenging crop improvement efforts. A deeper and more practical understanding of plant genetic resources is necessary if these assets are to be used effectively in developing improved varieties. In general, current varieties and potential varieties have a narrow genetic base, making them prone to suffer the consequences of new and different abiotic and biotic stresses that can reduce crop yield and quality. The deployment of genomic technologies and sophisticated statistical analysis procedures has generated a dramatic change in the way we characterize and access genetic diversity in crop plants, including barley. Various mapping strategies can be used to identify the genetic variants that lead to target phenotypes and these variants can be assigned coordinates in reference genomes. In this way, new genes and/or new alleles at known loci present in wild ancestors, germplasm accessions, land races, and un-adapted introductions can be located and targeted for introgression. In principle, the introgression process can now be streamlined and linkage drag reduced. In this review, we present an overview of (1) past and current efforts to identify diversity that can be tapped to improve barley yield and quality, and (2) case studies of our efforts to introgress resistance to stripe and stem rust from un-adapted germplasm. We conclude with a description of a modified Nested Association Mapping (NAM) population strategy that we are implementing for the development of multi-use naked barley for organic systems and share perspectives on the use of genome editing in introgression breeding.

Genome-Wide Association Study for Maize Leaf Cuticular Conductance Identifies Candidate Genes Involved in the Regulation of Cuticle Development

The cuticle, a hydrophobic layer of cutin and waxes synthesized by plant epidermal cells, is the major barrier to water loss when stomata are closed at night and under water-limited conditions. Elucidating the genetic architecture of natural variation for leaf cuticular conductance (g_c) is important for identifying genes relevant to improving crop productivity in drought-prone environments. To this end, we conducted a genome-wide association study of g_c of adult leaves in a maize inbred association panel that was evaluated in four environments (Maricopa, AZ, and San Diego, CA, in 2016 and 2017). Five genomic regions significantly associated with g_c were resolved to seven plausible candidate genes (ISTL1, two SEC14 homologs, cyclase-associated protein, a CER7 homolog, GDSL lipase, and β-D-XYLOSIDASE 4). These candidates are potentially involved in cuticle biosynthesis, trafficking and deposition of cuticle lipids, cutin polymerization, and cell wall modification. Laser microdissection RNA sequencing revealed that all these candidate genes, with the exception of the CER7 homolog, were expressed in the zone of the expanding adult maize leaf where cuticle maturation occurs. With direct application to genetic improvement, moderately high average predictive abilities were observed for whole-genome prediction of g_c in locations (0.46 and 0.45) and across all environments (0.52). The findings of this study provide novel insights into the genetic control of g_c and have the potential to help breeders more effectively develop drought-tolerant maize for target environments.

Whole-Genome Sequencing of the NARO World Rice Core Collection (WRC) as the Basis for Diversity and Association Studies

Genebanks provide access to diverse materials for crop improvement. To utilize and evaluate them effectively, core collections, such as the World Rice Core Collection (WRC) in the Genebank at the National Agriculture and Food Research Organization, have been developed. Because the WRC consists of 69 accessions with a high degree of genetic diversity, it has been used for >300 projects. To allow deeper investigation of existing WRC data and to further promote research using Genebank rice accessions, we performed whole-genome resequencing of these 69 accessions, examining their sequence variation by mapping against the Oryza sativa ssp. japonica Nipponbare genome. We obtained a total of 2,805,329 single nucleotide polymorphisms (SNPs) and 357,639 insertion-deletions. Based on the principal component analysis and population structure analysis of these data, the WRC can be classified into three major groups. We applied TASUKE, a multiple genome browser to visualize the different WRC genome sequences, and classified haplotype groups of genes affecting seed characteristics and heading date. TASUKE thus provides access to WRC genotypes as a tool for reverse genetics. We examined the suitability of the compact WRC population for genome-wide association studies (GWASs). Heading date, affected by a large number of quantitative trait loci (QTLs), was not associated with known genes, but several seed-related phenotypes were associated with known genes. Thus, for QTLs of strong effect, the compact WRC performed well in GWAS. This information enables us to understand genetic diversity in 37,000 rice accessions maintained in the Genebank and to find genes associated with different phenotypes. The sequence data have been deposited in DNA Data Bank of Japan Sequence Read Archive (DRA) (Supplementary Table S1).

Genomic Prediction Accuracy of Seven Breeding Selection Traits Improved by QTL Identification in Flax

Molecular markers are one of the major factors affecting genomic prediction accuracy and the cost of genomic selection (GS). Previous studies have indicated that the use of quantitative trait loci (QTL) as markers in GS significantly increases prediction accuracy compared with genome-wide random single nucleotide polymorphism (SNP) markers. To optimize the selection of QTL markers in GS, a set of 260 lines from bi-parental populations with 17,277 genome-wide SNPs were used to evaluate the prediction accuracy for seed yield (YLD), days to maturity (DTM), iodine value (IOD), protein (PRO), oil (OIL), linoleic acid (LIO), and linolenic acid (LIN) contents. These seven traits were phenotyped over four years at two locations. Identification of quantitative trait nucleotides (QTNs) for the seven traits was performed using three types of statistical models for genome-wide association study: two SNP-based single-locus (SS), seven SNP-based multi-locus (SM), and one haplotype-block-based multi-locus (BM) models. The identified QTNs were then grouped into QTL based on haplotype blocks. For all seven traits, 133, 355, and 1208 unique QTL were identified by SS, SM, and BM, respectively. A total of 1420 unique QTL were obtained by SS+SM+BM, ranging from 254 (OIL, LIO) to 361 (YLD) for individual traits, whereas a total of 427 unique QTL were achieved by SS+SM, ranging from 56 (YLD) to 128 (LIO). SS models alone did not identify sufficient QTL for GS. The highest prediction accuracies were obtained using single-trait QTL identified by SS+SM+BM for OIL (0.929 ± 0.016), PRO (0.893 ± 0.023), YLD (0.892 ± 0.030), and DTM (0.730 ± 0.062), and by SS+SM for LIN (0.837 ± 0.053), LIO (0.835 ± 0.049), and IOD (0.835 ± 0.041). In terms of the number of QTL markers and prediction accuracy, SS+SM outperformed other models or combinations thereof. The use of all SNPs or QTL of all seven traits significantly reduced the prediction accuracy of traits. The results further validated that QTL outperformed high-density genome-wide random markers, and demonstrated that the combined use of single and multi-locus models can effectively identify a comprehensive set of QTL that improve prediction accuracy, but further studies on detection and removal of redundant or false-positive QTL to maximize prediction accuracy and minimize the number of QTL markers in GS are warranted.

Resequencing 200 Flax Cultivated Accessions Identifies Candidate Genes Related to Seed Size and Weight and Reveals Signatures of Artificial Selection

Seed size and weight are key traits determining crop yield, which often undergo strongly artificial selection during crop domestication. Although seed sizes differ significantly between oil flax and fiber flax, the genetic basis of morphological differences and artificial selection characteristics in seed size remains largely unclear. Here we re-sequenced 200 flax cultivated accessions to generate a genome variation map based on chromosome assembly reference genomes. We provide evidence that oil flax group is the ancestor of cultivated flax, and the oil-fiber dual purpose group (OF) is the evolutionary intermediate transition state between oil and fiber flax. Genome-wide association studies (GWAS) were combined with LD Heatmap to identify candidate regions related to seed size and weight, then candidate genes were screened based on detailed functional annotations and estimation of nucleotide polymorphism effects. Using this strategy, we obtained 13 candidate genes related to seed size and weight. Selective sweeps analysis indicates human-involved selection of small seeds during the oil to fiber flax transition. Our study shows the existence of elite alleles for seed size and weight in flax germplasm and provides molecular insights into approaches for further improvement.

Exploring genetic architecture of grain yield and quality traits in a 16-way indica by japonica rice MAGIC global population

Identification of Quantitative Trait Loci (QTL) has been a challenge for complex traits due to the use of populations with narrow genetic base. Most of QTL mapping studies were carried out from crosses made within the subspecies, either indica × indica or japonica × japonica. In this study we report advantages of using Multi-parent Advanced Generation Inter-Crosses global population, derived from a combination of eight indica and eight japonica elite parents, in QTL discovery for yield and grain quality traits. Genome-wide association study and interval mapping identified 38 and 34 QTLs whereas Bayesian networking detected 60 QTLs with 22 marker-marker associations, 32 trait-trait associations and 65 marker-trait associations. Notably, nine known QTLs/genes qPH₁/OsGA20ox2, qDF₃/OsMADS50, PL, QDg1, qGW-5b, grb7-2, qGL₃/GS3, Amy6/Wx gene and OsNAS3 were consistently identified by all approaches for nine traits whereas qDF₃/OsMADS50 was co-located for both yield and days-to-flowering traits on chromosome 3. Moreover, we identified a number of candidate QTLs in either one or two analyses but further validations will be needed. The results indicate that this new population has enabled identifications of significant QTLs and interactions for 16 traits through multiple approaches. Pyramided recombinant inbred lines provide a valuable source for integration into future breeding programs.

Genome-wide association study of the seed transmission rate of soybean mosaic virus and associated traits using two diverse population panels

KEY MESSAGE

Genome-wide association analyses identified candidates for genes involved in restricting virus movement into embryonic tissues, suppressing virus-induced seed coat mottling and preserving yield in soybean plants infected with soybean mosaic virus. Soybean mosaic virus (SMV) causes significant reductions in soybean yield and seed quality. Because seedborne infections can serve as primary sources of inoculum for SMV infections, resistance to SMV seed transmission provides a means to limit the impacts of SMV. In this study, two diverse population panels, Pop1 and Pop2, composed of 409 and 199 soybean plant introductions, respectively, were evaluated for SMV seed transmission rate, seed coat mottling, and seed yield from SMV-infected plants. The phenotypic data and genotypic data from the SoySNP50K dataset were analyzed using GAPIT and rrBLUP. For SMV seed transmission rate, a single locus was identified on chromosome 9 in Pop1. For SMV-induced seed coat mottling, loci were identified on chromosome 9 in Pop1 and on chromosome 3 in Pop2. For seed yield from SMV-infected plants, a single locus was identified on chromosome 3 in Pop2 that was within the map interval of a previously described quantitative trait locus for seed number. The high linkage disequilibrium regions surrounding the markers on chromosomes 3 and 9 contained a predicted nonsense-mediated RNA decay gene, multiple pectin methylesterase inhibitor genes (involved in restricting virus movement), two chalcone synthase genes, and a homolog of the yeast Rtf1 gene (involved in RNA-mediated transcriptional gene silencing). The results of this study provided additional insight into the genetic architecture of these three important traits, suggested candidate genes for downstream functional validation, and suggested that genomic prediction would outperform marker-assisted selection for two of the four trait-marker associations.

Reference Genome Anchoring of High-Density Markers for Association Mapping and Genomic Prediction in European Winter Wheat

In this study, we anchored genotyping-by-sequencing data to the International Wheat Genome Sequencing Consortium Reference Sequence v1.0 assembly to generate over 40,000 high quality single nucleotide polymorphism markers on a panel of 376 elite European winter wheat varieties released between 1946 and 2007. We compared association mapping and genomic prediction accuracy for a range of productivity traits with previous results based on lower density dominant DArT markers. The results demonstrate that the availability of RefSeq v1.0 supports higher precision trait mapping and provides the density of markers required to obtain accurate predictions of traits controlled by multiple small effect loci, including grain yield.

Mapping of resistance to corn borers in a MAGIC population of maize

BACKGROUND

Corn borers constitute an important pest of maize around the world; in particular Sesamia nonagrioides Lefèbvre, named Mediterranean corn borer (MCB), causes important losses in Southern Europe. Methods of selection can be combined with transgenic approaches to increase the efficiency and durability of the resistance to corn borers. Previous studies of the genetic factors involved in resistance to MCB have been carried out using bi-parental populations that have low resolution or using association inbred panels that have a low power to detect rare alleles. We developed a Multi-parent Advanced Generation InterCrosses (MAGIC) population to map with high resolution the genetic determinants of resistance to MCB.

RESULTS

We detected multiple single nucleotide polymorphisms (SNPs) of low effect associated with resistance to stalk tunneling by MCB. We dissected a wide region related to stalk tunneling in multiple studies into three smaller regions (at ~ 150, ~ 155, and ~ 165 Mb in chromosome 6) that closely overlap with regions associated with cell wall composition. We also detected regions associated with kernel resistance and agronomic traits, although the co-localization of significant regions between traits was very low. This indicates that it is possible the concurrent improvement of resistance and agronomic traits.

CONCLUSIONS

We developed a mapping population which allowed a finer dissection of the genetics of maize resistance to corn borers and a solid nomination of candidate genes based on functional information. The population, given its large variability, was also adequate to map multiple traits and study the relationship between them.

Multivariate Genome-Wide Association Analyses Reveal the Genetic Basis of Seed Fatty Acid Composition in Oat (Avena sativa L.)

Oat (Avena sativa L.) has a high concentration of oils, comprised primarily of healthful unsaturated oleic and linoleic fatty acids. To accelerate oat plant breeding efforts, we sought to identify loci associated with variation in fatty acid composition, defined as the types and quantities of fatty acids. We genotyped a panel of 500 oat cultivars with genotyping-by-sequencing and measured the concentrations of ten fatty acids in these oat cultivars grown in two environments. Measurements of individual fatty acids were highly correlated across samples, consistent with fatty acids participating in shared biosynthetic pathways. We leveraged these phenotypic correlations in two multivariate genome-wide association study (GWAS) approaches. In the first analysis, we fitted a multivariate linear mixed model for all ten fatty acids simultaneously while accounting for population structure and relatedness among cultivars. In the second, we performed a univariate association test for each principal component (PC) derived from a singular value decomposition of the phenotypic data matrix. To aid interpretation of results from the multivariate analyses, we also conducted univariate association tests for each trait. The multivariate mixed model approach yielded 148 genome-wide significant single-nucleotide polymorphisms (SNPs) at a 10% false-discovery rate, compared to 129 and 73 significant SNPs in the PC and univariate analyses, respectively. Thus, explicit modeling of the correlation structure between fatty acids in a multivariate framework enabled identification of loci associated with variation in seed fatty acid concentration that were not detected in the univariate analyses. Ultimately, a detailed characterization of the loci underlying fatty acid variation can be used to enhance the nutritional profile of oats through breeding.

The Future of Digital Sequence Information for Plant Genetic Resources for Food and Agriculture

The recent debates on the legal status of "digital sequence information" (DSI) at the international level could have extensive consequences for the future of agriculture and food security. A large majority of recent advances in biology, medicine, or agriculture were achieved by sharing and mining of freely accessible sequencing data. It is most probably because of the tremendous success of modern genomics and advances of synthetic biology that concerns were raised about possible fair and equitable ways of sharing data. The DSI concept is relatively new, and all concerned parties agreed upon the need for a clear definition. For example, the extent to which DSI understanding is limited only to genetic sequence data has to be clarified. In this paper, I focus on a subset of DSI essential to humankind: the DSI originating from plant genetic resources for food and agriculture (PGRFA). Two international agreements shape the conservation and use of plant genetic resources: the Convention on Biodiversity and the International Treaty for Plant Genetic Resources for Food and Agriculture. In an attempt to mobilize DSI users and producers involved in research, breeding, and conservation, I describe here how the increasing amount of genomic data, information, and studies interact with the existing legal framework at the global level. Using possible scenarios, I will emphasize the complexity of the issues surrounding DSI for PGRFA and propose potential ways forward for developing an inclusive governance and fair use of these genetic resources.

Epistasis Detection and Modeling for Genomic Selection in Cowpea (Vigna unguiculata L. Walp.)

Genetic architecture reflects the pattern of effects and interaction of genes underlying phenotypic variation. Most mapping and breeding approaches generally consider the additive part of variation but offer limited knowledge on the benefits of epistasis which explains in part the variation observed in traits. In this study, the cowpea multiparent advanced generation inter-cross (MAGIC) population was used to characterize the epistatic genetic architecture of flowering time, maturity, and seed size. In addition, consideration for epistatic genetic architecture in genomic-enabled breeding (GEB) was investigated using parametric, semi-parametric, and non-parametric genomic selection (GS) models. Our results showed that large and moderate effect-sized two-way epistatic interactions underlie the traits examined. Flowering time QTL colocalized with cowpea putative orthologs of Arabidopsis thaliana and Glycine max genes like PHYTOCLOCK1 (PCL1 [Vigun11g157600]) and PHYTOCHROME A (PHY A [Vigun01g205500]). Flowering time adaptation to long and short photoperiod was found to be controlled by distinct and common main and epistatic loci. Parametric and semi-parametric GS models outperformed non-parametric GS model, while using known quantitative trait nucleotide(s) (QTNs) as fixed effects improved prediction accuracy when traits were controlled by large effect loci. In general, our study demonstrated that prior understanding of the genetic architecture of a trait can help make informed decisions in GEB.

Genomic signatures of seed mass adaptation to global precipitation gradients in sorghum

Seed mass is a key component of adaptation in plants and a determinant of yield in crops. The climatic drivers and genomic basis of seed mass variation remain poorly understood. In the cereal crop Sorghum bicolor, globally-distributed landraces harbor abundant variation in seed mass, which is associated with precipitation in their agroclimatic zones of origin. This study aimed to test the hypothesis that diversifying selection across precipitation gradients, acting on ancestral cereal grain size regulators, underlies seed mass variation in global sorghum germplasm. We tested this hypothesis in a set of 1901 georeferenced and genotyped sorghum landraces, 100-seed mass from common gardens, and bioclimatic precipitation variables. As predicted, 100-seed mass in global germplasm varies significantly among botanical races and is correlated to proxies of the precipitation gradients. With general and mixed linear model genome-wide associations, we identified 29 and 56 of 100 a priori candidate seed size genes with polymorphisms in the top 1% of seed mass association, respectively. Eleven of these genes harbor polymorphisms associated with the precipitation gradient, including orthologs of genes that regulate seed size in other cereals. With FarmCPU, 13 significant SNPs were identified, including one at an a priori candidate gene. Finally, we identified eleven colocalized outlier SNPs associated with seed mass and precipitation that also carry signatures of selection based on F_ST scans and PCAdapt, which represents a significant enrichment. Our findings suggest that seed mass in sorghum was shaped by diversifying selection on drought stress, and can inform genomics-enabled breeding for climate-resilient cereals.