Skip to main content

ORIGINAL RESEARCH article

Front. Plant Sci., 26 November 2018
Sec. Plant Breeding

Sequencing Analysis of Genetic Loci for Resistance for Late Leaf Spot and Rust in Peanut (Arachis hypogaea L.)

\r\nKenta Shirasawa*Kenta Shirasawa1*Ramesh S. BhatRamesh S. Bhat2Yogendra P. KhedikarYogendra P. Khedikar3Venkataswamy SujayVenkataswamy Sujay3Rohini M. KolekarRohini M. Kolekar2Sharanabasappa B. YeriSharanabasappa B. Yeri2Mallenahally SukruthMallenahally Sukruth2Sarvamangala CholinSarvamangala Cholin3Byregowda AshaByregowda Asha2Manish K. PandeyManish K. Pandey4Rajeev K. VarshneyRajeev K. Varshney4Makanahally V. C. GowdaMakanahally V. C. Gowda3
  • 1Department of Frontier Research and Development, Kazusa DNA Research Institute (KDRI), Chiba, Japan
  • 2Department of Biotechnology, University of Agricultural Sciences, Dharwad, India
  • 3Department of Genetics and Plant Breeding, University of Agricultural Sciences, Dharwad, India
  • 4Center of Excellence in Genomics and System Biology (CEGSB), International Crops Research Institute for the Semi-Arid Tropics (ICRISAT), Hyderabad, India

The aim of this study was to identify candidate resistance genes for late leaf spot (LLS) and rust diseases in peanut (Arachis hypogaea L.). We used a double-digest restriction-site associated DNA sequencing (ddRAD-Seq) technique based on next-generation sequencing (NGS) for genotyping analysis across the recombinant inbred lines (RILs) derived from a cross between a susceptible line, TAG 24, and a resistant line, GPBD 4. A total of 171 SNPs from the ddRAD-Seq together with 282 markers published in the previous studies were mapped on a genetic map covering 1510.1 cM. Subsequent quantitative trait locus (QTL) analysis revealed major genetic loci for LLS and rust resistance on chromosomes A02 and A03, respectively. Heterogeneous inbred family-derived near isogenic lines and the pedigree of the resistant gene donor, A. cardenasii Krapov. & W.C. Greg., including the resistant derivatives of ICGV 86855 and VG 9514 as well as GPBD 4, were employed for whole-genome resequencing analysis. The results indicated the QTL candidates for LLS and rust resistance were located in 1.4- and 2.7-Mb genome regions on A02 and A03, respectively. In these regions, four and six resistance-related genes with deleterious mutations were selected as candidates for LLS and rust resistance, respectively. These delimited genomic regions may be beneficial in breeding programs aimed at improving disease resistance and enhancing peanut productivity.

Introduction

Peanut (Arachis hypogaea L.), so-called groundnut, is an important legume crop widely cultivated for food, oil, and fodder productions. Peanut productivity in most areas is hampered by foliar diseases, particularly late leaf spots (LLS) caused by Mycosphaerella berkeleyi W.A. Jenkins, also known as Phaeoisariopsis personata (Berk. & M.A. Curtis) van Arx, and rust by Puccinia arachidis Speg. These diseases can reduce yield by up to 70% (Subrahmanyam et al., 1984) and can adversely affect kernel and fodder quality (Dwivedi et al., 2002). The production of disease resistant cultivars is limited by genetic bottlenecks, strong associations of disease resistance with poor productivity and undesirable pod characteristics, and inefficient phenotypic selection for resistance, all of which pose challenges for breeding foliar disease resistant varieties. Efficiency of selection can be enhanced using genomics-assisted breeding, which also considerably reduces the time needed for variety development.

The use of molecular breeding techniques is dependent on the availability of trait-linked markers. Efforts have been made to develop linkage maps and identify quantitative trait loci (QTLs) and flanking markers associated with LLS and rust resistance in peanut. The markers flanking these QTL were validated (Yeri et al., 2014; Sukruth et al., 2015) and employed for transferring the LLS and rust resistance regions to susceptible varieties (Varshney et al., 2014; Yeri and Bhat, 2016; Kolekar et al., 2017). Two mapping populations, TAG 24 × GPBD 4 and TG 26 × GPBD 4, were used to construct genetic maps (Khedikar et al., 2010; Sarvamangala et al., 2011; Sujay et al., 2012; Pandey et al., 2014; Kolekar et al., 2016). Initial efforts were concentrated on partial mapping with single sequence repeat (SSR) markers (Khedikar et al., 2010; Sarvamangala et al., 2011; Pandey et al., 2014), followed by map saturation with SSR and transposable element (TE) markers (Sujay et al., 2012; Kolekar et al., 2016). Then, the maps were used for QTL analysis to identify two major genomic regions associated with LLS and rust resistance (Sujay et al., 2012; Kolekar et al., 2016). The first region, on chromosome A03 (linkage group AhXV), governing LLS and rust resistance with very high percentage of variation explained (PVE; up to 82.96%). The second genomic region, on chromosome A02 (linkage group AhXII), contributed to LLS resistance only, with relatively low PVE (62.34%). In addition, another mapping population of VG 9514 × TAG 24 was also employed to detect QTLs for LLS and rust resistance on the A03 (Mondal et al., 2012). There results speculate that the genomic regions on A02 and A03 governing LLS and rust resistance are contributed by A. cardenasii Krapov. & W.C. Greg., since the disease resistant parents (GPBD 4 and VG 9514) in the three mapping populations commonly had diploid A. cardenasii Krapov. & W.C. Greg. in their pedigrees (Bertioli et al., 2016; Pandey et al., 2017).

