A combination of linkage mapping and GWAS brings new elements on the genetic basis of yield-related traits in maize across multiple environments

Comparative transcriptome analysis identified candidate genes associated with kernel row number in maize

Kernel row number (KRN) is a crucial trait in maize that has a high impact on yield. However, KRN is a typical quantitative trait with only a few genes being verified so far. Here, two maize inbred lines with contrasting KRN were used to perform transcriptome analysis at five early ear developmental stages. Pairwise differential gene expression analyses were performed, and a total of 11,897 line-specific differentially expressed genes (DEGs) were detected between the two lines across the five development stages. Clustering analysis of line-specific DEGs revealed that the trends of gene expression changed significantly in the five stages, thus the five stages were further divided into two development phases: Phase I (V6-V8) and Phase II (V9-V10). Gene ontology enrichment analysis revealed that different transcriptional pathways were activated in the two phases. DEGs in Phase I were significantly enriched in morphogenesis and differentiation processes and hormone regulation. Of the 5,850 line-specific DEGs in Phase I, 2,132 genes were in known quantitative trait loci (QTLs) or flanking regions of quantitative trait nucleotides (QTNs), of which 92 were repeatedly detected in QTLs where QTNs also exist. The 92 high-probability candidate genes included development-related transcription factors (SBP-box and AP2/EREBP TFs) as well as genes involved in hormone homeostasis and signaling. Our study provides genetic resources for the elucidation of the molecular mechanisms of KRN development and reference for the cloning of candidate genes.

Combined Genome-Wide Association Study and Linkage Analysis for Mining Candidate Genes for the Kernel Row Number in Maize (Zea mays L.)

Kernel row number (KRN) is one of the key traits that significantly affect maize yield and productivity. Therefore, investigating the candidate genes and their functions in regulating KRN provides a theoretical basis and practical direction for genetic improvement in maize breeding, which is vital for increasing maize yield and understanding domestication. In this study, three recombinant inbred line (RIL) populations were developed using the parental lines AN20, YML1218, CM395, and Ye107, resulting in a multiparent population comprising a total of 490 F9 RILs. Phenotypic evaluation of the RILs for KRN was performed in three distinct environments. The heritability estimates of the RILs ranged from 81.40% to 84.16%. Genotyping-by-sequencing (GBS) of RILs identified 569,529 high-quality single nucleotide polymorphisms (SNPs). Combined genome-wide association study (GWAS) and linkage analyses revealed 120 SNPs and 22 quantitative trait loci (QTLs) which were significantly associated with KRN in maize. Furthermore, two novel candidate genes, Zm00001d042733 and Zm00001d042735, regulating KRN in maize were identified, which were located in close proximity to the significant SNP3-178,487,003 and overlapping the interval of QTL qKRN3-1. Zm00001d042733 encodes ubiquitin carboxyl-terminal hydrolase and Zm00001d042735 encodes the Arabidopsis Tóxicos en Levadura family of proteins. This study identified novel candidate loci and established a theoretical foundation for further functional validation of candidate genes. These findings deepen our comprehension of the genetic mechanisms that underpin KRN and offer potential applications of KRN-related strategies in developing maize varieties with higher yield.

Genetic architecture of kernel-related traits in sweet and waxy maize revealed by genome-wide association analysis

Introduction

Maize (Zea mays L.) is one of the most important crops worldwide, the kernel size-related traits are the major components of maize grain yield.

Methods

To dissect the genetic architecture of four kernel-related traits of 100-kernel weight, kernel length, kernel width, and kernel diameter, a genome-wide association study (GWAS) was conducted in the waxy and sweet maize panel comprising of 447 maize inbred lines re-sequenced at the 5× coverage depth. GWAS analysis was carried out with the mixed linear model using 1,684,029 high-quality SNP markers.

Results

In total, 49 SNPs significantly associated with the four kernel-related traits were identified, including 46 SNPs on chromosome 3, two SNPs on chromosome 4, and one SNP on chromosome 7. Haplotype regression analysis identified 338 haplotypes that significantly affected these four kernel-related traits. Genomic selection (GS) results revealed that a set of 10,000 SNPs and a training population size of 30% are sufficient for the application of GS in waxy and sweet maize breeding for kernel weight and kernel size. Forty candidate genes associated with the four kernel-related traits were identified, including both Zm00001d000707 and Zm00001d044139 expressed in the kernel development tissues and stages with unknown functions.

Discussion

These significant SNPs and important haplotypes provide valuable information for developing functional markers for the implementation of marker-assisted selection in breeding. The molecular mechanism of Zm00001d000707 and Zm00001d044139 regulating these kernel-related traits needs to be investigated further.

Genome-wide association analysis of grain yield and Striga hermonthica and S. asiatica resistance in tropical and sub-tropical maize populations

BACKGROUND

Genetic improvement for Striga hermonthica (Sh) and S. asiatica (Sa) resistance is the most economical and effective control method to enhance the productivity of maize and other major cereal crops. Hence, identification of quantitative trait loci (QTL) associated with Striga resistance and economic traits will guide the pace and precision of resistance breeding in maize. The objective of this study was to undertake a genome-wide association analysis of grain yield and Sh and Sa resistance among tropical and sub-tropical maize populations to identify putative genetic markers and genes for resistance breeding. 126 maize genotypes were evaluated under controlled environment conditions using artificial infestation of Sh and Sa. The test genotypes were profiled for grain yield (GY), Striga emergence counts at 8 (SEC8) and 10 (SEC10) weeks after planting, and Striga damage rate scores at 8 (SDR8) and 10 (SDR10) weeks after planting. Population structure analysis and genome-wide association mapping were undertaken based on 16,000 single nucleotide polymorphism (SNP) markers.

