MetaQTL: a package of new computational methods for the meta-analysis of QTL mapping experiments

Dissection of QTLs underlying the genetic basis of drought resistance in wheat: a meta-analysis

Wheat (Triticum aestivum L.) is one of the most important cereal crops, with its grain serving as a predominant staple food source on a global scale. However, there are many biotic and abiotic stresses challenging the stability of wheat production. Among the abiotic stresses, drought is recognized as a significant stress and poses a substantial threat to food production and quality throughout the world. Raising drought tolerance of wheat varieties through genetic regulation is therefore considered as one of the most effective ways to combat the challenges caused by drought stress. Meta-QTL analysis has demonstrated its effectiveness in identifying consensus QTL regions in wheat drought resistance in numerous instances. In this study, we present a comprehensive meta-analysis aimed at unraveling the drought tolerance genetic basis associated with agronomic traits in bread wheat. Extracting data from 34 previously published studies, we aggregated a corpus of 1291 Quantitative Trait Loci (QTL) pertinent to wheat drought tolerance. Then, the translation of the consensus genetic map yielded a comprehensive compendium of 49 distinct MQTLs, each associated with diverse agronomic traits. Prominently featured among the MQTLs were MQTLs 1.1, 1.7, 1.8 (1D), 4.1 (4A), 4.6 (4D), 5.2 (5B), 6.6 (6B), and 7.2 (7B), distinguished as pivotal MQTLs offering significant potential for application in marker-assisted breeding endeavors. Altogether, a total of 66 putative candidate genes (CGs)-related drought tolerance were identified. This work illustrates a translational research approach in transferring information from published mapping studies to genomic regions hosting major QTLs governing key agronomical traits in wheat.

Meta-QTL mapping for wheat thousand kernel weight

Wheat domestication and subsequent genetic improvement have yielded cultivated species with larger seeds compared to wild ancestors. Increasing thousand kernel weight (TKW) remains a crucial goal in many wheat breeding programs. To identify genomic regions influencing TKW across diverse genetic populations, we performed a comprehensive meta-analysis of quantitative trait loci (MQTL), integrating 993 initial QTL from 120 independent mapping studies over recent decades. We refined 242 loci into 66 MQTL, with an average confidence interval (CI) 3.06 times smaller than that of the original QTL. In these 66 MQTL regions, a total of 4,913 candidate genes related to TKW were identified, involved in ubiquitination, phytohormones, G-proteins, photosynthesis, and microRNAs. Expression analysis of the candidate genes showed that 95 were specific to grain and might potentially affect TKW at different seed development stages. These findings enhance our understanding of the genetic factors associated with TKW in wheat, providing reliable MQTL and potential candidate genes for genetic improvement of this trait.

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.

Integrating RTM-GWAS and meta‑QTL data revealed genomic regions and candidate genes associated with the first fruit branch node and its height in upland cotton

KEY MESSAGE

Two genomic regions associated with FFBN and HFFBN and a potential regulatory gene (GhE6) of HFFBN were identified through the integration of RTM-GWAS and meta‑QTL analyses. Abstract The first fruit branch node (FFBN) and the height of the first fruit branch node (HFFBN) are two important traits that are related to plant architecture and early maturation in upland cotton. Several studies have been conducted to elucidate the genetic basis of these traits in cotton using biparental and natural populations. In this study, by using 9,244 SNP linkage disequilibrium block (SNPLDB) loci from 315 upland cotton accessions, we carried out restricted two-stage multilocus and multiallele genome-wide association studies (RTM-GWASs) and identified promising haplotypes/alleles of the four stable and true major SNPLDB loci that were significantly associated with FFBN and HFFBN. Additionally, a meta-quantitative trait locus (MQTL) analysis was conducted on 274 original QTLs that were reported in 27 studies, and 40 MQTLs associated with FFBN and HFFBN were identified. Through the integration of the RTM-GWAS and meta‑QTL analyses, two stable and true major SNPLDBs (LDB_5_15144433 and LDB_16_37952328) that were distributed in the two MQTLs were identified. Ultimately, 142 genes in the two genomic regions were annotated, and three candidate genes associated with FFBN and HFFBN were identified in the genomic region (A05:14.64-15.64 Mb) via RNA-Seq and qRT‒PCR. The results of virus-induced gene silencing (VIGS) experiments indicated that GhE6 was a key gene related to HFFBN and that GhDRM1 and GhGES were important genes associated with early flowering in upland cotton. These findings will aid in the future identification of molecular markers and genetic resources for developing elite early-maturing cultivars with ideal plant characteristics.

Refinement of rice blast disease resistance QTLs and gene networks through meta-QTL analysis

Rice blast disease is the most devastating disease constraining crop productivity. Vertical resistance to blast disease is widely studied despite its instability. Clusters of genes or QTLs conferring blast resistance that offer durable horizontal resistance are important in resistance breeding. In this study, we aimed to refine the reported QTLs and identify stable meta-QTLs (MQTLs) associated with rice blast resistance. A total of 435 QTLs were used to project 71 MQTLs across all the rice chromosomes. As many as 199 putative rice blast resistance genes were identified within 53 MQTL regions. The genes included 48 characterized resistance gene analogs and related proteins, such as NBS-LRR type, LRR receptor-like kinase, NB-ARC domain, pathogenesis-related TF/ERF domain, elicitor-induced defense and proteins involved in defense signaling. MQTL regions with clusters of RGA were also identified. Fifteen highly significant MQTLs included 29 candidate genes and genes characterized for blast resistance, such as Piz, Nbs-Pi9, pi55-1, pi55-2, Pi3/Pi5-1, Pi3/Pi5-2, Pikh, Pi54, Pik/Pikm/Pikp, Pb1 and Pb2. Furthermore, the candidate genes (42) were associated with differential expression (in silico) in compatible and incompatible reactions upon disease infection. Moreover, nearly half of the genes within the MQTL regions were orthologous to those in O. sativa indica, Z. mays and A. thaliana, which confirmed their significance. The peak markers within three significant MQTLs differentiated blast-resistant and susceptible lines and serve as potential surrogates for the selection of blast-resistant lines. These MQTLs are potential candidates for durable and broad-spectrum rice blast resistance and could be utilized in blast resistance breeding.

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.