Genome sequences are now available for the A- and B-genome diploid progenitors of peanut, namely A. duranensis Krapov. & W.C. Greg. and A. ipaeënsis Krapov. & W.C. Greg., respectively (Bertioli et al., 2016). As a result, mapping efforts are currently focused on genome sequence-based analysis (Pandey et al., 2017; Clevenger et al., 2018) and single-nucleotide polymorphism (SNP)-enriched mapping to delimit regions of interest and identify candidate genes. Pandey et al. (2017) reported 25 candidate genes for LLS resistance and nine candidates for rust on a 3.06-Mb region of the A03 by a QTL-Seq approach, in which bulks of resistance and susceptible lines of the TAG 24 × GPBD 4 population were employed. Recently, Mondal and Badigannavar (2018) have also identified five candidate resistance genes for rust on A03 using the gene function prediction/annotation data and domain search analysis.

In this study an effort was made to map the genomic regions governing LLS and rust resistance using an improved genetic map and extensive phenotypic data. To dissect the A02 and A03 chromosomal regions using genetic mapping, 190 SNPs newly developed in this study were added to the previously reported map of TAG 24 × GPBD 4 with the 263 published markers (Khedikar et al., 2010; Sarvamangala et al., 2011; Sujay et al., 2012). The updated high-density map was employed for high-resolution QTL mapping, in which phenotypic data for the LLS and rust across 11 seasons (2004–2014) were used. In addition, the chromosomal regions were validated for LLS and rust resistance by comparing SNPs from WGRS data in bulked samples of heterogeneous inbred family (HIF)-derived near isogenic lines (NILs), that genetic background is different from the RILs of TAG 24 × GPBD 4 as Pandey et al. (2017) were used for the QTL-Seq analysis. Also, co-segregation between the disease resistance phenotypes and genomic regions was examined by comparing SNPs from whole-genome resequencing (WGRS) data between resistant (GPBD 4, ICGV 86855, and VG 9514) and susceptible (TAG 24 and five Japanese lines) genotypes. This study integrated genetic and genomic analysis to identify the genic SNPs on A02 and A03 that co-segregated with LLS and rust resistance.

Materials and Methods

Plant Materials

The pedigree of the plant materials used in this study is shown in Supplementary Figure S1 and Supplementary Table S1. A TAG 24 × GPBD 4 RIL mapping population consisting of 266 RILs was used for QTL analysis. This mapping population was developed previously by crossing an LLS and rust susceptible variety, TAG 24 (Patil et al., 1995), with a resistant variety, GPBD 4 (Gowda et al., 2002), and advancing the generations by single seed decent (SSD) (Khedikar et al., 2010; Sarvamangala et al., 2011; Sujay et al., 2012). In addition, a total of 29 heterogeneous inbred family (HIF)-derived near isogenic lines (NILs) (Yeri et al., 2014), derived from crosses between TAG 24 and GPBD 4, and between another susceptible variety, TG 26 (Kale et al., 1997), and GPBD 4, were used for bulked segregant analysis by sequencing. The HIF-NILs consisted of 13 rust-resistant (1-1, 7-2, 9-2, 14-1, 50-2, 60-3, 77-1, 83-1, 89-1, 46-3, 53-1, 101-1, and 116-1) and 16 rust-susceptible (1-2, 7-1, 9-1, 9-3, 14-2, 46-1, 46-2, 50-1, 53-2, 60-1, 60-2, 77-2, 83-2, 89-2, 101-2, and 116-2) lines. The phenotyping scores of the resistance lines and susceptible were reported to be 4.55 (4.13–4.88) and 6.41 (5.38–7.75) on average, respectively (Yeri et al., 2014), where susceptible and resistance should have scores of >5 and <5.

SNP Genotyping by ddRAD-Seq Analysis

Genome-wide SNP analysis across the RIL population was conducted using ddRAD-Seq (Peterson et al., 2012). Genomic DNAs extracted from leaves of each RIL and their parental lines with DNeasy Plant Mini Kit (Qiagen, Germany) were double-digested with the restriction enzymes PstI and MspI (Thermo Fisher Scientific, MA, United States). The ddRAD-Seq libraries were constructed as described by Shirasawa et al. (2016a,b). The libraries were sequenced on HiSeq2000 (Illumina) in paired-end mode (93-base). Nucleotide sequence data were deposited in the DDBJ Sequence Read Archive (accession numbers DRA006500, DRA006514, and DRA006515).

Data processing including quality control of raw sequence reads, adaptor trimming, mapping of trimmed reads onto reference sequences, and SNP calling, was performed as described by Shirasawa et al. (2016a,b). Low-quality sequences were removed and adapters were trimmed using PRINSEQ and fastx_clipper in FASTX-Toolkit1. The filtered reads were mapped onto a genome sequence concatenated. from two diploid progenitorS of peanut, A. duranensis Krapov. & W.C. Greg. and A. ipaënsis Krapov. & W.C. Greg. (Bertioli et al., 2016), as a reference using Bowtie 2 (Langmead and Salzberg, 2012). The resultant sequence alignment–map format (SAM) files were converted to binary sequence alignment–map format (BAM) files and subjected to SNP calling using the mpileup option of SAMtools and the view option of BCFtools (Li et al., 2009). Lengths of genome regions covered with at least two reads were calculated with the genome Coverage option of BEDtools (Quinlan and Hall, 2010). High-confidence biallelic SNP candidates were selected using VCFtools (version 0.1.12b: Danecek et al., 2011) with the following criteria: (i) depth of coverage ≥5 for each data point, (ii) SNP quality score of ≥20 for each locus, and (iii) proportion of missing data of <50% for each locus.

