Mycologia
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS

This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Services
Right arrow Similar articles in this journal
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Ceresini, P. C.
Right arrow Articles by Cubeta, M. A.
Right arrow Search for Related Content
PubMed
Right arrow Articles by Ceresini, P. C.
Right arrow Articles by Cubeta, M. A.
Agricola
Right arrow Articles by Ceresini, P. C.
Right arrow Articles by Cubeta, M. A.
Mycologia, 94(3), 2002, pp. 450-460.
© 2002 by The Mycological Society of America

Genetic structure of populations of Rhizoctonia solani AG-3 on potato in eastern North Carolina


Paulo C. Ceresini

     Departamento de Biologia, Faculdade de Engenharia, Universidade Estadual Paulista "Júlio de Mesquita Filho" (UNESP), 15385-000, Ilha Solteira, São Paulo, Brazil

H. David Shew

     Department of Plant Pathology, North Carolina State University, Raleigh, North Carolina 27695

Rytas J. Vilgalys

     Department of Biology, Duke University, Durham, North Carolina 27708

U. Liane Rosewich

     USDA-ARS, Cereal Disease Laboratory, University of Minnesota, Saint Paul, Minnesota 55108

Marc A. Cubeta 1

     Department of Plant Pathology, North Carolina State University, Raleigh, North Carolina 27695

    ABSTRACT
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 LITERATURE CITED
 

A polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) method was developed to identify and differentiate genotypes of Rhizoctonia solani anastomosis group 3 subgroup PT (AG-3 PT), a fungal pathogen of potato. Polymorphic co-dominant single-locus PCR-RFLP markers were identified after sequencing of clones from a genomic library and digestion with restriction enzymes. Multilocus genotypes were determined by a combination of PCR product and digestion with a specific restriction enzyme for each of seven loci. A sample of 104 isolates from one commercial field in each of five counties in eastern North Carolina was analyzed, and evidence for high levels of gene flow between populations was revealed. When data were clone-corrected and samples pooled into one single North Carolina population, random associations of alleles were found for all loci or pairs of loci, indicating random mating. However, when all genotypes were analyzed, the observed genotypic diversity deviated from panmixia and alleles within and between loci were not randomly associated. These findings support a model of population structure for R. solani AG-3 PT on potato that includes both recombination and clonality.

Key words: migration, population subdivision, Thanatephorus cucumeris


    INTRODUCTION
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 LITERATURE CITED
 
The soilborne basidiomycete Rhizoctonia solani Kühn [Thanatephorus cucumeris (Frank) Donk] anastomosis group 3 subgroup PT (AG-3 PT) is a major pathogen infecting potato (Solanum tuberosum L.) worldwide. Symptoms and signs of Rhizoctonia disease of potato can be found on many plant parts including roots, sprouts, stems, stolons, and tubers. As potato plants reach maturity, the fungus can infest daughter tubers, producing black sclerotia on their surface (Carling et al 1989Citation ).

The ecology and pathology of R. solani AG-3 PT on potato has been extensively studied (Ogoshi 1987Citation ), but little information is available concerning the role of asexual and sexual reproduction in determining population structure of this important potato pathogen. In general, R. solani AG-3 PT is assumed to be asexual, surviving primarily as mycelium and sclerotia in soil and on potato seed tubers (Cubeta and Vilgalys 2000Citation ). However, the sexual stage (T. cucumeris) is frequently observed on potato stems in commercial production fields. Although basidiospores of T. cucumeris are assumed to have limited dispersal and little or no importance in the disease cycle (Carling et al 1989Citation ), they are an important aspect of fungal life history that contributes to the genetic diversity and structure of field populations. The mating system of T. cucumeris (anamorph = R. solani AG-3 PT) is not known. Likewise, it is not known whether recombination associated with sexual reproduction occurs among field isolates (Vilgalys and Cubeta 1994Citation , Cubeta and Vilgalys 1997Citation ).

Except for the studies of Rosewich et al (1999)Citation , no information about the genetic diversity and structure of field populations of R. solani exists. Previous studies with R. solani AG-3 have employed various methods to examine the genetic variation among isolates from culture collections obtained from distinct Solanaceous hosts, including eggplant, potato, tobacco and tomato (Laroche et al 1992Citation , Liu and Sinclair 1992Citation , Liu et al 1993Citation , Stevens Johnk et al 1993Citation ). However, basic questions remain about the population structure of R. solani AG-3 PT. For example: i) How diverse are populations of R. solani AG-3 PT from potato fields? and ii) Are populations of R. solani AG-3 PT on potato clonal or panmictic?

In this study two null hypotheses were tested: a) populations of R. solani AG-3 PT from potato in North Carolina (NC) have no genetic subdivision or differentiation; and b) populations of R. solani AG-3 PT from potato plants are randomly mating (recombining). To address these hypotheses, a strategy for rapid detection of individual genotypes of R. solani AG-3 PT using co-dominant polymerase chain reaction (PCR)-RFLP markers was employed. These markers were used to assess the genetic diversity and structure of population samples of R. solani AG-3 PT obtained from one field from each of five different counties in NC. Analysis of genetic variation (e.g., similarity) and random mating among these samples were conducted. Data were also analyzed to determine whether populations of R. solani AG-3 PT were subdivided.


    MATERIALS AND METHODS
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 LITERATURE CITED
 
