Development of a multiple-hybrid population for genome-wide association studies: theoretical consideration and genetic mapping of flowering traits in maize

Genome-wide association analysis of four yield-related traits using a maize (Zea mays L.) F1 population

Increasing the yield of maize F1 hybrid is one of the most important target for breeders. However, as a result of the genetic complexity and extremely low heritability, it is very difficult to directly dissect the genetic basis and molecular mechanisms of yield, and reports on genetic analysis of F1 hybrid yield are rare. Taking F1 hybrid as the research object and dividing the yield into different affect factors, this approach may be the best strategy for clarifying the genetic mechanism of yield. Therefore, in this study, a maize F1 population consisting of 300 hybrids with 17,652 single nucleotide polymorphisms (SNPs) markers was used for genome-wide association study (GWAS) to filtrate candidate genes associated with the four yield-related traits, i.e., kernel row number (KRN), kernel number per row (KNPR), ear tip-barrenness (ETB), and hundred kernel weight (HKW). Combined with the results of previous studies and functional annotation information of candidate genes, a total of six candidate genes were identified as being associated with the four traits, which were involved in plant growth and development, protein synthesis response, phytohormone biosynthesis and signal transduction. Our results improve the understanding of the genetic basis of the four yield-related traits and may be provide a new strategy for the genetic basis of maize yield.

Dual-Model GWAS Analysis and Genomic Selection of Maize Flowering Time-Related Traits

An appropriate flowering period is an important selection criterion in maize breeding. It plays a crucial role in the ecological adaptability of maize varieties. To explore the genetic basis of flowering time, GWAS and GS analyses were conducted using an associating panel consisting of 379 multi-parent DH lines. The DH population was phenotyped for days to tasseling (DTT), days to pollen-shedding (DTP), and days to silking (DTS) in different environments. The heritability was 82.75%, 86.09%, and 85.26% for DTT, DTP, and DTS, respectively. The GWAS analysis with the FarmCPU model identified 10 single-nucleotide polymorphisms (SNPs) distributed on chromosomes 3, 8, 9, and 10 that were significantly associated with flowering time-related traits. The GWAS analysis with the BLINK model identified seven SNPs distributed on chromosomes 1, 3, 8, 9, and 10 that were significantly associated with flowering time-related traits. Three SNPs 3_198946071, 9_146646966, and 9_152140631 showed a pleiotropic effect, indicating a significant genetic correlation between DTT, DTP, and DTS. A total of 24 candidate genes were detected. A relatively high prediction accuracy was achieved with 100 significantly associated SNPs detected from GWAS, and the optimal training population size was 70%. This study provides a better understanding of the genetic architecture of flowering time-related traits and provides an optimal strategy for GS.

Multi-locus genome-wide association studies reveal the dynamic genetic architecture of flowering time in chrysanthemum

KEY MESSAGE

The dynamic genetic architecture of flowering time in chrysanthemum was elucidated by GWAS. Thirty-six known genes and 14 candidate genes were identified around the stable QTNs and QEIs, among which ERF-1 was highlighted. Flowering time (FT) adaptation is one of the major breeding goals in chrysanthemum, a multipurpose ornamental plant. In order to reveal the dynamic genetic architecture of FT in chrysanthemum, phenotype investigation of ten FT-related traits was conducted on 169 entries in 2 environments. The broad-sense heritability of five non-conditional FT traits, i.e., budding (FBD), visible coloring (VC), early opening (EO), full-bloom (OF) and decay period (DP), ranged from 56.93 to 84.26%, which were higher than that of the five derived conditional FT traits (38.51-75.13%). The phenotypic variation coefficients of OF_EO and DP_OF were relatively large ranging from 30.59 to 36.17%. Based on 375,865 SNPs, the compressed variance component mixed linear model 3VmrMLM was applied for a multi-locus genome-wide association study (GWAS). As a result, 313 quantitative trait nucleotides (QTNs) were identified for the non-conditional FT traits in single-environment analysis, while 119 QTNs and 67 QTN-by-environment interactions (QEIs) were identified in multi-environment analysis. As for the conditional traits, 343 QTNs were detected in single-environment analysis, and 119 QTNs and 83 QEIs were identified in multi- environment analysis. Among the genes around stable QTNs and QEIs, 36 were orthologs of known FT genes in Arabidopsis and other plants; 14 candidates were mined by combining the transcriptomics data and functional annotation, including ERF-1, ACA10, and FOP1. Furthermore, the haplotype analysis of ERF-1 revealed six elite accessions with extreme FBD. Our findings contribute to the understanding of dynamic genetic architecture of FT and provide valuable resources for future chrysanthemum molecular breeding programs.

