Original Research ARTICLE
Fine Mapping of QTLs for Ascochyta Blight Resistance in Pea Using Heterogeneous Inbred Families
- Crop Development Centre – Department of Plant Sciences, University of Saskatchewan, Saskatoon, SK, Canada
Ascochyta blight (AB) is an important disease of pea which can cause severe grain yield loss under wet conditions. In our previous study, we identified two quantitative trait loci (QTLs) abIII-1 and abI-IV-2 for AB resistance and these QTLs were consistent across locations and/or years in an inter-specific pea population (PR-19) developed from a cross between Alfetta (Pisum sativum) and P651 (P. fulvum). The objectives of this study were to fine map the abIII-1 and abI-IV-2 QTLs using a high density single nucleotide polymorphism (SNP)-based genetic linkage map and analyze identified markers in heterogeneous inbred family (HIF) populations. Selective genotyping of 51 PR-19 recombinant inbred lines was performed using genotyping-by-sequencing (GBS) and the resulting high density genetic linkage map was used to identify eight new SNP markers within the abI-IV-2 QTL, whereas no additional SNPs were identified within the abIII-1 QTL. Two HIF populations HIF-224 (143 lines) and HIF-173 (126 lines) were developed from F6 RILs PR-19-224 and PR-19-173, respectively. The HIF populations evaluated under field conditions in 2015 and 2016 showed a wide range of variation for reaction to AB resistance. Lodging score had significant positive (P < 0.001) correlation with AB scores. HIFs were genotyped using SNP markers within targeted QTLs. The genotypic and phenotypic data of the HIFs were used to identify two new QTLs, abI-IV-2.1 and abI-IV-2.2 for AB resistance within the abI-IV-2 QTL. These QTLs individually explained 5.5 to 14% of the total phenotypic variation. Resistance to lodging was also associated with these two QTLs. Identified SNP markers will be useful in marker assisted selection for development of pea cultivars with improved AB resistance.
Ascochyta blight (AB), caused by Peyronellaea pinodes (Berk. & A. Bloxam) Aveskamp, Gruyter & Verkley (Aveskamp et al., 2010), is the most important pea (Pisum sativum) disease which can severely affect grain yield under wet conditions in most pea growing regions in the world (Lawyer, 1984; Xue et al., 1997; Kraft et al., 1998). The impact of the disease under field conditions is greatly affected by agronomic traits including lodging and plant height (Tar’an et al., 2003; Banniza et al., 2005; Le May et al., 2009; Jha et al., 2013, 2016). Genetic resistance is the optimal approach to reduce the disease impact (Zimmer and Sabourin, 1986). More than 3500 cultivated pea accessions were evaluated for their reaction to the disease resulting in the identification of a few lines with low to moderate levels of resistance (Kraft et al., 1998; Zhang et al., 2006). In contrast, a higher level of resistance was identified in wild pea (P. fulvum) accessions (Clulow et al., 1991; Wroth, 1998; Fondevilla et al., 2005; Jha et al., 2012). Further, Fondevilla et al. (2005) reported the highest level of resistance in accession P651 (P. fulvum) compared to other wild peas, P670 (P. sativum ssp. elatius) and P665 (P. sativum ssp. syriacum). Promising accessions (P. fulvum and P. sativum ssp. elatius) were identified upon evaluation of 44 wild pea accessions which had the potential for improvement of AB resistance (Jha et al., 2012). Among them, the most promising accession, P651 (P. fulvum) was utilized for resistance breeding (Sindhu et al., 2014; Jha et al., 2016).
Previously, more than 30 quantitative trait loci (QTLs) were identified for resistance to AB in P. sativum mapping populations on all seven linkage groups (LGs) under field or controlled conditions (Timmerman-Vaughan et al., 2002, 2004; Tar’an et al., 2003; Prioul et al., 2004). QTLs were also identified in a cross involving wild pea, P. sativum subsp. syriacum (Fondevilla et al., 2008, 2011; Carrillo et al., 2014). Co-localization of QTLs for disease resistance with candidate genes including RGAs (resistance gene analogs), PsDof1 (a putative transcription factor) and DRR230-b (a pea defensin) involved in defense responses to P. pinodes was reported in pea (Timmerman-Vaughan et al., 2002, 2016; Prioul-Gervais et al., 2007). Further, Jha et al. (2015) reported significant association of SNPs detected within candidate genes PsDof1 (PsDof1p308) and RGA-G3A (RGA-G3Ap103) with AB scores. Most recently, nine QTLs were identified for AB resistance in an inter-specific pea population (PR-19) developed from a cross between Alfetta (P. sativum) and wild pea accession P651 (P. fulvum) (Jha et al., 2016). These QTLs individually explained 7.5 to 28% of the phenotypic variation.
Quantitative trait loci mapping studies in several pea crosses have resulted in the identification of genomic regions associated with AB resistance, however, these QTLs cover large regions which may not be effective for marker-assisted selection (MAS). Though several markers linked to resistance genes have been identified, even the closest markers are not necessarily tightly linked to the gene of interest (reviewed by Michelmore, 1995). Recombination could occur between a marker and QTL if markers are not tightly linked to genes (Collard et al., 2005). High-resolution or fine mapping of QTLs can be used to identify more tightly-linked or perfect markers within the gene sequence that can be efficiently utilized for MAS (reviewed by Mohan et al., 1997). Development of an advanced population, such as near isogenic lines (NILs), is required for fine mapping. Conventional consecutive backcrossing method was the original method for NIL development. Tuinstra et al. (1997) proposed development of heterogeneous inbred family (HIF) populations, an alternative, more efficient method than the NILs. This approach has been widely used in several species including Arabidopsis, soybean and maize for fine mapping of QTLs (Meng et al., 2008; Bai et al., 2010; Todesco et al., 2010; Coles et al., 2011; Dwiyanti et al., 2011; Watanabe et al., 2011; Bouteillé et al., 2012).
Among the nine AB resistance QTLs identified in PR-19 population, two QTLs abIII-1 and abI-IV-2 were consistent across locations and/or years (Jha et al., 2016). The objectives of this research were to identify additional SNP markers within abIII-I and abI-IV-2 QTLs and to fine map them using HIF populations for identification of closely linked markers for AB resistance in pea.
Materials and Methods
Previously, PR-19 recombinant inbred line (RIL) population was generated from a cross between Alfetta (P. sativum) and P651 (P. fulvum) (Sindhu et al., 2014). P651 (original code IFPI3232) was first identified in Syria, then characterized by Consejo Superior de Investigaciones Científicas (Cordoba, Spain). For fine mapping of QTLs abI-IV-2 and abIII-1, HIF populations HIF-224 and HIF-173 were developed from F6 RILs of PR-19-224 and PR-19-173, respectively.
Selection of PR-19 Lines for HIF Populations
RILs PR-19-57, PR-19-132, PR-19-176, and PR-19-224 segregated for marker loci associated with the QTL abI-IV-2, and PR-19-04, PR-19-65, PR-19-115, and PR-19-173 segregated for marker loci associated with the QTL abIII-1. Three seeds of each of these RILs were sown in 2 gallon pots in a greenhouse with 22 + 3°C day/20 + 3°C night temperature under an 18-h photoperiod with approximately 60% relative humidity. Genomic DNA was extracted from freeze-dried leaf tissue collected from each plant using DNeasy Plant Mini Kit (QIAGEN Inc., Valencia, CA, USA) and used for Kompetitive Alelle Specific PCR (KASP) assays to validate heterozygous alleles for SNP loci within the QTLs. Allele-specific primers were designed for SNP loci PsC8780p118 (abIII-1) and PsC6805p316 (abI-IV-2) (Supplementary Table S1) using Primer-Picker software (LGC Genomics, Beverly, MA, USA). A total reaction volume of 10 μl was prepared by adding 20 ng of template DNA, 5 μl of KASP 2X Reaction Mix and 0.14 μl of KASP assay mixture (LGC Genomics, Beverly, MA, USA) in a 96-well plate format. Amplifications were performed using StepOnePlus Real-Time PCR system (Applied Biosystems, USA) according to the program described in Jha et al. (2015). Genotypic data were analyzed using SNPViewer software (LGC Genomics, Beverly, MA, USA).
Development of HIF-224 and HIF-173
Ten F6 seeds each for PR-19-224 and PR-19-173 were grown under greenhouse conditions and tested for heterogeneity by KASP assays as described earlier. Based on these assays, five seeds for PR-19-224 and seven seeds for PR-19-173 had heterozygous alleles for markers associated with abI-IV-2 and abIII-1, respectively. Seeds were bulked from five plants of PR-19-224 and seven plants of PR-19-173. Using single seed descent, self-pollination and bulking of seeds were conducted for F7 to F8 generation. Progenies at F8 were represented HIF-224 and HIF-173 for PR-19-224 and PR-19-173, respectively (Figure 1).
FIGURE 1. Development of heterogeneous inbred family (HIF) populations, HIF-224 and HIF-173 from lines PR-19-224 and PR-19-173, respectively.
Assessment of AB Resistance and Other Agronomic Traits Under Field Conditions
HIF-224 (143 lines) and HIF-173 (126 lines) along with parental checks (Alfetta and P651) were evaluated for reaction to AB and other agronomic traits including days to flower (DTF), plant height, lodging, days to maturity (DTM), and grain yield on a plot basis under field conditions in 2015 at Saskatoon with two replicates, and in 2016 at Saskatoon and Rosthern with three replicates at each location. The experimental design was a randomized complete block design with three-row plots of 1.0 m × 1.0 m, a plant density of 75 plants m-2 and row spacing of 0.25 m. Plants were inoculated at the start of the flowering stage with approximately 3 g per plot of pea straw that had been naturally infected by P. pinodes in the previous season, air dried, and chopped into approximately 2-cm pieces. HIFs were evaluated for AB severity at pod filling and physiological maturity stages (80% of pods in the plot turned brown) using a scale of 0 (no disease) to 9 (whole plant severely blighted) based on Xue et al. (1996). Lodging was assessed on a 1 (upright) to 9 (completely lodged) scale. Plant height was measured from the soil level to the tip of the central stem at physiological maturity. DTF and DTM were calculated as the number of days from planting to 50% bloom and physiological maturity, respectively.
Identification of Additional SNPs in QTLs
A high density genetic linkage map of PR-19 based on selective genotyping of the RIL population was developed for identification of additional SNP loci within the two targeted QTLs. Fifty-one F7 RILs of PR-19 including PR-19-224 and PR-19-173 along with the parents (Alfetta and P651) were genotyped using genotyping-by-sequencing (GBS) method as described by Elshire et al. (2011). Twenty ng/μL DNA of each RIL as quantified using picogreen was digested with restriction enzymes PstI and MspI. Digested DNA of individual RILs was ligated with a unique 4 to 8 base pair barcode adapter. At this stage the DNA samples were pooled for construction of a single library for sequencing. Paired-end sequencing of the library was done in a single lane of an Illumina Hiseq sequencer using V4 sequencing chemistry.
The raw sequence reads were assigned to individual RILs based on the ligated barcode adapter. Following this deconvolution, barcode sequences were removed from the sequence. The reads were then trimmed for quality with Trimmomatic-0.33, and mapped to the draft genome assembly provided through the pea genome sequencing consortium (Madoui et al., 2016) using Bowtie2-2.2.5. SNP variants were identified and converted to VCF format using Samtools-1.1 and BCFtools-1.1.
After filtering for missing values and heterozygosity, 6160 SNP markers were selected for linkage analysis. Segregation data of these markers were combined with 733 polymorphic SNP markers previously genotyped using Illumina GoldenGate 1536 SNP array (Jha et al., 2016). Combined SNP marker segregation data were used for linkage analysis using MstMap. SNP markers from the GoldenGate assay served as anchor markers to identify additional SNP loci within the targeted QTLs. All the SNP markers identified within QTLs were converted to KASP assays (Supplementary Table S1) and used for genotyping of the complete set of 144 RILs of PR-19 for cross-validation of their genetic linkage positions.
Genotyping of HIF-224 and HIF-173
Genomic DNA was extracted from freeze dried leaf tissues collected from single plants of HIF-224 (143 lines) and HIF-173 (126 lines) progenies using DNeasy Plant Mini Kit (QIAGEN Inc., Valencia, California, USA). HIF-224 lines were genotyped using 20 SNP markers (Supplementary Table S1) representing the QTL abI-IV-2 and the region adjacent to the QTL by KASP assays. HIF-173 population segregating for QTL abIII-1 was genotyped with three SNP markers, PsC22609p103, PsC8780p118, and PsC23317p284, each representing a unique locus within this QTL.
Linkage Mapping and QTL Analysis in HIF Populations
The linkage map was constructed separately for PR-19, HIF-224, and HIF-173 using MAPMAKER (Lander et al., 1987). QTL mapping was performed by composite interval mapping (CIM) using Windows QTL Cartographer 2.5 (Wang et al., 2012). The significance threshold (P < 0.05) was used to declare the presence of QTLs by performing 1000 permutations of the data (Churchill and Doerge, 1994). MapChart 2.2 was used for graphical presentation of linkage maps (Voorrips, 2002).
PROC MIXED implemented in SAS®9.3 (SAS Institute Inc. Cary, NC, USA) was used for data analysis. Line was treated as a fixed effect whereas replication was treated as a random effect across the HIFs. Homogeneity of variance test (HOVTEST) was used to assess the homogeneity of variance among replications.
Selection of PR-19 Lines for HIF Populations
Four RILs each tested for abI-IV-2 (PR-19-57, PR-19-132, PR-19-176, and PR-19-224) and abIII-1 (PR-19-04, PR-19-65, PR-19-115, and PR-19-173) had heterozygous alleles within QTLs, i.e., these lines were segregating for markers associated with AB, which is a prerequisite for HIF development. On the basis of KASP assays and AB scores of lines, PR-19-224 and PR-19-173 were selected for development of HIF-224 and HIF-173, respectively.
Assessment of AB Resistance and Other Agronomic Traits under Field Conditions
HIF-224 and HIF-173 showed a wide range of variation for reaction to AB, plant height, lodging, and grain yield under field conditions in 2015 at Saskatoon and in 2015 and 2016 at Saskatoon and Rosthern locations in Saskatchewan (Tables 1, 2 and Figures 2–4). Data from different station years could not be combined for analysis of variance due to significant effect of locations and years in the HOVTEST. In general, the effect of line was significant (P < 0.05) for AB scores, plant height, lodging, and grain yield. AB scores of HIF-224 ranged from 2 to 7 at pod filling, and 2 to 8 at physiological maturity (0–9 scale), whereas for HIF-173, scores ranged from 1 to 7 at pod filling, and 2 to 9 at physiological maturity. Alfetta had disease score of 3 to 4 at pod filling and 4 to 5 at physiological stage, whereas P651 had disease score of 2 to 3 and 3 to 4 at pod filling and physiological maturity stage, respectively. Lodging scores varied from 1 to 9 for HIF-224, whereas for HIF-173, scores varied from 1 to 7 on the 1–9 scale. Alfetta had 1 to 3 lodging score whereas P651 had 8 to 9 score. HIFs had a small range of variation for DTF and DTM at different station years, while plant height and grain yield had a wide range of variation among tested HIF lines. For both HIFs, AB scores were positively correlated with lodging (P < 0.001) and negatively correlated with plant height (P < 0.001) and grain yield (P < 0.01) (Tables 3, 4).
TABLE 1. F-values, coefficients of variations (CV) of statistical analyses, and means with standard deviations (SD) of ascochyta blight (AB) scores and other agronomic assessments for 143 lines of heterogeneous inbred family (HIF) population, HIF-224 evaluated under field conditions in 2015 at Saskatoon and in 2016 at Saskatoon and Rosthern, Saskatchewan.
TABLE 2. F-values, coefficients of variations (CV) of statistical analyses, and means with standard deviations (SD) of AB scores and other agronomic assessments for 126 lines of heterogeneous inbred family (HIF) population, HIF-173 evaluated under field conditions in 2015 at Saskatoon and in 2016 at Saskatoon and Rosthern, Saskatchewan.
FIGURE 2. Frequency distribution of 143 lines of HIF population, HIF-224 using least square means of Saskatoon 2015, Saskatoon 2016 and Rosthern 2016 for the reaction to ascochyta blight (AB) resistance at pod filling and physiological maturity stages under field conditions.
FIGURE 3. Frequency distribution of 126 lines of HIF population, HIF-173 using least square means of Saskatoon 2015, Saskatoon 2016 and Rosthern 2016 for the reaction to AB resistance at pod filling and physiological maturity stages under field conditions.
FIGURE 4. Frequency distribution of 143 lines of HIF population, HIF-224 (A) and 126 lines of HIF-173 (B) using least square means of Saskatoon 2015, Saskatoon 2016 and Rosthern 2016 for lodging at physiological maturity stage under field conditions.
TABLE 3. Pearson correlation coefficients for traits of 143 lines of heterogeneous inbred family (HIF) population, HIF-224 evaluated under field conditions in 2015 at Saskatoon and in 2016 at Saskatoon and Rosthern, Saskatchewan.
TABLE 4. Pearson correlation coefficients for traits of 126 lines of heterogeneous inbred family (HIF) population, HIF-173 evaluated under field conditions in 2015 at Saskatoon and in 2016 at Saskatoon and Rosthern, Saskatchewan.
Identification of Additional SNP Markers within QTLs
Overall, 10,985 SNPs were identified at a read depth of 10 by selective genotyping of 51 PR-19 RILs using GBS method. After filtering for allele distribution, 6160 SNPs along with 733 previously genotyped SNPs were used for construction of a high density genetic linkage map to identify markers within QTLs. Based on the high density genetic linkage map, 12 SNP markers were identified within abI-IV-2 QTL. Ten of the 12 markers along with previously identified SNP markers from an Illumina GoldenGate array were used for genotyping of a complete set of PR-19 RILs (144) to re-confirm their position and order within the QTL. Following linkage analysis of markers of this QTL, the eight SNP markers identified from the high density genetic linkage map were confirmed to localize within the existing QTL flanked by SNP markers PsC943p541/PsC4233p498 and PsC8970p349/PsC7884p449, whereas two SNPs were located to the region adjacent to the QTL. Mapping of eight additional SNPs within the QTL has increased the map distance of the QTL from 13.4 to 17.1 cM (Figures 5A,B).
FIGURE 5. Quantitative trait locus (QTL) abI-IV-2 based on SNP linkage map in PR-19 (Jha et al., 2016) (A), additional SNPs identified by fine mapping using genotyping-by-sequencing (GBS) method in PR-19 (B). Locations of QTLs for AB are shown by vertical bars. S14 and R14 associated with QTLs name denote 2014 Saskatoon and 2014 Rosthern, respectively.
In the case of abIII-1, based on the high density genetic linkage map, no additional SNP marker was identified within the QTL (Figure 6B). Two flanking markers of the QTL were converted to KASP assays and were used for genotyping the complete set of RILs. Linkage analysis of this region based on these two flanking markers and known existing markers within the QTL reconfirmed the order of SNP markers on the high density genetic linkage map.
FIGURE 6. Quantitative trait locus abIII-1 based on SNP linkage map in PR-19 (Jha et al., 2016) (A), additional SNPs identified by fine mapping using GBS method in PR-19 (B). Locations of QTLs for AB are shown by vertical bars. S13, S14, and R14 associated with QTLs name denote 2013 Saskatoon, 2014 Saskatoon, and 2014 Rosthern, respectively.
Fine Mapping of QTLs for AB Resistance
For fine mapping of abI-IV-2 QTL, 143 lines of HIF-224 segregating for this QTL were genotyped with 20 SNP markers using KASP assays. This set of 20 SNP markers included 10 previously known SNP markers and 10 markers currently identified through the high density genetic linkage map. Of the total genotyped, 17 SNP markers were used for linkage analysis to verify the marker order and distance in the HIF population. The 17 SNP markers represented a map distance of 86.3 cM in HIF-224 population (Figure 7). Based on the field evaluation of HIF-224 population in 2015 and 2016 trials, two new QTLs, abI-IV-2.1 and abI-IV-2.2 were identified for AB resistance within the abI-IV-2 QTL (Table 5 and Figure 7). QTL abI-IV-2.1 explained 5.5 to 14% of the total phenotypic variation, whereas abI-IV-2.2 explained 7 to 10% of the total variation. QTLs for lodging resistance were also associated with these two QTLs. Alfetta contributed alleles for AB resistance as well as for lodging resistance. Fine mapping with HIF lines has confirmed the occurrence of AB resistance QTLs within the previously reported QTL ab-IV-2, and provided additional markers for MAS of this QTL in breeding populations.
FIGURE 7. Quantitative trait locus abI-IV-2 based on GoldenGate and GBS markers in HIF population, HIF-224. Locations of QTLs for AB and lodging (lodg) are shown by vertical bars. S15, S16, and R16 associated with QTLs name denote 2015 Saskatoon, 2016 Saskatoon, and 2016 Rosthern, respectively.
TABLE 5. Quantitative trait loci (QTLs) detected for reaction to AB resistance and lodging in abI-IV-2 QTL in HIF population, HIF-224 evaluated under field conditions in 2015 at Saskatoon and in 2016 at Saskatoon and Rosthern.
Additional SNP makers within the abIII-1 QTL were not identified using the high density genetic linkage map. The extreme distortion of allele segregation determined based on the existing three SNP markers within this QTL did not allow for the determination of the linkage order of these markers in HIF-173 population. Additionally, significant recombination within this QTL was not identified in the HIF family to continue with other tests to determine the significance of these markers.
ABs are the most important diseases of pulse crops (Tivoli et al., 2006; Muehlbauer and Chen, 2007). Resistance breeding is considered the most effective method of control; however, this process is slow due to the complex nature of resistance (Muehlbauer and Chen, 2007; Rubiales and Fondevilla, 2012). Significant progress has been made in resistance breeding with the advancement of innovative tools including next generation sequencing. Several QTLs have been reported for AB resistance in pea (Timmerman-Vaughan et al., 2002, 2004; Tar’an et al., 2003; Prioul et al., 2004; Fondevilla et al., 2008, 2011; Jha et al., 2016), chickpea (Udupa and Baum, 2003; Lichtenzveig et al., 2006; Tar’an et al., 2007; Sabbavarapu et al., 2013), lentil (Sudheesh et al., 2016), and faba bean (Atienza et al., 2016).
Sudheesh et al. (2016) reported validation of previously reported QTLs for AB resistance in lentil on genetic maps based on SNP and SSR markers developed from three RIL populations. Further, they identified two common genomic regions for disease resistance in two out of three maps that could provide validated markers associated with disease for lentil improvement. Similarly, Atienza et al. (2016) studied validation and stability of major QTLs located on chromosomes II and III for AB resistance in faba bean under field and controlled conditions and reported that QTL Af2 located on chromosome II was the same QTL reported previously by other researchers. In chickpea, QTLs were identified for AB resistance on LGs 2, 3, 4, 5, 6, and 8 (Udupa and Baum, 2003; Lichtenzveig et al., 2006; Tar’an et al., 2007; Anbessa et al., 2009; Sabbavarapu et al., 2013). Among them, one major QTL on LG 4 was reported by several researchers under different conditions (Lichtenzveig et al., 2006; Tar’an et al., 2007; Anbessa et al., 2009; Sabbavarapu et al., 2013). Most recently, Li et al. (2017) identified a 100 kb genomic region containing 12 candidate genes for disease resistance associated with a major QTL on chromosome 4 of chickpea using Fst genome-scan and genome-wide association mapping.
Grain yield loss due to AB is a major cause for concern in pea growing regions. Several studies have been conducted to identify improved sources of resistance for pea breeding. Many QTLs were reported for AB resistance in pea (Timmerman-Vaughan et al., 2002, 2004; Tar’an et al., 2003; Prioul et al., 2004; Fondevilla et al., 2008, 2011; Carrillo et al., 2014). Under field conditions, Timmerman-Vaughan et al. (2002, 2004) reported several QTLs for resistance on LGs I, II, III, IV, V, VII, and Group A in two pea mapping populations, whereas Tar’an et al. (2003) identified three QTLs on LGs II, IV, and VI. Prioul et al. (2004) reported six QTLs on LGs III, Va, VI, and VII and 10 QTLs on LGs II, III, Va, and VII under controlled and field conditions, respectively. In P. sativum ssp. syriacum, six QTLs were reported on LGs II, III, IV and V by Fondevilla et al. (2008), whereas three additional QTLs were identified by Fondevilla et al. (2011) on LGs III and VI. Carrillo et al. (2014) identified four new QTLs on LGs II, III, and V controlling cellular mechanisms involved in AB resistance in P. sativum ssp. syriacum. A comparative analysis showed that QTL MpIII.1 (Fondevilla et al., 2008) was located on the same distal part of LG III where Prioul et al. (2004) identified mpIII-1. Fondevilla et al. (2011) indicated that QTLs MpIII.1, MpIII.3, and MpIII.2 detected in P. sativum ssp. syriacum corresponded to the QTLs mpIII-1, mpIII-3, and mpIII-5 identified in P. sativum by Prioul et al. (2004).
With the long-term objective to develop disease resistant pea cultivars, P651 (P. fulvum) a wild accession with improved resistance was identified and utilized for the development of an inter-specific pea population (PR-19) (Jha et al., 2012; Sindhu et al., 2014). Nine QTLs were identified for AB resistance in PR-19 and these QTLs individually explained 7.5 to 28% of phenotypic variation (Jha et al., 2016). Among these QTLs, abI-IV-2 and abIII-1 were consistent across locations and/or years with greater effects (16 to 28% of the total phenotypic variation) and P651 contributed alleles for disease resistance. Based on shared anchored markers, none of the identified QTLs were located in the regions of previously reported QTLs for AB resistance in pea (Jha et al., 2016).
In this research, abI-IV-2 and abIII-1 were selected for fine mapping to develop closely linked markers associated with AB resistance. For this purpose, four RILs each were identified in the abI-IV-2 and abIII-1 QTLs, for development of HIF populations. Among these RILs, lines PR-19-224 and PR-19-173 were selected for development of HIF-224 and HIF-173, respectively, on the basis of presence of heterozygous alleles as determined by KASP assay and AB scores. These HIFs served as segregating populations for fine mapping.
To find additional markers within QTLs, selective genotyping of 51 PR-19 RILs was performed using GBS. Based on linkage map construction from these RILs, 12 SNPs were identified in regions next to the highly linked markers within QTL abI-IV-2. Ten of the 12 markers were further genotyped on the complete set of PR-19 RILs (144) to determine the exact position and order of the tested markers in the QTL. Eight out of 10 SNPs from GBS were mapped within QTL abI-IV-2. Three markers (Sc34405_60551, Sc33468_44352, and Sc12023_67096) were located within the closest flanking markers (PsC6805p316 and PsC19558p107) located on either side of marker (PsC8031p219) having highest LOD in the QTL. The presence of QTL abI-IV-2 was validated on linkage map of PR-19 lines enriched with additional GBS markers. GBS marker Sc33287_25420 was the closest marker to the identified QTL with maximum LOD value and co-located with PsC6805p316. HIF-224 lines were genotyped using 20 SNPs including 10 GBS markers. A linkage map was constructed from 17 markers which covered 86.3 cM distance. The order and distance of markers were different compared to abI-IV-2 QTL obtained for PR-19. This could be due to recombination between the nearest markers within the QTL, or with markers near this QTL. In the case of PR-19, markers PsC20402p281 and PsC7497p542 were adjacent to the abI-IV-2 QTL and distant (around 7 cM) from the closest marker (PsC8031p219) to the QTL. However, in HIF-224, these markers were present within the QTL abI-IV-2 and covered more than 40 cM distance out of 86.3 cM. The larger map distance in HIF population compared to PR-19 RIL population could be due to the possibility that RIL PR-19-224 selected for HIF development was not heterozygous for the entire QTL. This RIL was fixed for alleles from Alfetta at several loci and was the best line that could be selected for maximum heterozygosity within this QTL based on genotyping of the F6 generation. Further, line PR-19-224 selected for HIF development might also contain positive alleles at other ascochyta resistance QTLs which might have affected determining the true effect of this QTL on disease resistance, thus there was no spike observed in LOD value in the HIF population.
Two new QTLs, abI-IV-2.1 and abI-IV-2.2 were identified within abI-IV-2 QTL due to additional SNP markers identified and these QTLs individually explained 5.5 to 14% of the total phenotypic variation. In general, improvement in LOD value was observed in comparison to previously identified QTL. QTLs for lodging resistance were co-located with QTLs associated with AB resistance. The parent Alfetta contributed the alleles for AB resistance as well as for lodging resistance. In this research, it was observed that the difference in AB score was relatively narrow between the parents under field conditions. On a 0-9 scale, Alfetta had 3 to 4 and 4 to 5 disease scores at pod filling and physiological maturity stage, respectively, whereas P651 had 2 to 3 at pod filling stage and 3 to 4 at physiological maturity stage. Further, Alfetta (1–3) had very low lodging score compared to P651 (8–9) on the 1–9 scale. Previous studies reported positive correlation between AB and lodging scores (Tar’an et al., 2003; Banniza et al., 2005; Jha et al., 2013, 2016). Under field conditions, lodging might play an important role in the disease progression and AB score could be related to disease avoidance rather than resistance per se. Our previous study (Jha et al., 2016) reported that three out of six QTLs identified under field conditions could account for disease avoidance as these loci were also associated with traits including lodging or plant height. Alternatively, resistance under field conditions could be due to physiological resistance (Khan et al., 2013) as canopy architecture features including branching, lodging resistance, and leaf area index could affect the impact of disease by affecting the microclimate within the canopy and splash dispersal of P. pinodes conidia (Schoeny et al., 2008; Le May et al., 2009). In case of abIII-1, no additional marker was identified within the QTL (2.1 cM). Five and three additional SNP markers identified in QTLs abI-IV-2.1 and abI-IV-2.2, respectively, by fine mapping can be used for marker assisted selection. Further, promising HIF lines harboring QTLs for disease resistance can be utilized as donors for development of cultivars with improved AB resistance.
AJ, BT, and TW conceived and designed the experiments. AJ and TW were involved in the development of HIFs and the multi-year field trials. GBS markers were developed by KG. AJ and KG were involved in the genotyping of PR-19 and HIFs and data analysis. AJ wrote the manuscript with input from KG, TW, and BT. All authors have read and approved the final manuscript.
We greatly acknowledge the financial support from the Saskatchewan Ministry of Agriculture, Saskatchewan Pulse Growers, and Western Grains Research Foundation.
Conflict of Interest Statement
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.
We are thankful to Kamal Bandara, Brent Barlow, Scott Ife, and other pulse crop field lab members for technical assistance.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/article/10.3389/fpls.2017.00765/full#supplementary-material
Anbessa, Y., Tar’an, B., Warkentin, T. D., Tullu, A., and Vandenberg, A. (2009). Genetic analyses and conservation of QTL for ascochyta blight resistance in chickpea (Cicer arietinum L.). Theor. Appl. Genet. 119, 757–765. doi: 10.1007/s00122-009-1086-2
Atienza, S. G., Palomino, C., Gutiérrez, N., Alfaro, C. M., Rubiales, D., Torres, A. M., et al. (2016). QTLs for ascochyta blight resistance in faba bean (Vicia faba L.): validation in field and controlled conditions. Crop Pasture Sci. 67, 216–224. doi: 10.1071/cp15227
Aveskamp, M. M., de Gruyter, J., Woudenberg, J. H. C., Verkley, G. J. M., and Crous, P. W. (2010). Highlights of the Didymellaceae: a polyphasic approach to characterise Phoma and related pleosporalean genera. Stud. Mycol. 65, 1–60. doi: 10.3114/sim.2010.65.01
Bai, X., Luo, L., Yan, W., Kovi, M. R., Zhan, W., and Xing, Y. (2010). Genetic dissection of rice grain shape using a recombinant inbred line population derived from two contrasting parents and fine mapping a pleiotropic quantitative trait locus qGL7. BMC Genet. 26:16. doi: 10.1186/1471-2156-11-16
Banniza, S., Hashemi, P., Warkentin, T. D., Vandenberg, A., and Davis, A. R. (2005). The relationship among lodging, stem anatomy, degree of lignification and susceptibility to mycosphaerella blight in field pea (Pisum sativum). Can. J. Bot. 83, 954–967. doi: 10.1139/b05-044
Bouteillé, M., Rolland, G., Balsera, C., Loudet, O., and Muller, B. (2012). Disentangling the intertwined genetic bases of root and shoot growth in Arabidopsis. PLoS ONE 7:e32319. doi: 10.1371/journal.pone.0032319
Carrillo, E., Satovic, Z., Aubert, G., Boucherot, K., Rubiales, D., and Fondevilla, S. (2014). Identification of quantitative trait loci and candidate genes for specific cellular resistance responses against Didymella pinodes in pea. Plant Cell Rep. 33, 1133–1145. doi: 10.1007/s00299-014-1603-x
Collard, B. C. Y., Jahufer, M. Z. Z., Brouwer, J. B., and Pang, E. C. K. (2005). An introduction to markers, quantitative trait loci (QTL) mapping and marker-assisted selection for crop improvement: the basic concepts. Euphytica 142, 169–196. doi: 10.1007/s10681-005-1681-5
Dwiyanti, M. S., Yamada, T., Sato, M., Abe, J., and Kitamura, K. (2011). Genetic variation of γ-tocopherol methyltransferase gene contributes to elevated a-tocopherol content in soybean seeds. BMC Plant Biol. 11:152. doi: 10.1186/1471-2229-11-152
Elshire, R. J., Glaubitz, J. C., Sun, Q., Poland, J. A., Kawamoto, K., Buckler, E. S., et al. (2011). A robust, simple genotyping-by-sequencing (GBS) approach for high diversity species. PLoS ONE 6:e19379. doi: 10.1371/journal.pone.0019379
Fondevilla, S., Almeida, N. F., Satovic, Z., Rubiales, D., Patto, M. C. V., Cubero, J. I., et al. (2011). Identification of common genomic regions controlling resistance to Mycosphaerella pinodes, earliness and architectural traits in different pea genetic backgrounds. Euphytica 182, 43–52. doi: 10.1007/s10681-011-0460-8
Fondevilla, S., Avila, C. M., Cubero, J. I., and Rubiales, D. (2005). Response to Ascochyta pinodes in a germplasm collection of Pisum spp. Plant Breed. 124, 313–315. doi: 10.1111/j.1439-0523.2005.01104.x
Fondevilla, S., Satovic, Z., Rubiales, D., Moreno, M. T., and Torres, A. M. (2008). Mapping of quantitative trait loci for resistance to Ascochyta pinodes in Pisum sativum subsp. syriacum. Mol. Breed. 21, 439–454. doi: 10.1007/s11032-007-9144-4
Jha, A. B., Arganosa, G., Tar’an, B., Diederichsen, A., and Warkentin, T. D. (2013). Characterization of 169 diverse pea germplasm accessions for agronomic performance, mycosphaerella blight resistance and nutritional profile. Genet. Resour. Crop Evol. 60, 747–761. doi: 10.1007/s10722-012-9871-1
Jha, A. B., Tar’an, B., Diapari, M., Sindhu, A., Shunmugam, A., Bett, K., et al. (2015). Allele diversity analysis to identify SNPs associated with ascochyta blight resistance in pea. Euphytica 202, 189–197. doi: 10.1007/s10681-014-1254-6
Jha, A. B., Tar’an, B., Stonehouse, R., and Warkentin, T. D. (2016). Identification of QTLs associated with improved resistance to ascochyta blight in an interspecific pea recombinant inbred line population. Crop Sci. 56, 2926–2939. doi: 10.2135/cropsci2016.01.0001
Jha, A. B., Warkentin, T. D., Gurusamy, V., Tar’an, B., and Banniza, S. (2012). Identification of ascochyta blight resistance in wild Pisum species for use in pea breeding. Crop Sci. 52, 2462–2468. doi: 10.2135/cropsci2012.04.0242
Khan, T. N., Timmerman-Vaughan, G. M., Rubiales, D., Warkentin, T. D., Siddique, K. H. M., Erskine, W., et al. (2013). Didymella pinodes and its management in field pea: challenges and opportunities. Field Crop Res. 148, 61–77. doi: 10.1016/j.fcr.2013.04.003
Lander, E., Green, P., Abrahamson, J., Barlow, A., Daley, M., Lincoln, S., et al. (1987). MAPMAKER: an interactive computer package for constructing primary genetic linkage maps of experimental and natural populations. Genomics 1, 174–181. doi: 10.1016/0888-7543(87)90010-3
Le May, C., Ney, B., Lemarchand, E., Schoeny, A., and Tivoli, B. (2009). Effect of pea plant architecture on spatiotemporal epidemic development of ascochyta blight (Mycosphaerella pinodes) in the field. Plant Pathol. 58, 332–343. doi: 10.1111/j.1365-3059.2008.01947.x
Li, Y., Ruperao, P., Batley, J., Edwards, D., Davidson, J., Hobson, K., et al. (2017). Genome analysis identified novel candidate genes for ascochyta blight resistance in chickpea using whole genome re-sequencing data. Front. Plant Sci. 8:359. doi: 10.3389/fpls.2017.00359
Lichtenzveig, J., Bonfil, D. J., Zhang, H.-B., Shtienberg, D., and Abbo, S. (2006). Mapping quantitative trait loci in chickpea associated with time to flowering and resistance to Didymella rabiei the causal agent of Ascochyta blight. Theor. Appl. Genet. 113, 1357–1369. doi: 10.1007/s00122-006-0390-3
Madoui, M., Labadie, K., Agata, L., Aury, J., Kreplak, J., Gali, K. K., et al. (2016). “Assembly of the pea genome by integration of high throughput sequencing (PacBio and Illumina) and whole genome profiling (WGPTM) data,” in Proceedings of the Plant and Animal Genome XXIV Conference 2016 (San Diego, CA: Plant and Animal Genome).
Meng, P. H., Macquet, A., Loudet, O., Marion-Poll, A., and North, H. M. (2008). Analysis of natural allelic variation controlling Arabidopsis thaliana seed germinability in response to cold and dark: identification of three major quantitative trait loci. Mol. Plant 1, 145–154. doi: 10.1093/mp/ssm014
Mohan, M., Nair, S., Bhagwat, A., Krishna, T. G., Yano, M., Bhatia, C. R., et al. (1997). Genome mapping, molecular markers and marker-assisted selection in crop plants. Mol. Breed. 3, 87–103. doi: 10.1023/A:1009651919792
Prioul, S., Frankewitz, A., Deniot, G., Morin, G., and Baranger, A. (2004). Mapping of quantitative trait loci for partial resistance to Ascochyta pinodes in pea (Pisum sativum L.) at the seedling and adult plant stages. Theor. Appl. Genet. 108, 1322–1334. doi: 10.1007/s00122-003-1543-2
Prioul-Gervais, S., Deniot, G., Receveur, E. M., Frankewitz, A., Fourmann, M., Rameau, C., et al. (2007). Candidate genes for quantitative resistance to Ascochyta pinodes in pea (Pisum sativum L.). Theor. Appl. Genet. 114, 971–984. doi: 10.1007/s00122-016-2669-3
Sabbavarapu, M. M., Sharma, M., Chamarthi, S. K., Swapna, N., Rathore, A., Thudi, M., et al. (2013). Molecular mapping of QTLs for resistance to Fusarium wilt (race 1) and Ascochyta blight in chickpea (Cicer arietinum L.). Euphytica 193, 121–133. doi: 10.1007/s10681-013-0959-2
Schoeny, A., Menat, J., Darsonval, A., Rouault, F., Jumel, S., and Tivoli, B. (2008). Effect of pea canopy architecture on splash dispersal of Mycosphaerella pinodes conidia. Plant Pathol. 57, 1073–1085. doi: 10.1111/j.1365-3059.2008.01888.x
Sindhu, A., Ramsay, L., Sanderson, L. A., Stonehouse, R., Li, R., Condie, J., et al. (2014). Gene-based SNP discovery and genetic mapping in pea. Theor. Appl. Genet. 127, 2225–2241. doi: 10.1007/s00122-014-2375-y
Sudheesh, S., Rodda, M. S., Davidson, J., Javid, M., Stephens, A., Slater, A. T., et al. (2016). SNP-based linkage mapping for validation of QTLs for resistance to ascochyta blight in lentil. Front. Plant Sci. 7:1604. doi: 10.3389/fpls.2016.01604
Tar’an, B., Warkentin, T., Somers, D. J., Miranda, D., Vandenberg, A., Balde, S., et al. (2003). Quantitative trait loci for lodging resistance, plant height and partial resistance to ascochyta blight in field pea (Pisum sativum L.). Theor. Appl. Genet. 107, 1482–1491. doi: 10.1007/s00122-003-1379-9
Tar’an, B., Warkentin, T. D., Tullu, A., and Vandenberg, A. (2007). Genetic mapping of ascochyta blight resistance in chickpea (Cicer arietinum L.) using a simple sequence repeat linkage map. Genome 50, 26–34. doi: 10.1139/g06-137
Timmerman-Vaughan, G. M., Frew, T. J., Butler, R., Murray, S., Gilpin, M., Falloon, K., et al. (2004). Validation of quantitative trait loci for Ascochyta blight resistance in pea (Pisum sativum L.), using populations from two crosses. Theor. Appl. Genet. 109, 1620–1631. doi: 10.1007/s00122-004-1779-5
Timmerman-Vaughan, G. M., Frew, T. J., Russell, A. C., Khan, T., Butler, R., Gilpin, M., et al. (2002). QTL mapping of partial resistance to field epidemics of ascochyta blight of pea. Crop Sci. 42, 2100–2111. doi: 10.2135/cropsci2002.2100
Timmerman-Vaughan, G. M., Moya, L., Frew, T. J., Murray, S. R., and Crowhurst, R. (2016). Ascochyta blight disease of pea (Pisum sativum L.): defence-related candidate genes associated with QTL regions and identification of epistatic QTL. Theor. Appl. Genet. 129, 879–896. doi: 10.1007/s00122-016-2669-3
Tivoli, B., Baranger, A., Avila, C. M., Banniza, S., Barbetti, M., Chen, W., et al. (2006). Screening techniques and sources of resistance to foliar diseases caused by major necrotrophic fungi in grain legumes. Euphytica 147, 223–253. doi: 10.1007/s10681-006-3131-4
Todesco, M., Balasubramanian, S., Hu, T. T., Traw, M. B., Horton, M., Epple, P., et al. (2010). Natural allelic variation underlying a major fitness trade-off in Arabidopsis thaliana. Nature 465, 632–636. doi: 10.1038/nature09083
Tuinstra, M. R., Ejeta, G., and Goldsbrough, P. B. (1997). Heterogeneous inbred family (HIF) analysis: a method for developing near-isogenic lines that differ at quantitative trait loci. Theor. Appl. Genet. 95, 1005–1011. doi: 10.1007/s001220050654
Udupa, S. M., and Baum, M. (2003). Genetic dissection of pathotype-specific resistance to ascochyta blight disease in chickpea (Cicer arietinum L.) using microsatellite markers. Theor. Appl. Genet. 106, 1196–1202. doi: 10.1007/s00122-002-1168-x
Watanabe, S., Xia, Z., Hideshima, R., Tsubokura, Y., Sato, S., Yamanaka, N., et al. (2011). A map-based cloning strategy employing a residual heterozygous line reveals that the GIGANTEA gene is involved in soybean maturity and flowering. Genetics 188, 395–407. doi: 10.1534/genetics.110.125062
Xue, A. G., Warkentin, T. D., Greeniaus, M. T., and Zimmer, R. C. (1996). Genotypic variability in seed borne infection of field pea by Ascochyta pinodes and its relation to foliar disease severity. Can. J. Plant Pathol. 18, 370–374. doi: 10.1080/07060669609500590
Xue, A. G., Warkentin, T. D., and Kenaschuk, E. O. (1997). Effect of timings of inoculation with Ascochyta pinodes on yield and seed infection on field pea. Can. J. Plant Sci. 77, 685–689. doi: 10.4141/P96-150
Zhang, R., Hwang, S. F., Chang, K. F., Gossen, B. D., Strelkov, S. E., Turnbull, G. D., et al. (2006). Genetic resistance to Ascochyta pinodes in 558 field pea accessions. Crop Sci. 46, 2409–2414. doi: 10.2135/cropsci2006.02.0089
Keywords: ascochyta blight, genotyping-by-sequencing, heterogeneous inbred family, quantitative trait loci, Pisum fulvum, P. sativum
Citation: Jha AB, Gali KK, Tar’an B and Warkentin TD (2017) Fine Mapping of QTLs for Ascochyta Blight Resistance in Pea Using Heterogeneous Inbred Families. Front. Plant Sci. 8:765. doi: 10.3389/fpls.2017.00765
Received: 08 March 2017; Accepted: 24 April 2017;
Published: 09 May 2017.
Edited by:Jennifer Davidson, South Australian Research and Development Institute, Australia
Copyright © 2017 Jha, Gali, Tar’an and Warkentin. 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) or licensor 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: Thomas D. Warkentin, firstname.lastname@example.org