Population sampling and DNA isolation – Isolates of R. solani AG-3 PT for this study were obtained from infected subterranean potato stems with cankers. Infected plants were obtained from one commercial field in each of 5 counties in eastern NC (Ceresini 2000Citation ). Ten plots (2 per field), each consisting of 6–8 rows of 9 m length were planted with either 240 or 320 potato seed pieces (cv. ‘Atlantic’). Seed potatoes used to establish the plots were disinfested before planting with 2% formaldehyde to eliminate seedborne R. solani AG-3 PT (Carling et al 1989Citation ). All plants from each plot were harvested and examined for symptoms of Rhizoctonia stem canker.

Samples were assayed on alkaline water agar (AWA, pH 8.5) to isolate R. solani AG-3 PT (Gutierrez et al 1997Citation ). A total of 104 pure cultures of R. solani were established by transferring mycelium from single hyphal tips to potato dextrose agar. Isolates were stored on sterilized rice grains at -20 C. Isolates were examined for nuclear condition by the DAPI staining procedure (Martin 1987Citation ). The anastomosis grouping of each isolate of R. solani was determined by standard protocols (Gutierrez et al 1997Citation , Herr and Roberts 1980Citation ). Genomic DNA was extracted by methods described previously (Rosewich et al 1999Citation ).

Development of PCR-RFLP markers – Anonymous single-copy DNA markers were developed using a strategy outlined by Karl and Avise (1993)Citation and Xu et al (1999)Citation . A genomic library was constructed as described previously (Rosewich et al 1999Citation ) with Hind III-digested DNA of a randomly chosen NC isolate of R. solani AG-3 PT (P-211) cloned into pUC18 (Rosewich et al 1999Citation ). Sizes of clones resulting from this library were determined by digestion with Hind III and separation of DNA on a 0.8% agarose gel with a 1-kb ladder. Clones with a single insert between 1 and 1.8 kb were selected (32 total). An extra clone with a 3.2 kb insert was also included in the sample. Both termini from 15 cloned DNA fragments were sequenced with forward and reverse primers for the multiple cloning site of pUC18 in an ABI Prism 377 automated sequencer (PE Applied Biosystems). Sequence chromatograms were compiled using Sequencher software (vers. 2.0, GeneCodes Corp.). Locus-specific PCR primers for 21–26 nucleotides were designed for further screening. PCR reactions were conduced in a Model 9600 DNA ThermoCycler (Perkin-Elmer) employing standard conditions described previously for R. solani (Vilgalys and Gonzalez 1990Citation ) with Amplitaq DNA polymerase (Perkin-Elmer). Restriction polymorphisms were detected initially by screening 21 random isolates of R. solani AG-3 PT from all five NC counties with restriction enzymes with 6-base recognition sites (Bam HI, Eco RI, Eco RV, Hind III, Hinc II, Hae II, Pst I, Sac I, Sty I, Xho I) or restriction enzymes with 4-base recognition sites (Hha I, Hae III, Mbo I, Msp I, Rsa I, Tru 9I). The entire sample of 104 isolates was then screened with combinations of PCR fragments and specific restriction enzymes that generated polymorphisms. Typical restriction enzyme reactions consisted of 5 µL of PCR product, 7.5 U of enzyme (Promega), 2.0 µL of 10x reaction buffer, 0.2 µL of 10x bovine serum albumin and H2O to a total volume of 20 µL. Reactions were incubated for 5 h at 37 C (or according to manufacturer's recommendations) then stopped by incubation at 70 C for 5 min. Restriction digests of PCR products were separated by electrophoresis in 1.2 to 2.4% MetaPhor (FMC BioProducts) agarose gels in 1x TBE, containing ethidium bromide, and viewed by UV transillumination.

Data analysis – Three data sets were constructed: (a) all isolates and retaining the five geographical populations (NC counties, uncorrected data); (b) clone-corrected data; isolates with the same multilocus PCR-RFLP genotype within each geographic population considered only once; and (c) either uncorrected or clone-corrected data of the combined isolates from all populations.

For analysis of genetic diversity within R. solani AG-3 PT populations, three population genetics measures were applied: (a) Stoddart's genotypic diversity (Stoddart and Taylor 1988Citation ) Go = 1/({Sigma}pi2), where pi is the frequency of a particular multilocus genotype, with Go ranging from 1 to N, where N is the sample size. The sample-corrected genotypic diversity, Gc = Go/N, was used to allow comparisons between different populations. The expected overall genotypic diversity for a sample of size N was calculated using Stoddart's expected genotypic diversity (Stoddart 1983Citation ) GE = 1/(d + p/N), where d = {Sigma}gi2 for all genotypes when (gi* N) >= 1 and p = {Sigma}gi for all the genotypes when gi* N) < 1; gi is the expected frequency of the multilocus genotype calculated by Hardy-Weinberg expectations; (b) observed heterozygosity (HO) which represents the percentage of heterozygosity at each locus (Hartl and Clark 1997Citation ). The mean observed heterozygosity (i.e., gene diversity) was calculated as the arithmetic mean of all loci sampled; and (c) unbiased estimate of expected heterozygosity (HE) or mean unbiased gene diversity, HE = {Sigma}rk=1 hk/r, where hk is the value of h for kth locus, and r is the total number of loci studied; h is an unbiased estimate of heterozygosity for a single locus given by h = 2n(1 - {Sigma}xi2)/(2n - 1), where n is the number of diploid individuals, and xi is the corresponding frequency of the ith allele at a locus in a sample from the population (Levene 1949Citation , Nei 1978, 1987Citation ). The expected heterozygosity was calculated with TFPGA version 1.3 (Miller 1997Citation ).