Unraveling the genetic basis of quantitative resistance to diseases in tomato: a meta-QTL analysis and mining of transcript profiles

KEY MESSAGE

Whole-genome QTL mining and meta-analysis in tomato for resistance to bacterial and fungal diseases identified 73 meta-QTL regions with significantly refined/reduced confidence intervals. Tomato production is affected by a range of biotic stressors, causing yield losses and quality reductions. While sources of genetic resistance to many tomato diseases have been identified and characterized, stability of the resistance genes or quantitative trait loci (QTLs) across the resources has not been determined. Here, we examined 491 QTLs previously reported for resistance to tomato diseases in 40 independent studies and 54 unique mapping populations. We identified 29 meta-QTLs (MQTLs) for resistance to bacterial pathogens and 44 MQTLs for resistance to fungal pathogens, and were able to reduce the average confidence interval (CI) of the QTLs by 4.1-fold and 6.7-fold, respectively, compared to the average CI of the original QTLs. The corresponding physical length of the CIs of MQTLs ranged from 56 kb to 6.37 Mb, with a median of 921 kb, of which 27% had a CI lower than 500 kb and 53% had a CI lower than 1 Mb. Comparison of defense responses between tomato and Arabidopsis highlighted 73 orthologous genes in the MQTL regions, which were putatively determined to be involved in defense against bacterial and fungal diseases. Intriguingly, multiple genes were identified in some MQTL regions that are implicated in plant defense responses, including PR-P2, NDR1, PDF1.2, Pip1, SNI1, PTI5, NSL1, DND1, CAD1, SlACO, DAD1, SlPAL, Ph-3, EDS5/SID1, CHI-B/PR-3, Ph-5, ETR1, WRKY29, and WRKY25. Further, we identified a number of candidate resistance genes in the MQTL regions that can be useful for both marker/gene-assisted breeding as well as cloning and genetic transformation.

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.

Genome-wide screening of meta-QTL and candidate genes controlling yield and yield-related traits in barley (Hordeum vulgare L.)

Increasing yield is an important goal of barley breeding. In this study, 54 papers published from 2001-2022 on QTL mapping for yield and yield-related traits in barley were collected, which contained 1080 QTLs mapped to the barley high-density consensus map for QTL meta-analysis. These initial QTLs were integrated into 85 meta-QTLs (MQTL) with a mean confidence interval (CI) of 2.76 cM, which was 7.86-fold narrower than the CI of the initial QTL. Among these 85 MQTLs, 68 MQTLs were validated in GWAS studies, and 25 breeder's MQTLs were screened from them. Seventeen barley orthologs of yield-related genes in rice and maize were identified within the hcMQTL region based on comparative genomics strategy and were presumed to be reliable candidates for controlling yield-related traits. The results of this study provide useful information for molecular marker-assisted breeding and candidate gene mining of yield-related traits in barley.

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.

Meta-analysis of identified genomic regions and candidate genes underlying salinity tolerance in rice (Oryza sativa L.)

Rice output has grown globally, yet abiotic factors are still a key cause for worry. Salinity stress seems to have the more impact on crop production out of all abiotic stresses. Currently one of the most significant challenges in paddy breeding for salinity tolerance with the help of QTLs, is to determine the QTLs having the best chance of improving salinity tolerance with the least amount of background noise from the tolerant parent. Minimizing the size of the QTL confidence interval (CI) is essential in order to primarily include the genes responsible for salinity stress tolerance. By considering that, a genome-wide meta-QTL analysis on 768 QTLs from 35 rice populations published from 2001 to 2022 was conducted to identify consensus regions and the candidate genes underlying those regions responsible for the salinity tolerance, as it reduces the confidence interval (CI) to many folds from the initial QTL studies. In the present investigation, a total of 65 MQTLs were extracted with an average CI reduced from 17.35 to 1.66 cM including the smallest of 0.01 cM. Identification of the MQTLs for individual traits and then classifying the target traits into correlated morphological, physiological and biochemical aspects, resulted in more efficient interpretation of the salinity tolerance, identifying the candidate genes and to understand the salinity tolerance mechanism as a whole. The results of this study have a huge potential to improve the rice genotypes for salinity tolerance with the help of MAS and MABC.

Uncovering the Genomic Regions Associated with Yield Maintenance in Rice Under Drought Stress Using an Integrated Meta-Analysis Approach

The complex trait of yield is controlled by several quantitative trait loci (QTLs). Given the global water deficit issue, the development of rice varieties suitable for non-flooded cultivation holds significant importance in breeding programs. The powerful approach of Meta-QTL (MQTL) analysis can be used for the genetic dissection of complicated quantitative traits. In the current study, a comprehensive MQTL analysis was conducted to identify consistent QTL regions associated with drought tolerance and yield-related traits under water deficit conditions in rice. In total, 1087 QTLs from 134 rice populations, published between 2000 to 2021, were utilized in the analysis. Distinct MQTL analysis of the relevant traits resulted in the identification of 213 stable MQTLs. The confidence interval (CI) for the detected MQTLs was between 0.12 and 19.7 cM. The average CI of the identified MQTLs (4.68 cM) was 2.74 times narrower compared to the average CI of the initial QTLs. Interestingly, 63 MQTLs coincided with SNP peak positions detected by genome-wide association studies for yield and drought tolerance-associated traits under water deficit conditions in rice. Considering the genes located both in the QTL-overview peaks and the SNP peak positions, 19 novel candidate genes were introduced, which are associated with drought response index, plant height, panicle number, biomass, and grain yield. Moreover, an inclusive MQTL analysis was performed on all the traits to obtain "Breeding MQTLs". This analysis resulted in the identification of 96 MQTLs with a CI ranging from 0.01 to 9.0 cM. The mean CI of the obtained MQTLs (2.33 cM) was 4.66 times less than the mean CI of the original QTLs. Thirteen MQTLs fulfilling the criteria of having more than 10 initial QTLs, CI < 1 cM, and an average phenotypic variance explained greater than 10%, were designated as "Breeding MQTLs". These findings hold promise for assisting breeders in enhancing rice yield under drought stress conditions.

