Prediction and association mapping of agronomic traits in maize using multiple omic data

Vegetable Signatures Derived from Human Urinary Metabolomic Data in Controlled Feeding Studies

Examination of changes in urinary metabolomic profiles after vegetable ingestion may lead to new methods of assessing plant food intake. To this regard, we developed a proof-of-principle methodology to identify urinary metabolomic signatures for spinach, celery, and onion. Three feeding studies were conducted. In the first study, healthy individuals were fed with spinach, celery, onion, and no vegetables in four separate experiments with pooled urinary samples for metabolite discovery. The same protocol was used to validate the finding at the individual level in the second study and when feeding all three vegetables simultaneously in the third study. An LC-MS-based metabolomics approach was adopted to search for indicative metabolites from urine samples collected during multiple time periods before and after the meal. Consequently, a total of 1, 9, and 3 nonoverlapping urinary metabolites were associated with the intake of spinach, celery, and onion, respectively. The PCA signature of these metabolites followed a similar "time cycle" pattern, which maximized at approximately 2-4 h after intake. In addition, the metabolite profiles for the same vegetable were consistent across samples, regardless of whether it was consumed individually or in combination. The developed methodology along with the identified urinary metabolomic signatures were potential tools for assessing plant food intake.

Modeling copy number variation in the genomic prediction of maize hybrids

Our study indicates that copy variants may play an essential role in the phenotypic variation of complex traits in maize hybrids. Moreover, predicting hybrid phenotypes by combining additive-dominance effects with copy variants has the potential to be a viable predictive model. Non-additive effects resulting from the actions of multiple loci may influence trait variation in single-cross hybrids. In addition, complementation of allelic variation could be a valuable contributor to hybrid genetic variation, especially when crossing inbred lines with higher contents of copy gains. With this in mind, we aimed (1) to study the association between copy number variation (CNV) and hybrid phenotype, and (2) to compare the predictive ability (PA) of additive and additive-dominance genomic best linear unbiased prediction model when combined with the effects of CNV in two datasets of maize hybrids (USP and HELIX). In the USP dataset, we observed a significant negative phenotypic correlation of low magnitude between copy number loss and plant height, revealing a tendency that more copy losses lead to lower plants. In the same set, when CNV was combined with the additive plus dominance effects, the PA significantly increased only for plant height under low nitrogen. In this case, CNV effects explicitly capture relatedness between individuals and add extra information to the model. In the HELIX dataset, we observed a pronounced difference in PA between additive (0.50) and additive-dominance (0.71) models for predicting grain yield, suggesting a significant contribution of dominance. We conclude that copy variants may play an essential role in the phenotypic variation of complex traits in maize hybrids, although the inclusion of CNVs into datasets does not return significant gains concerning PA.

Construction of a high-density genetic map and its application for leaf shape QTL mapping in poplar

High-quality and dense genetic maps were constructed, and leaf shape variation was dissected by QTL mapping in poplar. Species in the genus Populus, also known as poplars, are important woody species and considered model plants for perennial trees. High-density genetic maps are valuable genomic resources for population genetics. Here, we generated a high-quality and dense genetic map for an F1 poplar population using high-throughput NGS-based genotyping. A total of 92,097 high-quality SNP markers were developed by stringent filtering and identification. In total, 889 and 1650 SNPs formed the female and male genetic maps, respectively. To test the application of the genetic maps, QTL mapping of leaf shape was conducted for this F1 population. A total of nine parameters were scored for leaf shape variation in three different environments. Combining genetic maps and measurements of the nine leaf shape parameters, we mapped a total of 42 significant QTLs. The highest LOD score of all QTLs was 9.2, and that QTL explained the most (15.13%) trait variation. A total of nine QTLs could be detected in at least two environments, and they were located in two genomic regions. Within these two QTL regions, some candidate genes for regulating leaf shape were predicted through functional annotation. The successful mapping of leaf shape QTLs demonstrated the utility of our genetic maps. According to the performance of this study, we were able to provide high-quality and dense genetic maps and dissect the leaf shape variation in poplar.

Genomic selection using principal component regression

Many statistical methods are available for genomic selection (GS) through which genetic values of quantitative traits are predicted for plants and animals using whole-genome SNP data. A large number of predictors with much fewer subjects become a major computational challenge in GS. Principal components regression (PCR) and its derivative, i.e., partial least squares regression (PLSR), provide a solution through dimensionality reduction. In this study, we show that PCR can perform better than PLSR in cross validation. PCR often requires extracting more components to achieve the maximum predictive ability than PLSR and thus may be associated with a higher computational cost. However, application of the HAT method (a strategy of describing the relationship between the fitted and observed response variables with a hat matrix) to PCR circumvents conventional cross validation in testing predictive ability, resulting in substantially improved computational efficiency over PLSR where cross validation is mandatory. Advantages of PCR over PLSR are illustrated with a simulated trait of a hypothetical population and four agronomical traits of a rice population. The benefit of using PCR in genomic selection is further demonstrated in an effort to predict 1000 metabolomic traits and 24,973 transcriptomic traits in the same rice population.

