Nikwan Shariatipour
Bahram Heidari
Ahmad Tahmasebi
Christopher Richards2†- 1Department of Plant Production and Genetics, School of Agriculture, Shiraz University, Shiraz, Iran
- 2USDA ARS National Laboratory for Genetic Resources Preservation, Fort Collins, CO, United States
Comparative genomics and meta-quantitative trait loci (MQTLs) analysis are important tools for the identification of reliable and stable QTLs and functional genes controlling quantitative traits. We conducted a meta-analysis to identify the most stable QTLs for grain yield (GY), grain quality traits, and micronutrient contents in wheat. A total of 735 QTLs retrieved from 27 independent mapping populations reported in the last 13 years were used for the meta-analysis. The results showed that 449 QTLs were successfully projected onto the genetic consensus map which condensed to 100 MQTLs distributed on wheat chromosomes. This consolidation of MQTLs resulted in a three-fold reduction in the confidence interval (CI) compared with the CI for the initial QTLs. Projection of QTLs revealed that the majority of QTLs and MQTLs were in the non-telomeric regions of chromosomes. The majority of micronutrient MQTLs were located on the A and D genomes. The QTLs of thousand kernel weight (TKW) were frequently associated with QTLs for GY and grain protein content (GPC) with co-localization occurring at 55 and 63%, respectively. The co- localization of QTLs for GY and grain Fe was found to be 52% and for QTLs of grain Fe and Zn, it was found to be 66%. The genomic collinearity within Poaceae allowed us to identify 16 orthologous MQTLs (OrMQTLs) in wheat, rice, and maize. Annotation of promising candidate genes (CGs) located in the genomic intervals of the stable MQTLs indicated that several CGs (e.g., TraesCS2A02G141400, TraesCS3B02G040900, TraesCS4D02G323700, TraesCS3B02G077100, and TraesCS4D02G290900) had effects on micronutrients contents, yield, and yield-related traits. The mapping refinements leading to the identification of these CGs provide an opportunity to understand the genetic mechanisms driving quantitative variation for these traits and apply this information for crop improvement programs.
Introduction
Most agricultural systems focus on increasing crop productivity and grain yield (GY) and fewer efforts have been devoted to the grain yield–quality tradeoff. However, a shift from prioritizing yield to more emphasis on quality, such as nutrient content is gaining ground in breeding programs (Khush et al., 2012). Extending the existing concepts for a simultaneous selection of GY, quality traits, and micronutrient contents seems necessary to facilitate the development of varieties with an effective combination of yield potential and end-use quality (Michel et al., 2019). A rapid increase in micronutrient deficiency in food grains has resulted in micronutrient malnutrition among consumers. Fe and Zn deficiencies are serious and prevalent sources of malnutrition in developing countries with high consumption of cereals, such as wheat (Black et al., 2013; Shahzad et al., 2014; Kumar et al., 2019). Diets based on staple food crops with Zn and Fe deficiency have been widely recognized as a major global health problem that affects almost three billion people (Murray and Lopez, 2013). Breeding crops through biofortification is a practical approach to cope with Fe and Zn deficiencies by increasing the grain Fe content (GFeC) and grain Zn content (GZnC) within the edible parts of staple food crops, especially cereals (Stein, 2010; Liu et al., 2019; Shariatipour and Heidari, 2020).
Wheat biofortification through breeding methods is a promising strategy to ameliorate Fe and Zn deficiencies in developing countries (Liu et al., 2019). One of the most important challenges in breeding for micronutrients are negative genetic trade-offs between yield and micronutrient traits (Flatt and Heyland, 2011; Fabian and Flatt, 2012). At the genetic level, such trade-offs are thought to be caused by alleles with antagonistic pleiotropic effects or by linkage disequilibrium between loci (Fabian and Flatt, 2012). While GFeC and GZnC biofortification is an important objective in wheat breeding programs, other important traits, such as GY and grain protein content (GPC) typically cannot be compromised. The wheat quality is measured by its rheological traits, such as GPC (Goel et al., 2019) since wheat is a major source of protein accounting for 19% of human protein intake in the developing countries (Braun et al., 2010). The genetic control of quality traits and GY is complex (Zilic et al., 2011; Velu et al., 2018; Giancaspro et al., 2019; Liu et al., 2019). A high priority of breeders is to develop high-quality genotypes that balance acceptable yield potential while maintaining quality characteristics, both of which are highly dependent on the co-variance between GY and the major quality (i.e., GZnC, GFeC, and GPC) criteria (Michel et al., 2019). However, a negative correlation between the GY and quality traits is challenging (Simmonds, 1995; Velu et al., 2016; Michel et al., 2019). In addition to this, our research indicates that protein concentration in grain decreases under elevated air CO2 concentrations of 550 μmol/mol (Fernando et al., 2012). By the end of the twenty-first century, it is predicted that the global temperature will rise from 1.1 to −3.1°C and the atmospheric CO2 concentration will reach above 550 ppm under intermediate scenarios (Pachauri et al., 2014) and that heat waves will occur with a higher frequency and longer duration (Pachauri et al., 2014). Given these future climate scenarios, it is critical to anticipate the effects of future growing environments and focus on breeding strategies that compensate for changes in grain quality. Likewise, it is important to have an inclusive breeding objective that tracks a portfolio of indirect traits responsible for grain quality and productivity, such as the assessment of dry matter accumulation, photosynthesis, coleoptile growth, carbon isotope discrimination, plant senescence, and rheological properties (Rebetzke et al., 2008; Liang et al., 2010; Vijayalakshmi et al., 2010; Goel et al., 2019).
Because of the large genome size and limited genome sequence information in wheat, the typical mapping intervals are quite large in most studies especially for the complex quality traits (Li Q. et al., 2020) and so further refinement is needed to narrow down the QTL intervals. Developing a statistically derived catalog of relevant loci is critical for developing marker-assisted selection (MAS) approaches in breeding programs. These markers can be applied to the quantitative trait loci (QTLs) that regulate the accumulation of high mineral nutrient concentration in grain along with QTLs for GY and grain quality traits. The QTL mapping method involves creating a QTL continuity map to identify genomic regions associated with quantitative traits (Mohan et al., 1997). Although QTL mapping is a powerful approach for detecting the genomic regions associated with complex traits, the genetic effects of QTL identified in different studies may not be present or are simply not tested in different genetic backgrounds and environments (Zhang L. Y. et al., 2010). In addition, the number of traits that can be measured in any single study is always resource limited (Acuña-Galindo et al., 2015). Overall, biparental populations are strongly influenced by different factors consisting of parents, the size and type of population, the choice of marker sets, and environmental conditions (Li et al., 2013; Izquierdo et al., 2018; Lei et al., 2018; Zhao et al., 2018).
In the last decade, an efficient approach called meta-QTL (MQTL) analysis has emerged in order to circumvent these restrictions. The MQTL method was initially developed by Goffinet and Gerber (2000) and was then improved by Veyrieras et al. (2007) is a method that gathers QTL data from independent experiments, years, location, and genetic backgrounds to detect stable QTLs (Goffinet and Gerber, 2000; Arcade et al., 2004; Hanocq et al., 2007; Sosnowski et al., 2012). The meta-QTL analysis integrates the information of QTLs from different population types and sizes identified in different environmental conditions to find stable MQTLs in a narrower genomic region with small CI (Goffinet and Gerber, 2000; Hanocq et al., 2007; Li et al., 2013).
The MQTL analysis allows for the dissection of genetic correlation among different traits (Truntzler et al., 2010; Danan et al., 2011; Xiang et al., 2012; Badji et al., 2018; Delfino et al., 2019). Hence, the MQTL analysis helps to enquire co-location of QTLs relying on dense marker maps that are responsible for different desirable traits including micronutrient contents, GY, and quality traits (Delfino et al., 2019). Currently, the MQTL analysis has become popular research in most studies (Goffinet and Gerber, 2000) related to micronutrients content (Raza et al., 2019), yield (Zhang L. Y. et al., 2010; Avni et al., 2018), and quality traits (Quraishi et al., 2017) to overcome the inconsistent QTL information reported for GZnC, GFeC, GPC, and yield traits. Further, MQTL analysis helps to identify the candidate genes (CGs) and refine genomic regions for yield and quality traits (Raza et al., 2019). However, little is known about the relation of micronutrients and quality-related MQTLs with agronomic traits in wheat.
The transferability of QTLs between cereals based on the analysis of syntenic regions and genomic collinearity helps to identify stable and important QTLs for use in breeding programs. The aims of this study were to perform a QTL meta-analysis to (1) identify QTLs that are consistently associated with grain quality, yield traits, and micronutrient content (2) explore the co-localized QTLs controlling GY, GPC, GZnC, and GFeC in the wheat genome, and (3) assess the transferability of QTLs between wheat, rice, and maize based on the comparative genomics and the orthologous MQTL (OrMQTS) mining. The outcome of this study will aid plant breeders in refining micronutrients, GY, and quality traits for crop improvement through marker-assisted breeding. Refining our understanding of the genetic architecture of micronutrients, GY, and quality traits often leads to the interrelationship between the regions of the genome that may be more challenging to breed independently. The MQTL is an analytic procedure that helps refine these relationships between CGs by providing for precision and statistical power.
Materials and Methods
QTL Database Development
A database consisted of 735 quantitative trait loci (QTLs) derived from 27 independent mapping populations (assessed between 2006 and 2019) assigned to 70 traits (Table 1) was used for the meta-QTL (MQTL) analysis. The independent populations consisted of 20 recombinant inbred lines (RILs) and 7 double haploids (DH) populations, with population sizes ranging from 92 to 485 lines (Supplementary Table 1). The reported position, the proportion of the phenotypic variance (R2), and the logarithm of odds (LOD score) of the initial QTLs were used for the analysis of meta-QTLs. For QTLs with missing LOD or R2, the values were estimated by the following equation (Nagelkerke, 1991):
where n represents the size of the population.
Table 1. The list of assessed traits in meta-quantitative trait loci (MQTL) analysis.
Constructing Consensus Genetic Map and QTL Projection
The data files of the 27 maps were integrated with the Somers (Somers et al., 2004) reference map for the construction of a consensus genetic map. Attempts to use other mapping studies consisting of single-nucleotide polymorphism (SNP) were unsuccessful due to the lack of SNP density in the regions for the meta-QTL analysis. The constructed map file for each population consisted of information on cross-type, population size, map function, map units, and the position of different markers in different linkage groups. The individual QTLs derived from independent populations were projected onto the consensus genetic map consisted of 3,394 markers with a total length of 3,412.5 cM.
Meta-QTL Analysis
Meta-QTL analysis was performed in BioMercator v4.2 (Arcade et al., 2004; Sosnowski et al., 2012). For n QTLs, the BioMercator tests the most likely assumption based on Akaike information criterion (AIC), corrected AIC (AICc), AIC 3 candidate models (AIC3), Bayesian information criterion (BIC), and an average weight of evidence (AWE) criteria in which the prevalent value among them was considered as the best fit. The consensus QTL from the optimal model was reported as MQTL. Consequently, the MQTL position and distribution on each linkage group were presented as a heatmap using pheatmap R package (Kolde, 2013). Moreover, the initial QTLs with 95% confidence interval (CI), QTL density in the identified MQTL, and the distribution of MQTLs were drawn on the linkage groups using shinyCircos web tool based on the R program (Yu et al., 2017). The variation of QTL density for different traits toward centromeric and telomeric genomic regions was estimated following the approach by Martinez et al. (2016). The QTL density was determined by counting the number of QTLs for each trait on 50 cM intervals across the wheat genome, starting from the centromere region of a linkage group at position 0. The centromere position was retrieved from the study by Wan et al. (2017).
Functional Candidate Genes in MQTLs Intervals
The MQTLs containing more than five trait-QTLs from different experiments were considered as the most stable consensus regions and were analyzed for the detection of functional candidate genes (CGs). To identify the functional CGs, the sequences of the flanking markers for each MQTL were retrieved from “Grain Genes” database (https://wheat.pw.usda.gov/browse?class=probe;query=BARC%2A;begin=351) for the simple sequence repeat (SSR) flanking markers and “Diversity Array Technology” (https://www.diversityarrays.com/) (DArT) flanking markers. For flanking markers lacking a definite position on the wheat genome, the closest markers on the genetic consensus map were selected to determine the MQTL position. Additionally, for those flanking markers lacking sequence information in databases, the forward and reverse sequences were retrieved from the “Grain Genes” database and were used for the Basic Local Alignment Search Tool (BLAST) analysis against the newest wheat reference genome (IWGSC RefSeq v2.0) for detecting the genomic position of each MQTL. The annotation and gene ontology (GO) of genes lying at the MQTL interval were retrieved from EnsemblPlants (http://plants.ensembl.org/index.html) using the new wheat genome (IWGSC v2.0). Finally, the orthologous of genes located at each MQTL interval were investigated in rice to describe the functional CGs based on their reported functions in wheat or rice.
Identification of Traits Within MQTLs
To analyze traits within the MQTL regions, the MQTL results were converted into binary scores (0 or 1) on the basis of the absence/presence of an individual trait-QTL within an MQTL region. We tabulated the number of times a trait was present within an MQTL, the number of QTLs for a trait present within an MQTL (implying confirmation of the QTL), and the number of times the traits were able to be co-localized within an MQTL. A chi-squared test with one degree of freedom was performed to determine traits showing significant co-localization with grain protein content (GPC), grain zinc content (GZnC), grain Fe content (GFeC), and grain yield (GY) beyond what would be expected for a random distribution of QTL within MQTL throughout the genome. The expected number of MQTL associated with a trait and each GPC, GZnC, GFeC, and GY were separately calculated by multiplying the number of observed MQTL for a trait by the proportion of MQTL containing GPC, GZnC, GFeC, and GY QTL(s). Traits within the MQTL regions were also analyzed using IBM SPSS Statistics v.24. The simple regression analysis was performed using Minitab v. 18 to determine the effect of MQTLs on the association of GY, GPC, GFeC, and GZnC.
MQTLs and GWAS Comparison
The detected MQTLs were compared with the significant loci associated with different quantitative traits identified in wheat genome-wide association studies (GWAS). The mapped coordinates of the identified significant loci through GWAS were compared to those found with the MQTL analysis.
Orthologous MQTL
Due to the high synteny among genes in Poacaea, the most stable and promising wheat MQTLs were evaluated for the detection of the orthologous MQTLs (OrMQTL) in rice (Lei et al., 2018; Raza et al., 2019; Khahani et al., 2020), and maize (Semagn et al., 2013; Wang et al., 2013, 2016). A set of orthologous genes at the MQTLs regions was considered as a criterion of a syntenic region using EnsemblPlants (http://plants.ensembl.org/index.htmldatabase).
Results
Genomic Quantitative Trait Loci Distributions
The individual traits of quantitative trait loci (QTLs) used in this study are listed in Supplementary Table 2. The QTLs for thousand kernel weight (TKW) (18.03%), grain yield (GY), (13.11%), and a number of grains per spike (NS) (13.11%) were the most frequently reported agronomic QTLs identified in the tested mapping populations. The grain Fe content (GFeC) (33.33 %) and grain Zn content (GZnC) (28.21 %) of QTLs for micronutrient traits and grain protein content (GPC) (64.00 %) and sudden death syndrome (SDS) (10.00 %) of QTLs for quality traits were frequent.
Our results suggest a non-random distribution of QTLs within the wheat genome. Distribution of QTLs on the basis of physical size [χ2 (2) = 60.17, P = 8.58E−14] showed that 254, 326, and 155 QTLs were located on the A, B, and D genomes, respectively. The QTL distribution was significantly different among the seven chromosome groups [ = 47.12, P = 1.77E−8], ranging from as few as 76 QTLs on group 6 to as many as 165 QTLs on Group 2. Chromosome 3B with 75 QTLs had the highest number of QTLs, followed by chromosome 2B (62 QTLs) and 2A (61 QTLs), while chromosome 3A with 10 QTLs had the lowest QTL. The distribution of QTL over the genetic linkage map with respect to centromeric and telomeric regions was distinctly non-random (Figure 1). The non-telomeric region of each chromosome (−50 up to + 50 cM intervals) had the highest number of QTLs (Figure 1). We did not detect QTLs at 100 cM for the tested traits.
Figure 1. Distribution of quantitative trait loci (QTLs) for traits on all chromosomes represented as the number of QTLs per distance (50 cM), starting from the centromeric region of each chromosome where it was considered at the position 0 cM. Each dot represents the exact location of each QTL.
The distribution of meta QTLs (MQTLs) indicated that a cluster of MQTLs was mapped to the non-telomeric regions. Besides, MQTL_5B_4 and MQTL_6B_1 were located near the centromeric region of chromosomes 5B and 6B, respectively (Figure 2). There was a significant correlation between the number of initial QTLs and MQTLs (r = 0.46, P < 0.03). The number of MQTLs per chromosome varied from two (chromosomes 3A and 5D) to eight (chromosomes 1B, 2A, and 6B).
Figure 2. Distribution of meta-quantitative trait loci (MQTLs) for assessed traits on all chromosomes represented as a number of MQTLs per distance (50 cM), starting from the centromeric region of each chromosome where it was considered at position 0 cM. Each dot represents the exact location of each MQTL.
Meta-QTL Analysis
Of the 735 initial QTLs, 449 were successfully projected onto the genetics consensus map and used in the meta-QTL analysis (Figure 3). A total of 100 MQTLs were detected and the number of individual QTL per MQTL ranged from 1 to 43 (Figures 4, 5; Supplementary Figure 1). The number of traits present per MQTL region ranged from 1 to 18. Among the identified MQTLs, MQTL_3B_1 that contained 43 QTLs had the highest number of initial QTLs followed by MQTL_7A_3 with 29 initial QTLs (Supplementary Table 3). These two MQTLs can be considered as the most stable QTLs under different experimental conditions. The detailed information of MQTLs consisted of the chromosome number, position, confidence interval (CI), flanking markers, and traits which are shown in Table 2. The higher marker density of the consensus map compared with the lower marker density in the independent linkage maps helped to reduce the CI of QTLs up to three-fold with an average of 4.63 cM in MQTLs compared with the mean CI of 13.73 cM for the original QTLs. Among the detected MQTLs, the CI of 11 MQTLs was reduced up to <1 cM (Table 2).
Figure 3. Distribution of projected quantitative trait loci (QTLs) across the wheat (Triticum aestivum L.) chromosomes. The numbers inside each parenthesis represent the number of QTLs.
Figure 4. Position of detected meta-quantitative trait loci (MQTLs) on the wheat genome associated with micronutrient content, grain quality, and quantitative traits with 95% confidence interval (CI). Each color in a different linkage group indicates the number of initial QTLs involved in each MQTL. The flanking markers for each MQTL are presented on the left side of the linkage groups.
Figure 5. Circus plot showing distribution of QTLs and MQTLs on 21 linkage group of wheat. The outermost circle indicates the length of each linkage map on the consensus genetic map. The second circle indicates the initial QTLs with 95% confidence intervals. The third circle displays the density of QTLs in each MQTLs. The fourth circle displays the position of detected MQTL related to assessed traits based on heatmap illustration. The last inner circle demonstrated common markers links between different genome (A, B, and D) of each linkage groups.
Table 2. Description of detected meta-quantitative trait loci (MQTL).
Functional Identification of Candidate Genes
The genomic positions of the stable MQTLs and the number of functional candidate genes (CGs) in their intervals are reported in Table 3. The range for the number of the CGs annotated in the tested meta-QTLs was between 20 and 802. The MQTL_4D_1, MQTL_4B_3, MQTL_3B_1, and MQTL_5A_4 harbored the greatest number of CGs. Among the detected CGs on MQTL regions, several well-known genes including psbL (21.896318-21.896434 Mb), psbT (21.905522-21.905638 Mb), rpl33 (21.899696-21.899896 Mb), and rps4 (24.160684-24.161289 Mb) genes were located in the MQTL_3D_4 region. In addition, the miR166 gene was detected in the MQTL_4D_1 (450.196302-450.196403 Mb) and MQTL_7A_4 (561.152174-561.152339 Mb) regions. Furthermore, several CGs with unknown annotation in wheat and were orthologous to genes in rice were identified (Supplementary Tables 4, 5). Overall, some CGs, such as TraesCS2A02G141400, TraesCS3B02G040900, TraesCS4D02G323700, TraesCS3B02G077100, and TraesCS4D02G290900 were uncovered with a possible role in micronutrient contents, yield, and yield-related traits.
Table 3. The genomic position of most stable meta-quantitative trait loci (MQTLs) on the wheat genome and number of genes in their genomic intervals.
Traits Analysis Within MQTLs
Our data show the unequal contribution of individual QTL across the detected MQTLs (Figure 6 and Supplementary Table 3). Individual QTL for TKW were present in 41 of the 100 MQTL regions, the most for any agronomic trait, followed by GY (38 MQTL). Among quality traits, QTLs for GPC was the most distributed QTLs which were located among 19 MQTLs. The GFeC and GZnC QTLs were presented in 19 and 12 MQTLs, respectively (Supplementary Table 3).
Figure 6. Distribution of QTLs controlling different traits in detected MQTL regions.
Analysis for the co-localization of QTLs revealed that QTLs for GY and TKW were frequently co-localized with QTLs of the target traits (GPC, GFeC, and GZnC). The QTLs for TKW showed 52% co-localizations with QTLs for GPC (Supplementary Table 6). The GY QTLs showed 52% co-localization with the QTLs for GFeC. Results also indicated 66% co-localization between grain Fe and grain Zn QTLs. Co-localization frequency for GY (R2 = 80.81%) was strongly associated with the total number of MQTL for a trait. The association of co-localization frequency with the overall number of MQTL for a trait was relatively strong for GPC (R2 = 65.58%), GFeC (R2 = 58.74%), and GZnC (R2 = 56.00%). For most traits, association with target traits (GPC, GZnC, GFeC, and GY) did not differ from the expected on the basis of chi-squared analysis. However, the association of traits with GY, GPC, GFeC, and GZnC as target traits was more than expected (Supplementary Table 6).
Comparison of the Identified MQTLs and QTL Mapping in the Wheat GWAS
The comparison of the MQTL locations with genome-wide association studies (GWAS) QTL regions showed that 21 significant signals (SNPs-linked QTLs) of the available wheat GWAS map were co-located with MQTLs of seven of the 35 traits tested in our study (Table 4). The results indicated the co-localization of significant single-nucleotide polymorphisms (SNPs) for GY (6 SNPs), a number of grains per spike (NG) (1 SNPs), plant height (PH) (4 SNPs), spike length (SL) (2 SNPs), spike number (SN) (1 SNPs), spike number per spike (SNS) (4 SNPs), and thousand kernel weight (TKW) (3 SNPs) traits in the wheat GWAS with the identified MQTLs in our study. For instance, the MQTL_4A_3, MQTL_4A_4, and MQTL_4B_3 identified for TKW in our study were positioned in the genomic regions of the major signal reported for TKW in the wheat GWSA map. Overall, the co-located MQTLs and significant GWAS signals were distributed on chromosomes 1B, 2B, 3B, 4A, 4B, 4D, 5A, 7A, and 7B (Table 4).
Table 4. The collinear meta-quantitative trait loci (MQTLs) with the significant loci in wheat genome-wide association studies.
Orthologous MQTL Mining of Wheat in Rice and Maize
The comparative analysis for QTLs in wheat, rice, and maize resulted in the identification of orthologous MQTL (OrMQTL). Nine OrMQTLs were detected for wheat and rice including five OrMQTL for GY and two for PH, and GFeC/GZnC, respectively. Moreover, seven OrMQTL were identified in wheat and maize consisting of six and one OrMQTL for GY and PH, respectively. Among the uncovered OrMQTLs, the OrMQTL_10 was a cross-species QTL in wheat, rice, and maize (Figure 7; Table 5). The MQTL_7B_2, MQTL_3B_4, MQTL_4D_1, and MQTL_5A_4 for GY in wheat were in the co-linear regions for rice GY MQTLs on chromosome 6 (MQTL6-2), 1 (MQTL-YLD3), 3 (MQTL-YLD9), and 11 (MQTL-YLD19), respectively (Table 5). The wheat MQTL_4B_1 and MQTL_7B_5 were in the co-linear region of MQTLs for PH on chromosome 3 (MQTL-PH11) and 10 (MQTL-PH26) in rice (Table 5). In addition, wheat MQTL_2A_1, MQTL_4D_1 and MQTL_4D_1 were in the co-linear regions of MQTLs of GFeC and GZnC on chromosome 7 (rMQTL7.1), 6 (rMQTL6.3) and 7 (rMQTL7.2) in rice (Table 5). The MQTL_1B_1, MQTL_2A_1, MQTL_4A_3, MQTL_4B_2, and MQTL_4B_3 on wheat chromosomes 1B, 2A, 4A, and 4B were in located in the co-linear regions of the MQTLs of GY on chromosomes 1 (MQTL7), 5 (MQTL29), 1 (MQTL10), 8 (MQTL44), and 2 (MQTL23) in maize (Table 5). Furthermore, an OrMQTL for MQTL_3D_3 on wheat chromosome 3D was detected on chromosome 10 (MQTL107) harboring the MQTL for PH in maize (Table 5). One of the wheat MQTLs (MQTL_7B_3) on chromosome 7B was located in the syntenic region of maize MQTLs chromosome 5 (MQTL5.5) and 6 (MQTL66), respectively, as well as rice MQTLs on chromosome 6 (MQTL-YLD14) for GY (Figure 7; Table 5). The genes located at the OrMQTL regions and their annotations were reported in Supplementary Table 7. Interestingly, the well-known and proved orthologous genes including OsCYP72A18, OsFbox146, NAC22, SD37, BRD2, OsMAPK4, and ZmMPK5 were identified in rice and maize OrMQTLs.
Figure 7. The syntenic region of meta-quantitative trait loci (MQTLs) among the wheat, rice, and maize. Orthologous MQTL (OrMQTL)_10 indicates syntenic regions among identified MQTLs for grain yield (GY) control in wheat (MQTL_7B_3), rice (MQTL-YLD14), and maize (MQTL5.5 and MQTL66). The genomic position, chromosome number, and common genes among the wheat, rice, and maize are indicated.
Table 5. The OrMQTLs detected in rice and maize according to the syntenic region with MQTLs in wheat.
Discussion
Quantitative Trait Loci and Meta-Quantitative Trait Loci Distribution Over the Wheat Genome
Analysis of the genetic control of quantitative traits is a challenge in plant breeding that is due to the complexity and the difficulty of stacking numerous alleles for their control (Izquierdo et al., 2018). Meta-quantitative trait loci (MQTL) analysis made possible the consolidation of individual QTL into MQTL regions. The MQTL analysis and the comparative genomic approaches help to identify consensus QTLs and conserved genes across genomes and species. In the current MQTL analysis, we used information from over 60% of the individual QTLs to identify MQTLs that showed the statistical power of MQTL analysis for combining wheat QTLs that were in the range of 54–84% that have been reported in the previous studies of MQTL in wheat (Löffler et al., 2009; Acuña-Galindo et al., 2015; Soriano and Alvaro, 2019). In this study, we narrowed down the confidence interval (CI) of the detected MQTL regions compared to the initial QTLs used in our MQTL analysis. The efficacy of MQTL analysis in refining the CI for previously known QTLs as well as for validating their effects across different genetic backgrounds and environments is well-demonstrated (Goffinet and Gerber, 2000; Wu and Hu, 2012).
The results showed that QTLs and MQTLs had higher density in the non-telomeric and near-centromeric regions. The QTL result from the genetic segregation of sequence polymorphisms at functional elements, such as regulatory sequences upstream of genes and/or coding sequences (Flint and Mackay, 2009; Salvi and Tuberosa, 2015). Therefore, it is expected that QTL density on a genetic map is driven by gene density, polymorphism rate at functional sites in genetic regions, and the frequency of recombination. This finding was in line with the result of Martinez et al. (2016), which illustrated higher QTLs densely mapped to the near centromeres on the genetic maps.
In the quality traits category, QTLs for grain protein content (GPC), grain Fe content (GFeC), and grain Zn content (GZnC) traits were the most observed QTLs in the identified MQTLs (Supplementary Table 2). In agronomic traits, QTLs for thousand kernel weight (TKW), grain yield (GY), and number of grains per spike (NG) were the most frequent QTLs identified in the MQTLs regions. The outcome of a MQTL analysis in tetraploid wheat revealed that more than 10 loci were associated with TKW (Avni et al., 2018). Besides, in another MQTL analysis in wheat, individual QTLs for TKW, GY, and kernel number had the highest number of individual QTLs in MQTL regions. This result shows the importance of these QTLs for tested traits. The probable assumptions for a higher frequency of QTLs for agronomic traits are easy to measure and frequent data for these traits in different genetic mapping studies. On the other hand, the TKW, GY, NG, GPC, GFeC, and GZnC traits are multi-genic, highly heritable, and relatively insensitive to the environment (especially, TKW) which suggests a high likelihood of different populations carrying different suites of relevant alleles (Cooper et al., 1995; Bezant et al., 1997; Avni et al., 2018; Velu et al., 2018; Zhang et al., 2020).
MQTLs and Functional Candidate Genes
Gene annotation analysis for the MQTL regions helps clarify our understanding of their genetic architecture and refining the targets of breeding for these traits. The results of functional genomics for the identified MQTLs showed that several well-known genes consisted of psbL for electron transfer in photosystem II (PSII) (Ozawa et al., 1997), psbT encoding a PSII subunit for maintaining optimal PSII activity under adverse growth conditions (Monod et al., 1994), rpl33f or structural constituent of ribosome and translation which confers tolerance to cold stress (Rogalski et al., 2008; Moin et al., 2017), and the rps4 gene for the regulation of translational fidelity in wheat were located in the MQTl_3D_4 region. In addition, the miR166 gene was detected in the MQTL_4D_1 and MQTL_7A_4 regions. Comparative genomic approaches start with making some form of alignment of genome sequences and looking for orthologous sequences in the aligned genomes and checking to what extent those sequences are conserved. Based on these, genome and molecular evolution are inferred and this may, in turn, be put in the context of phenotypic evolution or population genetics. Analysis of the conservation and diversification of the miR166 family have shown that miR166 members play a wide and important regulatory role in seed development. More recently, a short-tandem target mimic (STTM) method was used to verify if miR166 regulates important agronomic traits in rice (Zhang et al., 2017). In a study, transgenic STTM165/166 plants showed significantly reduced seed number and sterile siliques in Arabidopsis, suggesting that miR166 plays a vital role in seed development and it might be useful evidence to improve inferior grain size in wheat (Wang et al., 2018).
Among the detected and annotated genes in the MQTLs regions in this study, the OsMED9 (TraesCS3B02G077100), which is an orthologous gene of rice, was identified in the MQTL_3B_1 region that was associated with GY and its components. A diverse array of MED genes has been identified for the regulation of GY and yield components in crop plants (Malik et al., 2016). The results indicated that the MQTL_4D_1 of our study was located in the regions of the rice orthologous BG1 (TraesCS4D02G290900), OsIDD1 (TraesCS4D02G262500), and OsPAO (TraesCS4D02G309000) genes in wheat. These genes are associated with metal ion binding, flowering time, days to heading, senescence, gravitropism, yield, and yield-related traits (Wu et al., 2008; Liu et al., 2015; Chen et al., 2016; Deng et al., 2017; Mishra et al., 2017). Overexpression of BG1 leads to larger grain size in rice (Liu et al., 2015) and manipulation of BG1 increases plant biomass, grain size, and GY in rice and Arabidopsis (Liu et al., 2015; Mishra et al., 2017). OsIDD1 could rescue the never-flowering phenotype of rid1 by a transition from vegetative to reproductive growth in rice (Deng et al., 2017), and the PAOs protein-encoding genes (Chen et al., 2016) regulate cellular polyamine levels which are critical for embryogenesis (Bertoldi et al., 2004; De-la-Pena et al., 2008), germination (Bethke et al., 2004; Liszkay et al., 2004), root growth (Cona et al., 2005), flowering and senescence (Kakkar and Sawhney, 2002), and mineral deficiency (Moschou et al., 2008, 2009). The orthologous Ghd7 (TraesCS5A02G541200) that was among 427 genes located in the MQTL_5A_4 region in our study involves photoperiodism, flowering, days to heading, plant height (PH), and yield traits. Enhanced expression of Ghd7 under long-day conditions delays heading and increases PH and panicle size in rice. The Ghd7 gene plays crucial role in increasing the productivity and adaptability of rice globally (Xue et al., 2008). The uncovered rice orthologous D27 (TraesCS7B02G319100) and BRD2 (TraesCS7B02G484200) genes in wheat in chromosome 7B region are known to control tillering, heading date, PH, and yield traits in rice (Hong et al., 2005; Lin et al., 2009; Liu et al., 2016).
The identified rice orthologous D18 (TraesCS2B02G570800), OsRPK1 (TraesCS5D02G034200), DRUS1, and DRUS2 (TraesCS4A02G133800) genes within the MQTL_2B_2, MQTL_4A_2, and MQTL_5D_2 regions of our MQTLs were associated with PH in wheat (Itoh et al., 2002; Zou et al., 2014; Pu et al., 2017). Manipulating these genes lead to varieties with dwarf and semi-dwarf phenotype and such phenotypes possess short and strong stalks, exhibit less lodging, and a greater proportion of assimilation partitioned into the grain, resulting in further yield increases (Hedden, 2003).
The identified rice orthologous OsRLCK189 (TraesCS1A02G317300), OsGH3-4 (TraesCS1A02G320200), FT-L (TraesCS2B02G511400), OsRLCK57 (TraesCS3B02G600300), OsHK1 (TraesCS2B02G495500), AIM1 (TraesCS3D02G077200), and WOX11 (TraesCS2A02G100700) genes in wheat was located in MQTL_1A_4, MQTL_2A_2, MQTL_2B_1, MQTL_3B_5, and MQTl_3D_3 regions show a role in root growth and inflorescence and floral development (Richmond and Bleecker, 1999; Izawa et al., 2002; Vij et al., 2008; Ogiso-Tanaka et al., 2013; Zhao et al., 2015, 2020; Cheng et al., 2016; Xu et al., 2017; Lehner et al., 2018; Kong et al., 2019). Inflorescence development in cereals directly affects grain number and size which are key determinants of yield (Yamburenko et al., 2017).
The rice orthologous OsMTP12 (TraesCS2A02G141400), OsMTP9 (TraesCS3B02G040900), and OsMTP1 (TraesCS4D02G323700) belonging to the metal tolerance protein (MTP) gene family were identified in the MQTL_2A_1, MQTL_3B_1, and MQTL_4D_1 interval of the wheat genome in the present study. The MTP gene family plays a critical role in metal transport, mainly in Zn, Mn, Fe, Cd, Co, and Ni, metal homeostasis, and tolerance (Gustin et al., 2011; Zhang and Liu, 2017; Ram et al., 2019). The increased expression of MTP genes during seed development and their potential role in metal homeostasis during the seed filling stage had been documented (Ram et al., 2019). The MTP1 and MTP12 genes share a characteristic histidine-rich loop toward the c-terminal, which is known to have a role in Zn-binding (Ram et al., 2019).
The MQTL_1A_4, MQTL_3B_5, MQTL_4B_1, MQTL_4B_2, MQTL_4B_3, MQTL_4B_4, MQTL_7A_1, and MQTL_7A_4 interval regions harbored rice orthologous VAL1 (TraesCS4B02G314600), PROG1 (TraesCS4B02G354000), D14 (TraesCS4B02G258200), LPL2 (TraesCS4B02G308000), OsINV3 (TraesCS7A02G009100), and OsRLCK218 (TraesCS7A02G385300) genes with a role in diverse development and growth traits including leaf development, PH, shoot branching, panicle and tiller number, grain number, grain and spikelet size, and GY (Vij et al., 2008; Zhou et al., 2016; Morey et al., 2018; Wu et al., 2018; Yao et al., 2018; Zhang et al., 2018; Deng et al., 2020). These candidate genes (CGs) at the detected MQTL regions can potentially have the same function as their orthologous varieties in rice and therefore regulate various developmental and growth-related traits in wheat. Identification of these confines stable chromosomal regions and CGs that influence economically important traits in wheat can help to expedite wheat improvement in future breeding programs.
The MQTLs detected in this study will help identify CGs in these regions responsible for desirable traits and generate allele-specific markers through allele mining (Leung et al., 2015; Ogawa et al., 2018) for marker-assisted selection (MAS) application in pre-breeding population. Allele mining is a promising approach to dissect naturally occurring allelic variation at QTLs/CGs controlling desirable traits (e.g., yield, quality, and micronutrient content) which has potential applications in crop improvement programs (Kumar et al., 2010; Gokidi et al., 2017; Kumari et al., 2018). The data raised from this MQTL study help to refine genomic targets for validation and subsequent development of haplotypic markers in breeding programs. Information of the identified MQTLs can also be used for genetic transformation or allele screening in germplasm collections (ecoTilling) for finding new alleles capable of improving yield, quality, and micronutrient traits (Izquierdo et al., 2018; Belzile et al., 2020). Additionally, a promising application of the variation within these MQTL regions might be their introduction as fixed effects in genomic selection (GS) models to increase the accuracy of the prediction models in their use in wheat breeding programs (Spindel et al., 2016; Izquierdo et al., 2018).
MQTL and Traits Analysis
The QTLs projection on a consensus map allows for the inspection of co-location across traits and categories, which is especially relevant for complex traits (Delfino et al., 2019). The association between trait classes by analyzing the co-localization frequency of individual trait-QTLs demonstrated that TKW with 55% and 63% co-localization frequencies was frequently associated with GY and GPC. In addition, GY and GFeC had the highest co-localization frequency with GFeC (52%) and GZnC (66%), respectively. The results of MQTL analysis of our study confirmed correlations identified for QTLs of TKW and GY (An-Ming et al., 2011; Azadi et al., 2014; Mahdi-Nezhad et al., 2019), TKW and GPC (Wang L. I. et al., 2012; Goel et al., 2019), and GFeC and GZn (Roshanzamir et al., 2013; Pu et al., 2014; Tiwari et al., 2016; Liu et al., 2019) in other studies. Co-localization of QTLs for correlated traits has been identified by Wang et al. (2018). Co-localization of QTLs could be due to the pleiotropy or the presence of different linked genes in the same regions that can partly explain the correlation that exists between traits (Bhatta et al., 2018). The tightly linked genomic regions or pleiotropy can partly explain the correlation that exists between traits (Crespo-Herrera et al., 2016). The result of co-localization frequency of target traits (GY, GPC, GFeC, and GZnC) with the detected MQTLs in our study showed interacting MQTLs which affects the association of traits at the genomic level. Acuña-Galindo et al. (2015) demonstrated that TKW was most frequently associated with GY in MQTL regions, with 57% co-localization. The co-localization of TKW and GPC QTLs has been observed in the genetic analysis of wheat (Goel et al., 2019). The positive and highly significant associations of GPC and TKW traits have been reported in different wheat populations in other studies (Peleg et al., 2009; Badakhshan et al., 2013; Krishnappa et al., 2017). According to Liu et al. (2019), mineral nutrient concentrations and yield components have shown common QTLs which is in line with co-localized QTLs for grain Zn, Fe QTLs in the tetraploid and hexaploid wheat in other studies (Peleg et al., 2009; Xu et al., 2012; Crespo-Herrera et al., 2016; Krishnappa et al., 2017; Velu et al., 2017). In a recent QTL analysis study on rice, grain Fe and grain yield QTLs were found to be co-localized (Dixit et al., 2019). Crespo-Herrera et al. (2016) suggested the possibility of simultaneous breeding for GFeC and GZnC traits due to their co-localized QTLs. The current co-localization analysis within MQTLs regions indicates the possibility of simultaneous breeding of micronutrient content, grain quality, and GY by pyramiding the QTL regions through marker-assisted selection (MAS). Using -MAS, the specific regions can be transferred to the elite wheat genotypes to simultaneously increase the contents of various traits (Saini et al., 2020). Therefore, an attempt to pyramid QTLs responsible for GY, TKW, GPC, GZnC, and GFeC may accelerate progress in wheat variety development. The QTL pyramiding strategy has been used for simultaneous improvement of traits through MAS in wheat (Wang P. et al., 2012; Feng et al., 2018; Gautam et al., 2020; Muthu et al., 2020). A clear understanding of the co-location of QTLs and their effect on target traits, such as grain Fe and Zn concentrations and yield is very important for using the major effect of QTLs in marker-assisted breeding (Swamy et al., 2018). The results of MQTL analysis in our study showed that the location of 18 MQTLs was in agreement with the position of the QTLs in the meta-QTL studies by Zhang A. et al. (2010), Acuña-Galindo et al. (2015), and Kumar et al. (2020). However, our work adds to this body of genomic mapping research by identifying new MQTL specific for GY, grain quality, and micronutrient content.
Comparative Genomic Analysis and Orthologous MQTL
The genome-wide association studies (GWAS) approach can help in gene discovery and in the analysis of the genetic basis of complex traits for the improvement of wheat. Comparative analysis of the wheat GWAS with the identified current MQTLs suggested the co-localization of 12 MQTLs (MQTL_1B_2, MQTL_2B_1, MQTL_3B_1, MQTL_4A_2, MQTL_4A_3, MQTL_4A_4, MQTL_4B_3, MQTL_4D_1, MQTL_5A_4, MQTL_7A_1, MQTL_7A_3, and MQTL_7B_2) and significant genomic regions of wheat traits that have been identified in the wheat GWAS studies (Jamil et al., 2019; Li et al., 2019; Mathew et al., 2019; Hu et al., 2020; Li X. et al., 2020). Some of the QTL hotspots suggested common genetic markers for wheat traits for further use in marker-assisted breeding that increase genetic gain in breeding programs. The identified OrMQTLs in this study can facilitate the detection of the underlying regulatory genes with evolutionary history and conservative function. The current MQTL analysis defines a genome-wide landscape on the most stable genetic markers and CGs related to micronutrient content, yield, and yield-related traits as the most economically important traits in wheat. The straight-up comparative genomics in our study indicated high synteny between wheat, rice, and maize QTLs especially for GY suggesting possible corroborating evidence of acting the same way in different species that was in line with the results of other MQTL studies (Ahn et al., 1993; Moore et al., 1995; Gale and Devos, 1998; Feuillet and Keller, 1999; Minx et al., 2005). Despite the high interest in the identification of genes involved in GY and yield-related traits in maize and wheat as two economically important crops, the responsible genes have largely remained unknown due to their complex genomes. Given a close evolutionary relation among grass genomes (Gaut, 2002), a synteny analysis of wheat, maize, and rice through the identification of OrMQTLs enabled us to broaden our genetic information and leads to uncover the possible function of unknown CGs among these species. Here we selected the most promising wheat MQTLs to explore their conserved syntenic regions reported in similar MQTLs studies on the same traits in rice and maize to identify OrMQTLs (Table 5). For the wheat MQTL_3B_4 of our study, there is an MQTL in the syntenic region in rice (Khahani et al., 2020) controlling OrMQTL_2 containing a rice OsCYP72A18 orthologous gene (TraesCS3B02G609400 and TraesCS3B02G609600) for grGY (Swamy and Sarla, 2011). Moreover, the wheat MQTL_4D_1 and MQTL_7B_2 of our study were located in the syntenic regions of the rice genome on chromosomes 3 (OrMQTL_3) and 6 (OrMQTL_1), respectively. These two OrMQTLs encompass OsFbox146 (TraesCS4D02G288900) and OsFbox297 (TraesCS7B02G405900) orthologous genes for genetic control of GY in rice (Swamy and Sarla, 2011; Lei et al., 2018). Furthermore, the syntenic region of wheat MQTL_4B_1 possessed the two orthologous rice genes, NAC22 (TraesCS4B02G328600) and CYP96B4/SD37 (TraesCS4B02G342400) located on chromosome 3 (OrMQTL_11) which control plant growth (Tamiru et al., 2015; Hong et al., 2016). In the syntenic regions of the MQTL_2A_1 and MQTL_4D_1 of our study, there was orthologous OsZIP1/OsZIP8 (TraesCS2A02G143400) and OsFerroportin (TraesCS4D02G323100) genes on chromosomes 7 (OrMQTL_14) and 6 (OrMQTL_15) of rice which are related to Zn and Fe transport/homeostasis (Morrissey and Guerinot, 2009; Bashir et al., 2011; Alagarasan et al., 2017) (Table 5).
More intriguingly, the syntenic regions of the wheat MQTL_7B_3 on both chromosome 6 of rice (OrMQTL_10) and chromosome 5 of maize, harbored the rice (OsMAPK4, Os06g0699400) and maize (ZmMPK5, Zm00001d014658) and orthologous gene in wheat (TraesCS7B02G322900). The mitogen-activated protein kinase (MAPK) cascades play important roles in regulating plant growth (PH), development, and stress responses (Zhu et al., 2020). The maize ZmMPK5 is induced by various stimuli, involved in defense signaling pathways in various abiotic/biotic stress, confers tolerance to stresses (Zhang A. et al., 2010; Zhang et al., 2014), and subsequently leads to enhancement of plant growth and yield. Due to the key role of OsMAPK4 in plant growth, grain development (Liu et al., 2018; Chen et al., 2021), and subsequently in GY, the results of the current syntenic analysis suggest the same function for TraesCS7B02G322900 and Zm00001d014658 genes located in OrMQTL_10 interval (Figure 7; Table 5; Supplementary Table 7).
Conclusion
The results of this study introduced several novel meta-quantitative trait loci (MQTLs) for improving wheat in multi-purpose breeding programs by identifying key genomic regions associated with agronomic performance, grain quality traits, and micronutrients content. The results of our MQTL analysis have significantly increased the power and precision of our ability to map wheat traits and will provide greater resolution for future fine mapping and marker development. Importantly, our data identify co-localization between grain yield (GY) QTLs with grain Zn content (GZnC), grain Fe content (GFeC), and grain protein content (GPC) QTLs that suggest an opportunity for simultaneous breeding for these traits. This study also provides an example for the value of comparative analysis between evolutionarily close cereal species for the identification of genomic regions and candidate genes (CGS) controlling quantitative traits. Our finding shows the utility of MQTL analysis for refining the location of genomic regions associated with a variety of traits and helps understand how their relative map positions can be exploited for crop improvement. Overall, these findings can lead to both increased selection efficiency and accuracy for breeding by providing the basis for marker development in a marker-assisted selection (MAS) program and for identifying a novel source of variation through allele mining efforts in genetic resource collections. Lastly, these refined MQTLs provide the basis for further focus on the genetic mechanisms controlling micronutrients, GY, and quality traits.
Data Availability Statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.
Author Contributions
NS conducted bioinformatics analysis and wrote the draft of the article. BH conceived and designed the project, the bioinformatics analysis, and complemented the writing of the article. AT helped in bioinformatics analysis. CR helped with reviewing and providing critical advice on the article. All authors have read and approved the manuscript.
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.
Acknowledgments
The authors gratefully acknowledge the supports from Shiraz University.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2021.709817/full#supplementary-material
Supplementary Figure 1. The location of detected MQTLs with 95% confidence interval associated with quantitative traits in wheat chromosomes. The Lines on the left side of the linkage groups indicate the confidence interval (CI) of QTLs. The colored boxes on each linkage groups represent MQTLs region and the colors inside the vertical lines on the left of illustrate the best model of MQTLs. The molecular markers and their genetic distance (cM) over linkage groups are shown on the right side.
Supplementary Table 1. The information of wheat population used for meta- quantitative trait loci (MQTL) analysis.
Supplementary Table 2. The initial and projected QTL data on the genetic consensus map and identified MQTLs for the studied traits.
Supplementary Table 3. The information of QTLs in the detected MQTLs.
Supplementary Table 4. The list of uncovered candidate genes (CGS) and all annotated genes anchoring at each MQTL interval.
Supplementary Table 5. The list of rice orthologous genes at each MQTL interval.
Supplementary Table 6. The information of the chi-square test for co-localization frequency analysis of the QTLs of the tested traits in the detected MQTLs regions.
Supplementary Table 7. The list of orthologous genes located at the syntenic regions of detected ortho-MQTLs between wheat, maize, and rice.
Abbreviations
AIC, akaike information criterion; AICc, corrected AIC; AWE, the average weight of evidence; BIC, Bayesian information criterion; CI, confidence interval; GCs, candidate genes; GO, gene ontology; GWAS, genome wide association studies; LOD, the logarithm of odds; MQTL, meta- quantitative trait loci; MTP, metal tolerance proteins; OrMQTL, orthologous MQTL; QTL, quantitative trait loci; SNP, single-nucleotide polymorphism.
References
Abdollahi-Sisi, N., Mohammadi, S. A., and Razeghi, J. (2018). Mapping QTLs for grain yield components in bread wheat under well-watered and rain-fed conditions. J. Bio. Env. Sci.13, 306–314.
Acuña-Galindo, M. A., Mason, R. E., Subramanian, N. K., and Hays, D. B. (2015). Meta-analysis of wheat QTL regions associated with adaptation to drought and heat stress. Crop. Sci. 55, 477–492. doi: 10.2135/cropsci2013.11.0793
Ahn, S., Anderson, J. A., Sorrells, M. E., and Tanksley, S. D. (1993). Homoeologous relationships of rice, wheat and maize chromosomes. Mol. Gen. Genet. 241, 483–490. doi: 10.1007/BF00279889
Alagarasan, G., Dubey, M., Aswathy, K. S., and Chandel, G. (2017). Genome wide identification of orthologous ZIP genes associated with zinc and Iron translocation in Setaria italica. Front. Plant. Sci. 8:775. doi: 10.3389/fpls.2017.00775
An-Ming, D. I., Jun, L. I., Fa, C. U., Chun-Hua, Z. H., Hang-Yun, M. A., and Hong-Gang, W. A. (2011). Mapping QTLs for yield related traits using two associated RIL populations of wheat. Zuo Wu Xue Bao 37, 1511–1524. doi: 10.3724/SP.J.1006.2011.01511
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. Bioinformatics 20, 2324–2326. doi: 10.1093/bioinformatics/bth230
Avni, R., Oren, L., Shabtay, G., Assili, S., Pozniak, C., Hale, I., et al. (2018). Genome based meta-QTL analysis of grain weight in tetraploid wheat identifies rare alleles of GRF4 associated with larger grains. Genes 9:636. doi: 10.3390/genes9120636
Azadi, A., Mardi, M., Hervan, E. M., Mohammadi, S. A., Moradi, F., Tabatabaee, M. T., et al. (2014). QTL mapping of yield and yield components under normal and salt-stress conditions in bread wheat (Triticum aestivum L.). Plant Mol. Biol. Rep. 33, 102–120. doi: 10.1007/s11105-014-0726-0
Badakhshan, H., Moradi, N., Mohammadzadeh, H., and Zakeri, M. R. (2013). Genetic variability analysis of grains Fe, Zn and beta-carotene concentration of prevalent wheat varieties in Iran. Int. J. Agri. Crop Sci. 6:57.
Badji, A., Otim, M., Machida, L., Odong, T., Kwemoi, D. B., Okii, D., et al. (2018). Maize combined insect resistance genomic regions and their co-localization with Cell Wall constituents revealed by tissue-specific QTL meta-analyses. Front. Plant Sci. 9:895. doi: 10.3389/fpls.2018.00895
Bashir, K., Ishimaru, Y., Shimo, H., Nagasaka, S., Fujimoto, M., Takanashi, H., et al. (2011). The rice mitochondrial iron transporter is essential for plant growth. Nat. Commun. 2, 1–7. doi: 10.1038/ncomms1326
Belzile, F., Abed, A., and Torkamaneh, D. (2020). Time for a paradigm shift in the use of plant genetic resources. Genome 63, 189–194. doi: 10.1139/gen-2019-0141
Bertoldi, D., Tassoni, A., Martinelli, L., and Bagni, N. (2004). Polyamines and somatic embryogenesis in two Vitis vinifera cultivars. Physiol. Plant. 120, 657–666. doi: 10.1111/j.0031-9317.2004.0282.x
Bethke, P. C., Gubler, F., Jacobsen, J. V., and Jones, R. L. (2004). Dormancy of Arabidopsis seeds and barley grains can be broken by nitric oxide. Planta 219, 847–855. doi: 10.1007/s00425-004-1282-x
Bezant, H., Laurie, D. A., Pratchett, N., Chojecki, J., and Kearsey, M. J. (1997). Mapping of QTL controlling NIR predicted hot water extract and grain nitrogen content in a spring barley cross using marker-regression. Plant Breed. 116, 141–145. doi: 10.1111/j.1439-0523.1997.tb02168.x
Bhatta, M., Baenziger, P., Waters, B., Poudel, R., Belamkar, V., Poland, J., et al. (2018). Genome-wide association study reveals novel genomic regions associated with 10 grain minerals in synthetic hexaploid wheat. Int. J. Mol. Sci. 19:3237. doi: 10.3390/ijms19103237
Bhusal, N., Sarial, A. K., Sharma, P., and Sareen, S. (2017). Mapping QTLs for grain yield components in wheat under heat stress. PLoS ONE 12:e0189594. doi: 10.1371/journal.pone.0189594
Black, R. E., Victora, C. G., Walker, S. P., Bhutta, Z. A., Christian, P., de Onis, M., et al. (2013). Maternal and child undernutrition and overweight in low-income and middle-income countries. Lancet 382, 427–451. doi: 10.1016/S0140-6736(13)60937-X
Braun, H. J., Atlin, G., and Payne, T. (2010). “Multi-location testing as a tool to identify plant response to global climate change,” in Climate change and crop production, ed. M. P. Reynolds (CABI Press, Oxford), 115–138. doi: 10.1079/9781845936334.0115
Chen, B. X., Li, W. Y., Gao, Y. T., Chen, Z. J., Zhang, W. N., Liu, Q. J., et al. (2016). Involvement of polyamine oxidase-produced hydrogen peroxide during coleorhiza-limited germination of rice seeds. Front. Plant. Sci. 7:1219. doi: 10.3389/fpls.2016.01219
Chen, J., Wang, L., and Yuan, M. (2021). Update on the roles of rice MAPK cascades. Int. J. Mol. Sci. 22:1679. doi: 10.3390/ijms22041679
Cheng, S., Zhou, D. X., and Zhao, Y. (2016). WUSCHEL-related homeobox gene WOX11 increases rice drought resistance by controlling root hair formation and root system development. Plant Signal Behav. 11:e1130198. doi: 10.1080/15592324.2015.1130198
Cona, A., Moreno, S., Cenci, F., Federico, R., and Angelini, R. (2005). Cellular re-distribution of flavin-containing polyamine oxidase in differentiating root and mesocotyl of Zea mays L. seedlings. Planta 221, 265–276. doi: 10.1007/s00425-004-1435-y
Cooper, M., Woodruff, D. R., Eisemann, R. L., Brennan, P. S., and DeLacy, I. H. (1995). A selection strategy to accommodate genotype-by-environment interaction for grain yield of wheat: managed-environments for selection among genotypes. Theor. Appl. Genet. 90, 492–502. doi: 10.1007/BF00221995
Crespo-Herrera, L. A., Velu, G., and Singh, R. P. (2016). Quantitative trait loci mapping reveals pleiotropic effect for grain iron and zinc concentrations in wheat. Ann. Appl. Biol. 169, 27–35. doi: 10.1111/aab.12276
Danan, S., Veyrieras, J. B., and Lefebvre, V. (2011). Construction of a potato consensus map and QTL meta-analysis offer new insights into the genetic architecture of late blight resistance and plant maturity traits. BMC Plant Biol. 11:16. doi: 10.1186/1471-2229-11-16
De-la-Pena, C., Galaz-Avalos, R. M., and Loyola-Vargas, V. M. (2008). Possible role of light and polyamines in the onset of somatic embryogenesis of Coffea canephora. Mol. Biotechnol. 39, 215–224. doi: 10.1007/s12033-008-9037-8
Delfino, P., Zenoni, S., Imanifard, Z., Tornielli, G. B., and Bellin, D. (2019). Selection of candidate genes controlling veraison time in grapevine through integration of meta-QTL and transcriptomic data. BMC Genomics 20, 1–19. doi: 10.1186/s12864-019-6124-0
Deng, L., Li, L., Zhang, S., Shen, J., Li, S., Hu, S., et al. (2017). Suppressor of rid1 (SID1) shares common targets with RID1 on florigen genes to initiate floral transition in rice. PLoS Genet. 13:e1006642. doi: 10.1371/journal.pgen.1006642
Deng, X., Han, X., Yu, S., Liu, Z., Guo, D., He, Y., et al. (2020). OsINV3 and its homolog, OsINV2, control grain size in rice. Int. J. Mol. Sci. 21:2199. doi: 10.3390/ijms21062199
Dixit, S., Singh, U. M., Abbai, R., Ram, T., Singh, V. K., Paul, A., et al. (2019). Identification of genomic region (s) responsible for high iron and zinc content in rice. Sci. Rep. 9:8136. doi: 10.1038/s41598-019-43888-y
El-Feki, W. M., Byrne, P. F., Reid, S. D., and Haley, S. D. (2018). Mapping quantitative trait loci for agronomic traits in winter wheat under different soil moisture levels. Agronomy 8:133. doi: 10.3390/agronomy8080133
Feng, B., Chen, K., Cui, Y., Wu, Z., Zheng, T., Zhu, Y., et al. (2018). Genetic dissection and simultaneous improvement of drought and low nitrogen tolerances by designed QTL pyramiding in rice. Front. Plant. Sci. 9:306. doi: 10.3389/fpls.2018.00306
Fernando, N., Panozzo, J., Tausz, M., Norton, R., Fitzgerald, G., and Seneweera, S. (2012). Rising atmospheric CO2 concentration affects mineral nutrient and protein concentration of wheat grain. Food. Chem. 133, 1307–1311. doi: 10.1016/j.foodchem.2012.01.105
Feuillet, C., and Keller, B. (1999). High gene density is conserved at syntenic loci of small and large grass genomes. Proc. Natl. Acad. Sci. USA. 96, 8265–8270. doi: 10.1073/pnas.96.14.8265
Flatt, T., and Heyland, A. (2011). Mechanisms of Life History Evolution: the Genetics and Physiology of Life History Traits and Trade-Offs. Oxford: Oxford University Press. doi: 10.1093/acprof:oso/9780199568765.001.0001
Flint, J., and Mackay, T. F. (2009). Genetic architecture of quantitative traits in mice, flies, and humans. Genome Res. 19, 23–733. doi: 10.1101/gr.086660.108
Gahlaut, V., Jaiswal, V., Tyagi, B. S., Singh, G., Sareen, S., Balyan, H. S., et al. (2017). QTL mapping for nine drought-responsive agronomic traits in bread wheat under irrigated and rain-fed environments. PLoS ONE 12:e0182857. doi: 10.1371/journal.pone.0182857
Gale, M. D., and Devos, K. M. (1998). Comparative genetics in the grasses. Proc. Natl. Acad. Sci. USA. 95, 1971–1974. doi: 10.1073/pnas.95.5.1971
Gaut, B. S. (2002). Evolutionary dynamics of grass genomes. New Phytol. 154, 15–28. doi: 10.1046/j.1469-8137.2002.00352.x
Gautam, T., Dhillon, G. S., Saripalli, G., Singh, V. P., Prasad, P., Kaur, S., et al. (2020). Marker-assisted pyramiding of genes/QTL for grain quality and rust resistance in wheat (Triticum aestivum L.). Mol. Breed. 40, 1–14. doi: 10.1007/s11032-020-01125-9
Giancaspro, A., Giove, S. L., Zacheo, S. A., Blanco, A., and Gadaleta, A. (2019). Genetic variation for protein content and yield-related traits in a durum population derived from an inter-specific cross between hexaploid and tetraploid wheat cultivars. Front. Plant Sci. 10:1509. doi: 10.3389/fpls.2019.01509
Goel, S., Singh, K., Singh, B., Grewal, S., Dwivedi, N., Alqarawi, A. A., et al. (2019). Analysis of genetic control and QTL mapping of essential wheat grain quality traits in a recombinant inbred population. PLoS ONE 14:e0200669. doi: 10.1371/journal.pone.0200669
Goffinet, B., and Gerber, S. (2000). Quantitative trait loci: a meta-analysis. Genetics 155, 463–473. doi: 10.1093/genetics/155.1.463
Gokidi, Y., Bhanu, A. N., Chandra, K., Singh, M. N., and Hemantaranjan, A. (2017). Allele mining–an approach to discover allelic variation in crops. J. Plant Sci. Res. 33, 167–180.
Guo, Y., Du, Z., Chen, J., and Zhang, Z. (2017). QTL mapping of wheat plant architectural characteristics and their genetic relationship with seven QTLs conferring resistance to sheath blight. PLoS ONE 12:e0174939. doi: 10.1371/journal.pone.0174939
Gustin, J. L., Zanis, M. J., and Salt, D. E. (2011). Structure and evolution of the plant cation diffusion facilitator family of ion transporters. BMC Evol. Biol. 11:76. doi: 10.1186/1471-2148-11-76
Hanocq, E., Laperche, A., Jaminon, O., Lainé, A. L., and Le Gouis, J. (2007). Most significant genome regions involved in the control of earliness traits in bread wheat, as revealed by QTL meta-analysis. Theor. Appl. Genet. 114, 569–584. doi: 10.1007/s00122-006-0459-z
Hedden, P. (2003). The genes of the Green Revolution. Trends Genet. 19, 5–9. doi: 10.1016/S0168-9525(02)00009-4
Hong, Y., Zhang, H., Huang, L., Li, D., and Song, F. (2016). Overexpression of a stress-responsive NAC transcription factor gene ONAC022 improves drought and salt tolerance in rice. Front. Plant Sci. 7:4. doi: 10.3389/fpls.2016.00004
Hong, Z., Ueguchi-Tanaka, M., Fujioka, S., Takatsuto, S., Yoshida, S., Hasegawa, Y., et al. (2005). The rice brassinosteroid-deficient dwarf2 mutant, defective in the rice homolog of Arabidopsis DIMINUTO/DWARF1, is rescued by the endogenously accumulated alternative bioactive brassinosteroid, dolichosterone. Plant Cell 17, 2243–2254. doi: 10.1105/tpc.105.030973
Hu, P., Zheng, Q., Luo, Q., Teng, W., Li, H., Li, B., et al. (2020). Genome-wide association study of yield and related traits in common wheat under salt-stress conditions. BMC Plant Biol. 21:27. doi: 10.1186/s12870-020-02799-1
Itoh, H., Ueguchi-Tanaka, M., Sakamoto, T., Kayano, T., Tanaka, H., Ashikari, M., et al. (2002). Modification of rice plant height by suppressing the height-controlling gene, D18, in rice. Breed. Sci. 52, 215–218. doi: 10.1270/jsbbs.52.215
Izawa, T., Oikawa, T., Sugiyama, N., Tanisaka, T., Yano, M., and Shimamoto, K. (2002). Phytochrome mediates the external light signal to repress FT orthologs in photoperiodic flowering of rice. Genes Dev. 16, 2006–2020. doi: 10.1101/gad.999202
Izquierdo, P., Astudillo, C., Blair, M. W., Iqbal, A. M., Raatz, B., and Cichy, K. A. (2018). Meta-QTL analysis of seed iron and zinc concentration and content in common bean (Phaseolus vulgaris L.). Theor. Appl. Genet. 131, 1645–1658. doi: 10.1007/s00122-018-3104-8
Jamil, M., Ali, A., Gul, A., Ghafoor, A., Napar, A. A., Ibrahim, A. M., et al. (2019). Genome-wide association studies of seven agronomic traits under two sowing conditions in bread wheat. BMC Plant Biol. 19:149. doi: 10.1186/s12870-019-1754-6
Kakkar, R. K., and Sawhney, V. K. (2002). Polyamine research in plants—a changing perspective. Physiol. Plant. 116, 281–292. doi: 10.1034/j.1399-3054.2002.1160302.x
Khahani, B., Tavakol, E., Shariati, V., and Fornara, F. (2020). Genome wide screening and comparative genome analysis for meta-QTLs, ortho-MQTLs and candidate genes controlling yield and yield-related traits in rice. BMC Genomics 21, 1–24. doi: 10.1186/s12864-020-6702-1
Khush, G. S., Lee, S., Cho, J. I., and Jeon, J. S. (2012). Biofortification of crops for reducing malnutrition. Plant Biotechnol. Rep. 6, 195–202. doi: 10.1007/s11816-012-0216-5
Kolde, R. (2013). pheatmap: Pretty Heatmaps. R package version 1.0.12. Available online at: http://CRAN.R-project.org/package=pheatmap (accessed December 26, 2018).
Kong, W., Zhang, Y., Deng, X., Zhang, C., and Li, Y. (2019). Comparative genomic and transcriptomic analyses suggests the evolutionary dynamic of gh3 genes in gramineae crops. Front. Plant Sci. 10:1297. doi: 10.3389/fpls.2019.01297
Krishnappa, G., Singh, A. M., Chaudhary, S., Ahlawat, A. K., Singh, S. K., Shukla, R. B., et al. (2017). Molecular mapping of the grain iron and zinc concentration, protein content and thousand kernel weight in wheat (Triticum aestivum L.). PLoS ONE 12:e0174972. doi: 10.1371/journal.pone.0174972
Kumar, A., Saripalli, G., Jan, I., Kumar, K., Sharma, P. K., Balyan, H. S., et al. (2020). Meta-QTL analysis and identification of candidate genes for drought tolerance in bread wheat (Triticum aestivum L.). Physiol. Mol. Biol. Plants 26, 1713–1725. doi: 10.1007/s12298-020-00847-6
Kumar, G. R., Sakthivel, K., Sundaram, R. M., Neeraja, C. N., Balachandran, S. M., Rani, N. S., et al. (2010). Allele mining in crops: prospects and potentials. Biotechnol. Adv. 28, 451–461. doi: 10.1016/j.biotechadv.2010.02.007
Kumar, S., Palve, A., Joshi, C., and Srivastava, R. K. (2019). Crop biofortification for iron (Fe), zinc (Zn) and vitamin A with transgenic approaches. Heliyon 5:e01914. doi: 10.1016/j.heliyon.2019.e01914
Kumari, R., Kumar, P., Sharma, V. K., and Kumar, H. (2018). Allele mining for crop improvement. Int. J. Pure Appl. Biosci. 6, 1456–1465. doi: 10.18782/2320-7051.6073
Lehner, K. R., Taylor, I., McCaskey, E. N., Jain, R., Ronald, P. C., Goldman, D. I., et al. (2018). A histidine kinase gene is required for large radius root tip circumnutation and surface exploration in rice. bioRxiv. 437012. doi: 10.1101/437012
Lei, L., Zheng, H. L., Wang, J. G., Liu, H. L., Sun, J., Zhao, H. W., et al. (2018). Genetic dissection of rice (Oryza sativa L.) tiller, plant height, and grain yield based on QTL mapping and metaanalysis. Euphytica 214, 1–7. doi: 10.1007/s10681-018-2187-2
Leung, H., Raghavan, C., Zhou, B., Oliva, R., Choi, I. R., Lacorte, V., et al. (2015). Allele mining and enhanced genetic recombination for rice breeding. Rice 8:34. doi: 10.1186/s12284-015-0069-y
Li, F., Wen, W., Liu, J., Zhang, Y., Cao, S., He, Z., et al. (2019). Genetic architecture of grain yield in bread wheat based on genome-wide association studies. BMC Plant Biol. 19:168. doi: 10.1186/s12870-019-1781-3
Li, Q., Pan, Z., Gao, Y., Li, T., Liang, J., Zhang, Z., et al. (2020). Quantitative trait locus (QTLs) mapping for quality traits of wheat based on high density genetic map combined with bulked segregant analysis RNA-seq (BSR-Seq) indicates that the basic 7S globulin gene is related to falling number. Front. Plant Sci. 11:1918. doi: 10.3389/fpls.2020.600788
Li, W. T., Liu, C., Liu, Y. X., Pu, Z. E., Dai, S. F., Wang, J. R., et al. (2013). Meta-analysis of QTL associated with tolerance to abiotic stresses in barley. Euphytica 189, 31–49. doi: 10.1007/s10681-012-0683-3
Li, X., Xu, X., Liu, W., Li, X., Yang, X., Ru, Z., et al. (2020). Dissection of superior alleles for yield-related traits and their distribution in important cultivars of wheat by association mapping. Front. Plant Sci. 11:175. doi: 10.3389/fpls.2020.00175
Liang, Y., Zhang, K., Zhao, L., Liu, B., and Meng, Q., Tian, J., et al. (2010). Identification of chromosome regions conferring dry matter accumulation and photosynthesis in wheat (Triticum aestivum L.). Euphytica 171, 145–156. doi: 10.1007/s10681-009-0024-3
Lin, H., Wang, R., Qian, Q., Yan, M., Meng, X., Fu, Z., et al. (2009). DWARF27, an iron-containing protein required for the biosynthesis of strigolactones, regulates rice tiller bud outgrowth. Plant Cell 21, 1512–1525. doi: 10.1105/tpc.109.065987
Liszkay, A., van der Zalm, E., and Schopfer, P. (2004). Production of reactive oxygen intermediates (O2-, H2O2, and OH) by maize roots and their role in wall loosening and elongation growth. Plant Physiol. 136, 3114–3123. doi: 10.1104/pp.104.044784
Liu, J., Wu, B., Singh, R. P., and Velu, G. (2019). QTL mapping for micronutrients concentration and yield component traits in a hexaploid wheat mapping population. J. Cereal Sci. 88, 57–64. doi: 10.1016/j.jcs.2019.05.008
Liu, L., Tong, H., Xiao, Y., Che, R., Xu, F., Hu, B., et al. (2015). Activation of Big Grain1 significantly improves grain size by regulating auxin transport in rice. Proc. Natl. Acad. Sci. USA. 112, 11102–11107. doi: 10.1073/pnas.1512748112
Liu, X., Feng, Z. M., Zhou, C. L., Ren, Y. K., Mou, C. L., et al. (2016). Brassinosteroid (BR) biosynthetic gene lhdd10 controls late heading and plant height in rice (Oryza sativa L.). Plant Cell Rep. 35, 357–368. doi: 10.1007/s00299-015-1889-3
Liu, X., Li, J., Xu, L., Wang, Q., and Lou, Y. (2018). Expressing OsMPK4 impairs plant growth but enhances the resistance of rice to the striped stem borer Chilo suppressalis. Int. J. Mol. Sci. 19:1182. doi: 10.3390/ijms19041182
Löffler, M., Schön, C. C., and Miedaner, T. (2009). Revealing the genetic architecture of FHB resistance in hexaploid wheat (Triticum aestivum L.) by QTL meta-analysis. Mol. Breed. 23, 473–488. doi: 10.1007/s11032-008-9250-y
Mahdi-Nezhad, N., Kamali, M. J., McIntyre, C. L., Fakheri, B. A., Omidi, M., and Masoudi, B. (2019). Mapping QTLs with main and epistatic effect on Seri ‘M82× Babax ‘wheat population under salt stress. Euphytica 215:130. doi: 10.1007/s10681-019-2450-1
Malik, N., Dwivedi, N., Singh, A. K., Parida, S. K., Agarwal, P., Thakur, J. K., et al. (2016). An integrated genomic strategy delineates candidate mediator genes regulating grain size and weight in rice. Sci. Rep. 6, 1–12. doi: 10.1038/srep23253
Martinez, A. K., Soriano, J. M., Tuberosa, R., Koumproglou, R., Jahrmann, T., and Salvi, S. (2016). Yield QTLome distribution correlates with gene density in maize. Plant Sci. 242, 300–309. doi: 10.1016/j.plantsci.2015.09.022
Marza, F., Bai, G. H., Carver, B. F., and Zhou, W. C. (2006). Quantitative trait loci for yield and related traits in the wheat population Ning7840× Clark. Theor. Appl. Genet. 112, 688–698. doi: 10.1007/s00122-005-0172-3
Mathew, I., Shimelis, H., Shayanowako, A. I. T., Laing, M., and Chaplot, V. (2019). Genome-wide association study of drought tolerance and biomass allocation in wheat. PLoS ONE 14:e0225383. doi: 10.1371/journal.pone.0225383
Michel, S., Löschenberger, F., Ametz, C., Pachler, B., Sparry, E., and Bürstmayr, H. (2019). Combining grain yield, protein content and protein quality by multi-trait genomic selection in bread wheat. Theor. Appl. Genet. 132, 2767–2780. doi: 10.1007/s00122-019-03386-1
Minx, P., Cordum, H., and Wilson, R. (2005). Sequence, annotation, and analysis of synteny between rice chromosome 3 and diverged grass species. Genome Res. 15, 1284–1291. doi: 10.1101/gr.3869505
Mishra, B. S., Jamsheer, K. M., Singh, D., Sharma, M., and Laxmi, A. (2017). Genome-wide identification and expression, protein–protein interaction and evolutionary analysis of the seed plant-specific BIG GRAIN and BIG GRAIN LIKE gene family. Front. Plant Sci. 8:1812. doi: 10.3389/fpls.2017.01812
Mohan, M., Nair, S., Bhagwat, A., Krishna, T. G., Yano, M., Bhatia, C. R., et al. (1997). Genome mapping, molecular markers and marker-assisted selection in crop plants. Mol. Breed. 3, 87–103. doi: 10.1023/A:1009651919792
Moin, M., Bakshi, A., Madhav, M. S., and Kirti, P. B. (2017). Expression profiling of ribosomal protein gene family in dehydration stress responses and characterization of transgenic rice plants overexpressing RPL23A for water-use efficiency and tolerance to drought and salt stresses. Front. Chem. 5:97. doi: 10.3389/fchem.2017.00097
Monod, C., Takahashi, Y., Goldschmidt-Clermont, M., and Rochaix, J. D. (1994). The chloroplast ycf8 open reading frame encodes a photosystem II polypeptide which maintains photosynthetic activity under adverse growth conditions. EMBO J. 13, 2747–2754. doi: 10.1002/j.1460-2075.1994.tb06568.x
Moore, G., Devos, K. M., Wang, Z., and Gale, M. D. (1995). Cereal genome evolution: grasses, line up and form a circle. Curr. Biol. 5, 737–739. doi: 10.1016/S0960-9822(95)00148-5
Morey, S. R., Hirose, T., Hashida, Y., Miyao, A., Hirochika, H., Ohsugi, R., et al. (2018). Genetic evidence for the role of a rice vacuolar invertase as a molecular sink strength determinant. Rice 11:6. doi: 10.1186/s12284-018-0201-x
Morrissey, J., and Guerinot, M. L. (2009). Iron uptake and transport in plants: the good, the bad, and the ionome. Chem. Rev. 109, 4553–4567. doi: 10.1021/cr900112r
Moschou, P. N., Sanmartin, M., Andriopoulou, A. H., Rojo, E., Sanchez-Serrano, J. J., and Roubelakis-Angelakis, K. A. (2008). Bridging the gap between plant and mammalian polyamine catabolism: a novel peroxisomal polyamine oxidase responsible for a full back-conversion pathway in Arabidopsis. Plant Physiol. 147, 1845–1857. doi: 10.1104/pp.108.123802
Moschou, P. N., Sarris, P. F., Skandalis, N., Andriopoulou, A. H., Paschalidis, K. A., Panopoulos, N. J., et al. (2009). Engineered polyamine catabolism preinduces tolerance of tobacco to bacteria and oomycetes. Plant Physiol. 149, 1970–1981. doi: 10.1104/pp.108.134932
Murray, C. J., and Lopez, A. D. (2013). Measuring the global burden of disease. N. Engl. J. Med. 369, 448–457. doi: 10.1056/NEJMra1201534
Muthu, V., Abbai, R., Nallathambi, J., Rahman, H., Ramasamy, S., Kambale, R., et al. (2020). Pyramiding QTLs controlling tolerance against drought, salinity, and submergence in rice through marker assisted breeding. PLoS ONE 15:e0227421. doi: 10.1371/journal.pone.0227421
Nagelkerke, N. J. (1991). A note on a general definition of the coefficient of determination. Biometrika 78, 691–692. doi: 10.1093/biomet/78.3.691
Ogawa, D., Yamamoto, E., Ohtani, T., Kanno, N., Tsunematsu, H., Nonoue, Y., et al. (2018). Haplotype-based allele mining in the Japan-MAGIC rice population. Sci. Rep. 8:4379. doi: 10.1038/s41598-018-22657-3
Ogiso-Tanaka, E., Matsubara, K., Yamamoto, S. I., Nonoue, Y., Wu, J., Fujisawa, H., et al. (2013). Natural variation of the RICE FLOWERING LOCUS T 1 contributes to flowering time divergence in rice. PLoS ONE 8:e75959. doi: 10.1371/journal.pone.0075959
Ozawa, S., Kobayashi, T., Sugiyama, R., Hoshida, H., Shiina, T., and Toyoshima, Y. (1997). Role of PSII-L protein (psbL gene product) on the electron transfer in photosystem II complex. 1. Over-production of wild-type and mutant versions of PSII-L protein and reconstitution into the PSII core complex. Plant Mol. Biol. 34, 151–161. doi: 10.1023/A:1005800909495
Pachauri, R. K., Allen, M. R., Barros, V. R., Broome, J., Cramer, W., Christ, R., et al. (2014). Climate change 2014: synthesis report. Contribution of Working Groups I, II and III to the fifth assessment report of the Intergovernmental Panel on Climate Change (p. 151). IPCC.
Peleg, Z., Cakmak, I., Ozturk, L., Yazici, A., Jun, Y., Budak, H., et al. (2009). Quantitative trait loci conferring grain mineral nutrient concentrations in durum wheat × wild emmer wheat RIL population. Theor. Appl. Genet. 119, 353–369. doi: 10.1007/s00122-009-1044-z
Pu, C. X., Han, Y. F., Zhu, S., Song, F. Y., Zhao, Y., Wang, C. Y., et al. (2017). The rice receptor-like kinases DWARF AND RUNTISH SPIKELET1 and 2 repress cell death and affect sugar utilization during reproductive development. Plant Cell 29, 70–89. doi: 10.1105/tpc.16.00218
Pu, Z. E., Ma, Y.u., He, Q. Y., Chen, G. Y., Wang, J. R., Liu, Y. X., et al. (2014). Quantitative trait loci associated with micronutrient concentrations in two recombinant inbred wheat lines. J. Integr. Agric. 13, 2322–2329. doi: 10.1016/S2095-3119(13)60640-1
Quraishi, U. M., Pont, C., Ain, Q. U., Flores, R., Burlot, L., Alaux, M., et al. (2017). Combined genomic and genetic data integration of major agronomical traits in bread wheat (Triticum aestivum L.). Front. Plant Sci. 8:1843. doi: 10.3389/fpls.2017.01843
Ram, H., Kaur, A., Gandass, N., Singh, S., Deshmukh, R., Sonah, H., et al. (2019). Molecular characterization and expression dynamics of MTP genes under various spatio-temporal stages and metal stress conditions in rice. PLoS ONE 14:e0217360. doi: 10.1371/journal.pone.0217360
Ramya, P., Chaubal, A., Kulkarni, K., Gupta, L., Kadoo, N., Dhaliwal, H. S., et al. (2010). QTL mapping of 1000-kernel weight, kernel length, and kernel width in bread wheat (Triticum aestivum L.). J. Appl. Genet. 51, 421–429. doi: 10.1007/BF03208872
Raza, Q., Riaz, A., Sabar, M., Atif, R. M., and Bashir, K. (2019). Meta-analysis of grain iron and zinc associated QTLs identified hotspot chromosomal regions and positional candidate genes for breeding biofortified rice. Plant Sci. 288:110214. doi: 10.1016/j.plantsci.2019.110214
Rebetzke, G. J., Condon, A. G., Farquhar, G. D., Appels, R., and Richards, R. A. (2008). Quantitative trait loci for carbon isotope discrimination are repeatable across environments and wheat mapping populations. Theor. Appl. Genet. 118, 123–137. doi: 10.1007/s00122-008-0882-4
Richmond, T. A., and Bleecker, A. B. (1999). A defect in β-oxidation causes abnormal inflorescence development in Arabidopsis. Plant Cell 11, 1911–1923. doi: 10.1105/tpc.11.10.1911
Rogalski, M., Schöttler, M. A., Thiele, W., Schulze, W. X., and Bock, R. (2008). Rpl33, a nonessential plastid-encoded ribosomal protein in tobacco, is required under cold stress conditions. Plant Cell 20, 2221–2237. doi: 10.1105/tpc.108.060392
Roshanzamir, H., Kordenaeej, A., and Bostani, A. (2013). Mapping QTLs related to Zn and Fe concentrations in bread wheat (Triticum aestivum) grain using microsatellite markers. Iranian J. Genet. Plant Breed. 2, 10–17.
Saini, D. K., Devi, P., and Kaushik, P. (2020). Advances in genomic interventions for wheat biofortification: a review. Agronomy 10:62. doi: 10.3390/agronomy10010062
Salvi, S., and Tuberosa, R. (2015). The crop QTLome comes of age. Curr. Opin. Biotechnol. 32, 179–185. doi: 10.1016/j.copbio.2015.01.001
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 Genomics 14, 1–6. doi: 10.1186/1471-2164-14-313
Shahzad, Z., Rouached, H., and Rakha, A. (2014). Combating mineral malnutrition through iron and zinc biofortification of cereals. Compr. Rev. Food Sci. Food Saf. 13, 329–346. doi: 10.1111/1541-4337.12063
Shariatipour, N., and Heidari, B. (2020). “Genetic-based biofortification of staple food crops to meet zinc and iron deficiency-related challenges,” in Plant Micronutrients, Deficiency and Toxicity Management, ed. T. Aftab, K. R. Hakeem (Cham: Springer), 173–223. doi: 10.1007/978-3-030-49856-6_8
Shukla, S., Singh, K., Patil, R. V., Kadam, S., Bharti, S., Prasad, P., et al. (2014). Genomic regions associated with grain yield under drought stress in wheat (Triticum aestivum L.). Euphytica 203, 449–467. doi: 10.1007/s10681-014-1314-y
Simmonds, N. W. (1995). The relation between yield and protein in cereal grain. J. Sci. Food Agric. 67, 309–315. doi: 10.1002/jsfa.2740670306
Somers, D. J., Isaac, P., and Edwards, K. (2004). A high-density microsatellite consensus map for bread wheat (Triticum aestivum L.). Theor. Appl. Genet. 109, 1105–1114. doi: 10.1007/s00122-004-1740-7
Soriano, J. M., and Alvaro, F. (2019). Discovering consensus genomic regions in wheat for root-related traits by QTL meta-analysis. Sci. Rep. 9:10537. doi: 10.1038/s41598-019-47038-2
Sosnowski, O., Charcosset, A., and Joets, J. (2012). BioMercator V3: an upgrade of genetic map compilation and quantitative trait loci meta-analysis algorithms. Bioinformatics 28, 2082–2083. doi: 10.1093/bioinformatics/bts313
Spindel, J. E., Begum, H., Akdemir, D., Collard, B., Redoña, E., Jannink, J. L., et al. (2016). Genome-wide prediction models that incorporate de novo GWAS are a powerful new tool for tropical rice improvement. Heredity 116, 395–408. doi: 10.1038/hdy.2015.113
Stein, A. J. (2010). Global impacts of human mineral malnutrition. Plant Soil 335, 133–154. doi: 10.1007/s11104-009-0228-2
Swamy, B. P. M., Kaladhar, K., Anuradha, K., Batchu, A. K., Longvah, T., and Sarla, N. (2018). QTL analysis for grain iron and zinc concentrations in two O. nivara derived backcross populations. Rice Sci. 25, 197–207. doi: 10.1016/j.rsci.2018.06.003
Swamy, B. P. M., and Sarla, N. (2011). Meta-analysis of yield QTLs derived from interspecific crosses of rice reveals consensus regions and candidate genes. Plant Mol. Biol. Rep. 29, 663–680. doi: 10.1007/s11105-010-0274-1
Tamiru, M., Undan, J. R., Takagi, H., Abe, A., Yoshida, K., Undan, J. Q., et al. (2015). A cytochrome P450, OsDSS1, is involved in growth and drought stress responses in rice (Oryza sativa L.). Plant Mol. Biol. 88, 85–99. doi: 10.1007/s11103-015-0310-5
Tiwari, C., Wallwork, H., Arun, B., Mishra, V. K., Velu, G., Stangoulis, J., et al. (2016). Molecular mapping of quantitative trait loci for zinc, iron and protein content in the grains of hexaploid wheat. Euphytica 207, 563–570. doi: 10.1007/s10681-015-1544-7
Tiwari, V. K., Rawat, N., Chhuneja, P., Neelam, K., Aggarwal, R., Randhawa, G. S., et al. (2009). Mapping of quantitative trait loci for grain iron and zinc concentration in diploid A genome wheat. J. Hered. 100, 771–776. doi: 10.1093/jhered/esp030
Truntzler, M., Barriere, Y., Sawkins, M. C., Lespinasse, D., Betran, J., Charcosset, A., et al. (2010). Meta-analysis of QTL involved in silage quality of maize and comparison with the position of candidate genes. Theor. Appl. Genet. 121, 1465–1482. doi: 10.1007/s00122-010-1402-x
Velu, G., Guzman, C., Mondal, S., Autrique, J. E., Huerta, J., and Singh, R. P. (2016). Effect of drought and elevated temperature on grain zinc and iron concentrations in CIMMYT spring wheat. J. Cereal Sci. 69, 182–186. doi: 10.1016/j.jcs.2016.03.006
Velu, G., Singh, R. P., Crespo-Herrera, L., Juliana, P., Dreisigacker, S., Valluru, R., et al. (2018). Genetic dissection of grain zinc concentration in spring wheat for mainstreaming biofortification in CIMMYT wheat breeding. Sci. Rep. 8:13526. doi: 10.1038/s41598-018-31951-z
Velu, G., Tutus, Y., Gomez-Becerra, H. F., Hao, Y., Demir, L., Kara, R., et al. (2017). QTL mapping for grain zinc and iron concentrations and zinc efficiency in a tetraploid and hexaploid wheat mapping populations. Plant Soil 411, 81–99. doi: 10.1007/s11104-016-3025-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
Vij, S., Giri, J., Dansana, P. K., Kapoor, S., and Tyagi, A. K. (2008). The receptor-like cytoplasmic kinase (OsRLCK) gene family in rice: organization, phylogenetic relationship, and expression during development and stress. Mol. Plant 1, 732–750. doi: 10.1093/mp/ssn047
Vijayalakshmi, K., Fritz, A. K., Paulsen, G. M., Bai, G., Pandravada, S., and Gill, B. S. (2010). Modeling and mapping QTL for senescence-related traits in winter wheat under high temperature. Mol. Breed. 26, 163–175. doi: 10.1007/s11032-009-9366-8
Wan, Y., King, R., Mitchell, R. A., Hassani-Pak, K., and Hawkesford, M. J. (2017). Spatiotemporal expression patterns of wheat amino acid transporters reveal their putative roles in nitrogen transport and responses to abiotic stress. Sci. Rep. 7:5461. doi: 10.1038/s41598-017-04473-3
Wang, L. I., Cui, F. A., Wang, J., Jun, L. I., Ding, A., Zhao, C., et al. (2012). Conditional QTL mapping of protein content in wheat with respect to grain yield and its components. J. Genet. 91, 303–312. doi: 10.1007/s12041-012-0190-2
Wang, P., Xing, Y., Li, Z., and Yu, S. (2012). Improving rice yield and quality by QTL pyramiding. Mol. Breed. 29, 903–913. doi: 10.1007/s11032-011-9679-2
Wang, Y., Huang, Z., Deng, D., Ding, H., Zhang, R., Wang, S., et al. (2013). Meta-analysis combined with syntenic metaQTL mining dissects candidate loci for maize yield. Mol. Breed. 31, 601–614. doi: 10.1007/s11032-012-9818-4
Wang, Y., Shi, C., Yang, T., Zhao, L., Chen, J., Zhang, N., et al. (2018). High-throughput sequencing revealed that microRNAs were involved in the development of superior and inferior grains in bread wheat. Sci. Rep. 8, 1–18. doi: 10.1038/s41598-018-31870-z
Wang, Y., Xu, J., Deng, D., Ding, H., Bian, Y., Yin, Z., et al. (2016). A comprehensive meta-analysis of plant morphology, yield, stay-green, and virus disease resistance QTL in maize (Zea mays L.). Planta 243, 459–471. doi: 10.1007/s00425-015-2419-9
Wu, C., You, C., Li, C., Long, T., Chen, G., Byrne, M. E., et al. (2008). RID1, encoding a Cys2/His2-type zinc finger transcription factor, acts as a master switch from vegetative to floral development in rice. Proc. Natl. Acad. Sci. USA. 105, 12915–12920. doi: 10.1073/pnas.0806019105
Wu, X. L., and Hu, Z. L. (2012). Meta-analysis of QTL mapping experiments. Methods Mol. Biol. 871, 145–171. doi: 10.1007/978-1-61779-785-9_8
Wu, Y., Zhao, S., Li, X., Zhang, B., Jiang, L., Tang, Y., et al. (2018). Deletions linked to PROG1 gene participate in plant architecture domestication in Asian and African rice. Nat. Commun. 9, 1–10. doi: 10.1038/s41467-018-06509-2
Xiang, K., Reid, L. M., Zhang, Z. M., Zhu, X. Y., and Pan, G. T. (2012). Characterization of correlation between grain moisture and ear rot resistance in maize by QTL meta-analysis. Euphytica 183, 185–195. doi: 10.1007/s10681-011-0440-z
Xu, L., Zhao, H., Ruan, W., Deng, M., Wang, F., Peng, J., et al. (2017). ABNORMAL INFLORESCENCE MERISTEM1 functions in salicylic acid biosynthesis to maintain proper reactive oxygen species levels for root meristem activity in rice. Plant Cell 29, 560–574. doi: 10.1105/tpc.16.00665
Xu, Y., An, D., Liu, D., Zhang, A., Xu, H., and Li, B. (2012). Molecular mapping of QTLs for grain zinc, iron and protein concentration of wheat across two environments. Field Crops Res. 138, 57–62. doi: 10.1016/j.fcr.2012.09.017
Xue, W., Xing, Y., Weng, X., Zhao, Y., Tang, W., Wang, L., et al. (2008). Natural variation in Ghd7 is an important regulator of heading date and yield potential in rice. Nat. Genet. 40, 761–767. doi: 10.1038/ng.143
Yamburenko, M. V., Kieber, J. J., and Schaller, G. E. (2017). Dynamic patterns of expression for genes regulating cytokinin metabolism and signaling during rice inflorescence development. PLoS ONE 12:e0176060. doi: 10.1371/journal.pone.0176060
Yao, R., Wang, L., Li, Y., Chen, L., Li, S., Du, X., et al. (2018). Rice DWARF14 acts as an unconventional hormone receptor for strigolactone. J. Exp. Bot. 69, 2355–2365. doi: 10.1093/jxb/ery014
Yu, Y., Ouyang, Y., and Yao, W. (2017). shinyCircos: an R/Shiny application for interactive creation of Circos plot. Bioinformatics 34, 1229–1231. doi: 10.1093/bioinformatics/btx763
Zhang, A., Zhang, J., Ye, N., Cao, J., Tan, M., Zhang, J., et al. (2010). ZmMPK5 is required for the NADPH oxidase-mediated self-propagation of apoplastic H2O2 in brassinosteroid-induced antioxidant defence in leaves of maize. J. Exp. Bot. 61, 4399–4411. doi: 10.1093/jxb/erq243
Zhang, D., Jiang, S., Pan, J., Kong, X., Zhou, Y., Liu, Y., et al. (2014). The overexpression of a maize mitogen-activated protein kinase gene (ZmMPK5) confers salt stress tolerance and induces defence responses in tobacco. Plant Biol. 16, 558–570. doi: 10.1111/plb.12084
Zhang, H., Zhang, J., Yan, J., Gou, F., Mao, Y., Tang, G., et al. (2017). Short tandem target mimic rice lines uncover functions of miRNAs in regulating important agronomic traits. Proc. Natl. Acad. Sci. USA. 114, 5277–5282. doi: 10.1073/pnas.1703752114
Zhang, L. Y., Liu, D. C., Guo, X. L., Yang, W. L., Sun, J. Z., Wang, D. W., et al. (2010). Genomic distribution of quantitative trait loci for yield and yield-related traits in common wheat. J. Integr. Plant. Biol. 52, 996–1007. doi: 10.1111/j.1744-7909.2010.00967.x
Zhang, M., Cui, Y., Liu, Y. H., Xu, W., Sze, S. H., Murray, S. C., et al. (2020). Accurate prediction of maize grain yield using its contributing genes for gene-based breeding. Genomics 112, 225–236. doi: 10.1016/j.ygeno.2019.02.001
Zhang, M., and Liu, B. (2017). Identification of a rice metal tolerance protein OsMTP11 as a manganese transporter. PLoS ONE 12:e0174987. doi: 10.1371/journal.pone.0174987
Zhang, T., Feng, P., Li, Y., Yu, P., Yu, G., Sang, X., et al. (2018). VIRESCENT-ALBINO LEAF 1 regulates leaf colour development and cell division in rice. J. Exp. Bot. 69, 4791–4804. doi: 10.1093/jxb/ery250
Zhao, H., Duan, K. X., Ma, B., Yin, C. C., Hu, Y., Tao, J. J., et al. (2020). Histidine kinase MHZ1/OsHK1 interacts with ethylene receptors to regulate root growth in rice. Nat. Commun. 11, 1–13. doi: 10.1038/s41467-020-14313-0
Zhao, X., Peng, Y., Zhang, J., Fang, P., and Wu, B. (2018). Identification of QTLs and meta-QTLs for seven agronomic traits in multiple maize populations under well-watered and water-stressed conditions. Crop Sci. 58, 507–520. doi: 10.2135/cropsci2016.12.0991
Zhao, Y., Cheng, S., Song, Y., Huang, Y., Zhou, S., Liu, X., et al. (2015). The interaction between rice ERF3 and WOX11 promotes crown root development by regulating gene expression involved in cytokinin signaling. Plant Cell 27, 2469–2483. doi: 10.1105/tpc.15.00227
Zhou, W., Wang, Y., Wu, Z., Luo, L., Liu, P., Yan, L., et al. (2016). Homologs of SCAR/WAVE complex components are required for epidermal cell morphogenesis in rice. J. Exp. Bot. 67, 4311–4323. doi: 10.1093/jxb/erw214
Zhu, D., Chang, Y., Pei, T., Zhang, X., Liu, L., Li, Y., et al. (2020). MAPK-like protein 1 positively regulates maize seedling drought sensitivity by suppressing ABA biosynthesis. Plant J. 102, 747–760. doi: 10.1111/tpj.14660
Zilic, S., Barac, M., Pešic, M., Dodig, D., and Ignjatovic-Micic, D. (2011). Characterization of proteins from grain of different bread and durum wheat genotypes. Int. J. Mol. Sci. 12, 5878–5894. doi: 10.3390/ijms12095878
Keywords: meta QTL analysis, comparative genomics, iron, zinc, candidate gene, GWAS
Citation: Shariatipour N, Heidari B, Tahmasebi A and Richards C (2021) Comparative Genomic Analysis of Quantitative Trait Loci Associated With Micronutrient Contents, Grain Quality, and Agronomic Traits in Wheat (Triticum aestivum L.). Front. Plant Sci. 12:709817. doi: 10.3389/fpls.2021.709817
Received: 14 May 2021; Accepted: 06 September 2021;
Published: 12 October 2021.
Edited by:
Leonardo Abdiel Crespo Herrera, International Maize and Wheat Improvement Center, MexicoReviewed by:
Velu Govindan, International Maize and Wheat Improvement Center, MexicoRoi Ben-David, Agricultural Research Organization (ARO), Israel
Copyright © 2021 Shariatipour, Heidari, Tahmasebi and Richards. 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: Bahram Heidari, YmhlaWRhcmkmI3gwMDA0MDtzaGlyYXp1LmFjLmly
†ORCID: Nikwan Shariatipour orcid.org/0000-0003-4174-4375
Bahram Heidari orcid.org/0000-0002-5856-4592
Christopher Richards orcid.org/0000-0002-9978-6079