RESULTS

A linkage disequilibrium (LD) analysis in 798,675 marker pairs revealed that 21.52% of pairs were in significant linkage (P < 0.001). Across the chromosomes, the LD between SNPs decayed below a critical level (r2 = 0.1) at a map distance of 0.19 Mbp. The genome-wide association study identified 50 significant loci associated with Sh resistance and 22 significant loci linked to Sa resistance, corresponding to 39 and 19 candidate genes, respectively.

CONCLUSION

The study found non-significant QTL associated with dual resistance to the two examined Striga species Some of the detected genes reportedly conditioned insect and pathogen resistance, plant cell development, variable senescence, and pollen fertility. The markers detected in the present study for Sa resistance were reported for the first time. The gene Zm00001eb219710 was pleiotropic, and conditioned GY and SEC10, while Zm00001eb165170 affected SDR8 and SDR10, and Zm00001eb112030 conditioned SDR8 and SDR10 associated with Sh resistance. The candidate genes may facilitate simultaneous selection for Sh and Sa resistance and grain yield in maize after further validation and introgression in breeding pipelines. Overall, we recommend breeding maize specifically for resistance to each Striga species using germplasm adapted to the endemic region of each parasite.

Quantitative trait locus analysis of gray leaf spot resistance in the maize IBM Syn10 DH population

KEY MESSAGE

The exploration and dissection of a set of QTLs and candidate genes for gray leaf spot disease resistance using two fully assembled parental genomes may help expedite maize resistance breeding. The fungal disease of maize known as gray leaf spot (GLS), caused by Cercospora zeae-maydis and Cercospora zeina, is a significant concern in China, Southern Africa, and the USA. Resistance to GLS is governed by multiple genes with an additive effect and is influenced by both genotype and environment. The most effective way to reduce the cost of production is to develop resistant hybrids. In this study, we utilized the IBM Syn 10 Doubled Haploid (IBM Syn10 DH) population to identify quantitative trait loci (QTLs) associated with resistance to gray leaf spot (GLS) in multiple locations. Analysis of seven distinct environments revealed a total of 58 QTLs, 49 of which formed 12 discrete clusters distributed across chromosomes 1, 2, 3, 4, 8 and 10. By comparing these findings with published research, we identified colocalized QTLs or GWAS loci within eleven clustering intervals. By integrating transcriptome data with genomic structural variations between parental individuals, we identified a total of 110 genes that exhibit both robust disparities in gene expression and structural alterations. Further analysis revealed 19 potential candidate genes encoding conserved resistance gene domains, including putative leucine-rich repeat receptors, NLP transcription factors, fucosyltransferases, and putative xyloglucan galactosyltransferases. Our results provide a valuable resource and linked loci for GLS marker resistance selection breeding in maize.

Combining genome-wide association study and linkage mapping in the genetic dissection of amylose content in maize (Zea mays L.)

KEY MESSAGE

Integrated linkage and association analysis revealed genetic basis across multiple environments. The genes Zm00001d003102 and Zm00001d015905 were further verified to influence amylose content using gene-based association study. Maize kernel amylose is an important source of human food and industrial raw material. However, the genetic basis underlying maize amylose content is still obscure. Herein, we used an intermated B73 × Mo17 (IBM) Syn10 doubled haploid population composed of 222 lines and a germplasm set including 305 inbred lines to uncover the genetic control for amylose content under four environments. Linkage mapping detected 16 unique QTL, among which four were individually repeatedly identified across multiple environments. Genome-wide association study revealed 17 significant (P = 2.24E-06) single-nucleotide polymorphisms, of which two (SYN19568 and PZE-105090500) were located in the intervals of the mapped QTL (qAC2 and qAC5-3), respectively. According to the two population co-localized loci, 20 genes were confirmed as the candidate genes for amylose content. Gene-based association analysis indicated that the variants in Zm00001d003102 (Beta-16-galactosyltransferase GALT29A) and Zm00001d015905 (Sugar transporter 4a) affected amylose content across multi-environment. Tissue expression analysis showed that the two genes were specifically highly expressed in the ear and stem, respectively, suggesting that they might participate in sugar transport from source to sink organs. Our study provides valuable genetic information for breeding maize varieties with high amylose.

A combination of QTL mapping and genome-wide association study revealed the key gene for husk number in maize

KEY MESSAGE

Two key genes Zm00001d021232 and Zm00001d048138 were identified by QTL mapping and GWAS. Additionally, they were verified to be significantly associated with maize husk number (HN) using gene-based association study. As a by-product of maize production, maize husk is an important industrial raw material. Husk layer number (HN) is an important trait that affects the yield of maize husk. However, the genetic mechanism underlying HN remains unclear. Herein, a total of 13 quantitative trait loci (QTL) controlling HN were identified in an IBM Syn 10 DH population across different locations. Among these, three QTL were individually repeatedly detected in at least two environments. Meanwhile, 26 unique single nucleotide polymorphisms (SNPs) were detected to be significantly (p < 2.15 × 10-6) associated with HN in an association pool. Of these SNPs, three were simultaneously detected across multiple environments or environments and best linear unbiased prediction (BLUP). We focused on these environment-stable and population-common genetic loci for excavating the candidate genes responsible for maize HN. Finally, 173 initial candidate genes were identified, of which 22 were involved in both multicellular organism development and single-multicellular organism process and thus confirmed as the candidate genes for HN. Gene-based association analyses revealed that the variants in four genes were significantly (p < 0.01/N) correlated with HN, of which Zm00001d021232 and Zm00001d048138 were highly expressed in husks and early developing ears among different maize tissues. Our study contributes to the understanding of genetic and molecular mechanisms of maize husk yield and industrial development in the future.

Genetic dissection of ten photosynthesis-related traits based on InDel- and SNP-GWAS in soybean

KEY MESSAGE

A total of 416 InDels and 112 SNPs were significantly associated with soybean photosynthesis-related traits. GmIWS1 and GmCDC48 might be related to chlorophyll fluorescence and gas-exchange parameters, respectively. Photosynthesis is one of the main factors determining crop yield. A better understanding of the genetic architecture for photosynthesis is of great significance for soybean yield improvement. Our previous studies identified 5,410,112 single nucleotide polymorphisms (SNPs) from the resequencing data of 219 natural soybean accessions. Here, we identified 634,106 insertions and deletions (InDels) from these 219 accessions and used these InDel variations to perform principal component and linkage disequilibrium analysis of this population. The genome-wide association study (GWAS) were conducted on six chlorophyll fluorescence parameters (chlorophyll content, light energy absorbed per reaction center, quantum yield for electron transport, probability that a trapped exciton moves an electron into the electron transport chain beyond primary quinone acceptor, maximum quantum yield of photosystem II primary photochemistry in the dark-adapted state, performance index on absorption basis) and four gas-exchange parameters (intercellular carbon dioxide concentration, stomatal conductance, net photosynthesis rate, transpiration rate) and revealed 416 significant InDels and 112 significant SNPs. Based on GWAS results, GmIWS1 (encoding a transcription elongation factor) and GmCDC48 (encoding a cell division cycle protein) with the highest expression in the mapping region were determined as the candidate genes responsible for chlorophyll fluorescence and gas-exchange parameters, respectively. Further identification of favorable haplotypes with higher photosynthesis, seed weight and seed yield were carried out for GmIWS1 and GmCDC48. Overall, this study revealed the natural variations and candidate genes underlying the photosynthesis-related traits based on abundant phenotypic and genetic data, providing valuable insights into the genetic mechanisms controlling photosynthesis and yield in soybean.

Multi-omics analysis reveals novel loci and a candidate regulatory gene of unsaturated fatty acids in soybean (Glycine max (L.) Merr)

BACKGROUND

Soybean is a major oil crop; the nutritional components of soybean oil are mainly controlled by unsaturated fatty acids (FA). Unsaturated FAs mainly include oleic acid (OA, 18:1), linoleic acid (LLA, 18:2), and linolenic acid (LNA, 18:3). The genetic architecture of unsaturated FAs in soybean seeds has not been fully elucidated, although many independent studies have been conducted. A 3 V multi-locus random single nucleotide polymorphism (SNP)-effect mixed linear model (3VmrMLM) was established to identify quantitative trait loci (QTLs) and QTL-by-environment interactions (QEIs) for complex traits.

RESULTS

In this study, 194 soybean accessions with 36,981 SNPs were calculated using the 3VmrMLM model. As a result, 94 quantitative trait nucleotides (QTNs) and 19 QEIs were detected using single-environment (QTN) and multi-environment (QEI) methods. Three significant QEIs, namely rs4633292, rs39216169, and rs14264702, overlapped with a significant single-environment QTN.

CONCLUSIONS

For QTNs and QEIs, further haplotype analysis of candidate genes revealed that the Glyma.03G040400 and Glyma.17G236700 genes were beneficial haplotypes that may be associated with unsaturated FAs. This result provides ideas for the identification of soybean lipid-related genes and provides insights for breeding high oil soybean varieties in the future.

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.

Genomic prediction in multi-environment trials in maize using statistical and machine learning methods

In the context of multi-environment trials (MET), genomic prediction is proposed as a tool that allows the prediction of the phenotype of single cross hybrids that were not tested in field trials. This approach saves time and costs compared to traditional breeding methods. Thus, this study aimed to evaluate the genomic prediction of single cross maize hybrids not tested in MET, grain yield and female flowering time. We also aimed to propose an application of machine learning methodologies in MET in the prediction of hybrids and compare their performance with Genomic best linear unbiased prediction (GBLUP) with non-additive effects. Our results highlight that both methodologies are efficient and can be used in maize breeding programs to accurately predict the performance of hybrids in specific environments. The best methodology is case-dependent, specifically, to explore the potential of GBLUP, it is important to perform accurate modeling of the variance components to optimize the prediction of new hybrids. On the other hand, machine learning methodologies can capture non-additive effects without making any assumptions at the outset of the model. Overall, predicting the performance of new hybrids that were not evaluated in any field trials was more challenging than predicting hybrids in sparse test designs.

GWAS and Meta-QTL Analysis of Yield-Related Ear Traits in Maize

Maize ear traits are an important component of yield, and the genetic basis of ear traits facilitates further yield improvement. In this study, a panel of 580 maize inbred lines were used as the study material, eight ear-related traits were measured through three years of planting, and whole genome sequencing was performed using the maize 40 K breeding chip based on genotyping by targeted sequencing (GBTS) technology. Five models were used to conduct a genome-wide association study (GWAS) on best linear unbiased estimate (BLUE) of ear traits to find the best model. The FarmCPU (Fixed and random model Circulating Probability Unification) model was the best model for this study; a total of 104 significant single nucleotide polymorphisms (SNPs) were detected, and 10 co-location SNPs were detected simultaneously in more than two environments. Through gene function annotation and prediction, a total of nine genes were identified as potentially associated with ear traits. Moreover, a total of 760 quantitative trait loci (QTL) associated with yield-related traits reported in 37 different articles were collected. Using the collected 760 QTL for meta-QTL analysis, a total of 41 MQTL (meta-QTL) associated with yield-related traits were identified, and 19 MQTL detected yield-related ear trait functional genes and candidate genes that have been reported in maize. Five significant SNPs detected by GWAS were located within these MQTL intervals, and another three significant SNPs were close to MQTL (less than 1 Mb). The results provide a theoretical reference for the analysis of the genetic basis of ear-related traits and the improvement of maize yield.