Detection of consensus genomic regions and candidate genes for quality traits in barley using QTL meta-analysis

Improving barley grain quality is a major goal in barley breeding. In this study, a total of 35 papers focusing on quantitative trait loci (QTLs) mapping for barley quality traits published since 2000 were collected. Among the 454 QTLs identified in these studies, 349 of them were mapped onto high-density consensus maps, which were used for QTL meta-analysis. Through QTL meta-analysis, the initial QTLs were integrated into 41 meta-QTLs (MQTLs) with an average confidence interval (CI) of 1. 66 cM, which is 88.9% narrower than that of the initial QTLs. Among the 41 identified MQTLs, 25 were subsequently validated in publications using genome-wide association study (GWAS). From these 25 validated MQTLs, ten breeder's MQTLs were selected. Synteny analysis comparing barley and wheat MQTLs revealed orthologous relationships between eight breeder's MQTLs and 45 wheat MQTLs. Additionally, 17 barley homologs associated with rice quality traits were identified within the regions of the breeder's MQTLs through comparative analysis. The findings of this study provide valuable insights for molecular marker-assisted breeding and the identification of candidate genes related to quality traits in barley.

Genetic Dissection of Panicle Morphology Traits in Super High-Yield Hybrid Rice Chaoyou 1000

The morphological characteristics of the rice panicle play a pivotal role in influencing yield. In our research, we employed F2 and F2:3 populations derived from the high-yielding hybrid rice variety Chaoyou 1000. We screened 123 pairs of molecular markers, which were available, to construct the genetic linkage map. Subsequently, we assessed the panicle morphology traits of F2 populations in Lingshui County, Hainan Province, in 2017, and F2:3 populations in Hangzhou City, Zhejiang Province, in 2018. These two locations represent two types of ecology. Hangzhou's climate is characterized by high temperatures and humidity, while Lingshui's climate is characterized by a tropical monsoon climate. In total, 33 QTLs were identified, with eight of these being newly discovered, and two of them were consistently detected in two distinct environments. We identified fourteen QTL-by-environment interactions (QEs), which collectively explained 4.93% to 59.95% of the phenotypic variation. While most of the detected QTLs are consistent with the results of previous tests, the novel-detected QTLs will lay the foundation for rice yield increase and molecular breeding.

Meta-QTL and ortho analysis unravels the genetic architecture and key candidate genes for cold tolerance at seedling stage in rice

Rice, a critical cereal crop, grapples with productivity challenges due to its inherent sensitivity to low temperatures, primarily during the seedling and booting stages. Recognizing the polygenic complexity of cold stress signaling in rice, a meta-analysis was undertaken, focusing on 20 physiological traits integral to cold tolerance. This initiative allowed the consolidation of genetic data from 242 QTLs into 58 meta-QTLs, thereby significantly constricting the genetic and physical intervals, with 84% of meta-QTLs (MQTLs) being reduced to less than 2 Mb. The list of 10,505 genes within these MQTLs, was further refined utilizing expression datasets to pinpoint 46 pivotal genes exhibiting noteworthy differential regulation during cold stress. The study underscored the presence of several TFs such as WRKY, NAC, CBF/DREB, MYB, and bHLH, known for their roles in cold stress response. Further, ortho-analysis involving maize, barley, and Arabidopsis identified OsWRKY71, among others, as a prospective candidate for enhancing cold tolerance in diverse crop plants. In conclusion, our study delineates the intricate genetic architecture underpinning cold tolerance in rice and propounds significant candidate genes, offering crucial insights for further research and breeding strategies focused on fortifying crops against cold stress, thereby bolstering global food resilience.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12298-024-01412-1.

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.

High confidence QTLs and key genes identified using Meta-QTL analysis for enhancing heat tolerance in chickpea (Cicer arietinum L.)

The rising global temperatures seriously threaten sustainable crop production, particularly the productivity and production of heat-sensitive crops like chickpeas. Multiple QTLs have been identified to enhance the heat stress tolerance in chickpeas, but their successful use in breeding programs remains limited. Towards this direction, we constructed a high-density genetic map spanning 2233.5 cM with 1069 markers. Using 138 QTLs reported earlier, we identified six Meta-QTL regions for heat tolerance whose confidence interval was reduced by 2.7-folds compared to the reported QTLs. Meta-QTLs identified on CaLG01 and CaLG06 harbor QTLs for important traits, including days to 50% flowering, days to maturity, days to flower initiation, days to pod initiation, number of filled pods, visual score, seed yield per plant, biological yield per plant, chlorophyll content, and harvest index. In addition, key genes identified in Meta-QTL regions like Pollen receptor-like kinase 3 (CaPRK3), Flowering-promoting factor 1 (CaFPF1), Flowering Locus C (CaFLC), Heat stress transcription factor A-5 (CaHsfsA5), and Pollen-specific leucine-rich repeat extensins (CaLRXs) play an important role in regulating the flowering time, pollen germination, and growth. The consensus genomic regions, and the key genes reported in this study can be used in genomics-assisted breeding for enhancing heat tolerance and developing heat-resilient chickpea cultivars.

A meta-QTL analysis highlights genomic hotspots associated with phosphorus use efficiency in rice (Oryza sativa L.)