Between-population analyses were conducted to test for geographic structure in populations of R. solani AG-3 PT. Three population genetics analyses were performed: (a) analysis of molecular variance (AMOVA), which estimates variance components considering the number of differences between molecular genotypes (Weir and Cockerham 1984Citation , Excoffier et al 1992Citation , Weir 1996Citation ); (b) pairwise {Phi} statistics comparisons (F statistics equivalents produced by AMOVA), which reveals similarity between the individual geographic populations. AMOVA analyses were conducted using ARLEQUIN ver 2.000 (Schneider et al 2000Citation ); and (c) similarity among the population samples calculated with F statistics (theta P or FWC(ST) as in Weir and Cockerham (1984)Citation ) and Nei's unbiased genetic distance (D) (Nei 1978, 1987Citation ). These tests (FWC(ST) and D), which are based on allele frequency differences between population samples, were calculated with GDA version 1.0 (d15) (Lewis and Zaykin 2000Citation ).

To compare the relative contribution of clonality and recombination in each sample of R. solani AG-3 PT, two population genetic measures were used: (a) the Hardy-Weinberg equilibrium (HWE) test, which examines the random association of alleles within a locus by testing the observed genotypic frequencies for goodness-of-fit with Hardy-Weinberg expectations (Hartl and Clark 1997Citation ) and (b) composite genotype disequilibrium or linkage disequilibrium (Weir 1996Citation ), which estimates the allelic association between pairs of polymorphic markers by an exact test for association through permutation (Zaykin et al 1995Citation ). To test for HWE, an exact test analogous to Fisher's exact test (Guo and Thompson 1992Citation ) was calculated with ARLEQUIN 2.000 (Schneider et al 2000Citation ). In addition, the likelihood ratio statistic G2 (recommended for genotypic data for which the haplotypic phase is unknown) was calculated with POPGENE (Yeh et al 1997Citation ). Linkage disequilibrium between all pairs of loci was calculated by an exact test analogous to Fisher's exact test (Slatkin 1994Citation , Slatkin and Excoffier 1996Citation ). Calculations were performed with ARLEQUIN 2.000 (Schneider et al 2000Citation ).


    RESULTS
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 LITERATURE CITED
 
Development of PCR-RFLP markers – From the total of 32 randomly cloned Hind III-digested R. solani AG-3 PT genomic DNA fragments selected based on size, 15 were partially sequenced and used to design locus-specific PCR primers. Eight fragments were found to be monomorphic, hypervariable or otherwise problematic (e.g., failing to amplify, producing multiple PCR fragments). The remaining seven DNA fragments could be re-amplified from all isolates and possessed unique restriction polymorphisms. Each of four fragments had one polymorphic restriction site, whereas three had two polymorphic sites. Each polymorphic restriction site, which corresponded to a unique combination of primer pair, restriction enzyme, and restriction site, was defined as a locus (Table I ). Scoring by this method is unambiguous as each locus has potentially two alleles: allele 1 represents the absence of the restriction site, and allele 2 its presence. There were two cases where a third allele was present (i.e., markers pP42-Hind III and pP83-Hind III). When the restriction enzyme Hind III was used for digestion, allele 1 corresponded to the original (undigested PCR product) fragment size. However, for other restriction enzymes, allele 1 or 2 was smaller than the original fragment. Examples of the PCR-RFLP genotyping and allele designation methods are presented in Fig. 1 .


View this table:
[in this window]
[in a new window]
 
TABLE I. Polymorphic PCR-RFLP markers for analysis of Rhizoctonia solani AG-3 PT from potato

 


View larger version (28K):
[in this window]
[in a new window]
 
  FIG. 1. PCR restriction fragment length polymorphisms (RFLPs) of Rhizoctonia solani AG-3 PT total DNA amplified with the locus-specific primers for pP09(F09/R12), pP42(F12/R09), pP45(F09/R03), pP46(F09/R01), pP47(F08/R07), pP83(F05/R04), and pP89(F02/R10) markers showing homozygous (1/1, 2/2, and 3/3) and heterozygous (1/2) genotypes for each locus. The approximate size of each PCR fragment is presented in base pairs (bp) and was determined by comparing each migrating fragment with a DNA standard (either 100 or 1000 bp ladder, Promega). Polymorphic restriction sites for each locus are indicated in bold

 
