Next Article in Journal
New Insights in Thyroid Cancer and p53 Family Proteins
Next Article in Special Issue
Olive Mill Waste Enhances α-Glucan Content in the Edible Mushroom Pleurotus eryngii
Previous Article in Journal
Identification and Validation of Reference Genes for Seashore Paspalum Response to Abiotic Stresses
Previous Article in Special Issue
Translational Development and Application of (1→3)-β-d-Glucan for Diagnosis and Therapeutic Monitoring of Invasive Mycoses
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 18
Issue 6
10.3390/ijms18061329
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Development of a High-Density SNP-Based Linkage Map and Detection of QTL for β-Glucans, Protein Content, Grain Yield per Spike and Heading Time in Durum Wheat

by
Ilaria Marcotuli
1,
Agata Gadaleta
1,*,
Giacomo Mangini
2,
Antonio Massimo Signorile
2,
Silvana Addolorata Zacheo
2,
Antonio Blanco
2,
Rosanna Simeone
2 and
Pasqualina Colasuonno
1
1
Department of Agricultural and Environmental Science, University of Bari “Aldo Moro”, Via G. Amendola 165/A, 70126 Bari, Italy
2
Department of Soil, Plant and Food Sciences, Section of Genetic and Plant Breeding, University of Bari “Aldo Moro”, Via G. Amendola 165/A, 70126 Bari, Italy
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2017, 18(6), 1329; https://doi.org/10.3390/ijms18061329
Submission received: 11 May 2017 / Revised: 16 June 2017 / Accepted: 17 June 2017 / Published: 21 June 2017
(This article belongs to the Special Issue Glucan: New Perspectives on Biochemistry and Application)
Browse Figures
Versions Notes

Abstract

:
High-density genetic linkage maps of crop species are particularly useful in detecting qualitative and quantitative trait loci for important agronomic traits and in improving the power of classical approaches to identify candidate genes. The aim of this study was to develop a high-density genetic linkage map in a durum wheat recombinant inbred lines population (RIL) derived from two elite wheat cultivars and to identify, characterize and correlate Quantitative Trait Loci (QTL) for β-glucan, protein content, grain yield per spike and heading time. A dense map constructed by genotyping the RIL population with the wheat 90K iSelect array included 5444 single nucleotide polymorphism (SNP) markers distributed in 36 linkage groups. Data for β-glucan and protein content, grain yield per spike and heading time were obtained from replicated trials conducted at two locations in southern Italy. A total of 19 QTL were detected in different chromosome regions. In particular, three QTL for β-glucan content were detected on chromosomes 2A and 2B (two loci); eight QTL controlling grain protein content were detected on chromosomes 1B, 2B, 3B (two loci), 4A, 5A, 7A and 7B; seven QTL for grain yield per spike were identified on chromosomes 1A, 2B, 3A (two loci), 3B (two loci) and 6B; and one marker-trait association was detected on chromosome 2A for heading time. The last was co-located with a β-glucan QTL, and the two QTL appeared to be negatively correlated. A genome scan for genomic regions controlling the traits and SNP annotated sequences identified five putative candidate genes involved in different biosynthesis pathways (β-glucosidase, GLU1a; APETALA2, TaAP2; gigantea 3, TaGI3; 14-3-3 protein, Ta14A; and photoperiod sensitivity, Ppd-A1). This study provides additional information on QTL for important agronomic traits that could be useful for marker-assisted breeding to obtain new genotypes with commercial and nutritional relevance.

Graphical Abstract

1. Introduction

Durum wheat (Triticum turgidum L. subsp. Durum) is the tenth most widely grown crop on the planet with production of more than 36 million metric tons. It is commonly cultivated on both south and north shores of the Mediterranean Sea, with the largest proportion of the global production focused there. This crop was domesticated at the origin of human civilization and Phoenicians traded it along these waters in historical times.
Durum is an excellent dietary source for its high levels of fibre, protein and carotenoids. It is also a significant source of rare minerals like magnesium and selenium. The grain quality of wheat varieties is of paramount importance in determining the nutritional value of the derived end products. The price of durum grain is linked to consumer appreciation and overall dietary quality; better quality is linked with higher prices. Considerable variation exists in grain quality of different varieties and this difference is often due to heritable genetic components, even though the environment and management practices (especially genotype × environment and genotype × agronomic management interactions) also have significant effects [1]. Many studies aimed at improving wheat quality have focused on the grain/flour value of breeding lines, environmental effects, and end-use processing quality [2]. Thus, there are different ways of considering the general words “wheat quality”. The commercial wheat quality value is primarily determined by the milling rate, the quantity of wheat flour derived from 100 kg of seeds. The maximum quantity of flour obtained is affected by seed size, thousand kernel weight, ash mass (or mineral content) and percentage of defective seeds (such as sprouted, small or white seeds, and darkness of grain caused by black point). However, at the technological level, wheat quality is determined by the strength of the flour protein when mixed with water to make dough, and the quantity of water required to produce workable dough. The protein quantity and quality influence the characteristics of the final food products [3]. For example, high protein content brings high water absorbance and increases the shelf life of the products. Another important technological parameter is fibre content, which is linked with many health benefits including those associated with immunomodulatory activity, blood cholesterol, faecal bulking effects, mineral absorption, prebiotic effects, and reduced type II diabetes [4]. The principal components of fibre are β-glucans and arabinoxylans. Even though β-glucans are a minor component of wheat cell walls they are among the most important portion of the fibre-soluble fraction with beneficial effects on human health. Furthermore, the hydrocolloid properties of β-glucan may positively affect bread quality by extending its shelf life in contributing to sustained firmness and prevention of dehydration during storage [5]. Mohamed et al. [6] reported a delay in the hardening of bread containing β-glucan during storage, whereas Kurek et al. [7] found a positive effect of β-glucan on product quality due to increased water retention in the bread. Moreover, they have potential use as functional ingredients to improve mouth-feel of low-fat dairy products [4]. β-glucan and grain protein contents, as well as yield components, are controlled by complex genetic systems that are influenced by environmental factors and management practices. During recent decades, breeding programs have been more focused on quality traits, such as protein and fibre content. Nevertheless, maintenance of yield remains a major objective.
Due to the importance of wheat as a major crop, its genome has been extensively studied. A high-density wheat single nucleotide polymorphism (SNP) iSelect array recently developed by Wang et al. [8] comprises approximately 90,000 gene-associated SNPs and provides a dense coverage of the wheat genome. Many wheat germplasm collections [9,10,11,12] and biparental cross populations [2,13] have also been used to achieve SNP-based maps and SNP consensus maps in both bread wheat and durum [7,14,15], which were also used in quantitative trait studies.
Many studies have reported Quantitative Trait Loci (QTL) for protein content [16], β-glucan content [10] and yield [2], and correlations between protein and yield components in durum wheat, largely investigated in segregating populations, germplasm collections and recurrent selection material [3,17]. Unfortunately, there is a lack of information concerning β-glucan content in relation to protein content and yield components. Therefore, our objectives were to: (a) evaluate the genetic variation in β-glucan and protein content, grain yield per spike and heading time in a recombinant inbred line (RIL) population derived from a cross between two elite durum wheat cultivars; (b) develop a high-density SNP-based genetic linkage map; (c) identify and characterise QTL for β-glucans and protein content, grain yield per spike and heading time; and (d) identify candidate genes associated with the analysed traits.

2. Results

2.1. Field Trait Analysis