The Genetic Structures and Molecular Mechanisms Underlying Ear Traits in Maize (Zea mays L.)

Maize (Zea mays L.) is one of the world's staple food crops. In order to feed the growing world population, improving maize yield is a top priority for breeding programs. Ear traits are important determinants of maize yield, and are mostly quantitatively inherited. To date, many studies relating to the genetic and molecular dissection of ear traits have been performed; therefore, we explored the genetic loci of the ear traits that were previously discovered in the genome-wide association study (GWAS) and quantitative trait locus (QTL) mapping studies, and refined 153 QTL and 85 quantitative trait nucleotide (QTN) clusters. Next, we shortlisted 19 common intervals (CIs) that can be detected simultaneously by both QTL mapping and GWAS, and 40 CIs that have pleiotropic effects on ear traits. Further, we predicted the best possible candidate genes from 71 QTL and 25 QTN clusters that could be valuable for maize yield improvement.

Mapping and Functional Analysis of QTL for Kernel Number per Row in Tropical and Temperate-Tropical Introgression Lines of Maize (Zea mays L.)

Kernel number per row (KNR) is an essential component of maize (Zea mays L.) grain yield (GY), and understanding its genetic mechanism is crucial to improve GY. In this study, two F₇ recombinant inbred line (RIL) populations were created using a temperate-tropical introgression line TML418 and a tropical inbred line CML312 as female parents and a backbone maize inbred line Ye107 as the common male parent. Bi-parental quantitative trait locus (QTL) mapping and genome-wide association analysis (GWAS) were then performed on 399 lines of the two maize RIL populations for KNR in two different environments using 4118 validated single nucleotide polymorphism (SNP) markers. This study aimed to: (1) detect molecular markers and/or the genomic regions associated with KNR; (2) identify the candidate genes controlling KNR; and (3) analyze whether the candidate genes are useful in improving GY. The authors reported a total of 7 QTLs tightly linked to KNR through bi-parental QTL mapping and identified 21 SNPs significantly associated with KNR through GWAS. Among these, a highly confident locus qKNR7-1 was detected at two locations, Dehong and Baoshan, with both mapping approaches. At this locus, three novel candidate genes (Zm00001d022202, Zm00001d022168, Zm00001d022169) were identified to be associated with KNR. These candidate genes were primarily involved in the processes related to compound metabolism, biosynthesis, protein modification, degradation, and denaturation, all of which were related to the inflorescence development affecting KNR. These three candidate genes were not reported previously and are considered new candidate genes for KNR. The progeny of the hybrid Ye107 × TML418 exhibited strong heterosis for KNR, which the authors believe might be related to qKNR7-1. This study provides a theoretical foundation for future research on the genetic mechanism underlying KNR in maize and the use of heterotic patterns to develop high-yielding hybrids.

GWAS across multiple environments and WGCNA suggest the involvement of ZmARF23 in embryonic callus induction from immature maize embryos

Combined GWAS, WGCNA, and gene-based association studies identified the co-expression network and hub genes for maize EC induction. ZmARF23 bound to ZmSAUR15 promoter and regulated its expression, affecting EC induction. Embryonic callus (EC) induction in immature maize embryos shows high genotype dependence, which limits the application of genetic transformation in transgenic breeding and gene function elucidation in maize. Herein, we conducted a genome-wide association mapping (GWAS) for four EC induction-related traits, namely rate of embryonic callus induction (REC), increased callus diameter (ICD), ratio of shoot formation (RSF), and length of shoot (LS) across different environments. A total of 77 SNPs were significantly associated these traits under three environments and using the averages (across environments). Among these significant SNPs, five were simultaneously detected under multiple environments and 11 had respective phenotypic variation explained > 10%. A total of 257 genes were located in the linkage disequilibrium decay of these REC- and ICD-associated SNPs, of which 178 were responsive to EC induction. According to the expression values of the 178 genes, we performed a weighted gene co-expression network analysis (WGCNA) and revealed an EC induction-associated module and five hub genes. Hub gene-based association studies uncovered that the intragenic variations in GRMZM2G105473 and ZmARF23 influenced EC induction efficiency among different maize lines. Dual-luciferase reporter assay indicated that ZmARF23 bound to the promoter of a known causal gene (ZmSAUR15) for EC induction and positively regulated its expression on the transcription level. Our study will deepen the understanding of genetic and molecular mechanisms underlying EC induction and contribute to the use of genetic transformation in maize.

Mapping quantitative trait loci and predicting candidate genes for Striga resistance in maize using resistance donor line derived from Zea diploperennis

The parasitic weed, Striga is a major biological constraint to cereal production in sub-Saharan Africa (SSA) and threatens food and nutrition security. Two hundred and twenty-three (223) F_2:3 mapping population involving individuals derived from TZdEI 352 x TZEI 916 were phenotyped for four Striga-adaptive traits and genotyped using the Diversity Arrays Technology (DArT) to determine the genomic regions responsible for Striga resistance in maize. After removing distorted SNP markers, a genetic linkage map was constructed using 1,918 DArTseq markers which covered 2092.1 cM. Using the inclusive composite interval mapping method in IciMapping, twenty-three QTLs influencing Striga resistance traits were identified across four Striga-infested environments with five stable QTLs (qGY4, qSC2.1, qSC2.2, qSC5, and qSC6) detected in more than one environment. The variations explained by the QTLs ranged from 4.1% (qSD2.3) to 14.4% (qSC7.1). Six QTLs each with significant additive × environment interactions were also identified for grain yield and Striga damage. Gene annotation revealed candidate genes underlying the QTLs, including the gene models GRMZM2G077002 and GRMZM2G404973 which encode the GATA transcription factors, GRMZM2G178998 and GRMZM2G134073 encoding the NAC transcription factors, GRMZM2G053868 and GRMZM2G157068 which encode the nitrate transporter protein and GRMZM2G371033 encoding the SBP-transcription factor. These candidate genes play crucial roles in plant growth and developmental processes and defense functions. This study provides further insights into the genetic mechanisms of resistance to Striga parasitism in maize. The QTL detected in more than one environment would be useful for further fine-mapping and marker-assisted selection for the development of Striga resistant and high-yielding maize cultivars.