Genotypic analysis – Genotypes and genotypic counts for each of the five samples of Rhizoctonia solani AG-3 PT are presented in Table II . Among the 104 R. solani isolates analyzed, 32 multilocus PCR-RFLP genotypes (MRG) were identified. On average, 8.8 MRG were found per locality (range 2–15). Most of the MRG (about 81%) were site specific (Tables II and III ). However, isolates with the same genotype were also commonly observed between field sites. Five multilocus genotypes (MRG 3, 2, 1, 48 and 50) represented 13.5, 9.6, 8.7, 8.7, and 6.7% of the sample, respectively. Altogether, these five genotypes represented isolates, which accounted for 47.1% of the total population sample. The observed count of 14 for MRG 3, and 9 for either MRG 1 or MRG 48 was higher than the expected counts of 5.28 for both MRG 3 and MRG 1, and of 3.82 for MRG 48, which were derived under the hypothesis of panmixia (Stoddart and Taylor 1988Citation ). The observed and expected counts of genotypes for MRG 2 and 50 did not differ statistically.


View this table:
[in this window]
[in a new window]
 
TABLE II. Frequency of multi-locus PCR-RFLP genotypes (MRG) of Rhizoctonia solani AG-3 PT from potato fields in North Carolina, and estimates of overall observed (GO) and expected (GE*) genotypic diversity

 

View this table:
[in this window]
[in a new window]
 
TABLE III. Genotypic diversity of five population samples of Rhizoctonia solani AG-3 PT in North Carolina

 
As a control, a sample of the same potato seed tubers used for establishing the plots was subjected to isolation on AWA, before disinfestation with 2% formaldehyde. Six additional MRG were identified on these infested seed potatoes: one MRG 16 (1/1 1/1 1/1 1/2 1/1 1/1 2/2), two MGR 17 (1/1 1/1 1/1 1/2 1/1 1/2 2/2), one MRG 28 (1/1 1/2 1/1 1/1 1/1 1/1 2/2), one MRG 37 (1/1 1/2 1/1 2/2 1/1 2/2 2/2), one MRG 38 (1/1 1/2 1/2 1/1 1/1 1/1 2/2), and two MRG 59 (1/2 1/1 1/2 1/1 1/1 1/1 1/2). No isolates of R. solani AG-3 PT were recovered from potato seed tubers after disinfestation.

Stoddart's genotypic diversity (Go) and the sample size-corrected genotypic diversity (Gc) of R. solani AG-3 PT are presented in Table III . The mean genotypic diversity (Gc) was 0.2945. Overall, there was no difference between any pairs of samples considering their Gc values. In addition, the observed overall genotypic diversity (Go for uncorrected data = 16.64) was smaller than the expected value under panmixia (GE* = 38.78) (Table II ).

Gene diversity, patterns of heterozygosity and deviations from HWE – Six of the seven loci surveyed from five geographical populations of R. solani AG-3 PT were polymorphic at all localities, with an average of 1.97 alleles per locus. The only exception was the locus pP45-Hind III, which was polymorphic only at Hyde County. Most alleles (11/16) were present in all five populations. A rare allele was present at low frequencies in one population (allele 2 at locus pP45-Hind III in the population sample from Hyde).

The observed and expected heterozygosity of all samples from NC are presented in Table IV . Heterozygosity values varied greatly among loci, from 0.010 to 0.481. Two of the seven loci had heterozygosity less than 0.100 (pP45-Hind III and pP47-Eco RV). The mean heterozygosity (over loci) of the samples was between 0.186 and 0.381 when all the isolates were included in the analysis and between 0.229 and 0.357 for clone-corrected data. Comparisons of mean heterozygosity of each sample indicated no difference between expected and observed values.


View this table:
[in this window]
[in a new window]
 
TABLE IV. Patterns of heterozygosity in different samples of Rhizoctonia solani AG-3 PT in North Carolina

 
Each sample of R. solani AG-3 PT from NC was also examined for the association of alleles within a locus (HWE tests) for each particular locus (Table IV ). For all samples, when data were not clone corrected, most loci (at least 5/7) had genotypic counts that did not deviate from Hardy-Weinberg expectations (HWE). Two loci (pP09-Hha I and pP42-Hind III) in the sample from Currituck, one locus (pP83-Hind III) in the sample from Hyde, and one locus (pP46-Hind III) in the sample from Camden showed significant deviation from HWE. After clone correction, all loci had genotypic counts that did not deviate from HWE.

Population structure – Analysis of population structure was conducted for both uncorrected and clone-corrected data. Using AMOVA (Analysis of Molecular Variance), overall {Phi}ST between geographic populations was 0.0745 for uncorrected data. When the geographical populations were subdivided into sampling units within geographic populations (counties), subdivision was found to be even slightly higher between sampling units than between geographic populations ({Phi}ST = 0.1076). For clone-corrected data, overall {Phi}ST was -0.0001 between geographic populations. With further hierarchical subdivision of geographical populations into sampling units, the {Phi}ST value was again observed to be slightly higher (0.0111). In general, most of the molecular variation was detected within populations, with very little variation between geographic populations.

Population differentiation and gene flow were also analyzed by estimates of theta P (FST) as described by Weir and Cockerham (1984)Citation (Table V ). Overall FST between populations was 0.0741 for uncorrected data, which was significantly greater than zero. Nevertheless, this value of FST also indicates little genetic differentiation among populations. For clone-corrected data overall FST was -0.0009, which was not significantly different from zero.


View this table:
[in this window]
[in a new window]
 