Linkage Map Construction and QTL Analysis

Linkage analysis was performed for the TAG 24 × GPBD 4 population with JoinMap 4.0 (Van Ooijen, 2006). The “Locus Genotype Frequency” function was used to calculate Chi-square (χ2) values for each marker to test for the expected 1:1 segregation ratio Markers were placed into linkage groups with the “LOD Groupings” and “Create Groups for Mapping” command using the Kosambi mapping function (Kosambi, 1943). Calculation parameters were set for a minimum LOD (logarithm of the odds) of 3 and recombination fraction of 0.45. Marker order in groups was established using the “Calculate Map” command. The linkage map was constructed using MapChart 2.2 software (Voorrips, 2002).

Detection of QTL associated with responses to LLS and rust among the TAG 24 × GPBD 4 population was performed using a composite interval mapping (CIM) approach (Zeng, 1994) using WinQTL Cartographer, version 2.5 (Wang et al., 2011). CIM was performed using Model 6, with scanning intervals of 2.0 cM between markers and putative QTL with a window size of 10.0 cM. The number of marker cofactors for the background control was set by forward–backward stepwise regression. The “Locate QTLs” option was used automatically with a minimum of 5 cM between QTL in order to define a QTL region. One thousand permutationsaq were employed for determining the QTL using the option “permutations times” with 0.05 significance level.

Bulked Segregant Analysis by Sequencing

For bulked segregant analysis by sequencing, genomic DNAs from the 13 resistant and 16 susceptible lines of the HIF-NILs (Yeri et al., 2014) were mixed into two pools (resistant and susceptible) and subjected to library preparation (paired-end libraries with insert size of 500 bp) as described by Shirasawa et al. (2016b). The libraries were sequenced on NextSeq 500 systems (Illumina) in paired-end mode (151-base). Nucleotide sequence data were deposited in the DDBJ Sequence Read Archive (accession number DRA006500). SNP calling was modified as follows. BAM files were subjected to SNP calling using the mpileup option of SAMtools (version 0.1.19) and the mpileup2snp option of VarScan 2 (version 2.3: Koboldt et al., 2012) to obtain a variant call format (VCF) file including SNP information. High-confidence biallelic SNP candidates were selected using VCFtools (version 0.1.12b: Danecek et al., 2011) with the following criteria: (i) depth of coverage ≥10 for each data point and (ii) no missing data for each locus. The effects of SNPs on gene function were predicted using SnpEff v4.1g (Cingolani et al., 2012). From the filtered VCF files, values for variant allele frequency (FREQ) for each position were extracted from genotype fields. FREQ values over the genome were plotted using R (R Core Team, 2013). Genetic loci significantly associated with SNPs (P < 0.01) were selected under the null hypothesis of no QTLs as reported by Takagi et al. (2013) and Pandey et al. (2017).

Whole Genome Resequencing Analysis

WGRS data for ICGV 86855 and VG 9514 (resistant lines) and Chiba-handachi, Satonoka, Kintoki, Nakateyutaka, and YI-0311 (susceptible lines) were obtained from the DNA Data Bank of Japan (DDBJ) Sequence Read Archive (Shirasawa et al., 2016b; Gayathri et al., 2018; accession numbers DRA004503 and DRA0062392). Data for TAG 24 (a susceptible line) and GPBD 4 (a resistant line) were provided by The International Crops Research Institute for the Semi-Arid Tropics (Pandey et al., 2017). Data processing of WGRS reads was performed as above except for SNP filtering. High-confidence biallelic SNP candidates were selected using VCFtools (version 0.1.12b: Danecek et al., 2011) with the following criteria: (i) depth of coverage ≥5 for each data point, (ii) SNP quality score of ≥999 for each locus, and (iii) no missing data for each locus. The effects of SNPs on gene function were predicted using SnpEff v4.1g (Cingolani et al., 2012).

Results

QTL Analysis of Rust and LLS Resistance

An average of 1.7 million (M) high-quality ddRAD-Seq reads per sample were obtained from a TAG 24 × GPBD 4 RIL population (Supplementary Table S2). Reads were mapped onto genome reference sequences of A. duranensis Krapov. & W.C. Greg. and A. ipaënsis Krapov. & W.C. Greg., with an average map rate of 85.4% for TAG 24 × GPBD 4. Approximately 3.2 Mb (0.13%) of the genome was covered with at least five reads. From the read alignments, 190 SNPs were detected with high confidence. The SNPs were named as a combination of names of chromosomes and positions by linking with an underscore. These SNPs together with the 326 SSR (Sujay et al., 2012) and A. hypogaea L. transposable element (AhTE) (Kolekar et al., 2016) markers were used for a linkage analysis. The resultant map consisted of 29 linkage groups with 453 loci (171 SNPs, 89 transposons, and 193 SSRs) covering a total length of 1510.1 cM (Supplementary Table S3 and SupplementaryFigure S2). Clusters of SNPs were observed in AhXIIb (A02) and AhXV (A03). In accordance with the SNPs, 23 linkage groups were assigned to 18 chromosomes expect for A07 and A10, while the other six groups had no SNP loci. Among the 171 mapped SNPs, only three were unexpectedly located on different chromosomes (but homoeologous chromosomes), e.g., Aradu.A06_106278061 and Aradu.A06_108579782 were on AhX (B06) and Araip.B03_125319254 was on AhIIIa (A03).