Advancements and Prospects of Genome-Wide Association Studies (GWAS) in Maize

Genome-wide association studies (GWAS) have emerged as a powerful tool for unraveling intricate genotype-phenotype association across various species. Maize (Zea mays L.), renowned for its extensive genetic diversity and rapid linkage disequilibrium (LD), stands as an exemplary candidate for GWAS. In maize, GWAS has made significant advancements by pinpointing numerous genetic loci and potential genes associated with complex traits, including responses to both abiotic and biotic stress. These discoveries hold the promise of enhancing adaptability and yield through effective breeding strategies. Nevertheless, the impact of environmental stress on crop growth and yield is evident in various agronomic traits. Therefore, understanding the complex genetic basis of these traits becomes paramount. This review delves into current and future prospectives aimed at yield, quality, and environmental stress resilience in maize and also addresses the challenges encountered during genomic selection and molecular breeding, all facilitated by the utilization of GWAS. Furthermore, the integration of omics, including genomics, transcriptomics, proteomics, metabolomics, epigenomics, and phenomics has enriched our understanding of intricate traits in maize, thereby enhancing environmental stress tolerance and boosting maize production. Collectively, these insights not only advance our understanding of the genetic mechanism regulating complex traits but also propel the utilization of marker-assisted selection in maize molecular breeding programs, where GWAS plays a pivotal role. Therefore, GWAS provides robust support for delving into the genetic mechanism underlying complex traits in maize and enhancing breeding strategies.

Genome-Wide Association Study for Maize Hybrid Performance in a Typical Breeder Population

Maize is one of the major crops that has demonstrated success in the utilization of heterosis. Developing high-yield hybrids is a crucial part of plant breeding to secure global food demand. In this study, we conducted a genome-wide association study (GWAS) for 10 agronomic traits using a typical breeder population comprised 442 single-cross hybrids by evaluating additive, dominance, and epistatic effects. A total of 49 significant single nucleotide polymorphisms (SNPs) and 69 significant pairs of epistasis were identified, explaining 26.2% to 64.3% of the phenotypic variation across the 10 traits. The enrichment of favorable genotypes is significantly correlated to the corresponding phenotype. In the confident region of the associated site, 532 protein-coding genes were discovered. Among these genes, the Zm00001d044211 candidate gene was found to negatively regulate starch synthesis and potentially impact yield. This typical breeding population provided a valuable resource for dissecting the genetic architecture of yield-related traits. We proposed a novel mating strategy to increase the GWAS efficiency without utilizing more resources. Finally, we analyzed the enrichment of favorable alleles in the Shaan A and Shaan B groups, as well as in each inbred line. Our breeding practice led to consistent results. Not only does this study demonstrate the feasibility of GWAS in F1 hybrid populations, it also provides a valuable basis for further molecular biology and breeding research.

Genome-wide association study of multiyear dynamic growth traits in hybrid Liriodendron identifies robust genetic loci associated with growth trajectories

The genetic factors underlying growth traits differ over time points or stages. However, most current studies of phenotypes at single time points do not capture all loci or explain the genetic differences underlying growth trajectories. Hybrid Liriodendron exhibits obvious heterosis and is widely cultivated, although its complex genetic mechanism underlying growth traits remains unknown. A genome-wide association study (GWAS) is an effective method for elucidating the genetic architecture by identifying genetic loci underlying complex quantitative traits. In the present study, using a GWAS, we identified robust loci associated with growth trajectories in hybrid Liriodendron populations. We selected 233 hybrid progenies derived from 25 crosses for resequencing, and measured their tree height (H) and diameter at breast height (DBH) for 11 consecutive years; 192 972 high-quality single nucleotide polymorphisms (SNPs) were obtained. The dynamics of the multiyear single-trait GWAS showed that year-specific SNPs predominated, and only five robust SNPs for DBH were identified in at least three different years. Multitrait GWAS analysis with model parameters as latent variables also revealed 62 SNPs for H and 52 for DBH associated with the growth trajectory, displaying different biomass accumulation patterns, among which four SNPs exerted pleiotropic effects. All identified SNPs also exhibited temporal variations in effect sizes and inheritance patterns potentially related to different growth and developmental stages. The haplotypes resulting from these significant SNPs might pyramid favorable loci, benefitting the selection of superior genotypes. The present study provides insights into the genetic architecture of dynamic growth traits and lays a basis for future molecular-assisted breeding.

