Meta-QTLs, ortho-MQTLs, and candidate genes for thermotolerance in wheat (Triticum aestivum L.)

Physical map of QTL for eleven agronomic traits across fifteen environments, identification of related candidate genes, and development of KASP markers with emphasis on terminal heat stress tolerance in common wheat

KEY MESSAGE

Key message This study identified stable QTL, promising candidate genes and developed novel KASP markers for heat tolerance, providing genomic resources to assist breeding for the development of high-yielding and heat-tolerant wheat germplasm and varieties. To understand the genetic architecture of eleven agronomic traits under heat stress, we used a doubled-haploid population (177 lines) derived from a heat-sensitive cultivar (PBW343) and a heat-tolerant genotype (KSG1203). This population was evaluated under timely, late and very late sown conditions over locations and years comprising fifteen environments. Best linear unbiased estimates and a genetic map (5,710 SNPs) developed using sequencing-based genotyping were used for QTL mapping. The identified 66 QTL (20 novel) were integrated into wheat physical map (14,263.4 Mb). These QTL explained 5.3% (QDth.ccsu-4A for days to heading and QDtm.ccsu-5B for days to maturity) to 24.9% (QGfd.ccsu-7D for grain filling duration) phenotypic variation. Thirteen stable QTL explaining high phenotypic variation were recommended for marker-assisted recurrent selection (MARS) for optimum/heat stress environments. Selected QTL were validated by their presence in high-yielding doubled-haploid lines. Some QTL for 1000-grain weight (TaERF3-3B, TaFER-5B, and TaZIM-A1), grain yield (TaCol-B5), and developmental traits (TaVRT-2) were co-localized with known genes. Specific known genes for traits like abiotic/biotic stress, grain quality and yield were co-located with 26 other QTL. Furthermore, 209 differentially expressed candidate genes for heat tolerance in plants that encode 28 different proteins were identified. KASP markers for three major/stable QTL, namely QGfd.ccsu-7A for grain filling duration on chromosome 7A (timely sown), QNgs.ccsu-3A for number of grains per spike on 3A, and QDth.ccsu-7A for days to heading on 7A (late and very late sown) environments were developed for MARS focusing on the development of heat-tolerant wheat varieties/germplasm.

Unveiling the role of the ERD15 gene in wheat's tolerance to combined drought and salinity stress: a meta-analysis of QTL and RNA-Seq data

The coexistence of drought and salinity stresses in field conditions significantly hinders wheat (Triticum aestivum L.) productivity. Understanding the molecular mechanisms governing response and tolerance to these stresses is crucial for developing resilient wheat varieties. Our research, employing a combination of meta-QTL and meta-RNA-Seq transcriptome analyses, has uncovered the genome functional landscape of wheat in response to drought and salinity. We identified 118 meta-QTLs (MQTLs) distributed across all 21 wheat chromosomes, with ten designated as the most promising. Additionally, we found 690 meta-differentially expressed genes (mDEGs) shared between drought and salinity stress. Notably, our findings highlight the Early Responsive to Dehydration 15 (ERD15) gene, located in one of the most promising MQTLs, as a key gene in the shared gene network of drought and salinity stress. ERD15, differentially expressed between contrasting wheat genotypes under combined stress conditions, significantly regulates water relations, photosynthetic activity, antioxidant activity, and ion homeostasis. These findings not only provide valuable insights into the molecular genetic mechanisms underlying combined stress tolerance in wheat but also hold the potential to contribute significantly to the development of stress-resilient wheat varieties.

Insights into maydis leaf blight resistance in maize: a comprehensive genome-wide association study in sub-tropics of India

BACKGROUND

In the face of contemporary climatic vulnerabilities and escalating global temperatures, the prevalence of maydis leaf blight (MLB) poses a potential threat to maize production. This study endeavours to discern marker-trait associations and elucidate the candidate genes that underlie resistance to MLB in maize by employing a diverse panel comprising 336 lines. The panel was screening for MLB across four environments, employing standard artificial inoculation techniques. Genome-wide association studies (GWAS) and haplotype analysis were conducted utilizing a total of 128,490 SNPs obtained from genotyping-by-sequencing (GBS).

RESULTS

GWAS identified 26 highly significant SNPs associated with MLB resistance, among the markers examined. Seven of these SNPs, reported in novel chromosomal bins (9.06, 5.01, 9.01, 7.04, 4.06, 1.04, and 6.05) were associated with genes: bzip23, NAGS1, CDPK7, aspartic proteinase NEP-2, VQ4, and Wun1, which were characterized for their roles in diminishing fungal activity, fortifying defence mechanisms against necrotrophic pathogens, modulating phyto-hormone signalling, and orchestrating oxidative burst responses. Gene mining approach identified 22 potential candidate genes associated with SNPs due to their functional relevance to resistance against necrotrophic pathogens. Notably, bin 8.06, which hosts five SNPs, showed a connection to defense-regulating genes against MLB, indicating the potential formation of a functional gene cluster that triggers a cascade of reactions against MLB. In silico studies revealed gene expression levels exceeding ten fragments per kilobase million (FPKM) for most genes and demonstrated coexpression among all candidate genes in the coexpression network. Haplotype regression analysis revealed the association of 13 common significant haplotypes at Bonferroni ≤ 0.05. The phenotypic variance explained by these significant haplotypes ranged from low to moderate, suggesting a breeding strategy that combines multiple resistance alleles to enhance resistance to MLB. Additionally, one particular haplotype block (Hap_8.3) was found to consist of two SNPs (S8_152715134, S8_152460815) identified in GWAS with 9.45% variation explained (PVE).

CONCLUSION