This genetic map was used for QTL analysis using the phenotypic data on the LLS and rust reaction collected over 11 seasons (2004–2014) (Khedikar et al., 2010; Sujay et al., 2012; Kolekar et al., 2016). The phenotypic data has suggested high heritability and independent nature of inheritance between both the diseases (Khedikar et al., 2010). For LLS resistance, QTLs were detected on 12 linkage groups (Supplementary Table S4 and Supplementary Figure S2). Of these, two QTL regions with high LOD/PVE were detected on AhXIIb (A02) and AhXV (A03). The QTL on A02 were stable across eight seasons with a maximum additive effect 1.81 for the favorable allele contributed from GPBD 4, and PVE of up to 91.51%. For rust resistance, a single QTL region was found on AhXV (A03) with a maximum PVE of 35.72% and the additive effect of 2.06 for the favorable allele of GPBD 4 (Supplementary Table S4 and Supplementary Figure S2). The genotypes of the markers on the QTLs were highly associated with the phenotypes (Supplementary Figure S3).

Identification of Genome Fragment Conferring Rust Disease Resistance

WGRS analysis produced 315 and 319 million high-quality 151 nucleotide reads from the bulk samples of resistant and susceptible lines, respectively (Supplementary Table S2). Reads were mapped onto the reference sequence with an average alignment rate of 94.0%. The sequence reads for each sample covered approximately 1.7 Gb of the genomic regions (72.0% of the reference sequences) with at least 10 reads. Average sequence read depth was 28× for each sample. A total of 173,995 high-quality SNPs, with 67,251 and 106,744 SNPs on the A and B genomes, respectively, were detected between the two bulk samples. The expected FREQ values of the GPBD 4 alleles in the resistant and susceptible bulk-sample were 1 and 0, respectively, at the resistance loci. These predicted values were observed in a 5.2 Mb region of chromosome A03 (position 129.7–134.9 Mb) with a significant level (P < 0.01) (Figure 1), suggesting that the candidate loci for rust resistance were located in this region.

FIGURE 1
www.frontiersin.org

FIGURE 1. SNP signals detected by bulked segregant analysis by sequencing. Frequencies of GPBD 4 alleles of genome-wide SNPs in bulk samples of rust-disease susceptible lines, bulk S (A), and those of resistant lines, bulk R (B), are plotted. Differences between the bulks R and S (bulk R – bulk S) are shown in (C) with an arrow indicating a prominent peak at the end of the chromosome A03. Twenty chromosomes were painted by different colors.

High-quality WGRS reads, approximately 20–40× genome coverage, were produced for the three resistant (GPBD 4, ICGV 86855, and VG 9514) and six susceptible (TAG 24, Chiba-handachi, Nakateyutaka, YI-0311, Satonoka, and Kintoki) lines and were mapped onto the reference sequence with an alignment rate of 96.1% (Supplementary Table S2). The reads covered 2.0 Gb (83.8%) of the genome with at least five reads. From the alignment data, 3,532,384 and 168,116 SNPs between the tested lines and the reference sequences were detected on the A and B genomes, respectively. Conversely, 189,559 A-genome SNPs and 252,983 B-genome SNPs were found among the tested lines. With respect to TAG 24, there were three SNP clusters in ICGV 86855 (0–4.5 Mb on A02; 72.9–85.9 Mb on A02; and 131.6–134.6 Mb on A03), two SNP clusters in GPBD 4 (0–1.4 Mb on A02 and 131.6–134.6 Mb on A03), and two SNP clusters in VG 9514 (0–3.1 Mb on A02 and 131.9–134.6 Mb on A03) (Figure 2). No prominent signals were observed in the other five susceptible lines (Supplementary Figure S4). The positions of the SNP “hot-spots” in the resistance lines overlapped, indicating that the three resistant lines, descendants of A. cardenasii Krapov. & W.C. Greg. (GKP10017, PI262141), shared the same haplotypes in these regions. Subsequently, the candidate region for rust resistance was refined to an approximately 2.7 Mb genome fragment (131.9–134.6 Mb on A03) (Figure 3). The region contributing for LLS resistance was found at 0–1.4 Mb on A02.

FIGURE 2
www.frontiersin.org

FIGURE 2. Density of SNPs in resistant lines with respect to TAG 24. Numbers of SNPs every 100-kb length over the genome between TAG 24 and either ICGV 86855 (A), GPBD 4 (B), or VG 9514 (C) are plotted.

FIGURE 3
www.frontiersin.org

FIGURE 3. Graphical genotypes of A03 for resistant lines. Black bars show chromosome segments derived from the resistant donor, GPBD 4. Numbers on the top are chromosome positions in a mega-base scale. An arrow indicates an overlapping region of the GPDB 4 genome corresponding to the candidate region for rust resistance.

Potential Candidate Genes for Rust and LLS Resistance