Phosphorus use efficiency (PUE) is a complex trait, governed by many minor quantitative trait loci (QTLs) with small effects. Advances in molecular marker technology have led to the identification of QTLs underlying PUE. However, their practical use in breeding programs remains challenging due to the unstable effects in different genetic backgrounds and environments, interaction with soil status, and linkage drag. Here, we compiled PUE QTL information from 16 independent studies. A total of 192 QTLs were subjected to meta-QTL (MQTL) analysis and were projected into a high-density SNP consensus map. A total of 60 MQTLs, with significantly reduced number of initial QTLs and confidence intervals (CI), were identified across the rice genome. Candidate gene (CG) mining was carried out for the 38 MQTLs supported by multiple QTLs from at least two independent studies. Genes related to amino and organic acid transport and auxin response were found to be abundant in the MQTLs linked to PUE. CGs were cross validated using a root transcriptome database (RiceXPro) and haplotype analysis. This led to the identification of the eight CGs (OsARF8, OsSPX-MFS3, OsRING141, OsMIOX, HsfC2b, OsFER2, OsWRKY64, and OsYUCCA11) modulating PUE. Potential donors for superior PUE CG haplotypes were identified through haplotype analysis. The distribution of superior haplotypes varied among subspecies being mostly found in indica but were largely scarce in japonica. Our study offers an insight on the complex genetic networks that modulate PUE in rice. The MQTLs, CGs, and superior CG haplotypes identified in our study are useful in the combination of beneficial alleles for PUE in rice.

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.

Comprehensive meta-analysis of QTL and gene expression studies identify candidate genes associated with Aspergillus flavus resistance in maize

Aflatoxin (AF) contamination, caused by Aspergillus flavus, compromises the food safety and marketability of commodities, such as maize, cotton, peanuts, and tree nuts. Multigenic inheritance of AF resistance impedes conventional introgression of resistance traits into high-yielding commercial maize varieties. Several AF resistance-associated quantitative trait loci (QTLs) and markers have been reported from multiple biparental mapping and genome-wide association studies (GWAS) in maize. However, QTLs with large confidence intervals (CI) explaining inconsistent phenotypic variance limit their use in marker-assisted selection. Meta-analysis of published QTLs can identify significant meta-QTLs (MQTLs) with a narrower CI for reliable identification of genes and linked markers for AF resistance. Using 276 out of 356 reported QTLs controlling resistance to A. flavus infection and AF contamination in maize, we identified 58 MQTLs on all 10 chromosomes with a 66.5% reduction in the average CI. Similarly, a meta-analysis of maize genes differentially expressed in response to (a)biotic stresses from the to-date published literature identified 591 genes putatively responding to only A. flavus infection, of which 14 were significantly differentially expressed (-1.0 ≤ Log2Fc ≥ 1.0; p ≤ 0.05). Eight MQTLs were validated by their colocalization with 14 A. flavus resistance-associated SNPs identified from GWAS in maize. A total of 15 genes were physically close between the MQTL intervals and SNPs. Assessment of 12 MQTL-linked SSR markers identified three markers that could discriminate 14 and eight cultivars with resistance and susceptible responses, respectively. A comprehensive meta-analysis of QTLs and differentially expressed genes led to the identification of genes and makers for their potential application in marker-assisted breeding of A. flavus-resistant maize varieties.

Meta-analysis of seed weight QTLome using a consensus and highly dense genetic map in Brassica napus L

KEY MESSAGE

We report here the discovery of high-confidence MQTL regions and of putative candidate genes associated with seed weight in B. napus using a highly dense consensus genetic map and by comparing various large-scale multiomics datasets. Seed weight (SW) is a direct determinant of seed yield in Brassica napus and is controlled by many loci. To unravel the main genomic regions associated with this complex trait, we used 13 available genetic maps to construct a consensus and highly dense map, comprising 40,401 polymorphic markers and 9191 genetic bins, harboring a cumulative length of 3047.8 cM. Then, we performed a meta-analysis using 639 projected SW quantitative trait loci (QTLs) obtained from studies conducted since 1999, enabling the identification of 57 meta-QTLS (MQTLs). The confidence intervals of our MQTLs were 9.8 and 4.3 times lower than the average CIs of the original QTLs for the A and C subgenomes, respectively, resulting in the detection of some key genes and several putative novel candidate genes associated with SW. By comparing the genes identified in MQTL intervals with multiomics datasets and coexpression analyses of common genes, we defined a more reliable and shorter list of putative candidate genes potentially involved in the regulation of seed maturation and SW. As an example, we provide a list of promising genes with high expression levels in seeds and embryos (e.g., BnaA03g04230D, BnaC03g08840D, BnaA10g29580D and BnaA03g27410D) that can be more finely studied through functional genetics experiments or that may be useful for MQTL-assisted breeding for SW. The high-density genetic consensus map and the single nucleotide polymorphism (SNP) physical map generated from the latest B. napus cv. Darmor-bzh v10 assembly will be a valuable resource for further mapping and map-based cloning of other important traits.

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-Analysis of Rose Rosette Disease-Resistant Quantitative Trait Loci and a Search for Candidate Genes

Rose rosette disease (RRD), caused by the rose rosette emaravirus (RRV), is a major viral disease in roses (Rosa sp.) that threatens the rose industry. Recent studies have revealed quantitative trait loci (QTL) for reduced susceptibility to RRD in the linkage groups (LGs) 1, 5, 6, and 7 in tetraploid populations and the LGs 1, 3, 5, and 6 in diploid populations. In this study, we seek to better localize and understand the relationship between QTL identified in both diploid and tetraploid populations. We do so by remapping the populations found in these studies and performing a meta-analysis. This analysis reveals that the peaks and intervals for QTL using diploid and tetraploid populations co-localized on LG 1, suggesting that these are the same QTL. The same was seen on LG 3. Three meta-QTL were identified on LG 5, and two were discovered on LG 6. The meta-QTL on LG 1, MetaRRD1.1, had a confidence interval (CI) of 10.53 cM. On LG 3, MetaRRD3.1 had a CI of 5.94 cM. MetaRRD5.1 had a CI of 17.37 cM, MetaRRD5.2 had a CI of 4.33 cM, and MetaRRD5.3 had a CI of 21.95 cM. For LG 6, MetaRRD6.1 and MetaRRD6.2 had CIs of 9.81 and 8.81 cM, respectively. The analysis also led to the identification of potential disease resistance genes, with a primary interest in genes localized in meta-QTL intervals on LG 5 as this LG was found to explain the greatest proportion of phenotypic variance for RRD resistance. The results from this study may be used in the design of more robust marker-based selection tools to track and use a given QTL in a plant breeding context.

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 association study to identify genomic loci associated with early vigor in bread wheat under simulated water deficit complemented with quantitative trait loci meta-analysis