The identified SNPs/ haplotypes associated with the trait of interest contribute to the enrichment of allelic diversity and hold direct applicability in Genomics Assisted Breeding for enhancing MLB resistance in maize.

Exploring the genetic architecture of powdery mildew resistance in wheat through QTL meta-analysis

Powdery mildew (PM), caused by Blumeria graminis f. sp. tritici, poses a significant threat to wheat production, necessitating the development of genetically resistant varieties for long-term control. Therefore, exploring genetic architecture of PM in wheat to uncover important genomic regions is an important area of wheat research. In recent years, the utilization of meta-QTL (MQTL) analysis has gained prominence as an essential tool for unraveling the complex genetic architecture underlying complex quantitative traits. The aim of this research was to conduct a QTL meta-analysis to pinpoint the specific genomic regions in wheat responsible for governing PM resistance. This study integrated 222 QTLs from 33 linkage-based studies using a consensus map with 54,672 markers. The analysis revealed 39 MQTLs, refined to 9 high-confidence MQTLs (hcMQTLs) with confidence intervals of 0.49 to 12.94 cM. The MQTLs had an average physical interval of 41.00 Mb, ranging from 0.000048 Mb to 380.71 Mb per MQTL. Importantly, 18 MQTLs co-localized with known resistance genes like Pm2, Pm3, Pm8, Pm21, Pm38, and Pm41. The study identified 256 gene models within hcMQTLs, providing potential targets for marker-assisted breeding and genomic prediction programs to enhance PM resistance. These MQTLs would serve as a foundation for fine mapping, gene isolation, and functional genomics studies, facilitating a deeper understanding of molecular mechanisms. The identification of candidate genes opens up exciting possibilities for the development of PM-resistant wheat varieties after validation.

Meta-QTL analysis and identification of candidate genes for multiple-traits associated with spot blotch resistance in bread wheat

In bread wheat, a literature search gave 228 QTLs for six traits, including resistance against spot blotch and the following five other related traits: (i) stay green; (ii) flag leaf senescence; (iii) green leaf area duration; (iv) green leaf area of the main stem; and (v) black point resistance. These QTLs were used for metaQTL (MQTL) analysis. For this purpose, a consensus map with 72,788 markers was prepared; 69 of the above 228 QTLs, which were suitable for MQTL analysis, were projected on the consensus map. This exercise resulted in the identification of 16 meta-QTLs (MQTLs) located on 11 chromosomes, with the PVE ranging from 5.4% (MQTL7) to 21.8% (MQTL5), and the confidence intervals ranging from 1.5 to 20.7 cM (except five MQTLs with a range of 36.1-57.8 cM). The number of QTLs associated with individual MQTLs ranged from a maximum of 17 in MQTL3 to 8 each in MQTL5 and MQTL8 and 5 each in MQTL7 and MQTL14. The 16 MQTLs, included 12 multi-trait MQTLs; one of the MQTL also overlapped a genomic region carrying the major spot blotch resistance gene Sb1. Of the total 16 MQTLs, 12 MQTLs were also validated through marker-trait associations that were available from earlier genome-wide association studies. The genomic regions associated with MQTLs were also used for the identification of candidate genes (CGs) and led to the identification of 516 CGs encoding 508 proteins; 411 of these proteins are known to be associated with resistance against several biotic stresses. In silico expression analysis of CGs using transcriptome data allowed the identification of 71 differentially expressed CGs, which were examined for further possible studies. The findings of the present study should facilitate fine-mapping and cloning of genes, enabling Marker Assisted Selection.

Meta QTL analysis for dissecting abiotic stress tolerance in chickpea

BACKGROUND

Chickpea is prone to many abiotic stresses such as heat, drought, salinity, etc. which cause severe loss in yield. Tolerance towards these stresses is quantitative in nature and many studies have been done to map the loci influencing these traits in different populations using different markers. This study is an attempt to meta-analyse those reported loci projected over a high-density consensus map to provide a more accurate information on the regions influencing heat, drought, cold and salinity tolerance in chickpea.

RESULTS

A meta-analysis of QTL reported to be responsible for tolerance to drought, heat, cold and salinity stress tolerance in chickpeas was done. A total of 1512 QTL responsible for the concerned abiotic stress tolerance were collected from literature, of which 1189 were projected on a chickpea consensus genetic map. The QTL meta-analysis predicted 59 MQTL spread over all 8 chromosomes, responsible for these 4 kinds of abiotic stress tolerance in chickpea. The physical locations of 23 MQTL were validated by various marker-trait associations and genome-wide association studies. Out of these reported MQTL, CaMQAST1.1, CaMQAST4.1, CaMQAST4.4, CaMQAST7.8, and CaMQAST8.2 were suggested to be useful for different breeding approaches as they were responsible for high per cent variance explained (PVE), had small intervals and encompassed a large number of originally reported QTL. Many putative candidate genes that might be responsible for directly or indirectly conferring abiotic stress tolerance were identified in the region covered by 4 major MQTL- CaMQAST1.1, CaMQAST4.4, CaMQAST7.7, and CaMQAST6.4, such as heat shock proteins, auxin and gibberellin response factors, etc. CONCLUSION: The results of this study should be useful for the breeders and researchers to develop new chickpea varieties which are tolerant to drought, heat, cold, and salinity stresses.

Integrated meta-analysis and transcriptomics pinpoint genomic loci and novel candidate genes associated with submergence tolerance in rice

BACKGROUND

Due to rising costs, water shortages, and labour shortages, farmers across the globe now prefer a direct seeding approach. However, submergence stress remains a major bottleneck limiting the success of this approach in rice cultivation. The merger of accumulated rice genetic resources provides an opportunity to detect key genomic loci and candidate genes that influence the flooding tolerance of rice.