TABLE V. Population differentiation and gene flow of Rhizoctonia solani AG-3 PT in North Carolina

 
Pairwise comparison between geographic populations was performed using estimates of FST equivalents (pairwise {Phi}ST values) and Nei's unbiased genetic distances (D), for both uncorrected and clone-corrected data (Table VI ). Pairwise {Phi}ST values were significantly greater than zero, which would indicate population differentiation, only for uncorrected data. For example, the population from Washington was significantly different from Tyrrell and Currituck and very similar to Hyde and Camden county populations. For clone-corrected data, however, there was no indication of population differentiation considering both pairwise {Phi}ST values and D.


View this table:
[in this window]
[in a new window]
 
TABLE VI. Population pairwise {Phi}STa of Rhizoctonia solani AG-3 PT isolates collected from five commercial potato fields in North Carolina, for clone-corrected data (above diagonal) and total data (below diagonal)

 
HWE test for association of alleles within loci and exact test for linkage disequilibrium between pairs of loci – Based on the lack of geographic subdivision, all MRG data were pooled into a single data set (for both uncorrected and clone-corrected samples) to test for HWE (Table VII ) and composite linkage disequilibrium (Table VIII ) between all pairs of loci.


View this table:
[in this window]
[in a new window]
 
TABLE VII. Test for Hardy-Weinberg equilibrium for multilocus PCR-RFLP genotypes of Rhizoctonia solani AG-3 PT

 

View this table:
[in this window]
[in a new window]
 
TABLE VIII. Exact tests for linkage disequilibriuma among pairs of individual PCR-RFLP loci based on 32 multilocus genotypes

 
The results from the likelihood ratio test (G2) and the exact test analogous to Fisher's indicated genotypic proportions similar to HWE proportions at all seven loci (for clone-corrected data) and for six loci considering the uncorrected sample. For uncorrected data, only the locus pP09-Hha I showed proportions of genotypic counts that deviated from HWE. Overall, for uncorrected data, most pairs of loci were in genotypic disequilibrium (13/21). However, for clone-corrected data, the exact test for linkage disequilibrium was not significant for all pairs of loci.


    DISCUSSION
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 LITERATURE CITED
 
This study describes a PCR-RFLP method for genotyping isolates of R. solani AG-3 PT. We employed a combination of two strategies for development of co-dominant markers as proposed by Rosewich et al (1999)Citation and Xu et al (1999)Citation . The method established here enables unambiguous scoring of genotypes of isolates of R. solani AG-3 PT by allowing the distinction between homozygotes and heterozygotes. The PCR-RFLP method is also efficient, as only small amounts of DNA are required as starting material. Consistent with the heterokaryotic nature of R. solani, which has been confirmed by a similar molecular typing method (Rosewich et al 1999Citation ), every isolate had one or more heterozygous loci (Table II ).

During the process of screening, four of the seven PCR-RFLP markers used for genotyping R. solani AG-3 PT (pP42, pP46, pP47 and pP89) also amplified genomic DNA from R. solani AG-3 TB, a pathogen of tobacco. Except for marker pP42, none of the alleles present in the sample of 20 tobacco isolates of R. solani AG-3 TB were shared with the potato sample (data not shown). Our findings corroborate previous information that populations of R. solani AG-3 are subdivided by host specialization (Ceresini et al 2001Citation , Kuninaga et al 1997Citation , Kuninaga et al 2000Citation , Ogoshi 1987Citation , Shew and Melton 1995Citation , Stevens Johnk et al 1993Citation ). These four genetic markers could have potential applicability for phylogenetic analysis within R. solani AG-3 as they were shown to be conserved enough to amplify genes across AG-3 PT and AG-3 TB, two genetically distinct subgroups of AG-3.

The first objective of this study was to determine the extent of genetic differentiation among populations of R. solani AG-3 PT and make inferences about gene flow and population subdivision. The allelic information from the seven single locus PCR-RFLP markers allowed us to group 104 isolates of R. solani AG-3 PT into 32 MRG. Gene flow was high, even though relatively few MRG were shared among localities. The overall {Phi}ST between geographic populations was 0.0745 for total data. When data were subdivided into sampling locations (plots) within a field, subdivision was found to be even slightly higher between sampling units than between geographic populations ({Phi}ST = 0.1076). However, using clone-corrected data, {Phi}ST and FST estimates were essentially zero. In addition, for clone-corrected data, all pairwise {Phi}ST values were not significant at 5%, and Nei's genetic distance (D) indicated no differentiation between populations.

Our results support the existence of a mechanism that ensures efficient long distance dispersal of genotypes between geographic populations. Recent analysis of the contribution of tuber-borne genotypes to the overall population of R. solani AG-3 PT present in potato fields of North Carolina provided evidence on how tuber-borne inoculum introduced every year affects population structure of the pathogen. High levels of genetic diversity were detected in introduced potato tubers from Canada (New Brunswick), Maine, and Wisconsin (50 different genotypes among a total sample of 58 isolates) (Ceresini 2000Citation ). Nevertheless, seven of these genotypes were among the most frequent genotypes observed in North Carolina fields (Ceresini 2000Citation ).