A genome-wide association study (GWAS) was used to identify associated loci with early vigor under simulated water deficit and grain yield under field drought in a diverse collection of Iranian bread wheat landraces. In addition, a meta-quantitative trait loci (MQTL) analysis was used to further expand our approach by retrieving already published quantitative trait loci (QTL) from recombinant inbred lines, double haploids, back-crosses, and F2 mapping populations. In the current study, around 16%, 14%, and 16% of SNPs were in significant linkage disequilibrium (LD) in the A, B, and D genomes, respectively, and varied between 5.44% (4A) and 21.85% (6A). Three main subgroups were identified among the landraces with different degrees of admixture, and population structure was further explored through principal component analysis. Our GWAS identified 54 marker-trait associations (MTAs) that were located across the wheat genome but with the highest number found in the B sub-genome. The gene ontology (GO) analysis of MTAs revealed that around 75% were located within or closed to protein-coding genes. In the MQTL analysis, 23 MQTLs, from a total of 215 QTLs, were identified and successfully projected onto the reference map. MQT-YLD4, MQT-YLD9, MQT-YLD13, MQT-YLD17, MQT-YLD18, MQT-YLD19, and MQTL-RL1 contributed to the highest number of projected QTLs and were therefore regarded as the most reliable and stable QTLs under water deficit conditions. These MQTLs greatly facilitate the identification of putative candidate genes underlying at each MQTL interval due to the reduced confidence of intervals associated with MQTLs. These findings provide important information on the genetic basis of early vigor traits and grain yield under water deficit conditions and set the foundation for future investigations into adaptation to water deficit in bread wheat.

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.

Meta-QTL Analysis for Yield Components in Common Bean (Phaseolus vulgaris L.)

Common bean is one of the most important legumes produced and consumed worldwide because it is a highly valuable food for the human diet. However, its production is mainly carried out by small farmers, who obtain average grain yields below the potential yield of the species. In this sense, numerous mapping studies have been conducted to identify quantitative trait loci (QTL) associated with yield components in common bean. Meta-QTL (MQTL) analysis is a useful approach to combine data sets and for creating consensus positions for the QTL detected in independent studies. Consequently, the objective of this study was to perform a MQTL analysis to identify the most reliable and stable genomic regions associated with yield-related traits of common bean. A total of 667 QTL associated with yield-related traits reported in 21 different studies were collected. A total of 42 MQTL associated with yield-related traits were identified, in which the average confidence interval (CI) of the MQTL was 3.41 times lower than the CIs of the original QTL. Most of the MQTL (28) identified in this study contain QTL associated with yield and phenological traits; therefore, these MQTL can be useful in common bean breeding programs. Finally, a total of 18 candidate genes were identified and associated with grain yield within these MQTL, with functions related to ubiquitin ligase complex, response to auxin, and translation elongation factor activity.

Identification of quantitative trait loci (QTL) and meta-QTL analysis for kernel size-related traits in wheat (Triticum aestivum L.)

BACKGROUND

Kernel size-related traits, including kernel length (KL), kernel width (KW), kernel diameter ratio (KDR) and kernel thickness (KT), are critical determinants for wheat kernel weight and yield and highly governed by a type of quantitative genetic basis. Genome-wide identification of major and stable quantitative trait loci (QTLs) and functional genes are urgently required for genetic improvement in wheat kernel yield. A hexaploid wheat population consisting of 120 recombinant inbred lines was developed to identify QTLs for kernel size-related traits under different water environments. The meta-analysis and transcriptome evaluation were further integrated to identify major genomic regions and putative candidate genes.

RESULTS

The analysis of variance (ANOVA) revealed more significant genotypic effects for kernel size-related traits, indicating the moderate to high heritability of 0.61-0.89. Thirty-two QTLs for kernel size-related traits were identified, explaining 3.06%-14.2% of the phenotypic variation. Eleven stable QTLs were detected in more than three water environments. The 1103 original QTLs from the 34 previous studies and the present study were employed for the MQTL analysis and refined into 58 MQTLs. The average confidence interval of the MQTLs was 3.26-fold less than that of the original QTLs. The 1864 putative candidate genes were mined within the regions of 12 core MQTLs, where 70 candidate genes were highly expressed in spikes and kernels by comprehensive analysis of wheat transcriptome data. They were involved in various metabolic pathways, such as carbon fixation in photosynthetic organisms, carbon metabolism, mRNA surveillance pathway, RNA transport and biosynthesis of secondary metabolites.

CONCLUSIONS

Major genomic regions and putative candidate genes for kernel size-related traits in wheat have been revealed by an integrative strategy with QTL linkage mapping, meta-analysis and transcriptomic assessment. The findings provide a novel insight into understanding the genetic determinants of kernel size-related traits and will be useful for the marker-assisted selection of high yield in wheat breeding.

Genetic structure and molecular mechanism underlying the stalk lodging traits in maize (Zea mays L.)