The 1.4 Mb candidate region for LLS on A02 included 147 predicted genes (Supplementary Table S5). Of these, nine genes with “leucine-rich repeat (LRR)” and “mildew resistance locus O (MLO)” annotations were considered to be candidates for LLS resistance. Three of these nine genes (Aradu.0G2IC, Aradu.M8RMT, and Aradu.F918E) had sequence variations in the coding regions that caused nonsense- or missense mutations (Table 1). Nonsense mutations were seen across susceptible-line specific alleles of Aradu.0G2IC (encoding LRR and NB-ARC domain disease resistance proteins) and in resistance-line alleles of Aradu.M8RMT (disease resistance protein). Aradu.F918E (LRR and NB-ARC domain disease resistance protein) had a missense mutation at amino acid 66 that produced glutamate and lysine, respectively, in the alleles of the susceptible and resistance lines (hereafter termed GluR66LysS, with superscripts indicating alleles of resistant (R) and susceptible (S) lines).

TABLE 1
www.frontiersin.org

TABLE 1. Candidate genes for LLS and rust resistance.

The 2.7 Mb region of A03 associated with rust resistance contained 221 genes (Supplementary Table S6). Seven candidate resistance genes were selected using the annotation terms “leucine-rich repeat (LRR),” “beta-glucan,” “lipoxygenase (LOX),” and “two-component system.” Of these, while no deleterious mutation was found in one gene, 11 missense variations were predicted in the remaining six genes (Table 1): IleS27ValR in Aradu.Z87JB [putative disease resistance protein (TIR-NBS-LRR class)]; HisS307ArgR, GlyR49ArgS, and LeuR34SerS in Aradu.C88Z1 (seed linoleate 9S-lipoxygenase); CysR8TyrS in Aradu.1WV86 (glucan endo-1,3-beta-glucosidase-like protein 2-like); GluS344AlaR in Aradu.RW91L (Lipase/lipooxygenase, PLAT/LH2 family protein); LysS127GluR, ProR116LeuS, SerS72CysR, and GlyR69CysS in Aradu.NG5IQ (glucan endo-1,3-beta-glucosidase 4-like); and ArgR47SerS in Aradu.YL3ZN (receptor-like kinase 1).

Discussion

There were three points that we achieved in this study. The first is establishment of a high-density genetic map of TAG 24 × GPBD 4 (Supplementary Table S3 and Supplementary Figure S2) by employing 171 SNPs newly developed from this study, for which the ddRAD-Seq technology was indeed helpful. The number of the SNPs was lower and higher than those detected by Zhou et al. (2014) and Shirasawa et al. (2016b), respectively. This might due to differences of genetic diversity of the tested lines as reported by Shirasawa et al. (2016a), as well as data processing conditions such as criteria for filtering low quality SNPs. Out of the 171 SNPs, three (1.8%) were assigned to unexpected linkage groups, homoeologous chromosomes, as reported by Agarwa et al. (2018) where 739 of 8869 SNPs (8.3%) were mis-assigned due to the similarity of the sequences between the A and B genomes. Totally, however, this mapping result contributed to increase the number of markers on the map up to 1.5 times rather than the previous version with 289 marker loci (Kolekar et al., 2016), and marker density was also improved from 6.0 cM/locus (Kolekar et al., 2016) to 3.3 cM/locus in this study. Subsequent QTL analysis indicated the results of the previous studies (Khedikar et al., 2010; Sarvamangala et al., 2011; Sujay et al., 2012) were reproducible and promised the SNPs clustered in the QTL regions could be a source for DNA markers to select LLS and rust resistance lines in future breeding programs (Supplementary Figure S3).

Second point is that, to the best of our knowledge, this study is the first detailed genomic exploration of the QTLs on A02 and A03. NGS-based strategies were used to refine the candidate regions to a 1.4 Mb genome segment of A02 for LLS resistance and a 2.7 Mb segment of A03 for rust resistance (Figure 3 and Supplementary Table S5, S6). These segment sizes were shorter than those reported by Pandey et al. (2017), suggesting that it would be effective to use materials with novel genetic background that are expected to have new recombination breakpoints as in VG 9514 (Supplementary Table S1 and Supplementary Figure S1). Furthermore, WGRS analysis indicated that there were SNP hotspots in the QTL regions (Figure 2), that haplotype was conserved in the three resistant lines, ICGV 86855, GPBD 4, and VG 9514. Whereas we postulated that these regions might be derived from A. cardenasii Krapov. & W.C. Greg. in accordance with the pedigree (Supplementary Table S1 and Supplementary Figure S1) and the result of IntroMap (Clevenger et al., 2017), further study providing direct experimental evidences are required to confirm our speculation.

In addition, as third point, WGRS data supported the presumed effects of SNPs on gene function, which information was helpful to identify four and six candidate genes for LLS and rust resistance, respectively (Table 1). Certainly, gene annotation terms were also useful to delimit the candidates as reported by Pandey et al. (2017) and Mondal and Badigannavar (2018). Out of the candidates, three genes (Aradu.Z87JB, Aradu.1WV86, and Aradu.NG5IQ) were overlapped those suggested in the previous studies (Pandey et al., 2017; Mondal and Badigannavar, 2018). Obviously, advanced genetic analysis is required to identify the responsible sequence genes/variations for rust resistance. Because we identified a single RIL with a recombination breakpoint in the A03 candidate region (Figure 3), progeny tests of this line would be one possibility to further delimit the candidate region.

Conclusion