Hybrid performance evaluation and genome-wide association analysis of root system architecture in a maize association population

KEY MESSAGE

The genetic architecture of RSA traits was dissected by GWAS and coexpression networks analysis in a maize association population. Root system architecture (RSA) is a crucial determinant of water and nutrient uptake efficiency in crops. However, the maize genetic architecture of RSA is still poorly understood due to the challenges in quantifying root traits and the lack of dense molecular markers. Here, an association mapping panel including 356 inbred lines were crossed with a common tester, Zheng58, and the test crosses were phenotyped for 12 RSA traits in three locations. We observed a 1.3 ~ sixfold phenotypic variation for measured RSA in the association panel. The association panel consisted of four subpopulations, non-stiff stalk (NSS) lines, stiff stalk (SS), tropical/subtropical (TST), and mixed. Zheng58 × TST has a 2.1% higher crown root number (CRN) and 8.6% less brace root number (BRN) than Zheng58 × NSS and Zheng58 × SS, respectively. Using a genome-wide association study (GWAS) with 1.25 million SNPs and correction for population structure, 191 significant SNPs were identified for root traits. Ninety (47%) of the significant SNPs showed positive allelic effects, and 101 (53%) showed negative effects. Each locus could explain 0.39% to 11.8% of phenotypic variation. By integrating GWAS results and comparing coexpression networks, 26 high-priority candidate genes were identified. Gene GRMZM2G377215, which belongs to the COBRA-like gene family, affected root growth and development. Gene GRMZM2G468657 encodes the aspartic proteinase nepenthesin-1, related to root development and N-deficient response. Collectively, our research provides progress in the genetic dissection of root system architecture. These findings present the further possibility for the genetic improvement of root traits in maize.

Epistasis Activation Contributes Substantially to Heterosis in Temperate by Tropical Maize Hybrids

Epistasis strongly affects the performance of superior maize hybrids. In this study, a multiple-hybrid population, consisting of three hybrid maize sets with varied interparental divergence, was generated by crossing 28 temperate and 23 tropical inbred lines with diverse genetic backgrounds. We obtained 1,154 tested hybrids. Among these tested hybrids, heterosis increased steadily as the heterotic genetic distance increased. Mid-parent heterosis was significantly higher in the temperate by tropical hybrids than in the temperate by temperate hybrids. Genome-wide prediction and association mapping was performed for grain weight per plant (GWPP) and days to silking (DTS) using 20K high-quality SNPs, showing that epistatic effects played a more prominent role than dominance effects in temperate by tropical maize hybrids. A total of 33 and 420 epistatic QTL were identified for GWPP and DTS, respectively, in the temperate by tropical hybrids. Protein-protein interaction network and gene-set enrichment analyses showed that epistatic genes were involved in protein interactions, which play an important role in photosynthesis, biological transcription pathways, and protein synthesis. We showed that the interaction of many minor-effect genes in the hybrids could activate the transcription activators of epistatic genes, resulting in a cascade of amplified yield heterosis. The multiple-hybrid population design enhanced our understanding of heterosis in maize, providing an insight into the acceleration of hybrid maize breeding by activating epistatic effects.

Nutritional improvement of cereal crops to combat hidden hunger during COVID-19 pandemic: Progress and prospects