Genetic variation in β-glucan content (BG), grain protein content (GPC), grain yield per spike (GYS) and heading time (HT) was assessed in a Duilio × Avonlea (D × A) RIL population grown at Policoro (Matera, Italy) and Valenzano (Bari, Italy) in 2014 (Figure 1).
All traits were determined in three field replicates in order to evaluate the environmental influence. Descriptive statistics (general mean for each parent and RIL, range of the RIL population, genetic variance, and broad-sense heritability) are reported in Table 1.
The phenotypic data for yield and adaptive components highlighted considerable differences between the parents and significant genetic variation in the RIL population in both environments. Avonlea was later heading and had lower GYS than Duilio. Qualitative trait data showed differences in BG, with cv. Avonlea having higher content, whereas GPC was similar for both parents. For all traits, the range of the RIL population was much larger than the range of the parental lines, suggesting transgressive segregation and that favourable alleles were present in both parents.
The means of the RIL population for BG were near Duilio, but with a big range of variability (0.30–0.63% and 0.22–0.68%) in both the environments with a typical pattern of a quantitative trait. The broad-sense heritability estimates (genotype mean basis) ranged from 0.80 to 0.82, indicating the stability of the trait and the phenotypic expression mainly due to genotypic effects.
Protein content in the RIL population ranged from 10.2 to 17.6% with mean values corresponding to the parental values. Broad-sense heritability (0.52 and 0.44) indicated the degree to which phenotypic expression was influenced by genetic versus environmental effects.
Similarly low estimates of heritability were observed for GYS. Moreover, GYS data showed a wide range of variability in the D × A RIL population (0.94–2.64 g and 1.58–3.47 g), typical of a quantitative trait and with mean values corresponding to the parental means. There was wide variation and transgressive segregation for HT with a range of values from 9 to 41. The heritability was 0.99–0.98, indicative of a high degree of stability for this trait. Contrary to the other traits analysed, heading date was bimodally distributed.
Correlation analyses (Table 2) between the various traits indicated significant and positive relationships for GPC and HT in Valenzano 2014 (p ≤ 0.01). A significant and negative correlation was detected between GPC and BG (p ≤ 0.01) in Valenzano 2014. Negative correlations between GPC and GYS were observed at both locations (p ≤ 0.001).

2.2. A Duilio × Avonlea Linkage Map

The D × A RIL population was genotyped with the wheat 90K SNP iSelect assay [7]. Of 81,587 assays, 5577 (6.8%) failed and 69,692 (85.4%) were monomorphic across the mapping population. The remaining 6831 (8.37%) were polymorphic; however, 513 had more than 10% missing data and were excluded from further analyses. Hence, 6318 markers were used for map construction.
Of 6318 loci, 5444 (6.7% from the original SNP dataset) assembled in 36 linkage groups when using a LOD score ≥ 6, and 874 remained unlinked or assembled in small linkage groups (LGs) without a chromosomal anchor locus. Linkage groups were assigned to the A and B genome chromosomes using the physical position of SNP reported in Colasuonno et al. [13] and the SNP anchor loci reported in the durum consensus map of Maccaferri et al. [15] (Figure 2, Table S1).
A total of 2904 (53.3%) markers were localised on the B genome with a total length of 1001.3 cM, whereas 2540 (46.6%) were mapped on the A genome (total length 961.2 cM). The entire map covered 1962.5 cM with an average chromosome length of 140.2 cM. The lengths of individual chromosomes varied from 80.6 cM (chromosome 4B) to 214.7 cM (chromosome 5A). Chromosomes 2A and 7A showed the highest number of markers per cM (4.2 and 4.6, respectively), and chromosome 3B had the lowest marker density (1.6) (Table 3).
Linkage analysis showed a high number of co-segregating SNP (1596). The SNP markers were generally well distributed throughout the genome, although some chromosomes exhibited higher densities.
The overall SNP density was 2.7 markers/cM, with a maximum of 4.6 for chromosome 7B and a minimum of 1.6 for chromosome 3B. The recombinant frequency was 49.7%. The percentage of mapped SNPs showing segregation distortion at p ≥ 0.05 was 5.6%.

2.3. Detection of QTL for β-Glucan Contents, Protein Content, Grain Yield per Spike and Heading Time

Identification of QTL associated with BG, GPC, GYS and HT was carried out based on the genetic map obtained as previously described. Composite interval mapping (CIM) as proposed by Zeng [18] identified three QTL for β-glucans with LOD values ≥3.0 on chromosomes 2A and 2B, consistent in the mean of the two environments and in one location (Policoro), respectively (Figure 2, Table 4). For additive QTL effects, positive and negative signs indicate the relative trait contributions from Avonlea and Duilio, respectively. The QTL on chromosome 2A had a negative effect due to the Duilio alleles, whereas both QTL on chromosome 2B showed positive effects by Avonlea alleles.
Eight QTL for GPC were located on chromosomes 1B, 2B, 3B (two loci), 4A, 5A, 7A and 7B. Positive alleles in Avonlea were on chromosomes 1B, 3B (both loci) and 7B, and those on 2B, 4A and 5A were from Duilio. The 7B QTL (QGpc.mgb-7B.1) was detected in both two environments and the mean.
QTL for GYS were detected on chromosomes 1A, 2B, 3A (two loci), 3B (two loci) and 6B with QGys.mgb-1A.1, QGys.mgb-3A.1 and QGys.mgb-3B.1 having positive effects from Avonlea and a negative effect for the rest of the loci.
Regression analysis allowed identification of a single region on chromosome 2A associated with HT that was consistent in both environments with a high LOD score of 25.5.

2.4. Candidate Genes Related to All the Detected QTL

Sequences of SNP markers in Wang et al. [7] were used to identify potential candidate genes associated with QTL for BG, GPC, GYS and HT detected in the D × A population. All the sequences of the SNP markers located in the QTL regions were analysed for annotated genes or proteins involved in each considered pathway in wheat. For BG, the marker IWB23783 located on chromosome 2B corresponded to β-glucosidase (GLU1a) (Table 5).
For GYS, candidate genes were detected on chromosomes 2B, 3A and 3B. Since the low number of SNPs in the 2B QTL region limited this study, we considered the same region in the consensus durum map [15]. This allowed us to detect the APETALA2 (TaAP2) gene corresponding to the wmc213-2B and wmc243-2B markers. Both SSRs mapped at 4.0 cM from the closest associated (IWA8152) marker. The gigantea 3 (TaGI3) gene showed sequence similarity with the IWB65703 marker (contig_3AS_200_3324175) at 0.6 cM from IWB48828 in the 3A QTL region. The 14-3-3 protein (Ta14A) gene, identified by the SNP marker IWB66165, mapped in the 3B QTL region at 7.0 cM from the closest SNP marker IWB64877 (Table 5).
The heading time gene QHt.mgb-2A.1 appeared to be associated with the photoperiod sensitivity (Ppd-A) gene through SNP marker IWB54033 in the chromosome 2A QTL region.

3. Discussion