In conclusion, we demonstrated integration of genetic and physical mapping could facilitate identification and validation of genomic regions governing traits, and use of genotypes with different genomic backgrounds would assist in rapid identification of candidate genes. However, the physical mapping in this study was based on the reference genome sequences of two wild diploids, A. duranensis Krapov. & W.C. Greg. and A. ipaensis Krapov. & W.C. Greg. Once the genome sequence of the cultivated peanut is disclosed and used for physical mapping, any discrepancies due to chromosome-level structural variations such as translocations, inversions, copy number variations, and present/absent variations between the wild and cultivated genome sequences can be resolved. The expected results from the cultivated peanut genome as well as the current knowledge based on the wild progenitor genomes would prove valuable for transferring resistance to LLS and rust to elite varieties to enhance peanut productivity.

Author Contributions

KS and RSB conceived the study and designed the experiments. MVCG developed and provided the mapping population and GPBD 4. YPK, SC, VS, RMK, BA, and SBY conducted the phenotyping. MS made the genetic map and performed the QTL mapping. KS produced and processed the NGS data. MKP and RKV provided the sequence data for TAG 24 and GPBD 4. KS and RSB analyzed and interpreted the data and wrote the manuscript. All authors read and approved the final manuscript.

Funding

This work was supported by the DST-JSPS Bilateral Program and the Kazusa DNA Research Institute 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.

Acknowledgments

We are grateful to S. Sasamoto and C. Mimani (Kazusa DNA Research Institute) for technical assistance.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2018.01727/full#supplementary-material

Footnotes

  1. ^ http://hannonlab.cshl.edu/fastx_toolkit
  2. ^ www.ddbj.nig.ac.jp

References

Agarwa, G., Clevenger, J., Pandey, M. K., Wang, H., Shasidhar, Y., Chu, Y., et al. (2018). High-density genetic map using whole-genome resequencing for fine mapping and candidate gene discovery for disease resistance in peanut. Plant Biotechnol. J. 16, 1954–1967. doi: 10.1111/pbi.12930

PubMed Abstract | CrossRef Full Text | Google Scholar

Bertioli, D. J., Cannon, S. B., Froenicke, L., Huang, G., Farmer, A. D., Cannon, E. K., et al. (2016). The genome sequences of Arachis duranensis and Arachis ipaensis, the diploid ancestors of cultivated peanut. Nat. Genet. 48, 438–446. doi: 10.1038/ng.3517

PubMed Abstract | CrossRef Full Text | Google Scholar

Cingolani, P., Platts, A., Wang, L. L., Coon, M., Nguyen, T., Wang, L., et al. (2012). A program for annotating and predicting the effects of single nucleotide polymorphisms, SnpEff: SNPs in the genome of Drosophila melanogaster strain w1118; iso-2; iso-3. Fly 6, 80–92. doi: 10.4161/fly.19695

PubMed Abstract | CrossRef Full Text | Google Scholar

Clevenger, J., Bertioli, D. J., Leal-Bertioli, S. C. M., Chu, Y., Stalker, H. T., and Ozias-Akins, P. (2017). IntroMap: a pipeline and set of diagnostic diploid Arachis SNPs as a tool for mapping alien introgressions in Arachis hypogaea. Peanut Sci. 44, 66–73. doi: 10.3146/PS17-5.1

CrossRef Full Text | Google Scholar

Clevenger, J., Chu, Y., Chavarro, C., Botton, S., Culbreath, A., Isleib, T. G., et al. (2018). Mapping late leaf spot resistance in peanut (Arachis hypogaea) using QTL-seq reveals markers for marker-assisted selection. Front. Plant Sci. 9:83. doi: 10.3389/fpls.2018.00083

PubMed Abstract | CrossRef Full Text | Google Scholar

Danecek, P., Auton, A., Abecasis, G., Albers, C. A., Banks, E., DePristo, M. A., et al. (2011). The variant call format and VCFtools. Bioinformatics 27, 2156–2158. doi: 10.1093/bioinformatics/btr330

PubMed Abstract | CrossRef Full Text | Google Scholar

Dwivedi, S. L., Pande, S., Rao, J. N., and Nigam, S. N. (2002). Components of resistance to late leaf spot and rust among interspecific derivatives and their significance in a foliar disease resistance breeding in groundnut (Arachis hypogaea L.). Euphytica 125, 81–88. doi: 10.1023/A:1015707301659

CrossRef Full Text | Google Scholar

Gayathri, G., Shirasawa, K., Varshney, R. K., Pandey, M. K., and Bhat, R. S. (2018). Development of AhMITE1 markers through genome-wide analysis in peanut (Arachis hypogaea L.). BMC Res. Notes 11:10. doi: 10.1186/s13104-017-3121-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Gowda, M. V. C., Motagi, B. N., Naidu, G. K., Diddimani, S. B., and Sheshagiri, R. (2002). GPBD 4: a spanish bunch groundnut genotype resistant to rust and late leaf spot. Int. Arachis News Lett. 22, 29–32. doi: 10.1111/pbi.12686

PubMed Abstract | CrossRef Full Text | Google Scholar

Kale, D., Mouli, C., Murty, G., and Rao, M. (1997). Development of a new groundnut variety, TG-26 by using induced mutations in cross breeding. Mutat. Breed. Newslet. 43, 25–27.

