Ilaria Marcotuli
Jose Miguel Soriano
Agata Gadaleta- 1Department of Agricultural and Environmental Science, University of Bari Aldo Moro, Bari, Italy
- 2Sustainable Field Crops Programme, IRTA (Institute for Food and Agricultural Research and Technology), Lleida, Spain
The present work focused on the identification of durum wheat QTL hotspots from a collection of genome-wide association studies, for quality traits, such as grain protein content and composition, yellow color, fiber, grain microelement content (iron, magnesium, potassium, selenium, sulfur, calcium, cadmium), kernel vitreousness, semolina, and dough quality test. For the first time a total of 10 GWAS studies, comprising 395 marker-trait associations (MTA) on 57 quality traits, with more than 1,500 genotypes from 9 association panels, were used to investigate consensus QTL hotspots representative of a wide durum wheat genetic variation. MTA were found distributed on all the A and B genomes chromosomes with minimum number of MTA observed on chromosome 5B (15) and a maximum of 45 on chromosome 7A, with an average of 28 MTA per chromosome. The MTA were equally distributed on A (48%) and B (52%) genomes and allowed the identification of 94 QTL hotspots. Synteny maps for QTL were also performed in Zea mays, Brachypodium, and Oryza sativa, and candidate gene identification allowed the association of genes involved in biological processes playing a major role in the control of quality traits.
Introduction
Cereals are the main species in Mediterranean cropping systems, well adapted to semi-arid climate conditions and able to give a stable economic sustain to the farmers. Wheat represents a target crop for Mediterranean agriculture and provides 20% of calories to the world population, highlighting the relevance of this crop for current and future strategic cultivation (Royo et al., 2017). In this context, we must intensify efforts toward crop improvement and yield stability under conditions of sustainable agricultural production.
Durum wheat (Triticum turgidum ssp. durum) is largely produced in the Mediterranean basin being used for human nutrition, prevalently transformed into semolina for pasta and couscous, but it can also be used to obtain flour for bread (Kadkol and Sissons, 2016). Due to the ongoing climate changes and steady increase in average temperatures, the flour quality of commercial genotypes may no longer be able to perform well in the coming years. The importance of grain quality parameters for durum wheat end-products in the food chain makes it a crucial tool in maintaining or increasing durum wheat production under disease pressure and adverse climatic conditions to preserve the grain quality (Beres et al., 2020).
Since crop quality is a complex trait, it can be either related to end-use properties or to nutritional content and these parameters are regulated by different compounds, the genetic determination of quality traits is rather complicated. Due to the high environmental influence on crop quality traits, the search for genes related to them is more complicated such as the transfer of these characters, despite the presence of novel genomic tools (Yang et al., 2022).
The research activities for wheat quality improvement have been focused in the last decade on the valorization of wheat germplasm collections including old varieties and wild relatives besides the obtainment and evaluation of new breeding lines. Several studies have been conducted on the environmental effects on quality, development of evaluation methods, and processing for end-users (Pour-Aboughadareh et al., 2021). Therefore, the nutritional improvement of crops has a positive impact on millions of people around the world without the need to alter their eating habits. Currently, the aspects related to human health have taken great relevance in cereal breeding programs to develop biofortified crops (Swamy et al., 2021).
The increasing awareness of the use of cereal-based products in a healthy diet is currently becoming more and more evident, and so plenty of studies focus on the identification and exploitation of natural variations of bioactive components in the grain. Wheat has a protein content of about 13% and is the leading source of vegetable protein in human food, and it is also an important source of carbohydrates. When wheat is eaten as a whole grain, it is an excellent source of dietary fiber and nutrients (Marcotuli et al., 2020). Many of the abovementioned traits are polygenic traits associated with quantitative trait loci (QTL) located on all the tetraploid wheat chromosomes. To identify main effective genes for wheat quality is a main target of wheat breeders worldwide, however, the efficiency of selection is constrained by the following: 1) the decline of genetic diversity in elite germplasm (by the pursuit of elite high-performing cultivars) that leads to scarcity or even absence of suitable loci in modern breeding lines; 2) limited knowledge on quality traits that are often complex QTL influenced by numerous genes and environmental conditions; 3) lack of adequate molecular knowledge to lay the foundation for molecular breeding. These limitations are more severe in durum wheat with respect to bread wheat. The underlined problem can be solved with the use of high-density genetic maps, new molecular markers, or a wide collection summarized in GWAS analyses.
GWAS detects the association between genotype and trait of interest using conserved linkage disequilibrium (LD) present in a selected panel of accessions (Myles et al., 2009). Recently, association or LD mapping, utilizing genome-wide markers, has been adopted in wheat because of two main advantages: 1) association mapping does not require the cost and time associated with the population and genetic map development, and 2) GWAS provides high mapping resolution as it efficiently uses the multiple historical crossover events occurred in the diverse association panel used (Saini et al., 2022).
GWAS give, also, information about MTAs that could be utilized for candidate gene discovery and characterization adoptable in breeding programs. A study from Saini et al. (2022) underlines that among 86,122 wheat lines studied under various GWAS analyses 46,940 loci associated with traits were reported, but further utilization of these markers was largely limited. To solve this situation QTL hotspots (genomic locations enriched in QTL) are a common and notable feature when collecting many QTL for various traits in many areas of biological studies (Wu et al., 2021). This approach is a good instrument to study at the same time many traits, trying to find the consensus and most robust QTL using the information reported in multiple studies for the reliability of their location and effect across different genetic backgrounds and environments, as well as to refine QTL positions on a consensus map (Goffinet and Gerber, 2000). The approach developed by these authors is called QTL meta-analysis and it allows the identification of the genome regions most involved in trait variation, being a suitable approach to narrow down QTL regions, to tackle map-based cloning strategies, and to identify candidate genes (Soriano and Alvaro, 2019; Wu et al., 2021). Regarding GWAS, a statistical framework for QTL hotspot detection was reported by Wu et al. (2021), allowing the integration of independent GWAS studies in a reference map for durum wheat. The use of MTA in meta-analysis represents an important and complementary tool for the identification of new stable Genome-wide QTL hotspots and the discovery of new alleles for durum wheat quality. The most powerful hotspots have also been used for the detection of the so-called ortho-Meta-QTL, which are conserved among species and hence will be more reliable and can be used in different cereals.
In the present manuscript, MTA detected by GWAS were projected for the first time on the durum wheat consensus map from Maccaferri et al. (2019) allowing the identification of QTL hotspots related to quality traits in durum wheat, genes directly associated and controlling different quality traits and regulatory genes, and finally conserved orthologous regions among bread wheat, durum wheat, rice, maize, and Brachypodium.
Materials and methods
Quality traits and marker-trait association database
Genome-wide association studies for quality traits were retrieved using the keywords “durum wheat quality GWAS” from the Web of Science server (http://apps.webofknowledge.com). The marker-trait association (MTA) for the GWAS studies database was created from 10 studies published from 2015 to 2020 (Table 1). The database reported information on the diversity panel used, number of genotypes, number of MTA, and traits analyzed in each study. A total of 57 different traits were collected from the GWAS studies (Table 2). The MTA database reported information on the name of the chromosome and position of the MTA, the confidence interval (CI), and the phenotypic variance explained (PVE) by each MTA.
TABLE 1. Summary of QTL studies included in the meta-analysis.
TABLE 2. Traits analyzed for QTL meta-analysis.
Projection on the consensus map
To represent all the MTA in the same linkage map, the durum wheat consensus map developed by Maccaferri et al. (2015) was used for projection following the homothetic approach described by Chardon et al. (2004) as previously described in Colasuonno et al. (2021). When the CIs were not reported in the original studies, they were estimated according to the linkage disequilibrium (LD) decay for each chromosome.
Identification of quantitative trait loci hotspots
To simplify the MTA information, the associations were grouped into QTL hotspots. First, the CIs were standardized using the formula described by Chardon et al. (2004):
where CI corresponded with the original CI or the LD decay for each chromosome. To define a hotspot, the density of MTAs along the chromosome was calculated as the QTL overview index (Chardon et al., 2004) for each cM of the genetic map reported by Maccaferri et al. (2015):
where nbQTL is the number of MTA and nbE is the total number of studies.
Breeding QTL hotspots were selected based on three criteria according to Löffler et al. (2009): small supporting intervals of the mQTL, high number of initial QTL, and high phenotypic variance explained (PVE) by the initial QTL. In our case, we considered for QTL hotspots a CI lower than 20Mb, and subsequently, a PVE > 0.04 when 5 or more MTAs per hotspot were reported or PVE > 0.1 when MTAs per hotspots was 3. QTL hotspot distribution was compared among the genetic and physical maps of durum wheat (Maccaferri et al., 2015; Maccaferri et al., 2019). Chromosomes were equally distributed in five bins based on their genetic and physical length. Subsequently, the total number of QTL hotspots per bin was counted.
Synteny analysis and identification of candidate genes underlying the MQTL
A total of 12,606 markers associated with QTL for quality from Triticum ssp. were palmed on the genomes of O. sativa, Z. mays, and B. distachyon to extract the syntenic positions of the markers in the three genomes using a mapping approach and draw a circus diagram for data view for each analysis. The following reference genomes were used: B. distachyon version 1 (http://www.plantgdb.org/BdGDB/); O. sativa Japonica Group version IRGSP-1.0 (http://rice.uga.edu/); and Z. mays version AGPv3 (https://www.maizegdb.org/).
An ad-hoc mapping pipeline was developed to map the markers and highlight all the syntenic positions between Triticum spp. and three different monocot genomes. The core algorithm of this pipeline was the bwa-mem aligner, and the workflow was divided into three parts: 1) perform three different mappings for each genome with different quality parameters, from higher to lower stringency; 2) merge the results of the five mapping steps considering several quality parameters including the redundancy of the alignments, the quality of the markers, and the quality of the alignments; 3) filter by mapping length >60%, which means to keep all the hits with a mapping length >60% with respect to the marker sizes, and a spread-scaled quality >10, which measures the posterior probability that the mapping position is wrong.
Gene models within syntenic QTL hotspot were identified using the high-confidence genes reported for the Svevo durum wheat reference sequence, available at https://wheat.pw.usda.gov/GG3/jbrowse_Durum_Svevo, and the bread wheat Chinese spring reference sequence, available at https://wheat-urgi.versailles.inra.fr/Seq-Repository. Homologous genes from “Chinese Spring” and “Svevo” were subsequently identified in the syntheny region of Maize (https://www.maizegdb.org/), Rice (http://rice.uga.edu/), and Brachypodium (http://www.plantgdb.org/BdGDB/) databases.
Circular figures of the QTL hotspots and chromosome synteny were created using the online software Clico FS (Cheong et al., 2015) available at http://clicofs.codoncloud.com.
Results
Marker-trait association distribution and quantitative trait loci hotspots identification
Ten studies published from 2015 to 2020 reporting 395 MTA for quality traits were collected (Table 1; Supplementary Table S1). The studies covered a total of 1,598 genotypes (including elite cultivars, breeding lines, and landraces). MTA was distributed throughout the 14 chromosomes (A and B genomes) of durum wheat. The number of MTA per chromosome ranged from 15 on chromosome 5B to 45 on chromosome 7A, with an average of 28 MTA per chromosome. Forty-eight percent of the MTA were identified in genome A, and 52% in genome B (Figure 1A). When the trait was considered, out of 39 traits, yellow pigment content (YPC) presented the highest number of MTA (24%). Confidence intervals (CI) ranged from 0.1 to 43 cm with an average of 3.3 cm. For 203 MTA not reporting a CI in the original study, the intra-chromosomal LD decay was used to establish their CI. Most of the MTA (97%) showed a CI lower or equal to 10 cm, whereas the CI for 83% of the MTAs was lower or equal to 5 cm (Figure 1B). The proportion of the phenotypic variance explained (PVE) followed a typical quantitative inheritance as previously reported in other QTL meta-analysis studies (Soriano and Alvaro, 2019; Soriano et al., 2021). 79% of the MTA showed a PVE lower than 0.1, whereas it increases to 95% for a PVE < 0.2 (Figure 1C).
FIGURE 1. MTA statistics. (A) Number of MTA per chromosome. (B) Confidence interval from original MTA. (C) PVE from original MTA. (D) Percentage of QTL hotspots per bin in the genetic and physical maps.
In order to simplify the MTA information, QTL hotspots were defined using the QTL overview density index defined by Chardon et al. (2004) for each centimeter of the durum wheat consensus map. A total of 564 peaks were identified using the mean of the overview index across the 14 chromosomes (0.15) as the threshold, whereas using a high threshold (0.75), a total of 158 peaks were detected (Figure 2). These 158 peaks were reduced to 92 QTL hotspots (Supplementary Table S2), 47 in genome A and 45 in genome B. QTL hotspots were selected based on three criteria for candidate gene identification. First, hotspots with a physical distance lower than 20 Mb were identified, and later only those containing three or more MTA with a PVE higher or equal to 0.04 were chosen. As a result, 20 QTL hotspots including 121 MTA with a PVE mean of 0.11 were selected (Table 3). These hotspots were considered as breeding QTL.
FIGURE 2. QTL hotspots defined by the overview density index (red line). Blue boxes represent the projected MTAs on the durum wheat consensus map.
TABLE 3. Selected QTL hotspots. PVE refers to the mean of all the MTA in the hotspot.
The comparison of the breeding QTL coverage in the genetic and physical maps resulted in different distribution along chromosomes. Consensus genetic and physical chromosomes were defined based on the differentiation of five bins per chromosome according to its genetic and physical distance respectively, and each bin corresponding to 1/5 of the total length. Whereas for the genetic chromosome the QTL hotspots were equally distributed in each one of the bins (with 16% on bin 3, 21% on bins 1, 2, and 5, and 22% on bin 4), for the physical chromosome and unequal distribution was observed, with the higher number of QTL hotspots on telomeric bins (83% on bins 1 and 5, 12% on bin 4, and 2 and 3% on bins 2 and 3, respectively) (Figure 1D).
Synteny among species and MQTL candidate genes identification
Thanks to the pipeline described in Materials and Method, it was possible to identify syntenic genome regions between the Triticum spp., markers associated with the breeding QTL on three genomes (Brachypodium, Z. mays, and O. sativa) (Figure 3; Table 4; Supplementary Table S3).
FIGURE 3. Synteny among genomes for the selected quality QTL hotspots. For better visualization of the chromosome links the size of the Brachypodium and rice chromosomes is multiplied by 10.
TABLE 4. Synteny among species in the QTL hotspot regions. bd, Brachypodium distachyum; os, Oryza sativa; zm, Zea mays.
A total of 134 markers from the breeding QTL were identified in syntenic regions (ortho-MQTL) of the Brachypodium (78), rice (49), and maize (7) genomes. Chromosome 6A showed the highest number of markers with syntenic regions with other genomes, whereas chromosomes 3A, 4A, 5B, 6B, and 7B did not show any homologous marker. Genome A represented 70% of the syntenic markers.
Using the genomic sequences located in the breeding QTL from Svevo, Z. mays, O. sativa, and B. distachyum, candidate genes correlated to ortho-MQTL were identified (Table 5).
TABLE 5. Genes identified through meta-QTL analysis in each qhotspot and detected in the following databases: bread wheat cultivar Chinese Spring (CS) genome sequence, durum wheat cultivar Svevo genome sequence, Brachypodium db, Rice db, and MaizeGDB.
The investigation of collinear regions within the five genomes resulted in the identification of 17 synthenic regions containing paralogous genes with similar functions that can be considered as promising candidate genes controlling the quality trait considered (Table 5). On qhotspot1, which was correlated to the dough extensibility, firmness, grain protein content, and gluten index and strength, the genes gamma gliadin-A1, A3, A4, and LMW-A2, HMW glutenin A gene, and LMW glutenin subunit (GLU-3) were identified on Chinese Spring, Svevo and Brachypodium, respectively confirming the correlation between the genomic region considered and the traits control. Instead, two different genes were identified in the same qhotspot1 for rice and maize (ATPase, AAA family domain, and paired amphipathic helix protein Sin3-like 3, respectively).
On qhotspot4.qhotspot10, qhotspot15, qhotspot35, qhotspot47, qhotspot52, qhotspot64, and qhotspot81 were identified in Svevo, maize, rice, and Brachypodium genes for biological processes, such as DNA replication, cell wall development, and secondary metabolite production, with no direct correlation with the traits corresponding to the QTL regions (Table 5). No colinear genes were detected on qhotspot5 between the species analyzed, but candidate genes correlated to the dough extensibility and tenacity, grain protein content, gluten index, and strength (gamma gliadin-B1, B2, B4, B6, delta gliadin-B1, omega gliadin-B3, B6, LMW-B2, and LMW-B3 genes) were detected in Svevo.
Qhotspot26, qhotspot73, qhotspot79, qhotspot85, and qhotspot92 showed the same genes for all the genomes considered, involved in biological process, but not directly correlated with the traits considered, except for qhotspot92, which showed a correlation between trait IP and sucrose synthase genes identified in Chinese Spring, Svevo and Brachypodium. The other two correlations were detected between traits and genes on qhotspot33 and qhotspot57 in Brachypudium and rice, respectively, but no collinearity was found between species.
Discussion
The detection of major loci for quantitative traits (QTL) is a key tool of modern genetics to identify candidate genes controlling polygenic traits. Kernel quality traits are often quantitative, so controlled by multiple genes/QTL, strongly influencing the end-product quality of commercial wheat varieties, and determining the type of products that can be produced. Thus, it is possible to distinguish a commercial value at the beginning of wheat production and a technological value linked to the worldwide market requirements for end-product uses (Colasuonno et al., 2021). The development of new genotypes with high-quality value is an important tool for breeders besides industrial and consumer requirements. (Colasuonno et al., 2021).
The study of quantitative traits is complex and required the identification of molecular markers and/or genes tightly linked to quality traits to be used in Marker Assisted Selection (MAS) programs. A modern approach based on GWAS, and meta-QTL analysis allowed the evaluation of wide collections identifying numerous new genes/QTL and the comparison of results through meta-analysis. In the present study, we provide a wide study for the evaluation of genetic loci controlling kernel quality traits in durum wheat through the identification of new stable genome-wide QTL hotspots, and correlation through a synteny analysis between different species including Z. mays, B. distachyon, and O. sativa.
The study has been focused on 10 GWAS studies reporting 57 quality traits recorded in wide collections and representative of durum wheat genetic variability for protein content and composition (HMW, LMW), yellow color, fiber, grain microelement content such as iron, magnesium, potassium, selenium, sulfur, calcium, cadmium, kernel vitreousness, semolina, and dough quality test (all reported in Table 1). A total of 395 MTA, equally distributed on the A and B genomes, were used in genome-wide QTL hotspot detection. In order to include all the MTA in a single map they were projected on the durum wheat consensus map developed by Maccaferri et al. (2015) as described by Chardon et al. (2004) and Colasuonno et al. (2021), estimating the CIs, not reported in the original studies, according to the linkage disequilibrium (LD) decay for each chromosome (as reported in materials and method). 97% of the MTAs reported a CI lower than 10 cm, with a phenotypic variance explained (PVE) lower than 0.1, these values are in accordance with those reported by Soriano and Alvaro (2019), Soriano et al. (2021) for quantitative inheritance.
In order to reduce the number of MTA and simplify the information, 92 QTL hotspots (47 in genome A and 45 in genome B) were identified using the QTL overview density index defined by Chardon et al. (2004). Some features were observed when comparing the distribution pattern of QTL hotspots in the genetic and physical maps. Chromosomes were divided into five bins (genetic and physical) and percentage of QTL hotspots per bin was calculated (as it is reported in the Material and Methods section). While QTL hotspots were equally distributed along the genetic chromosomes, a different distribution was observed for the physical chromosomes with higher number of QTL hotspots on telomeric bins (83% on bins 1 and 5). This is in accordance with what was reported by several authors in QTL mapping, which observed that QTL are highly abundant in some genomic regions, and, frequently, correlated traits are clustered closely together in some specific loci as compared to others (Goffinet and Gerber, 2000; Schadt et al., 2003; Chardon et al., 2004; West et al., 2007; Breitling et al., 2008; Wu et al., 2008; Wang et al., 2014; Basnet et al., 2015; Martinez et al., 2016; Yang et al., 2019).
A clustering between QTL and traits has also been reported in classical MQTL analysis by Soriano et al. (2021) that grouped 318 QTL including quality, biotic and abiotic traits in 85 MQTL with number of traits involved in each MQTL ranging from 1 to 12, and the number of QTL per MQTL ranged from 2 to 11. This phenomenon may have several causes, including QTL with high allelic polymorphisms and interesting pleiotropic effect or closely linked QTL controlling correlated traits and frequently co-localized in the same regions (Zhao et al., 2011; Vuong et al., 2015; Mengistu et al., 2016; Zhang et al., 2020). As the QTL hotspots can lead to the identification of genes that affect the traits of interest, and further help to build networks among QTL hotspots, genes, and traits, the qhotspot detection analysis at genome-wide level has been a key step toward deciphering the genetic architectures of quantitative traits in genes, genomes, and genetics studies (Breitling et al., 2008; Fu et al., 2009; Neto et al., 2012; Wang et al., 2014; Yang et al., 2019; Wu et al., 2021).
Previous studies identified ortho-MQTL for yield-related traits and nitrogen use efficiency (Saini et al., 2021; Singh et al., 2022). According to these authors, the conserved nature of the ortho-MQTL suggests that may be associated with regulatory elements affecting many genes.
In this study, we used the breeding QTL for the investigation of collinear regions between durum and bread wheat, with Brachypodium, maize, and rice.
The candidate gene investigation allowed the correlation between some ortho-MQTL with genes, such as for qhotspot1 which showed collinearity between Chinese Spring, Svevo, and Brachypodium, a clear association with DOE, FIRM, GPC, GS, and genes for gliadin and glutenin. A similar association was investigated on qhotspot5 where it was detected on the Chinese Spring genome again genes for gliadin and glutenin associated with DOE, DOTE, GPS, and GS. Both these associations are in accordance with data reported in the literature which located genes encoding glutenin and gliadin on chromosome groups 1 (Gli-A1, Gli-B1, Gli-D1, Glu-A1, Glu-B1, Glu-D1, Glu-A3, Glu-B3, and Glu-D3 loci) and 6 (Gli-A2, Gli-B2, and Gli-D2 loci) (Zaefizadeh et al., 2010; Dubois et al., 2016; Utebayev et al., 2019; Asri et al., 2021).
Another association was detected on qhotspot26, which was mainly related to protein content traits and showed co-location with potassium ion transport in wheat. In fact, the influence of potassium transport efficiency in kernel protein content has been largely studied and reported in the literature (Beaton and Sekhon, 1985). Analyzing the collinearity with the other genomes used in this study, the syntheny with rice, maize, and Brachypodium was also detected.
An important correlation was determined on qhotspot92 between the gliadins (IP traits) and the gene S. synthase, which was found in Chinese Spring, Svevo, and Bachypodium genome sequences. The association between this trait and the gene was already highlighted by Marín-Sanz et al. (2022) who reported that genes implicated in starch synthesis, including S. synthase, were up-or down-regulated in RNAi lines characterized by a strong reduction in α/β-, ω-, and γ-gliadins. All the other association qhotspot/genes corresponded to genes involved in biological processes, mainly regulatory or transport genes.
Conclusion
This study aimed to a better understanding of the genetic architecture controlling quality traits in durum wheat. GWAS is considered a powerful and complementary tool to QTL analysis for the detection of genes directly correlated to phenotypic traits. However, a meta-analysis of GWAS permits to identify consensus regions among different germplasm controlling different traits and regulatory genes which for some of them play a major role in the control of the trait per se. Comparative genomics analysis among bread and durum wheat, rice, maize, and Brachypodium allowed the identification of the conserved orthologous set of DNA sequences for candidate genes underpinning quantitative traits. The strict relation between comparative genomics and gene identification emphasizes the importance of the orthologous genes identified, opening the possibility to transfer information across species.
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material. Further inquiries can be directed to the corresponding authors.
Author contributions
All authors listed have made a substantial, direct, and intellectual contribution to the work and approved it for publication.
Funding
This research was funded by the projects PRIMA 2019 “CEREALMED” (Italy) and PID 2019-109089RB-C31 (Ministerio de Ciencia e Innovación, Spain).
Acknowledgments
JS acknowledges the contribution of the CERCA Program (Generalitat de Catalunya) (http://cerca.cat/).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors, and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2022.982418/full#supplementary-material
References
Alemu, A., Feyissa, T., Tuberosa, R., Maccaferri, M., Sciara, G., Letta, T., et al. (2020). Genome-wide association mapping for grain shape and color traits in Ethiopian durum wheat (Triticum turgidum ssp. durum). Crop J. 8, 757–768. doi:10.1016/j.cj.2020.01.001
Asri, N., Rostami-Nejad, M., Anderson, R. P., and Rostami, K. (2021). The gluten gene: unlocking the understanding of gluten sensitivity and intolerance. Appl. Clin. Genet. 14, 37–50. doi:10.2147/TACG.S276596
Basnet, R. K., Duwal, A., Tiwari, D. N., Xiao, D., Monakhos, S., Bucher, J., et al. (2015). Quantitative trait locus analysis of seed germination and seedling vigor in Brassica rapa reveals OTL hotspots and epistatic interactions. Front. Plant Sci. 6, 1032. doi:10.3389/fpls.2015.01032
Beaton, J., and Sekhon, G. (1985). “Potassium nutrition of wheat and other small grains,” in Potassium in agriculture. Editor R. D. Munson (Hoboken, NJ: Willey), 701–752.
Beres, B. L., Rahmani, E., Clarke, J. M., Grassini, P., Pozniak, C. J., Geddes, C. M., et al. (2020). A systematic review of durum wheat: Enhancing production systems by exploring genotype, environment, and management (G × E × M) synergies. Front. Plant Sci. 11, 568657. doi:10.3389/fpls.2020.568657
Breitling, R., Li, Y., Tesson, B. M., Fu, J., Wu, C., Wiltshire, T., et al. (2008). Genetical genomics: spotlight on QTL hotspots. PLoS Genet. 4, e1000232. doi:10.1371/journal.pgen.1000232
Chardon, F., Virlon, B., Moreau, L., Falque, M., Joets, J., Decousset, L., et al. (2004). Genetic architecture of flowering time in maize as inferred from quantitative trait loci meta-analysis and synteny conservation with the rice genome. Genetics 168, 2169–2185. doi:10.1534/genetics.104.032375
Cheong, W. H., Tan, Y. C., Yap, S. J., and Ng, K. P. (2015). ClicO FS: an interactive web-based service of circos. Bioinformatics 31, 3685–3687. doi:10.1093/bioinformatics/btv433
Colasuonno, P., Lozito, M. L., Marcotuli, I., Nigro, D., Giancaspro, A., Mangini, G., et al. (2017). The carotenoid biosynthetic and catabolic genes in wheat and their association with yellow pigments. BMC Genomics 18, 122. doi:10.1186/s12864-016-3395-6
Colasuonno, P., Marcotuli, I., Gadaleta, A., and Soriano, J. M. (2021). From genetic maps to QTL cloning: an overview for durum wheat. Plants 10, 315. doi:10.3390/plants10020315
Dubois, B., Bertin, P., and Mingeot, D. (2016). Molecular diversity of α-gliadin expressed genes in genetically contrasted spelt (Triticum aestivum ssp. spelta) accessions and comparison with bread wheat (T. aestivum ssp. aestivum) and related diploid Triticum and Aegilops species. Mol. Breed. 36, 152. doi:10.1007/s11032-016-0569-5
Fu, J., Keurentjes, J. J. B., Bouwmeester, H., America, T., Verstappen, F. W. A., Ward, J. L., et al. (2009). System-wide molecular evidence for phenotypic buffering in Arabidopsis. Nat. Genet. 41, 166–167. doi:10.1038/ng.308
Goffinet, B., and Gerber, S. (2000). Quantitative trait loci: a meta-analysis. Genetics 155, 463–473. doi:10.1093/genetics/155.1.463
Johnson, M., Kumar, A., Oladzad-Abbasabadi, A., Salsman, E., Aoun, M., Manthey, F. A., et al. (2019). Association mapping for 24 traits related to protein content, gluten strength, color, cooking, and milling quality using balanced and unbalanced data in durum wheat [Triticum turgidum L. var. durum (desf).]. Front. Genet. 10, 717. doi:10.3389/fgene.2019.00717
Kadkol, G. P., and Sissons, M. (2016). “Durum wheat: overview,” in Encyclopedia of food grains. Editors C. Wrigley, H. Corke, K. K. Seetharaman, and J. Faubion (Amsterdam: Academic Press), 117–124.
Löffler, M., Schoen, C. C., and Miedaner, T. (2009). Revealing the genetic architecture of FHB resistance in hexaploid wheat (Triticum aestivum L.) by QTL meta-analysis. Mol. Breed. 23, 473–488. doi:10.1007/s11032-008-9250-y
Maccaferri, M., Harris, N. S., Twardziok, S. O., Pasam, R. K., Gundlach, H., Spannagl, M., et al. (2019). Durum wheat genome highlights past domestication signatures and future improvement targets. Nat. Genet. 51, 885–895. doi:10.1038/s41588-019-0381-3
Maccaferri, M., Ricci, A., Salvi, S., Milner, S. G., Noli, E., Martelli, P. L., et al. (2015). A high-density, SNP-based consensus map of tetraploid wheat as a bridge to integrate durum and bread wheat genomics and breeding. Plant Biotechnol. J. 13, 648–663. doi:10.1111/pbi.12288
Marcotuli, I., Colasuonno, P., Hsieh, Y. S. Y., Fincher, G. B., and Gadaleta, A. (2020). Non-starch polysaccharides in durum wheat: a review. Int. J. Mol. Sci. 21, 2933. doi:10.3390/ijms21082933
Marcotuli, I., Houston, K., Schwerdt, J. G., Waugh, R., Fincher, G. B., Burton, R. A., et al. (2016). Genetic diversity and genome wide association study of β-glucan content in tetraploid wheat grains. PLoS One 11, e0152590. doi:10.1371/journal.pone.0152590
Marcotuli, I., Houston, K., Waugh, R., Fincher, G. B., Burton, R. A., Blanco, A., et al. (2015). Genome wide association mapping for arabinoxylan content in a collection of tetraploid wheats. PLoS One 10, e0132787. doi:10.1371/journal.pone.0132787
Marín-Sanz, M., Masaru Iehisa, J. C., and Barro, F. (2022). New transcriptomic insights in two RNAi wheat lines with the gliadins strongly down-regulated by two endosperm specific promoters. Crop J. 10, 194–203. doi:10.1016/j.cj.2021.04.009
Martinez, A. K., Soriano, J. M., Tuberosa, R., Koumproglou, R., Jahrmann, T., and Salvi, S. (2016). Yield QTLome distribution correlates with gene density in maize. Plant Sci. 242, 300–309. doi:10.1016/j.plantsci.2015.09.022
Mengistu, D. K., Kidane, Y. G., Catellani, M., Frascaroli, E., Fadda, C., Pe, M. E., et al. (2016). High-density molecular characterization and association mapping in Ethiopian durum wheat landraces reveals high diversity and potential for wheat breeding. Plant Biotechnol. J. 14, 1800–1812. doi:10.1111/pbi.12538
Myles, S., Peiffer, J., Brown, P. J., Ersoz, E. S., Zhang, Z., Costich, D. E., et al. (2009). Association mapping: critical considerations shift from genotyping to experimental design. Plant Cell 21, 2194–2202. doi:10.1105/tpc.109.068437
N’Diaye, A., Haile, J. K., Cory, A. T., Clarke, F. R., Clarke, J. M., Knox, R. E., et al. (2017). Single marker and haplotype-based association analysis of semolina and pasta colour in elite durum wheat breeding lines using a high-density consensus map. PLoS One 12, e0170941. doi:10.1371/journal.pone.0170941
N’Diaye, A., Haile, J. K., Nilsen, K. T., Walkowiak, S., Ruan, Y., Singh, A. K., et al. (2018). Haplotype loci under selection in Canadian durum wheat germplasm over 60 years of breeding: association with grain yield, quality traits, protein loss, and plant height. Front. Plant Sci. 9 (1589), 1589. doi:10.3389/fpls.2018.01589
Neto, E. C., Keller, M. P., Broman, A. F., Attie, A. D., Jansen, R. C., Broman, K. W., et al. (2012). Quantile-based permutation thresholds for quantitative trait loci hotspots. Genetics 191, 1355–1365. doi:10.1534/genetics.112.139451
Pour-Aboughadareh, A., Kianersi, F., Poczai, P., and Moradkhani, H. (2021). Potential of wild relatives of wheat: ideal genetic resources for future breeding programs. Agronomy 11, 1656. doi:10.3390/agronomy11081656
Reimer, S., Pozniak, C. J., Clarke, F. R., Clarke, J. M., Somers, D. J., Knox, R. E., et al. (2008). Association mapping of yellow pigment in an elite collection of durum wheat cultivars and breeding lines. Genome 51, 1016–1025. doi:10.1139/G08-083
Rosello, M., Royo, C., Alvaro, F., Villegas, D., Nazco, R., and Soriano, J. M. (2018). Pasta-making quality QTLome from Mediterranean durum wheat landraces. Front. Plant Sci. 9, 1512. doi:10.3389/fpls.2018.01512
Royo, C., Soriano, J. M., and Alvaro, F. (2017). “Wheat: a crop in the bottom of the mediterranean diet pyramid,” in Mediterranean identities - environment, society, culture. Editor B. Fuerst-Bjelis (London: IntechOpen). doi:10.5772/intechopen.69184
Saini, D. K., Chopra, Y., Pal, N., Chahal, A., Srivastava, P., and Gupta, P. K. (2021). Meta-QTL, ortho-MQTL and candidate genes for nitrogen use efficiency and root system architecture in bread wheat (Triticum aestivum L.) Physiol. Mol. Biol. Plants. 27, 2245–2267. doi:10.1007/s12298-021-01085-0
Saini, D. K., Chopra, Y., Singh, J., Sandhu, K. S., Kumar, A., Bazzer, S., et al. (2022). Comprehensive evaluation of mapping complex traits in wheat using genome-wide association studies. Mol. Breed. 42, 1. doi:10.1007/s11032-021-01272-7
Schadt, E. E., Monks, S. A., Drake, T. A., Lusis, A. J., Che, N., Colinayo, V., et al. (2003). Genetics of gene expression surveyed in maize, mouse and man. Nature 422, 297–302. doi:10.1038/nature01434
Singh, R., Saripalli, G., Gautam, T., Kumar, A., Batra, R., Kumar, J., et al. (2022). Meta-QTL, ortho-MetaQTL and candidate genes for grain Fe and Zn contents in wheat (Triticum aestivum L.) Physiol. Mol. Biol. Plants. 28, 637–650. doi:10.1007/s12298-022-01149-9
Soriano, J. M., and Alvaro, F. (2019). Discovering consensus genomic regions in wheat for root-related traits by QTL meta-analysis. Sci. Rep. 9, 10537. doi:10.1038/s41598-019-47038-2
Soriano, J. M., Colasuonno, P., Marcotuli, I., and Gadaleta, A. (2021). Meta-QTL analysis and identification of candidate genes for quality, abiotic and biotic stress in durum wheat. Sci. Rep. 11, 11877. doi:10.1038/s41598-021-91446-2
Swamy, B. P. M., Marathi, B., Ribeiro-Barros, A. I. F., and Ricachenevsky, F. K. (2021). Editorial: development of healthy and nutritious cereals: recent insights on molecular advances in breeding. Front. Genet. 12, 635006. doi:10.3389/fgene.2021.635006
Taranto, F., D’Agostino, N., Catellani, M., Laviano, L., Ronga, D., Milc, J., et al. (2020). Characterization of celiac disease-related epitopes and gluten fractions, and identification of associated loci in durum wheat. Agronomy 10, 1231. doi:10.3390/agronomy10091231
Utebayev, M., Dashkevich, S., Bome, N., Bulatova, K., and Shavrukov, Y. (2019). Genetic diversity of gliadin-coding alleles in bread wheat (Triticum aestivum L.) from Northern Kazakhstan. PeerJ 7, e7082. doi:10.7717/peerj.7082
Vuong, T. D., Sonah, H., Meinhardt, C. G., Deshmukh, R., Kadam, S., Nelson, R. L., et al. (2015). Genetic architecture of cyst nematode resistance revealed by genome-wide association study in soybean. BMC Genomics 16, 593. doi:10.1186/s12864-015-1811-y
Wang, J., Yu, H., Weng, X., Xie, W., Xu, C., Li, X., et al. (2014). An expression quantitative trait loci-guided co-expression analysis for constructing regulatory network using a rice recombinant inbred line population. J. Exp. Bot. 65, 1069–1079. doi:10.1093/jxb/ert464
West, M. A. L., Kim, K., Kliebenstein, D. J., van Leeuwen, H., Michelmore, R. W., Doerge, R. W., et al. (2007). Global eQTL mapping reveals the complex genetic architecture of transcript-level variation in Arabidopsis. Genetics 175, 1441–1450. doi:10.1534/genetics.106.064972
Wu, C. L., Delano, D. L., Mitro, N., Su, S. V., Janes, J., McClurg, P., et al. (2008). Gene set enrichment in eQTL data identifies novel annotations and pathway regulators. PLoS Genet. 4, e1000070. doi:10.1371/journal.pgen.1000070
Wu, P. Y., Yang, M. H., and Kao, C. H. (2021). A statistical framework for QTL hotspot detection. G3 (Bethesda) 11, jkab056. doi:10.1093/g3journal/jkab056
Yang, M. H., Wu, D. H., and Kao, C. H. (2019). A statistical procedure for genome-wide detection of QTL hotspots using public databases with application to rice. G3 (Bethesda) 9, 439–452. doi:10.1534/g3.118.200922
Yang, Y., Xu, C., Shen, Z., and Yan, C. (2022). Crop quality improvement through genome editing strategy. Front. Genome Ed. 3, 819687. doi:10.3389/fgeed.2021.819687
Zaefizadeh, M., Jamaati-e-Somarin, S., Ojaghi, J., Seyedi, S. M., Zabihi-e-Mahmoodabad, R., and Ochi, M. (2010). Genetic diversity for gliadin patterns of durum wheat landraces in the Northwest of Iran and Azerbaijan. Pesq. Agropec. Bras. 45, 1425–1432. doi:10.1590/s0100-204x2010001200013
Zhang, Z., Li, J., Jamshed, M., Shi, Y., Liu, A., Gong, J., et al. (2020). Genome-wide quantitative trait loci reveal the genetic basis of cotton fibre quality and yield-related traits in a Gossypium hirsutum recombinant inbred line population. Plant Biotechnol. J. 18, 239–253. doi:10.1111/pbi.13191
Keywords: candidate genes, grain quality, synteny, Brachypodium, rice, maize
Citation: Marcotuli I, Soriano JM and Gadaleta A (2022) A consensus map for quality traits in durum wheat based on genome-wide association studies and detection of ortho-meta QTL across cereal species. Front. Genet. 13:982418. doi: 10.3389/fgene.2022.982418
Received: 30 June 2022; Accepted: 21 July 2022;
Published: 30 August 2022.
Edited by:
Awais Rasheed, Quaid-i-Azam University, PakistanReviewed by:
Zhiwu Zhang, Washington State University, United StatesLiangliang Gao, Kansas State University, United States
Copyright © 2022 Marcotuli, Soriano and Gadaleta. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Ilaria Marcotuli, ilaria.marcotuli@uniba.it; Jose Miguel Soriano, josemiguel.soriano@irta.cat
†These authors share first authorship