The recent emphasis on grain quality over grain quantity has shifted the major objective of wheat breeding programs to increasing nutrient properties while maintaining or increasing grain yield in lines to be released for commercial production. For a long time the major quality objective in durum breeding has been increased GPC [19] and yellow pigment content [9], but recently, more attention has been given to fibre [11], for which the main objective has been improved BG [10]. Different components contribute to grain yield, and among them the fundamental one is grain yield per spike, which is generally negatively correlated to GPC [2]. Another trait correlated with yield is HT, which in cereal crops is determined by vernalisation requirement, photoperiod sensitivity, and narrow-sense earliness or earliness per se [20]. Thanks to the recently released consensus linkage maps for bread wheat [8] and durum [15], development and saturation of genetic maps for detection of QTL underlying agronomic traits has rapidly increased in genetics studies.
In the present work, we evaluated BG, GPC, GYS and HT in an Italian cv. Duilio × Canadian cv. Avonlea RIL population grown in two field trials and developed a high-density linkage map. The total map length was 1962.5 cM, and included 5444 SNP markers, of which 374 SNPs were used to anchor the linkage groups (LGs) to chromosomes. A good coverage and representation of all 14 durum chromosomes was obtained, but the presence of gaps (over 20 cM) led to 36 LGs. Despite the high density of polymorphic markers detected using the 90K SNP assay, linkage groups ranging from 20 to 33 cM have been reported in previous SNP maps [13,21,22,23]. This fragmentation observed for the D × A map could be dependent upon the strict conditions used in determining the marker order (LOD score > 6 and SNP markers with less than 10% missing data). Despite the numerous LGs, our map represented a useful basis for studies on marker-assisted selection and positional cloning. The D × A map had a higher length (1962.5 cM) than previous maps [2,3,13], but a lower marker density (2.7 cM). The average number of mapped markers per chromosome was 389, with a range from 230 SNPs on chromosome 4B to 536 on chromosome 7B. More markers were mapped on the B genome as observed in other linkage studies [2,3,13]. Chromosome 1A was represented by only one LG [15].
A total of 19 QTL for the examined traits were detected, including three for BG, eight for GPC, seven for GYS, and one for HT.
BG QTL were mapped on chromosomes 2A and 2B, as previously reported by Marcotuli et al. [10], but in different chromosome positions. QGbg.mgb-2A.1 was co-located with QHt.mgb-2A.1 controlling HT and with a strong and negative effect on the trait (−9.96) due to Duilio contributing the negative effect. Previous work on yield components reported a major QTL for HT on chromosomes 2BL with negative alleles on the trait [24,25]. We found no significant correlation between BG and HT. Similar results were reported in oat [26,27] and barley [28,29]. Due to the paucity of studies on BG in wheat at any ploidy level, the present data provide some knowledge and understanding of the trait in durum wheat.
A large number of studies have identified major and minor QTL for GPC in wheat [16,30,31]. We mapped eight QTL on chromosomes 1B, 2B, 3B, 4A, 5A, 7A and 7B in the present RIL population. The major stable QTL on chromosome 7B was repeatedly detected in common wheat and durum [16,31,32]. The other QTL on chromosomes 1B [24,33], 3B [30,33,34], 4A [30,31] and 7A [30] were not stable across environments, but were already reported in the literature.
QTL for GYS have been reported on almost all wheat chromosomes [1,33,35,36,37,38,39,40,41,42,43,44]. Most of these studies identified several grain yield QTL; however, the majority were detected in only single environments, and QTL detected in more than one environment showed considerable variation in magnitude of effect [1,36,37,41,42,43]. Our study detected QTL on chromosomes 1A, 2B, 3A, 3B and 6B as already reported in literature [2,45,46,47], of which QGys.mbg.3A.1 was consistent in both environments.
To identify coding sequences in regions associated with QTL, BLAST analyses on the NCBI database [48] were carried out. No gene co-located with the BG QTL QGbg.mgb-2A.1 and QGbg.mgb-2B.1, whereas the mapped SNP IWB23783 associated with QGbg.mgb-2B.2 corresponded to β-glucosidase 1a (GLU1a), a member of the glycoside hydrolase family 1 (GH1). This enzyme catalyses hydrolysis of the glycosidic linkage in glycosides, leading to formation of a hemiacetal or hemiketal sugar and a corresponding free aglycon [49]. We also detected the photoperiod sensitivity (Ppd-A1) gene in the QTL QHt.mgb-2A.1 region. The contig sequence of IWB54033 located in the 2A QTL corresponded to the Ppd-1A sequence, confirming previous results on the relationship between photoperiod response and HT [50,51]. No marker located in gene sequences for GPC was identified, but previous studies reported genes involved in the nitrogen pathway located in the same position of QTL reported in the present work [2,52,53,54,55].
Bioinformatics analysis of markers located in QTL region QGys.mgb-2B.1 for GYS detected no candidate gene associated with the trait. Because previous studies reported a gene (APETALA2) associated with a QTL for grain yield in different species [56,57,58,59], we decided to use all the markers (SNPs and Simple Sequence Repeat, SSR) from the durum wheat consensus map [15] co-locating with all the SNP mapped in the 2B chromosome interval 19.9–37.9 cM of the D × A map. Due to the higher number of considered markers, we identified the APETALA2 (TaAP2) gene associated with the SNP marker IWA8152.
In Arabidopsis, Aukerman and Sakai [56] and Jung et al. [57] demonstrated that overexpression or knockout of the APETALA 2-like (AP2) gene delayed flowering time or caused earlier flowering, respectively. Similar studies were conducted in rice [58], barley [59] and wheat [59]. Another gene involved in flowering promotion and associated with QGys.mgb-3A.2 was gigantea 3 (TaGI3) [60]. A third gene detected in the region controlling GYS was 14-3-3 protein (Ta14A) associated with two SNPs (IWB3714 and IWB11529). Ta14A was strongly expressed in leaf and stem tissues, but was undetectable in roots, suggesting tissue-specific functions [61].

4. Materials and Methods

4.1. Genetic Materials and Field Experiments

A set of 134 F2:7 RILs were developed from a cross between Italian cultivar Duilio and Canadian cultivar Avonlea. After the last generation of selfing, each line was harvested in bulk to provide seeds for replicated trials and DNA extraction. The parental lines and RILs were evaluated for β-glucan (BG), grain protein content (GPC), grain yield per spike (GYS), and heading time (HT) in replicated trials at Policoro (metropolitan city of Matera) and Valenzano (metropolitan city of Bari) in southern Italy in 2014.
A randomised complete block design with three replications was used in field experiments. Plots consisted of 1 m rows spaced 30 cm apart; each plot had 80 germinated seedlings. During the growing season, 10 g of nitrogen per m2 were applied at flowering stage and standard cultivation practices were adopted.
Plots were hand-harvested at maturity and GYS was determined by dividing grain yield by the number of spikes (about 70–80 spikes) per row.
β-glucan (percentage w/w) content in the whole grain was assayed using a Mixed-Linkage β-Glucan Assay Kit (Megazyme International Ireland Ltd, Wicklow, Ireland) based on the method by McCleary and Codd [62] and included the industrial standard for barley (4.1% β-glucan). GPC, expressed as a percentage of protein on a dry weight basis, was determined on a 2 g sample of whole-meal flour using an InfraAlyzer spectrophotometer (Technicon Industrial Systems, Tarrytown, NY, USA) (that is, using near-infrared reflectance spectroscopy).
HT was recorded as the number of days from 1 April 2014, to 50% ear-emergence, corresponding to stage 55 on the Zadoks et al., scale [63].

4.2. SNP Genotyping