Khedikar, Y. P., Gowda, M. V., Sarvamangala, C., Patgar, K. V., Upadhyaya, H. D., and Varshney, R. K. (2010). A QTL study on late leaf spot and rust revealed one major QTL for molecular breeding for rust resistance in groundnut (Arachis hypogaea L.). Theor. Appl. Genet. 121, 971–984. doi: 10.1007/s00122-010-1366-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Koboldt, D. C., Zhang, Q., Larson, D. E., Shen, D., McLellan, M. D., Lin, L., et al. (2012). VarScan 2: somatic mutation and copy number alteration discovery in cancer by exome sequencing. Genome Res. 22, 568–576. doi: 10.1101/gr.129684.111

PubMed Abstract | CrossRef Full Text | Google Scholar

Kolekar, R. M., Sujay, V., Shirasawa, K., Sukruth, M., Khedikar, Y. P., Gowda, M. V. C., et al. (2016). QTL mapping for late leaf spot and rust resistance using an improved genetic map and extensive phenotypic data on a recombinant inbred line population in peanut (Arachis hypogaea L.). Euphytica 209, 147–156. doi: 10.1007/s10681-016-1651-0

CrossRef Full Text | Google Scholar

Kolekar, R. M., Sukruth, M., Shirasawa, K., Nadaf, H. L., Motagi, B. N., Lingaraju, S., et al. (2017). Marker-assisted backcrossing to develop foliar disease-resistant genotypes in TMV 2 variety of peanut (Arachis hypogaea L.). Plant Breed 136, 948–953. doi: 10.1111/pbr.12549

CrossRef Full Text | Google Scholar

Kosambi, D. D. (1943). The estimation of map distances from recombination values. Ann. Eugen. 12, 172–175. doi: 10.1111/j.1469-1809.1943.tb02321.x

CrossRef Full Text | Google Scholar

Langmead, B., and Salzberg, S. L. (2012). Fast gapped-read alignment with bowtie 2. Nat. Methods 9, 357–359. doi: 10.1038/nmeth.1923

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, H., Handsaker, B., Wysoker, A., Fennell, T., Ruan, J., Homer, N., et al. (2009). The sequence alignment/map format and SAMtools. Bioinformatics 25, 2078–2079. doi: 10.1093/bioinformatics/btp352

PubMed Abstract | CrossRef Full Text | Google Scholar

Mondal, S., and Badigannavar, A. M. (2018). Mapping of a dominant rust resistance gene revealed two R genes around the major Rust_QTL in cultivated peanut (Arachis hypogaea L.). Theor. Appl. Genet. doi: 10.1007/s00122-018-3106-6 [Epub ahead of print].

PubMed Abstract | CrossRef Full Text | Google Scholar

Mondal, S., Badigannavar, A. M., and D’Souza, S. F. (2012). Molecular tagging of a rust resistance gene in cultivated groundnut (Arachis hypogaea L.) introgressed from Arachis cardenasii. Mol. Breed. 29, 467–476. doi: 10.1007/s11032-011-9564-z

CrossRef Full Text | Google Scholar

Pandey, M. K., Khan, A. W., Singh, V. K., Vishwakarma, M. K., Shasidhar, Y., Kumar, V., et al. (2017). QTL-seq approach identified genomic regions and diagnostic markers for rust and late leaf spot resistance in groundnut (Arachis hypogaea L.). Plant Biotechnol. J. 15, 927–941. doi: 10.1111/pbi.12686

PubMed Abstract | CrossRef Full Text | Google Scholar

Pandey, M. K., Upadhyaya, H. D., Rathore, A., Vadez, V., Sheshshayee, M. S., Sriswathi, M., et al. (2014). Genomewide association studies for 50 agronomic traits in peanut using the reference set comprising 300 genotypes from 48 countries of semi-arid tropics of the world. PLoS One 9:e105228. doi: 10.1371/journal.pone.0105228

PubMed Abstract | CrossRef Full Text | Google Scholar

Patil, S., Kale, D., Deshmukh, S., Fulzele, G., and Weginwar, B. (1995). Semi-dwarf, early maturing and high yielding new groundnut variety. TAG-24. J. Oilseeds Res. 12, 254–257.

Google Scholar

Peterson, B. K., Weber, J. K., Kay, E. H., Fisher, H. S., and Hoekstra, H. E. (2012). Double digest RADseq: an inexpensive method for de novo SNP discovery and genotyping in model and non-model species. PLoS One 7:e37135. doi: 10.1371/journal.pone.0037135

PubMed Abstract | CrossRef Full Text | Google Scholar

Quinlan, A. R., and Hall, I. M. (2010). BEDTools: a flexible suite of utilities for comparing genomic features. Bioinformatics 26, 841–842. doi: 10.1093/bioinformatics/btq033

PubMed Abstract | CrossRef Full Text | Google Scholar

R Core Team (2013). R: A Language and Environment for Statistical Computing. Vienna: R Foundation for Statistical Computing.

Google Scholar

Sarvamangala, C., Gowda, M. V. C., and Varshney, R. K. (2011). Identification of quantitative trait loci for protein content, oil content and oil quality for groundnut (Arachis hypogaea L.). Field Crops Res. 122, 49–59. doi: 10.3389/fpls.2013.00023

PubMed Abstract | CrossRef Full Text | Google Scholar

Shirasawa, K., Hirakawa, H., and Isobe, S. (2016a). Analytical workflow of double-digest restriction site-associated DNA sequencing based on empirical and in silico optimization in tomato. DNA Res. 23, 145–153. doi: 10.1093/dnares/dsw004