RESULTS

In the present study, a whole-genome meta-analysis was conducted on 120 quantitative trait loci (QTL) obtained from 16 independent QTL studies reported from 2004 to 2023. These QTL were confined to 18 meta-QTL (MQTL), and ten MQTL were successfully validated by independent genome-wide association studies from diverse natural populations. The mean confidence interval (CI) of the identified MQTL was 3.44 times narrower than the mean CI of the initial QTL. Moreover, four core MQTL loci with genetic distance less than 2 cM were obtained. By combining differentially expressed genes (DEG) from two transcriptome datasets with 858 candidate genes identified in the core MQTL regions, we found 38 common differentially expressed candidate genes (DECGs). In silico expression analysis of these DECGs led to the identification of 21 genes with high expression in embryo and coleoptile under submerged conditions. These DECGs encode proteins with known functions involved in submergence tolerance including WRKY, F-box, zinc fingers, glycosyltransferase, protein kinase, cytochrome P450, PP2C, hypoxia-responsive family, and DUF domain. By haplotype analysis, the 21 DECGs demonstrated distinct genetic differentiation and substantial genetic distance mainly between indica and japonica subspecies. Further, the MQTL7.1 was successfully validated using flanked marker S2329 on a set of genotypes with phenotypic variation.

CONCLUSION

This study provides a new perspective on understanding the genetic basis of submergence tolerance in rice. The identified MQTL and novel candidate genes lay the foundation for marker-assisted breeding/engineering of flooding-tolerant cultivars conducive to direct seeding.

Association study of morpho-phenological traits in quinoa (Chenopodium quinoa Willd.) using SSR markers

In this study, the genetic and molecular diversity of 60 quinoa accessions was assessed using agronomically important traits related to grain yield as well as microsatellite (SSR) markers, and informative markers linked to the studied traits were identified using association study. The results showed that most of the studied traits had a relatively high diversity, but grain saponin and protein content showed the highest diversity. High diversity was also observed in all SSR markers, but KAAT023, KAAT027, KAAT036, and KCAA014 showed the highest values for most of the diversity indices and can be introduced as the informative markers to assess genetic diversity in quinoa. Population structure analysis showed that the studied population probably includes two subclusters, so that out of 60 quinoa accessions, 29 (48%) and 23 (38%) accessions were assigned to the first and second subclusters, respectively, and eight (13%) accessions were considered as the mixed genotypes. The study of the population structure using Structure software showed two possible subgroups (K = 2) in the studied population and the results of the bar plot confirmed it. Association study using the general linear model (GLM) and mixed linear model (MLM) identified the number of 35 and 32 significant marker-trait associations (MTAs) for the first year (2019) and 37 and 35 significant MTAs for the second year (2020), respectively. Among the significant MTAs identified for different traits, the highest number of significant MTAs were obtained for grain yield and 1000-grain weight with six and five MTAs, respectively.

Meta-QTL analysis in wheat: progress, challenges and opportunities

Wheat, an important cereal crop globally, faces major challenges due to increasing global population and changing climates. The production and productivity are challenged by several biotic and abiotic stresses. There is also a pressing demand to enhance grain yield and quality/nutrition to ensure global food and nutritional security. To address these multifaceted concerns, researchers have conducted numerous meta-QTL (MQTL) studies in wheat, resulting in the identification of candidate genes that govern these complex quantitative traits. MQTL analysis has successfully unraveled the complex genetic architecture of polygenic quantitative traits in wheat. Candidate genes associated with stress adaptation have been pinpointed for abiotic and biotic traits, facilitating targeted breeding efforts to enhance stress tolerance. Furthermore, high-confidence candidate genes (CGs) and flanking markers to MQTLs will help in marker-assisted breeding programs aimed at enhancing stress tolerance, yield, quality and nutrition. Functional analysis of these CGs can enhance our understanding of intricate trait-related genetics. The discovery of orthologous MQTLs shared between wheat and other crops sheds light on common evolutionary pathways governing these traits. Breeders can leverage the most promising MQTLs and CGs associated with multiple traits to develop superior next-generation wheat cultivars with improved trait performance. This review provides a comprehensive overview of MQTL analysis in wheat, highlighting progress, challenges, validation methods and future opportunities in wheat genetics and breeding, contributing to global food security and sustainable agriculture.

Unravelling the genetic framework associated with grain quality and yield-related traits in maize (Zea mays L.)

Maize serves as a crucial nutrient reservoir for a significant portion of the global population. However, to effectively address the growing world population's hidden hunger, it is essential to focus on two key aspects: biofortification of maize and improving its yield potential through advanced breeding techniques. Moreover, the coordination of multiple targets within a single breeding program poses a complex challenge. This study compiled mapping studies conducted over the past decade, identifying quantitative trait loci associated with grain quality and yield related traits in maize. Meta-QTL analysis of 2,974 QTLs for 169 component traits (associated with quality and yield related traits) revealed 68 MQTLs across different genetic backgrounds and environments. Most of these MQTLs were further validated using the data from genome-wide association studies (GWAS). Further, ten MQTLs, referred to as breeding-friendly MQTLs (BF-MQTLs), with a significant phenotypic variation explained over 10% and confidence interval less than 2 Mb, were shortlisted. BF-MQTLs were further used to identify potential candidate genes, including 59 genes encoding important proteins/products involved in essential metabolic pathways. Five BF-MQTLs associated with both quality and yield traits were also recommended to be utilized in future breeding programs. Synteny analysis with wheat and rice genomes revealed conserved regions across the genomes, indicating these hotspot regions as validated targets for developing biofortified, high-yielding maize varieties in future breeding programs. After validation, the identified candidate genes can also be utilized to effectively model the plant architecture and enhance desirable quality traits through various approaches such as marker-assisted breeding, genetic engineering, and genome editing.