Fifty ng/µL of genomic DNA from each RIL and parental line (Duilio and Avonlea) was analysed with the wheat 90K iSelect array [8]. Genotyping was performed at TraitGenetics GmbH (http://www.traitgenetics.de [64]) following the manufacturer’s recommendations as described in Akhunov et al. [65]. The genotyping assays were carried out using the Illumina iScan reader and performed using GenomeStudio software version 2011.1 (Illumina, San Diego, CA, USA).

4.3. Segregation Analysis and Map Construction

Chi-squared tests were used to determine the goodness-of-fit at p > 0.001 of segregation ratios to expected 1:1 ratios for each SNP. All markers with more than 10% missing data, or segregating presence/absence in the mapping population were excluded from further analysis. Linkage analysis between markers and determination of the linear order of loci was performed by JoinMap 4.0 (Kyazma, Wageningen, Netherlands) [66] using the regression mapping algorithm. Grouping was performed using the independence LOD parameter, with groups presenting a LOD 6. The Haldane mapping function was used to calculate map distances [67]. Physically mapped SNP markers [13] and SNP data from the durum consensus maps [15] were used as anchor loci and for assigning linkage groups to specific chromosomes. Linkage groups were named according to the wheat chromosome nomenclature followed by a number linkage groups within chromosomes.

4.4. Statistical Analysis, and QTL and Genes Detections

Analysis of variance was carried out for each trait using GenStat (version 18, VSN International Ltd, Hemel Hempstead, UK,). QTL detection was performed in QGene 8.3.16 using composite interval mapping [18]. A marker-trait association was considered significant when one or more markers was associated with a trait at a threshold of −log10(p) ≥ 3.0, determined by modified Bonferroni correction. For additive effects, positive and negative signs for estimates indicated the contributions of Duilio and Avonlea, respectively, toward higher trait values. Regression analysis implemented in QGene was used for HT. Graphical representations of linkage groups and QTL were carried out using MapChart 2.2 software (SolarWinds, Austin, TX, USA) [68].
Putative candidate genes associated with QTL for BG, GPC, GYS and HT were identified using SNP sequences from Wang et al. [8] blasted against Triticum sequences annotated in NCBI [69] and contig sequences from the Unité de Recherche Génomique Info (URGI) website [70].

5. Conclusions

In the present study, we report a high-density SNP-based genetic map for the Duilio × Avonlea durum wheat RIL population. It comprised 5444 SNP loci. This high-resolution map and the data on β-glucan content, protein content, grain yield per spike, and heading time evaluated in two environments provided increased knowledge of BG in durum wheat and its interaction with other quantitative traits. The three QTL for BG were influenced by environmental differences. The QTL on chromosome 2B was coincident with the β-glucosidase 1a gene involved in the BG pathway. The QTL for BG on chromosome 2A was coincident with HT providing additional information on trait correlations. The marker-trait associations for GPC and GYS confirmed previous data reported in the literature on the polygenic nature of the traits and correspondent QTL in both common and durum wheat. The overall results provide a basis for marker-assisted selection for BG in breeding programs.

Supplementary Materials

Supplementary materials can be found at www.mdpi.com/1422-0067/18/6/1329/s1.

Acknowledgments

The research project was supported by grants from Ministero dell’Istruzione, dell’Università e della Ricerca, project “PON-01_01145 ‘ISCOCEM’ and ‘PRIN 2010–2011’”. The funding organizations had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Author Contributions

Ilaria Marcotuli, Pasqualina Colasuonno, Agata Gadaleta, Antonio Blanco and Rosanna Simeone conceived and designed the experiments; Ilaria Marcotuli, Pasqualina Colasuonno, Giacomo Mangini, Antonio Massimo Signorile and Silvana Addolorata Zacheo performed the experiments; Ilaria Marcotuli and Pasqualina Colasuonno analysed the data; Agata Gadaleta contributed reagents/materials/analysis tools; Ilaria Marcotuli and Pasqualina Colasuonno wrote the paper.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Cobb, J.N.; DeClerck, G.; Greenberg, A.; Clark, R.; McCouch, S. Next-generation phenotyping: Requirements and strategies for enhancing our understanding of genotype–phenotype relationships and its relevance to crop improvement. TAG 2013, 126, 867–887. [Google Scholar] [CrossRef] [PubMed]
  2. Nam, J.H.; Ph, H.Y.W.; Yang, M.H. Achievement and prospect of good quality breeding in winter cereals. Korean J. Crop Sci. 1998, 43, 11–15. [Google Scholar]
  3. Blanco, A.; Mangini, G.; Giancaspro, A.; Giove, S.; Colasuonno, P.; Simeone, R.; Signorile, A.; De Vita, P.; Mastrangelo, A.M.; Cattivelli, L.; et al. Relationships between grain protein content and grain yield components through quantitative trait locus analyses in a recombinant inbred line population derived from two elite durum wheat cultivars. Mol. Breed. 2012, 30, 79–92. [Google Scholar] [CrossRef]
  4. Collins, H.M.; Burton, R.A.; Topping, D.L.; Liao, M.-L.; Bacic, A.; Fincher, G.B. Variability in Fine structures of noncellulosic cell wall polysaccharides from cereal grains: Potential importance in human health and nutrition. Cereal Chem. 2010, 87, 272–282. [Google Scholar] [CrossRef]
  5. Guarda, A.; Rosell, C.M.; Benedito, C.; Galotto, M.J. Different hydrocolloids as bread improvers and antistaling agents. Food Hydrocoll. 2004, 18, 241–247. [Google Scholar] [CrossRef]
  6. Mohamed, A.; Rayas-Duarte, P.; Xu, J. Hard Red Spring wheat/C-TRIM 20 bread: Formulation, processing and texture analysis. Food Chem. 2008, 107, 516–524. [Google Scholar] [CrossRef]
  7. Kurek, M.A.; Wyrwisz, J.; Piwińska, M.; Wierzbicka, A. Influence of the wheat flour extraction degree in the quality of bread made with high proportions of β-glucan. Food Sci. Technol. 2015, 35, 273–278. [Google Scholar] [CrossRef]
  8. Wang, S.W.; Forrest, D.; Allen, K.; Chao, A.; Huang, S.; Huang, B.E.; Maccaferri, M.; Salvi, S.; Milner, S.G.; Cattivelli, L.; et al. Characterization of polyploid wheat genomic diversity using a high-density 90,000 single nucleotide polymorphismarray. Plant Biotechnol. J. 2014, 12, 787–796. [Google Scholar] [CrossRef] [PubMed]
  9. Colasuonno, P.; Lozito, M.L.; Marcotuli, I.; Nigro, D.; Giancaspro, A.; Mangini, G.; de Vita, P.; Mastrangelo, A.M.; Pecchioni, N.; Houston, K.; et al. The carotenoid biosynthetic and catabolic genes in wheat and their association with yellow pigments. BMC Genom. 2017, 18, 122–140. [Google Scholar] [CrossRef] [PubMed]
  10. Marcotuli, I.; Houston, K.; Schwerdt, J.G.; Waugh, R.; Fincher, G.B.; Burton, R.A.; Blanco, A.; Gataleda, A. Genetic diversity and genome wide association study of β-glucan content in tetraploid wheat grains. PLoS ONE 2016, 11, e0152590. [Google Scholar] [CrossRef] [PubMed]
  11. Marcotuli, I.; Houston, K.; Waugh, R.; Fincher, G.B.; Burton, R.A.; Blanco, A.; Gataleda, A. Genome Wide Association Mapping for Arabinoxylan Content in a Collection of Tetraploid Wheats. PLoS ONE 2015, 10, e0132787. [Google Scholar] [CrossRef] [PubMed]
  12. Laidò, G.; Mangini, G.; Taranto, F.; Gadaleta, A.; Blanco, A.; Cattivelli, L.; Marone, D.; Mastrangelo, A.M.; Papa, R.; de Vita, R. Genetic diversity and population structure of tetraploid wheats (L.) estimated by SSR, DArT and pedigree data. PLoS ONE 2013, 8, e67280. [Google Scholar] [CrossRef] [PubMed]
  13. Colasuonno, P.; Gadaleta, A.; Giancaspro, A.; Nigro, D.; Giove, S.; Incerti, O.; Mangini, G.; Signorile, A.; Simeone, R.; Blanco, A. Development of a high-density SNP-based linkage map and detection of yellow pigment content QTLs in durum wheat. Mol. Breed. 2014, 34, 1563–1578. [Google Scholar] [CrossRef]
  14. Cavanagh, C.R.; Chao, S.; Wang, S.; Huang, B.E.; Stephen, S.; Kiani, S.; Forrest, K.; Saintenac, C.; Brown-Guedira, G.L.; Akhunova, A.; et al. Genome-wide comparative diversity uncovers multiple targets of selection for improvement in hexaploid wheat landraces and cultivars. Proc. Natl. Acad. Sci. USA 2014, 110, 8057–8062. [Google Scholar] [CrossRef] [PubMed]
  15. Maccaferri, M.; Ricci, A.; Salvi, S.; Milner, S.G.; Noli, E.; Martelli, P.L.; Casadio, R.; Akhunov, E.; Scalabrin, S.; Vendramin, V.; et al. 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. 2015, 13, 648–663. [Google Scholar] [CrossRef] [PubMed]
  16. Blanco, A.; Simeone, R.; Gadaleta, A. Detection of QTL for grain protein content in durum wheat. Theor. Appl. Genet. 2006, 112, 1195–1204. [Google Scholar] [CrossRef] [PubMed]
  17. Oury, F.X.; Berard, P.; Brancourt-Hulmel, M.; Depatureaux, C.; Doussinault, G.; Galic, N.; Giraud, A.; Heumez, E.; Lecomte, C.; Pluchard, P.; et al. Yield and grain protein concentration in bread wheat: A review and a study of multi-annual data from a French breeding program. J. Genet. Breed. 2003, 57, 59–68. [Google Scholar]
  18. Zeng, Z. Precision mapping of quantitative trait loci. Genetics 1994, 136, 1457–1468. [Google Scholar] [PubMed]
  19. Tabbita, F.; Pearce, S.; Barneix, A.J. Breeding for increased grain protein and micronutrient content in wheat: Ten years of the GPC-B1 gene. J. Cereal Sci. 2017, 73, 183–191. [Google Scholar] [CrossRef]
  20. Kiseleva, A.A.; Shcherban, A.B.; Leonova, I.N.; Frenkel, Z.; Salina, E.A. Identification of new heading date determinants in wheat 5B chromosome. BMC Plant Biol. 2016, 16, 8. [Google Scholar] [CrossRef] [PubMed]
  21. Gadaleta, A.; Giancaspro, A.; Nigro, D.; Giove, S.L.; Incerti, O.; Simeone, R.; Piarulli, L.; Colasuonno, P.; Valè, G.; Cattivelli, L.; et al. A new genetic and deletion map of wheat chromosome 5A to detect candidate genes for quantitative traits. Mol. Breed. 2014, 34, 1599. [Google Scholar] [CrossRef]
  22. Russo, M.A.; Ficco, D.B.M.; Laidò, G.; Marone, D.; Papa, R.; Blanco, A.; Gadaleta, A.; De Vita, P.; Mastrangelo, A.M. A dense durum wheat × T. dicoccum linkage map based on SNP markers for the study of seed morphology. Mol. Breed. 2014, 34, 1579–1597. [Google Scholar] [CrossRef]
  23. Giancaspro, A.; Giove, S.L.; Zito, D.; Blanco, A.; Gadaleta, A. Mapping QTLs for Fusarium head blight resistance in an interspecific wheat population. Front. Plant Sci. 2016, 7, 1–13. [Google Scholar] [CrossRef] [PubMed]
  24. Sourdille, P.; Cadalen, T.; Guyomarc’h, H.; Snape, J.W.; Perretant, M.R.; Charmet, G.; Boeuf, C.; Bernard, S.; Bernard, M. An update of the Courtot × Chinese Spring intervarietal molecular marker linkage map for the QTL detection of agronomic traits in wheat. Theor. Appl. Genet. 2003, 106, 530–538. [Google Scholar] [CrossRef] [PubMed]
  25. Hanocq, E.; Niarquin, M.; Heumez, E.; Rousset, M.; Le Gouis, J. Detection and mapping of QTL for earliness components in a bread wheat recombinant inbred lines population. Theor. Appl. Genet. 2004, 110, 106–115. [Google Scholar] [CrossRef] [PubMed]
  26. Peterson, D.M.; Wesenberg, D.M.; Burrup, D.E. β-glucan content and its relationship to agronomic characteristics in elite oat germplasm. Crop Sci. 1995, 35, 965–970. [Google Scholar] [CrossRef]
  27. Peterson, D.M.; Wesenberg, D.M.; Burrup, D.E.; Erickson, C.A. Relationships among agronomic traits and grain composition in oat genotypes grown in different environments. Crop Sci. 2005, 45, 1249–1255. [Google Scholar] [CrossRef]
  28. Thomas, W.T.B.; Powell, W.; Swanston, J.S. The effects of major genes on quantitatively varying characters in barley. 4. The GPert and denso loci and quality characters. Heredity 1991, 66, 381–389. [Google Scholar] [CrossRef]
  29. Eticha, F.; Grausgruber, H.; Berghoffer, E. Multivariate analysis of agronomic and quality traits of hull-less spring barley (Hordeum vulgare L.). J. Plant Breed. Crop Sci. 2010, 2, 81–95. [Google Scholar]
  30. Prasad, M.; Kumar, N.; Kulwal, P.L.; Roder, M.S.; Balyan, H.S.; Dhaliwal, H.S.; Gupta, P.K. QTL analysis for grain protein content using SSR markers and validation studies using NILs in bread wheat. Theor. Appl. Genet. 2003, 106, 659–667. [Google Scholar] [CrossRef] [PubMed]
  31. Kunert, A.; Naz, A.A.; Dedeck, O.; Pillen, K.; Leon, J. AB-QTL analysis in winter wheat: I. Synthetic hexaploid wheat (T. turgidum ssp. dicoccoides × T. tauschii) as a source of favourable alleles for milling and baking quality traits. Theor. Appl. Genet. 2007, 115, 683–695. [Google Scholar]
  32. Huang, X.Q.; Cloutier, S.; Lycar, L.; Radovanovic, N.; Humphreys, D.G.; Noll, J.S.; Somers, D.J.; Brown, P.D. Molecular detection of QTLs for agronomic and quality traits in a doubled haploid population derived from two Canadian wheats (Triticum aestivum L.). Theor. Appl. Genet. 2006, 113, 753–766. [Google Scholar] [CrossRef] [PubMed]
  33. Kulwal, P.; Kumar, N.; Kumar, A.; Gupta, R.K.; Balyan, H.S.; Gupta, P.K. Gene networks in hexaploid wheat: Interacting quantitative trait loci for grain protein content. Funct. Integr. Genom. 2005, 5, 254–259. [Google Scholar] [CrossRef] [PubMed]
  34. Groos, C.; Robert, N.; Bervas, E.; Charmet, G. Genetic analysis of grain protein content, grain yield and thousand-kernel weight in bread wheat. Theor. Appl. Genet. 2003, 106, 1032–1040. [Google Scholar] [CrossRef] [PubMed]
  35. Bӧrner, A.; Schumann, E.; Furste, A.; Coster, H.; Leithold, B.; Rӧder, M.S.; Weber, W.E. Mapping of quantitative trait loci determining agronomic important characters in hexaploid wheat (Triticum aestivum L.). Theor. Appl. Genet. 2002, 105, 921–936. [Google Scholar]
  36. Huang, X.Q.; Cӧster, H.; Ganal, M.W.; Rӧder, M.S. Advanced backcross QTL analysis for the identification of quantitative trait loci alleles from wild relatives of wheat (Triticum aestivum L.). Theor. Appl. Genet. 2003, 106, 1379–1389. [Google Scholar] [CrossRef] [PubMed]
  37. Huang, X.Q.; Kempf, H.; Ganal, M.W.; Rӧder, M.S. Advanced backcross QTL analysis in progenies derived from a cross between a German elite winter wheat variety and a synthetic wheat (Triticum aestivum L.). Theor. Appl. Genet. 2004, 109, 933–943. [Google Scholar] [CrossRef] [PubMed]
  38. Marza, F.; Bai, G.H.; Carver, B.F.; Zhou, W.C. Quantitative trait loci for yield and related traits in the wheat population Ning7840 × Clark. Theor. Appl. Genet. 2005, 21, 1–11. [Google Scholar] [CrossRef] [PubMed]
  39. McCartney, C.A.; Somers, D.J.; Humphreys, D.G.; Lukow, O.; Ames, N.; Noll, J.; Cloutier, S.; McCallum, B.D. Mapping quantitative trait loci controlling agronomic traits in the spring wheat cross RL4452 × ‘AC Domain’. Genome 2005, 48, 870–883. [Google Scholar] [CrossRef] [PubMed]
  40. Quarrie, S.A.; Steed, A.; Calestani, C.; Semikhodskii, A.; Lebreton, C.; Chinoy, C.; Steele, N.; Pljevljakusić, D.; Waterman, E.; Weyen, J.; et al. A high-density genetic map of hexaploid wheat (Triticum aestivum L.) from the cross Chinese Spring × SQ1 and its use to compare QTLs for grain yield across a range of environments. Theor. Appl. Genet. 2005, 110, 865–880. [Google Scholar] [CrossRef] [PubMed]
  41. Kuchel, H.; Williams, K.J.; Langridge, P.; Eagles, H.A.; Jefferies, S.P. Genetic dissection of grain yield in bread wheat: I. QTL analysis. Theor. Appl. Genet. 2007, 115, 1029–1041. [Google Scholar] [CrossRef] [PubMed]
  42. Kumar, N.; Kulwal, P.L.; Balyan, H.S.; Gupta, P.K. QTL mapping for yield and yield contributing traits in two mapping populations of bread wheat. Mol. Breed. 2007, 19, 163–177. [Google Scholar] [CrossRef]
  43. McIntyre, C.L.; Mathews, K.L.; Rattey, A.; Chapman, S.C.; Drenth, J.; Ghaderi, M.; Reynolds, M.; Shorter, R. Molecular detection of genomic regions associated with grain yield and yield-related components in an elite bread wheat cross evaluated under irrigated and rainfed conditions. Theor. Appl. Genet. 2010, 120, 527–541. [Google Scholar] [CrossRef] [PubMed]
  44. Zheng, B.S.; Le Gouis, J.; Leflon, M.; Rong, W.Y.; Laperche, A.; Bran- court-Hulmel, M. Using probe genotypes to dissect QTL × environment interactions for grain yield components in winter wheat. Theor. Appl. Genet. 2010, 121, 1501–1517. [Google Scholar] [CrossRef] [PubMed]
  45. Munkvold, J.D.; Greene, R.A.; Bermudez-Kandianis, C.E.; La Rota, C.M.; Edwards, H.; Sorrells, S.F.; Dake, T.; Benscher, D.; Kantety, R.; Linkiewicz, A.M.; et al. Group 3 chromosome bin maps of wheat and their relationship to rice chromosome 1. Genetics 2004, 168, 639–650. [Google Scholar] [CrossRef] [PubMed]
  46. Dilbirligi, M.; Erayman, M.; Campbell, B.T.; Randhawa, H.S.; Baenziger, P.S.; Dweikat, I.; Gill, K.S. High-density mapping and comparative analysis of agronomically important traits on wheat chromosome 3A. Genomics 2006, 88, 74–87. [Google Scholar] [CrossRef] [PubMed]
  47. Campbell, B.T.; Baezinger, P.S.; Gill, K.S.; Eskridge, K.M.; Budak, H.; Erayman, M.; Dweikat, I.; Yen, Y. Identification of QTLs and environmental interactions associated with agronomic traits on chromosome 3A of wheat. Crop Sci. 2003, 43, 1493–1505. [Google Scholar] [CrossRef]
  48. Basic Local Alignment Search Tool (BLAST). Available online: https://blast.ncbi.nlm.nih.gov/Blast.cgi (accessed on 20 June 2016).
  49. Lombard, V.; Golaconda, R.H.; Drula, E.; Coutinho, P.M.; Henrissat, B. The Carbohydrate-active enzymes database (CAZy) in 2013. Nucleic Acids Res. 2014, 42, D490–D495. [Google Scholar] [CrossRef] [PubMed]
  50. Beales, J.; Turner, A.; Griffiths, S.; Snape, J.W.; Laurie, D.A. A pseudo-response regulator is misexpressed in the photoperiod insensitive Ppd-D1a mutant of wheat (Triticum aestivum L.). Theor. Appl. Genet. 2007, 115, 721–733. [Google Scholar] [CrossRef] [PubMed]
  51. Scarth, R.; Law, C.N. The control of the day-length response in wheat by the group 2 chromosomes. Z. Pflanzenzüchtung 1984, 92, 140–150. [Google Scholar]
  52. Bernard, S.M.; BlomMøller, A.L.; Dionisio, G.; Kichey, T.; Jahn, T.P.; Dubois, F.; Baudo, M.; Lopes, M.S.; Tercé-Laforgue, T.; Foyer, C.H.; et al. Gene expression, cellular localisation and function of glutamine synthetase isozymes in wheat (Triticum aestivum L.). Plant Mol. Biol. 2008, 67, 89–105. [Google Scholar] [CrossRef] [PubMed]
  53. Fontaine, J.; Ravel, C.; Pageau, K.; Heumez, E.; Dubois, F.; Hirel, B.; Le Gouis, J. A quantitative genetic study for elucidating the contribution of glutamine synthetase; glutamate dehydrogenase and other nitrogen-related physiological traits to the agronomic performance of common wheat. Theor. Appl. Genet. 2009, 119, 645–662. [Google Scholar] [CrossRef] [PubMed]
  54. Gadaleta, A.; Nigro, D.; Giancaspro, A.; Blanco, A. The glutamine synthetase (GS2) genes in relation to grain protein content of durum wheat. Funct. Integr. Genom. 2011, 11, 665–670. [Google Scholar] [CrossRef] [PubMed]
  55. Bordes, J.; Ravel, C.; Jaubertie, J.P.; Duperrier, B.; Gardet, O.; Balfourier, F. Use of a global wheat core collection for association analysis of flour and dough quality traits. J. Cereal Sci. 2011, 54, 137–147. [Google Scholar] [CrossRef]
  56. Aukerman, M.J.; Sakai, H. Regulation of flowering time and floral organ identity by a microRNA and its APETALA2-like target genes. Plant Cell 2003, 15, 2730–2741. [Google Scholar] [CrossRef] [PubMed]
  57. Jung, J.-H.; Seo, Y.-H.; Seo, P.J.; Reyes, J.L.; Yun, J.; Chua, N.-H.; Park, C.M. The GIGANTEA-regulated microRNA172 mediates photoperiodic flowering independent of CONSTANS in Arabidopsis. Plant Cell 2007, 19, 2736–2748. [Google Scholar] [CrossRef] [PubMed]
  58. Chen, Y.C.; Lin, W.S.; Chen, R.K.; Chao, Y.Y.; Chin, S.W.; Chen, F.C.; Lee, C.Y. Establishment of a prediction model for the miRNA-based heading date characteristics of rice in the booting stage. Genet. Mol. Res. 2015, 14, 4381–4390. [Google Scholar] [CrossRef] [PubMed]
  59. Gil-Humanes, J.; Pistón, F.; Martín, A.; Barro, F. Comparative genomic analysis and expression of the APETALA2-like genes from barley, wheat, and barley-wheat amphiploids. BMC Plant Biol. 2009, 9, 66. [Google Scholar] [CrossRef] [PubMed]
  60. Zhao, X.; Zhang, X. Rhythmic expression of wheat TaGI genes introduced by photoperiod and ectopic expression of TaGI1 promoting flowering in transgenic Arabidopsis. Unpublished.
  61. Wang, C.; Ma, Q.H.; Lin, Z.B.; He, P.; Liu, J.Y. Cloning and characterization of a cDNA encoding 14–3-3 protein with leaf and stem-specific expression from wheat. DNA Seq. 2008, 19, 130–136. [Google Scholar] [CrossRef] [PubMed]
  62. McCleary, B.V.; Codd, R. Measurement of (1–3),(1–4)-β-d-glucan in barley and oats: A streamlined enzymatic procedure. J. Sci. Food Agric. 1991, 55, 303–312. [Google Scholar] [CrossRef]
  63. Zadoks, J.C.; Chang, T.T.; Konzak, C.F. A decimal code for the growth stages of cereals. Weed Res. 1974, 14, 415–421. [Google Scholar] [CrossRef]
  64. TraitGenetics. Available online: http://www.traitgenetics.com/en/ (accessed on 20 June 2001).
  65. Akhunov, E.; Nicolet, C.; Dvorak, J. Single nucleotide polymorphism genotyping in polyploid wheat with the Illumina GoldenGate assay. Theor. Appl. Genet. 2009, 119, 507–517. [Google Scholar] [CrossRef] [PubMed]
  66. Van Ooijen, J.W. JoinMap, version 4.0; Software for the Calculation of Genetic Linkage Maps in Experimental Populations; Kyazma B.V.: Wageningen, The Netherlands, 2006.
  67. Haldane, J.B.S. The combination of linkage values, and the calculation of distance between linked factors. J. Genet. 1919, 8, 299–309. [Google Scholar]
  68. Voorrips, R.E. MapChart: Software for the graphical presentation of linkage maps and QTLs. J. Hered. 2002, 93, 77–78. [Google Scholar] [CrossRef] [PubMed]
  69. National Center for Biotechnology Information (NCBI). Available online: http://www.ncbi.nlm.nih.gov/ (accessed on 20 June 1988).
  70. URGI Website. Unité de Recherche Génomique Info. Available online: https://urgi.versailles.inra.fr/ (accessed on 21 September 2015).
Ijms 18 01329 g001 550
Figure 1. Frequency distribution of β-glucan content (BG), grain protein content (GPC), grain yield per spike (GYS) and heading time (HT) in the Duilio × Avonlea recombinant inbred lines (RIL) population grown at Policoro and Valenzano in 2014. Reported values are means of three biological replicates.
Figure 1. Frequency distribution of β-glucan content (BG), grain protein content (GPC), grain yield per spike (GYS) and heading time (HT) in the Duilio × Avonlea recombinant inbred lines (RIL) population grown at Policoro and Valenzano in 2014. Reported values are means of three biological replicates.
Ijms 18 01329 g001
Ijms 18 01329 g002a 550Ijms 18 01329 g002b 550
Figure 2. Schematic representation of the durum wheat chromosomes and Quantitative Trait Loci (QTL) for BG, GPC, GYS and HT detected in 134 RILs derived from Duilio × Avonlea. Map positions are given in cM (on the left). Each chromosome is represented by single nucleotide polymorphism (SNP) markers (on the right) mapped on durum wheat consensus maps [13,15]. QTL are represented by bars on the right of each chromosome bar. QTL names indicate the trait (BG, GPC, GYS and HT) and the environment in which the QTL was detected (Mean14, mean in 2014; POL14, Policoro in 2014; VAL14, Valenzano in 2014); the closest SNP marker is indicated in colour and the candidate gene is reported next to the corresponding marker. Solid bars represent QTL declared at LOD score ≥ 3.0.
Figure 2. Schematic representation of the durum wheat chromosomes and Quantitative Trait Loci (QTL) for BG, GPC, GYS and HT detected in 134 RILs derived from Duilio × Avonlea. Map positions are given in cM (on the left). Each chromosome is represented by single nucleotide polymorphism (SNP) markers (on the right) mapped on durum wheat consensus maps [13,15]. QTL are represented by bars on the right of each chromosome bar. QTL names indicate the trait (BG, GPC, GYS and HT) and the environment in which the QTL was detected (Mean14, mean in 2014; POL14, Policoro in 2014; VAL14, Valenzano in 2014); the closest SNP marker is indicated in colour and the candidate gene is reported next to the corresponding marker. Solid bars represent QTL declared at LOD score ≥ 3.0.
Ijms 18 01329 g002aIjms 18 01329 g002b
Table
Table 1. Descriptive statistics of BG, GYS, GPC and HT for the two parental lines (Avonlea and Duilio) and the RIL population D × A grown at Policoro (Matera) and Valenzano (Bari) in 2014. The genetic variance (s2G) and broad-sense heritability (h2) of the traits are reported for the RIL population.
Table 1. Descriptive statistics of BG, GYS, GPC and HT for the two parental lines (Avonlea and Duilio) and the RIL population D × A grown at Policoro (Matera) and Valenzano (Bari) in 2014. The genetic variance (s2G) and broad-sense heritability (h2) of the traits are reported for the RIL population.
TraitEnvironment
Policoro 2014Valenzano 2014
β-glucan (%)
Duilio0.480.48
Avonlea0.560.59
RIL mean0.470.49
Range0.30–0.630.22–0.68
s2G0.0040.004
h20.800.82
Grain protein content (%)
Duilio14.312.0
Avonlea13.912.1
RIL mean14.6312.7
Range12.6–17.610.2–16.5
s2G0.6110.550
h20.520.44
Grain yield per spike (g)
Duilio1.712.37
Avonlea1.932.55
RIL mean1.672.31
Range0.94–2.641.58–3.47
s2G0.0680.091
h20.470.53
Heading time (days)
Duilio2012
Avonlea3826
RIL mean2919
Range20–419–30
s2G35.16344.044
h20.990.98
Table
Table 2. Correlations between BG, GPC, GYS and HT for the Duilio × Avonlea RIL population grown at Policoro and Valenzano in 2014.
Table 2. Correlations between BG, GPC, GYS and HT for the Duilio × Avonlea RIL population grown at Policoro and Valenzano in 2014.
TraitEnvironmentGPCGYSHT
BGPolicoro−0.040.12−0.10
Valenzano−0.24 **−0.06−0.02
GPCPolicoro −0.36 ***−0.11
Valenzano −0.38 ***0.20 **
GYSPolicoro 0.16
Valenzano 0.05
** significant at p = 0.01; *** significant at p = 0.001.
Table
Table 3. Number and distribution of SNP markers in the Duilio × Avonlea map.
Table 3. Number and distribution of SNP markers in the Duilio × Avonlea map.
ChromosomeNumber of Linkage GroupsNumber of SNPsChromosome Length (cM)Number of SNPs/cM
1A1302149.22.0
1B3522137.93.8
2A3422101.44.2
2B4492169.62.9
3A3283128.52.2
3B2294180.61.6
4A2273116.82.3
4B323080.62.8
5A2473214.72.2
5B4454151.23.0
6A2287143.72.0
6B2376123.53.0
7A4500106.94.6
7B2536157.93.4
Total genome A172540961.22.6
Total genome B1929041001.32.9
Total genomes AB3654441962.52.7
Table
Table 4. Marker-trait associations for BG, GPC, GYS, and HT determined for each location (Policoro and Valenzano) and for the mean values of the two environments. QTL, closest marker, chromosome locations, interval (cM), marker effects (additive effect with positive and negative value indicating contribution of parent Avonlea and Duilio, respectively), LOD value (logarithm-of-odds) ≥3.0, and R2 (% of phenotypic variation explained by QTL) are reported for each SNP associated with each trait.
Table 4. Marker-trait associations for BG, GPC, GYS, and HT determined for each location (Policoro and Valenzano) and for the mean values of the two environments. QTL, closest marker, chromosome locations, interval (cM), marker effects (additive effect with positive and negative value indicating contribution of parent Avonlea and Duilio, respectively), LOD value (logarithm-of-odds) ≥3.0, and R2 (% of phenotypic variation explained by QTL) are reported for each SNP associated with each trait.
QTLClosest MarkerChromQTL IntervalMean Across EnvironmentsPolicoroValenzano
cMEffectLODR2EffectLODR2EffectLODR2
β-glucan content
QGbg.mgb-2A.1IWB12802A(35.8–48.0)−0.034.50.14−0.026.10.19---
QGbg.mgb-2B.1IWB301152B(0.1–3.9)0.024.70.15---0.036.10.19
QGbg.mgb-2B.2IWB237832B(29.9–47.9)0.023.80.120.033.10.10---
Protein content
QGpc.mgb-1B.1IWB419241B(81.9–84.5)0.643.00.10------
QGpc.mgb-2B.1IWA5442B(0.1–5.9)---−0.534.10.13---
QGpc.mgb-3B.1IWB668423B(11.7–14.0)---0.224.80.15---
QGpc.mgb-3B.2IWB718423B(126.0–130.9)0.183.60.12------
QGpc.mgb-4A.1IWB711804A(57.1–62.0)−0.193.70.12------
QGpc.mgb-5A.1IWB283505A(63.4–80.5)−0.434.30.14−0.385.40.17---
QGpc.mgb-7A.1IWB203817A(1.9–15.9)---−0.315.60.17---
QGpc.mgb-7B.1IWB714997B(53.9–61.9)0.273.60.120.566.30.290.363.30.12
Grain yield per spike
QGys.mgb-1A.1IWB476511A(8.3–16.0)0.157.10.22---0.249.50.28
QGys.mgb-2B.1IWA81522B(19.9–37.9)---0.073.00.10---
QGys.mgb-3A.1IWB696013A(5.5–12.5)0.099.30.280.105.30.170.094.70.15
QGys.mgb-3A.2IWB488283A(42.0–50.0)−0.064.50.14−0.128.10.24---
QGys.mgb-3B.1IWB648773B(20.6–28.5)0.083.60.120.134.10.13---
QGys.mgb-3B.2IWB254953B(72.0–84.0)−0.064.10.13------
QGys.mgb-6B.1IWB236596B(26.0–32.1)−0.087.00.21---−0.095.50.17
Heading time
QHt.mgb-2A.1IWB540332A(27.4–49.3)−9.9625.50.58−10.4826.00.59−9.3123.80.56
Table
Table 5. Putative candidate genes associated with QTL for BG, GYS and HT detected in Duilio × Avonlea RIL population. Gene name, enzyme, contig, marker name and ID, type of mutation and localization on D × A consensus map and the durum consensus map (Maccaferri et al., [15]) are listed for each associated gene trait.
Table 5. Putative candidate genes associated with QTL for BG, GYS and HT detected in Duilio × Avonlea RIL population. Gene name, enzyme, contig, marker name and ID, type of mutation and localization on D × A consensus map and the durum consensus map (Maccaferri et al., [15]) are listed for each associated gene trait.
TraitGeneEnzymeContigMarker IDSNPWheat Map Position
ChromD × A MapDurum Map
β-glucan contentGLU1aβ-glucosidase 1acontig_2BL_8006334IWB23783T/C2B42.345.9
Grain yield per spikeTaAP2APETALA2contig_2BS_5227696wmc213-2B-2B-31
contig_2BS_5208541wmc243-2B-2B-31.1
TaGI3gigantea 3contitg_3AS_200_3324175IWB65703T/C3A50.651.5
Ta14A14-3-3 proteincontig_3B_10411870IWB3714A/G3B31.533.2
contig_3B_10411871IWB11529A/G3B31.633.5
Heading timePpd-A1Photoperiod sensitivitycontig_2AS_5262553IWB54033A/C2A47.846.2

Share and Cite

MDPI and ACS Style

Marcotuli, I.; Gadaleta, A.; Mangini, G.; Signorile, A.M.; Zacheo, S.A.; Blanco, A.; Simeone, R.; Colasuonno, P. Development of a High-Density SNP-Based Linkage Map and Detection of QTL for β-Glucans, Protein Content, Grain Yield per Spike and Heading Time in Durum Wheat. Int. J. Mol. Sci. 2017, 18, 1329. https://doi.org/10.3390/ijms18061329

AMA Style

Marcotuli I, Gadaleta A, Mangini G, Signorile AM, Zacheo SA, Blanco A, Simeone R, Colasuonno P. Development of a High-Density SNP-Based Linkage Map and Detection of QTL for β-Glucans, Protein Content, Grain Yield per Spike and Heading Time in Durum Wheat. International Journal of Molecular Sciences. 2017; 18(6):1329. https://doi.org/10.3390/ijms18061329

Chicago/Turabian Style

Marcotuli, Ilaria, Agata Gadaleta, Giacomo Mangini, Antonio Massimo Signorile, Silvana Addolorata Zacheo, Antonio Blanco, Rosanna Simeone, and Pasqualina Colasuonno. 2017. "Development of a High-Density SNP-Based Linkage Map and Detection of QTL for β-Glucans, Protein Content, Grain Yield per Spike and Heading Time in Durum Wheat" International Journal of Molecular Sciences 18, no. 6: 1329. https://doi.org/10.3390/ijms18061329

APA Style

Marcotuli, I., Gadaleta, A., Mangini, G., Signorile, A. M., Zacheo, S. A., Blanco, A., Simeone, R., & Colasuonno, P. (2017). Development of a High-Density SNP-Based Linkage Map and Detection of QTL for β-Glucans, Protein Content, Grain Yield per Spike and Heading Time in Durum Wheat. International Journal of Molecular Sciences, 18(6), 1329. https://doi.org/10.3390/ijms18061329

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop