Published online 30 July 2007
Published in Crop Sci 47:1336-1343 (2007)
© 2007 Crop Science Society of America
677 S. Segoe Rd., Madison, WI 53711 USA
CROP BREEDING & GENETICS
SSR Mapping and Confirmation of Soybean QTL from PI 437654 Conditioning Resistance to Reniform Nematode
Bo-Keun Haa,*,
Robert T. Robbinsb,
Feng Hanc,
Richard S. Husseyd,
John F. Soperc and
H. Roger Boermaa
a Dep. of Crop and Soil Science, Univ. of Georgia, Center for Applied Genetic Technologies, 111 Riverbend Rd., Athens, GA 30602
b Dep. of Plant Pathology, Univ. of Arkansas, Fayetteville, AR 72701
c Pioneer Hi-Bred International, Inc., Johnston, IA 50131
d Dep. of Plant Pathology, Univ. of Georgia, Athens, GA 30602
* Corresponding author (ha{at}uga.edu).
 |
ABSTRACT
|
|---|
Recently, reniform nematode (Rotylenchulus reniformis Linford and Oliveira) has been recognized to be a major problem on soybean [Glycine max (L.) Merr.] in the southern USA. Planting reniform nematode–resistant soybean cultivars is the most cost effective control method. With restriction fragment length polymorphism (RFLP) markers, Pioneer Hi-Bred International, Inc. had previously identified quantitative trait loci (QTL) on linkage group B1 (LG-B1) and another on LG-L associated with reniform nematode resistance in a cross between BSR101 and PI 437654. In this study, we refined the QTL locations and identified additional QTL conditioning resistance to reniform nematode in a population of 228 recombinant inbred lines (RILs) from a cross of BSR101 x PI 437654. A QTL (R2 = 21%) conditioning reniform reproductive index (RI) was found on LG-L and was flanked by Sat_184 and Satt513. Two other QTL were identified, one on LG-B1 (R2 = 16%) and the other on LG-G (R2 = 8%). These two QTL acted in an epistatic manner with lines homozygous for PI 437654 alleles at both QTL providing the lowest reniform RI. The 31 RILs from BSR101 x PI 437654 that were homozygous for the PI 437654 alleles at Satt513, Satt359, and Sat_168 averaged 1077 reniform eggs and juveniles (RI = 0.63), while the 33 RILs homozygous for the BSR101 alleles averaged 10606 (RI = 6.24). By screening the RILs population of Prichard x Anand we confirmed that the QTL on LG-G and-B1 conditioned reniform RI.
Abbreviations: LG, linkage group MAS, marker-assisted selection NIL, near isogenic line PCR, polymerase chain reaction Pf, final population Pi, initial inoculum level QTL, quantitative trait loci RFLP, restriction fragment length polymorphism RI, reproductive index RIL, recombinant inbred line SCN, soybean cyst nematode SSR, simple sequence repeat
 |
INTRODUCTION
|
|---|
THE RENIFORM nematode (Rotylenchulus reniformis Linford and Oliveira), a semi-endoparasite of roots (Linford and Oliveira, 1940; Robinson et al., 1997), is recognized to be a major problem on soybean [Glycine max (L.) Merr.] in the southern USA (Robbins et al., 1994a). The reniform nematode induces a feeding site, called a syncytium, and causes root decay, unthrifty growth, empty pods, and significant yield reduction in susceptible soybean cultivars (Birchfield et al., 1971; Williams and Birchfield, 1974). Parasitized plants exhibit varying degrees of stunting and chlorosis and initial field infestations may exhibit patchy areas of poor growth. Moderate to heavy soybean yield loss was observed from damage caused by reniform nematode in both Louisiana (Birchfield et al., 1971) and South Carolina (Fassuliotis and Rau, 1967). Rebois and Johnson (1973) reported that an initial population of 6579 reniform juveniles and pre-adult females per liter of soil consistently reduced yields on both moderately resistant and susceptible cultivars by an average of 33%, and reduced the content of P and K in harvested seeds.
Since nematicides and/or rotation to a nonhost crop are usually not economical for soybean producers, genetic resistance is the most desirable control method (Harville et al., 1985). Rebois et al. (1968) found that the soybean cultivars Pickett and Dyer were resistant to the reniform nematode. Also Rebois et al. (1970) reported that resistance to reniform nematode could be expected in genotypes with resistance to the soybean cyst nematode (SCN) (Heterodera glycines Ichinohe) due to the similarities in the feeding sites. However, Caviness and Riggs (1976) reported that Mack was resistant to SCN and susceptible to reniform nematode. Recent studies have shown that SCN-resistant soybean cultivars that derive resistance from Peking (e.g., Forrest and Centennial), PI 437654 (e.g., Hartwig), or PI 90763 (e.g., Cordell) are potentially resistant to R. reniformis, whereas those that derive SCN resistance from PI 88788 are not resistant to R. reniformis (Robbins et al., 1994a, 1994b; Davis et al., 1996; Robbins and Rakes, 1996). The phenotypic relationship suggests that common genes and/or linked genes are responsible for resistance to both SCN and reniform nematode in soybean.
The genetics of resistance to R. reniformis has been studied in soybean (Williams et al., 1981; Harville et al., 1985). Williams et al. (1981) reported that in a cross of Forrest (resistant) and Ransom (susceptible), reniform nematode resistance was recessive and controlled by alleles at one locus. In other soybean populations, resistance in soybean was controlled by two loci with unequal effects (Harville et al., 1985). The proposed genotypes of the four parental cultivars were (i) susceptible soybean cultivar Davis Rn1 Rn1 Rn2 Rn2, (ii) moderately susceptible soybean cultivar Bragg Rn1 Rn1 rn2 rn2, (iii) moderately resistant soybean cultivar Dare rn1 rn1 Rn2 Rn2, and (iv) resistant soybean cultivar Pickett 71 rn1 rn1 rn2 rn2. In an unpublished study (2000), Pioneer Hi-Bred International, Inc. identified two quantitative trait loci (QTL) for reniform nematode resistance in a cross between BSR101 and PI 437654 with restriction fragment length polymorphism (RFLP) markers covering whole linkage groups. This study basically employed the same 298 recombinant inbred lines (RILs) and 355 RFLP markers reported by Webb et al. (1995). One QTL on LG-L explained 20% of the total variation in reniform reproductive index (RI), and another QTL on LG-B1 accounted for 11% of variation.
Many studies using DNA markers have concentrated on QTL that control parasitic nematodes in soybean. Southern root-knot nematode [M. incognita (Kofoid and White) Chitwood] (Tamulonis et al., 1997c; Li et al., 2001), peanut root-knot nematode [M. arenaria (Neal) Chitwood] (Tamulonis et al., 1997b), and javanese root-knot nematode [M. javanica (Treub) Chitwood] resistance QTL (Tamulonis et al., 1997a) have been identified. In addition, DNA markers have been used to dissect quantitatively inherited resistance to the SCN (Concibido et al., 1994, 2004; Webb et al., 1995; Cregan et al., 1999). Webb et al. (1995) mapped three major QTL to LG-A2, -G, and -M, respectively, that confer resistance to SCN race 3 in a cross between PI 437654, resistant to SCN and reniform nematode (Robbins and Rakes, 1996), and BSR101. Cregan et al. (1999) reported that marker-assisted selection (MAS) using Satt309, located 1 to 2 cM from rhg1 on LG-G, should be effective to select for SCN resistance at rhg1 derived from Peking, PI 90763, and PI 437654 in crosses with most SCN-susceptible soybean genotypes. In a direct comparison, genotypic selection with Satt309 was 99% accurate in predicting lines that were susceptible in subsequent greenhouse assays (Mudge et al., 1997; Cregan et al., 1999).
The objectives of this study were to use simple sequence repeat (SSR) markers to: (i) refine the location of existing, but unpublished QTL and identify additional QTL conferring resistance to reniform nematode in the population of BSR101 x PI 437654 previously used in reniform and SCN mapping study, (ii) determine the magnitude of resistance of each QTL, and (iii) evaluate the potential interaction of identified QTL.
 |