A similar lack of population structure was found for R. solani AG-1-IA from rice in Texas (Rosewich et al 1999Citation ). The authors characterized seven co-dominant RFLP loci for 182 AG-1 IA isolates in Texas rice fields and found 36 different multilocus genotypes. As there was no evidence for population subdivision, all isolates from the rice growing areas in Texas were considered to be part of one single population of R. solani AG-1 IA. Considering the similar observations of lack of population structure within R. solani AG-3 PT, we also pooled all multilocus PCR-RFLP genotypes into one single data set for further analysis of HWE and linkage disequilibrium.

The other specific question addressed in this study was: Does multilocus genetic diversity analysis support a sexual structure in local populations of R. solani AG-3 PT? Several studies have provided information about the structure of R. solani populations on wheat in Australia and rice in Texas. For example, MacNish et al (1997)Citation found an association of anastomosis reaction type with allozymes and an overrepresentation of certain phenotypes, which was considered evidence for clonality. Evidence for clonality was also demonstrated for R. solani AG-1-IA through findings of repeated recovery of genotypes in the same and different rice field (Rosewich et al 1999Citation ).

To answer the question on population structure of R. solani AG-3 PT from potato, we considered the main effects of asexual reproduction on the pathogen population. The first major effect of asexual reproduction is the production of fewer recombinant genotypes, which decreases genotypic diversity. In this study the hypothesis of clonality was partially supported in that at least three of the population samples from potato shared the most common multilocus genotypes (MRG 1, 3 and 48) (Table II ). In addition, observed counts for MRG 1, MRG 3, or MRG 48 were much higher than the expected counts derived under the hypothesis of panmixia (Stoddart and Taylor 1988Citation ). The higher frequency of these genotypes in potato samples may be the consequence of asexual reproduction, since all three MRGs were also found on seed potato tubers introduced from Maine into NC (Ceresini 2000Citation ). In addition, if populations of R. solani AG-3 PT were totally panmictic, the observed genotypic diversity would be similar to the expected value under panmixia (Stoddart 1983Citation ). Our results indicated that the overall genotypic diversity (Go for uncorrected data = 16.64) was considerably smaller than the expected value under panmixia (GE* = 38.78) (Table II ), supporting clonality. The clonal reproduction of R. solani AG-3 PT via mycelia or sclerotia provides a mechanism for a local collection of clones to predominate within a much larger recombining population of the pathogen on potato. However, the biological features of these asexual structures (mycelia and scletoria), especially the mode of dissemination, will influence the scale and intensity of clonal spread (Taylor et al 1999Citation ). A clonal structure for a population of R. solani AG-3 PT in South Australia could not be rejected because one DNA fingerprint pattern represented 30% of the isolates tested (Balali et al 1996Citation ).

The repeated occurrence of genotypes of R. solani AG-3 PT clearly suggests that populations are not only randomly mating (Milgroom 1996Citation ). However, the multiple occurrence of genotypes can be a confounding factor in the analyses of population structure of this potato pathogen. Therefore, the approach of analyzing the population structure of R. solani AG-3 PT considering only one representative of each multilocus genotype is appropriate. This approach has been recommended for making inferences about recombination in other populations of plant pathogenic fungi (Liu et al 1996Citation , McDonald et al 1994Citation , Rosewich et al 1999Citation ).

The second major effect of recombination on population structure is the random association of alleles both within (HWE) and between loci (linkage equilibrium). When all genotypes were included in the analysis, only one locus (pP09 Hha I) was not in HWE. In addition, linkage disequilibrium was detected for eight pairs of loci. However, for clone-corrected data, both tests for HWE indicated that all loci were in equilibrium and no linkage disequilibrium was detected for any pair of loci. This emphasizes the importance of using clone-corrected data to better understand aspects of reproduction that contribute to population structure (McDonald et al 1994Citation , Rosewich et al 1999Citation ). Comparatively, analyses of the population structure of R. solani AG-1-IA from rice indicated that four of seven RFLP loci were in HWE and only one of all possible pairs of loci was in linkage disequilibrium, which indicated that recombination was occurring (Rosewich et al 1999Citation ).

At the extremes of models of population structure, populations may be panmictic or strictly clonal. An intermediate model was proposed by Maynard Smith et al (1993)Citation for bacterial pathogens. This model may also be appropriate for many fungal plant pathogens that have annual sexual cycles and asexual phases. This model is characterized by frequent recombination and by the occurrence of one or a few successful individual genotypes that reproduce clonally and increase to high frequencies. Differences in fitness among genotypes expressed in the asexual phase could be responsible for the prevalence of certain genotypes. Our observations support a model of population structure for R. solani AG-3 PT that includes both clonality and recombination.

This study provided initial information on population genetics of R. solani AG-3 PT in NC. Further studies are needed to elucidate the magnitude of recombination in populations of R. solani AG-3 PT. In addition, increased sampling is needed to maximize genotypic diversity and enhance the power of statistical tests to infer population genetic processes. Increased resolution to infer population structure also may be achieved by the examination of additional genetic loci to ascertain that MRGs identified in this study actually represent clones.