Combined linkage mapping and association analysis uncovers candidate genes for 25 leaf-related traits across three environments in maize

Combined linkage and association analysis revealed five co-localized genetic loci across multiple environments. The key gene Zm00001d026491 was further verified to influence leaf length by candidate gene association analysis. Leaf morphology and number determine the canopy structure and thus affect crop yield. Herein, the genetic basis and key genes for 25 leaf-related traits, including leaf lengths (LL), leaf widths (LW), and leaf areas (LA) of eight continuous leaves under the tassel, and the number of leaves above the primary ear (LAE), were dissected by using an association panel and a biparental population. Using an intermated B73 × Mo17 (IBM) Syn10 doubled haploid (DH) population, 290 quantitative trait loci (QTL) controlling these traits were detected across different locations, among which 115 QTL were individually repeatedly identified in at least two environments. Using the association panel, 165 unique significant single-nucleotide polymorphisms (SNPs) were associated with target traits (P < 2.15E-06), of which 35 were separately detected across multiple environments. In total, 42 pleiotropic QTL/SNPs (pQTL/SNPs) were responsible for at least two of the LL, LW, LA, and LAE traits across multiple environments. Combining the QTL mapping and association study, five unique SNPs were located within the confidence intervals of seven QTL, and 77 genes were identified based on the linkage disequilibrium regions of co-localized SNP loci. Gene-based association studies verified that the intragenic variants in the candidate gene Zm00001d026491 influenced LL of the third leaf counted from the top node. These findings will provide vital information to understanding the genetic basis of leaf-related traits and help to cultivate maize varieties with ideal plant architecture.

CRISPR-Cas technology opens a new era for the creation of novel maize germplasms

Maize (Zea mays) is one of the most important food crops in the world with the greatest global production, and contributes to satiating the demands for human food, animal feed, and biofuels. With population growth and deteriorating environment, efficient and innovative breeding strategies to develop maize varieties with high yield and stress resistance are urgently needed to augment global food security and sustainable agriculture. CRISPR-Cas-mediated genome-editing technology (clustered regularly interspaced short palindromic repeats (CRISPR)-Cas (CRISPR-associated)) has emerged as an effective and powerful tool for plant science and crop improvement, and is likely to accelerate crop breeding in ways dissimilar to crossbreeding and transgenic technologies. In this review, we summarize the current applications and prospects of CRISPR-Cas technology in maize gene-function studies and the generation of new germplasm for increased yield, specialty corns, plant architecture, stress response, haploid induction, and male sterility. Optimization of gene editing and genetic transformation systems for maize is also briefly reviewed. Lastly, the challenges and new opportunities that arise with the use of the CRISPR-Cas technology for maize genetic improvement are discussed.

Integrating GWAS, linkage mapping and gene expression analyses reveal the genetic control of first branch height in Brassica napus L

Rapeseed (Brassica napus L.) is a crucial oil crop cultivated worldwide. First branch height, an essential component of rapeseed plant architecture, has an important effect on yield and mechanized harvesting; however, the underlying genetic mechanism remains unclear. In this study, based on the 60K single nucleotide polymorphism array and a recombinant inbred lines population derived from M083 and 888-5, a total of 19 QTLs were detected in five environments, distributed on linkage groups A02, A09, A10, C06, and C07, which explained phenotypic variation ranging from 4.87 to 29.87%. Furthermore, 26 significant SNPs were discovered on Chr.A02 by genome-wide association study in a diversity panel of 324 re-sequencing accessions. The major QTL of the first branch height trait was co-located on Chr.A02 by integrating linkage mapping and association mapping. Eleven candidate genes were screened via allelic variation analysis, inter-subgenomic synteny analysis, and differential expression of genes in parental shoot apical meristem tissues. Among these genes, BnaA02g13010D, which encodes a TCP transcription factor, was confirmed as the target gene according to gene function annotation, haplotype analysis, and full-length gene sequencing, which revealed that TATA insertion/deletion in the promoter region was closely linked to significantly phenotypic differences BnaA02.TCP1 ^M083 overexpression resulted in decreased branch height and increased branch number in Arabidopsis. These results provide a genetic basis for first branch height and the ideal architecture of B. napus.

Genome-wide association study of agronomic traits related to nitrogen use efficiency in wheat

GWAS identified 347 QTLs associated with eight traits related to nitrogen use efficiency in a 389-count wheat panel. Four novel candidate transcription factor genes were verified using qRT-PCR. Nitrogen is an essential nutrient for plants that determines crop yield. Improving nitrogen use efficiency (NUE) should considerably increase wheat yield and reduce the use of nitrogen fertilisers. However, knowledge on the genetic basis of NUE during wheat maturity is limited. In this study, a diversity panel incorporating 389 wheat accessions was phenotyped for eight NUE-related agronomic traits across five different environments. A total of 347 quantitative trait loci (QTLs) for low nitrogen tolerance indices (ratio of agronomic characters under low and high nitrogen conditions) were identified through a genome-wide association study utilising 397,384 single nucleotide polymorphisms (SNPs) within the MLM (Q + K) model, including 11 stable QTLs. Furthermore, 69 candidate genes were predicted for low nitrogen tolerance indices of best linear unbiased predictions values of the eight studied agronomic traits, and four novel candidate transcription factors (TraesCS5A02G237500 for qFsnR5A.2, TraesCS5B02G384500 and TraesCS5B02G384600 for qSLR5B.1, and TraesCS3B02G068800 for qTKWR3B.1) showed differing expression patterns in contrasting low-nitrogen-tolerant wheat genotypes. Moreover, the number of favourable marker alleles calculated using NUE that were significantly related to SNP in accessions decreased over the decades, indicating a decline in the NUE of the 389 wheat varieties. These findings denote promising NUE markers that could be useful in breeding high-NUE wheat varieties, and the candidate genes could further detail the NUE-related regulation network in wheat.

