Identification of major QTLs underlying tomato spotted wilt virus resistance in peanut cultivar Florida-EPTM ‘113’
BMC Genetics volume 17, Article number: 128 (2016)
Spotted wilt caused by tomato spotted wilt virus (TSWV) is one of the major peanut (Arachis hypogaea L.) diseases in the southeastern United States. Occurrence, severity, and symptoms of spotted wilt disease are highly variable from season to season, making it difficult to efficiently evaluate breeding populations for resistance. Molecular markers linked to spotted wilt resistance could overcome this problem and allow selection of resistant lines regardless of environmental conditions. Florida-EPTM ‘113’ is a spotted wilt resistant cultivar with a significantly lower infection frequency. However, the genetic basis is still unknown. The objective of this study is to map the major quantitative trait loci (QTLs) linked to spotted wilt resistance in Florida-EPTM ‘113’.
Among 2,431 SSR markers located across the whole peanut genome screened between the two parental lines, 329 were polymorphic. Those polymorphic markers were used to further genotype a representative set of individuals in a segregating population. Only polymorphic markers on chromosome A01 showed co-segregation between genotype and phenotype. Genotyping by sequencing (GBS) of the representative set of individuals in the segregating population also depicted a strong association between several SNPs on chromosome A01 and the trait, indicating a major QTL on chromosome A01. Therefore marker density was enriched on the A01 chromosome. A linkage map with 23 makers on chromosome A01 was constructed, showing collinearity with the physical map. Combined with phenotypic data, a major QTL flanked by marker AHGS4584 and GM672 was identified on chromosome A01, with up to 22.7 % PVE and 9.0 LOD value.
A major QTL controlling the spotted wilt resistance in Florida-EPTM ‘113’ was identified. The resistance is most likely contributed by PI 576638, a hirsuta botanical-type line, introduced from Mexico with spotted wilt resistance. The flanking markers of this QTL can be used for further fine mapping and marker assisted selection in peanut breeding programs.
Cultivated Peanut (Arachis hypogaea L, 2n = 4x = 40) is a member of the legume family and is cultivated mainly in semi-arid tropic and sub-tropic regions. Worldwide, more than 42 million tons of peanuts were produced in 2013. Peanut seeds have high concentrations of oil (45–52 %) and protein (about 25 %), serving as a principal nutrition source in some developing countries . Cultivated peanut is an allotetraploid with a genome composition of AABB. The genomes A and B were most likely derived from two wild diploid Arachis species, A. duranensis (A genome, 2n = 20) and A. ipaensis (B genome, 2n = 20) respectively . Peanut production is significantly affected by spotted wilt disease caused by Tomato spotted wilt virus (TSWV) (genus Tospovirus, family Bunyaviridae), specifically in the United States. In 1997, the production losses due to spotted wilt in Georgia were estimated to be approximately $40 million . The typical symptoms of spotted wilt on peanuts are yellowing, stunting, concentric ringspots, chlorosis, and necrosis of various sizes and shapes on leaflets .
Many factors affect the severity of spotted wilt including peanut variety, planting date, plant population, row pattern, crop rotation, and tillage. Host resistance (peanut variety) is the most important factor to reduce disease risk. Hence, the development of spotted wilt resistance varieties has become a major breeding objective in peanut breeding programs in the United States. Traditionally, peanut breeders select resistant plants in the field under natural conditions. However, expression of spotted wilt disease is highly variable from season to season, reducing selection efficiency. The recent advances in genetic and genomic tools and resources for peanut increased the potential of using molecular markers for selection to accelerate the peanut cultivar improvement . Much progress was made in the past few years. Specifically, the two ancestor genomes, the A genome from A. duranensis and the B genome from A. ipaensis, have been sequenced and annotated , which provided a fundamental resource for molecular marker development.
Simple sequence repeat (SSR) marker is a valuable type of marker. Its abundance in the genome, co-dominance, multiple alleles, high polymorphism, PCR-based simple analysis, and transferability from other related species made SSR marker a great choice [7, 8] for peanut genetic studies. Single nucleotide polymorphism (SNP) is another type of molecular marker, which is very abundant and distributed throughout the whole genome. Those molecular markers can be used to help construct peanut linkage maps [9–11]. Linkage maps provide the basic framework for genetic and genomics studies, such as quantitative trait locus (QTL) analysis, marker assisted selection (MAS), comparative genomics, and genome assembly. QTL analysis based on the linkage map is critical in identifying markers linked to agronomically important traits. Several QTLs for important traits have been identified in peanut, such as drought tolerance, disease resistance, and nutritional quality [9, 12–15].
Florida-EPTM ‘113’ is a new runner-type variety with superior spotted wilt resistance . It has been tested under earlier planting date (April) and reduced seed density (13.1 seed per meter). Both conditions favor for spotted wilt epidemics. However, Florida-EPTM ‘113’proved to have resistance sufficient to obviate high risk situations presented by earlier planting date and lower seed density . One of the parental lines of Florida-EPTM ‘113’ is NC94022. High level of field resistance has been reported for NC94022  and the crosses utilizing NC94022 as a parental line in peanut breeding programs were initiated . In 2012, Qin et al.  reported the first genetic linkage map based on NC94022-derived population and a major QTL asssociated with resistance to TSWV on linkage group ‘A1’ was reported, which had a phenotypic variation explained (PVE) of 35.8 %. Recently Khera et al.  reported an improved genetic map for the same NC94022-derived population and QTLs on chromosome ‘A01’ associated with multi-year TSWV phenotypic data were identified. However, whether genetic basis of the spotted wilt resistance in Florida-EPTM ‘113’ is the same as the NC94022 was unknown. The objective of this study was to understand the genetic basis of the supurior spotted wilt resistance in Florida-EPTM ‘113’ through maping and QTL analysis using Florida-EPTM ‘113’-derived populations. If the major QTLs of the resistance in Florida-EPTM-‘113’ can be detected, MAS can be applied in the breeding program to avoid uncertain environmental impacts on selection of new cultivars.
An F2 population was derived from cross between Florida-EPTM ‘113’ and Georgia Valencia made in 2009. Georgia Valencia is a large-podded Valencia market type peanut susceptible to spotted wilt . The F2 segregating population comprised of 200 lines was planted at the North Florida Research and Education Center (NFREC) near Marianna, FL, in 2011. All F2 plants were self-pollinated to generate next generation seeds. Since seeds harvested from a few F2 lines were mixed up during seeds processing and thus discarded, only 163 pure F2:3 (F2 -derived in F3) families were kept and planted in 2012. Planting of F2:3 lines took place at Plant Science Research and Education Unit (PSREU) located near Citra, Florida in early April and at NFREC in the middle of April 2012. An augmented experimental design with two parental lines as controls was used at each location. Each plot had two rows, 0.9 m wide and 4.5 m long, and a single family was planted in each plot. Seed planting density was one seed per 0.3 m.
The plants from each F2:3 family allowed self-pollinating. The F2:3 seeds were bulk harvested and seeds from each family were randomly chosen to generate 163 F2:4 (F2 -derived in F4) families in 2013. After the harvest of the F2:4 families, 163 F2:5 (F2 -derived in F5) families were planted subsequently in 2014. The experimental design for the F2:5 was the same as the F2:3 and F2:4 but the planting took place only at the NFREC location (Additional file 1). A single plant was harvested from each F2:5 families and seeds were randomly chosen to plant the F6 generation as RILs (163 lines) with a high level of homozygosity. F2:3, F2:4 and F2:5 generations were included in the research for phenotyping.
Rating for disease resistance
Two different disease evaluation methods were conducted to assess the severity of spotted wilt. One was a visual rating on a scale ranging from 1 to 10, and the other was a form of immunoassay (immunostrip testing), which was used to test for the presence of TSWV in the root crown. A visual rating was conducted on a whole plot basis prior to digging. Each plot was assessed for typical symptoms of spotted wilt such as stunting and foliar symptoms of ringspot, leaf necrosis, and chlorosis (yellowing) . The 1 to 10 scale represented a percentage of infected plants (1, 1.5, 2, 3, 4, 5, 6, 7, 8, 9, 10 equals 0 %, 1–10 %, 11–20 %, 21–30 %, 31–40 %, 41–50 %, 51–60 %, 61–70 %71–80 %, 81–90 %, and 91–100 %, respectively).
The immunostrip testing was conducted by using the ImmunoStrip Kit (Agdia Inc., Elkhart, IN, USA). Ten individual plants were randomly harvested from each plot and root crowns of each plant were collected and air-dried after digging. Approximately 0.4 g root crown sample was trimmed and subjected to immunostrip test following the manufacturer’s instructions. The percentage of infected plants out of the 10 randomly selected plants from each plot represented the TSWV infection frequency of each plot, thus each plot had the possibility of TSWV infection percentage ranging from 0 to 100 %.
The visual disease evaluation method was conducted in F2:3 and F2:4 families at both NFREC and PSREU. For F2:5 families, only at NFREC, the evaluation was applied. Immunostrip testing was utilized in F2:3 and F2:4 generations at only NFREC (Additional file 1). Phenotypic correlations were analyzed using SPSS 22.0 to calculate Spearman’s rank correlation coefficient (Spearman’s rho).
The genomic DNA of each F2 individual was extracted from approximately 500 mg young leaf tissues using the method described by Dellaporta et al. , with a modification that 1 % polyvinylpyrrolidone (PVP) was added to extraction buffer to remove phenolic compounds. DNA quality and quantity were evaluated by using 1 % agarose gel electrophoresis and NanoDrop (Thermo Scientific, United States). The isolated DNA was diluted to 5 to 30 ng/ul for further polymerase chain reaction (PCR) process.
Publicly available SSR markers were selected based on polymorphism information content (PIC) and linkage group information from related literatures [7–9, 11, 13, 14, 20, 24–39] (Additional file 2). The primers of selected SSR markers were synthesized by InvitrogenTM, Life Technologies. PCR was performed in 10 ul volumes containing 1 μl of 10 × PCR buffer, 1 μl of Magnesium Chloride (25 mM), 1 μl of dNTP (2 mM), 0.3 μl of Taq enzyme (home-made), 2 μl of forward and reverse primers (2 mM), 2 μl of DNA template (10 ng/μl), and 2.8 μl of distilled deionized water. The PCR program was operated using a touchdown program with an initial denaturation at 94 °C for 3 min; 10 cycles of amplification at 94 °C for 30 s, 65 to 55 °C for 20 s (every cycle drops one degree until 55 °C), 72 °C for 40 s; 30 cycles of amplification at 94 °C for 30 s, 55 °C for 20 s, 72 °C for 40 s; and a final extension at 72 °C for 7 min. The PCR products were separated on 6 % non-denatured polyacrylamide gel electrophoresis (PAGE) under 150 volts for 2 h in 1X TBE buffer with DYCZ- 30B gel rigs system (Beijing, China) . The gels were stained by ethidium bromide and visualized under UV light.
Genotyping by sequencing (GBS)
Two parental lines (Florida-EPTM ‘113’ and Georgia Valencia) plus 10 F6 lines samples were selected for GBS. Among the 10 F6 lines, five lines consistently showed disease resistance for 3 years and the other five lines showed susceptibility. The 12 DNA samples with high quality and a concentration higher than 50 ng/ul were delivered to Cornell University, Institute of biotechnology (Genomic Diversity Facility) for 96-plex GBS library preparation (in combination with other 83 samples unrelated with this experiment) and sequencing with Illumina HiSeq 2000 platform following the optimized protocol . The restriction enzyme, ApeKI was used for reduced representation library preparation.
The Tassel-GBS pipeline was utilized for SNP calling . Trimmed and cleaned sequence tags were aligned to peanut A and B genomes respectively  using Bowtie2 . The minor allele frequency (MAF) cutoff was set at 0.01. Raw SNPs were called and the resulting HapMap.hmp.txt files were input into TASSEL software  for further filtering. SNPs with more than 25 % missing data were removed. The remaining SNPs were applied for single marker analysis first. Then association mapping analysis was conducted. The general linear model (GLM) was performed and a Manhattan plot was generated by Tassel v5.0 after inputting phenotyping data of 12 samples [44, 45].
Linkage and QTL analysis
Linkage analysis was performed using software QTL IciMapping V4.0  in combination with JoinMap 4.0 . To construct linkage groups, a minimum log-of-odds (LOD) threshold of 3.0 was applied and map distances were converted to centi morgans (cMs) using Kosambi mapping function . The 3-year phenotyping data were incorporated with linkage map information for QTL mapping using the software QTL IciMapping V4.0 . The Inclusive Composite Interval Mapping (ICIM-ADD) method [49, 50] was applied for QTL analysis with a minimum 3.0 LOD, 0.001 probability in stepwise regression, and a scanning interval of 1.0 cM/step.
Disease rating distribution in the segregating populations
Out of 163 F2:3 lines, more than 68 (40 %) at the PSREU in 2012 had a disease rating of “1”, indicating no visible symptoms. Only 25 lines showed more than 20 % disease symptoms (disease rating more than 2) (Fig. 1a). The two checks (two parental lines) showed a distinct difference at PSREU in 2012. Out of the 10 Florida-EPTM ‘113’ plots, eight did not show symptoms, but all Georgia Valencia plots showed various levels of disease symptoms with ratings ranging from 2 to 8 (Fig. 1b). Compared to the rating data collected from PSREU, data from NFREC had more plots with ratings of 1.5 and less plots with ratings of 1 (Fig. 1c). Florida-EPTM ‘113’ checks had mostly ratings of 1, while Georgia Valencia had a wide rating range (Fig. 1d). In general, the spotted wilt epidemics detected by visual rating was higher at NFREC than at PSREU.
The distribution of Immunostrip test results in 2012 was relatively even in contrast to the skewed distribution in visual ratings (Fig. 2a). The disease rating distributions of Florida-EPTM ‘113” and Georgia Valencia were distinct in 2012 (Fig. 2b).
The distribution of visual ratings of the F2:4 population at PSREU in 2013 were skewed toward low disease rating. There were 132 plots, whose rating range were from 1, 1.5 to 2, 60 plots were scored 1.5; 34 plots were scored 1, and 38 plots were scored 2 (Additional file 3). Florida-EPTM ‘113’ received all ratings below 2; however, Georgia Valencia also received many rating scores below 2 (Additional file 3). At NFREC, the visual rating resulted in most plots being scored 2 (45 plots), followed by rating score of 3 (38 plots). Only one plot was free of visual symptoms, whereas more plots had a high rating score (>3) (Additional file 3). Georgia Valencia plots had a wide range of disease ratings ranging from 1 to 9 indicating the uneven infection or disease development in the field (Additional file 3). In general, the disease pressure was higher in 2013 than in 2012, and the disease rating was higher at NFREC than at PSREU.
The distribution of immunostrip test results of the population in 2013 was nearly a normal distribution with a majority of plots falling between 40 and 70 % (Fig. 2c). As for the disease severity distribution of the two checks, it was bimodal, which represented their resistant and susceptible feature respectively (Fig. 2d).
The F2:5 population was only rated visually at NFREC, 2014. More lines (38, 35, and 33 lines, respectively) received rating scores of 2, 3 and 4 than other lower or higher rating scores (Additional file 4). Most plots of Florida-EPTM ‘113’ received rating scores of 2 and lower, and the rating scores of Georgia Valencia plots were all higher than 2 (Additional file 4).
A total of seven phenotypic datasets representing different years, locations, generations of the population, and measurement methods were recorded. The seven datasets were entitled as 2012NF-VR (2012: year 2012; NF: NFREC; VR: visual rating), 2013NF-VR (2013: year 2013; NF: NFREC; VR: visual rating), 2014NF-VR (2014: year 2014; NF: NFREC; VR: visual rating), 2012PS-VR (2012: year 2012; PS: PSREU VR: visual rating), 2013PS-VR (2013: year 2013; PS: PSREU; VR: visual rating), 2012NF-IS (2012: year 2012; NF: NFREC; IS: immunostrip), and 2013NF-IS (2013: year 2013; NF: NFREC; IS: immunostrip). There were a total of 21 possible pair-wise combinations among the seven datasets. The correlations of all the combinations were significant, except the combination between 2012NF-VR and 2013PS-VR (Table 1).
The highest correlation was found between 2013NF-IS and 2013NF-VR with a coefficient of 0.78 (p < 0.01). The second highest correlation coefficient was 0.59 (p < 0.01) between 2012NF-IS and 2012NF-VR. In general, the datasets from the same location tended to have higher correlations than those from different locations, indicating environmental variation between the two locations. Even though 2012NF-IS and 2013NF-IS were datasets of different years and generations, they were still correlated with a coefficient of 0.498 (p < 0.01). However, 2012NF-IS and 2013NF-VR were also correlated with a coefficient of 0.503 (p < 0.01), but they were datasets from different years and using different measurement methods. Overall, the majority of the correlations had coefficients in a range between 0.2 and 0.4. While comparing the correlations within the year of 2013 (0.781, 0.401, and 0.316, p < 0.01) and within the year of 2012 (0.591, 0.270, and 0.266, p < 0.01), the correlations of datasets collected in 2013 were higher than those in 2012.
Within the correlations of datasets between 2012 and 2013 that used the same measurement method (visual rating), the dataset generated from NFREC had greater correlation coefficient (0.355, p < 0.01) than that from PSREU (0.236, p < 0.01). While comparing the correlations between datasets in 2012 and 2013 at the same location (NFREC), the correlation coefficient of the IS (0.498, p < 0.01) was higher than that of the VR (0.355, p < 0.01).
SSR marker screening
A total of 2,431 markers across the whole peanut genomes (Table 2) were screened between the two parental lines, Florida-EPTM ‘113’ and Georgia Valencia. Each linkage group (LG) had an average of 88.4 markers screened, ranging from 62 to 128 markers/LG. The average number of markers screened was 93.9 on each A chromosome and 82.9 on each B chromosome, with additional 663 markers having no linkage group information. The sequences from 2,431 markers were aligned to A and B reference genomes  and 1,637 (67.34 %) could be mapped to the genomes and 794 SSR primers (32.66 %) could not be mapped, which may due to the genome rearrangement of cultivated peanut or incompleteness of the reference genomes.
The number of amplifiable markers was 2,221 with an amplification ratio of 91.36 %. Polymorphism was detected between the two parental lines at 329 SSR marker loci with a polymorphic ratio of 13.51 %. The lowest and highest polymorphic ratio was found on A04 (2.94 %) and B02 (26.39 %), respectively (Table 2).
Major QTL identification
To identify the chromosome location of the major QTLs for TSWV resistance, the 329 polymorphic SSR markers were further screened using 12 selected individuals, consisting of two parental lines, five susceptible and five resistant lines. Only 19 markers on chromosome A01 displayed co-segregation with phenotypes (Additional file 5). Markers on the other 19 chromosomes showed random band patterns and did not correspond to their phenotypic data (Additional file 5) indicating that potential major QTLs of spotted wilt disease resistance in Florida-EPTM ‘113’were on the chromosome A01. In addition, the GBS experiment generated a total of 1,856,429 sequences reads and 7,972 raw SNPs were called by Tassel-GBS pipeline. After removing the SNPs with high level of missing data (>25 %), a total of 2,670 SNPs were identified. Single marker analysis indicated that 18 SNP markers co-segregated with phenotypes of the 12 individuals including 8 SNPs on A01 chromosome, two on B7 chromosome one marker on A3, A4, A6, A7, A8, A9, B8 and B9 chromosome respectively. Association mapping further confirmed some of the markers associated with the TSWV resistance as showed on the manhattan plot (Fig. 3). The results further indicated that there were major loci on chromosome A01 associated with spotted wilt disease resistance (Fig. 3). Chromosomes A3, A9, B8 and B9 also showed a few SNPs associated with TSWV resistance. However, compared to the number of associated SNPs on chromosome A01, the numbers of SNPs on other chromosomes were relatively few.
To increase the mapping resolution on chromosome A01, an additional 154 SSR markers from latest literature [9, 11, 34] on chromosome A01 were utilized for polymorphism screening using the parental lines. Overall, 2,583 markers have been screened (Additional file 2) and 29 polymorphic SSR markers on A01 were used to genotype the whole F2 population of 163 individuals.
Linkage analysis revealed that out of 29 polymorphic SSR markers 23 markers were mapped on one linkage group. The genetic distance was 157.80 cM (Fig. 4b). Out of 23 linked SSR markers, 19 markers can be aligned to the reference genome on chromosome A01. The physical positions of seven SNP markers from GBS were also indicated on the physical map (Fig. 4a). On the physical map, the top marker aligned to A01 was AhTE0369 and the position was at 5.7 megabase. The bottom marker was AHGS1351 and the position was at 105.1 megabase. A high collinearity between the genetic and physical maps was observed (Fig. 4, Table 3).
The seven phenotypic datasets and the genotypic results were used for QTL analysis. Two QTLs representing two locations, NFREC and PSREU, were detected on chromosome A01 (Fig. 5). The same QTL was identified using phenotypic datasets of 2012NF-IS, 2013NF-VS, 2013NF-IS, and 2014NF-VS with flanking markers, AHGS4584 (80.73 cM) and GM672 (83.28 cM). A similar QTL was identified using 2012NF-VS dataset with flanking markers, AHGS1646 (78.83 cM) and AHGS4584 (80.73 cM). The QTL using 2013NF-IS dataset showed the highest LOD score (9.00) and PVE (22.7 %) (Table 4). QTL identified using datasets of 2012PS-VS and 2013PS-VS had same flanking markers, AHGS1713 (90.00 cM) and AHGS1760 (92.57 cM). However, the QTL identified using 2013PS-VS dataset showed the lowest LOD score (3.76) and PVE (10.02 %) (Table 4).
Phenotypic variability in population
Mechanical transmission of TSWV is difficult to achieve. Several factors affect the transmission efficiency . In addition, it is time consuming to inoculate every single plant in large-scale breeding experiment. In this study, the phenotypic data were collected based on natural inoculation. Under natural inoculation conditions, single plant rating is not reliable for spotted wilt disease evaluation, so the plot-based evaluation with multiple years at two different locations was applied in this study. Each line in the F2:4 and F2:3 families can be traced back to a single F2 plant. The family mean of phenotypic value from multiple individuals can represent the phenotype of the single F2 plant .
Visual rating was a comprehensive plot-based method including the assessment of spotted wilt severity and incidence. However, compared to immunostrip test, visual rating had a narrow disease scoring scale (1 to 10). When conducting QTL analysis, the discrete scales might not be able to represent the diverse variability of the populations. In addition, visual rating relies on experience and thus is relative subjective and may produce rating biases by a single person.
Immunostrip results were derived from 10 random plants in a plot and can be considered as the average of a plot. Immunostrip testing can capture TSWV reaction more accurately than visual scoring, since it can detect the virus in non-symptomatic plants. Approximately 50 % of non-symptomatic plants of some cultivars were found to be infected with TSWV . This means that visual rating does not capture the potential of disease expression. Visual rating resulted more ratings in the low range (rating 1, 1.5 and 2) and much fewer in the high range (rating 7, 8 and 9) mostly due to the non-symptomatic issue and non-synchronized infection. Varying disease pressure over seasons and locations can affect the spotted wilt epidemic, making visual ratings less reliable. Even the homozygous parental lines exhibited inconsistent disease levels across the field, specifically when disease pressure was not sufficient; the susceptible and resistant checks were not easily distinguished. However, the distribution of immunostrip testing results of the two parents showed a distinction. According to immunostrip results, more plots had severe infection compared to those with no or minor infection.
Major QTL underlying the TSWV resistance in Florida-EPTM ‘113’
Bulk segregant analysis based on polymorphic SSR markers indicated that only markers located on A01 chromosome showed co-segregation with resistant/susceptible plants (Additional file 5). In addition, according to association analysis based on GBS data (Fig. 5), 11 SNPs showed a high degree of association with spotted wilt resistance (−log10 (p-value) = 8 as cutoff) and 7 of 11 SNPs were located on chromosome A01. The other four were located at chromosome A3, A9, B8, and B9 respectively, most likely representing minor QTLs, which will be investigated in the future. In this study, we focused on the major QTLs located at chromosome A01.
The QTL analysis in this study (Table 4 and Fig. 5) revealed two putative QTLs on chromosome A01 next to each other. One was associated with NFREC phenotypic data and the other was associated with PSREU phenotypic data, which most likely are the same QTL with a slight shift to different positions due to the environment effects . In addition, genotype-by-environment interaction, population size, marker density, and genotyping errors  can also influence the QTL position. Therefore, most likely one major QTL on chromosome A01 with 23 % PVE controls the spotted wilt disease resistance in Florida-EPTM ‘113’. Association analysis based on GBS data revealed seven SNPs on chromosome A01 associated with TSWV resistance with three of them between flanking markers AHGS 4584 and GM 672. Two are between markers AHGS 1646 and AHGS4584 (Fig. 3). The resutls were in accordance with the QTL mapping (Fig. 4 and Table 3). In this study, the map was consructed based on a F2 popualtion, which might lead to overestimated phentoypic effects of QTL due to large number of heterozygous loci in the F2 individuals and transgressive segreation . Further validation of the effects of this major QTL with advanced RIL population may be necessory.
Qin et al.  and Khera et al.  also reported major QTLs associated to TSWV resistance on chromosome A01. Due to large number of different polymorphic markers between the maps, we noticed that one (ARS721) of the two flanking markers (Seq12F7 and ARS721) of the QTL identified by Qin et al.  and two markers (Ah21 and Ah126) of the four flanking markers (Ah21, Ah126, Seq13A10 and GNB842) of QTLs identified by Khera et al.  were also mapped to A01 linkage group in our study. Based on the physical locations of these common SSR markers (Table 3), we could tell that our QTL location flanked by AHGS4584 and GM 672 were between the QTLs identified by Qin et al.  and Khera et al. , not quite the same as either of them. It could be either that the major QTL related to TSWV resistance in Florida-EPTM ‘113’ was different from that in the NC94022 or most likely that these reported QTLs earlier were the same but with shifted location due to environment effects . There was another possible TSWV resistance QTL located on chromosome A09 . Five markers showed polymorphism between two parental lines on chromosome A09 were further used to genotype 163 F2 plants in our experiment, however, these markers were not linked together based on this mapping population. Our association analysis using GBS data indicated one SNP on chromosome A09 significantly associated with TSWV resistance, which might be a minor QTL and will be investigated in the future.
Spotted wilt resistance QTLs
Spotted wilt dominant resistance genes have been found in tomato (Solanum lycopersicum) and pepper (Capsicum annuum), named Sw-5 and Tsw, respectively with the hypersensitive response (HR) resistant mechanism [56, 57]. In this study, the spotted wilt resistance in Florida-EPTM ‘113’ does not appear to be a HR mechanism, although one major QTL was identified. Sw-5 has an ortholog gene in peanut, Ahsw with 37 % amino acid identity to Sw-5 . The gene was characterized as peanut oxalate oxidase  and two gene-specific SSR markers (Seq2F10 and TC7G10) were identified. However, none of them displayed polymorphism between Florida-EPTM ‘113’ and Georgia Valencia, thus were not used to validate the linkage within our segregating population.
Two mapping populations have been utilized for spotted wilt resistance QTL mapping previously. They were derived from the cross between Tifrunner and GT-C20 (referred as T-population) and the cross between SunOleic 97R and NC94022 (referred as S-population) . Two resistant QTLs (qTSWV1 and qTSWV2) were first reported on S populations . qTSWV1 was located on linkage group 15 (LGJ15) on T-population and qTSWV2 was on A01 chromosome on S-population. Further 15 and 9 QTLs were identified using different generations (F2 and F5) on T-population , respectively. With S-population, there were six spotted wilt related QTLs were identified . The major QTL located on LG A01 identified based on S population was at the similar location as the QTLs identified in the study here. Most likely the resistance segregating in both populations was derived from same genetic source, PI 576638, known as hirsuta botanical-type line introduced from the highlands of Mexico.
Peanut is not a native crop to the United States, so plant introduction (PI) played an important role in peanut cultivar development. Two PIs, PI 203396 and PI 576638, have provided extensive disease resistant sources [60, 61]. It was reported that several breeding lines containing PI 576638 origin had better TSWV resistance than the lines derived from PI 203396. The two PI accessions may contain different resistant genes and have different resistance mechanisms . At the similar QTL region on A01, QTL for early leaf spot resistance was also detected , indicating that wild Arachis species displayed a high level resistance to several diseases (early leaf spot, late leaf spot, and stem rot) .
Marker assisted selection (MAS)
In this study, the spotted wilt resistant QTL with 22.7 % PVE was identified and the two flanking makers, AHGS4584 and GM672 can be applied to conduct MAS. The map distance between two markers was 2.55 cM on linkage map with 14.4 megabase (Mb) distance on physical map. Although one major QTL was identified , the linkage map distance on entire A01 chromosome was only 34.0 cM constructed by using the S segregating population reported by Qin et al. previously  and the QTL region (LOD >3.0) was quite wide, covering approximately one-third of the A01 chromosome. In our current study, the whole linkage group is much larger (157.85 cM) and the identified QTL was located in a much narrow region (2.55 cM), which can be used for MAS with increased confidence. The physical distance was estimated by A genome (A. duranensis) sequencing, not directly from the genomes of cultivated peanut, however, the distance was still far. Crossovers could happen between traits (the resistant QTL) and markers, causing recombination, which may lead to certain false positive results and decreased selection accuracy. More closely linked markers should be obtained to increase the efficiency of MAS.
GBS is feasible for large genome species with a low cost . In this study, 12 samples were selected to conduct GBS and the results indicated there were high associations between spotted wilt resistant loci and SNPs on chromosome A01. The SNPs identified at the region should be further validated in order to develop more potential markers, which are closely linked to disease resistance. Besides the disease resistant trait, other traits, for examples, seed maturity, seed size, seed number per pod, seed coat color, and leaf spot resistance, also showed differences in this mapping population. GBS can help to develop SNP markers associated to the traits mentioned above and will be available for constructing a high-density linkage map, QTL analysis, and MAS.
Expression of spotted wilt disease in peanut is highly variable depending on the disease pressure which varies by location and season. Phenotypic distribution results in this study indicated that disease pressure at NFREC was higher than PSREU. By using immunostrip testing, Florida-EPTM ‘113’ proved to be a good variety to efficiently separate resistant and susceptible genotypes. Since visual rating cannot detect asymptomatic infection, many “resistant” looking plants were actually susceptible genotypes. By screening polymorphic makers between Florida-EPTM ‘113’ and Georgia Valencia, only makers located on A01 chromosome co-segregated with spotted wilt resistance. The GBS method also indicated the strongest association between SNP markers and disease resistant loci on A01 chromosome. Other minor SNPs are worthy of investigating further in the future. The A01 linkage map constructed using the Florida-EPTM ‘113’ derived population had good marker collinearity with the physical map. Two QTLs have been identified on chromosome A01. One QTL was PSREU-specific and another was NFREC-specific. The latter one had higher LOD and PVE values than the former one. Historically, the spotted wilt pressure was lower at PSREU, so QTL regions identified based on phenotypic data at NFREC are more reliable with two flanking markers, AHGS4584 and GM672. This major QTL is most likely contributed by PI 576638, a hirsuta botanical-type line, introduced from Mexico. Marker enrichment in the region needs to be conducted to perform fine mapping to refine the QTL region. Next generation sequencing technology could be utilized in the future.
Genotyping by sequencing
General linear model
Inclusive Composite Interval Mapping
Minor allele frequency
North Florida Research and Education Center
Polyacrylamide gel electrophoresis
Polymerase chain reaction
Polymorphism information content
Plant Science Research and Education Unit
Phenotypic variation explained
Quantitative trait locus
Recombinant inbred line
Single nucleotide polymorphism
Simple sequence repeat
Tomato spotted wilt virus
Savage GP, Keenan JI: The composition and nutritive value of groundnut kernels. In The Groundnut Crop. Springer; 1994:173–213. http://link.springer.com/chapter/10.1007%2F978-94-011-0733-4_6.
Kochert G, Halward T, Branch WD, Simpson CE. RFLP variability in peanut (Arachis hypogaea L.) cultivars and wild species. Theor Appl Genet. 1991;81:565–70.
Bertrand PF. Georgia plant disease loss estimates. Univ Ga Coop Ext Pub Pathol. 1997;1998(81):98–107.
Culbreath AK, Todd JW, Brown SL. Epidemiology and management of tomato spotted wilt in peanut. Annu Rev Phytopathol. 2003;41:53–75.
Pandey MK, Monyo E, Ozias-Akins P, Liang X, Guimarães P, Nigam SN, Upadhyaya HD, Janila P, Zhang X, Guo B. Advances in Arachis genomics for peanut improvement. Biotechnol Adv. 2012;30:639–51.
Bertioli DJ, Cannon SB, Froenicke L, Huang G, Farmer AD, Cannon EKS, Liu X, Gao D, Clevenger J, Dash S. The genome sequences of Arachis duranensis and Arachis ipaensis, the diploid ancestors of cultivated peanut. Nat Genet. 2016;48:438–46.
Ferguson ME, Burow MD, Schulze SR, Bramel PJ, Paterson AH, Kresovich S, Mitchell S. Microsatellite identification and characterization in peanut (A. hypogaea L.). Theor Appl Genet. 2004;108:1064–70.
He G, Meng R, Newman M, Gao G, Pittman RN, Prakash CS. Microsatellites as DNA markers in cultivated peanut (Arachis hypogaea L.). BMC Plant Biol. 2003;3:3.
Gautami B, Foncéka D, Pandey MK, Moretzsohn MC, Sujay V, Qin H, Hong Y, Faye I, Chen X, BhanuPrakash A. An international reference consensus genetic map with 897 marker loci based on 11 mapping populations for tetraploid groundnut (Arachis hypogaea L.). PLoS One. 2012;7:e41213.
Nagy ED, Guo Y, Tang S, Bowers JE, Okashah RA, Taylor CA, Zhang D, Khanal S, Heesacker AF, Khalilian N. A high-density genetic map of Arachis duranensis, a diploid ancestor of cultivated peanut. BMC Genomics. 2012;13:469.
Shirasawa K, Bertioli DJ, Varshney RK, Moretzsohn MC, Leal-Bertioli SCM, Thudi M, Pandey MK, Rami JF, Ka DF, Gowda MVC, Qin H, Guo B, Hong Y, Liang X, Hirakawa H, Tabata S, Isobe S. Integrated consensus map of cultivated peanut and wild relatives reveals structures of the a and b genomes of arachis and divergence of the legume genomes. DNA Res. 2013;20:173–84.
Khedikar YP, Gowda MVC, Sarvamangala C, Patgar KV, Upadhyaya HD, Varshney RK. 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. 2010;121:971–84.
Nagy ED, Chu Y, Guo Y, Khanal S, Tang S, Li Y, Dong WB, Timper P, Taylor C, Ozias-Akins P. Recombination is suppressed in an alien introgression in peanut harboring Rma, a dominant root-knot nematode resistance gene. Mol Breed. 2010;26:357–70.
Sujay V, Gowda MVC, Pandey MK, Bhat RS, Khedikar YP, Nadaf HL, Gautami B, Sarvamangala C, Lingaraju S, Radhakrishan T, Knapp SJ, Varshney RK. 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. 2012;30:773–88.
Sarvamangala C, Gowda MVC, Varshney RK. Identification of quantitative trait loci for protein content, oil content and oil quality for groundnut (Arachis hypogaea L.). F Crop Res. 2011;122:49–59.
Tillman BL, Gorbet DW. Peanut cultivar UFT113, United States Patent 8178752. 2012. https://www.google.com/patents/US8178752. Accessed 15 May 2012.
Mckinney JL. Influence Of Planting Date, Plant Population, And Cultivar On Management Of Spotted Wilt In Peanut (Arachis Hypogaea L.). 2013. http://cms.uflib.ufl.edu/. Accessed Dec 2013.
Culbreath AK, Gorbet DW, Martinez-Ochoa N, Holbrook CC, Todd JW, Isleib TG, Tillman B. High Levels of Field Resistance to Tomato spotted wilt virus in Peanut Breeding Lines Derived from hypogaea and hirsuta Botanical Varieties. Peanut Sci. 2005;32:20–4.
Culbreath AK, Srinivasan R. Epidemiology of spotted wilt disease of peanut caused by Tomato spotted wilt virus in the southeastern US. Virus Res. 2011;159:101–9.
Qin H, Feng S, Chen C, Guo Y, Knapp S, Culbreath A, He G, Wang ML, Zhang X, Holbrook CC, Ozias-Akins P, Guo B. An integrated genetic linkage map of cultivated peanut (Arachis hypogaea L.) constructed from two RIL populations. Theor Appl Genet. 2012;124:653–64.
Khera P, Pandey MK, Wang H, Feng S, Qiao L, Culbreath AK, Kale S, Wang J, Holbrook CC, Zhuang W. Mapping Quantitative Trait Loci of Resistance to Tomato Spotted Wilt Virus and Leaf Spots in a Recombinant Inbred Line Population of Peanut (Arachis hypogaea L.) from SunOleic 97R and NC94022. PLoS One. 2016;11:e0158452.
Branch WD, Branch WD. Registration of “Georgia Valencia”peanut. Crop Sci. 2001;41:2002–a.
Dellaporta SL, Wood J, Hicks JB. A plant DNA minipreparation: version II. Plant Mol Biol Report. 1983;1:19–21.
Hopkins MS, Casa AM, Wang T, Mitchell SE, Dean RE, Kochert GD, Kresovich S. Discovery and characterization of polymorphic simple sequence repeats (SSRs) in peanut. Crop Sci. 1999;39:1243–7.
Liu Z, Feng S, Pandey MK, Chen X, Culbreath AK, Varshney RK, Guo B. Identification of expressed resistance gene analogs from peanut (Arachis hypogaea L.) expressed sequence tags. J Integr Plant Biol. 2013;55:453–61.
Gimenes MA, Hoshino AA, Barbosa AVG, Palmieri DA, Lopes CR. Characterization and transferability of microsatellite markers of the cultivated peanut (Arachis hypogaea). BMC Plant Biol. 2007;7:9.
Leal-Bertioli SCM, José ACVF, Alves-Freitas DMT, Moretzsohn MC, Guimarães PM, Nielen S, Vidigal BS, Pereira RW, Pike J, Fávero AP, Parniske M, Varshney RK, Bertioli DJ. Identification of candidate genome regions controlling disease resistance in Arachis. BMC Plant Biol. 2009;9:112.
Moretzsohn MC, Hopkins MS, Mitchell SE, Kresovich S, Valls JFM, Ferreira ME. Genetic diversity of peanut (Arachis hypogaea L.) and its wild relatives based on the analysis of hypervariable regions of the genome. BMC Plant Biol. 2004;4:11.
Moretzsohn MC, Barbosa AVG, Alves-Freitas DMT, Teixeira C, Leal-Bertioli SCM, Guimarães PM, Pereira RW, Lopes CR, Cavallari MM, Valls JFM, Bertioli DJ, Gimenes MA. A linkage map for the B-genome of Arachis (Fabaceae) and its synteny to the A-genome. BMC Plant Biol. 2009;9:40.
He G, Meng R, Gao H, Guo B, Gao G, Newman M, Pittman RN, Prakash CS. Simple sequence repeat markers for botanical varieties of cultivated peanut (Arachis hypogaea L.). Euphytica. 2005;142:131–6.
Liang X, Chen X, Hong Y, Liu H, Zhou G, Li S, Guo B. Utility of EST-derived SSR in cultivated peanut (Arachis hypogaea L.) and Arachis wild species. BMC Plant Biol. 2009;9:35.
Bertioli DJ, Moretzsohn MC, Madsen LH, Sandal N, Leal-Bertioli SCM, Guimarães PM, Hougaard BK, Fredslund J, Schauser L, Nielsen AM. An analysis of synteny of Arachis with Lotus and Medicago sheds new light on the structure, stability and evolution of legume genomes. BMC Genomics. 2009;10:45.
Cuc LM, Mace ES, Crouch JH, Quang VD, Long TD, Varshney RK. Isolation and characterization of novel microsatellite markers and their application for diversity assessment in cultivated groundnut (Arachis hypogaea). BMC Plant Biol. 2008;8:55.
Peng Z, Gallo M, Tillman BL, Rowland D, Wang J. Molecular marker development from transcript sequences and germplasm evaluation for cultivated peanut (Arachis hypogaea L.). Mol Genet Genomics. 2016;291(1):363–81.
Tang Wang C, Dao Yang X, Xu Chen D, Lin Yu S, Zhen Liu G, Yi Tang Y, Zhi Xu J. Isolation of simple sequence repeats from groundnut. Electron J Biotechnol. 2007;10:473–80.
Macedo SE, Moretzsohn MC, Leal-Bertioli SCM, Alves DMT, Gouvea EG, Azevedo VCR, Bertioli DJ. Development and characterization of highly polymorphic long TC repeat microsatellite markers for genetic analysis of peanut. BMC Res Notes. 2012;5:86.
Moretzsohn MC, Leoi L, Proite K, Guimarães PM, Leal-Bertioli SCM, Gimenes MA, Martins WS, Valls JFM, Grattapaglia D, Bertioli DJ. A microsatellite-based, gene-rich linkage map for the AA genome of Arachis (Fabaceae). Theor Appl Genet. 2005;111:1060–71.
Wang H, Penmetsa RV, Yuan M, Gong L, Zhao Y, Guo B, Farmer AD, Rosen BD, Gao J, Isobe S, Bertioli DJ, Varshney RK, Cook DR, He G. Development and characterization of BAC-end sequence derived SSRs, and their incorporation into a new higher density genetic map for cultivated peanut (Arachis hypogaea L.). BMC Plant Biol. 2012;12:10.
Koilkonda P, Sato S, Tabata S, Shirasawa K, Hirakawa H, Sakai H, Sasamoto S, Watanabe A, Wada T, Kishida Y. Large-scale development of expressed sequence tag-derived simple sequence repeat markers and diversity analysis in Arachis spp. Mol Breed. 2012;30:125–38.
Fountain J, Qin H, Chen C, Dang P, Wang ML, Guo B. A Note on Development of a Low-cost and High-throughput SSR-based Genotyping Method in Peanut (Arachis hypogaea L.). Peanut Sci. 2011;38:122–7.
Elshire RJ, Glaubitz JC, Sun Q, Poland JA, Kawamoto K, Buckler ES, Mitchell SE. A robust, simple genotyping-by-sequencing (GBS) approach for high diversity species. PLoS One. 2011;6:e19379.
Glaubitz JC, Casstevens TM, Lu F, Harriman J, Elshire RJ, Sun Q, Buckler ES. TASSEL-GBS: A high capacity genotyping by sequencing analysis pipeline. PLoS One. 2014;9(2):e90346.
Langmead B, Salzberg SL. Fast gapped-read alignment with Bowtie 2. Nat Methods. 2012;9:357–9.
Bradbury PJ, Zhang Z, Kroon DE, Casstevens TM, Ramdoss Y, Buckler ES. TASSEL: software for association mapping of complex traits in diverse samples. Bioinformatics. 2007;23:2633–5.
Yu J, Pressoir G, Briggs WH, Bi IV, Yamasaki M, Doebley JF, McMullen MD, Gaut BS, Nielsen DM, Holland JB. A unified mixed-model method for association mapping that accounts for multiple levels of relatedness. Nat Genet. 2006;38:203–8.
Wang J, Li H, Zhang L, Meng L: Users’ manual of QTL IciMapping version 3.2. Quant Genet Group, Inst Crop Sci Chinese Acad Agric Sci (CAAS), Beijing 2012, 100081. http://www.sciencedirect.com/science/article/pii/S2214514115000161.
Van Ooijen JW. JoinMap 4. Wageningen: Softw Calc Genet Link maps Exp Popul Kyazma BV; 2006.
Kosambi DD. The estimation of map distances from recombination values. Ann Eugen. 1943;12:172–5.
Von Bargen S, Salchert K, Paape M, Piechulla B, Kellmann JW. Interactions between the tomato spotted wilt virus movement protein and plant proteins showing homologies to myosin, kinesin and DnaJ-like chaperones. Plant Physiol Biochem. 2001;39:1083–93.
Li H, Ye G, Wang J. A modified algorithm for the improvement of composite interval mapping. Genetics. 2007;175:361–74.
Mandal B, Pappu HR, Culbreath a K, Pathology P: Factors Affecting Mechanical Transmission of Tomato spotted wilt virus to Peanut (Arachis hypogaea). 2001(December):1259–1263.
Singh BD, Singh AK. Marker-Assisted Plant Breeding: Principles and Practices. 2015.
Weinig C, Schmitt J. Environmental effects on the expression of quantitative trait loci and implications for phenotypic evolution. Bioscience. 2004;54:627–35.
Collard BCY, Jahufer MZZ, Brouwer JB, Pang ECK. An introduction to markers, quantitative trait loci (QTL) mapping and marker-assisted selection for crop improvement: The basic concepts. Euphytica. 2005;142:169–96.
Wang H, Pandey M, Qiao L, Qin H, Culbreath A, He G, Varshney R, Scully B. Genetic mapping of quantitative trait loci analysis for disease resistance Using F2 and F5 Generation-based Genetic Maps Derived from “Tifrunner” × “GT-C20” in Peanut. Plant Genome. 2013;6:1771–4.
Stevens MR, Scott SJ, Gergerich RC. Inheritance of a gene for resistance to tomato spotted wilt virus (TSWV) from Lycopersicon peruvianum Mill. Euphytica. 1991;59:9–17.
Boiteux LS, De Avila AC. Inheritance of a resistance specific to tomato spotted wilt tospovirus in Capsicum chinense “PI 159236.”. Euphytica. 1994;75:139–42.
Chen X, Culbreath A, Brenneman T, Holbrook Jr C, Guo B: Identification and cloning of TSWV resistance gene (s) in cultivated peanuts and development of markers for breeding selection. In American Phytopathological Society Abstracts; 2008. https://www.ars.usda.gov/research/publications/publication/?seqNo115=226033.
Chen X, Wang ML, Holbrook C, Culbreath A, Liang X, Brenneman T, Guo B. Identification and characterization of a multigene family encoding germin-like proteins in cultivated peanut (Arachis hypogaea L.). Plant Mol Biol Report. 2011;29:389–403.
Barrientos-Priego L, Isleib TG, Pattee HE. Variation in Oil Content Among Mexican and Peruvian hirsuta Peanut Landraces and Virginia-Type hypogaea Lines 1. Peanut Sci. 2002;29:72–7.
Isleib TG, Holbrook CC, Gorbet DW. Use of plant introductions in peanut cultivar development. Peanut Sci. 2001;28:96–113.
Holbrook CC, Stalker HT. Peanut breeding and genetic resources. Plant Breed Rev. 2003;22:297–356.
The research presented in this article was sponsored by Plant Molecular Breeding initiative (PMBI) at University of Florida and Southeastern Peanut Research Initiative (SPRI). The field technical support was provided by the peanut breeding group at Marianna and Citra, FL. We are grateful to Spurthi Nayak, Andrea Villa, Liping Wang, Lubin Tan, and Dev Paudel at University of Florida for technical support. We highly appreciate the support from Dr. Baozhou Guo at USDA ARS in providing the gel electrophoresis units and involving in project discussion.
Availability of data and materials
All relevant data are available within the manuscript or in the additional files.
JW and BT conceived and designed the experiments. BT developed the plant materials. YT conducted the experiment. YT and ZP analyzed the data. YT drafted the manuscript. JW and BT revised the manuscript. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
The number of plot, replication, check at different sites and phenotyping methods in F3, F4 and F5 populations. (DOCX 62 kb)
Literature reference of sources of SSR primers screened. (DOCX 49 kb)
TSWV infection results on F2:4 population and parental lines (checks) by visual rating in 2013. a) F2:4 population at PSREU; b) checks of F2:4 population at PSREU; c) F2:4 population at NFREC; d) checks of F2:4 population at NFREC. (DOCX 66 kb)
TSWV infection results on F2:5 population and parental lines (checks) by visual rating in 2014 at NFREC. a) F2:5 population; b) checks of F2:5 population. (DOCX 48 kb)
PAGE gel images showing genotypes of 12 samples by using polymorphic SSR markers located on A01, A09 and A10 chromosomes. a) GA110 markers located on A01 chromosome; b) AHGS1910 marker located on A01 chromosome; c) AHGS1647 markers located on A09 chromosome; d) AHGS1390 marker located on A10 chromosome. 1: Georgia Valencia, 2: Florida-EPTM ‘113’, 3 to 7: Susceptible plants, 8–12: resistant plants. (DOCX 763 kb)
About this article
Cite this article
Tseng, YC., Tillman, B.L., Peng, Z. et al. Identification of major QTLs underlying tomato spotted wilt virus resistance in peanut cultivar Florida-EPTM ‘113’. BMC Genet 17, 128 (2016). https://doi.org/10.1186/s12863-016-0435-9