Consideration also should be given to the elucidation of the role of other evolutionary processes that could cause linkage disequilibrium in populations of R. solani AG-3 PT, such as gene flow or migration. Gene flow can cause gametic disequilibrium by the immigration of individuals from populations with different allele frequencies (Milgroom 1996Citation ). Future research also should examine populations of R. solani AG-3 PT from the northern United States (Wisconsin and Maine) and eastern Canada, as these areas are responsible for the production of most of the seed potatoes imported into eastern North Carolina.


    ACKNOWLEDGMENTS
 
We thank Dr. H. Corby Kistler (USDA/ARS, Cereal Diseases Lab, St. Paul, Minnesota), Steven Leath, (USDA/ARS, Dept. of Plant Pathology, NCSU), Tim James, and Daniel Snyder (Dept. of Biology, Duke University,). Special thanks for Pamela Puryear (Tobacco Literature Service, NCSU) for locating literature. This research was supported in part with a research assistantship from CNPq (Conselho Nacional deDesenvolvimento Científico e Tecnológico)—Brasília, DF, Brazil (200873/92–8) and with a graduate student fellowship and E.E. Butler travel award from The Mycological Society of America awarded to PC. A leave of absence was also granted to PC by UNESP (Universidade Estadual Paulista "Júlio de Mesquita Filho"), Ilha Solteira, São Paulo, Brazil, for attending the Ph.D. program on Plant Pathology at North Carolina State University. Special thanks to Drs. Ana Maria R. Cassiolato and Marli F. Stradioto Papa for assuming the teaching responsibilities of PC at UNESP during his absence.


    FOOTNOTES
 
1 Corresponding author, marc_cubeta{at}ncsu.edu Back

Accepted for publication October 26, 2001.


    LITERATURE CITED
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 LITERATURE CITED
 
Balali GR, Whisson DL, Scott ES, Neate SM., 1996 DNA fingerprinting probe specific to isolates of Rhizoctonia solani AG-3 Mycol Res 100:467-470

Carling DE, Leiner RH, Westphale PC., 1989 Symptoms, signs, and yield reduction associated with rhizoctonia disease of potato induced by tuberborne inoculum of Rhizoctonia solani AG-3 Am Potato J 66:693-702

Ceresini PC., 2000 Population biology and genetics of Rhizoctonia solani anastomosis group 3 (AG-3) Dept. of Plant Pathology, North Carolina State University. [PhD Dissertation]. 89 p

Ceresini PC., Shew HD, Vilgalys R, Cubeta MA., 2001 Genetic diversity of Rhizoctonia solani anastomosis AG-3 from potato and tobacco in North Carolina Mycologia 94: (in press)

Cubeta MA, Vilgalys R., 1997 Population biology of the Rhizoctonia solani complex Phytopathology 87:480-484

Cubeta MA, 2000 Rhizoctonia In. Lederberg J, ed. Encyclopedia of microbiology Vol. 4. San Diego: Academic Press. pp. 109–116

Excoffier L, Smouse P, Quattro J., 1992 Analysis of molecular variance inferred from metric distances among DNA haplotypes: application to human mitochondrial restriction data Genetics 131:479-491[Abstract]

Guo SW, Thompson EA., 1992 Performing the exact test of Hardy-Weinberg proportions for multiple alleles Biometrics 48:361-372[Medline]

Gutierrez WA, Shew HD, Melton TA., 1997 Sources of inoculum and management for Rhizoctonia solani damping-off on tobacco transplants under greenhouse conditions Plant Disease 81:64-66

Hartl DL, Clark AG., 1997 Principles of population genetics. 3rd ed Sunderland, Massachusetts: Sinauer. 542 p

Herr LJ, Roberts DL., 1980 Characterization of Rhizoctonia solani populations obtained from sugarbeet fields with differing soil textures Phytopathology 70:476-480

Karl SA, Avise JC., 1993 PCR-based assays of Mendelian polymorphisms from anonymous single copy nuclear DNA: techniques and applications for population genetics Mol Biol and Evol 10:342-361

Kuninaga S, Natsuaki T, Takeuchi T, Yokosawa R., 1997 Sequence variation of the rDNA ITS regions within and between anastomosis groups in Rhizoctonia solani Curr Genetics 32:237-243[Medline]

Kuninaga S, Carling DE, Takeuchi T, Yokosawa R., 2000 Comparison of rDNA-ITS sequences from potato and tobacco strains in Rhizoctonia solani AG-3 J Gen Plant Pathol 66:2-11

Laroche JP, Jabaji-Hare SH, Charest PM., 1992 Differentiation of two anastomosis groups of Rhizoctonia solani by isozyme analysis Phytopathology 82:1387-1393

Levene H., 1949 On a matching problem arising in Genetics Ann Math Stat 20:91-94

Lewis PO, Zaykin D., 2000 Genetic data analysis: computer program for the analysis of allelic data. Version 1.0 (d15) Free program distributed by the authors over the internet from the GDA Home Page at http://alleyn.eeb.uconn.edu/gda

Liu Y-C, Cortesi P, Double ML, McDonald WL, Milgroom MG., 1996 Diversity and multilocus genetic structure in populations of Cryphonectria parasitica Phytopathology 86:1344-1351

Liu ZL, Domier LL, Sinclair JB., 1993 ISG-specific ribosomal DNA polymorphism of the Rhizoctonia solani species complex Mycologia 85:795-800