Stalk lodging seriously affects yield and quality of crops, and it can be caused by several factors, such as environments, developmental stages, and internal chemical components of plant stalks. Breeding of stalk lodging-resistant varieties is thus an important task for maize breeders. To better understand the genetic basis underlying stalk lodging resistance, several methods such as quantitative trait locus (QTL) mapping and genome-wide association study (GWAS) have been used to mine potential gene resources. Based on different types of genetic populations and mapping methods, many significant loci associated with stalk lodging resistance have been identified so far. However, few work has been performed to compare and integrate these reported genetic loci. In this study, we first collected hundreds of QTLs and quantitative trait nucleotides (QTNs) related to stalk lodging traits in maize. Then we mapped and integrated the QTLs and QTNs in maize genome to identify overlapped hotspot regions. Based on the genomic confidence intervals harboring these overlapped hotspot regions, we predicted candidate genes related to stalk lodging traits. Meanwhile, we mapped reported genes to these hotspot regions. Finally, we constructed molecular regulatory networks underlying stalk lodging resistance in maize. Collectively, this study provides not only useful genetic loci for deeply exploring molecular mechanisms of stalk lodging resistance traits, but also potential candidate genes and targeted strategies for improving stalk lodging resistance to increase crop yields in future.

A meta-quantitative trait loci analysis identified consensus genomic regions and candidate genes associated with grain yield in rice

Improving grain yield potential in rice is an important step toward addressing global food security challenges. The meta-QTL analysis offers stable and robust QTLs irrespective of the genetic background of mapping populations and phenotype environment and effectively narrows confidence intervals (CI) for candidate gene (CG) mining and marker-assisted selection improvement. To achieve these aims, a comprehensive bibliographic search for grain yield traits (spikelet fertility, number of grains per panicle, panicles number per plant, and 1000-grain weight) QTLs was conducted, and 462 QTLs were retrieved from 47 independent QTL research published between 2002 and 2022. QTL projection was performed using a reference map with a cumulative length of 2,945.67 cM, and MQTL analysis was conducted on 313 QTLs. Consequently, a total of 62 MQTLs were identified with reduced mean CI (up to 3.40 fold) compared to the mean CI of original QTLs. However, 10 of these MQTLs harbored at least six of the initial QTLs from diverse genetic backgrounds and environments and were considered the most stable and robust MQTLs. Also, MQTLs were compared with GWAS studies and resulted in the identification of 16 common significant loci modulating the evaluated traits. Gene annotation, gene ontology (GO) enrichment, and RNA-seq analyses of chromosome regions of the stable MQTLs detected 52 potential CGs including those that have been cloned in previous studies. These genes encode proteins known to be involved in regulating grain yield including cytochrome P450, zinc fingers, MADs-box, AP2/ERF domain, F-box, ubiquitin ligase domain protein, homeobox domain, DEAD-box ATP domain, and U-box domain. This study provides the framework for molecular dissection of grain yield in rice. Moreover, the MQTLs and CGs identified could be useful for fine mapping, gene cloning, and marker-assisted selection to improve rice productivity.

Candidate loci for leaf angle in maize revealed by a combination of genome-wide association study and meta-analysis

Leaf angle (LA) is a key component of maize plant architecture that can simultaneously govern planting density and improve final yield. However, the genetic mechanisms underlying LA have not been fully addressed. To broaden our understanding of its genetic basis, we scored three LA-related traits on upper, middle, and low leaves of 492 maize inbred lines in five environments. Phenotypic data revealed that the three LA-related traits were normally distributed, and significant variation was observed among environments and genotypes. A genome-wide association study (GWAS) was then performed to dissect the genetic factors that control natural variation in maize LA. In total, 85 significant SNPs (involving 32 non-redundant QTLs) were detected (p ≤ 2.04 × 10-6), and individual QTL explained 4.80%-24.09% of the phenotypic variation. Five co-located QTL were detected in at least two environments, and two QTLs were co-located with multiple LA-related traits. Forty-seven meta-QTLs were identified based on meta-analysis combing 294 LA-related QTLs extracted from 18 previously published studies, 816 genes were identified within these meta-QTLs, and seven co-located QTLs were jointly identified by both GWAS and meta-analysis. ZmULA1 was located in one of the co-located QTLs, qLA7, and its haplotypes, hap1 and hap2, differed significantly in LA-related traits. Interestingly, the temperate materials with hap2 had smallest LA. Finally, we also performed haplotype analysis using the reported genes that regulate LA, and identified a lot of maize germplasms that aggregated favorable haplotypes. These results will be helpful for elucidating the genetic basis of LA and breeding new maize varieties with ideal plant architecture.

Mega Meta-QTLs: A Strategy for the Production of Golden Barley (Hordeum vulgare L.) Tolerant to Abiotic Stresses

Abiotic stresses cause a significant decrease in productivity and growth in agricultural products, especially barley. Breeding has been considered to create resistance against abiotic stresses. Pyramiding genes for tolerance to abiotic stresses through selection based on molecular markers connected to Mega MQTLs of abiotic tolerance can be one of the ways to reach Golden Barley. In this study, 1162 original QTLs controlling 116 traits tolerant to abiotic stresses were gathered from previous research and mapped from various populations. A consensus genetic map was made, including AFLP, SSR, RFLP, RAPD, SAP, DArT, EST, CAPS, STS, RGA, IFLP, and SNP markers based on two genetic linkage maps and 26 individual linkage maps. Individual genetic maps were created by integrating individual QTL studies into the pre-consensus map. The consensus map covered a total length of 2124.43 cM with an average distance of 0.25 cM between markers. In this study, 585 QTLs and 191 effective genes related to tolerance to abiotic stresses were identified in MQTLs. The most overlapping QTLs related to tolerance to abiotic stresses were observed in MQTL6.3. Furthermore, three MegaMQTL were identified, which explained more than 30% of the phenotypic variation. MQTLs, candidate genes, and linked molecular markers identified are essential in barley breeding and breeding programs to develop produce cultivars resistant to abiotic stresses.

Meta-analysis and co-expression analysis revealed stable QTL and candidate genes conferring resistances to Fusarium and Gibberella ear rots while reducing mycotoxin contamination in maize