Single-Locus and Multi-Locus Genome-Wide Association Studies in the Genetic Dissection of Fiber Quality Traits in Upland Cotton (Gossypium hirsutum L.)

A major breeding target in Upland cotton (Gossypium hirsutum L.) is to improve the fiber quality. To address this issue, 169 diverse accessions, genotyped by 53,848 high-quality single-nucleotide polymorphisms (SNPs) and phenotyped in four environments, were used to conduct genome-wide association studies (GWASs) for fiber quality traits using three single-locus and three multi-locus models. As a result, 342 quantitative trait nucleotides (QTNs) controlling fiber quality traits were detected. Of the 342 QTNs, 84 were simultaneously detected in at least two environments or by at least two models, which include 29 for fiber length, 22 for fiber strength, 11 for fiber micronaire, 12 for fiber uniformity, and 10 for fiber elongation. Meanwhile, nine QTNs with 10% greater sizes (R²) were simultaneously detected in at least two environments and between single- and multi-locus models, which include TM80185 (D13) for fiber length, TM1386 (A1) and TM14462 (A6) for fiber strength, TM18616 (A7), TM54735 (D3), and TM79518 (D12) for fiber micronaire, TM77489 (D12) and TM81448 (D13) for fiber uniformity, and TM47772 (D1) for fiber elongation. This indicates the possibility of marker-assisted selection in future breeding programs. Among 455 genes within the linkage disequilibrium regions of the nine QTNs, 113 are potential candidate genes and four are promising candidate genes. These findings reveal the genetic control underlying fiber quality traits and provide insights into possible genetic improvements in Upland cotton fiber quality.

Genome-Wide Association Mapping of Starch Pasting Properties in Maize Using Single-Locus and Multi-Locus Models

Maize starch plays a critical role in food processing and industrial application. The pasting properties, the most important starch characteristics, have enormous influence on fabrication property, flavor characteristics, storage, cooking, and baking. Understanding the genetic basis of starch pasting properties will be beneficial for manipulation of starch properties for a given purpose. Genome-wide association studies (GWAS) are becoming a powerful tool for dissecting the complex traits. Here, we carried out GWAS for seven pasting properties of maize starch with a panel of 230 inbred lines and 145,232 SNPs using one single-locus method, genome-wide efficient mixed model association (GEMMA), and three multi-locus methods, FASTmrEMMA, FarmCPU, and LASSO. We totally identified 60 quantitative trait nucleotides (QTNs) for starch pasting properties with these four GWAS methods. FASTmrEMMA detected the most QTNs (29), followed by FarmCPU (19) and LASSO (12), GEMMA detected the least QTNs (7). Of these QTNs, seven QTNs were identified by more than one method simultaneously. We further investigated locations of these significantly associated QTNs for possible candidate genes. These candidate genes and significant QTNs provide the guidance for further understanding of molecular mechanisms of starch pasting properties. We also compared the statistical powers and Type I errors of the four GWAS methods using Monte Carlo simulations. The results suggest that the multi-locus method is more powerful than the single-locus method and a combination of these multi-locus methods could help improve the detection power of GWAS.

Metabolome-wide association studies for agronomic traits of rice

Identification of trait-associated metabolites will advance the knowledge and understanding of the biosynthetic and catabolic pathways that are relevant to the complex traits of interest. In the past, the association between metabolites (treated as quantitative traits) and genetic variants (e.g., SNPs) has been extensively studied using metabolomic quantitative trait locus (mQTL) mapping. Nevertheless, the research on the association between metabolites with agronomic traits has been inadequate. In practice, the regular approaches for QTL mapping analysis may be adopted for metabolites-phenotypes association analysis due to the similarity in data structure of these two types of researches. In the study, we compared four regular QTL mapping approaches, i.e., simple linear regression (LR), linear mixed model (LMM), Bayesian analysis with spike-slab priors (Bayes B) and least absolute shrinkage and selection operator (LASSO), by testing their performances on the analysis of metabolome-phenotype associations. Simulation studies showed that LASSO had the higher power and lower false positive rate than the other three methods. We investigated the associations of 839 metobolites with five agronomic traits in a collection of 533 rice varieties. The results implied that a total of 25 metabolites were significantly associated with five agronomic traits. Literature search and bioinformatics analysis indicated that the identified 25 metabolites are significantly involved in some growth and development processes potentially related to agronomic traits. We also explored the predictability of agronomic traits based on the 839 metabolites through cross-validation, which showed that metabolomic prediction was efficient and its application in plant breeding has been justified.