PubMed Abstract | CrossRef Full Text | Google Scholar

Shirasawa, K., Kuwata, C., Watanabe, M., Fukami, M., Hirakawa, H., and Isobe, S. (2016b). Target amplicon sequencing for genotyping genome-wide single nucleotide polymorphisms identified by whole-genome resequencing in peanut. Plant Genome 9. doi: 10.3835/plantgenome2016.06.0052

PubMed Abstract | CrossRef Full Text | Google Scholar

Subrahmanyam, P., Williams, J., McDonald, D., and Gibbons, R. (1984). The influence of foliar diseases and their control by selective fungicides on a range of groundnut (Arachis hypogaea L.) genotypes. Ann. Appl. Biol. 104, 467–476. doi: 10.1111/j.1744-7348.1984.tb03029.x

CrossRef Full Text | Google Scholar

Sujay, V., Gowda, M. V., Pandey, M. K., Bhat, R. S., Khedikar, Y. P., Nadaf, H. L., et al. (2012). Quantitative trait locus analysis and construction of consensus genetic map for foliar disease resistance based on two recombinant inbred line populations in cultivated groundnut (Arachis hypogaea L.). Mol. Breed. 30, 773–788. doi: 10.1007/s11032-011-9661-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Sukruth, M., Paratwagh, S. A., Sujay, V., Kumari, V., Gowda, M. V. C., Nadaf, H. L., et al. (2015). Validation of markers linked to late leaf spot and rust resistance, and selection of superior genotypes among diverse recombinant inbred lines and backcross lines in peanut (Arachis hypogaea L.). Euphytica 204, 343–351. doi: 10.1007/s10681-014-1339-2

CrossRef Full Text | Google Scholar

Takagi, H., Abe, A., Yoshida, K., Kosugi, S., Natsume, S., Mitsuoka, C., et al. (2013). QTL-seq: rapid mapping of quantitative trait loci in rice by whole genome resequencing of DNA from two bulked populations. Plant J. 74, 174–183. doi: 10.1111/tpj.12105

PubMed Abstract | CrossRef Full Text | Google Scholar

Van Ooijen, J. W. (2006). JoinMap®4, Software for the Calculation of Genetic Linkage Maps in Experimental Populations. Wageningen: Kyazma BV.

Google Scholar

Varshney, R. K., Pandey, M. K., Janila, P., Nigam, S. N., Sudini, H., Gowda, M. V., et al. (2014). Marker-assisted introgression of a QTL region to improve rust resistance in three elite and popular varieties of peanut (Arachis hypogaea L.). Theor. Appl. Genet. 127, 1771–1781. doi: 10.1007/s00122-014-2338-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Voorrips, R. E. (2002). MapChart: software for the graphical presentation of linkage maps and QTLs. J. Hered. 93, 77–78. doi: 10.1093/jhered/93.1.77

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, S., Basten, C. J., and Zeng, Z. B. (2011). Windows QTL Cartographer 2.5. Raleigh, NC: North Carolina State University.

Google Scholar

Yeri, S. B., and Bhat, R. S. (2016). Development of late leaf spot and rust resistant backcross lines in Jl 24 variety of groundnut (Arachis hypogaea L.). Electron J. Plant Breed. 7, 37–41. doi: 10.5958/0975-928X.2016.00005.3

CrossRef Full Text | Google Scholar

Yeri, S. B., Shirasawa, K., Pandey, M. K., Gowda, M. V. C., Sujay, V., Shriswathi, M., et al. (2014). Development of NILs from heterogeneous inbred families for validating the rust resistance QTL in peanut (Arachis hypogaea L.). Plant Breed. 133, 80–85. doi: 10.1111/pbr.12130

CrossRef Full Text | Google Scholar

Zeng, Z. B. (1994). Precision mapping of quantitative trait loci. Genetics 136, 1457–1468.

Google Scholar

Zhou, X., Xia, Y., Ren, X., Chen, Y., Huang, L., Huang, S., et al. (2014). Construction of a SNP-based genetic linkage map in cultivated peanut based on large scale marker development using next-generation double-digest restriction-site-associated DNA sequencing (ddRADseq). BMC Genomics 15:351. doi: 10.1186/1471-2164-15-351

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: late leaf spot and rust diseases, peanut, quantitative trait locus, restriction-site associated DNA sequencing, whole-genome resequencing analysis

Citation: Shirasawa K, Bhat RS, Khedikar YP, Sujay V, Kolekar RM, Yeri SB, Sukruth M, Cholin S, Asha B, Pandey MK, Varshney RK and Gowda MVC (2018) Sequencing Analysis of Genetic Loci for Resistance for Late Leaf Spot and Rust in Peanut (Arachis hypogaea L.). Front. Plant Sci. 9:1727. doi: 10.3389/fpls.2018.01727

Received: 03 August 2018; Accepted: 07 November 2018;
Published: 26 November 2018.

Edited by:

Petr Smýkal, Palacký University Olomouc, Czechia

Reviewed by:

Harsh Raman, New South Wales Department of Primary Industries, Australia
Amit A. Deokar, University of Saskatchewan, Canada

Copyright © 2018 Shirasawa, Bhat, Khedikar, Sujay, Kolekar, Yeri, Sukruth, Cholin, Asha, Pandey, Varshney and Gowda. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Kenta Shirasawa, shirasaw@kazusa.or.jp

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.