Fusarium (FER) and Gibberella ear rots (GER) are the two most devastating diseases of maize (Zea mays L.) which reduce yield and affect grain quality worldwide, especially by contamination with mycotoxins. Genetic improvement of host resistance to effectively tackle FER and GER diseases requires the identification of stable quantitative trait loci (QTL) to facilitate the application of genomics-assisted breeding for improving selection efficiency in breeding programs. We applied improved meta-analysis algorithms to re-analyze 224 QTL identified in 15 studies based on dense genome-wide single nucleotide polymorphisms (SNP) in order to identify meta-QTL (MQTL) and colocalized genomic loci for fumonisin (FUM) and deoxynivalenol (DON) accumulation, silk (SR) and kernel (KR) resistances of both FER and GER, kernel dry-down rate (KDD) and husk coverage (HC). A high-resolution genetic consensus map with 36,243 loci was constructed and enabled the projection of 164 of the 224 collected QTL. Candidate genes (CG) mining was performed within the most refined MQTL, and identified CG were cross-validated using publicly available transcriptomic data of maize under Fusarium graminearum infection. The meta-analysis revealed 40 MQTL, of which 29 were associated each with 2-5 FER- and/or GER-related traits. Twenty-eight of the 40 MQTL were common to both FER and GER resistances and 19 MQTL were common to silk and kernel resistances. Fourteen most refined MQTL on chromosomes 1, 2, 3, 4, 7 and 9 harbored a total of 2,272 CG. Cross-validation identified 59 of these CG as responsive to FER and/or GER diseases. MQTL ZmMQTL2.2, ZmMQTL9.2 and ZmMQTL9.4 harbored promising resistance genes, of which GRMZM2G011151 and GRMZM2G093092 were specific to the resistant line for both diseases and encoded "terpene synthase21 (tps21)" and "flavonoid O-methyltransferase2 (fomt2)", respectively. Our findings revealed stable refined MQTL harboring promising candidate genes for use in breeding programs for improving FER and GER resistances with reduced mycotoxin accumulation. These candidate genes can be transferred into elite cultivars by integrating refined MQTL into genomics-assisted backcross breeding strategies.

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

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

Genome-wide meta-QTL analyses provide novel insight into disease resistance repertoires in common bean

BACKGROUND

Common bean (Phaseolus vulgaris) is considered a staple food in a number of developing countries. Several diseases attack the crop leading to substantial economic losses around the globe. However, the crop has rarely been investigated for multiple disease resistance traits using Meta-analysis approach.

RESULTS AND CONCLUSIONS

In this study, in order to identify the most reliable and stable quantitative trait loci (QTL) conveying disease resistance in common bean, we carried out a meta-QTL (MQTL) analysis using 152 QTLs belonging to 44 populations reported in 33 publications within the past 20 years. These QTLs were decreased into nine MQTLs and the average of confidence interval (CI) was reduced by 2.64 folds with an average of 5.12 cM in MQTLs. Uneven distribution of MQTLs across common bean genome was noted where sub-telomeric regions carry most of the corresponding genes and MQTLs. One MQTL was identified to be specifically associated with resistance to halo blight disease caused by the bacterial pathogen Pseudomonas savastanoi pv. phaseolicola, while three and one MQTLs were specifically associated with resistance to white mold and anthracnose caused by the fungal pathogens Sclerotinia sclerotiorum and Colletotrichum lindemuthianum, respectively. Furthermore, two MQTLs were detected governing resistance to halo blight and anthracnose, while two MQTLs were detected for resistance against anthracnose and white mold, suggesting putative genes governing resistance against these diseases at a shared locus. Comparative genomics and synteny analyses provide a valuable strategy to identify a number of well‑known functionally described genes as well as numerous putative novels candidate genes in common bean, Arabidopsis and soybean genomes.

A comprehensive meta-QTL analysis for yield-related traits of durum wheat (Triticum turgidum L. var. durum) grown under different water regimes

Abiotic stress strongly affects yield-related traits in durum wheat, in particular drought is one of the main environmental factors that have effect on grain yield and plant architecture. In order to obtain new genotypes well adapted to stress conditions, the highest number of desirable traits needs to be combined in the same genotype. In this context, hundreds of quantitative trait loci (QTL) have been identified for yield-related traits in different genetic backgrounds and environments. Meta-QTL (MQTL) analysis is a useful approach to combine data sets and for creating consensus positions for the QTL detected in independent studies for the reliability of their location and effects. MQTL analysis is a useful method to dissect the genetic architecture of complex traits, which provide an extensive allelic coverage, a higher mapping resolution and allow the identification of putative molecular markers useful for marker-assisted selection (MAS). In the present study, a complete and comprehensive MQTL analysis was carried out to identify genomic regions associated with grain-yield related traits in durum wheat under different water regimes. A total of 724 QTL on all 14 chromosomes (genomes A and B) were collected for the 19 yield-related traits selected, of which 468 were reported under rainfed conditions, and 256 under irrigated conditions. Out of the 590 QTL projected on the consensus map, 421 were grouped into 76 MQTL associated with yield components under both irrigated and rainfed conditions, 12 genomic regions containing stable MQTL on all chromosomes except 1A, 4A, 5A, and 6B. Candidate genes associated to MQTL were identified and an in-silico expression analysis was carried out for 15 genes selected among those that were differentially expressed under drought. These results can be used to increase durum wheat grain yields under different water regimes and to obtain new genotypes adapted to climate change.

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.

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-analysis of QTLs and candidate genes associated with seed germination in rice (Oryza sativa L.)