COVID-19 has posed a severe challenge on food security by limiting access to food for the marginally placed population. While access to food is a challenge, access to nutritional food is a greater challenge to the population. The present-day foods are not sufficient to meet the nutritional requirements of the human body. In a pandemic condition, providing nutritious food to the population is imperative to ensure the health and well-being of humankind. Exploiting the existing biodiversity of crop species and deploying classical and modern tools to improve the nutritional potential of these species holds the key to addressing the above challenge. Breeding has been a classical tool of crop improvement that relied predominantly on genetic diversity. Collecting and conserving diverse germplasms and characterizing their diversity using molecular markers is essential to preserve diversity and use them in genetic improvement programs. These markers are also valuable for association mapping analyses to identify the genetic determinants of traits-of-interest in crop species. Association mapping identifies the quantitative trait loci (QTL) underlying the trait-of-interest by exploring marker-trait associations, and these QTLs can further be exploited for the genetic improvement of cultivated species through genomics-assisted breeding. Conventional breeding and genomics approaches are also being applied to develop biofortified cereal crops to reduce nutritional deficiencies in consumers. In this context, chapter explains the prerequisites for association mapping, population structure, genetic diversity, different approaches of performing association mapping to dissect nutritional traits, use the information for genomics-assisted breeding for nutrient-rich cereal crops, and application of genomics strategies in crop biofortification. These approaches will ensure food and nutrition security for all amidst the current COVID-19 crisis.

Genome-wide association mapping in maize: status and prospects

Genome-wide association study (GWAS) provides a robust and potent tool to retrieve complex phenotypic traits back to their underlying genetics. Maize is an excellent crop for performing GWAS due to diverse genetic variability, rapid decay of linkage disequilibrium, availability of distinct sub-populations and abundant SNP information. The application of GWAS in maize has resulted in successful identification of thousands of genomic regions associated with many abiotic and biotic stresses. Many agronomic and quality traits of maize are severely affected by such stresses and, significantly affecting its growth and productivity. To improve productivity of maize crop in countries like India which contribute only 2% to the world's total production in 2019-2020, it is essential to understand genetic complexity of underlying traits. Various DNA markers and trait associations have been revealed using conventional linkage mapping methods. However, it has achieved limited success in improving polygenic complex traits due to lower resolution of trait mapping. The present review explores the prospects of GWAS in improving yield, quality and stress tolerance in maize besides, strengths and challenges of using GWAS for molecular breeding and genomic selection. The information gathered will facilitate elucidation of genetic mechanisms of complex traits and improve efficiency of marker-assisted selection in maize breeding.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s13205-021-02799-4.

Large-Scale Analysis of Combining Ability and Heterosis for Development of Hybrid Maize Breeding Strategies Using Diverse Germplasm Resources

Understanding combining ability and heterosis among diverse maize germplasm resources is important for breeding hybrid maize (Zea mays L.). Using 28 temperate and 23 tropical maize inbreds that represent different ecotypes and worldwide diversity of maize germplasm, we first developed a large-scale multiple-hybrid population (MHP) with 724 hybrids, which could be divided into three subsets, 325 temperate diallel hybrids and 136 tropical diallel hybrids generated in Griffing IV, and 263 temperate by tropical hybrids generated in NCD II. All the parental lines and hybrids were evaluated for 11 traits in replicated tests across two locations and three years. Several widely used inbreds showed strong general combining ability (GCA), and their derived hybrids showed strong specific combining ability (SCA). Heterosis is a quantifiable, trait-dependent and environment-specific phenotype, and the response of parental lines and their hybrids to environments resulted in various levels of heterosis. For all the tested traits except plant height and hundred grain weight (HGW), NCD II (temperate × tropical) hybrids showed higher average heterosis than the temperate and tropical diallel hybrids, with higher hybrid performance for ear length, ear diameter, and HGW. Tropical maize germplasm can be used to improve the yield potential for temperate lines. Grain number per row and grain number per ear were two most important traits that determined yield heterosis, which can be used as direct selection criteria for yield heterosis. The hybrids from heterotic groups, Reid × SPT, Reid × LRC, SPT × PA, and Lancaster × LRC, contributed highly significant positive SCA effects and strong heterosis to yield-related traits, and the heterotic patterns identified in this study were potentially useful for commercial maize breeding. Heterosis was more significantly and positively correlated with SCA than GCA, indicating that SCA can be used in heterosis prediction to develop potential hybrids in commercial maize breeding. The results of the present study not only contribute to developing breeding strategies, but also improve targeted breeding efficiency by using both temperate and tropical maize to broaden genetic basis. Large sets of parental lines with available genotypic information can be shared and used in worldwide hybrid breeding programs through an open-source breeding strategy. Potential applications of the reported results in developing hybrid maize breeding strategies were also discussed.