An Integration of Linkage Mapping and GWAS Reveals the Key Genes for Ear Shank Length in Maize

Ear shank length (ESL) has significant effects on grain yield and kernel dehydration rate in maize. Herein, linkage mapping and genome-wide association study were combined to reveal the genetic architecture of maize ESL. Sixteen quantitative trait loci (QTL) were identified in the segregation population, among which five were repeatedly detected across multiple environments. Meanwhile, 23 single nucleotide polymorphisms were associated with the ESL in the association panel, of which four were located in the QTL identified by linkage mapping and were designated as the population-common loci. A total of 42 genes residing in the linkage disequilibrium regions of these common variants and 12 of them were responsive to ear shank elongation. Of the 12 genes, five encode leucine-rich repeat receptor-like protein kinases, proline-rich proteins, and cyclin11, respectively, which were previously shown to regulate cell division, expansion, and elongation. Gene-based association analyses revealed that the variant located in Cyclin11 promoter affected the ESL among different lines. Cyclin11 showed the highest expression in the ear shank 15 days after silking among diverse tissues of maize, suggesting its role in modulating ESL. Our study contributes to the understanding of the genetic mechanism underlying maize ESL and genetic modification of maize dehydration rate and kernel yield.

Identification of quantitative trait loci for related traits of stalk lodging resistance using genome-wide association studies in maize (Zea mays L.)

BACKGROUND

Stalk lodging is one of the main factors affecting maize (Zea mays L.) yield and limiting mechanized harvesting. Developing maize varieties with high stalk lodging resistance requires exploring the genetic basis of lodging resistance-associated agronomic traits. Stalk strength is an important indicator to evaluate maize lodging and can be evaluated by measuring stalk rind penetrometer resistance (RPR) and stalk buckling strength (SBS). Along with morphological traits of the stalk for the third internodes length (TIL), fourth internode length (FIL), third internode diameter (TID), and the fourth internode diameter (FID) traits are associated with stalk lodging resistance.

RESULTS

In this study, a natural population containing 248 diverse maize inbred lines genotyped with 83,057 single nucleotide polymorphism (SNP) markers was used for genome-wide association study (GWAS) for six stalk lodging resistance-related traits. The heritability of all traits ranged from 0.59 to 0.72 in the association mapping panel. A total of 85 significant SNPs were identified for the association mapping panel using best linear unbiased prediction (BLUP) values of all traits. Additionally, five candidate genes were associated with stalk strength traits, which were either directly or indirectly associated with cell wall components.

CONCLUSIONS

These findings contribute to our understanding of the genetic basis of maize stalk lodging and provide valuable theoretical guidance for lodging resistance in maize breeding in the future.

Genome-wide meta-analysis of QTL for morphological related traits of flag leaf in bread wheat

Flag leaf is an important organ for photosynthesis of wheat plants, and a key factor affecting wheat yield. In this study, quantitative trait loci (QTL) for flag leaf morphological traits in wheat reported since 2010 were collected to investigate the genetic mechanism of these traits. Integration of 304 QTLs from various mapping populations into a high-density consensus map composed of various types of molecular markers as well as QTL meta-analysis discovered 55 meta-QTLs (MQTL) controlling morphological traits of flag leaves, of which 10 MQTLs were confirmed by GWAS. Four high-confidence MQTLs (MQTL-1, MQTL-11, MQTL-13, and MQTL-52) were screened out from 55 MQTLs, with an average confidence interval of 0.82 cM and a physical distance of 9.4 Mb, according to the definition of hcMQTL. Ten wheat orthologs from rice (7) and Arabidopsis (3) that regulated leaf angle, development and morphogenesis traits were identified in the hcMQTL region using comparative genomics, and were speculated to be potential candidate genes regulating flag leaf morphological traits in wheat. The results from this study provides valuable information for fine mapping and molecular markers assisted selection to improve morphological characters in wheat flag leaf.

Genome wide association analysis for yield related traits in maize

BACKGROUND

Understanding the genetic basis of yield related traits contributes to the improvement of grain yield in maize.

RESULTS

Using 291 excellent maize inbred lines as materials, six yield related traits of maize, including grain yield per plant (GYP), grain length (GL), grain width (GW), kernel number per row (KNR), 100 kernel weight (HKW) and tassel branch number (TBN) were investigated in Jinan, in 2017, 2018 and 2019. The average values of three environments were taken as the phenotypic data of yield related traits, and they were statistically analyzed. Based on 38,683 high-quality SNP markers in the whole genome of the association panel, the MLM with PCA model was used for genome-wide association analysis (GWAS) to obtain 59 significantly associated SNP sites. Moreover, 59 significantly associated SNPs (P < 0.0001) referring to GYP, GL, GW, KNR, HKW and TBN, of which 14 SNPs located in yield related QTLs/QTNs previously reported. A total of 66 candidate genes were identified based on the 59 significantly associated SNPs, of which 58 had functional annotation.