MATERIALS AND METHODS
|
|---|
Mapping Population
A population of 228 F6–derived RILs from a cross of BSR101 x PI 437654 was used in this study. This population was maintained and increased in Iowa by Pioneer Hi-Bred International, Inc. PI 437654 is a plant introduction that originated from China but it was received in the USA from the Russian Federation in 1979 and placed in the USDA Soybean Germplasm Collection (Nelson et al., 1988). It is Maturity Group III and is resistant to reniform nematode (Davis et al., 1996; Robbins and Rakes, 1996). BSR101 is a Maturity Group I cultivar developed at Iowa State University and is susceptible to reniform nematode (Tachibana et al., 1987).
Reniform Nematode Screening
Reniform nematode resistance screening was conducted in a greenhouse at the University of Arkansas. Seeds of each genotype were germinated in vermiculite and as the cotyledons where starting to open the seedlings were transplanted into 10-cm-diam. clay pots (one seedling per pot) containing 500 cm3 of pasteurized fine sandy loam soil (ca. 91% sand, 5% silt, 4% clay, <1% organic matter). One day after transplanting each pot was inoculated with 1700 vermiform reniform nematodes (i.e., initial population; Pi). The reniform population used was from a culture originating in Jefferson County, Arkansas, and maintained on Braxton soybean. The experimental design was a randomized complete block with five replications. All pots were watered at least twice daily and other times, if needed, and fertilized each week with 12 mg N, 6 mg P, and 10 mg K.
Thirteen weeks after inoculation, the number of reniform nematode eggs and vermiform nematodes contained in egg masses on the roots and the number of vermiform nematodes in the soil of each pot were determined. The eggs and vermiform nematodes in the egg masses on roots were collected with 0.525% NaOCl (Hussey and Barker, 1973) and counted. To calculate the final reniform nematode soil population, a 100-cm3 aliquot of well-mixed soil from each pot was suspended in water and poured through nested 850- and 38-µm-pore sieves to remove plant debris and collect the nematodes. Nematodes caught on the 38-µm-pore sieve were separated from soil by sucrose centrifugal-flotation (Jenkins, 1964), counted, and multiplied by 5 to give the total soil number per pot. The total number of reniform nematode eggs and vermiform nematodes per pot was calculated by adding the number from the soil to the number from the roots. A reproductive index (RI), defined as the number of eggs + vermiform nematodes at test termination (Pf)/initial infestation level (Pi), was calculated for each pot. Because of an association between the mean and variance of a genotype, the RI of each pot was transformed to log ratio [log10 (RI + 1)] values for statistical analysis and presented as antilogs.
SSR Marker Data Collection
Soybean DNA was extracted from seeds of the 228 RILs according to modified procedures of Kang et al. (1998), quantified using a PerkinElmer Lambda 2 UV/VIS Spectrometer (PerkinElmer Life and Analytical Sciences, Shelton, CT), and diluted to 20 ng µL–1 for polymerase chain reaction (PCR). Ten seeds from each line were ground with coffee grinder (Braun KSM2, Boston, MA), and 0.1 g of each homogenate was transferred to a new 1.5-mL tube containing 500 µL extraction buffer (200 mM Tris pH 8, 200 mM NaCl, 25 mM ethylenediamine tetraacetic acid [EDTA], 0.5% sodium dodecyl sulfate) and 10 µL proteinase K (20 mg m L–1). The samples were incubated in a water bath at 50°C for 1 h and then 500 µL CTAB solution (2% cetyltrimethylammonium bromide, 100 mM Tris-HCl pH 8, 20 mM EDTA pH 8, 1.4 M NaCl) was added. After shaking, the samples were spun at 10000 g (Beckman Microfuge E, Beckman Instruments, Carlsbad, CA) for 10 min. The supernatant was transferred to a new 1.5-mL tube, and then one volume chloroform/isoamyl alcohol (24:1, v/v) was added. After shaking for 1 min at room temperature, the samples were spun at 10000 g for 10 min. The supernatant was transferred to a new 1.5-mL tube, and then 80% volume of isopropyl alcohol was added to precipitate DNA. The supernatant was decanted and the pellets were washed with 70% (v/v) ethanol. The DNA pellets were then dried and dissolved in 100 µL of Tris-EDTA buffer.
On the basis of relationship with SCN and reniform nematode resistance, 96 SSR markers were chosen. These markers were near previously reported QTL for SCN and reniform resistance and were located on five linkage groups (A2, B1, G, L, and M) (Concibido et al., 2004). The primer sequences for each SSR were obtained from SoyBase, a USDA-sponsored genome database (http://soybase.org/). Fluorescent-labeled forward primers and nonlabeled reverse primers were obtained from PE-ABI (Foster City, CA). Polymerase chain reactions were prepared based on the protocol of Diwan and Cregan (1997) with slight modifications. The 10-µL reaction mix contained 2 µL of 40 ng template DNA, 1.0x PCR buffer, 2.5 mM MgCl2, 100 µM of each dNTP, 0.2 µM each of forward and reverse primers, and 0.5 unit of Taq DNA polymerase. The reactions were performed in a MJ Research PTC-225 thermal cycler (MJ Research, Inc., Reno, NV). The PCR amplicons were analyzed on an ABI-Prism 377 DNA sequencer (PE-ABI) with a 4.8% acrylamide to bisacrylamide (19:1) gel at 750 V for 2 h. The marker fragments were analyzed with GeneScan software v.3.0 and scored with Genotyper software v.2.1 (PE-ABI).
Linkage Map and QTL Analysis
Initially, 134 randomly selected RILs were evaluated for the association of SSR marker alleles and reniform nematode RI. To further increase precision, 94 additional RILs were included. A broad sense heritability estimate based on variance component was calculated on entry-mean basis with the following formula: h2=
RILs2/[
RILs2+(
e2/r)] (Fehr, 1987), where
RILs2 is the genotypic variance component for reniform RI per plant among RILs,
e2 is the error variance, and r is the number of replications.
Single factor analysis of variance (SF-ANOVA) was used to determine the significance (P = 0.05) among SSR genotypic class means by use of General Linear Model (GLM) (SAS Institute, 2001). A multiple regression model with a FORWARD option was used for identifying the independent markers linked to the QTL within and among linkage groups at the 5% significance level.
Genetic linkage maps were constructed in MAP MANAGER QTX Version b16 with the Kosambi map function (Manly et al., 2001). Quantitative trait loci analyses of RI were conducted using interval mapping implemented in the program MAP MANAGER QTX (Manly et al., 2001). Genome-wide threshold values for declaring the presence of QTL were determined by 1000 permutations for reniform RI, as implemented by MAP MANAGER QTX. The additive effects (
) for individual QTL were estimated by MAP MANAGER QTX.
Multilocus QTL analyses were performed with three genetic marker loci in mixed linear models (Littel et al., 1996; Tang et al., 2006). Statistical analyses were performed by SAS PROC MIXED, where genotype (G) effects were fixed and RIL nested within genotype (RIL/G), and replications (R) were random. Variance components were estimated by the REML method. Type III sums of squares and F statistics were estimated for genotype effects, and least square means were estimated for genotypes. Additive (
), additive x additive (
x
), and additive x additive x additive (
x
x
) effects were estimated by ESTIMATE statements.
QTL Confirmation
To evaluate the effect of identified QTL in a different genetic background, F7–derived RILs from the cross of Prichard x Anand were evaluated for reniform RI. Anand originated from the cross Holladay x Hartwig and was reported as resistant to reniform nematode (Anand et al., 2001; Robbins et al., 1999). Hartwig was developed from the cross of Forrest (3) x PI 437654 (Anand, 1992). Prichard originated from the cross of Co82–622 x Howard and was susceptible to reniform nematode (Boerma et al., 2001). Both Howard and Co82–622 include Forrest in their pedigree. A total of 353 RILs were genotyped using the SSR markers Satt513 on LG-L, Satt359 on LG-B1, and Sat_168 on LG-G. Our goal was to identify eight RILs for each of the eight possible homozygous marker classes. However, one genotypic class contained only two RILs and another class had only one RIL out of 353 RILs. Therefore, a total of 51 RILs were evaluated for reniform nematode resistance in a greenhouse as described earlier. The screening was conducted in a randomized complete block experimental design with five replications. Each plant was inoculated with 1644 vermiform reniform nematode (Pi). Final nematode population (Pf) for each pot was determined only from the soil.
 |
RESULTS AND DISCUSSION
|
|---|
Phenotypic Variation and Heritability
The mean reniform RI and total number of eggs + vermiform nematode (in parentheses) was 0.20 (336) for the resistant parent PI 437654 and 3.51 (5971) for the susceptible parent BSR101. The mean RI for the 228 RILs ranged from 0.14 to 15.05 per plant (Fig. 1
). Log transformed reniform RI data for the 228 RILs were used to determine heritability. The heritability on entry-mean basis (mean of five replications) was 0.78, suggesting that there was a major genetic component conditioning reniform nematode resistance in this population.

View larger version (14K):
[in this window]
[in a new window]
|
Figure 1. Distribution of reniform reproductive index (RI) for 228 recombinant inbred lines from BSR101 x PI 437654.
|
|
QTL Conditioning Reniform Nematode Resistance
To map reniform nematode resistance, we examined five linkage groups that had been previously associated with resistance to reniform and SCN. Ninety-six previously mapped SSR markers on LG-A2, -B1, -G, -L, and -M (Song et al., 2004) were screened for polymorphism between the parental genotypes. Sixty-one of these 96 markers (64%) were polymorphic. Among these, 53 SSR markers were used to genotype 134 of the 228 RILs. Based on the association of the reniform RI phenotypes and the SSR marker alleles for those 53 markers by SF-ANOVA, we putatively identified three QTL associated with reniform nematode RI on LG-B1, -G, and -L (data not shown).
To further increase the precision for QTL mapping, 94 additional RILs were genotyped with 28 SSR markers on LG-B1, -L, and -G. Data from these 94 RILs were combined with the 134 original RILs. When the marker and phenotypic data from the 228 RILs were analyzed by interval mapping with Map Manager QTX, a QTL mapped to a 25.3-cM interval between Sat_184 and Satt513 on LG-L and explained 21% of the variation in reniform RI (Fig. 2
). The QTL peak was 8 cM from Satt513 (LOD = 11.8). Another QTL mapped to a 14.5-cM interval between Satt359 and Satt484 on LG-B1 and explained 16% of the phenotypic variation (LOD = 8.8). The QTL peak was 5 cM from Satt359. A third QTL, accounting for 8% of the phenotypic variation in reniform RI (LOD = 4.1), was linked to marker Sat_168 on LG-G (Fig. 2).

View larger version (21K):
[in this window]
[in a new window]
|
Figure 2. Maximum likelihood plots indicating genomic regions of quantitative trait loci associated with reniform reproductive index (RI) in 228 recombinant inbred lines from BSR101 x PI 437654. Lines parallel to the linkage groups indicate the genome-wide = 0.001 significance threshold (LOD = 3.8).
|
|
The three single-linkage group multiple regression analyses that included the SSR markers found significantly associated with reniform RI based on SF-ANOVA on LG-B1 (seven markers), -G (six markers), and -L (four markers) retained Satt359 on LG-B1, Sat_168 on LG-G, and Satt513 and Sat_184 on LG-L (Table 1). Satt513 on LG-L accounted for the most variation (18.1%) in reniform RI across linkage groups. Marker Sat_184 on LG-L explained an additional 3.2% of the variation in reniform RI. Given a single QTL peak was identified by interval mapping of LG-L, it is likely Satt513 and Sat_184 are flanking a single LG-L reniform RI QTL.
View this table:
[in this window]
[in a new window]
|
Table 1. Soybean simple sequence repeat (SSR) markers associated with reniform nematode reproductive index (RI) in 228 recombinant inbred lines from BSR101 x PI 437654.
|
|
Satt359 on LG-B1 accounted for the 17.3% and Sat_168 on LG-G accounted for 8.5% of phenotypic variation in reniform RI. RILs homozygous for PI 437654 alleles had significantly lower reniform RI than RILs homozygous for BSR101 alleles. When all four significant markers from LG-B1, -L, and -G were included in a multiple regression analysis, all four SSR markers were retained in the model, and collectively explained 36% of the phenotypic variation in reniform RI (Table 1).
QTL Effects on Reniform Nematode Resistance
Three QTL conferring reniform RI showed a significant (P < 0.05) additive main effect (Table 1). Also, significant (P < 0.01) additive x additive epistasis was detected between Satt359 and Sat_168 (Table 2). The LG-B1 QTL effect was greater than the LG-G QTL effect (Table 1 and 2). However, RILs homozygous for PI 437654 alleles at Satt359 and Sat_168 had a reniform RI of 0.79 which was significantly lower than the mean RI (RI = 4.5) for the other three genotypes (Table 2).
View this table:
[in this window]
[in a new window]
|
Table 2. Main effects and interaction of two reniform reproductive index (RI) quantitative trait loci in the recombinant inbred lines from BSR101 x PI 437654.
|
|
The mean reniform RI for RILs homozygous for PI 437654 alleles at Satt513, Satt359, and Sat_168 was 0.63 and was not significantly different from the resistant parent PI 437654 (RI = 0.20), whereas RILs homozygous for BSR101 alleles had a mean RI of 6.24, which was significantly greater than susceptible parent BSR 101 (RI = 3.51) (Table 3). The LG-L QTL had a larger effect on RI than the LG-B1 or LG-G QTL. The combined effect on RI of the two LG-B1 and LG-G QTL was very similar to the combined effect of the three LG-L, LG-B1, and LG-G QTL. Recombinant inbred lines homozygous for PI 437654 alleles at the QTL on LG-B1 and LG-G averaged a RI of 1.49 compared to a RI of 0.63 for the RILs that were homozygous for the PI 437654 alleles at all three QTL. These data suggest that having a resistance allele at only one of the three QTL would not produce a high level of resistance in a commercial soybean cultivar.
View this table:
[in this window]
[in a new window]
|
Table 3. The effect on reniform nematode reproductive index (RI) for quantitative trait loci at Satt513 (LG-L), Satt359 (LG-B1), and Sat_168 (LG-G) in 228 recombinant inbred lines from BSR101 x PI 437654.
|
|
The three individual markers, Satt359 (LG-B1), Sat_168 (LG-G), and Satt513 (LG-L) associated with QTL for reniform RI, the three two-way interaction, and one three-way interaction among marker loci were entered into multilocus QTL analyses performed by SAS PROC MIXED. The main effects of Satt359, Satt513, Sat_168, and the Satt359 x Sat_168 interaction were significant (P = 0.01) in the model and explained 60% of the phenotypic variation in reniform RI. With the heritability estimate of 78% for reniform RI in this population, the multilocus model explained 76% (60/78) of the genotypic variation.
Of the 228 RILs in this population, 206 RILs were homozygous and without missing data for the SSR markers Satt513 on LG-L, Satt359 on LG-B1, and Sat_168 on LG-G. For the eight possible homozygous classes of the three markers, we expected an equal distribution of lines in each class (approx. 26 lines per class). The actual number of lines in these marker classes was significantly different (
2 = 28.25, P < 0.01) from expectations (Table 3). One genotypic class that had the BSR101 (B) allele at Satt513 and PI 437654 (P) alleles at the Satt359 and Sat_168 included only seven lines. Also, the B/B/P and B/P/B classes for Satt513, Satt359, and Sat_168, respectively, contained fewer than the 26 expected lines. Apparently, the combination of alleles from BSR101 in the region of Satt513 on LG-L with alleles from PI 437654 in the region of Satt359 on LG-B1 and the combination of alleles from BSR101 in the region of Satt513 on LG-L with alleles from PI 437654 in the region of Sat_168 on LG-G were deleterious to survival or were unintentionally selected against during population development.
A distortion in marker classes also had been reported in another study that mapped SCN resistance using the same BSR101 x PI 437654 population that was also maintained and increased in Iowa by Pioneer Hi-Bred International, Inc. (Webb et al., 1995). A recent study with three near isogenic line (NIL) populations segregating at the SCN resistance locus rhg1 on LG-G reported that SCN-resistant NILs had significantly less field emergence than susceptible NILs (Kopisch-Obuch and Diers, 2006). Results of that study suggested that the segregation distortion was caused by pleiotropic effects of rhg1 or genetic linkage with unfavorable genes causing reduced field emergence.
Confirmation of the QTL in a Different Genetic Background
Before the application of MAS it is desirable to confirm QTL in the same population with a different sample of lines or in a different genetic background (Fasoula et al., 2004). A population of 353 RILs derived from the cross of Prichard x Anand was evaluated to confirm the previously described QTL for reniform RI.
Anand and its ancestors (Holladay, Hartwig, Forrest, and PI 437654) were screened with SSR markers flanking QTL on LG-L, -B1, and -G (Table 4). Anand inherited the allele at the reniform QTL on LG-B1 from Hartwig (via Forrest), the allele at the reniform QTL on LG-G from Hartwig (via PI 437654 or Forrest), and the allele at the QTL on LG-L from Holladay. Putatively, Anand appeared to possess resistance alleles at the LG-B1 and LG-G QTL, but not at the LG-L QTL.
View this table:
[in this window]
[in a new window]
|
Table 4. Comparison of band sizes of simple sequence repeat (SSR) markers linked to quantitative trait loci conditioning reniform nematode reproductive index (RI) for Anand, its ancestors, and Prichard.
|
|
In greenhouse screening, Anand averaged a reniform RI of 5.9 while Prichard averaged 19.7. Resistant checks, PI 437654, Forrest, and Hartwig, averaged an RI of 1.7, 6.4, and 4.7, respectively, and susceptible check, Braxton, averaged an RI of 38.5. These data indicate that Hartwig did not inherit all the reniform resistance alleles present in PI 437654 and that Prichard may possess partial reniform nematode resistance. A group of 51 RILs that were homozygous for Prichard or Anand alleles at Satt513 (LG-L), Satt359 (LG-B1), and Sat_168 (LG-G) were selected (Table 5). Eight RILs homozygous for the Prichard alleles at Satt513, Satt359, and Sat_168 averaged a RI of 26.8, while two RILs homozygous for the Anand alleles averaged 5.5. Only Sat_168, near the reniform QTL on LG-G, showed a significant (P < 0.01) additive main effect and appeared to impact reniform RI. Given the presence of the susceptible Holladay allele in Anand at LG-L QTL, the lack of difference between RILs homozygous for the Anand and Prichard alleles at this QTL was predicted based on the QTL mapping study. The lack of a difference in reniform RI for RILs homozygous for either Anand or Prichard alleles at the reniform QTL on LG-B1 and the large impact on RI of Sat_168 suggest that both Anand and Prichard possess a reniform resistance allele at this QTL. If correct, this would result in all the RILs in this population possessing a reniform resistance allele at the LG-B1 QTL and explain the lower reniform RI in Prichard (RI = 19.7) compared to Braxton (RI = 38.5). The fixation of a reniform resistance allele at the LG-B1 QTL would also explain the large reduction in reniform RI shown by RILs homozygous for the Anand allele compared to the Prichard allele at the LG-G QTL (Table 5). Prichard could have inherited an allele for reniform resistance from either of its parents, Co82-622 or Howard, since the parents of Prichard include Forrest in their pedigrees. While the lack of segregation of a resistance allele at the reniform QTL on LG-L precluded its confirmation in this population, our results serve to confirm the reniform resistance conditioned by the combination of reniform resistance alleles at QTL on LG-G and LG-B1.
View this table:
[in this window]
[in a new window]
|
Table 5. The effect on reniform nematode reproductive index (RI) of three quantitative trait loci linked to markers Satt513 (LG-L), Satt359 (LG-B1), and Sat_168 (LG-G) in the recombinant inbred lines from Prichard x Anand.
|
|
The confirmation population also showed segregation distortion in marker classes. One genotypic class that was homozygous for the Anand alleles at Satt513, Satt359, and Sat_168 contained only two lines out of 353 RILs evaluated (Table 5). Also, one class that had the Anand alleles at Satt513 and Satt359 and the Prichard allele at Sat_168 had only one line out of 353 RILs evaluated. The lines in this population were selected for late maturity and there are known maturity QTL in the same region of the genome as the three reniform QTL. Maturity gene E3, a major locus conferring photoperiod insensitivity in soybean, is located between Satt006 and Satt373 on LG-L (Tasma and Shoemaker, 2003). Pod maturity QTL are also located near Satt359 on LG-B1 and near Satt309 on LG-G, respectively (Lee et al., 1996; Brucker et al., 2005). Therefore due to genetic linkage of maturity and reniform RI QTL, it is likely the Prichard alleles would be more frequent in the lines from this population near markers Satt513 and Satt359.
This study was initiated based on the observation that there was a genetic linkage between SCN and reniform nematode resistance. Many studies reported that SCN-resistant soybean cultivars that derive resistance from Peking, PI 437654, or PI 90763 are also resistant to reniform nematode (Robbins et al., 1994a, 1994b; Davis et al., 1996; Robbins and Rakes, 1996). The SCN-resistant genotype PI 437654 has major SCN QTL on LG-G (rhg1) and LG-A2 (Rhg4). The SSR marker Satt309, located 1 to 2cM from rhg1 has been shown to be useful for MAS for SCN resistance (Cregan et al., 1999). Concibido et al. (2004) reported that QTL on LG-A1, -B1, -C1, -J, and -M conditioned SCN resistance in PI 437654. Interestingly, the QTL for SCN resistance linked to RFLP marker A006 on LG-B1 was reported by Vierling et al. (1996) to explain more than 90% of the phenotypic variation in a cross between Williams 82 and Hartwig. Hartwig was derived from backcrossing SCN resistance from PI 437654 into Forrest (Anand, 1992). Also, a SCN-resistant QTL linked to an RFLP marker (A567a) on LG-B1 was found in a cross from BSR101 x PI 437654, according to a patent claim by Webb (2003).
Harville et al. (1985) suggested that reniform nematode resistance was quantitatively controlled by two genes with unequal effects. By contrast, we found a QTL on LG-L and two epistatically interacting QTL on LG-B1 and LG-G. The reniform resistance QTL on LG-G with its resistance allele derived from PI 437654 was located near SSR marker Sat_168 which is closely linked to SCN resistance gene rhg1. Also, the reniform QTL on LG-B1 located between Satt359 and Satt484 is positioned between RFLP markers A006 and A567a, which were reported to be linked to SCN QTL (Vierling et al., 1996; Webb, 2003). Therefore, our results support linkage between certain SCN and reniform resistance QTL, provide confirmation of the QTL on LG-G, and evidence of an interaction of QTL for reniform RI on LG-G and LG-B1.
 |
NOTES
|
|---|
All rights reserved. No part of this periodical may be reproduced or transmitted in any form or by any means, electronic or mechanical, including photocopying, recording, or any information storage and retrieval system, without permission in writing from the publisher. Permission for printing and for reprinting the material contained herein has been obtained by the publisher.
Received for publication October 9, 2006.
 |
REFERENCES
|
|---|
- Anand, S.C. 1992. Registration of Hartwig soybean. Crop Sci. 32:1069–1070.[Free Full Text]
- Anand, S.C., T. Newman, and J. Fisher. 2001. Registration of Anand soybean. Crop Sci. 41:919–920.[Free Full Text]
- Birchfield, W., C. Williams, E.E. Hartwig, and L.R. Brister. 1971. Reniform nematode resistance in soybeans. Plant Dis. Rep. 55:1043–1045.
- Boerma, H.R., R.S. Hussey, D.V. Phillips, E.D. Wood, G.B. Rowan, S.L. Finnerty, and J.T. Griner. 2001. Registration of Prichard soybean. Crop Sci. 41:920–921.[Free Full Text]
- Brucker, E., T. Niblack, F.J. Kopisch-Obuch, and B.W. Diers. 2005. The effect of rhg1 on reproduction of Heterodera glycines in the field and greenhouse and associated effects on agronomic traits. Crop Sci. 45:1721–1727.[Abstract/Free Full Text]
- Caviness, C.E., and R.D. Riggs. 1976. Breeding for nematode resistance. p. 594–601. In L.D. Hill (ed.) World soybean research. Proc. of the World Soybean Research Conf., Champaign, IL. 3–8 Aug. 1975. The Interstate Printers & Publishers, Danville, IL.
- Concibido, V.C., R.L. Denny, S.R. Boutin, R. Hautea, J.H. Orf, and N.D. Young. 1994. DNA marker analysis of loci underlying resistance to soybean cyst nematode (Heterodera glycines Ichinohe). Crop Sci. 34:240–246.[Abstract/Free Full Text]
- Concibido, V.C., B.W. Diers, and P.R. Arelli. 2004. A decade of QTL mapping for cyst nematode resistance in soybean. Crop Sci. 44:1121–1131.[Abstract/Free Full Text]
- Cregan, P.B., J. Mudge, E.W. Fickus, D. Danesh, R. Denny, and N.D. Young. 1999. Two simple sequence repeat markers to select for soybean cyst nematode resistance conditioned by the rhg1 locus. Theor. Appl. Genet. 99:811–818.[CrossRef][Web of Science]
- Davis, E.L., S.R. Koenning, J.W. Burton, and K.R. Barker. 1996. Greenhouse evaluation of selected soybean germplasm for resistance to North Carolina populations of Heterodera glycines, Rotylenchulus reniformis, and Meloidogyne species. J. Nematol. 28:590–598.[Medline]
- Diwan, N., and P.B. Cregan. 1997. Automated sizing of fluorescent-labeled simple sequence repeat (SSR) markers to assay genetic variation in soybean. Theor. Appl. Genet. 95:723–733.[CrossRef][Web of Science]
- Fasoula, V.A., D.K. Harris, and H.R. Boerma. 2004. Validation and designation of quantitative trait loci for seed protein, seed oil, and seed weight from two soybean populations. Crop Sci. 44:1218–1225.[Abstract/Free Full Text]
- Fassuliotis, G., and G.J. Rau. 1967. The reniform nematode in South Carolina. Plant Dis. Rep. 54(7):557.
- Fehr, W.R. 1987. Heritability. p. 95–105. In Principles of cultivar development. Vol. I. Theory and technique. Macmillan, New York.
- Harville, B.G., A. Green, and W. Birchfield. 1985. Genetic resistance to reniform nematodes in soybeans. Plant Dis. 69:587–589.
- Hussey, R.S., and K.R. Barker. 1973. A comparison of methods of collecting inocula of Meloidogyne spp., including a new technique. Plant Dis. Rep. 57:1025–1028.
- Jenkins, W.R. 1964. A rapid centrifugal-flotation technique for separating nematodes from soil. Plant Dis. Rep. 48:692.
- Kang, H.W., Y.G. Cho, U.H. Yoon, and M.Y. Eun. 1998. A rapid DNA extraction method for RFLP and PCR analysis from a single dry seed. Plant Mol. Biol. Rep. 16:1–9.[CrossRef]
- Kopisch-Obuch, F.J., and B.W. Diers. 2006. Segregation at the SCN resistance locus rhg1 in soybean is distorted by an association between the resistance allele and reduced field emergence. Theor. Appl. Genet. 112:199–207.[CrossRef][Web of Science][Medline]
- Lee, S.H., M.A. Bailey, M.A.R. Mian, T.E. Carter, Jr., D.A. Ashley, R.S. Hussey, W.A. Parrott, and H.R. Boerma. 1996. Molecular markers associated with soybean plant height, lodging, and maturity across locations. Crop Sci. 36:728–735.[Abstract/Free Full Text]
- Li, Z., L. Jakkula, R.S. Hussey, J.P. Tamulonis, and H.R. Boerma. 2001. SSR mapping and confirmation of the QTL from PI 96354 conditioning soybean resistance to southern root-knot nematode. Theor. Appl. Genet. 103:1167–1173.[CrossRef][Web of Science]
- Linford, M.B., and J.M. Oliveira. 1940. Rotylenchulus reniformis, nov. gen., n. sp., a nematode parasite of roots. Proc. Helminthological Soc. Washington 7:35–42.
- Littel, R.C., G.A. Milliken, W.W. Stroup, and R.D. Wolfinger. 1996. SAS system for mixed models. SAS Inst., Cary, NC.
- Manly, K.F., R.H. Cudmore, Jr., and J.M. Meer. 2001. Map manager QTX, cross-platform software for genetic mapping. Mamm. Genome 12:930–932.[CrossRef][Web of Science][Medline]
- Mudge, J., P.B. Cregan, J.P. Kenworthy, W.J. Kenworthy, J.H. Orf, and N.D. Young. 1997. Two microsatellite markers that flank the major soybean cyst nematode resistance locus. Crop Sci. 37:1611–1615.[Abstract/Free Full Text]
- Nelson, R.L., P.J. Amdor, J.H. Orf, and J.F. Cavins. 1988. Evaluation of the USDA soybean germplasm collection: Maturity groups 000 to IV (PI 427.136 to PI 445.845). USDA-ARS Tech. Bull. 1726. USDA, Washington, DC.
- Rebois, R.V., J.M. Epps, and E.E. Hartwig. 1970. Correlation of resistance in soybeans to Heterodera glycines and Rotylenchulus reniformis. Phytopathology 60:695–700.[Web of Science]
- Rebois, R.V., and W.C. Johnson. 1973. Effect of Rotylenchulus reniformis on yield and nitrogen, potassium, phosphorus and amino acid content of seed of Glycine max. J. Nematol. 5:1–6.[Medline]
- Rebois, R.V., W.C. Johnson, and E.J. Cairns. 1968. Resistance in soybeans, Glycine max L. Merr., to the reniform nematode. Crop Sci. 8:394–395.[Abstract/Free Full Text]
- Robbins, R.T., and L. Rakes. 1996. Resistance to the reniform nematode in selected soybean cultivars and germplasm lines. J. Nematol. 28:612–615.[Medline]
- Robbins, R.T., L. Rakes, and C.R. Elkins. 1994a. Reniform nematode reproduction and soybean yield of four soybean cultivars in Arkansas. J. Nematol. 26:656–658.[Medline]
- Robbins, R.T., L. Rakes, and C.R. Elkins. 1994b. Reproduction of the reniform nematode on thirty soybean cultivars. J. Nematol. 26:659–664.[Medline]
- Robbins, R.T., L. Rakes, L.E. Jackson, and D.G. Dombek. 1999. Reniform nematode resistance in selected soybean cultivars. J. Nematol. 31(4S):667–677.
- Robinson, A.F., R.N. Inserra, E.P. Caswell-Chen, N. Vovlas, and A. Troccoli. 1997. Rotylenchus species: Identification, distribution, host ranges, and crop plant resistance. Nematropica 27:127–180.[Web of Science]
- SAS Institute. 2001. Version 8. SAS Inst., Cary, NC.
- Song, Q.J., L.F. Marek, R.C. Shoemaker, K.G. Lark, V.C. Concibido, X. Delannay, J.E. Specht, and P.B. Cregan. 2004. A new integrated genetic linkage map of soybean. Theor. Appl. Genet. 109:122–128.[CrossRef][Web of Science][Medline]
- Tachibana, H., B.K. Voss, and W.R. Fehr. 1987. Registration of BSR101 soybean. Crop Sci. 27:612.
- Tamulonis, J.P., B.M. Luzzi, R.S. Hussey, W.A. Parrott, and H.R. Boerma. 1997a. DNA markers associated with resistance to Javanese root-knot nematode in soybean. Crop Sci. 37:783–788.[Abstract/Free Full Text]
- Tamulonis, J.P., B.M. Luzzi, R.S. Hussey, W.A. Parrott, and H.R. Boerma. 1997b. DNA markers analysis of loci conferring resistance to peanut root-knot nematode in soybean. Theor. Appl. Genet. 95:664–670.[CrossRef][Web of Science]
- Tamulonis, J.P., B.M. Luzzi, R.S. Hussey, W.A. Parrott, and H.R. Boerma. 1997c. RFLP mapping or resistance to southern root-knot nematode in soybean. Crop Sci. 37:1903–1909.[Abstract/Free Full Text]
- Tang, S., A. Leon, W.C. Bridges, and S.J. Knapp. 2006. Quantitative trait loci for genetically correlated seed traits are tightly linked to branching and pericarp pigment loci in sunflower. Crop Sci. 46:721–734.[Abstract/Free Full Text]
- Tasma, I.M., and R.C. Shoemaker. 2003. Mapping flowering time gene homologs in soybean and their association with maturity (E) loci. Crop Sci. 43:319–328.[Abstract/Free Full Text]
- Vierling, R.A., J. Faghihi, V.R. Ferris, and J.M. Ferris. 1996. Association of RFLP markers conferring broad-based resistance to the soybean cyst nematode (Heterodera glycines). Theor. Appl. Genet. 92:83–86.[CrossRef][Web of Science]
- Webb, D.M. 2003. Quantitative trait loci associated with cyst nematode resistance and uses thereof. U.S. Patent 6538,175. Date issued: 25 March.
- Webb, D.M., B.M. Baltazar, A.P. Rao-Arelli, J. Schupp, K. Clayton, P. Keim, and W.D. Beavis. 1995. Genetic mapping of soybean cyst nematode race-3 resistance loci in the soybean PI 437.654. Theor. Appl. Genet. 91:574–581.[Web of Science]
- Williams, C., and W. Birchfield. 1974. Controlling reniform nematodes in soybeans. La. Agric. 17:14–15.
- Williams, C., D.F. Gilman, D.S. Fontenot, and W. Birchfield. 1981. Inheritance of reaction to the reniform nematode in soybean. Crop Sci. 21:93–94.[Abstract/Free Full Text]
This article has been cited by other articles:

|
 |

|
 |
 
T. Ishii, T. Hayashi, and K. Yonezawa
Optimization of the Marker-based Procedures for Pyramiding Genes from Multiple Donor Lines: III. Multiple-Gene Assemblage Using Background Marker Selection
Crop Sci.,
November 24, 2008;
48(6):
2123 - 2131.
[Abstract]
[Full Text]
[PDF]
|
 |
|

|
 |

|
 |
 
Y. Gao, R.M. Biyashev, M.A. S. Maroof, N.M. Glover, D.M. Tucker, and G.R. Buss
Validation of Low-Phytate QTLs and Evaluation of Seedling Emergence of Low-Phytate Soybeans
Crop Sci.,
July 1, 2008;
48(4):
1355 - 1364.
[Abstract]
[Full Text]
[PDF]
|
 |
|