Enhancing Genetic Gain through Genomic Selection: From Livestock to Plants

Although long-term genetic gain has been achieved through increasing use of modern breeding methods and technologies, the rate of genetic gain needs to be accelerated to meet humanity's demand for agricultural products. In this regard, genomic selection (GS) has been considered most promising for genetic improvement of the complex traits controlled by many genes each with minor effects. Livestock scientists pioneered GS application largely due to livestock's significantly higher individual values and the greater reduction in generation interval that can be achieved in GS. Large-scale application of GS in plants can be achieved by refining field management to improve heritability estimation and prediction accuracy and developing optimum GS models with the consideration of genotype-by-environment interaction and non-additive effects, along with significant cost reduction. Moreover, it would be more effective to integrate GS with other breeding tools and platforms for accelerating the breeding process and thereby further enhancing genetic gain. In addition, establishing an open-source breeding network and developing transdisciplinary approaches would be essential in enhancing breeding efficiency for small- and medium-sized enterprises and agricultural research systems in developing countries. New strategies centered on GS for enhancing genetic gain need to be developed.

Additive and heterozygous (dis)advantage GWAS models reveal candidate genes involved in the genotypic variation of maize hybrids to Azospirillum brasilense

Maize genotypes can show different responsiveness to inoculation with Azospirillum brasilense and an intriguing issue is which genes of the plant are involved in the recognition and growth promotion by these Plant Growth-Promoting Bacteria (PGPB). We conducted Genome-Wide Association Studies (GWAS) using additive and heterozygous (dis)advantage models to find candidate genes for root and shoot traits under nitrogen (N) stress and N stress plus A. brasilense. A total of 52,215 Single Nucleotide Polymorphism (SNP) markers were used for GWAS analyses. For the six root traits with significant inoculation effect, the GWAS analyses revealed 25 significant SNPs for the N stress plus A. brasilense treatment, in which only two were overlapped with the 22 found for N stress only. Most were found by the heterozygous (dis)advantage model and were more related to exclusive gene ontology terms. Interestingly, the candidate genes around the significant SNPs found for the maize-A. brasilense association were involved in different functions previously described for PGPB in plants (e.g. signaling pathways of the plant's defense system and phytohormone biosynthesis). Our findings are a benchmark in the understanding of the genetic variation among maize hybrids for the association with A. brasilense and reveal the potential for further enhancement of maize through this association.

Genome-wide association study dissects yield components associated with low-phosphorus stress tolerance in maize

Phosphorus deficiency in soil is a worldwide constraint threatening maize production. Through a genome-wide association study, we identified molecular markers and associated candidate genes and molecular pathways for low-phosphorus stress tolerance. Phosphorus deficiency in soils will severely affect maize (Zea mays L.) growth and development, thus decreasing the final yield. Deciphering the genetic basis of yield-related traits can benefit our understanding of maize tolerance to low-phosphorus stress. However, considering that yield-related traits should be evaluated under field condition with large populations rather than under hydroponic condition at a single-plant level, searching for appropriate field experimental sites and target traits for low-phosphorus stress tolerance is still very challenging. In this study, a genome-wide association analysis using two natural populations was performed to detect candidate genes in response to low-phosphorus stress at two experimental sites representative of different climate and soil types. In total, 259 candidate genes were identified and these candidate genes are mainly involved in four major pathways: transcriptional regulation, reactive oxygen scavenging, hormone regulation, and remodeling of cell wall. Among these candidate genes, 98 showed differential expression by transcriptome data. Based on a haplotype analysis of grain number under phosphorus deficiency condition, the positive haplotypes with favorable alleles across five loci increased grain number by 42% than those without favorable alleles. For further verifying the feasibility of genomic selection for improving maize low-phosphorus tolerance, we also validated the predictive ability of five genomic selection methods and suggested that moderate-density SNPs were sufficient to make accurate predictions for low-phosphorus tolerance traits. All these results will facilitate elucidating genetic basis of maize tolerance to low-phosphorus stress and improving marker-assisted selection efficiency in breeding process.

Genetic Loci Governing Androgenic Capacity in Perennial Ryegrass (Lolium perenne L.)

Immature pollen can be induced to switch developmental pathways from gametogenesis to embryogenesis and subsequently regenerate into homozygous, diploid plants. Such androgenic production of doubled haploids is particularly useful for species where inbreeding is hampered by effective self-incompatibility systems. Therefore, increasing the generally low androgenic capacity of perennial ryegrass (Lolium perenne L.) germplasm would enable the efficient production of homozygous plant material, so that a more effective exploitation of heterosis through hybrid breeding schemes can be realized. Here, we present the results of a genome-wide association study in a heterozygous, multiparental population of perennial ryegrass (n = 391) segregating for androgenic capacity. Genotyping-by-sequencing was used to interrogate gene- dense genomic regions and revealed over 1,100 polymorphic sites. Between one and 10 quantitative trait loci (QTL) were identified for anther response, embryo and total plant production, green and albino plant production and regeneration. Most traits were under polygenic control, although a major QTL on linkage group 5 was associated with green plant regeneration. Distinct genetic factors seem to affect green and albino plant recovery. Two intriguing candidate genes, encoding chromatin binding domains of the developmental phase transition regulator, Polycomb Repressive Complex 2, were identified. Our results shed the first light on the molecular mechanisms behind perennial ryegrass microspore embryogenesis and enable marker-assisted introgression of androgenic capacity into recalcitrant germplasm of this forage crop of global significance.

QTL Mapping in Three Connected Populations Reveals a Set of Consensus Genomic Regions for Low Temperature Germination Ability in Zea mays L

Improving seed vigor in response to cold stress is an important breeding objective in maize that allows early sowing. Using two cold tolerant inbred lines 220 and P9-10 and two susceptible lines Y1518 and PH4CV, three connected F_2:3 populations were generated for detecting quantitative trait locus (QTL) related to seed low-temperature germination ability. At 10°C, two germination traits (emergence rate and germination index) were collected from a sand bed and three seedling traits (seedling root length, shoot length, and total length) were extracted from paper rolls. Significant correlations were found among all traits in all populations. Via single-population analysis, 43 QTL were detected with explained phenotypic variance of 0.62%∼39.44%. Seventeen QTL explained more than 10% phenotypic variance; of them sixteen (94.12%) inherited favorable alleles from the tolerant lines. After constructing a consensus map, three meta-QTL (mQTL) were identified to include at least two initial QTL from different populations. mQTL1-1 included seven initial QTL for both germination and seedling traits; with three explaining more than 30% phenotypic variance. mQTL2-1 and mQTL9-1 covered two to three initial QTL. The favorable alleles of the QTL within these three mQTL regions were all inherited from the tolerant line 220 and P9-10. These results provided a basis for cloning of genes underlying the mQTL regions to uncover the molecular mechanisms of maize cold tolerance during germination.

Genome-wide association study reveals genetic loci and candidate genes for average daily gain in Duroc pigs

Objective

Average daily gain (ADG) is an important target trait of pig breeding programs. We aimed to identify single nucleotide polymorphisms (SNPs) and genomic regions that are associated with ADG in the Duroc pig population.

Methods

We performed a genome-wide association study involving 390 Duroc boars and by using the PorcineSNP60K Beadchip and two linear models.

Results

After quality control, we detected 3,5971 SNPs, which included seven SNPs that are significantly associated with the ADG of pigs. We identified six quantitative trait loci (QTL) regions for ADG. These QTLs included four previously reported QTLs on Sus scrofa chromosome (SSC) 1, SSC5, SSC9, and SSC13, as well as two novel QTLs on SSC6 and SSC16. In addition, we selected six candidate genes (general transcription factor 3C polypeptide 5, high mobility group AT-hook 2, nicotinamide phosphoribosyltransferase, oligodendrocyte transcription factor 1, pleckstrin homology and RhoGEF domain containing G4B, and ENSSSCG00000031548) associated with ADG on the basis of their physiological roles and positional information. These candidate genes are involved in skeletal muscle cell differentiation, diet-induced obesity, and nervous system development.

Conclusion

This study contributes to the identification of the casual mutation that underlies QTLs associated with ADG and to future pig breeding programs based on marker-assisted selection. Further studies are needed to elucidate the role of the identified candidate genes in the physiological processes involved in ADG regulation.