Seed germination is one of the critical stages of plant life, and many quantitative trait loci (QTLs) control this complex trait. Meta-analysis of QTLs is a powerful computational technique for estimating the most stable QTLs regardless of the population's genetic background. Besides, this analysis effectively narrows down the confidence interval (CI) to identify candidate genes (CGs) and marker development. In the current study, a comprehensive genome-wide meta-analysis was performed on QTLs associated with germination in rice. This analysis was conducted based on the data reported over the last two decades. In this case, various analyses were performed, including seed germination rate, plumule length, radicle length, germination percentage, coleoptile length, coleorhiza length, radicle fresh weight, germination potential, and germination index. A total of 67 QTLs were projected onto a reference map for these traits and then integrated into 32 meta-QTLs (MQTLs) to provide a genetic framework for seed germination. The average CI of MQTLs was considerably reduced from 15.125 to 8.73 cM compared to the initial QTLs. This situation identified 728 well-known functionally characterized genes and novel putative CGs for investigated traits. The fold change calculation demonstrated that 155 CGs had significant changes in expression analysis. In this case, 112 and 43 CGs were up-regulated and down-regulated during germination, respectively. This study provides an overview and compares genetic loci controlling traits related to seed germination in rice. The findings can bridge the gap between QTLs and CGs for seed germination.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12298-022-01232-1.

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-QTL analysis and candidate genes identification for various abiotic stresses in maize (Zea mays L.) and their implications in breeding programs

Global climate change leads to the concurrence of a number of abiotic stresses including moisture stress (drought, waterlogging), temperature stress (heat, cold), and salinity stress, which are the major factors affecting maize production. To develop abiotic stress tolerance in maize, many quantitative trait loci (QTL) have been identified, but very few of them have been utilized successfully in breeding programs. In this context, the meta-QTL analysis of the reported QTL will enable the identification of stable/real QTL which will pave a reliable way to introgress these QTL into elite cultivars through marker-assisted selection. In this study, a total of 542 QTL were summarized from 33 published studies for tolerance to different abiotic stresses in maize to conduct meta-QTL analysis using BiomercatorV4.2.3. Among those, only 244 major QTL with more than 10% phenotypic variance were preferably utilised to carry out meta-QTL analysis. In total, 32 meta-QTL possessing 1907 candidate genes were detected for different abiotic stresses over diverse genetic and environmental backgrounds. The MQTL2.1, 5.1, 5.2, 5.6, 7.1, 9.1, and 9.2 control different stress-related traits for combined abiotic stress tolerance. The candidate genes for important transcription factor families such as ERF, MYB, bZIP, bHLH, NAC, LRR, ZF, MAPK, HSP, peroxidase, and WRKY have been detected for different stress tolerances. The identified meta-QTL are valuable for future climate-resilient maize breeding programs and functional validation of candidate genes studies, which will help to deepen our understanding of the complexity of these abiotic stresses.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-022-01294-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.

Unravelling consensus genomic regions conferring leaf rust resistance in wheat via meta-QTL analysis

Leaf rust, caused by the fungus Puccinia triticina Erikss (Pt), is a destructive disease affecting wheat (Triticum aestivum L.) and a threat to food security. Developing resistant cultivars represents a useful method of disease control, and thus, understanding the genetic basis for leaf rust resistance is required. To this end, a comprehensive bibliographic search for leaf rust resistance quantitative trait loci (QTL) was performed, and 393 QTL were collected from 50 QTL mapping studies. Afterward, a consensus map with a total length of 4,567 cM consisting of different types of markers (simple sequence repeat [SSR], diversity arrays technology [DArT], chip-based single-nucleotide polymorphism [SNP] markers, and SNP markers from genotyping-by-sequencing) was used for QTL projection, and meta-QTL (MQTL) analysis was performed on 320 QTL. A total of 75 MQTL were discovered and refined to 15 high-confidence MQTL (hcmQTL). The candidate genes discovered within the hcmQTL interval were then checked for differential expression using data from three transcriptome studies, resulting in 92 differentially expressed genes (DEGs). The expression of these genes in various leaf tissues during wheat development was explored. This study provides insight into leaf rust resistance in wheat and thereby provides an avenue for developing resistant cultivars by incorporating the most important hcmQTL.

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.

Genome-Wide Expression Analysis of Root Tips in Contrasting Rice Genotypes Revealed Novel Candidate Genes for Water Stress Adaptation

Root system architecture (RSA) is an important agronomic trait with vital roles in plant productivity under water stress conditions. A deep and branched root system may help plants to avoid water stress by enabling them to acquire more water and nutrient resources. Nevertheless, our knowledge of the genetics and molecular control mechanisms of RSA is still relatively limited. In this study, we analyzed the transcriptome response of root tips to water stress in two well-known genotypes of rice: IR64, a high-yielding lowland genotype, which represents a drought-susceptible and shallow-rooting genotype; and Azucena, a traditional, upland, drought-tolerant and deep-rooting genotype. We collected samples from three zones (Z) of root tip: two consecutive 5 mm sections (Z1 and Z2) and the following next 10 mm section (Z3), which mainly includes meristematic and maturation regions. Our results showed that Z1 of Azucena was enriched for genes involved in cell cycle and division and root growth and development whereas in IR64 root, responses to oxidative stress were strongly enriched. While the expansion of the lateral root system was used as a strategy by both genotypes when facing water shortage, it was more pronounced in Azucena. Our results also suggested that by enhancing meristematic cell wall thickening for insulation purposes as a means of confronting stress, the sensitive IR64 genotype may have reduced its capacity for root elongation to extract water from deeper layers of the soil. Furthermore, several members of gene families such as NAC, AP2/ERF, AUX/IAA, EXPANSIN, WRKY, and MYB emerged as main players in RSA and drought adaptation. We also found that HSP and HSF gene families participated in oxidative stress inhibition in IR64 root tip. Meta-quantitative trait loci (QTL) analysis revealed that 288 differentially expressed genes were colocalized with RSA QTLs previously reported under drought and normal conditions. This finding warrants further research into their possible roles in drought adaptation. Overall, our analyses presented several major molecular differences between Azucena and IR64, which may partly explain their differential root growth responses to water stress. It appears that Azucena avoided water stress through enhancing growth and root exploration to access water, whereas IR64 might mainly rely on cell insulation to maintain water and antioxidant system to withstand stress. We identified a large number of novel RSA and drought associated candidate genes, which should encourage further exploration of their potential to enhance drought adaptation in rice.