CONCLUSIONS

Using genome-wide association analysis strategy to identify genetic loci related to maize yield, a total of 59 significantly associated SNP were detected. Those results aid in our understanding of the genetic architecture of maize yield and provide useful SNPs for genetic improvement of maize.

Combining datasets for maize root seedling traits increases the power of GWAS and genomic prediction accuracies

The identification of genomic regions associated with root traits and the genomic prediction of untested genotypes can increase the rate of genetic gain in maize breeding programs targeting roots traits. Here, we combined two maize association panels with different genetic backgrounds to identify single nucleotide polymorphisms (SNPs) associated with root traits, and used a genome-wide association study (GWAS) and to assess the potential of genomic prediction for these traits in maize. For this, we evaluated 377 lines from the Ames panel and 302 from the Backcrossed Germplasm Enhancement of Maize (BGEM) panel in a combined panel of 679 lines. The lines were genotyped with 232 460 SNPs, and four root traits were collected from 14-day-old seedlings. We identified 30 SNPs significantly associated with root traits in the combined panel, whereas only two and six SNPs were detected in the Ames and BGEM panels, respectively. Those 38 SNPs were in linkage disequilibrium with 35 candidate genes. In addition, we found higher prediction accuracy in the combined panel than in the Ames or BGEM panel. We conclude that combining association panels appears to be a useful strategy to identify candidate genes associated with root traits in maize and improve the efficiency of genomic prediction.

A Combination of a Genome-Wide Association Study and a Transcriptome Analysis Reveals circRNAs as New Regulators Involved in the Response to Salt Stress in Maize

Salinization seriously threatens the normal growth of maize, especially at the seedling stage. Recent studies have demonstrated that circular RNAs (circRNAs) play vital roles in the regulation of plant stress resistance. Here, we performed a genome-wide association study (GWAS) on the survival rate of 300 maize accessions under a salt stress treatment. A total of 5 trait-associated SNPs and 86 candidate genes were obtained by the GWAS. We performed RNA sequencing for 28 transcriptome libraries derived from 2 maize lines with contrasting salt tolerance under normal and salt treatment conditions. A total of 1217 highly expressed circRNAs were identified, of which 371 were responsive to a salt treatment. Using PCR and Sanger sequencing, we verified the reliability of these differentially expressed circRNAs. An integration of the GWAS and RNA-Seq analyses uncovered two differentially expressed hub genes (Zm00001eb013650 and Zm00001eb198930), which were regulated by four circRNAs. Based on these results, we constructed a regulation model of circRNA/miRNA/mRNA that mediated salt stress tolerance in maize. By conducting hub gene-based association analyses, we detected a favorable haplotype in Zm00001eb198930, which was responsible for high salt tolerance. These results help to clarify the regulatory relationship between circRNAs and their target genes as well as to develop salt-tolerant lines for maize breeding.

Genome-Wide Association Study Reveals the Genetic Basis of Kernel and Cob Moisture Changes in Maize at Physiological Maturity Stage

Low moisture content (MC) and high dehydration rate (DR) at physiological maturity affect grain mechanical harvest, transport, and storage. In this study, we used an association panel composed of 241 maize inbred lines to analyze ear moisture changes at physiological maturity stage. A genome-wide association study revealed nine significant SNPs and 91 candidate genes. One SNP (SYN38588) was repeatedly detected for two traits, and 15 candidate genes were scanned in the linkage disequilibrium regions of this SNP. Of these, genes Zm00001d020615 and Zm00001d020623 were individually annotated as a polygalacturonase (PG) and a copper transporter 5.1 (COPT5.1), respectively. Candidate gene association analysis showed that three SNPs located in the exons of Zm00001d020615 were significantly associated with the dehydration rate, and AATTAA was determined as the superior haplotype. All these findings suggested that Zm00001d020615 was a key gene affecting moisture changes of maize at the physiological maturity stage. These results have demonstrated the genetic basis of ear moisture changes in maize and indicated a superior haplotype for cultivating maize varieties with low moisture content and high dehydration rates.

Combined QTL mapping and RNA-Seq pro-filing reveal candidate genes related to low-temperature tolerance in maize

Maize (Zea mays L.) is the most important food crop in the world, with significant acreage and production across the globe. However, it is affected by low temperatures throughout its growth process, especially during germination. Therefore, it is important to identify more QTLs or genes associated with germination under low-temperature conditions. For the QTL analysis of traits related to low-temperature germination, we used a high-res genetic map of 213 lines of the intermated B73 × Mo17 (IBM) Syn10 doubled haploid (DH) population, which had 6618 bin markers. We detected 28 QTLs of eight phenotypic characteristics associated with low-temperature germination, while they explained the phenotypic contribution rate of 5.4 ~ 13.34%. Additionally, 14 overlapping QTLs produced six QTL clusters on every chromosome, except for 8 and 10. RNA-Seq found six genes related to low-temperature tolerance in these QTLs, while qRT-PCR analysis demonstrated that the expression trends of the Zm00001d045568 gene in the LT_BvsLT_M group and the CK_BvsCK_M group were highly significantly different at all four-time points (P < 0.01), and encoded the RING zinc finger protein. It was located on qRTL9-2 and qRSVI9-1 and is related to the total length and simple vitality index. These results provided potential candidate genes for further gene cloning and improving the low-temperature tolerance of maize.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-022-01297-6.

Association mapping uncovers maize ZmbZIP107 regulating root system architecture and lead absorption under lead stress

Lead (Pb) is a highly toxic contaminant to living organisms and the environment. Excessive Pb in soils affects crop yield and quality, thus threatening human health via the food chain. Herein, we investigated Pb tolerance among a maize association panel using root bushiness (BSH) under Pb treatment as an indicator. Through a genome-wide association study of relative BSH, we identified four single nucleotide polymorphisms (SNPs) and 30 candidate genes associated with Pb tolerance in maize seedlings. Transcriptome analysis showed that four of the 30 genes were differentially responsive to Pb treatment between two maize lines with contrasting Pb tolerance. Among these, the ZmbZIP107 transcription factor was confirmed as the key gene controlling maize tolerance to Pb by using gene-based association studies. Two 5' UTR_variants in ZmbZIP107 affected its expression level and Pb tolerance among different maize lines. ZmbZIP107 protein was specifically targeted to the nucleus and ZmbZIP107 mRNA showed the highest expression in maize seedling roots among different tissues. Heterologous expression of ZmbZIP107 enhanced rice tolerance to Pb stress and decreased Pb absorption in the roots. Our study provided the basis for revelation of the molecular mechanism underlying Pb tolerance and contributed to cultivation of Pb-tolerant varieties in maize.

Identification of yield-related genes through genome-wide association: case study of weeping forsythia, an emerging medicinal crop

KEY MESSAGE

This study identified candidate genes related to fruit yield for an emerging medicinal crop, weeping forsythia.

BACKGROUND

The genetic basis of crop yield is an agricultural research hotspot. Identifying the genes related to yield traits is the key to increase the yield. Weeping forsythia is an emerging medicinal crop that currently lacks excellent varieties. The genes related to fruit yield in weeping forsythia have not been identified.

OBJECTIVE

Thus, we aimed to screen the candidate genes related to fruit yield of weeping forsythia by using genome-wide association analysis.

METHODS

Here, 60 samples from the same field and source of weeping forsythia were collected to identify its yield-related candidate genes. Association analysis was performed on the variant loci and the traits related to yield, i.e., fruit length, width, thickness, and weight.

RESULTS

Results from admixture, neighbor-joining, and kinship matrix analyses supported the non-significant genetic differentiation of these samples. Significant association was found between 2 variant loci and fruit length, 8 loci and fruit width, 24 loci and fruit thickness, and 13 loci and fruit weight. Further search on the 20 kb up/downstream of these variant loci revealed 1 gene related to fruit length, 16 genes related to fruit width, 12 genes related to fruit thickness, and 13 genes related to fruit weight. Among which, 4 genes, namely, WRKY transcription factor 35, salicylic acid-binding protein, auxin response factor 6, and alpha-mannosidase were highly related to the fruit development of weeping forsythia.

CONCLUSION

This study identify four candidate genes related to fruit development, which will provide useful information for the subsequent molecular-assisted and genetic breeding of weeping forsythia.

GWAS and WGCNA uncover hub genes controlling salt tolerance in maize (Zea mays L.) seedlings

KEYMESSAGE

Two hub genes GRMZM2G075104 and GRMZM2G333183 involved in salt tolerance were identified by GWAS and WGCNA. Furthermore, they were verified to affect salt tolerance by candidate gene association analysis. Salt stress influences maize growth and development. To decode the genetic basis and hub genes controlling salt tolerance is a meaningful exploration for cultivating salt-tolerant maize varieties. Herein, we used an association panel consisting of 305 lines to identify the genetic loci responsible for Na+- and K+-related traits in maize seedlings. Under the salt stress, seven significant single nucleotide polymorphisms were identified using a genome-wide association study, and 120 genes were obtained by scanning the linkage disequilibrium regions of these loci. According to the transcriptome data of the above 120 genes under salinity treatment, we conducted a weighted gene co-expression network analysis. Combined the gene annotations, two SNaC/SKC (shoot Na+ content/shoot K+ content)-associated genes GRMZM2G075104 and GRMZM2G333183 were finally identified as the hub genes involved in salt tolerance. Subsequently, these two genes were verified to affect salt tolerance of maize seedlings by candidate gene association analysis. Haplotypes TTGTCCG-CT and CTT were determined as favorable/salt-tolerance haplotypes for GRMZM2G075104 and GRMZM2G333183, respectively. These findings provide novel insights into genetic architectures underlying maize salt tolerance and contribute to the cultivation of salt-tolerant varieties in maize.

Genetic Dissection of Grain Yield of Maize and Yield-Related Traits Through Association Mapping and Genomic Prediction

High yield is the primary objective of maize breeding. Genomic dissection of grain yield and yield-related traits contribute to understanding the yield formation and improving the yield of maize. In this study, two genome-wide association study (GWAS) methods and genomic prediction were made on an association panel of 309 inbred lines. GWAS analyses revealed 22 significant trait-marker associations for grain yield per plant (GYP) and yield-related traits. Genomic prediction analyses showed that reproducing kernel Hilbert space (RKHS) outperformed the other four models based on GWAS-derived markers for GYP, ear weight, kernel number per ear and row, ear length, and ear diameter, whereas genomic best linear unbiased prediction (GBLUP) showed a slight superiority over other modes in most subsets of the trait-associated marker (TAM) for thousand kernel weight and kernel row number. The prediction accuracy could be improved when significant single-nucleotide polymorphisms were fitted as the fixed effects. Integrating information on population structure into the fixed model did not improve the prediction performance. For GYP, the prediction accuracy of TAMs derived from fixed and random model Circulating Probability Unification (FarmCPU) was comparable to that of the compressed mixed linear model (CMLM). For yield-related traits, CMLM-derived markers provided better accuracies than FarmCPU-derived markers in most scenarios. Compared with all markers, TAMs could effectively improve the prediction accuracies for GYP and yield-related traits. For eight traits, moderate- and high-prediction accuracies were achieved using TAMs. Taken together, genomic prediction incorporating prior information detected by GWAS could be a promising strategy to improve the grain yield of maize.