Identification and characterization of RuvBL DNA helicase genes for tolerance against abiotic stresses in bread wheat (Triticum aestivum L.) and related species

Recombination UVB (sensitivity) like (RuvBL) helicase genes represent a conserved family of genes, which are known to be involved in providing tolerance against abiotic stresses like heat and drought. We identified nine wheat RuvBL genes, one each on nine different chromosomes, belonging to homoeologous groups 2, 3, and 4. The lengths of genes ranged from 1647 to 2197 bp and exhibited synteny with corresponding genes in related species including Ae. tauschii, Z. mays, O. sativa, H. vulgare, and B. distachyon. The gene sequences were associated with regulatory cis-elements and transposable elements. Two genes, namely TaRuvBL1a-4A and TaRuvBL1a-4B, also carried targets for a widely known miRNA, tae-miR164. Gene ontology revealed that these genes were closely associated with ATP-dependent formation of histone acetyltransferase complex. Analysis of the structure and function of RuvBL proteins revealed that the proteins were localized mainly in the cytoplasm. A representative gene, namely TaRuvBL1a-4A, was also shown to be involved in protein-protein interactions with ten other proteins. On the basis of phylogeny, RuvBL proteins were placed in two sub-divisions, namely RuvBL1 and RuvBL2, which were further classified into clusters and sub-clusters. In silico studies suggested that these genes were differentially expressed under heat/drought. The qRT-PCR analysis confirmed that expression of TaRuvBL genes differed among wheat cultivars, which differed in the level of thermotolerance. The present study advances our understanding of the biological role of wheat RuvBL genes and should help in planning future studies on RuvBL genes in wheat including use of RuvBL genes in breeding thermotolerant wheat cultivars.

Genome-wide association mapping of genomic regions associated with drought stress tolerance at seedling and reproductive stages in bread wheat

Understanding the genetic architecture of drought stress tolerance in bread wheat at seedling and reproductive stages is crucial for developing drought-tolerant varieties. In the present study, 192 diverse wheat genotypes, a subset from the Wheat Associated Mapping Initiative (WAMI) panel, were evaluated at the seedling stage in a hydroponics system for chlorophyll content (CL), shoot length (SLT), shoot weight (SWT), root length (RLT), and root weight (RWT) under both drought and optimum conditions. Following that, a genome-wide association study (GWAS) was carried out using the phenotypic data recorded during the hydroponics experiment as well as data available from previously conducted multi-location field trials under optimal and drought stress conditions. The panel had previously been genotyped using the Infinium iSelect 90K SNP array with 26,814 polymorphic markers. Using single as well as multi-locus models, GWAS identified 94 significant marker-trait associations (MTAs) or SNPs associated with traits recorded at the seedling stage and 451 for traits recorded at the reproductive stage. The significant SNPs included several novel, significant, and promising MTAs for different traits. The average LD decay distance for the whole genome was approximately 0.48 Mbp, ranging from 0.07 Mbp (chromosome 6D) to 4.14 Mbp (chromosome 2A). Furthermore, several promising SNPs revealed significant differences among haplotypes for traits such as RLT, RWT, SLT, SWT, and GY under drought stress. Functional annotation and in silico expression analysis revealed important putative candidate genes underlying the identified stable genomic regions such as protein kinases, O-methyltransferases, GroES-like superfamily proteins, NAD-dependent dehydratases, etc. The findings of the present study may be useful for improving yield potential, and stability under drought stress conditions.

Comprehensive meta-QTL analysis for dissecting the genetic architecture of stripe rust resistance in bread wheat

BACKGROUND

Yellow or stripe rust, caused by the fungus Puccinia striiformis f. sp. tritici (Pst) is an important disease of wheat that threatens wheat production. Since developing resistant cultivars offers a viable solution for disease management, it is essential to understand the genetic basis of stripe rust resistance. In recent years, meta-QTL analysis of identified QTLs has gained popularity as a way to dissect the genetic architecture underpinning quantitative traits, including disease resistance.

RESULTS

Systematic meta-QTL analysis involving 505 QTLs from 101 linkage-based interval mapping studies was conducted for stripe rust resistance in wheat. For this purpose, publicly available high-quality genetic maps were used to create a consensus linkage map involving 138,574 markers. This map was used to project the QTLs and conduct meta-QTL analysis. A total of 67 important meta-QTLs (MQTLs) were identified which were refined to 29 high-confidence MQTLs. The confidence interval (CI) of MQTLs ranged from 0 to 11.68 cM with a mean of 1.97 cM. The mean physical CI of MQTLs was 24.01 Mb, ranging from 0.0749 to 216.23 Mb per MQTL. As many as 44 MQTLs colocalized with marker-trait associations or SNP peaks associated with stripe rust resistance in wheat. Some MQTLs also included the following major genes- Yr5, Yr7, Yr16, Yr26, Yr30, Yr43, Yr44, Yr64, YrCH52, and YrH52. Candidate gene mining in high-confidence MQTLs identified 1,562 gene models. Examining these gene models for differential expressions yielded 123 differentially expressed genes, including the 59 most promising CGs. We also studied how these genes were expressed in wheat tissues at different phases of development.

CONCLUSION

The most promising MQTLs identified in this study may facilitate marker-assisted breeding for stripe rust resistance in wheat. Information on markers flanking the MQTLs can be utilized in genomic selection models to increase the prediction accuracy for stripe rust resistance. The candidate genes identified can also be utilized for enhancing the wheat resistance against stripe rust after in vivo confirmation/validation using one or more of the following methods: gene cloning, reverse genetic methods, and omics approaches.

Meta-QTLs, ortho-meta QTLs and related candidate genes for yield and its component traits under water stress in wheat (Triticum aestivum L.)

Meta-QTLs (MQTLs), ortho-MQTLs, and related candidate genes (CGs) for yield and its seven component traits evaluated under water deficit conditions were identified in wheat. For this purpose, a high density consensus map and 318 known QTLs were used for identification of 56 MQTLs. Confidence intervals (CIs) of the MQTLs were narrower (0.7-21 cM; mean = 5.95 cM) than the CIs of the known QTLs (0.4-66.6 cM; mean = 12.72 cM). Forty-seven MQTLs were co-located with marker trait associations reported in previous genome-wide association studies. Nine selected MQTLs were declared as 'breeders MQTLs' for use in marker-assisted breeding (MAB). Utilizing known MQTLs and synteny/collinearity among wheat, rice and maize, 12 ortho-MQTLs were also identified. A total of 1497 CGs underlying MQTLs were also identified, which were subjected to in-silico expression analysis, leading to identification of 64 differentially expressed CGs (DECGs) under normal and water deficit conditions. These DECGs encoded a variety of proteins, including the following: zinc finger, cytochrome P450, AP2/ERF domain-containing proteins, plant peroxidase, glycosyl transferase, glycoside hydrolase. The expression of 12 CGs at seedling stage (3 h stress) was validated using qRT-PCR in two wheat genotypes, namely Excalibur (drought tolerant) and PBW343 (drought sensitive). Nine of the 12 CGs were up-regulated and three down-regulated in Excalibur. The results of the present study should prove useful for MAB, for fine mapping of promising MQTLs and for cloning of genes across the three cereals studied.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12298-023-01301-z.

Genome-Wide Meta-Analysis of QTLs Associated with Root Traits and Implications for Maize Breeding

Root system architecture (RSA), also known as root morphology, is critical in plant acquisition of soil resources, plant growth, and yield formation. Many QTLs associated with RSA or root traits in maize have been identified using several bi-parental populations, particularly in response to various environmental factors. In the present study, a meta-analysis of QTLs associated with root traits was performed in maize using 917 QTLs retrieved from 43 mapping studies published from 1998 to 2020. A total of 631 QTLs were projected onto a consensus map involving 19,714 markers, which led to the prediction of 68 meta-QTLs (MQTLs). Among these 68 MQTLs, 36 MQTLs were validated with the marker-trait associations available from previous genome-wide association studies for root traits. The use of comparative genomics approaches revealed several gene models conserved among the maize, sorghum, and rice genomes. Among the conserved genomic regions, the ortho-MQTL analysis uncovered 20 maize MQTLs syntenic to 27 rice MQTLs for root traits. Functional analysis of some high-confidence MQTL regions revealed 442 gene models, which were then subjected to in silico expression analysis, yielding 235 gene models with significant expression in various tissues. Furthermore, 16 known genes viz., DXS2, PHT, RTP1, TUA4, YUC3, YUC6, RTCS1, NSA1, EIN2, NHX1, CPPS4, BIGE1, RCP1, SKUS13, YUC5, and AW330564 associated with various root traits were present within or near the MQTL regions. These results could aid in QTL cloning and pyramiding in developing new maize varieties with specific root architecture for proper plant growth and development under optimum and abiotic stress conditions.

Deciphering key genomic regions controlling flag leaf size in wheat via integration of meta-QTL and in silico transcriptome assessment

BACKGROUND

Grain yield is a complex and polygenic trait influenced by the photosynthetic source-sink relationship in wheat. The top three leaves, especially the flag leaf, are considered the major sources of photo-assimilates accumulated in the grain. Determination of significant genomic regions and candidate genes affecting flag leaf size can be used in breeding for grain yield improvement.

RESULTS

With the final purpose of understanding key genomic regions for flag leaf size, a meta-analysis of 521 initial quantitative trait loci (QTLs) from 31 independent QTL mapping studies over the past decades was performed, where 333 loci eventually were refined into 64 meta-QTLs (MQTLs). The average confidence interval (CI) of these MQTLs was 5.28 times less than that of the initial QTLs. Thirty-three MQTLs overlapped the marker trait associations (MTAs) previously reported in genome-wide association studies (GWAS) for flag leaf traits in wheat. A total of 2262 candidate genes for flag leaf size, which were involved in the peroxisome, basal transcription factor, and tyrosine metabolism pathways were identified in MQTL regions by the in silico transcriptome assessment. Of these, the expression analysis of the available genes revealed that 134 genes with > 2 transcripts per million (TPM) were highly and specifically expressed in the leaf. These candidate genes could be critical to affect flag leaf size in wheat.

CONCLUSIONS

The findings will make further insight into the genetic determinants of flag leaf size and provide some reliable MQTLs and putative candidate genes for the genetic improvement of flag leaf size in wheat.

Consensus genomic regions associated with grain protein content in hexaploid and tetraploid wheat

A meta-analysis of QTLs associated with grain protein content (GPC) was conducted in hexaploid and tetraploid wheat to identify robust and stable meta-QTLs (MQTLs). For this purpose, as many as 459 GPC-related QTLs retrieved from 48 linkage-based QTL mapping studies were projected onto the newly developed wheat consensus map. The analysis resulted in the prediction of 57 MQTLs and 7 QTL hotspots located on all wheat chromosomes (except chromosomes 1D and 4D) and the average confidence interval reduced 2.71-fold in the MQTLs and QTL hotspots compared to the initial QTLs. The physical regions occupied by the MQTLs ranged from 140 bp to 224.02 Mb with an average of 15.2 Mb, whereas the physical regions occupied by QTL hotspots ranged from 1.81 Mb to 36.03 Mb with a mean of 8.82 Mb. Nineteen MQTLs and two QTL hotspots were also found to be co-localized with 45 significant SNPs identified in 16 previously published genome-wide association studies in wheat. Candidate gene (CG) investigation within some selected MQTLs led to the identification of 705 gene models which also included 96 high-confidence CGs showing significant expressions in different grain-related tissues and having probable roles in GPC regulation. These significantly expressed CGs mainly involved the genes/gene families encoding for the following proteins: aminotransferases, early nodulin 93, glutamine synthetases, invertase/pectin methylesterase inhibitors, protein BIG GRAIN 1-like, cytochrome P450, glycosyl transferases, hexokinases, small GTPases, UDP-glucuronosyl/UDP-glucosyltransferases, and EamA, SANT/Myb, GNAT, thioredoxin, phytocyanin, and homeobox domains containing proteins. Further, eight genes including GPC-B1, Glu-B1-1b, Glu-1By9, TaBiP1, GSr, TaNAC019-A, TaNAC019-D, and bZIP-TF SPA already known to be associated with GPC were also detected within some of the MQTL regions confirming the efficacy of MQTLs predicted during the current study.

Delineating meta-quantitative trait loci for anthracnose resistance in common bean (Phaseolus vulgaris L.)

Anthracnose, caused by the fungus Colletotrichum lindemuthianum, is one of the devastating disease affecting common bean production and productivity worldwide. Several quantitative trait loci (QTLs) for anthracnose resistance have been identified. In order to make use of these QTLs in common bean breeding programs, a detailed meta-QTL (MQTL) analysis has been conducted. For the MQTL analysis, 92 QTLs related to anthracnose disease reported in 18 different earlier studies involving 16 mapping populations were compiled and projected on to the consensus map. This meta-analysis led to the identification of 11 MQTLs (each involving QTLs from at least two different studies) on 06 bean chromosomes and 10 QTL hotspots each involving multiple QTLs from an individual study on 07 chromosomes. The confidence interval (CI) of the identified MQTLs was found 3.51 times lower than the CI of initial QTLs. Marker-trait associations (MTAs) reported in published genome-wide association studies (GWAS) were used to validate nine of the 11 identified MQTLs, with MQTL4.1 overlapping with as many as 40 MTAs. Functional annotation of the 11 MQTL regions revealed 1,251 genes including several R genes (such as those encoding for NBS-LRR domain-containing proteins, protein kinases, etc.) and other defense related genes. The MQTLs, QTL hotspots and the potential candidate genes identified during the present study will prove useful in common bean marker-assisted breeding programs and in basic studies involving fine mapping and cloning of genomic regions associated with anthracnose resistance in common beans.

Understanding complex genetic architecture of rice grain weight through QTL-meta analysis and candidate gene identification

Quantitative trait loci (QTL) for rice grain weight identified using bi-parental populations in various environments were found inconsistent and have a modest role in marker assisted breeding and map-based cloning programs. Thus, the identification of a consistent consensus QTL region across populations is critical to deploy in marker aided breeding programs. Using the QTL meta-analysis technique, we collated rice grain weight QTL information from numerous studies done across populations and in diverse environments to find constitutive QTL for grain weight. Using information from 114 original QTL in meta-analysis, we discovered three significant Meta-QTL (MQTL) for grain weight on chromosome 3. According to gene ontology, these three MQTL have 179 genes, 25 of which have roles in developmental functions. Amino acid sequence BLAST of these genes indicated their orthologue conservation among core cereals with similar functions. MQTL3.1 includes the OsAPX1, PDIL, SAUR, and OsASN1 genes, which are involved in grain development and have been discovered to play a key role in asparagine biosynthesis and metabolism, which is crucial for source-sink regulation. Five potential candidate genes were identified and their expression analysis indicated a significant role in early grain development. The gene sequence information retrieved from the 3 K rice genome project revealed the deletion of six bases coding for serine and alanine in the last exon of OsASN1 led to an interruption in the synthesis of α-helix of the protein, which negatively affected the asparagine biosynthesis pathway in the low grain weight genotypes. Further, the MQTL3.1 was validated using linked marker RM7197 on a set of genotypes with extreme phenotypes. MQTL that have been identified and validated in our study have significant scope in MAS breeding and map-based cloning programs for improving rice grain weight.

Consensus genomic regions associated with multiple abiotic stress tolerance in wheat and implications for wheat breeding

In wheat, a meta-analysis was performed using previously identified QTLs associated with drought stress (DS), heat stress (HS), salinity stress (SS), water-logging stress (WS), pre-harvest sprouting (PHS), and aluminium stress (AS) which predicted a total of 134 meta-QTLs (MQTLs) that involved at least 28 consistent and stable MQTLs conferring tolerance to five or all six abiotic stresses under study. Seventy-six MQTLs out of the 132 physically anchored MQTLs were also verified with genome-wide association studies. Around 43% of MQTLs had genetic and physical confidence intervals of less than 1 cM and 5 Mb, respectively. Consequently, 539 genes were identified in some selected MQTLs providing tolerance to 5 or all 6 abiotic stresses. Comparative analysis of genes underlying MQTLs with four RNA-seq based transcriptomic datasets unravelled a total of 189 differentially expressed genes which also included at least 11 most promising candidate genes common among different datasets. The promoter analysis showed that the promoters of these genes include many stress responsiveness cis-regulatory elements, such as ARE, MBS, TC-rich repeats, As-1 element, STRE, LTR, WRE3, and WUN-motif among others. Further, some MQTLs also overlapped with as many as 34 known abiotic stress tolerance genes. In addition, numerous ortho-MQTLs among the wheat, maize, and rice genomes were discovered. These findings could help with fine mapping and gene cloning, as well as marker-assisted breeding for multiple abiotic stress tolerances in wheat.

Meta-QTLs for multiple disease resistance involving three rusts in common wheat (Triticum aestivum L.)

In wheat, multiple disease resistance meta-QTLs (MDR-MQTLs) and underlying candidate genes for the three rusts were identified which may prove useful for development of resistant cultivars. Rust diseases in wheat are a major threat to global food security. Therefore, development of multiple disease-resistant cultivars (resistant to all three rusts) is a major goal in all wheat breeding programs worldwide. In the present study, meta-QTLs and candidate genes for multiple disease resistance (MDR) involving all three rusts were identified using 152 individual QTL mapping studies for resistance to leaf rust (LR), stem rust (SR), and yellow rust (YR). From these 152 studies, a total of 1,146 QTLs for resistance to three rusts were retrieved, which included 368 QTLs for LR, 291 QTLs for SR, and 487 QTLs for YR. Of these 1,146 QTLs, only 718 QTLs could be projected onto the consensus map saturated with 2, 34,619 markers. Meta-analysis of the projected QTLs resulted in the identification of 86 MQTLs, which included 71 MDR-MQTLs. Ten of these MDR-MQTLs were referred to as the 'Breeders' MQTLs'. Seventy-eight of the 86 MQTLs could also be anchored to the physical map of the wheat genome, and 54 MQTLs were validated by marker-trait associations identified during earlier genome-wide association studies. Twenty MQTLs (including 17 MDR-MQTLs) identified in the present study were co-localized with 44 known R genes. In silico expression analysis allowed identification of several differentially expressed candidate genes (DECGs) encoding proteins carrying different domains including the following: NBS-LRR, WRKY domains, F-box domains, sugar transporters, transferases, etc. The introgression of these MDR loci into high-yielding cultivars should prove useful for developing high yielding cultivars with resistance to all the three rusts.

Unravelling consensus genomic regions associated with quality traits in wheat using meta-analysis of quantitative trait loci

Meta-analysis in wheat for three major quality traits identified 110 meta-QTL (MQTL) with reduced confidence interval (CI). Five GWAS validated MQTL (viz., 1A.1, 1B.2, 3B.4, 5B.2, and 6B.2), each involving more than 20 initial QTL and reduced CI (95%) (< 2 cM), were selected for quality breeding programmes. Functional characterization including candidate gene mining and expression analysis discovered 44 high confidence candidate genes associated with quality traits. A meta-analysis of quantitative trait loci (QTL) associated with dough rheology properties, nutritional traits, and processing quality traits was conducted in wheat. For this purpose, as many as 2458 QTL were collected from 50 interval mapping studies published during 2013-2020. Of the total QTL, 1126 QTL were projected onto the consensus map saturated with 249,603 markers which led to the identification of 110 meta-QTL (MQTL). These MQTL exhibited an 18.84-fold reduction in the average CI compared to the average CI of the initial QTL (ranging from 14.87 to 95.55 cM with an average of 40.35 cM). Of the 110, 108 MQTL were physically anchored to the wheat reference genome, including 51 MQTL verified with marker-trait associations (MTAs) reported from earlier genome-wide association studies. Candidate gene (CG) mining allowed the identification of 2533 unique gene models from the MQTL regions. In-silico expression analysis discovered 439 differentially expressed gene models with > 2 transcripts per million expressions in grains and related tissues, which also included 44 high-confidence CGs involved in the various cellular and biochemical processes related to quality traits. Nine functionally characterized wheat genes associated with grain protein content, high-molecular-weight glutenin, and starch synthase enzymes were also found to be co-localized with some of the MQTL. Synteny analysis between wheat and rice MQTL regions identified 23 wheat MQTL syntenic to 16 rice MQTL associated with quality traits. Furthermore, 64 wheat orthologues of 30 known rice genes were detected in 44 MQTL regions. Markers flanking the MQTL identified in the present study can be used for marker-assisted breeding and as fixed effects in the genomic selection models for improving the prediction accuracy during quality breeding. Wheat orthologues of rice genes and other CGs available from MQTLs can be promising targets for further functional validation and to better understand the molecular mechanism underlying the quality traits in wheat.

Meta-QTLs, ortho-MetaQTLs and candidate genes for grain Fe and Zn contents in wheat (Triticum aestivum L.)

Majority of cereals are deficient in essential micronutrients including grain iron (GFe) and grain zinc (GZn), which are therefore the subject of research involving biofortification. In the present study, 11 meta-QTLs (MQTLs) including nine novel MQTLs for GFe and GZn contents were identified in wheat. Eight of these 11 MQTLs controlled both GFe and GZn. The confidence intervals of the MQTLs were narrower (0.51-15.75 cM) relative to those of the corresponding QTLs (0.6 to 55.1 cM). Two ortho-MQTLs involving three cereals (wheat, rice and maize) were also identified. Results of MQTLs were also compared with the results of earlier genome wide association studies (GWAS). As many as 101 candidate genes (CGs) underlying MQTLs were also identified. Twelve of these CGs were prioritized; these CGs encoded proteins with important domains (zinc finger, RING/FYVE/PHD type, flavin adenine dinucleotide linked oxidase, etc.) that are involved in metal ion binding, heme binding, iron binding, etc. qRT-PCR analysis was conducted for four of these 12 prioritized CGs using genotypes which have differed for GFe and GZn. Significant differential expression in these genotypes was observed at 14 and 28 days after anthesis. The MQTLs/CGs identified in the present study may be utilized in marker-assisted selection (MAS) for improvement of GFe/GZn contents and also for understanding the molecular basis of GFe/GZn homeostasis in wheat.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12298-022-01149-9.

Meta-QTLs, ortho-meta-QTLs and candidate genes for grain yield and associated traits in wheat (Triticum aestivum L.)

In wheat, 2852 major QTLs of 8998 QTLs available for yield and related traits were used for meta-analysis; 141 meta-QTLs were identified, which included 13 breeder's MQTLs and 24 ortho-MQTLs; 1202 candidate genes and 50 homologues of genes for yield from other cereals were also identified. Meta-QTL analysis was conducted using 2852 of the 8998 known QTLs, retrieved from 230 reports published during 1999-2020 (including 19 studies on tetraploid wheat) for grain yield (GY) and the following ten component traits: (i) grain weight (GWei), (ii) grain morphology-related traits (GMRTs), (iii) grain number (GN), (iv) spikes-related traits (SRTs), (v) plant height (PH), (vi) tiller number (TN), (vii) harvest index (HI), (viii) biomass yield (BY), (ix) days to heading/flowering and maturity (DTH/F/M), and (x) grain filling duration (GFD). The study resulted in the identification of 141 meta-QTLs (MQTLs), with an average confidence interval (CI) of 1.4 cM as against a CI of > 12.1 cM (8.8 fold reduction) in the QTLs that were used. The corresponding physical length of CI ranged from 0.01 Mb to 661.9 Mb (mean, 31.5 Mb). Seventy-seven (77) of these 141 MQTLs overlapped marker-trait associations (MTAs) reported in genome-wide association studies. Also, 63 MQTLs (each based on at least 10 QTLs) were considered stable and robust, with 13 MQTLs described as breeder's MQTLs (selected based on small CI, large LOD, and high level of phenotypic variation explained). Thirty-five yield-related genes from rice, barley, and maize were also utilized to identify 50 wheat homologues in MQTLs. Further, the use of synteny and collinearity allowed the identification of 24 ortho-MQTLs which were common among the wheat, barley, rice, and maize. The results of the present study should prove useful for wheat breeding and future basic research in cereals including wheat, barley, rice, and maize.

Meta-analysis reveals consensus genomic regions associated with multiple disease resistance in wheat (Triticum aestivum L.)

In wheat, meta-QTLs (MQTLs) and candidate genes (CGs) were identified for multiple disease resistance (MDR). For this purpose, information was collected from 58 studies for mapping QTLs for resistance to one or more of the five diseases. As many as 493 QTLs were available from these studies, which were distributed in five diseases as follows: septoria tritici blotch (STB) 126 QTLs; septoria nodorum blotch (SNB), 103 QTLs; fusarium head blight (FHB), 184 QTLs; karnal bunt (KB), 66 QTLs; and loose smut (LS), 14 QTLs. Of these 493 QTLs, only 291 QTLs could be projected onto a consensus genetic map, giving 63 MQTLs. The CI of the MQTLs ranged from 0.04 to 15.31 cM with an average of 3.09 cM per MQTL. This is a ~ 4.39 fold reduction from the CI of QTLs, which ranged from 0 to 197.6 cM, with a mean of 13.57 cM. Of 63 MQTLs, 60 were anchored to the reference physical map of wheat (the physical interval of these MQTLs ranged from 0.30 to 726.01 Mb with an average of 74.09 Mb). Thirty-eight (38) of these MQTLs were verified using marker-trait associations (MTAs) derived from genome-wide association studies. As many as 874 CGs were also identified which were further investigated for differential expression using data from five transcriptome studies, resulting in 194 differentially expressed candidate genes (DECGs). Among the DECGs, 85 genes had functions previously reported to be associated with disease resistance. These results should prove useful for fine mapping and cloning of MDR genes and marker-assisted breeding.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-022-01282-z.

Meta-QTLs, ortho-MQTLs and candidate genes for the traits contributing to salinity stress tolerance in common wheat (Triticum aestivum L.)

A meta-analysis of QTLs associated with the traits contributing to salinity tolerance was undertaken in wheat to detect consensus and robust meta-QTLs (MQTLs) using 844 known QTLs retrieved from 26 earlier studies. A consensus map with a total length of 4621.56 cM including 7710 markers was constructed using 21 individual linkage maps and three previously published integrated genetic maps. Out of 844 QTLs, 571 QTLs were projected on the consensus map which gave origin to 100 MQTLs. Interestingly, 49 MQTLs were co-located with marker-trait associations reported in wheat genome-wide association studies for the traits contributing to salinity stress tolerance. Five potential MQTLs associated with the major salinity-responsive traits were also identified to be utilized in the breeding programme. In the resulted MQTLs, the average confidence interval (CI, 3.58 cM) was reduced up to 4.16 folds compared to the mean CI of the initial QTLs. Furthermore, as many as 617 gene models including 81 most likely candidate genes (CGs) were identified in the high confidence MQTL regions. These most likely CGs encoded proteins mainly belonging to the following families: B-box-type zinc finger, cytochrome P450 protein, pentatricopeptide repeat, phospholipid/glycerol acyltransferase, F-box protein, small auxin-up RNA, UDP-glucosyltransferase, glutathione S-transferase protein, etc. In addition, ortho-MQTL analysis based on synteny among wheat, rice and barley was also performed which permitted the identification of six ortho-MQTLs among these three cereals. This meta-analysis defines a genome-wide landscape on the most stable and consistent loci associated with reliable molecular markers and candidate genes for salinity tolerance in wheat.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s12298-021-01112-0.