Original Research ARTICLE
Identification of Quantitative Trait Loci Controlling Root and Shoot Traits Associated with Drought Tolerance in a Lentil (Lens culinaris Medik.) Recombinant Inbred Line Population
- 1Department of Plant Production, Faculty of Bioscience Engineering, Ghent University, Ghent, Belgium
- 2Institut National de la Recherche Agronomique du Maroc (INRA), Centre Régional de Settat, Settat, Morocco
- 3International Center for Agricultural Research in the Dry Areas, Institut National de la Recherche Agronomique Morocco Cooperative Research Project, Rabat, Morocco
- 4Plant Sciences Unit, Applied Genetics and Breeding, The Institute for Agricultural and Fisheries Research (ILVO), Melle, Belgium
- 5United States Department of Agriculture, Agricultural Research Service Grain Legume Genetics and Physiology Research, Pullman, WA, USA
- 6United States Department of Agriculture, Agricultural Research Service Western Regional Plant Introduction, Washington State University, Pullman, WA, USA
- 7Brotherton Seed Company, Washington, DC, USA
- 8Faculty of Tropical AgriSciences, Czech University of Life Sciences, Prague, Czech Republic
Drought is one of the major abiotic stresses limiting lentil productivity in rainfed production systems. Specific rooting patterns can be associated with drought avoidance mechanisms that can be used in lentil breeding programs. In all, 252 co-dominant and dominant markers were used for Quantitative Trait Loci (QTL) analysis on 132 lentil recombinant inbred lines based on greenhouse experiments for root and shoot traits during two seasons under progressive drought-stressed conditions. Eighteen QTLs controlling a total of 14 root and shoot traits were identified. A QTL-hotspot genomic region related to a number of root and shoot characteristics associated with drought tolerance such as dry root biomass, root surface area, lateral root number, dry shoot biomass and shoot length was identified. Interestingly, a QTL (QRSratioIX-2.30) related to root-shoot ratio, an important trait for drought avoidance, explaining the highest phenotypic variance of 27.6 and 28.9% for the two consecutive seasons, respectively, was detected. This QTL was closed to the co-dominant SNP marker TP6337 and also flanked by the two SNP TP518 and TP1280. An important QTL (QLRNIII-98.64) related to lateral root number was found close to TP3371 and flanked by TP5093 and TP6072 SNP markers. Also, a QTL (QSRLIV-61.63) associated with specific root length was identified close to TP1873 and flanked by F7XEM6b SRAP marker and TP1035 SNP marker. These two QTLs were detected in both seasons. Our results could be used for marker-assisted selection in lentil breeding programs targeting root and shoot characteristics conferring drought avoidance as an efficient alternative to slow and labor-intensive conventional breeding methods.
Lentil (Lens culinaris Medik.) is an important grain legume crop that is often grown in sustainable farming systems and for nutrition in the world. Its ability to enhance soil fertility through atmospheric nitrogen fixation allows substantial reduction in fertilizer use and significant production improvement in cereal-based cropping systems thanks to the benefits of rotation.
Lentil grains are a rich source of proteins and some important micronutrients such as iron and zinc (Grusak and Coyne, 2009; Thavarajah et al., 2011). Consumed as staple food in developing countries and as vegetarian dishes elsewhere, lentil grains are considered a very healthy food. The United Nations, in its 68th General Assembly, declared 2016 as the International Year of Pulses (annual leguminous crops harvested for dry grains) in order to highlight the nutritional benefits of pulses as part of sustainable food production aimed towards food and nutrition security (FAO, 2015).
In the arid and semi-arid areas and also in the context of climate change and global warming, drought is one of the major constraints that can limit lentil production and cause substantial yield losses (Malhotra et al., 2004; Stoddard et al., 2006; Sarker et al., 2009). Developing cultivars with enhanced drought tolerance by conventional breeding often has limited success due to the complexity of this trait and the difficulties with finding reliable and suitable phenotyping methods. For example, well-developed root systems have been shown to be linked to drought tolerance as an avoidance mechanism guaranteeing plant productivity under water-limited conditions (Kashiwagi et al., 2005; Sarker et al., 2005; Verslues et al., 2006; Gaur et al., 2008; Vadez et al., 2008; Aswaf and Blair, 2012; Comas et al., 2013; Idrissi et al., 2015a,b). However, it is difficult to screen large numbers of accessions for these root traits using conventional methods. Thus, applying a marker-assisted selection for these traits would offer an interesting alternative in breeding programs targeting drought tolerance. As such, identifying and mapping DNA markers linked to genes controlling rooting patterns associated with drought tolerance will assist in reliable and efficient identification and development of tolerant cultivars. Several studies have shown that root traits are polygenically controlled, whereas they also identified related quantitative trait loci (QTLs) for different species such as maize (Ruta, 2008), common bean (Cichy et al., 2009; Aswaf and Blair, 2012), barley (Sayed, 2011), soybean (Brensha et al., 2012) and chickpea (Kashiwagi et al., 2014).
Lentil has a genome size of about 4 Gbp (Arumuganathan and Earle, 1991); several kinds of DNA markers have been developed and mapped, including RAPDs, ISSRs, AFLPs, SRAPs, SSRs, and SNPs (Eujayl et al., 1998; Rubeena et al., 2003; Hamwieh et al., 2005; Saha et al., 2010; Sharpe et al., 2013). Idrissi et al. (2015a) confirmed evidence of high genetic variability, high heritability and polygenic control of root and shoot characteristics. However, to our knowledge, no QTLs related to root traits have been reported for lentil to date. Thus, the objective of this study was to identify and map QTLs related to root and shoot traits associated with drought tolerance in a lentil recombinant inbred line population (RIL) as a promising step towards a marker-assisted selection approach. It also aimed to investigate the stability of detected QTLs by performing the analysis on two consecutive seasons.
Materials and Methods
A recombinant inbred line (RIL) population developed from a cross between two contrasting parents, ILL6002 and ILL5888 (Saha et al., 2010), obtained from Fred J. Muehlbauer, USDA-ARS, Washington State University, Pullman, USA, was used in this study. The RIL population consisted of the two parents and 132 F6–8 lines. The lines were advanced to the F6–8 generation from individual F2using single seed descent. The ILL6002 parent is a vigorous line reported as drought tolerant and with a well-developed root system (Sarker et al., 2005; Singh et al., 2013). On the other hand, ILL5888 is a drought sensitive line and has a less-developed root system and vegetative biomass. The two parents also differ in disease resistance (Stemphylium blight), flowering and maturity time, seed diameter, 100-seed weight, growth habit and plant height (Saha et al., 2010, 2013).
RIL Root and Shoot Traits Phenotyping and Drought Tolerance Evaluation
This F6–8 population was previously characterized for root and shoot traits related to drought tolerance (Idrissi et al., 2015a). Briefly, the population was evaluated under greenhouse conditions for root and shoot traits associated with drought tolerance under two contrasting watering regimes (well-watered and progressive drought-stressed) using the standard nutrition solution EEG MESTSTOF 19-8-16 (4) for two consecutive growing seasons (2013 and 2014). A completely randomized block design with three replications was used. Four uniformly germinated seeds were planted in plastic pots (H 35 × D 24 cm) filled with fine perlite (diameter ≤ 2 mm) in order to be able to extract intact roots without damage (Day, 1991; Rabah Nasser, 2009). The initial moisture in all the pots of both watering regimes was 75% of field capacity. It decreased to about 22% for the drought-stressed regime where plants were watered only once in the beginning of the experiment, while it was maintained at 75% for the well-watered treatment by watering plants twice a week as described in Idrissi et al. (2015a). At 38 days after sowing, plants were carefully extracted without damage to the roots, then shoots and roots were separated into plastic bags. Washed roots were preserved in a refrigerator (4°C, 90% relative humidity) to avoid drying before being scanned using an EPSON Scan scanner. The images were then analyzed using Image J software (Abramoff et al., 2004) combined with Smart Roots software (Lobet et al., 2011). From the scanned images, taproot length (TRL; cm plant−1), average taproot diameter (TRD; mm plant−1), root surface area (RSA; cm2 plant−1) and lateral root number (LRN) were measured. Dry root and shoot biomass (DRW, DSW; mg plant−1) were measured after oven-drying at 72°C for 48 h. Chlorophyll content was estimated according to the SPAD values measured at 32 days after sowing using a SPAD-502Plus chlorophyll meter (Konica Minolta, Japan), four measures were taken in fully expanded leaves per plant. The wilting score (WS) corresponding to the degree of wilting severity was used to estimate drought tolerance using the following 0–4 score scale (Singh et al., 2013): 0 = healthy plants with no visible symptoms of drought stress; 1 = green plants with slight wilting; 2 = leaves turning yellowish green with moderate wilting; 3 = leaves yellow–brown with severe wilting; and 4 = completely dried leaves and/or stems. Seedling vigor (SV) was recorded following the 1–5 IBPGR and ICARDA (1985) scale: 1 = very poor; 2 = poor; 3 = average; 4 = good; 5 = excellent. Root–shoot ratio (RS ratio) was calculated by dividing the dry root weight by the dry shoot weight. Growth rate (GR; cm) was estimated as the gain of length between 12 (SL12DAS; cm) and 22 days after sowing (SL22DAS; cm; GR = SL22DAS–SL12DAS). Specific root length (SRL) and specific root surface area (SRSA) were estimated by dividing root length and root surface area, respectively, by dry root weight. All the measures were recorded as the mean value based on the four plants per individual genotype in each pot. A summary of genetic variation and heritability of these traits is provided in Table 4 (Supplementary Material).
The previously developed linkage map of Saha et al. (2010) was created based on the same RIL population used in this study (ILL6002 × ILL5888). The initial mapping data of Saha et al. (2010) which consisted of 23 SSR, 108 SRAP, and 30 RAPD markers, were kindly provided by the authors. The map was further enhanced using 220 polymorphic Single Nucleotide Polymorphism (SNP) markers developed using the Genotyping By Sequencing (GBS) technique and 180 polymorphic Amplified Fragment Length Polymorphism (AFLP) markers.
Genotyping by Sequencing for SNP Identification
SNP data were obtained from 92 (out of 132) RILs using GBS. The GBS procedure of Poland et al. (2012) was used, including their 48 bar-coded adapters with a Pst I overhang; genomic DNA was digested with the enzymes Pst I and Msp I. The ligation reaction was completed using bar-coded Adapter 1 and the common Y-adapter in a master mix of buffer, ATP and T4-ligase. Ligated samples were pooled and PCR-amplified in a single tube, producing libraries of 48 samples each. The libraries were sequenced on two lanes of Illumina HiSeq2000 (University of California Berkeley V.C. Genomic Sequencing Lab). The sequencing data were processed to remove low quality data using in-house scripts and analyzed using Stacks software (Catchen et al., 2011, 2013). Two hundred twenty SNPs that proved to be polymorphic between both parents of the RIL population ILL6002 × ILL5888 (Wong et al., 2015) were analyzed.
The AFLP protocol of Vos et al. (1995) with minor modifications (De Riek et al., 2001) was performed as described in Idrissi et al. (2015c). Out of 12 primer combinations tested, seven (EcoRI-ACA + MseI-CAG, EcoRI-ACA + MseI-CTG, EcoRI-ACA + MseI-CTT, EcoRI-ACG + MseI-CAA, EcoRI-AGC + MseI-CAA, EcoRI-AGC + MseI-CAG, EcoRI-AGC + MseI-CTG) were selected and used for genotyping the RIL population.
Linkage Analysis and Map Construction
Five hundred sixty-one molecular markers on 132 RILs were used for linkage analysis and construction of a linkage map using JoinMap®4 (Van Ooijen, 2006; Table 1). First, segregation according to Mendelian expectation ratio of 1:1 was tested using the chi-square test at a significance level of 0.05, markers with distorted segregation were removed prior to further analysis. The grouping tree of the JoinMap® program was calculated using independent LOD (Logarithm of odds) as grouping parameter with threshold ranges of 6 for start and 30 for end, and 1 for step. Stable sets of markers at higher LOD values were selected. After initial creation of groups, the Strongest Cross Link (SCL) information from the output results was used for inspecting assignment of markers to groups, those with SCL-values larger than 5, indicating that they have strong linkage outside their respective groups, were assigned to the corresponding groups. This was repeated until all markers of each group had SCL-values smaller than 5. Linkage groups were calculated using the maximum likelihood mapping algorithm with default values as in the software. Map order in each linkage group was verified using the regression mapping algorithm with the following parameters: LOD threshold larger than 4, recombination frequency smaller than 0.25, Kosambi function as mapping function for genetic distance calculation and the second round map of the algorithm. The final linkage map was generated using MapChart© 2.3 program (Voorrips, 2002).
QTL analysis was performed for each season separately for drought-stressed treatments in order to check the stability of detected QTLs using MapQTL® 5 program (Van Ooijen, 2004). First, Kruskal-Wallis test was performed to determine a set of markers linked to each quantitative trait. Simple Interval Mapping was performed to identify linkage groups and positions with significant LOD scores. For each trait, LOD score threshold was determined based on a permutation test using 1000 iterations at a P-value of 0.05; LOD scores above these values were considered as significant. Co-factor selection was performed based on automatic co-factor selection implemented in the software for each linkage group and on manual selection of individual markers with significant LOD scores from Simple Interval Mapping output before applying Multiple-QTL Models (MQM) mapping (also called Composite Interval Mapping). Performing MQM mapping with markers close to significant LOD score positions as co-factors allows reduction of residual variance, thus enhancing the power of QTL detection. For each quantitative trait, co-factor selection and MQM mapping were repeated until no further enhancement was obtained (no more QTLs detected, increase in LOD scores and explained variances). From the MQM mapping output, closest marker, flanking markers, additive affect and percentage of explained variance for each detected QTL and for each quantitative trait were determined for seasons. Final results, with significant LOD scores and intervals, for each detected QTL per linkage group, were generated using MapChart© 2.3 program (Voorrips, 2002).
All detected QTLs were named as follows: Q‘Trait name abreviation'“linkage group number” − “position in cM”. For example: QLRNIII−98.64 is a QTL associated with LRN identified in linkage group III at position 98.64 cM.
GBS for SNP Identification
Selection of the genotyping-by-sequencing two enzyme method of Poland et al. (2012) and the enzymes Msp I and Pst I was based on the results of Wong et al. (2015) lentil SNP discovery across the lentil species. Using GBS, 220 polymorphic SNPs were deemed high quality for mapping, after satisfying quality control filtering based on deleting low quality and redundant SNPs using haplotype information for read depth (3), lack of redundancy and segregation in the parents. Genome coverage was reasonable, but incomplete, across six linkage groups (LG I, LG II, LG III, LG IV, LG VI, and LG IX; Figure 1).
Figure 1. Genetic linkage map of lentil developed at LOD score > 6 using maximum likelihood mapping algorithm of JoinMap® 4 program and drawn using MapChart© program. LGI-LGIX correspond to linkage groups, marker names are presented right to the linkage group and the genetic positions in CentiMorgans (cM) left. SNP markers are denoted by *TP_, SSR by *GLLC_, SRAP starting by *M_ or F_, RAPD by *UBC_ and AFLP by PC_.
Linkage Analysis and Map Construction
Marker distortion tested by Chi-square test (P < 0.05) revealed that 35.4% of SNPs, 43% of SSRs, 18% of SRAPs, 52.7% of AFLPs and 20% of RAPDs did not segregate according to the expected 1:1 ratio and were removed from the analysis. Out of 17 stable groups selected from the grouping tree, a total of 252 out of the 561 polymorphic markers were finally mapped in nine linkage groups spanning a total length of 2022.8 cM (Tables 1, 2). Final linkage groups were established using the SCL information. Linkage group length ranged from 71.7 to 531.1 cM whereas average distance between two markers ranged from 5.12 (LG V) to 9.8 cM (LG II) (Table 2; Figure 1). Seven linkage groups had a length of more than 100 cM (LG I, LG II, LG III, LG IV, LG VI, LG VII, and LG IX). Both co-dominant (SNP and SSR) and dominant (SRAP, AFLP, and RAPD) markers were present in six linkage groups, while three linkage groups (LG V, LG VII, and LG VIII) were composed only out of dominant markers.
A total number of 18 QTLs associated with 14 root and shoot traits were detected under drought-stressed conditions during two seasons (Table 3; Figures 2, 3). LOD score, percentage of explained phenotypic variance and additive effect of detected QTLs ranged from 2.75 (TRL) to 8.14 (DSW), from 4.3 (QRSratiIX-77.72) to 28.9% (QRSratioIX-2.30) and from −5.17 (LRN) to 8.10 (DRW), respectively.
Table 3. Characteristics of quantitative trait loci (QTL) identified under progressive drought stress in the RIL population (ILL6002 × ILL5888) for the 2013 and 2014 seasons.
Figure 2. Linkage groups with identified QTLs related to root and shoot traits under progressive drought stress for the 2013 season detected by MapQTL® 5 program and drawn using MapChart© program. LOD score curves are presented right to linkage groups and significant thresholds are presented by dotted lines (when more than one QTL for different quantitative traits are detected in same position, the dotted line correspond to the smallest threshold value). WS_DR_2013, wilting score; SPAD_DR_2013, chlorophyll content as estimated by SPAD value; LRN_DR_2013, lateral root number; SL12DAS_DR_2013, shoot length at 12 days after sowing; SL22DAS_DR_2013, shoot length at 22 days after sowing; SRL_DR_2013, specific root length; DRW_DR_2013, dry root weight; DSW_DR_2013, dry shoot weight; SV_DR_2013, seedling vigor; RSA_DR_2013, root surface area; GR_DR_2013, growth rate; RSratio_DR_2013, root-shoot ratio.
Figure 3. Linkage groups with identified QTLs related to root and shoot traits under progressive drought stress for the 2014 season detected by MapQTL® 5 program and drawn using MapChart© program. LOD score curves are presented right to linkage groups and significant thresholds are presented by dotted lines (when more than one QTL for different quantitative traits are detected in same position, the dotted line correspond to the smallest threshold value). RSratio_DR_2014, root-shoot ratio; LRN_DR_2014, lateral root number; SRL_DR_2014, specific root length; TRL_DR_2014, taproot length; SL12DAS_DR_2014, shoot length at 12 days after sowing; DRW_DR_2014, dry root weight; DSW_DR_2014, dry shoot weight; SL22DAS_DR_2014, shoot length at 22 days after sowing; RSA_DR_2014, root surface area; SPAD_DR_2014, Chlorophyll content as estimated by SPAD value; RSratio_DR_2014, root-shoot ratio; TRD_DR_2014, average taproot diameter.
Seven of the detected QTLs were co-located on LG VII at position 21–22 cM, with UBC34 as the closest marker and ME5XR10—UBC1 as the two flanking markers: QDRWVII-21.94, QLRNVII-21.94, QRSAVII-21.94, QDSWVII-21.94, QSL12VII-20.75, QSL22VII-21.75, and QGRVII-21.94.
Among the 18 detected QTLs, 12 were evidenced for the drought-stressed treatment for both seasons: QDRWVII-21.93, QRSAVII-21.94, QDSWVII-22.94, QRSratioIX-2.30, QSL12IV-103.83, QSL12VI-170.87, QSL12VII-19.71, QSL22VII-21.94, QLRNIII-98.64, QLRNVII-21.94, QSRLIV-61.63, and QSPADVIII-72.15. Interestingly, among these stable QTLs,QRSratioIX-2.30, located at 2.30 cM on LG IX, is associated with a high root-shoot ratio and had LOD scores of 6.20 and 5.11 for 2013 and 2014 seasons, respectively. The explained phenotypic variance of this QTL was the highest with 27.6 and 28.9% and an additive effect of 1.23 and 1.84 for 2103 and 2014 seasons, respectively. The closest marker to this QTL is SNP marker TP6337 located at 2.3 cM whereby the two flanking markers are TP518 and TP1280, located respectively at 0 and 2.9 cM.
Two QTLs were identified for dry root biomass, QDRWVII-21.93, accounted for 22.2% (with a LOD score of 7.21) and 21.3% (with a LOD score of 6.88) of the phenotypic variance with additive effects of 8.10 and 7.47 for 2013 and 2014 seasons, respectively.
Among the three QTLs detected for LRN, QLRNIII-98.64, was located at 98.64 cM position on LG III close to TP3371 SNP marker and flanked by the two SNP markers TP5093–TP6072. The LOD scores, percentage of explained phenotypic variances and additive effects were 2.94, 23.5% and −5.17, and 3.31, 24%, and −5.15 for 2013 and 2014 seasons, respectively. An important QTL was also identified for SRL, namely QSRLIV-61.63, that was detected for both seasons with LOD scores, percentage of explained phenotypic variances and additive effects of 3.84, 16.8% and 0.83 and 3.63, 16.2% and 0.32, respectively, for 2013 and 2014.
Three QTLs were identified to be linked to chlorophyll content in which one was common for both seasons. The latter is the QTL QSPADVIII-72.15, which was detected with LOD scores, percentage of explained phenotypic variances and additive effects of respectively 3.98, 10.7% and −2.20 for 2013 season and 4.25, 13.1% and −2.20 for 2014 season.
Also, a QTL related to early vegetative vigor estimated by SV was detected for the 2013 experiment. This QTL, QSVVII4, was located on LG VII at position 4 cM, had a LOD score of 3.46, an additive effect of 0.29 and explained 14.9% of total phenotypic variance.
A QTL QWSI-22.53, related to drought tolerance as estimated by WS, is located at 22.53 cM position on LG I with a LOD score of 3.08 and 18.8% as percentage of explained phenotypic variance.
The genetic linkage map of lentil initially developed by Saha et al. (2010) using a ILL6002 × ILL5888 RIL population containing 139 markers and 14 linkage groups was enhanced by adding SNP co-dominant markers and AFLP dominant markers, thereby increasing marker density and total spanned length. The number of linkage groups was reduced to nine with a total number of 252 mapped markers covering 2022.8 cM compared to 1565.2 cM in the previous genetic map. Average distance between markers was reduced from 11.3 to 8 cM. Sharpe et al. (2013) reported a lentil map with seven linkage groups using SNP and SSR markers. 62.77% of markers from the Saha et al. (2010) linkage map were also mapped in the genetic map developed in our study. Several sets of markers from the previous genetic map were confirmed to be linked to each other in our map. For instance, all markers from LG 1 from the map of Saha et al. (2010) were also mapped in LG I of our map. Thirteen markers out of a total of 19 mapped in LG 2 were mapped in LG III and four in LG IV of our map. Nine markers from LG 3 were mapped in LG V of our map and all those from LG 4 except for two that ended up in LG VII of our enhanced map. All markers from LG 11 of the previous map (except two) were mapped in LG II. All markers of LG 13 and LG 14 were mapped in LG IX and LG VIII of our map, respectively. Our linkage groups could not be assigned per Sharpe et al. (2013), the best lentil linkage map with seven linkage groups likely corresponding to the seven chromosomes of the genome developed to date, due to lack of common markers. We used a combination of dominant and co-dominant markers to develop a linkage map with reduced gaps. Since SNP data were not available for the whole population, we added also dominant AFLP markers for map construction. In other studies, dominant markers were also used together with co-dominant ones for the development of linkage maps and QTL analysis to overcome different limits such as genetic marker availability and large gaps in linkage groups (Gaudet et al., 2007; De Keyser et al., 2010; Kaur et al., 2014; Muys et al., 2014; Ting et al., 2014). Although, maximum likelihood mapping algorithm often results in increased map length, it is considered to be more robust with missing data, genotyping errors and the use of markers with low information content (Lincoln and Lander, 1992; Van Ooijen, 2006; Cartwright et al., 2007; De Keyser et al., 2010). This algorithm uses multipoint analysis to approximate missing genotypes using nearby markers (De Keyser et al., 2010). Genetic linkage maps based on this approach giving the most likely marker order (De Keyser et al., 2010) were reported to be suitable for QTL mapping (Kim, 2007; De Keyser et al., 2010). Thus, we adopted this approach as the main objective of our study was to identify QTLs related to root and shoot traits. Furthermore, although we used dominant markers such as AFLPs known to result in longer map, our linkage groups did not have extreme lengths and the total map length of 2022.8 cM is among common reported values in similar studies on lentil. Duran and Perez De La Vega (2004) reported a genetic linkage map of 2172 cM length using SSR, AFLP, ISSR and RAPD markers. Gupta et al. (2012) used SSR, ISSR and RAPD markers to construct a map of 3843.4 cM length. Also, Kaur et al. (2014) used SSR and SNP markers to develop a map of 1178 cM length. Using predominantly SNP markers and few SSRs, Sharpe et al. (2013) constructed a shorter map of 834.7 cM length. More recently, Ates et al. (2016) developed a map spanning a total length of 4060.6 cM and composed of seven linkage groups using SSR and SNP markers to identify QTLs controlling genes for Selenium uptake in lentil.
High genetic variability, quantitative, continuous and normally distributed variation as well as high heritability estimate values of all studied traits were reported in Idrissi et al. (2015a).
In all, 18 QTLs were identified for root and shoot traits for both seasons under progressive drought-stressed treatments in the lentil RIL population ILL6002 × ILL5888. Among these QTLs, 12 were evidenced for both seasons. Aswaf and Blair (2012) reported a total of 15 putative QTLs for seven rooting pattern traits and four shoot traits under drought-stressed treatments in common bean (Phaseolus vulgaris L.). Varshney et al. (2014) reported drought tolerance-related root trait QTLs in chickpea (Cicer arietinum L.). In soybean (Glycine max L.), Manavalan et al. (2015) identified a QTL region controlling a number of root and shoot architectural traits. In lentil, to our knowledge, this is the first report on QTLs related to root and shoot traits associated with drought tolerance. Interestingly, QTL QRSratioIX−2.30 related to root-shoot ratio, an important trait for drought avoidance (Verslues et al., 2006), was confirmed to be present on LG IX at 2.30 cM position during the two seasons. Among detected QTLs, this QTL explained the highest percentage of phenotypic variance and was close to the co-dominant SNP marker TP6337 (C/T) and furthermore was flanked by the two SNP markers TP518 (A/G) and TP1280 (G/T). These markers are potentially important for their practical use for marker-assisted selection in breeding programs targeting drought tolerance. It should be pointed out that the same SNP markers were confirmed as being linked to root-shoot ratio when using only SNP markers on 92 RILs for linkage map construction and QTL mapping (data not shown).
A QTL-“hotspot” genomic region was identified on LG VII close to UBC34 RAPD marker and ME4XR16c SRAP marker, and was identified to be linked to the genetic control of a number of root and shoot traits for both seasons: DRW, LRN, RSA, DSW, and SL at 12 and 22 days after sowing. These traits were shown to be significantly correlated (Idrissi et al., 2015a). Similarly, a QTL-“hotspot” related to 12 root traits was reported by Varshney et al. (2014) in chickpea (Cicer arietinum L.). Although practical efficient use of the identified genomic region in the ILL6002 × ILL5888 lentil population for marker-assisted selection could be limited by the dominant character of the closest RAPD marker, SRAP markers identified close to this genomic region could be used for assisting in the selection for linked traits. SRAP markers targeting the coding regions of open-reading frames of the genome, considered as better than RAPDs and technically less challenging than AFLPs, are of potential interest for QTL mapping (Chen et al., 2007; Yuan et al., 2008; Zhang et al., 2009; Saha et al., 2010, 2013; Robarts and Wolfe, 2014). Furthermore, up to 20% of SRAP markers were found to be co-dominant (Li and Quiros, 2001). Dry root weight, reported to be associated with drought tolerance by Idrissi et al. (2015a) in lentil, and other root and shoot traits such as root surface area and dry shoot weight also associated with drought tolerance are linked to this “hotspot” genomic region.
QTL QLRNIII−98.64, related to LRN located at 98.64 cM position on LG III, was identified during both seasons explained 23.5 and 24% variations for 2013 and 2014, respectively. This QTL was close to SNP marker TP3371 (C/T) whereas its significant interval is between TP5093 (C/T) and TP6072 (A/G) SNP markers. Thus, the efficient use of these markers in breeding programs is possible for screening for higher LRN. High lateral root number was previously reported to be associated with drought tolerance and yield in lentil under drought stress (Sarker et al., 2005). Similarly, QTL QSRLIV−61.63, located at 61.63 cM on LG IV and related to SRL, was detected in both seasons with fairly high LOD scores of 3.84 and 3.63 for 2013 and 2014 respectively. This QTL, explaining 16.8% of phenotypic variance, was close to TP1873 (A/C) SNP marker (61.6 cM) and flanked by the two F7XEM6b SRAP (52.1 cM) and TP1035 (A/T) SNP markers (65 cM). Specific root length is considered an important root trait that can contribute to plant productivity under drought (Comas et al., 2013). Therefore, the use of these linked markers to screen lines with longer root length should be of potential interest. Three QTLs were identified for chlorophyll content as estimated by the SPAD value. Among them, QSPADI−158.76 is located at 158.76 cM on LG I close to AFLP marker PC3_208 and flanked by the two co-dominant SNP markers TP1954 (A/T) and TP5642 (A/T) that could be efficiently used in marker-assisted selection. Idrissi et al. (2015a) reported correlations of SPAD value of 0.46 and 0.48 with dry root biomass and drought tolerance, respectively, in the same mapping population used here.
A QTL QWSI−22.53 related to drought tolerance as estimated by the WS, located at 22.53 cM on LG I and explaining 18.8% of total phenotypic variance, is close to TP5779 (A/T) and flanked by TP6354 (C/T) and TP1655 (C/T) SNP markers, was identified during the 2013 season. After validation, these markers could be used for screening for drought tolerance. Wilting severity due to drought stress was reported to be correlated with relative water content in lentil (Idrissi et al., 2015b) indicating the importance of this parameter for the identification of drought tolerant cultivars. QTLs for drought tolerance as estimated by relative water content were reported for pea (Pisum sativum) by Iglesias-García et al. (2015).
A drought tolerance breeding strategy could be first based on laboratory screening of large collections of genetic material for the presence of the identified markers. Then, lines carrying alleles linked to QTLs of targeted traits could be evaluated under field conditions to finally identify drought-tolerant individuals. More focus should be on QTLs related to root-shoot ratio, LRN, SRL, and WS shown to be flanked by SNPs markers. However, the QTL related to WS needs further evaluation under different watering conditions and drought intensity to determine environments of expression of this QTL. This will allow to determine whether it is an adaptive or constitutive QTL.
It should be pointed out that results of QTL analysis using the second round map of JoinMap®4 program (Van Ooijen, 2006) obtained from regression mapping algorithm were closely similar to those obtained using maximum likelihood algorithm, although total lengths of the two maps were different (data not shown).
In this study, a total of 18 QTLs related to root and shoot traits associated with drought tolerance such as dry root biomass, LRN, root-shoot ratio, and specific root length were identified under progressive drought-stressed treatment. Interestingly, 12 of these QTLs were detected for both seasons, confirming their potential importance in conveying drought tolerance. DNA markers linked to these QTLs could be used for marker-assisted selection, thus making subsequent breeding efforts more reliable and efficient as the respective phenotyping-based methods are slow and labor-intensive, and affected by environment. Although, root characteristics are difficult to study as many environmental effects (especially soil characteristics) interact with genetic factors, our results provide significant information about QTLs related to root and shoot traits that could be used in marker-assisted breeding after validation.
OI designed the study, analyzed data, interpreted results and wrote the paper. OI, JDR, PVD, EDK, SU contributed to design the study, analyze data, interpret results and wrote the paper. OI, CC, RM, GS contributed to data acquisition. SU, FM contributed to critically revising the paper.
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.
OI and SU thank OCP-Foundation-INRA-ICARDA-IAV Hassan II project on India-Morocco Food Legume Initiative for the support.
Aswaf, A., and Blair, M. (2012). Quantitative trait loci for rooting pattern traits of common beans grown under drought stress versus non-stress conditions. Mol. Breed. 30, 681–695. doi: 10.1007/s11032-011-9654-y
Ates, D., Sever, T., Aldemir, S., Yagmur, B., Temel, H. Y., Kaya, H. B., et al. (2016). Identification of QTLs controlling genes for Se uptake in lentil seeds. PLoS ONE 11:e0149210. doi: 10.1371/journal.pone.0149210
Brensha, W., Kantartzi, S. K., Meksem, K., Grier, R. L., Barakat, A., Lightfoot, D. A., et al. (2012). Genetic analysis of root and shoot traits in the “Essex” by “Forrest” recombinant inbred line (RIL) population of soybean (Glycine max (L.) Merr.). J. Plant Genome Sci. 1, 1–9. doi: 10.5147/jpgs.2012.0051
Catchen, J. M., Amores, A., Hohenlohe, P., Cresko, W., and Postlethwait, J. H. (2011). Stacks: building and genotyping loci de novo from short-read sequences. G3: Genes Genom. Genet. 1, 171–182. doi: 10.1534/g3.111.000240
Chen, W., Zhang, Y., Liu, X., Chen, B., Tu, J., and Tingdong, F. (2007). Detection of QTL for six yield-related traits in oilseed rape (Brassica napus) using DH and immortalized F(2) populations. Theor. Appl. Genet. 115, 849–858. doi: 10.1007/s00122-007-0613-2
Cichy, K., Blair, M. W., Mendoza, C. H. G., Snapp, S. S., and Kelly, J. D. (2009). QTL analysis of root architecture traits and low phosphorus tolerance in an Andean bean population. Crop Sci. 49, 59–68. doi: 10.2135/cropsci2008.03.0142
Comas, L. H., Becker, S. R., Cruz, V. M. V., Byrne, P. F., and Dierig, D. A. (2013). Root traits contributing to plant productivity under drought. Front. Plant Sci. 4:442. doi: 10.3389/fpls.2013.00442
De Keyser, E., Yan Shu, Q., Van Bockstaele, E., and De Riek, J. (2010). Multipoint-likelihood maximization mapping on 4 segregating populations to achieve an integrated framework map for QTL analysis in pot azalea (Rhododendron simsii hybrids). BMC Mol. Biol. 11:1. doi: 10.1186/1471-2199-11-1
De Riek, J., Calsyn, E., Everaert, I., Van Bockstaele, E., and De Loose, M. (2001). AFLP based alternatives for the assessment of distinctness, uniformity and stability of sugar beet varieties. Theor. Appl. Genet. 103, 1254–1265. doi: 10.1007/s001220100710
Duran, Y., and Perez De La Vega, M. (2004). Assessment of genetic variation and species relationships in a collection of Lens using RAPD and ISSR. Span J. Agric. Res. 2, 538–544. doi: 10.5424/sjar/2004024-110
Eujayl, I., Baum, M., Powell, W., Erskine, W., and Pehu, E. (1998). A genetic linkage map of lentil (Lens sp.) based on RAPD and AFLP markers using recombinant inbred lines. Theor. Appl. Genet. 97, 83–89. doi: 10.1007/s001220050869
FAO (2015). International Year of Pulses 2016. Available online at: http://www.fao.org/pulses-2016/en/ (Accessed on: 24 Nov 2015).
Gaudet, M., Jorge, V., Paolucci, I., Beritognolo, I., Scarascia Mugnozza, G., and Sabatti, M. (2007). Genetic linkage maps of Populus nigra L. including AFLPs, SSRs, SNPs, and sex trait. Tree Genet. Genom. 4, 25–36. doi: 10.1007/s11295-007-0085-1
Gaur, P. M., Krishnamurthy, L., and Kashiwagi, J. (2008). Improving drought avoidance root traits in chickpea (Cicer arietinum)—current status of research at ICRISAT. Plant Prod. Sci. 11, 3–11. doi: 10.1626/pps.11.3
Grusak, M. A., and Coyne, C. J. (2009). “Variation for seed minerals and protein concentrations in diverse germplasm of lentilm” in Paper Presented at North America Pulse Improvement Association, 20th Biennial Meeting (Fort Collins, CO), 11.
Gupta, M., Verma, B., Kumar, N., Chahota, R. K., Rathour, R., Sharma, S. K., et al. (2012). Construction of intersubspecific molecular genetic map of lentil based on ISSR, RAPD and SSR markers. J. Genet. 91, 279–287. doi: 10.1007/s12041-012-0180-4
Hamwieh, A., Udupa, S. M., Choumane, W., Sarker, A., Dreyer, F., Jung, C., et al. (2005). A genetic linkage map of Lens sp. based on microsatellite and AFLP markers and the localization of fusarium vascular wilt resistance. Theor. Appl. Genet. 110, 669–677. doi: 10.1007/s00122-004-1892-5
Idrissi, O., Houasli, C. H., Udupa, S. M., De Keyser, E., Van Damme, P., and De Riek, J. (2015a). Genetic variability for root and shoot traits in a lentil (Lens culinaris Medik.) recombinant inbred line population and their association with drought tolerance. Euphytica 204, 693–709. doi: 10.1007/s10681-015-1373-8
Idrissi, O., Udupa, S. M., De Keyser, E., Van Damme, P., and De Riek, J. (2015b). Functional genetic diversity analysis and identification of associated simple sequence repeats and amplified fragment length polymorphism markers to drought tolerance in lentil (Lens culinaris ssp. culinaris Medicus) Landraces. Plant Mol. Biol. Rep. 34, 659–680. doi: 10.1007/s11105-015-0940-4
Idrissi, O., Udupa, S. M., Houasli, C. H., De Keyser, E., Van Damme, P., and De Riek, J. (2015c). Genetic diversity analysis of Moroccan lentil (Lens culinaris Medik.) landraces using simple sequence repeat and amplified fragment length polymorphisms reveals functional adaptation toward agro-environmental origins. Plant Breed. 134, 322–332. doi: 10.1111/pbr.12261
Iglesias-García, R., Prats, E., Fondevilla, S., Satovic, Z., and Rubiales, D. (2015). Quantitative trait loci associated to drought adaptation in pea (Pisum sativum L.). Plant Mol. Biol. Rep. 33, 1768–1778. doi: 10.1007/s11105-015-0872-z
Kashiwagi, J., Krishnamurthy, L., Gaur, P. M., Shandra, S., and Upadhyaya, H. D. (2014). Estimation of gene effects of the drought avoidance root characteristics in chickpea (C. arietinum L.). Field Crop Res. 105, 64–69. doi: 10.1007/s00122-013-2230-6
Kashiwagi, J., Krishnamurthy, L., Upadhyaya, H. D., Krishna, H., Chandra, S., Vadez, V., et al. (2005). Genetic variability of drought-avoidance root traits in the mini-core germplasm collection of chickpea (Cicer arietinum L.). Euphytica 146, 213–222. doi: 10.1007/s10681-005-9007-1
Kaur, S., Cogan, N. O. I., Stephens, A., Noy, D., Butsch, M., Forster, J. W., et al. (2014). EST-SNP discovery and dense genetic mapping in lentil (Lens culinaris Medik.) enable candidate gene selection for boron tolerance. Theor. Appl. Genet. 127, 703–713. doi: 10.1007/s00122-013-2252-0
Li, G., and Quiros, C. F. (2001). Sequence-related amplified polymorphism (SRAP), a new marker system based on a simple PCR reaction: its application to mapping and gene tagging in Brassica. Theor. Appl. Genet. 103, 455–461. doi: 10.1007/s001220100570
Malhotra, R. S., Sarker, A., and Saxena, M. C. (2004). “Drought tolerance in chickpea and lentil-present status and future strategies,” in Challenges and Strategies for Dryland Agriculture, eds S. C. Rao and J. Ryan (Wisconsin, EI: CSSA Special Publication no. 32, Crop Science Society of America and American Society of Agronomy), 257–273.
Manavalan, L. P., Prince, S. J., Musket, T. A., Chaky, J., Deshmukh, R., Vuong, T. D., et al. (2015). Identification of novel QTL governing root architectural traits in an interspecific soybean population. PLoS ONE 10:e0120490. doi: 10.1371/journal.pone.0120490
Muys, C., Noelle, C., Hienpont, T., Dauchot, N., Maudoux, O., Draye, X., et al. (2014). Integration of AFLPs, SSRs and SNPs markers into a new genetic map of industrial chicory (Cichorium Intybus L. var. sativum). Plant Breed. 133, 130–137. doi: 10.1111/pbr.12113
Poland, J. A., Brown, P. J., Sorrells, M. E., and Jannink, J. (2012). Development of high density maps for barley and wheat using a novel two-enzyme genotyping by sequencing approach. PLoS ONE 7:e32253. doi: 10.1371/journal.pone.0032253
Rabah Nasser, R. (2009). Physiological Aspects of the Response to Elevated CO2 in Lentils (Lens culinaris Medic.). Dissertation thesis: School of Biological Sciences of the University of Plymouth, Devon.
Robarts, D. W. H., and Wolfe, A. D. (2014). Sequence-Related Amplified Polymorphism (SRAP) Markers: a potential resource for studies in plant molecular biology. Appl. Plant Sci. 2:1400017. doi: 10.3732/apps.1400017
Saha, G. C., Sarker, A., Chen, W., Vandemark, G. J., and Muehlbauer, F. J., (2010). Inheritance and linkage map positions of genes conferring resistance to Stemphylium Blight in lentil. Crop Sci. 50, 1831–1839. doi: 10.2135/cropsci2009.12.0709
Saha, G. C., Sarker, A., Chen, W., Vandemark, G. J., and Muehlbauer, F. J. (2013). Inheritance and linkage map positions of genes conferring agromorphological traits in Lens culinaris Medik. Int. J. Agron. 2013. doi: 10.1155/2013/618926
Sarker, A., Aydogan, A., Chandra, S., Kharrat, M., and Sabaghpour, S. (2009). “Genetic enhancement for yield and yield stability,” in The Lentil: Botany, Production and Uses, eds W. Erskine, F. J. Muehlbauer, A. Sarker, and B. Sharma (Oxfordshire: CAB International), 102–120.
Sarker, A., Erskine, W., and Singh, M. (2005). Variation in shoot and root characteristics and their association with drought tolerance in lentil landraces. Genet. Resour. Crop. Evol. 52, 89–97. doi: 10.1007/s10722-005-0289-x
Sayed, M. A. A. (2011). QTL Analysis for Drought Tolerance Related to Root and Shoot Traits in Barley (Hordeum vulgare L.). Dissertation thesis, Institut für Nutzpflanzenwissenschaften und Ressourcenschutz, Hohen Landwirtschaftlichen Fakultät der Rheinischen Friedrich-Wilhelms-Universität Bonn, Bonn.
Sharpe, A. G., Ramsay, L., Sanderson, L. A., Fedoruk, M. J., Clarke, W. E., Li, R., et al. (2013). Ancient orphan crop joins modern era: gene-based SNP discovery and mapping in lentil. BMC Genomics 14:192. doi: 10.1186/1471-2164-14-192
Singh, D. H., Dikshit, H. K., and Singh, R. A. (2013). A new phenotyping technique for screening for drought tolerance in lentil (Lens culinaris Medik.). Plant Breed. 132, 185–190. doi: 10.1111/pbr.12033
Stoddard, F. L., Balko, C., Erskine, W., Khan, H. R., Link, W., and Sarker, A. (2006). Screening techniques and sources of resistance to abiotic stresses in cool-season food legumes. Euphytica 147, 167–186. doi: 10.1007/s10681-006-4723-8
Thavarajah, D., Thavarajah, P., Sarker, A., Materne, M., Vandemark, G., Shrestha, R., et al. (2011). A global survey of effects of genotype and environment on selenium concentration in lentils (Lens culinaris L.): implications for nutritional fortification strategies. Food Chem. 125, 72–76. doi: 10.1016/j.foodchem.2010.08.038
Ting, N. C., Jansen, J., Mayes, S., Massawe, F., Sambanthamurthi, R., Ooi, L. C. L., et al. (2014). High density SNP and SSR-based genetic maps of two independent oil palm hybrids. BMC Genomics 15:309. doi: 10.1186/1471-2164-15-309
Varshney, R. K., Thudi, M., Nayak, S. N., Gaur, P. M., Kashiwagi, J., Krishnamurthy, L., et al. (2014). Genetic dissection of drought tolerance in chickpea (Cicer arietinum L.). Theor. Appl. Genet. 127, 445–462. doi: 10.1007/s00122-013-2230-6
Verslues, P. E., Agarwal, M., Agarwal, S. K., Zhu, J., and Zhu, J. K. (2006). Methods and concepts in quantifying resistance to drought, salt and freezing, abiotic stresses that affect plant water status. Plant J. 45, 523–539. doi: 10.1111/j.1365-313X.2005.02593.x
Wong, M. M., Gujaria-Verma, N., Ramsay, L., Yuan, H. Y., Caron, C., Diapari, M., et al. (2015). Classification and characterization of species within the genus lens using genotyping-by-sequencing (GBS). PLoS ONE 27:e0122025. doi: 10.1371/journal.pone.0122025
Yuan, X. J., Li, X. Z., Pan, J. S., Wang, G., Jiang, S., Li, X. H., et al. (2008). Genetic linkage map construction and location of QTLs for fruit-related traits in cucumber. Plant Breed. 127, 180–188. doi: 10.1111/j.1439-0523.2007.01426.x
Zhang, Z., Hu, M., Zhang, J., Liu, D., Zheng, J., Zhang, K., et al. (2009). Construction of a comprehensive PCR-based marker linkage map and QTL mapping for fiber quality traits in upland cotton (Gossypium hirsutum L.). Mol. Breed. 24, 49–61. doi: 10.1007/s11032-009-9271-1
Keywords: lentil, ecophysiology, drought tolerance, breeding, plant, QTL, marker-assisted selection
Citation: Idrissi O, Udupa SM, De Keyser E, McGee RJ, Coyne CJ, Saha GC, Muehlbauer FJ, Van Damme P and De Riek J (2016) Identification of Quantitative Trait Loci Controlling Root and Shoot Traits Associated with Drought Tolerance in a Lentil (Lens culinaris Medik.) Recombinant Inbred Line Population. Front. Plant Sci. 7:1174. doi: 10.3389/fpls.2016.01174
Received: 31 March 2016; Accepted: 21 July 2016;
Published: 23 August 2016.
Edited by:Maria Carlota Vaz Patto, Universidade Nova de Lisboa, Portugal
Reviewed by:Maoteng Li, Huazhong University of Science and Technology, China
Rebeca Iglesias-Garcia, Nebrija University, Spain
Copyright © 2016 Idrissi, Udupa, De Keyser, McGee, Coyne, Saha, Muehlbauer, Van Damme and De Riek. 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.