A correction has been applied to this article in:
Corrigendum: Unravelling the genetic framework associated with grain quality and yield-related traits in maize (Zea mays L.)
Mehak Sethi1‡
Dinesh Kumar Saini2†‡
Veena Devi1
Charanjeet Kaur3
Mohini Prabha Singh4
Jasneet Singh5Gomsie Pruthi6Amanpreet Kaur1Alla Singh1
Dharam Paul Chaudhary1*- 1Division of Biochemistry, Indian Institute of Maize Research, Ludhiana, Punjab, India
- 2Department of Plant Breeding and Genetics, Punjab Agricultural University, Ludhiana, Punjab, India
- 3Department of Basic Sciences and Humanities, Punjab Agricultural University, Ludhiana, Punjab, India
- 4Department of Floriculture and Landscaping, Punjab Agricultural University, Ludhiana, Punjab, India
- 5Agricultural and Environmental Sciences, Macdonald Campus, McGill University, Montreal, QC, Canada
- 6Department of Biotechnology, Punjab Agricultural University, Ludhiana, Punjab, India
Maize serves as a crucial nutrient reservoir for a significant portion of the global population. However, to effectively address the growing world population’s hidden hunger, it is essential to focus on two key aspects: biofortification of maize and improving its yield potential through advanced breeding techniques. Moreover, the coordination of multiple targets within a single breeding program poses a complex challenge. This study compiled mapping studies conducted over the past decade, identifying quantitative trait loci associated with grain quality and yield related traits in maize. Meta-QTL analysis of 2,974 QTLs for 169 component traits (associated with quality and yield related traits) revealed 68 MQTLs across different genetic backgrounds and environments. Most of these MQTLs were further validated using the data from genome-wide association studies (GWAS). Further, ten MQTLs, referred to as breeding-friendly MQTLs (BF-MQTLs), with a significant phenotypic variation explained over 10% and confidence interval less than 2 Mb, were shortlisted. BF-MQTLs were further used to identify potential candidate genes, including 59 genes encoding important proteins/products involved in essential metabolic pathways. Five BF-MQTLs associated with both quality and yield traits were also recommended to be utilized in future breeding programs. Synteny analysis with wheat and rice genomes revealed conserved regions across the genomes, indicating these hotspot regions as validated targets for developing biofortified, high-yielding maize varieties in future breeding programs. After validation, the identified candidate genes can also be utilized to effectively model the plant architecture and enhance desirable quality traits through various approaches such as marker-assisted breeding, genetic engineering, and genome editing.
1 Introduction
Maize (Zea mays L.), the third most prevalent cereal crop after wheat and rice, is consumed as a staple food in many developed and developing countries worldwide (Prasanna et al., 2001). Being a C4 and day-neutral plant, maize has less water requirement and wider adaptability under a range of agro-climatic conditions (Kaur et al., 2022). Combating the ever-increasing demands for food due to the continuously increasing population is an extreme challenge for developing countries. With the emergence of the green revolution, the yield-related parameters took the whole attention of maize breeders, thus the nutritional quality remained a complex problem (Maqbool et al., 2021). Later, improvement in breeding strategies to develop nutritionally superior maize resulted in poor agronomic characteristics leading to huge productivity loss. Since the nutritional value is determined by both physiological traits and nutrient content, knowledge of the inheritance of grain yield, nutrient content (macro and micro-nutrients), as well as the genetic relationships among these traits, can facilitate maize breeders working on biofortified varieties with enhanced yield potential (Das et al., 2021).
Traits contributing to yield and grain quality are usually antagonist to each other. Moreover, such traits are quantitatively inherited and controlled by number of genetic loci. Numerous mapping studies have been conducted in maize which have led to the identification of hundreds of QTLs associated with nutritional quality traits such as protein, starch, fat, carotenoids, tocopherol, phosphorus contents, etc., and yield-related traits including grain yield per se, kernel length, kernel width, kernel weight, etc. (Chen et al., 2008; Li et al., 2012; Zdunic et al., 2014; Venado et al., 2017; Fei et al., 2022; Xu et al., 2022; Banerjee et al., 2023).
The validity of QTL mapping findings is influenced by various factors, such as the specific experimental conditions, size and composition of the mapping population, density of genetic markers employed, statistical methodologies utilized, and the presence of genotype-environment and epistatic interactions (Goffinet and Gerber, 2000; Mackay, 2001; Arcade et al., 2004; Halladakeri et al., 2023). Additionally, it is common for QTLs to span relatively large genetic intervals, which can pose challenges when transferring desired QTLs through marker-assisted breeding (MAB). This phenomenon is often referred to as “linkage drag,” where undesirable genetic regions linked to the target QTL are unintentionally transferred along with it. To overcome this limitation, it is essential to precisely localize QTLs within narrow genetic intervals. This not only minimizes the potential for linkage drag but also enhances the efficiency and effectiveness of MAB and QTL cloning efforts. To enhance the reliability and robustness of QTL studies, researchers often employ multiple mapping populations evaluated across diverse locations and over multiple years. This approach helps account for genotype-environment interactions, ensuring that the identified QTLs exhibit consistent effects across different genetic backgrounds and environmental conditions. Consensus QTLs that consistently demonstrate significant associations with the target trait across multiple populations and environments are regarded as particularly suitable candidates for integration into MAB programs (Goffinet and Gerber, 2000; Halladakeri et al., 2023).
Meta-QTL analysis represents a conceptual approach that combines QTL discoveries from various studies and refines the positions of QTLs on a consensus map (Goffinet and Gerber, 2000). When numerous QTLs associated with the specific trait are analysed together through meta-analysis, the resulting consensus QTLs are referred to as “meta-QTLs (MQTLs).” This approach can not only identify redundant QTLs and candidate genes (CGs), but also facilitate the selection of promising QTLs for breeding applications involving MAB (Goffinet and Gerber, 2000; Sandhu et al., 2021; Saini et al., 2022; Sheoran et al., 2022; Tanin et al., 2022; Kumar A. et al., 2023; Halladakeri et al., 2023; Karnatam et al., 2023).
In maize, significant advancements have been made in the field of MQTL analysis, which involves identifying genetic regions associated with various traits such as fungal disease resistance (Gupta et al., 2023), multiple abiotic stress tolerance (Sheoran et al., 2022), popping traits (Kaur et al., 2021), and root-related traits (Karnatam et al., 2023). Additionally, specific MQTL studies have focused on yield-related traits (Semagn et al., 2013; Wang Y. et al., 2016; Wang et al., 2020; Chen et al., 2017; Zhou et al., 2020) or quality-related traits (Jin et al., 2013; Dong et al., 2015). However, it is important to note that none of the previous studies have conducted MQTL analysis for QTLs associated with both yield-related and quality traits in maize. Additionally, several new studies on QTL analysis for both yield-related and quality traits in maize have been conducted since the publication of earlier studies, but they have not been utilized for predicting MQTLs or identifying promising candidate genes (CGs). In our present study, we aimed to address this gap by simultaneously investigating the genetic regions influencing both yield-related and quality traits in maize. As a result, we can contribute to the development of improved maize varieties with enhanced yield and quality traits, benefiting both farmers and consumers alike.
The current study utilized QTL data from 56 studies associated with various yield-related and quality traits to perform MQTL analysis. The primary objectives of this study were as follows: i) constructing a high-density consensus genetic map, ii) identifying robust and consistent MQTLs with narrower confidence intervals compared to initial QTLs, iii) determining precise flanking markers for MQTLs to facilitate MAB, iv) validating MQTLs through genome-wide association studies (GWAS), v) identifying CGs within promising MQTL regions, and vi) conducting expression analysis of the identified CGs. It is anticipated that these efforts will greatly enhance the selection efficiency for various yield-related and quality traits in maize.
2 Materials and methods
In the present analysis, a MQTL analysis was conducted to identify the genomic regions associated with yield and quality-related traits in maize. The analysis involved following five key steps: i) bibliographic search and compilation of QTL mapping studies related to quality and yield-associated traits, ii) integration of high quality linkage maps and markers from the individual studies to create a consensus map, iii) prediction of MQTLs and their validation using GWAS, iv) mining genes within potential MQTL regions and in silico expression analysis, v) synteny analysis among the maize, wheat, and rice for the key genomic regions. By employing these steps, the study aimed to gain insights into the genetic basis underlying both yield and quality-related traits in maize and related cereals.
2.1 Bibliographic search and compilation of QTL mapping studies
The various traits that impact both maize quality and yield are illustrated in Supplementary Table S1. Several QTL mapping studies have highlighted the significance of these traits in defining maize quality and enhancing yield. To gather relevant information, these studies were collected from online platforms such as Google Scholar (https://scholar.google.com/) and PubMed Central (https://www.ncbi.nlm.nih.gov/pmc/), using appropriate keywords. The collected studies provided essential data for the meta-QTL analysis, including the size and type of the mapping population, the specific map used in each study, the marker positions associated with QTLs, the traits linked to the QTLs, the confidence interval of the markers, the LOD score, and the phenotypic variation explained by individual QTLs (represented as R2 or PVE). However, certain studies were not included in the meta-analysis due to the lack of important data. Further, to facilitate clarity and uniqueness, the QTLs used in the present study were assigned distinctive names, which included the trait name, followed by the chromosome number, and finally, the number of the QTL on the chromosome.
2.2 Consensus map development and QTL projection
To construct a highly dense consensus map, the present study incorporated two high-quality, high resolution, mixed marker (involving RFLP, SSR, SNP and InDels) linkage maps, i.e., “ISU Integrated IBM 2009” and IBM2 (intermated B73/Mo17) (https://www.maizegdb.org/data_center/map) These maps were carefully compiled and enriched to ensure accuracy and reliability. Furthermore, all the markers flanking the QTLs extracted from different studies were integrated into this consensus map, facilitating the projection of a greater number of QTLs. The construction of the consensus map was accomplished using the LPMerge package within the R environment as described by Kumar R et al (2022).
The collected QTLs underwent a rigorous screening process to identify reliable QTLs with complete information essential for projection and meta-analysis. Our initial selection criteria involved identifying stable and consistent QTLs across diverse environments within each study, while excluding those that were specific to particular environments (i.e., environment-specific QTLs). Additionally, QTLs without available LOD or PVE values were excluded, as well as those lacking information on flanking markers and genetic positions. Subsequently, the selected QTLs were projected onto the consensus map using BioMercator v4.2.3 software, which served as the projection tool (ChardonVirlon et al., 2004). For QTL projection optimization, QTLProj, a dynamic strategy, was employed. This strategy utilizes pairs of common markers flanking the QTLs in the original map, along with an estimate of the distance between the initial and consensus maps. The QTL projection procedure and map distance calculation are influenced by the minimum value of the flanking marker distance ratio and the minimum p-value. These parameters were carefully adjusted to ensure the similarity of flanking marker interval distances between the initial and consensus maps, thereby enhancing the accuracy of QTL projection.
2.3 Prediction of MQTLs
In this study, maize quality and yield-related traits were considered as meta-traits. A meta-QTL refers to a genomic region that contains two or more QTLs for a trait, obtained from at least two different mapping studies. The prediction of meta-QTLs was carried out using a two-step process described by Veyrieras et al. (2007). In the first step, the best MQTL model was selected by comparing five different models: i) AIC (Akaike information criterion), ii) AICc (AIC correction), iii) AIC3 (AIC 3 candidate models), iv) BIC (Bayesian information criterion), and v) AWE (average weight of evidence). The model with the lowest value among three out of the five models was chosen. The second step involved determining the total number of MQTLs on a chromosome, their weightage according to the selected model, and the confidence interval (CI).
To calculate the LOD and PVE values of the predicted MQTLs, the average of the LOD and PVE values of the initial QTLs comprising the meta-QTL was taken. In cases where multiple MQTLs were predicted on a single chromosome; a standardized naming procedure was followed. Each meta-QTL was assigned a name starting with “MQTL,” followed by the numeric representation of the chromosome number, a period, and another numeric value indicating the sequence of the MQTL on that specific chromosome (e.g., MQTL1.1, MQTL1.2). To obtain the physical position of each meta-QTL, the flanking markers associated with them were used. Published studies or online databases such as MaizeGDB (https://www.maizegdb.org/), GrainGenes (https://wheat.pw.usda.gov/GG3/), and Gramene (https://www.gramene.org/) were utilized to obtain information on the nucleotide sequences of the flanking markers. The obtained sequence information was then subjected to BLAST analysis against the reference genome of maize available on the Ensembl Plants database (https://plants.ensembl.org/index.html).
2.4 Validation of MQTLs using GWAS
To validate the projected MQTLs, a validation process was conducted using 10 independent GWAS studies focused on traits related to grain yield and quality (viz., Li et al., 2013; Liu et al., 2016; Fang et al., 2021; Ma and Cao, 2021; Ndlovu et al., 2022; Qu et al., 2022; Suwarno et al., 2015; Zheng et al., 2021; Zeng et al., 2022; Zhang et al., 2022). These GWAS studies encompassed maize populations consisting of 180–410 genotypes/accessions. To determine the physical position of the identified significant markers (marker-trait associations; MTAs) from the GWAS studies, the database MaizeGDB was utilized. The retrieved physical positions were then compared to the physical positions of the MQTLs present on individual chromosomes. MQTLs that co-localized with markers from the significant GWAS studies were considered as GWAS-verified MQTLs, further enhancing their credibility and relevance in the context of the study.
2.5 Selection of breeder friendly MQTLs
The MQTLs investigated in this study cover broad genomic regions impacting yield and quality related traits in maize. To enable practical implementation by plant breeders, “breeder friendly MQTLs” were selected based on the following criteria: PVE >10%, genetic confidence interval <2 cM, physical confidence interval <1 Mb, and involvement of at least 10 initial QTLs. These precise MQTLs provide targeted regions of interest for plant breeders to enhance grain yield and quality in maize varieties.
2.6 Candidate gene mining within breeder friendly MQTL regions
The regions identified as BF-MQTLs were further investigated for candidate gene mining and expression analysis related to maize quality and yield traits. MaizeGDB, utilizing the MaizeMine tool, was employed for candidate gene mining. The genomic regions corresponding to each MQTL were analysed to identify genes of interest within the Zm-B73-REFERENCE-NAM-5.0 assembly. This approach provided valuable information about the various genes encompassed within the selected genomic regions.
2.7 Expression analysis of candidate genes and their syntenic regions in wheat and rice genomes
Each gene identified within the MQTL region was analysed for expression patterns using the qTeller tool available on MaizeGDB. qTeller provides information on gene expression across different tissues throughout maize development. To gain insights into the functional characteristics of the identified genes, as well as to explore potential orthologous relationships, functional descriptions and orthologous studies in wheat (Triticum aestivum) and rice (Oryza sativa) were conducted using the MaizeMine tool on MaizeGDB. The wheat and rice homologous collected through MaizeMine were assessed through Ensembl Plants to retrieve their chromosomal location on respective genomes. MQTLs containing similar gene models located on conserved genomic regions across wheat, rice, and maize were considered as ortho-MQTLs.
3 Results
3.1 Bibliographic search and collected QTL information
A comprehensive analysis was conducted, encompassing data from 56 interval mapping studies, yielding a total of 2,974 QTLs associated with 169 component traits, with 48 yield-associated traits, 23 oil-associated traits, 28 protein-associated traits, 13 starch-associated traits, 25 carotenoid-associated traits, 16 metal ion associated traits, 10 tocopherol associated traits and 6 traits deciphering associated macromolecules (Supplementary Table S1). Among these studies, 37 studies focused on quality-associated traits, 10 on yield-related traits and 9 studies had a combined focus on both yield and maize grain quality (Supplementary Table S2). The studies employed 67 different bi-parental populations, including doubled haploid (DH), recombinant inbred lines (RILs), F2, and backcross populations. It is worth noting that certain populations were assessed for various traits in multiple studies. The range of population sizes utilized in the study varied from 10 to 4,699 individuals (Supplementary Table S2).
The number of QTLs varied across different categories of traits, with 1,877 QTLs identified for quality-related traits, 387 QTLs for yield-related traits, and 710 QTLs for combined quality and yield-related traits (Figure 1). The distribution of collected QTLs across the 10 maize chromosomes was uneven, ranging from 190 QTLs on chromosome 10 to 377 QTLs on chromosome 1, with an average of approximately 284 QTLs per chromosome. The individual QTLs exhibited a wide range of LOD scores, spanning from 1.88 to 75.3, and averaging at 4.71. Notably, over 80% of the QTLs had LOD scores falling within the 3 to 75.3 range. As for the PVE, it varied from 0.01% to 49.3%, with an average of 9.01%. These QTLs demonstrated a typical L-shaped distribution, with approximately 70% of them showing a PVE of 10% or less. Regarding the confidence intervals (CIs) of the individual QTLs, they ranged from 0 to 119.5 cM, with an average of 11.31 cM. Interestingly, a considerable portion of the QTLs (more than 38%) had CIs equal to or greater than 10 cM.
FIGURE 1
FIGURE 1. Distribution of QTLs associated with quality and yield related traits.
3.2 Features of consensus map
The consensus map integrated an extensive set of 57,523 markers, encompassing a diverse range of SSR, RFLP, and SNP markers. These markers were strategically distributed across a genetic span of 3,959.9 cM. Across the individual chromosomes, the lengths exhibited notable variation, ranging from 157.4 cM for chromosome 10 to 1,135.4 cM for chromosome 1, with an overall average of 395.99 cM. Remarkably, an average of approximately 14 markers were identified per cM, emphasizing the marker-rich nature of the map. Examining the marker density across chromosomes, the maximum density of 34.3 markers per cM was observed on chromosome 5, signifying a concentrated marker distribution. Conversely, chromosome 7 displayed the minimum marker density of 8.9 markers per cM, indicating a relatively sparser marker distribution on that particular chromosome.
3.3 QTLs projected on consensus map
Among the QTLs selected from the collected QTLs, as many as 965 were successfully projected onto the consensus map. The 965 projected QTLs were grouped into 75 hotspot regions. However, seven hotspot regions consisting of only a single QTL from an individual study were not considered true MQTLs. As a result, a total of 68 MQTLs were obtained from the present study, distributed across the 10 chromosomes. The highest number of MQTLs, nine in total, were identified on chromosomes 2 and 7, while the minimum of two MQTLs were found on chromosome 9 (Figure 2A). The number of QTLs within each MQTL varied, ranging from 2 (MQTL 3.9) to 64 (MQTL 1.1). Sixteen MQTLs comprised 20 or more QTLs from different studies (Figure 2B). The proportion of MQTLs and the number of QTLs varied on each chromosome (Supplementary Table S3).
FIGURE 2
FIGURE 2. Some characteristic features of QTLs and MQTLs (A) Distribution of MQTLs on different maize chromosomes; (B) Number of QTLs under each MQTL; (C) Average CI of QTLs and MQTLs available on different maize chromosomes.
All the MQTLs identified in this study exhibited LOD scores ranging from 3 to 12.25 and phenotypic variation explained (PVE) ranging from 3.86% to 21.41%. Among the reported MQTLs, 19 had high LOD scores (≥5) and PVE values (≥10%), namely,: MQTL1.4, MQTL2.4, MQTL2.8, MQTL3.5, MQTL5.3, MQTL6.4, MQTL6.8, MQTL7.2, MQTL7.3, MQTL7.4, MQTL7.5, MQTL7.7, MQTL7.8, MQTL8.1, MQTL8.3, MQTL8.4, MQTL8.6, MQTL10.7, MQTL10.8 (Figure 2B). The average CI of the MQTLs was 1.59 cM, which was 6.12 times smaller than the average CI of the initial QTLs (11.32 cM). Significant reductions in CI were observed on different chromosomes, with the highest fold change observed on chromosome 1 (14.21-times) and the lowest on chromosome 9 (2.4-times) (refer to Figure 2C).
All the MQTLs were physically mapped onto the maize reference genome Zm-B73-REFERENCE-NAM-5.0. The physical CI of these MQTLs ranged from 2,049 bp (MQTL7.3) to 39.95 Mb (MQTL 4.4), with a mean physical CI of 3.30 Mb. Collectively, these MQTLs occupied a physical length of 224.52 Mb across the genome (refer to Figure 3). Among the 68 MQTLs, two were exclusively associated with yield related traits (MQTL10.7, MQTL10.8), and 14 were exclusively associated with quality traits (MQTL1.5, 1.6, 2.9, 3.3, 3.5, 4.1, 4.2, 6.1, 7.6, 7.7, 7.8, 8.1, 8.2, 8.3, 8.4). Notably, MQTL6.1 was found to be associated with carotenoid synthesis, while MQTL4.2 was exclusively associated with amino acid biosynthesis (Figure 3).
FIGURE 3
FIGURE 3. Distribution of MQTLs on different maize chromosomes. Common: MQTLs associated with both quality and yield-associated traits, BFQ: breeder-friendly quality trait MQTLs, BFY: breeder-friendly yield trait MQTLs, BFC: breeder-friendly common MQTLs (involving both quality and yield-related traits).
3.4 Breeders friendly MQTLs
Out of the 68 MQTLs reported, a subset of 10 MQTLs was selected and designated as “breeder-friendly MQTLs” (BF-MQTLs). The selection criteria for identifying BF-MQTLs were as follows: a PVE value of more than 10%, a genetic CI of less than 2 cM, a physical CI of less than 1 Mb, and involvement of at least 10 initial QTLs from different studies. Detailed information on these BF-MQTLs can be found in Table 1. Among the 10 BF-MQTLs, eight MQTLs encompassed QTLs associated with both quality and yield-related traits, while one MQTL included QTLs solely for quality-related traits, and another MQTL comprised QTLs exclusively associated with yield-related traits (refer to Figure 3).
TABLE 1
TABLE 1. Detailed information of Breeder Friendly MQTLs. These MQTLs were selected based on the following criteria: PVE>10%, CI < 2 cM, and number of initial QTLs involved> 10.
3.5 Validation of MQTLs with the GWAS based MTAs and the availability of known genes within MQTL regions
The physical positions of MQTLs identified through meta-analysis were compared with the physical positions of significant markers (MTAs) identified for the same traits using the GWAS approach. Data obtained from 10 different GWAS studies provided information on 250 MTAs, that overlapped with 47 MQTLs out of 68 reported MQTLs. The count of overlapping MTAs or adjacent MTAs (within a 1 Mb region surrounding the MQTL regions) varied for each MQTL, with MQTL9.3 exhibiting the highest number of overlaps with 36 MTAs. Notably, among the overlapped MTAs, 35 were associated with oil-related traits, 26 were associated with yield-related traits, 17 were associated with carotenoid-related traits, 2 were associated with sugar-related traits, and 1 was associated with protein-related traits (see Supplementary Table S4 and Supplementary Figure S1). It is important to mention that these GWAS-validated MQTLs also included all the breeder-friendly MQTLs (BF-MQTLs).
GWAS studies provided valuable insights into the genes associated with the traits under investigation (see Supplementary Table S4). Specifically, a total of 13 genes were found to be co-localized with MQTLs affecting yield. These included one gene associated with cob weight (i.e., Zm00001d002000), two genes associated with kernel length (viz., Zm00001d036024, Zm00001d000184), three genes associated with kernel number per row (viz., GRMZM2G049091, GRMZM2G138067, GRMZM2G446921), one gene associated with kernel thickness (i.e., Zm00001d041972), two genes associated with kernel weight (viz., GRMZM2G010933, GRMZM2G464985), one gene associated with moisture content (i.e., GRMZM2G069024), and four genes associated with grain yield (viz., GRMZM2G391473, GRMZM2G125557, GRMZM5G845736, GRMZM2G151649). Interestingly, three breeder-friendly MQTLs (BF-MQTL2.8, BF-MQTL7.3, BF-MQTL7.4) were found to be associated with these reported genes that influence grain yield.
Furthermore, 13 genes were identified to be co-localized with some MQTLs affecting grain quality. Among them, 10 genes were associated with grain oil content (viz., GRMZM2G059138, GRMZM5G828253, GRMZM2G118423, GRMZM2G176542, GRMZM2G125268, GRMZM2G122767, GRMZM2G417435, GRMZM2G169089, GRMZM2G077789, GRMZM2G136072), two genes were associated with starch content (viz., GRMZM2G404453, GRMZM2G104325), and one gene was validated for grain protein content (i.e., GRMZM2G049681). Two breeder-friendly MQTLs (viz., BF-MQTL6.4, BF-MQTL6.8) were found to be associated with these reported genes affecting grain quality. Notably, these genes encode essential enzymes (e.g., cytochrome c oxidase copper chaperone 1, serine/threonine-protein kinase D6PKL1, beta-glucosidase 11) involved in significant metabolic pathways and transcription factors (e.g., transcription factor MYBS3, transcription initiation factor IIF beta subunit), as well as proteins (e.g., protein LURP-one-related 5, F-box domain containing protein) that are believed to play crucial roles in influencing grain yield and quality (see Supplementary Table S4).
3.6 In silico gene mining and tissue-specific expression analysis
The genes residing within the genomic regions of BF-MQTLs were investigated using maizeGDB database. A total of 59 genes were identified within these BF-MQTL regions, which are located on chromosomes 1, 2, 6, 8, and 10 (see Figure 4). Among the analysed genes, 15 genes were found to encode enzymes involved in crucial metabolic processes such as lipid, amino acid, protein, sugar, nucleotide and secondary metabolite metabolism. Some enzymes were also involved in isomerization, endomembrane homeostasis and gene regulation, 10 genes encoded versatile groups of transcription factors (e.g., C2H2-type domain-containing protein, eukaryotic translation initiation factor 3 subunit E, Mediator of RNA polymerase II transcription subunit 17, NAC domain-containing protein, RING-H2 finger protein ATL80, transcription factor ICE1, transcription factor PHYTOCHROME INTERACTING FACTOR-LIKE 13, transposition associated factors, zinc finger protein WIP2), 6 genes encoded transporter proteins (e.g., magnesium transporter (2), protein transport protein Sec61 subunit gamma, HMA domain-containing protein, adenine/guanine permease AZG1, putative VHS/GAT domain containing family protein, 6 genes were associated with signalling (e.g., ras-related protein RABC1, calcium mediated signalling (3), auxin-responsive protein SAUR22, tr-type G domain-containing protein), 2 genes each were related to apoptosis (e.g., senescence associated gene 20, ubiquitin carboxyl-terminal hydrolase 19), ribonuclease (e.g., U5 small nuclear ribonucleoprotein 40 kDa protein, ribonuclease P protein subunit p29), and chloroplast (e.g., multiple chloroplast division site 1, protein EXECUTER 2 chloroplastic) proteins, and one gene was reported to function in peroxisomes (i.e., peroxisome biogenesis protein 22) (see Figure 5).
FIGURE 4
FIGURE 4. Diagram deciphering the salient characteristics of MQTLs, BF-MQTLs, candidate genes. Outermost circle represents the location of MQTLs predicted during the present study. Second circle (blue) provides an overview of density of MQTLs on different chromosomes. Third circle represents locations and densities of BF-MQTL along with putative gene density. Lines within innermost circle represents phylogenetic relationships between identified genes.
FIGURE 5
FIGURE 5. Biological functions of various classes of genes annotated under BF-MQTLs region.
Additionally, the expression analysis of these genes was performed using the expression analysis tool within the MaizeMine, revealing tissue-specific expression patterns (Figure 6; Supplementary Table S5). Among the eight different tissues studied, a single gene showed expression in anthers (i.e., Zm00001eb034640), four genes were expressed in shoot axis internodes, 41 genes were expressed in stem and shoot apical meristems, 49 genes were expressed in reproductive tissues, 52 genes were expressed in roots and internodes, 53 genes were expressed in leaves, and 54 genes were expressed in seeds.
FIGURE 6
FIGURE 6. Tissue-specific expression of different genes available from MQTL regions (A: anther, IN: internode, L: leaf, RP: reproductive, R: root, S: shoot, SA: shoot axis internode, SSAM: stem and shoot apical meristem).
3.7 Synteny of genes identified in maize MQTL regions with wheat and rice genomes
Synteny analysis was conducted to identify conserved regions of the reported maize MQTLs in the genomes of wheat and rice. In total, 59 genes identified in the maize genome within the 10 BF-MQTL regions corresponded to 42 genes in the rice genome (Figure 7A) and 150 genes in the wheat genome (see Figure 7B). The genes obtained from the BF-MQTLs showed variation in the number of orthologs in the wheat genome, with the maximum number of orthologs (34) observed for genes derived from MQTL10.8, while the minimum (12) was observed for genes from MQTL2.4. Similarly, in the rice genome, the maximum number of orthologous genes was reported for genes derived from MQTL1.4, while the minimum (3) was observed for genes from MQTL2.8. Overall, the synteny analysis demonstrated the presence of conserved genomic regions between maize and rice, as well as between maize and wheat. These regions can be considered as ortho-MQTL regions (see Figure 7).
FIGURE 7
FIGURE 7. Syntenic relationships of genes available from the maize MQTL regions with the (A) rice and (B) wheat genomes. The physical lengths of the chromosomes are indicated by the rulers drawn above on each chromosome.
4 Discussion
Grain quality and grain yield are critical traits in maize breeding, as they determine the nutritional value and productivity of the crop. Maize is a major source of macronutrients such as protein, starch, and fat, as well as micronutrients and phytochemicals including tocopherols, carotenoids, metal ions (Fe, Zn, etc.), and phytic acid, which are essential for the maize-consuming population (Owens et al., 2014; Tanumihardjo et al., 2019). Numerous studies have been conducted over the years, identifying a large number of QTLs associated with various traits related to maize grain quality and yield (Supplementary Table S1). However, most of the QTLs reported in previous studies have low PVE values, wide CIs and other potential limitations, limiting their potential for utilization in breeding programs aimed at developing nutrient-rich, high-yielding maize varieties. Moreover, several studies have demonstrated the co-localization of QTLs linked to various aspects of grain quality and yield-related characteristics in maize (Yang et al., 2013). These findings suggest that a single genetic locus may exert influence over multiple facets of grain development. Consequently, investigating these genomic regions and analysing the hotspot areas that affect both grain quality and yield can significantly enhance future maize breeding programs.
Meta-analyses of QTLs associated with a variety of traits have been recently conducted in different crops such as wheat (Kumar et al., 2021; Kumar et al., 2022 A. C.; Kumar et al., 2023 S.; Saini et al., 2021; Saini et al., 2022; Tanin et al., 2022), rice (Sandhu et al., 2021; Kumari et al., 2023), barley (Akbari et al., 2022), common bean (Shafi et al., 2022), pigeon pea (Halladakeri et al., 2023), including maize (Kaur et al., 2021; Sheoran et al., 2022; Wang et al., 2022; Gupta et al., 2023; Karnatam et al., 2023), for diverse traits, including both yield-related traits (Semagn et al., 2013; Wang Y. et al., 2016; 2020; Chen et al., 2017; Zhou et al., 2020) and quality traits (Jin et al., 2013; Dong et al., 2015). However, there is currently no comprehensive study on the genomic regions influencing both grain quality and yield in maize. This study aims to fill this gap by identifying QTL clusters in defined regions using a large dataset of yield-related and quality trait data from multiple studies.
The findings of this study can significantly contribute to advancing maize breeding programs. Further, our study represents a significant advancement over previous meta-analyses conducted on various yield and quality-related traits. It offers a more up-to-date and comprehensive analysis, distinguishing itself in several key aspects (Supplementary Table S6). The present study offers a broad perspective on recent mapping studies (published up to 2022) related to maize grain yield and quality. The consensus map used in this study is comprehensive and densely populated, covering the entire maize genome with high precision. A total of 57,523 markers, including SSR, SNP, and RFLP markers, were employed in the consensus maps, distributed across a length of 3,959.9 cM. The highest number of markers was reported on chromosome 1 (10,341), while the lowest was observed on chromosome 10 (4,566). The distribution of markers and QTL density on maize chromosomes observed in this study aligns with findings from previous meta-QTL studies on maize (Sheoran et al., 2022; Wang et al., 2022).
A total of 75 hotspot regions were distributed across the 10 chromosomes in the present study. The majority of these hotspot regions were found in the upper terminal regions of the chromosomes, which have been previously identified as regions with high gene density (Ghaffari et al., 2013). These regions exhibited multiple yield and quality-associated QTLs within a single MQTL region, indicating a strong correlation of these genomic regions in enhancing maize grain quality and yield. The hotspot regions established in present study can be considered in future breeding programs to obtain high yield and high-quality hybrids. Several biochemical and mapping studies have been conducted previously to establish correlation between yield and quality associated traits (Zheng et al., 2021; Amegbor et al., 2022).
In the context of maize, grain nutrient quality is primarily determined by the levels of tryptophan and lysine in the maize protein. Standard maize varieties generally lack sufficient amounts of these essential amino acids, resulting in reduced overall protein quality. This deficiency has led to a significant focus on quality protein maize breeding over the past decade (Prasanna et al., 2001; Prasanna et al., 2020). Maize grain quality-associated traits are complex, and they are regulated by multiple major genes and lots of minor genes. It is reported that there is a negative association between grain yield and protein quality. Tryptophan content showed a significant positive genotypic association with moisture, oil, and fiber content, but a strong negative correlation with protein and starch content (Amegbor et al., 2022). Given the intricate interplay between various quality and yield-related traits, gene pyramiding is strongly recommended to develop biofortified, high-yielding hybrids (Zheng et al., 2021). Despite several mapping studies, the goal of achieving both high yield and grain quality has proven challenging, with fluctuations observed under varying environmental conditions. The common hotspot regions for both grain yield and grain quality traits (BF-MQTL1.4, 2.4, 2.8, 6.4, 6.8, 7.2, 7.3, 7.4) can be beneficial for future breeding programs for targeting associations between yield and quality traits. Although, the significant inverse correlations observed between various traits along with the presence of shared or linked QTLs and epistasis with opposing effects, may pose potential challenges for plant breeders aiming to simultaneously enhance these traits in maize.
The present study also refined the CIs of different QTLs for various quality and yield traits, resulting in an average reduction of 6.12 times compared to the initial QTLs. Further, GWAS was employed as a precision technique to validate the reported MQTLs, with approximately 62% of the MQTLs being validated through previous GWA studies (Suwarno et al., 2015; Ma and Cao, 2021; Zheng et al., 2021; Ndlovu et al., 2022; Zeng et al., 2022; Zhang et al., 2022). MQTL9.3 showed the highest overlap with a maximum number of MTAs in GWA studies, with 36 MTAs identified. Additionally, 27 MQTLs overlapped with only one MTA. Similar results of MQTL validation through MTAs have been reported in previous studies, with validation rates of 54.6% and 63% (Saini et al., 2021; Saini et al., 2022).
The MQTLs validated through GWAS were further characterized as breeder-friendly based on their high PVE and low CI values. Among the 10 validated BF- MQTLs in the present study, a total of 59 CGs were reported. Supplementary Table S7 presents the orthologs of these 59 CGs in model crops, including Arabidopsis and rice, along with their respective function descriptions. Notably, some of these orthologs have been extensively studied and their associations with the traits under study have been well-documented. Among these BF-MQTLs, five (BF-MQTL 1.4, 2.4, 2.8, 6.4, 6.8) exhibited QTLs associated with both grain quality and yield. Among these promising BF-MQTLs, MQTL 6.8 had the highest number of associated genes (12), followed by MQTL1.4 with 10 genes, and 7 genes each for MQTL2.4, MQTL2.8, and MQTL6.4. Based on these findings, BF-MQTL1.4, 2.4, 2.8, 6.4, and 6.8 appear to be the most promising candidates for improving both grain yield and quality. However, further investigations are required to identify functional variants of these genes within the MQTL regions.
Several genes identified within the validated hotspot regions encode enzymes involved in various metabolic processes, such as terpenoid biosynthesis (1-deoxy-D-xylulose-5-phosphate reductoisomerase) (Zhang et al., 2020), auxin dependent endomembrane homeostasis (1-phosphatidylinositol-3-phosphate 5-kinase (Lo et al., 2022), oil biosynthesis (acyltransferase) (Yan et al., 2018), photorespiration (alanine--glyoxylate transaminase which is known to facilitate the transamination between L-alanine and glyoxylate to produce pyruvate and glycine) (Liepman et al., 2019), lipid metabolism (aspartic proteinase oryzasin-1) (Alexandrov et al., 2009), carbohydrate metabolism (endoglucanase) (Jiao et al., 2017), sugar isomerization (epimerase domain-containing protein), phosphorus remobilization (glycerophosphodiester phosphodiesterase) (Wang J. et al., 2021), etc.
In addition to enzymes, some genes also encode for potent transcription factors. An example of such transcription factors is the C2H2-type domain-containing proteins, known as zinc finger proteins, which play a pivotal role in regulating plant growth and development (Li et al., 2022). Eukaryotic translation initiation factor 3 subunit E is involved in translating specific mRNA subsets related to cell proliferation. Mediator of RNA polymerase II transcription subunit 17 acts as a coactivator, regulating the transcription of almost all RNA polymerase II-dependent genes. NAC domain-containing proteins have diverse functions, including nuclear localization, DNA binding, and the formation of homodimers or heterodimers with other NAC domain-containing proteins (Nuruzzaman et al., 2013). Moreover, certain NAC-TFs, such as ZmNAC128 and ZmNAC130, have been identified to play crucial roles in the accumulation of starch and protein in maize, as highlighted in a study by Zhang et al. (2019). Pentatricopeptide repeat-containing proteins participate in RNA regulation and metabolism within plant organelles (Wang X. et al., 2021), and they are recognized for their multifaceted roles in plant growth and development (Li et al., 2022). The RING-H2 finger protein ATL80, featuring zinc finger domains, plays a role in stress response (Gao et al., 2012) and has been associated with influencing plant architecture and grain yield in rice (Yan et al., 2022). Another transcription factor, ICE1 MYC-like bHLH, acts as a transcriptional activator during cold stress response (Waititu et al., 2021) and is also known to enhance yield in rice. TF PHYTOCHROME INTERACTING FACTOR-LIKE 13 participates in light-induced transcription and plays a role in plant growth and internode elongation (Wu et al., 2019). Additionally, there are transposition-associated genes and the Zinc finger protein WIP2, which have been linked to plant growth (Wang et al., 2010; Wang B. et al., 2016; Li et al., 2022).
A few genes reported within the MQTL regions were categorized as peroxisome biogenesis protein 22 that participates in peroxisome assembly (Jiao et al., 2021), multiple chloroplast division sites 1 that leads to natural variation in chloroplast size (Kadirjan-Kalbach et al., 2019), protein EXECUTER 2 chloroplastic (Zhang et al., 2019), U5 small nuclear ribonucleoprotein (part of spliceosome), ribonuclease P protein subunit p29 which is involved in Mg2 -dependent hydrolysis (Ohtani, 2018). Several genes are implicated in transport across bio-membranes, for instance, the putative VHS/GAT domain-containing family protein is known to be involved in protein transport and energy metabolism, and its overexpression has been reported to enhance biomass (Nunes et al., 2019). Additionally, the magnesium transporter is recognized for its role in chlorophyll synthesis (Li et al., 2018; Li, 2020), while the protein transport protein Sec61 subunit gamma is responsible for mediating endoplasmic reticulum translocation (Zimmermann et al., 2011). Furthermore, the HMA domain-containing protein contains a conserved domain found in several heavy metal transport or detoxification proteins (Li et al., 2018; Cao et al., 2019), and the Adenine/guanine permease AZG1 plays a crucial role in facilitating adenine and guanine transport across membranes.
Two genes were associated with apoptosis involving senescence-associated gene 20, and ubiquitin carboxyl-terminal hydrolase 19 (Zhang et al., 2014). Several signalling associated proteins were encoded by above mentioned genes such as ras-related protein RABC1 it is a GTP-binding protein that regulates stomatal movements and drought stress responses by mediating the interaction with ABA (Wu et al., 2015), calcium mediated signalling acquisition of stress response (Jiao et al., 2017; Li, 2020), auxin-responsive protein SAUR22 that is known to promote cell expansion (Ren and Gray, 2015), calcium-binding EF hand family protein which actively bind to Ca2+ and chelate the cytosolic Ca to regulate Ca homeostasis (Mohanta et al., 2015). Tr-type G domain-containing protein which control a multitude of biological processes, ranging from cell division, cell cycling, and signal transduction, to ribosome assembly and protein synthesis (Jin et al., 2021). The above-mentioned genes exhibited tissue-specific expression, with maximum expression observed during the seed-filling stage, while fewer genes were expressed during the anther development.
During the synteny analysis, for MQTL1.4, the maximum number of orthologs were present on chromosome 9 in rice and on chromosomes 5 and 7 in wheat. These conserved genes encode major transcription factors and enzymes that play important roles in major metabolic pathways. For MQTL2.4 and 2.8, rice orthologs were found on chromosome 4 and on chromosome 2 in wheat, respectively. MQTL6.4 showed conserved sequences on chromosome 6 in rice and chromosome 7 in wheat. In MQTL6.8, orthologs were present on chromosome 5 in rice and on chromosomes 1 and 7 in wheat. In MQTL8.1, orthologs accumulated on chromosomes 3, 4, 5, and 7 in wheat, and on chromosomes 6, 3, and 1 in rice. In MQTL10.8, conserved sequences were reported on chromosome 4 in rice and chromosome 2 in wheat. Overall, chromosomes 1, 2, 3, 4, 5, and 7 exhibited conserved regions in wheat for quality and yield-associated traits, with the maximum number of genes on chromosome 7 (Liu et al., 2020; Soriano et al., 2021; Arriagada et al., 2022; Gudi et al., 2022; Khan et al., 2022). In rice, conserved genes were reported on chromosomes 1, 3, 4, 5, 6, and 9, with maximum synteny on chromosomes 4 and 6 (Kulkarni et al., 2020; Khahani, et al., 2021; Zhang et al., 2021; Aloryi et al., 2022; Kumar R. et al., 2022). Therefore, these validated BF-MQTLs, showing synteny relations between different cereal crops, are believed to be promising for improving yield potential and developing high-quality maize varieties through marker-assisted breeding programs.
5 Conclusion
The present study successfully compiled the results of various mapping studies on yield and quality-related traits. A total of 68 MQTLs and 7 singletons were identified in the current study. More than half of the MQTLs were validated through GWAS, and a total of 10 MQTLs were defined as BF-MQTLs based on their CI and PVE values. A comprehensive analysis revealed nearly 60 different genes within the genomic regions of these BF-MQTLs. The putative genes analysed in this study encode essential enzymes involved in multiple metabolic pathways and transcription factors that impact the regulation of key genes. Additionally, some genes were found to be expressed in chloroplasts and peroxisomes, while others function as transporters and signal proteins. Four potential BF-MQTLs were identified to significantly influence both maize yield and grain nutrient quality, making them highly recommended for further utilization in marker-assisted breeding. The overall findings from this study provide valuable insights into robust genomic regions that can be targeted in future investigations to enhance maize yield and nutrient quality.
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.
Author contributions
MeS, VD, and CK collected the data. MeS, VD, MoS, GP, and CK prepared input files. MeS, DS, and AK conducted analysis. MeS, DS, and AS results interpretation. MeS, DS, CK, JS, VD, and DC manuscript writing. DS and DC review and manuscript correction. All authors contributed to the article and approved the submitted version.
Funding
We acknowledge the funds received from the DST funded Women Scientist Project (DST/WOS-A/LS-197/2021).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2023.1248697/full#supplementary-material
SUPPLEMENTARY FIGURE S1 | An overview of number of GWAS validated MQTLs associated with different traits.
References
Amegbor, I. K., Van Biljon, A., Shargie, N., Tarekegne, A., and Labuschagne, M. T. (2022). Heritability and associations among grain yield and quality traits in quality protein maize (QPM) and non-QPM hybrids. Plants 11, 713. doi:10.3390/plants11060713
Akbari, M., Sabouri, H., Sajadi, S. J., Yarahmadi, S., Ahangar, L., Abedi, A., et al. (2022). Mega meta-QTLs: A strategy for the production of golden barley (hordeum vulgare L) tolerant to abiotic stresses. Genes (Basel) 13 (11), 2087. doi:10.3390/genes13112087
Alexandrov, N. N., Brover, V. V., Freidin, S., Troukhan, M. E., Tatarinova, T. V., Zhang, H., et al. (2009). Insights into corn genes derived from large-scale cDNA sequencing. Plant Mol. Biol. 69 (1-2), 179–194. doi:10.1007/s11103-008-9415-4
Aloryi, K. D., Okpala, N. E., Amo, A., Bello, S. F., Akaba, S., and Tian, X. (2022). A meta-quantitative trait loci analysis identified consensus genomic regions and candidate genes associated with grain yield in rice. Front. Plant Sci. 13, 1035851. doi:10.3389/fpls.2022.1035851
Arcade, A., Labourdette, A., Falque, M., Mangin, B., Chardon, F., Charcosset, A., et al. (2004). BioMercator: Integrating genetic maps and QTL towards discovery of candidate genes. Bioinform 20 (14), 2324–2326. doi:10.1093/bioinformatics/bth230
Arriagada, O., Gadaleta, A., Marcotuli, I., Maccaferri, M., Campana, M., Reveco, S., et al. (2022). A comprehensive meta-QTL analysis for yield-related traits of durum wheat (Triticum turgidum L. var. durum) grown under different water regimes. Front. Plant Sci. 13, 984269. doi:10.3389/fpls.2022.984269
Banerjee, S., Roy, P., Nandi, S., and Roy, S. (2023). Advanced biotechnological strategies towards the development of crops with enhanced micronutrient content. Plant Growth Regul. 100, 355–371. doi:10.1007/s10725-023-00968-4
Cao, Y., Zhao, X., Liu, Y., Wang, Y., Wu, W., and Jiang, Y. (2019). Genome-wide identification of ZmHMAs and association of natural variation in ZmHMA2 and ZmHMA3 with leaf cadmium accumulation in maize. Peer. J. 7, e7877. doi:10.7717/peerj.7877
ChardonVirlon, F. B., Moreau, L., Falque, M., Joets, J., and Decousset, L. (2004). Genetic architecture of flowering time in maize as inferred from quantitative trait loci meta-analysis and synteny conservation with the rice genome. Genetics 168 (4), 2169–2185. doi:10.1534/genetics.104.032375
Chen, J., Xu, L., Cai, Y., and Xu, J. (2008). QTL mapping of phosphorus efficiency and relative biologic characteristics in maize (Zea mays L) at two sites. Plant Soil 313 (1-2), 251–266. doi:10.1007/s11104-008-9698-x
Chen, L., An, Y., Li, Y. X., Li, C., Shi, Y., Song, Y., et al. (2017). Candidate loci for yield-related traits in maize revealed by a combination of MetaQTL analysis and regional association mapping. Front. Plant Sci. 8, 2190. doi:10.3389/fpls.2017.02190
Das, A. K., Gowda, M. M., Muthusamy, V., Zunjare, R. U., Chauhan, H. S., Baveja, A., et al. (2021). Development of maize hybrids with enhanced vitamin-E, vitamin-A, lysine, and tryptophan through molecular breeding. Front. Plant Sci. 12, 659381. doi:10.3389/fpls.2021.659381
Dong, Y., Zhang, Z., Shi, Q., Wang, Q., Zhou, Q., and Li, Y. (2015). QTL identification and meta-analysis for kernel composition traits across three generations in popcorn. Euphytica 204, 649–660. doi:10.1007/s10681-015-1360-0
Fang, H., Fu, X., Ge, H., Zhang, A., Shan, T., Wang, Y., et al. (2021). Genetic basis of maize kernel oil-related traits revealed by high-density SNP markers in a recombinant inbred line population. BMC Plant Biol. 21 (1), 344. doi:10.1186/s12870-021-03089-0
Fei, J., Lu, J., Jiang, Q., Liu, Z., Yao, D., Qu, J., et al. (2022). Maize plant architecture trait QTL mapping and candidate gene identification based on multiple environments and double populations. BMC Plant Biol. 22 (1), 110. doi:10.1186/s12870-022-03470-7
Gao, Y., Li, H., DengChen, D. S., Jiang, W., and Chen, J. (2012). Characterization and expression analysis of the maize RING-H2 finger protein gene ZmXERICO responsive to plant hormones and abiotic stresses. Acta Physiol. Plant 34, 1529–1535. doi:10.1007/s11738-012-0951-9
Ghaffari, R., Cannon, E. K. S., Kanizay, L. B., Lawrence, C. J., and Dawe, R. K. (2013). Maize chromosomal knobs are located in gene-dense areas and suppress local recombination. Chromosoma 122 (1-2), 67–75. doi:10.1007/s00412-012-0391-8
Goffinet, B., and Gerber, S. (2000). Quantitative trait loci: A meta-analysis. Genet 155 (1), 463–473. doi:10.1093/genetics/155.1.463
Gudi, S., Saini, D. K., Singh, G., Halladakeri, P., Kumar, P., Shamshad, M., et al. (2022). Unravelling consensus genomic regions associated with quality traits in wheat using meta-analysis of quantitative trait loci. Planta 255 (6), 115. doi:10.1007/s00425-022-03904-4
Gupta, M., Choudhary, M., Singh, A., Sheoran, S., Singla, D., and Rakshit, S. (2023). Meta-QTL analysis for mining of candidate genes and constitutive gene network development for fungal disease resistance in maize (Zea mays L). Crop J. 11 (2), 511–522. doi:10.1016/j.cj.2022.07.020
Halladakeri, P., Gudi, S., Akhtar, S., Singh, G., Saini, D. K., Sakure, A., et al. (2023). Meta-analysis of the quantitative trait loci associated with agronomic traits, fertility restoration, disease resistance, and seed quality traits in pigeonpea (Cajanus cajan L). Plant Genome 2023, e20342. doi:10.1002/tpg2.20342
Jiao, Y., Peluso, P., Shi, J., Liang, T., Stitzer, M. C., Wang, B., et al. (2017). Improved maize reference genome with single-molecule technologies. Nature 546 (7659), 524–527. doi:10.1038/nature22971
Jiao, Z., Wang, J., TianWang, Y. S., Sun, X., and Li, S. (2021). Maize catalases localized in peroxisomes support the replication of maize chlorotic mottle virus. Phytopathol. Res. 3, 17. doi:10.1186/s42483-021-00094-x
Jin, T., Zhou, J., Chen, J., Zhu, L., Zhao, Y., and Huang, Y. (2013). The genetic architecture of zinc and iron content in maize grains as revealed by QTL mapping and meta-analysis. Breed. Sci. 63 (3), 317–324. doi:10.1270/jsbbs.63.317
Jin, Y. N., Cui, Z. H., Ma, K., Yao, J. L., Ruan, Y. Y., and Guo, Z. F. (2021). Characterization of ZmCOLD1, novel GPCR-Type G Protein genes involved in cold stress from Zea mays L. and the evolution analysis with those from other species. Physiol. Mol. Biol. Plants 27 (3), 619–632. doi:10.1007/s12298-021-00966-8
Kadirjan-Kalbach, D. K., Turmo, A., Wang, J., Smith, B. C., Chen, C., Porter, K. J., et al. (2019). Allelic variation in the chloroplast division gene FtsZ2-2 leads to natural variation in chloroplast size. Plant Physiol. 181 (3), 1059–1074. doi:10.1104/pp.19.00841
Karnatam, K. S., Chhabra, G., Saini, D. K., Singh, R., Kaur, G., Praba, U. P., et al. (2023). Genome-wide meta-analysis of QTLs associated with root traits and implications for maize breeding. Int. J. Mol. Sci. 24 (7), 6135. doi:10.3390/ijms24076135
Kaur, C., Sethi, M., Devi, V., Chaudhary, D. P., Phagna, R. K., Singh, A., et al. (2022). Optimization of protein quality assay in normal, opaque-2 and quality protein maize. Front. Sustain. Food Syst. 6, 743019. doi:10.3389/fsufs.2022.743019
Kaur, S., Rakshit, S., Choudhary, M., Das, A. K., and Kumar, R. R. (2021). Meta-analysis of QTLs associated with popping traits in maize (Zea mays L). PLoS One 16 (8), e0256389. doi:10.1371/journal.pone.0256389
Khahani, B., Tavakol, E., Shariati, V., and Rossini, L. (2021). Meta-QTL and ortho-MQTL analyses identified genomic regions controlling rice yield, yield-related traits and root architecture under water deficit conditions. Sci. Rep. 11, 6942–7018. doi:10.1038/s41598-021-86259-2
Khan, H., Krishnappa, G., Kumar, S., Mishra, C. N., Krishna, H., Devate, N. B., et al. (2022). Genome-wide association study for grain yield and component traits in bread wheat (Triticum aestivum L). Front. Genet. 13, 982589. doi:10.3389/fgene.2022.982589
Kulkarni, S. R., Balachandran, S. M., UlaganathanBalakrishnan, K. D., Praveen, M., and Prasad, A. S. (2020). Molecular mapping of QTLs for yield related traits in recombinant inbred line (RIL) population derived from the popular rice hybrid KRH-2 and their validation through SNP genotyping. Sci. Rep. 10 (1), 13695–13721. doi:10.1038/s41598-020-70637-3
Kumar, A. C., Sah, R. P., Muhammed Azharudheen, T. P., Behera, S., Singh, N., Prakash, N. R., et al. (2022a). Understanding complex genetic architecture of rice grain weight through QTL-meta analysis and candidate gene identification. Sci. Rep. 12 (1), 1–13. doi:10.1038/s41598-022-17402-w
Kumar, A., Saini, D. K., Saripalli, G., Sharma, P. K., Balyan, H. S., and Gupta, P. K. (2023a). Meta-QTLs, ortho-meta QTLs and related candidate genes for yield and its component traits under water stress in wheat (Triticum aestivum L). Physiol. Mol. Biol. Plants 29 (4), 525–542. doi:10.1007/s12298-023-01301-z
Kumar, R., Saini, D. K., Kumar, M., Priyanka, V., Akhatar, J., Kaushik, D., et al. (2022b). Revealing the genetic architecture of yield-related and quality traits in Indian mustard [Brassica juncea (L) czern. And coss] using meta-QTL analysis. Agron 12 (10), 2442. doi:10.3390/agronomy12102442
Kumar, S., Saini, D. K., Jan, F., Jan, S., Tahir, M., Djalovic, I., et al. (2023b). Comprehensive meta-QTL analysis for dissecting the genetic architecture of stripe rust resistance in bread wheat. BMC Genom 24 (1), 259–319. doi:10.1186/s12864-023-09336-y
Kumar, S., Singh, V. P., Saini, D. K., Sharma, H., Saripalli, G., Kumar, S., et al. (2021). Meta-QTLs, ortho-MQTLs, and candidate genes for thermotolerance in wheat (Triticum aestivum L). Mol. Breed. 41 (11), 69–22. doi:10.1007/s11032-021-01264-7
Kumari, A., Sharma, D., Sharma, P., Wang, C., Verma, V., Patil, A., et al. (2023). Meta-QTL and haplo-pheno analysis reveal superior haplotype combinations associated with low grain chalkiness under high temperature in rice. Front. Plant Sci. 14, 1133115. doi:10.3389/fpls.2023.1133115
Li, H., Liu, C., Zhou, L., Zhao, Z., Li, Y., Qu, M., et al. (2018). Molecular and functional characterization of the magnesium transporter gene ZmMGT12 in maize. Gene 665, 167–173. doi:10.1016/j.gene.2018.04.068
Li, H., Peng, Z., Yang, X., Wang, W., Fu, J., Wang, J., et al. (2013). Genome-wide association study dissects the genetic architecture of oil biosynthesis in maize kernels. Nat. Genet. 45 (1), 43–50. doi:10.1038/ng.2484
Li, Q., Yang, X., Xu, S., Cai, Y., Zhang, D., Han, Y., et al. (2012). Genome-wide association studies identified three independent polymorphisms associated with α-tocopherol content in maize kernels. PLoS One 7 (5), e36807. doi:10.1371/journal.pone.0036807
Li, S., Li, Y., Cai, Q., Li, X., Sun, Y., Yu, T., et al. (2022). Analysis of the C2H2 gene family in maize (Zea mays L) under cold stress: Identification and expression. Life (Basel) 13 (1), 122. doi:10.3390/life13010122
Li, Z. G. (2020). Regulative role of calcium signaling on methylglyoxal-improved heat tolerance in maize (Zea mays L) seedlings. Plant Signal Behav. 15 (9), 1788303. doi:10.1080/15592324.2020.1788303
Liepman, A. H., Vijayalakshmi, J., Peisach, D., Hulsebus, B., Olsen, L. J., and Saper, M. A. (2019). Crystal structure of photorespiratory alanine: Glyoxylate aminotransferase 1 (AGT1) from Arabidopsis thaliana. Front. Plant Sci. 10, 1229. doi:10.3389/fpls.2019.01229
Liu, H., Mullan, D., Zhang, C., Zhao, S., Li, X., Zhang, A., et al. (2020). Major genomic regions responsible for wheat yield and its components as revealed by meta-QTL and genotype-phenotype association analyses. Planta 252 (4), 65. doi:10.1007/s00425-020-03466-3
Liu, N. A., Xue, Y., Guo, Z., Li, W., and Tang, J. (2016). Genome-wide association study identifies candidate genes for starch content regulation in maize kernels. Front. Plant Sci. 7, 1046. doi:10.3389/fpls.2016.01046
Lo, W. T., Zhang, Y., VadasRoske, O. Y., Gulluni, F., and De Santis, M. C. (2022). Structural basis of phosphatidylinositol 3-kinase C2α function. Nat. Struct. Mol. Biol. 29, 218–228. doi:10.1038/s41594-022-00730-w
Ma, J., and Cao, Y. (2021). Genetic dissection of grain yield of maize and yield-related traits through association mapping and genomic prediction. Front. Plant Sci. 12, 690059. doi:10.3389/fpls.2021.690059
Mackay, T. F. (2001). The genetic architecture of quantitative traits. Annu. Rev. Genet. 35 (1), 303–339. doi:10.1146/annurev.genet.35.102401.090633
Maqbool, M. A., Beshir Issa, A., and Khokhar, E. S. (2021). Quality protein maize (QPM): Importance, genetics, timeline of different events, breeding strategies and varietal adoption. Plant Breed. 140 (3), 375–399. doi:10.1111/pbr.12923
Mohanta, T. K., Mohanta, N., Mohanta, Y. K., and Bae, H. (2015). Genome-wide identification of calcium dependent protein kinase gene family in plant lineage shows presence of novel D-x-D and D-E-L motifs in EF-hand domain. Front. Plant Sci. 6, 1146. doi:10.3389/fpls.2015.01146
Ndlovu, N., Spillane, C., McKeown, P. C., Cairns, J. E., Das, B., and Gowda, M. (2022). Genome-wide association studies of grain yield and quality traits under optimum and low-nitrogen stress in tropical maize (Zea mays L). Theor. Appl. Genet. 135, 4351–4370. doi:10.1007/s00122-022-04224-7
Nunes, R. O., Domiciano, G. A., Alves, W. S., Melo, A. C. A., Nogueira, F. C. S., Canellas, L. P., et al. (2019). Evaluation of the effects of humic acids on maize root architecture by label-free proteomics analysis. Sci. Rep. 9 (1), 12019. doi:10.1038/s41598-019-48509-2
Nuruzzaman, M., Sharoni, A. M., and Kikuchi, S. (2013). Roles of NAC transcription factors in the regulation of biotic and abiotic stress responses in plants. Front. Microbiol. 4, 248. doi:10.3389/fmicb.2013.00248
Ohtani, M. (2018). Plant snRNP biogenesis: A perspective from the nucleolus and cajal bodies. Front. Plant Sci. 8, 2184. doi:10.3389/fpls.2017.02184
Owens, B. F., Lipka, A. E., Magallanes-Lundback, M., Tiede, T., Diepenbrock, C. H., Kandianis, C. B., et al. (2014). A foundation for provitamin A biofortification of maize: Genome-wide association and genomic prediction models of carotenoid levels. Genet 198 (4), 1699–1716. doi:10.1534/genetics.114.169979
Prasanna, B. M., Palacios-Rojas, N., Hossain, F., Muthusamy, V., Menkir, A., Dhliwayo, T., et al. (2020). Molecular breeding for nutritionally enriched maize: Status and prospects. Front. Genet. 10, 1392. doi:10.3389/fgene.2019.01392
Prasanna, B. M., Vasal, S. K., Kassahun, B., and Singh, N. N. (2001). Quality protein maize. Curr. Sci. 81, 1308–1319.
Qu, Z., Wu, Y., Hu, D., Li, T., Liang, H., Ye, F., et al. (2022). Genome-wide association analysis for candidate genes contributing to kernel-related traits in maize. Front. Plant Sci. 13, 872292. doi:10.3389/fpls.2022.872292
Ren, H., and Gray, W. M. (2015). SAUR proteins as effectors of hormonal and environmental signals in plant growth. Mol. Plants. 8, 1153–1164. doi:10.1016/j.molp.2015.05.003
Saini, D. K., Chopra, Y., Pal, N., Chahal, A., Srivastava, P., and Gupta, P. K. (2021). Meta-QTLs, ortho-MQTLs and candidate genes for nitrogen use efficiency and root system architecture in bread wheat (Triticum aestivum L). Physiol. Mol. Biol. Plants 27 (10), 2245–2267. doi:10.1007/s12298-021-01085-0
Saini, D. K., Srivastava, P., Pal, N., and Gupta, P. K. (2022). Meta-QTLs, ortho-meta-QTLs and candidate genes for grain yield and associated traits in wheat (Triticum aestivum L). Theor. App. Genet. 135, 1049–1081. doi:10.1007/s00122-021-04018-3
Sandhu, N., Pruthi, G., Raigar, O. P., Singh, M. P., Phagna, K., Kumar, A., et al. (2021). MetaQTL analysis in rice and cross-genome talk of the genomic regions controlling nitrogen use efficiency in cereal crops revealing phylogenetic relationship. Front. Genet. 12, 807210. doi:10.3389/fgene.2021.807210
Semagn, K., Beyene, Y., Warburton, M. L., Tarekegne, A., Mugo, S., Meisel, B., et al. (2013). Meta-analyses of QTL for grain yield and anthesis silking interval in 18 maize populations evaluated under water-stressed and well-watered environments. BMC Genom 14 (1), 313. doi:10.1186/1471-2164-14-313
Shafi, S., Saini, D. K., Khan, M. A., Bawa, V., Choudhary, N., Dar, W. A., et al. (2022). Delineating meta-quantitative trait loci for anthracnose resistance in common bean (Phaseolus vulgaris L). Front. Plant Sci. 13, 966339. doi:10.3389/fpls.2022.966339
Sheoran, S., Gupta, M., Kumari, S., Kumar, S., and Rakshit, S. (2022). Meta-QTL analysis and candidate genes identification for various abiotic stresses in maize (Zea mays L) and their implications in breeding programs. Mol. Breed. 42 (5), 26. doi:10.1007/s11032-022-01294-9
Soriano, J. M., Colasuonno, P., Marcotuli, I., and Gadaleta, A. (2021). Meta-QTL analysis and identification of candidate genes for quality, abiotic and biotic stress in durum wheat. Sci. Rep. 11 (1), 11877–11915. doi:10.1038/s41598-021-91446-2
Suwarno, W. B., Pixley, K. V., Palacios-Rojas, N., Kaeppler, S. M., and Babu, R. (2015). Genome-wide association analysis reveals new targets for carotenoid biofortification in maize. Theor. Appl. Genet. 128 (5), 851–864. doi:10.1007/s00122-015-2475-3
Tanin, M. J., Saini, D. K., Sandhu, K. S., Pal, N., Gudi, S., Chaudhary, J., et al. (2022). Consensus genomic regions associated with multiple abiotic stress tolerance in wheat and implications for wheat breeding. Sci. Rep. 12 (1), 13680–13717. doi:10.1038/s41598-022-18149-0
Tanumihardjo, S., McCulley, L., Roh, R., Lopez-Ridaura, S., Palacios-Rojas, N., and Gunaratna, N. (2019). Maize agro-food systems to ensure food and nutrition security in reference to the sustainable development goals. Glob. Food Sec 25, 100327. doi:10.1016/j.gfs.2019.100327
Venado, R. E., Owens, B. F., Ortiz, D., Lawson, T., Mateos-Hernandez, M., Ferruzzi, M. G., et al. (2017). Genetic analysis of provitamin A carotenoid β-cryptoxanthin concentration and relationship with other carotenoids in maize grain (Zea mays L). Mol. Breed. 37, 127. doi:10.1007/s11032-017-0723-8
Veyrieras, J. B., Goffinet, B., and Charcosset, A. (2007). MetaQTL: A package of new computational methods for the meta-analysis of QTL mapping experiments. BMC Bioinform 8, 49. doi:10.1186/1471-2105-8-49
Waititu, J. K., Cai, Q., Sun, Y., Sun, Y., Li, C., Zhang, C., et al. (2021). Transcriptome profiling of maize (Zea mays L) leaves reveals key cold-responsive genes, transcription factors, and metabolic pathways regulating cold stress tolerance at the seedling stage. Genes (Basel) 12 (10), 1638. doi:10.3390/genes12101638
Wang, B., Tseng, E., Regulski, M., Clark, T. A., Hon, T., Jiao, Y., et al. (2016a). Unveiling the complexity of the maize transcriptome by single-molecule long-read sequencing. Nat. Commun. 7 (1), 11708–11713. doi:10.1038/ncomms11708
Wang, J., Pan, W., Nikiforov, A., King, W., Hong, W., Li, W., et al. (2021a). Identification of two glycerophosphodiester phosphodiesterase genes in maize leaf phosphorus remobilization. Crop J. 9 (1), 95–108. doi:10.1016/j.cj.2020.05.004
Wang, W., Ren, Z., Li, L., Du, Y., Zhou, Y., Zhang, M., et al. (2022). Meta-QTL analysis explores the key genes, especially hormone related genes, involved in the regulation of grain water content and grain dehydration rate in maize. BMC Plant Biol. 22, 346. doi:10.1186/s12870-022-03738-y
Wang, X., An, Y., Xu, P., and Xiao, J. (2021b). Functioning of PPR proteins in organelle RNA metabolism and chloroplast biogenesis. Front. Plant Sci. 12, 627501. doi:10.3389/fpls.2021.627501
Wang, Y., Deng, D., Bian, Y., Lv, Y., and Xie, Q. (2010). Genome-wide analysis of primary auxin-responsive Aux/IAA gene family in maize (Zea mays. L). Mol. Biol. Rep. 8, 3991–4001. doi:10.1007/s11033-010-0058-6
Wang, Y., Wang, Y., Wang, X., and Deng, D. (2020). Integrated meta-QTL and genome-wide association study analyses reveal candidate genes for maize yield. J. Plant Growth Reg. 39, 229–238. doi:10.1007/s00344-019-09977-y
Wang, Y., Xu, J., Deng, D., Ding, H., Bian, Y., Yin, Z., et al. (2016b). A comprehensive meta-analysis of plant morphology, yield, stay-green, and virus disease resistance QTL in maize (Zea mays L). Planta 243 (2), 459–471. doi:10.1007/s00425-015-2419-9
Wu, G., Zhao, Y., Shen, R., Wang, B., Xie, Y., Ma, X., et al. (2019). Characterization of maize phytochrome-interacting factors in light signaling and photomorphogenesis. Plant Physiol. 181 (2), 789–803. doi:10.1104/pp.19.00239
Wu, X., Cai, G., Gong, F., An, S., Cresti, M., and Wang, W. (2015). Proteome profiling of maize pollen coats reveals novel protein components. Plant Mol. Biol. Rep. 33, 975–986. doi:10.1007/s11105-014-0812-3
Xu, J., Qin, X., Ni, Z., Chen, F., Fu, X., and Yu, F. (2022). Identification of zinc efficiency-associated loci (ZEAL s) and candidate genes for zn deficiency tolerance of two recombination inbred line populations in maize. Int. J. Mol. Sci. 23 (9), 4852. doi:10.3390/ijms23094852
Yan, B., Xu, X., Gu, Y., Zhao, Y., Zhao, X., He, L., et al. (2018). Genome-wide characterization and expression profiling of diacylglycerol acyltransferase genes from maize. Genome 61 (10), 735–743. doi:10.1139/gen-2018-0029
Yan, P., Zhu, Y., Wang, Y., Ma, F., Lan, D., Niu, F., et al. (2022). A new RING finger protein, PLANT ARCHITECTURE and GRAIN NUMBER 1, affects plant architecture and grain yield in rice. Int. J. Mol. Sci. 23, 824. doi:10.3390/ijms23020824
Yang, G., Dong, Y., Li, Y., Wang, Q., Shi, Q., and Zhou, Q. (2013). Verification of QTL for grain starch content and its genetic correlation with oil content using two connected RIL populations in high-oil maize. PLoS One 8 (1), e53770. doi:10.1371/journal.pone.0053770
Zdunic, Z., Grljusic, S., Ledencan, T., Duvnjak, T., and Simic, D. (2014). Quantitative trait loci mapping of metal concentrations in leaves of the maize IBM population. Hereditas 151 (2-3), 55–60. doi:10.1111/hrd2.00048
Zeng, T., Meng, Z., Yue, R., Lu, S., Li, W., Li, W., et al. (2022). Genome wide association analysis for yield-related traits in maize. BMC Plant Biol. 22 (1), 449. doi:10.1186/s12870-022-03812-5
Zhang, G., Wang, R., MaGao, J. H., Deng, L., and Wang, N. (2021). Genome-wide association studies of yield-related traits in high-latitude japonica rice. BMC Genom. Data 22, 39. doi:10.1186/s12863-021-00995-y
Zhang, J., Wu, S., Boehlein, S. K., McCarty, D. R., Song, G., Walley, J. W., et al. (2019). Maize defective kernel5 is a bacterial TamB homologue required for chloroplast envelope biogenesis. J. Cell Biol. 218 (8), 2638–2658. doi:10.1083/jcb.201807166
Zhang, S., Ding, G., He, W., Liu, K., Luo, Y., Tang, J., et al. (2020). Functional characterization of the 1-deoxy-D-xylulose 5-phosphate synthase genes in Morus notabilis. Front. Plant Sci. 11, 1142. doi:10.3389/fpls.2020.01142
Zhang, X., Ren, Z., Luo, B., Zhong, H., Ma, P., Zhang, H., et al. (2022). Cramér-Rao bound-informed training of neural networks for quantitative MRI. Crop J. 10 (2), 436–448. doi:10.1002/mrm.29206
Zhang, Z., Liu, Z., Hu, Y., Li, W., Fu, Z., Ding, D., et al. (2014). QTL analysis of kernel-related traits in maize using an immortalized F2 population. PLoS One 9 (2), e89645. doi:10.1371/journal.pone.0089645
Zheng, Y., Yuan, F., Huang, Y., Zhao, Y., Jia, X., Zhu, L., et al. (2021). Genome-wide association studies of grain quality traits in maize. Sci. Rep. 11 (1), 9797–9812. doi:10.1038/s41598-021-89276-3
Zhou, Z., Li, G., Tan, S., Li, D., Weiß, T. M., Wang, X., et al. (2020). A QTL atlas for grain yield and its component traits in maize (Zea mays). Plant Breed. 139 (14), 562–574. doi:10.1111/pbr.12809
Keywords: maize, yield, quality, meta-QTLs, breeder-friendly, candidate genes
Citation: Sethi M, Saini DK, Devi V, Kaur C, Singh MP, Singh J, Pruthi G, Kaur A, Singh A and Chaudhary DP (2023) Unravelling the genetic framework associated with grain quality and yield-related traits in maize (Zea mays L.). Front. Genet. 14:1248697. doi: 10.3389/fgene.2023.1248697
Received: 27 June 2023; Accepted: 26 July 2023;
Published: 07 August 2023.
Edited by:
Ertugrul Filiz, Duzce University, TürkiyeReviewed by:
Jiban Shrestha, Nepal Agricultural Research Council, NepalShulei Guo, Henan Academy of Agricultural Sciences, China
Copyright © 2023 Sethi, Saini, Devi, Kaur, Singh, Singh, Pruthi, Kaur, Singh and Chaudhary. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Dharam Paul Chaudhary, chaudharydp@gmail.com
†Present Address: Dinesh Kumar Saini, Department of Plant and Soil Science, Texas Tech University, Lubbock, TX, United States
‡These authors have contributed equally to this work