Liu ZL, Sinclair JB., 1992 Genetic diversity of Rhizoctonia solani anastomosis group 2 Phytopathology 82:778-787

McDonald BA, Miles J, Nelson LR, Pettway RE., 1994 Genetic variability in nuclear DNA in field populations of Stagonospora nodorum Phytopathology 84:250-255

MacNish GC, Carling DE, Brainard KA., 1997 Relationship of microscopic and macroscopic vegetative reactions in Rhizoctonia solani and the occurrence of vegetatively compatible populations (VCPs) in AG-8 Mycol Res 101:61-68

Martin B., 1987 Rapid tentative identification of Rhizoctonia spp. associated with diseased turfgrasses Plant Dis 71:47-49

Maynard Smith J, Smith NH, O'Rourke M, Spratt BG., 1993 How clonal are bacteria? Proc Natl Acad Sci USA 90:4384-4388[Abstract/Free Full Text]

Milgroom MG., 1996 Recombination and the multilocus structure of fungal populations Annu Rev Phytopathol 34:457-477[Medline]

Miller MP., 1997 Tools for Population Genetic Analyses (TFPGA): a Windows® program for analysis of allozyme and molecular population genetic data. Version 1.3 Free program distributed by the author over the Internet from the website http://herb.bio.nau.edu/~miller

Nei M., 1978 Estimation of average heterozygosity and genetic distance from a smaller number of individuals Genetics 89:583-590[Abstract/Free Full Text]

Nei M., 1987 Molecular evolutionary genetics New York: Columbia University Press. 512 p

Ogoshi A., 1987 ECOLOGY and pathogenicity of anastomosis and intraspecific groups of Rhizoctonia solani Kühn Annu Rev Phytopathol 25:125-143

Rosewich UL, Pettway RE, McDonald BA, Kistler HC., 1999 High levels of gene flow and heterozygote excess characterize Rhizoctonia solani AG-1 IA (Thanatephorus cucumeris) from Texas Fungal Genet Biol 28:148-159[Medline]

Schneider S, Roeslli D, Excoffier L., 2000 ARLEQUIN: a software for population genetics data analysis. Version 2.000 Switzerland: University of Geneva. Free program distributed by the authors over the Internet from the website http://anthro.unige.ch/arlequin

Shew HD, Melton TA., 1995 Target spot of tobacco Plant Dis 79:6-11

Slatkin M., 1994 Linkage disequilibrium in growing and stable populations Genetics 137:331-336[Abstract]

Slatkin M., Excoffier L., 1996 Testing for linkage disequilibrium in genotypic data using EM algorithm Heredity 76:377-383

Stevens Johnk J, Jones RK, Shew HD, Carling DE., 1993 Characterization of populations of Rhizoctonia solani AG-3 from potato and tobacco Phytopathology 83:854-858

Stoddart JA., 1983 A genotypic diversity measure J Hered 74:489-490[Abstract/Free Full Text]

Stoddart JA., Taylor JF., 1988 Genotypic diversity: estimation and prediction in samples Genetics 118:705-711[Abstract/Free Full Text]

Taylor JW, Jacobson DJ, Fisher MC., 1999 The evolution of asexual fungi: reproduction, speciation and classification Annu Rev Phytopathol 37:197-246[Medline]

Vilgalys R, Cubeta MA., 1994 Molecular systematics and population biology of Rhizoctonia Annu Rev Phytopathol 32:135-155

Vilgalys R, Gonzalez D., 1990 Ribosomal DNA restriction fragment length polymorphisms in Rhizoctonia solani, Phytopathology 80:151-158

Weir BS., 1996 Genetic data analysis. 2nd. ed Sunderland, Massachusetts: Sinauer. 445 p

Weir BS., Cockerham CC., 1984 Estimating F-statistics for the analysis of population structure Evolution 38:1358-1370

Xu J, Mitchell T, Vilgalys R., 1999 PCR-restriction fragment length polymorphism (RFLP) analyses reveal both extensive clonality and local genetic differences in Candida albicans Mol Ecol 8:59-73[Medline]

Yeh FC, Yang R-C, Boyle TBJ, Ye ZH, Mao JX., 1997 POPGENE: microsoft Windows®-based freeware for population genetic analysis. Version 1.32 Canada: University of Alberta and Centre for International Forestry Research. Free program distributed by the authors over the Internet from the website at http://www.ualberta.ca/~fyeh/fyeh

Zaykin D, Zhivotovsky L, Weir BS., 1995 Exact tests for association between alleles at arbitrary numbers of loci Genetica 96:169-178[Medline]




This article has been cited by other articles:


Home page
MycologiaHome page
N. D. Charlton and M. A. Cubeta
Transmission of the M2 double-stranded RNA in Rhizoctonia solani anastomosis group 3 (AG-3).
Mycologia, November 1, 2007; 99(6): 859 - 867.
[Abstract] [Full Text] [PDF]


This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Services
Right arrow Similar articles in this journal
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Ceresini, P. C.
Right arrow Articles by Cubeta, M. A.
Right arrow Search for Related Content
PubMed
Right arrow Articles by Ceresini, P. C.
Right arrow Articles by Cubeta, M. A.
Agricola
Right arrow Articles by Ceresini, P. C.
Right arrow Articles by Cubeta, M. A.


HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS