- Research article
- Open Access
Phylogeography of the Solanaceae-infecting Basidiomycota fungus Rhizoctonia solani AG-3 based on sequence analysis of two nuclear DNA loci
© Ceresini et al; licensee BioMed Central Ltd. 2007
- Received: 16 February 2007
- Accepted: 13 September 2007
- Published: 13 September 2007
The soil fungus Rhizoctonia solani anastomosis group 3 (AG-3) is an important pathogen of cultivated plants in the family Solanaceae. Isolates of R. solani AG-3 are taxonomically related based on the composition of cellular fatty acids, phylogenetic analysis of nuclear ribosomal DNA (rDNA) and beta-tubulin gene sequences, and somatic hyphal interactions. Despite the close genetic relationship among isolates of R. solani AG-3, field populations from potato and tobacco exhibit comparative differences in their disease biology, dispersal ecology, host specialization, genetic diversity and population structure. However, little information is available on how field populations of R. solani AG-3 on potato and tobacco are shaped by population genetic processes. In this study, two field populations of R. solani AG-3 from potato in North Carolina (NC) and the Northern USA; and two field populations from tobacco in NC and Southern Brazil were examined using sequence analysis of two cloned regions of nuclear DNA (pP42F and pP89).
Populations of R. solani AG-3 from potato were genetically diverse with a high frequency of heterozygosity, while limited or no genetic diversity was observed within the highly homozygous tobacco populations from NC and Brazil. Except for one isolate (TBR24), all NC and Brazilian isolates from tobacco shared the same alleles. No alleles were shared between potato and tobacco populations of R. solani AG-3, indicating no gene flow between them. To infer historical events that influenced current geographical patterns observed for populations of R. solani AG-3 from potato, we performed an analysis of molecular variance (AMOVA) and a nested clade analysis (NCA). Population differentiation was detected for locus pP89 (Φ ST = 0.257, significant at P < 0.05) but not for locus pP42F (Φ ST = 0.034, not significant). Results based on NCA of the pP89 locus suggest that historical restricted gene flow is a plausible explanation for the geographical association of clades. Coalescent-based simulations of genealogical relationships between populations of R. solani AG-3 from potato and tobacco were used to estimate the amount and directionality of historical migration patterns in time, and the ages of mutations of populations. Low rates of historical movement of genes were observed between the potato and tobacco populations of R. solani AG-3.
The two sisters populations of the basidiomycete fungus R. solani AG-3 from potato and tobacco represent two genetically distinct and historically divergent lineages that have probably evolved within the range of their particular related Solanaceae hosts as sympatric species.
- Haplotype Network
- Anastomosis Group
- North Carolina
- Restricted Gene Flow
- Chromosome Region Maintenance
Rhizoctonia solani Kühn is a species complex composed of genetically distinct groups of fungi in the cantherelloid clade of the phylum Basidiomycota . Current classification within the R. solani species complex is based largely on grouping of isolates into anastomosis groups (AG) based on their hyphal interactions. At least 13 AG have been described within the R. solani species complex including AG-1 to AG-13 [2–4]. Current knowledge suggests that AG and their subgroups represent independent evolutionary units within R. solani [5, 6]. Anastomosis group 3 (AG-3) of Rhizoctonia solani Kühn is associated primarily with diseases of Solanaceous plants. Two well defined phylogenetic (sister) groups of R. solani AG-3 that cause diseases associated with potato and tobacco have been identified recently by analysis of sequence variation in ribosomal DNA (rDNA) [6–8] and beta-tubulin genes : the AG-3 PT and TB. Isolates from the two hosts are closely related, but differ in their dispersal, epidemiology, fatty acid composition, and AFLP patterns [10–16]. Moreover, the host range of isolates of R. solani from potato and tobacco does not overlap, suggesting that they may represent genetically subdivided populations that have evolved a high level of specificity on different Solanaceous plant hosts [3, 8, 16, 17].
The ecology and epidemiology of R. solani AG-3 on potato and tobacco has been extensively studied . In general, isolates of R. solani AG-3 from potato are predominantly asexual and survive as mycelium and sclerotia in soil and on potato seed tubers . When the sexual stage (teleomorph=Thanatephorus cucumeris Frank Donk) is formed during periods of cool and moist weather, meiospores (i.e., basidiospores) have limited dispersal and do not contribute directly to disease epidemics on potato . However, sexual spores produced as a result of heterothallic mating, might constitute an important aspect of fungal life history that contributes to the genetic diversity and structure of field populations of R. solani AG-3 on potato. Ceresini et al  have recently suggested a model of population structure that includes both recombination and clonality for R. solani AG-3 on potato that provides experimental support for this assertion.
In contrast, isolates of R. solani AG-3 from tobacco are predominantly sexual and basidiospores of T. cucumeris serve as infectious propagules that initiate disease on tobacco leaves. However, the mating system of T. cucumeris from tobacco is not known. Likewise, it is not known whether recombination associated with sexual reproduction occurs among field isolates of R. solani AG-3 from tobacco [5, 20].
Despite the evidence for genetic division among populations of R. solani AG-3 from potato and tobacco, many questions remain regarding their population biology and genetics. For example, although recombination had been detected within populations of R. solani AG-3 from potato, the role of asexual and sexual reproduction in determining population structure of R. solani AG-3 from potato is not known. Similarly, questions remain about the extent of gene flow among populations of R. solani AG-3 from tobacco and between populations of AG-3 from potato and tobacco.
For R. solani AG-3 from potato, a recent study supports the concept of migration of the pathogen on potato seed tubers from source populations from Northern US (Maine and Wisconsin) and Canada into NC, which provides evidence (based on classical F statistics) for a low level of genetic differentiation between source and recipient populations . However, the causal role of contemporary gene flow on the observed population structure has not been determined. Traditional F statistics do not use temporal information on allelic variation, which would allow for inferences of evolutionary relationships. For example, there is no information revealing how current geographical patterns of genetic diversity in R. solani AG-3 from potato is influenced by population structure, history, and by a combination of structural and historical events.
In this study, four hypotheses were tested: 1) populations of R. solani AG-3 from potato and tobacco have different levels of genetic diversity; 2)populations of AG-3 from potato and tobacco differ in their extent of recombination; 3) there is no gene flow between populations of R. solani AG-3 from potato and tobacco; and 4) there is no significant association between haplotypes of AG-3 from potato and geographic location.
The tobacco population of R. solani AG-3 is hypothesized to have a higher level of genetic diversity and a predominant recombining structure due to the importance of sexual spores on disease epidemics. In contrast, populations of R. solani AG-3 from potato are hypothesized to have lower level of genetic diversity and a non-recombining structure, probably due to the predominant clonal reproduction system. The complete absence of gene flow between AG-3 from potato and tobacco is hypothesized as the reason for these populations being genetically unconnected. The hypothesis of no geographical association implies that populations of R. solani AG-3 from potato constitutes one single panmictic population with no genetic subdivision.
Therefore our goal was to analyze and further elucidate the phylogeography of populations of R. solani AG-3 from potato and tobacco. To address these phylogeographical hypotheses, we have focused on the observation, description and analysis of the spatial distribution of genotypes of R. solani AG-3 from potato and tobacco and the inference of historical scenarios based on coalescent gene genealogies [22, 23]. In this study, we examined the sequence variation of two cloned nuclear DNA fragments from R. solani AG-3 from potato and tobacco. These regions of the genome were selected based on a priori information from seven PCR-RFLP markers previously used for genotyping R. solani AG-3 from potato and tobacco .
Measures of nucleotide diversity and intragenic recombination
Descriptive analysis of molecular variation within two cloned nuclear DNA fragments from samples of Rhizoctonia solani AG-3 isolates from potato and tobacco
Sample of isolates
Total number of mutations, Eta
Number of polymorphic (segregating) sites
Number of haplotypes, NHap
Haplotype (gene) diversity ± standard deviation
Nucleotide diversity, Pi
Average number of nucleotide differences, k
Sampling variance of k, Vs(k) b
Total variance of k, V(k) b
Total sample: potato and tobacco
0.856 ± 0.048
Eastern NC, Maine and Wisconsin
0.942 ± 0.026
Central NC and Southern Brazil
0.133 ± 0.112
Total Sample: potato and tobacco
0.869 ± 0.048
Eastern NC, Maine and Wisconsin
0.963 ± 0.029
Central NC and Southern Brazil
0.133 ± 0.112
For the potato sample of R. solani AG-3, 14 unique haplotypes were identified for locus pP42F and 20 for locus pP89 (Table 1); whereas, for AG-3 from tobacco only two haplotypes were identified at both pP42F and pP89. Consequently, higher haplotype diversity was observed for AG-3 from potato (0.942 ± 0.026 at pP42F and 0.963 ± 0.029 at pP89). In contrast, for AG-3 from tobacco, a considerably smaller value of haplotype diversity (0.133 ± 0.112) was observed at both pP42F and pP89. The average π value between two random samples was 0.01340 at pP42F and 0.00706 at pP89 for AG-3 from potato, whereas for R. solani AG-3 from tobacco smaller π values ranging from 0.00196 (pP89) to 0.00366 (pP42F) were observed. For R. solani AG-3 from potato, the average number of nucleotide differences (k) was 5.373 [V(k) = 1.947] at pP42F and 7.64 [V(k) = 2.759] at pP89. The average k values for AG-3 from tobacco were smaller and ranged from 1.467 [V(k) = 0.559] at pP42F to 2.133 [V(k) = 0.813] at pP89. Of the total aligned 401 positions at pP42F, 20 (or 4.99%) were polymorphic for AG-3 from potato and 11 (or 2.74%) for AG-3 from tobacco. At pP89, from a total of 1086 positions, 32 (or 2.95%) were variable for R. solani AG-3 from potato and 16 (or 1.47%) for AG-3 from tobacco.
Intraspecific evolution of R. solani AG-3 inferred by analysis of haplotype networks
Analysis of the population structure using AMOVA was conducted independently for both pP42F and pP89 loci. By grouping molecular sequencing haplotypes of R. solani based on their host of origin (either from potato or tobacco), most of the molecular variation was detected between potato and tobacco populations (87.7 to 92.0% of the total variance), with very little variation among geographical populations within each group (0.4 to 1.7%). The within-population variance corresponded to 11.9 and 6.3% of the total variance for pP42F and pP89 respectively. The overall ΦST was 0.88 for pP42F and 0.94 for pP89 (P ≤ 0.001).
Population pairwise ΦST a of Rhizoctonia solani AG-3 PT and TB collected from commercial potato or tobacco fields in Northern US, Eastern and Central North Carolina and Southern Brazil for loci pP42F (above diagonal) and pP89 (below diagonal)
Southern Brazil (PR, SC)
Southern Brazil (PR, SC)
Nested clade analysis
For the NCA of pP89, high sequence variation was observed, and such variation has caused a few cases of ambiguity (probably due to recombination) in the network estimation. As the NCA does not accommodate recombination, haplotypes P0471, P1511 and P1512 were removed from the reticulate network presented in Figure 3-B. A significant geographical association pattern indicating restricted gene flow with isolation by distance was found (Figure 4). The interior haplotypes P0062 and P08310 from Northern America contained in the one-step clade 1-1 showed significantly large (P < 0.05) D c . In addition, these internal haplotypes showed larger D c than the average clade distance for the tip haplotype P0061 that has presumably descended from it [i.e. (Int - Tip) c significantly large at P < 0.05]. At the higher level, within the entire cladogram nesting, the tip clade 3-2 (that included the two-step level clades 2–3 and 2–4, containing only haplotypes from Camden, NC), showed significantly small D c , indicating restricted gene flow as well. In contrast, clade 3-1 showed a significantly large D c , which would indicate a long distance dispersal event, probably from Wisconsin to Maine, which contained the tip haplotypes. Clade 3-1 includes the two-step clades 2-1 (composed by one-step clades 1-1 and 1–2) and 2-2 (which contained the one-step clade 1–3). This 3-step clade contained haplotypes only observed in Northern America.
Estimates from neutrality tests for each locus and populations of Rhizoctonia solani AG-3 from potato and tobacco
Loci and populations
Fu and Li's D
Fu and Li's F*
Potato (Northern US + Eastern NC)
Tobacco (US + Brazil)
Potato (Northern US)
Potato (Eastern NC)
Tobacco (US + Brazil)
Estimates of coalescent parameters from the divergence between potato and tobacco populations of Rhizoctonia solani AG-3
Loci and populations
Effective sample size
Population growth rate a
Theta (4Ne μ) b
Migration rate c 4Nm (+ = receiving population)
Potato (Northern US + Eastern NC)
Tobacco (US + Brazil)
Potato (Northern US)
Potato (Eastern NC)
Tobacco (US + Brazil)
The population parameter Theta was used as a relative measure of effective population size. The estimates of Theta are summarized in Table 4. The pattern observed for the pP42F locus suggests a larger effective population size for the pooled R. solani AG-3 populations (Northern US + Eastern NC) from potato (Theta = 0.01412) when compared to the pooled AG-3 populations (US + Brazil) from tobacco (Theta = 0.00422). However, for the pP89 locus the Theta values ranged from 0.00155 to 0.00206 for the populations of AG-3 from potato from Eastern NC and Northern US, respectively, and 0.00176 for the pooled population samples of R. solani AG-3 from tobacco.
Estimates of directional gene flow between populations of R. solani AG-3 from potato and tobacco were consistent with relatively low rates of migration between them and ranged from 0 to 0.0994. When the pooled populations of R. solani AG-3 from potato were compared with the pooled populations from tobacco for the pP42F locus, the estimates of directional gene flow (since the divergence of R. solani AG-3) were very low, with a historical movement of genes from the pooled potato populations (source) to the pooled tobacco populations (recipient) of R. solani AG-3 (4Nm1→2 = 0.0994). For the pP89 locus, there was no evidence of migration between either of the two populations from potato to the pooled populations from tobacco. However, the migration parameter was estimated to be 0.0417 (4Nm3→2) from the pooled tobacco populations to the Eastern NC potato population. While the highest historical migration contribution was observed from the Eastern NC to the Northern US potato population of AG-3 (migration rate 4Nm2→1 = 5.1115), a much smaller migration value was observed in the opposite direction (4Nm1→2 = 0.1258).
We have studied two closely related sister groups within the basidiomycete fungus R. solani AG-3 on taxonomically related hosts and this system provided us with a unique opportunity to examine aspects of their population genetics in a phylogeographic framework.
During the process of screening of markers for population genetics, four of seven PCR-RFLP markers previously used for genotyping isolates of R. solani AG-3 from potato (pP42, pP46, pP47 and pP89) also amplified genomic DNA from R. solani AG-3 from tobacco. Except for marker pP42F, none of the alleles present in the sample of 20 tobacco isolates of R. solani AG-3 from tobacco were shared with the potato sample . As these markers were shown to be conserved enough to amplify genes across potato and tobacco populations of AG-3, we explored their potential appropriateness for phylogenetic and population genetics analyses within R. solani AG-3. For this study, we selected two (pP42 and pP89) of these four genetic markers. An important outcome from the sequencing analysis of these two markers (pP42 mad pP89) was the increasing power for discriminating alleles within a locus. Genotyping by PCR-RFLPs indicated the occurrence of only two or three alleles (for pP42F and pP89 loci, respectively) within a potato population of R. solani AG-3 from NC . However, genotyping by sequencing analysis indicated that the actual number of alleles present within that population was much higher, with 14 at pP42F and 20 at pP89.
The first objective of this study was to determine the extent of genetic divergence between populations of R. solani AG-3 from potato and tobacco. The potato sample was more genetically diverse showing a higher number of nucleotide differences among haplotypes, which resulted in high haplotype (gene) diversity. In contrast with the sample from potato, limited to no genetic diversity was observed within the tobacco sample from NC and Brazil. The number of nucleotide differences was considerably smaller for tobacco than for potato isolates of R. solani AG-3. Only two haplotypes were found for the tobacco population in Brazil, with a single haplotype predominating in North America. Our results support the first hypothesis that populations of R. solani AG-3 from potato and tobacco have different levels of molecular diversity. Consistent with the heterokaryotic nature of R. solani, every isolate of AG-3 from potato had at least one heterozygous locus. However, all tobacco isolates from both NC and Brazil were homozygous.
The other specific objective of this study was to determine how populations of R. solani AG-3 from potato and tobacco are shaped by recombination. We addressed this objective by analyzing the intragenic recombination at pP42 and pP89 loci. This analysis was done by determining both the recombination parameter R within a locus  and the minimum number of recombination events between two adjacent sites at the same locus .
The tobacco population of R. solani AG-3 was hypothesized to have a predominantly recombining structure due to the importance of sexual spores in disease epidemics. In theory, progeny generated via this type of reproduction would be genetically different than parental individuals as a result of recombination of genetically different nuclei [5, 30–32]. In contrast, the potato population of AG-3 was hypothesized to have a non-recombining structure, probably because of the predominant clonal reproduction system.
Higher values of R were detected for the population of R. solani AG-3 from potato. In addition, a minimum of four recombination events were detected between sites along the pP42 locus and six recombination events at pP89 for the isolates of AG-3 from potato. Putative recombinants were also identified through observations of homoplasy in the networks of both pP42F and pP89 loci (Fig.4. 4A–B). In contrast, no recombination events were detected for AG-3 isolates from tobacco. Our data support the alternative hypotheses of a sexual recombining structure for R. solani AG-3 from potato. This is in agreement with the common occurrence and observation of the sexual stage in potato fields and previous research results that suggested recombination was occurring in population of AG-3 from potato based on analyses of the population structure of the fungus with seven co-dominant PCR-RFLP markers . All seven loci were found in Hardy Weinberg Equilibrium and no linkage disequilibrium was detected for any pair of loci .
In contrast with the potato population, we found no evidence of recombination for the tobacco population. Fungal progeny generated via asexual or homothallic sexual (self-fertile, inbreeding) reproduction would be nearly genetically identical to the parental individuals (and all parts of the genome will have similar evolutionary history) in the absence of outcrossing. No or limited outcrossing would contribute to the association between independent characters (genetic markers), a population with limited genetic diversity and the repeated recovery of similar genotypes (e.g. clonal population structure) and new genotype in the population would arise mainly as a result of mutation .
Our previous observation indicated a clonal population structure for AG-3 from tobacco in NC based on the one-to-one (strict) association between two independent criteria (somatic incompatibility and AFLP markers) . However, our current data is not sufficient for supporting the alternative hypothesis of clonality for AG-3 from tobacco because only two haplotypes were detected for each locus in the tobacco populations. The limited amount of genetic variation in this sample would not be sufficient to accurately assess the relative contribution of recombination in the pathogen population.
The third hypothesis tested in this study was that no gene flow occurs between populations of R. solani AG-3 from potato and tobacco, indicating that they constitute genetically distinct groups. The objective was to describe this pair of populations of R. solani AG-3 by the proportion of loci that show unilateral or reciprocal fixation, and by standard estimates of statistics based on gene flow.
In recently evolved species pairs it is expected to identify variable nucleotide positions where polymorphisms are shared by both sibling species, other nucleotide positions where only one species is polymorphic and the other is fixed for one allele, and where both species have fixed loci [11, 33]. An analysis of the DNA polymorphisms at each locus, suggested that most of the nucleotide diversity occurred between samples of potato and tobacco R. solani AG-3. As genetic isolation is considered to precede the loss of shared polymorphism, discovering a single locus or multiple loci that show fixation in one or the other of the phylogenetic species would be evidence of genetic isolation [11, 33]. No alleles were shared between the potato and tobacco populations of R. solani AG-3 used in this study. In addition, analysis of population structure indicated no gene flow between populations of R. solani AG-3 from potato and tobacco. Most of the molecular variation and differentiation was detected between potato and tobacco populations of R. solani AG-3. In contrast, a high level of gene flow was detected for tobacco group of R. solani AG-3. With only one exception (isolate TBR24), all of the isolates from NC and the Brazilian population of AG-3 from tobacco shared the same alleles at both pP42 and pP89 loci.
From the Bayesian phylogenetics analysis of cloned nuclear DNA fragments pP42F and pP89, two major clades in R. solani AG-3 were inferred, splitting haplotypes from potato and tobacco into two groups. In addition, the reconstruction of the haplotype network of both pP42F and pP89 loci by statistical parsimony indicated no mutational connection linking any of the haplotypes of R. solani AG-3 from potato and tobacco in this network. These results support the hypothesis of divergent evolution of genes from populations of R. solani AG-3 from potato and tobacco.
The fourth hypothesis tested was that there is no significant association between haplotypes of AG-3 from potato and geographic origin, which implies that there is a single panmictic population with no genetic subdivision.
Based on F statistics, an important observation about gene flow was inferred from the within-group analysis of population subdivision of the sample of R. solani AG-3 from potato. While gene flow (Φ ST = 0.03414, not significant) was detected between populations of R. solani AG-3 from potato at the pP42F locus, population differentiation was detected for the locus pP89.
To infer historical events influencing the current geographical patterns observed on populations of the pathogen, we performed NCA analysis on a reticulated phylogeny of AG-3 haplotypes from potato. The main benefit of using haplotype tree information is its qualitative advantage of discriminating among various biological explanations for any detected geographical association. By using a haplotype network analysis it was possible to examine spatial and temporal patterns of genetic variation in populations of R. solani AG-3 from potato whereas the F-statistics analysis could only examine the current spatial pattern.
Based on NCA, the null hypothesis of no association between haplotype variation and geography was rejected for pP89 locus. For this locus, the first prediction for the population structure hypothesis was that restricted gene flow caused the geographical associations. Within clade 1-1 from Northern America, clade distance equal to zero was associated with the tip pP89 haplotype P0061 whereas significantly large clade distance was associated with the interior haplotypes P0062 and P08310. In addition, the average interior clade distance minus the average tip clade distance was significantly large. These observations fulfilled the prediction that under restricted gene flow younger or tip clades should be less widespread relatively to older clades interior to them [34, 35]. This is because restricted gene flow implies only limited movement by individuals during any giving generation as it takes time for a new haplotype to spread geographically. Another aspect of restricted gene flow is that the new haplotypes initially reside within the range of its ancestor for many generations, which is predicted by the isolation by distance model. In contrast, the ancestor would have a wider geographic distribution .
At the higher clade level the tip clade 3-2 (that included the two-step level clades 2–3 and 2–4, containing only haplotypes from Camden, NC), showed significantly small D n , indicating restricted gene flow. In contrast, clade 3-1 showed a significantly large D n , which would indicate a long distance dispersal event, probably from Wisconsin to Maine, which contained the tip haplotypes. However, conclusions concerning a long distance dispersal event should be interpreted cautiously since our sampling procedure was geographically limited. Therefore, it may not be possible to discriminate between isolation by distance and long distance gene flow models .
The patterns of spatial structure and differentiation between potato populations of R. solani AG-3 revealed by the two distinct markers were conflicting and deserve further discussion. Discordant patterns of divergence in molecular markers and loci are commonly interpreted as differences in mutation rate, neutrality or linkage disequilibrium with other loci that are subject to selection, and the differential effect of natural selection on those markers and loci [37–40]. Although we do not have ancillary information on whether the genetic loci examined in this study have distinct mutation rates, our data suggest a model of neutral evolution for both loci (Table 3). A plausible explanation for the conflicting patterns of spatial structure is due to the distinct resolution capabilities associated with difference in per-locus sequence length. For pP42 we detected 20 polymorphic sites in 401 bp and for pP89 we detected 32 segregating sites in 1090 bp. The paucity of informative sites on pP42F resulting from the analyses of less sequence length might have introduced a downward bias in the estimate of population differentiation. Support for this assertion is based on simulation data on how many loci are needed and how much sequence is required from each locus to provide better estimates of population genetic parameters .
In this study, we also conducted coalescent-based simulations of genealogical relationships between populations of R. solani AG-3 from potato and tobacco to estimate the amount and directionality of historical migration patterns in time, the ages of mutations and of populations (Table 4; Figure 5). This is the first time that a genealogical approach (based on the coalescent model) was used to study the population genetics and evolution of a Rhizoctonia pathosystem.
Tests of non-neutral evolution could not be rejected for the two loci, considering only the populations from tobacco. We believe that the observed patterns of gene diversity is most likely due to changes in population size and not to a intrinsic feature of the loci. Given the contemporary Theta values  it was possible to extrapolate the relative increase or decrease in population size since the time of divergence. Growth rate estimates  were consistent with a significant decrease in size for the tobacco populations of R. solani AG-3. We hypothesize that this decrease in size reflects a split where only a small fraction of the ancestral pathogen population founded a new population during the domestication of tobacco from the wild habitats into agricultural fields. In contrast, results from Lamarc analyses  were consistent with expansion of the potato populations of R. solani AG-3. The increase in R. solani AG-3 effective population size on potato was most likely caused by a corresponding increase in host population size following the expansion of agriculture.
Although characterized as being younger than the tobacco population of AG-3, the potato population has experienced a very recent diversification in the US. In contrast, the two tobacco populations sampled from Southern Brazil and Central NC (US), had a small effective population size. Older ancestor populations with higher haplotype diversity of R. solani (and consequently higher effective population sizes) might exist in the center of origin of tobacco and other related Solanaceae.
Our findings corroborate previous information that populations of R. solani AG-3 are subdivided by host specialization, constituting two distinct phylotypes [3, 7, 8, 11, 16, 17]. R. solani AG-3 represents at least two phylogenetic species that display high phenotypic differences in disease biology and ecology [3, 16, 17].
Low rates of historical movement of genes were observed between the potato and tobacco populations of AG-3. These two populations of R. solani AG-3 represented two distinctly and historically divergent lineages, which have probably evolved within the range of their particular hosts as sympatric species. Reinforcement is considered as an important mechanism of speciation between sympatrically divergent populations . Reinforcement is the active mechanism of reproductive isolation due to reduced fitness of the interspecific hybrids, thereby favoring intraspecific mating and resulting in genetic differentiation between species . Reproductive success between R. solani AG-3 from potato and tobacco was not tested in our study, but this testing could provide knowledge about the mechanisms involved in sympatric speciation of R. solani AG-3.
The contrasting ecological context of these two host specific populations makes our study a good example of ecological speciation between lineages of plant pathogenic fungal populations. Host range and disease causing ability are clearly important factors in fungal speciation. Consistent with our study, evidence for lineage divergence strongly associated with closely-related host species has been previously reported for other fungi, such as Magnaporthe grisea and Mycosphaerella graminicola [43, 45, 46]. Several studies have also reported extensive genetic differentiation between divergent fungal species [47–51].
The results presented in this paper provide important insight into our understanding of the phylogeography of populations of R. solani AG-3 from potato and tobacco. The next step of our research will involve a more in depth study to reveal the origins and patterns of domestication of these two phylogenetically distinct populations of R. solani AG-3. Our experimental approach will involve a global sampling strategy for isolation of R. solani AG-3 from Solanum tuberosum (potato) and Nicotianna tabacum (tobacco) and including additional host species of Solanum and Nicotiana.
A global phylogeographical framework is critical to unravel the answers on large scale-hypotheses on the origins and evolution of R. solani AG-3 as a pathogen of domesticated crops within the Solanaceae. A sampling initiative towards the understanding of the phylogeography of populations from potato and tobacco should include the main potato and tobacco growing areas in the world, besides the centers of origin in the Andes mountain valleys of Bolivia, Peru, Columbia and Ecuador, in South America (for hosts on both Nicotiana and Solanum genera) . It also should include Africa and Australia for the Nicotiana, which is part of a truly widespread tribe within the Solanaceae with a Gondwanan distribution (the southern supercontinent, which included most of the landmasses in today's southern hemisphere) .
The two sister populations of the basidiomycete fungus R. solani AG-3 from potato and tobacco represent two distinct and historically divergent lineages that have probably evolved within the range of their particular related Solanaceae hosts as sympatric species.
Isolates and populations
Description and origin of pP42F and pP89 haplotypes of Rhizoctonia solani AG-3 from potato and tobacco
Origin (region and State)
Haplotypes designation for phylogenetic analysis
Haplotype designation for coalescent analysis*
Sequence code at GenBank® NCBI
Haplotypes designation for phylogenetic analysis
Haplotype designation for coalescent analysis*
Sequence code at GenBan® NCBI
Parana and Santa Catarina States
DNA technique for genotyping
Primers for PCR amplification and sequencing reactions of two nuclear DNA fragments from Rhizoctonia solani AG-3 PT and TB
Sequence (5' to 3')
PCR reactions were conduced in a Model 9600 DNA ThermoCycler (Perkin-Elmer) employing standard conditions described previously for R. solani  and using the Expand High Fidelity PCR System (Roche Molecular Biochemicals), which contains a thermostable, exonuclease-free, N-terminal deletion mutant of Taq DNA polymerase and a proofreading polymerase (Pfu) exhibiting 3'-exonuclease activity. After purification by ultrafiltration with QIAquick PCR Purification Kit (QIAGEN®), the PCR products were sequenced in both directions by the dideoxy chain terminating method using the Big Dye Terminator Cycle Sequencing Ready Reaction Kit (Applied Biosystems, Inc.) and analyzed on an ABI 3700 automatic sequencer (Perkin Elmer Corp., Norwalk, CT). Sequence chromatograms were compiled using Sequencher software (vers. 2.0, GeneCodes Corp.).
To separate different pP42F or pP89 alleles within heterogenous PCR reactions, amplicons were cloned into the PCR2.1-TOPO® vector (Invitrogen, San Diego). Plasmids from selected recombinant One Shot® DH5a™-T1R Escherichia coli (Invitrogen, San Diego) were extracted from each sample and purified using QIAprep Spin Miniprep Kit (QIAGEN®). Selected primers were used to reamplify and sequence the clones from PCR2.1-TOPO® (Table 6). Sequences were aligned using the software Sequencher (vers. 2.0, GeneCodes Corp.).
Measures of nucleotide diversity and intragenic recombination
Nucleotide diversity, or the average number of differences per site between two homologous sequences (π), was calculated using the program DNASP  according to the equation 10.5 of Nei . Diversities were calculated for the population and separately for the two major clades. The spectrum of polymorphism along each locus was measured and represented as π along a sliding window of 25 positions. The minimum number of recombination events between two adjacent sites (Rm) was calculated according to Hudson and Kaplan [24, 25]. Recombination parameters were estimated according to Hudson .
The phylogenetic analyses were conducted using maximum likelihood (ML) and Bayesian methods. For ML, the software PAUP* 4.0b10 was used . The Bayesian analysis was conducted using the Monte Carlo Markov Chain method coupled with the Metropolis algorithm (MCMCMC), implemented by the software MrBayes v.2.10 .
For both ML and MCMCMC analyses, MODELTEST 3.7  was used to determine the DNA bases substitution model more fit to the data. A hierarchical likelihood ratio test implemented by MODELTEST selected the following models for the markers: a) pP42F: K80 (K2P)+G model (number of substitution types = 2; transition/transversion ratio = 3.1683 (kappa = 6.3366); assumed proportion of invariable sites = none; assumed nucleotides frequencies = equal; gamma distribution of rates at variable sites, shape parameter (alpha) = 0.1147; under this evolution model for base substitution, the ML value for the pP42F tree was -lnL = 886.9029); b) pP89: K80 (K2P)+G+I model (number of substitution types = 2; transition/transversion ratio = 15.4892 (kappa = 30.9784); assumed proportion of invariable sites = 0.8057; assumed nucleotides frequencies = equal; gamma distribution of rates at variable sites, shape parameter (alpha) = 0.5569; under this evolution model for base substitution, the ML value for the pP89 tree was -lnL = 2271.1220).
The phylogenetic analysis by MCMCMC was performed using, as a priori hypothesis, the model of DNA bases substitution described above for the ML analysis. We have searched for 3,000,000 generations using four chains of search, eliminating the first 500 trees, and storing one tree every 500. Posterior probability supporting values (PP) were generated for the branches of the consensus tree considering the rule of a minimum of 50% of the 11002 trees showing a particular partition indicated by the branches.
Intraspecific evolution of R. solani AG-3 inferred by analysis of haplotype networks
For inferring intraspecific evolution, networks of haplotypes sampled from populations of R. solani AG-3 were built using the algorithm recommended by Posada & Crandall . This method starts by estimating the maximum number or differences between haplotypes as a result of single substitutions (i.e., those not resulting from multiple substitutions in a single site) with the statistical significance level of 95% . This is called the limit of parsimony, or the limit of parsimony connection. Haplotypes differing by only one change are connected, and then those differing by two, three and so on until all the haplotypes are included in a single network, or until the limit of parsimony connection is reached. This method also allows for the identification of probable recombinants through observations of spatial distribution of the sequence of homoplasies defined by the network . The estimation of the phylogeny from DNA sequences was implemented by the computer program TCS .
Between-population analyses were conducted to test for geographic structure in populations of R. solani AG-3. Pairwise Φ statistics comparisons (F statistics equivalents ) were used to describe the proportion of the total genetic variance at each locus due to differences among populations. The sequences were grouped by sample origin into four populations of R. solani AG-3: potato tubers from the Northern US; potato stems from Eastern NC; tobacco leaves from Central NC; and tobacco leaves from Southern Brazil. FST was calculated via analysis of molecular variance (AMOVA) using the software package ARLEQUIN [63, 64]. AMOVA estimates variance components considering the number of differences between molecular genotypes [63, 65, 66]. FST ranges from 0.0, in which all populations appears homogenous, to 1.0, in which all the variation is among populations.
Nested clade analysis
The nested clade analysis (NCA) on haplotypes of R. solani AG-3 from potato was performed with the objective of discriminating between phylogenetic patterns that result from historical events (for example, past fragmentation, expansion, colonization) and from those due to recurrent gene flow [34, 35]. This analysis utilizes both a phylogeny and the geographical distribution of haplotypes or haplotype clades. The NCA was completed using the computer program GeoDis v.2.0  considering the haplotype network (phylogeny) obtained by statistical parsimony analysis as described above. This analysis divides the phylogeny into nested clades of n steps, where n is equal to the number of mutations connecting the haplotypes of a clade. For each clade there is a geographic point around which all the individuals are centered. Two quantitative measures of how within-clade haplotypes are geographically dispersed were used: the clade distance = D c (X), which describes the average distance of each haplotype within clade X from the geographic center of clade X; and the nested clade distance = D n (X), which similarly describes the average distance of each haplotype within clade X from the geographic center of the next more inclusive clade within which the clade X is nested . The clade distance D c (X) measures the geographical range of a particular clade X (average distance that a haplotype lies from the geographical center of all haplotypes from the same clade). The nested clade distance D n (X) measures how a particular clade X is geographically distributed relatively to its closest evolutionary sister clades (i.e., clades in the same higher level nesting category). Accordingly, within each nested category and for both types of distances (D c and D n ), we calculated and determined the significance of the average interior distance minus the average tip distance [(Int - Tip) c and (Int - Tip) n ]. Contrasts of interior versus tip clades within a nested clade strongly indicate contrasts of older versus younger clades. Haplotypes were assigned to the following populations: Maine (46.51562 N 68.36468 W), Wisconsin (43.07980 N 89.38751 W), and the North Carolina counties of Washington (35.83681 N 76.56470 W), Tyrrell (35.87039 N 76.16513 W), Hyde (35.40617 N 76.15316 W) and Camden (36.33743 N 76.16263 W). Clade distances were calculated using a distance matrix among populations based on geographic coordinates. A test of 1000 permutations of population assignments to haplotypes was used to determine the significance of larger or smaller values of D c , D n , (Int - Tip) c and (Int - Tip) n than would be expected for null hypothesis of random geographical association of haplotypes .
SNAP Workbench Java program package was used to analyze gene genealogies and population parameters . SNAP Map was used to collapse sequences into haplotypes, removing indels and infinite sites violations among the mutations . Coalescent methods make strict assumptions, such as neutrality and the lack of recombination (identified based on the detection of homoplasious sites), which must be verified previously. To test for departures from neutrality, Tajima's D  and Fu and Li's D and F* [59, 71] test statistics were calculated using the statistical tests of neutrality of mutations against excess of recent mutations (rare alleles), written by Yun-Xin Fu . We hypothesized that any deviation from neutrality was due to change in effective population size. This hypothesis was tested using the population growth estimates inferred by the program Lamarc 2.0 [28, 29]. The Lamarc algorithm permits intragenic recombination and therefore allowed us to use all polymorphic sites to assess this population parameter. To estimate the population growth parameter, we used 10 initial chains with 2,000 genealogies sampled and 2 final chains with 10,000 genealogies sampled. The population growth rates were inferred using Bayesian analyses, assessing the 95% credibility intervals from the Bayesian search.
Based on the evidence for no geographic population structure, we pooled data for the two tobacco populations of R. solani AG-3 (US + Brazil) for both loci. Data for the potato populations (Northern US + Eastern NC) were also pooled for the pP42F marker, but kept separate for the pP89. Coalescent parameters were estimated for each locus separately, due to the lack of information on how the alleles at each locus are linked in the haploid phase of the fungus. To determine the order of coalescent events for haplotypes in time it is necessary to determine the amount of migration that occurred between populations. Haplotypes from populations linked by migration are postulated to coalesce before haplotypes from unlinked populations . Using the program Migrate , incorporated in the SNAP Workbench program package, migration matrices were constructed for each locus, indicating the number of migrants exchanged between populations. These backward migration matrices were used for ancestral inference in Genetree version 9.0 [74, 75], also incorporated in the SNAP Workbench. Subsequently, we reconstructed the genealogy with the highest root probability, the ages of mutation and the TMRCA (time for the most recent common ancestor) of the samples using 100,000 coalescent simulations with five runs with distinct starting random number seeds, considering population subdivision and distinct population sizes. The program allows for the estimation of the ancestral history of each haplotype. It also shows the distribution of mutations on a coalescent scale, allowing for comparison of the divergence of haplotypes between and within each population.
The first author wish to acknowledge the North Carolina State University, Department of Plant Pathology, for supporting his post-doctoral academic training by a research assistantship grant.
- Binder M, Hibbett D, Larsson K, Larsson E, Langer E, Langer G: The phylogenetic distribution of resupinate forms across the major clades of mushroom-forming fungi (Homobasidiomycetes). Systematics and Biodiversity. 2005, 3 (2): 1-45. 10.1017/S1477200005001623.View ArticleGoogle Scholar
- Carling DE: Grouping in Rhizoctonia solani by hyphal anastomosis reaction. Rhizoctonia Species: Taxonomy, Molecular Biology, Ecology, Pathology and Disease Control. Edited by: Sneh B, Jabaji-Hare S, Neate SM, Dijst G. 1996, Dordrecht: Kluwer Academic Publishers, 37-47.View ArticleGoogle Scholar
- Ogoshi A: Ecology and pathogenicity of anastomosis and intraspecific groups of Rhizoctonia solani Kühn. Annual Review of Phytopathology. 1987, 25: 125-143.View ArticleGoogle Scholar
- Sneh B, Jabaji-Hare S, Neate S, Dijst G: Rhizoctonia species: Taxonomy, Molecular Biology, Ecology, Pathology and Disease Control. 1996, Dordrecht: Kluwer Academic PublishersView ArticleGoogle Scholar
- Cubeta MA, Vilgalys R: Population biology of the Rhizoctonia solani complex. Phytopathology. 1997, 87: 480-484. 10.1094/PHYTO.19184.108.40.2060.View ArticlePubMedGoogle Scholar
- Gonzalez D, Carling DE, Kuninaga S, Vilgalys R, Cubeta MA: Ribosomal DNA systematics of Ceratobasidium and Thanatephorus with Rhizoctonia anamorphs. Mycologia. 2001, 93: 1138-1150. 10.2307/3761674.View ArticleGoogle Scholar
- Kuninaga S, Natsuaki T, Takeuchi T, Yokosawa R: Sequence variation of the rDNA ITS regions within and between anastomosis groups in Rhizoctonia solani. Current Genetics. 1997, 32 (3): 237-243. 10.1007/s002940050272.View ArticlePubMedGoogle Scholar
- Kuninaga S, Carling DE, Takeuchi T, Yokosawa R: Comparison of rDNA-ITS Sequences between Potato and Tobacco Strains in Rhizoctonia solani AG-3. Journal of General Plant Pathology. 2000, 66 (1): 2-11. 10.1007/PL00012917.View ArticleGoogle Scholar
- Gonzalez D, Cubeta MA, Vilgalys R: Phylogenetic utility of indels within ribosomal DNA and [beta]-tubulin sequences from fungi in the Rhizoctonia solani species complex. Molecular Phylogenetics and Evolution. 2006, 40 (2): 459-470. 10.1016/j.ympev.2006.03.022.View ArticlePubMedGoogle Scholar
- Balali GR, Whisson DL, Scott ES, Neate SM: DNA fingerprinting probe specific to isolates of Rhizoctonia solani AG-3. Mycological Reseach. 1996, 100 (4): 467-470.View ArticleGoogle Scholar
- Ceresini PC, Shew HD, Vilgalys RJ, Cubeta MA: Genetic diversity of Rhizoctonia solani AG-3 from potato and tobacco in North Carolina. Mycologia. 2002, 94 (3): 437-449. 10.2307/3761778.View ArticlePubMedGoogle Scholar
- Ceresini PC, Shew HD, Vilgalys RJ, Rosewich UL, Cubeta MA: Genetic structure of populations of Rhizoctonia solani AG-3 on potato in eastern North Carolina. Mycologia. 2002, 94 (3): 450-460. 10.2307/3761779.View ArticlePubMedGoogle Scholar
- Laroche JP, Jabaji-Hare SH, Charest PM: Differentiation of two anastomosis groups of Rhizoctonia solani by isozyme analysis. Phytopathology. 1992, 82: 1387-1392. 10.1094/Phyto-82-1387.View ArticleGoogle Scholar
- Liu ZL, Domier LL, Sinclair JB: ISG-specific ribosomal DNA polymorphism of the Rhizoctonia solani species complex. Mycologia. 1993, 85: 795-800. 10.2307/3760611.View ArticleGoogle Scholar
- Liu ZL, Sinclair JB: Genetic diversity of Rhizoctonia solani anastomosis group 2. Phytopathology. 1992, 82: 778-787. 10.1094/Phyto-82-778.View ArticleGoogle Scholar
- Stevens Johnk J, Jones RK, Shew HD, Carling DE: Characterization of populations of Rhizoctonia solani AG-3 from potato and tobacco. Phytopathology. 1993, 83: 854-858. 10.1094/Phyto-83-854.View ArticleGoogle Scholar
- Shew HD, Melton TA: Target spot of tobacco. Plant Disease. 1995, 79: 6-11.View ArticleGoogle Scholar
- Cubeta MA, Vilgalys R: Rhizoctonia. Encyclopedia of Microbiology. Edited by: Lederberg J. 2000, San Diego: Academic Press, 4: 109-116.Google Scholar
- Carling DE, Leiner RH, Westphale PC: Symptoms, signs, and yield reduction associated with rhizoctonia disease of potato induced by tuberborne inoculum of Rhizoctonia solani. American Potato Journal. 1989, 66: 693-702.View ArticleGoogle Scholar
- Vilgalys R, Cubeta MA: Molecular systematics and population biology of Rhizoctonia. Annual Review of Phytopathology. 1994, 32: 135-155. 10.1146/annurev.py.32.090194.001031.View ArticleGoogle Scholar
- Ceresini PC, Shew HD, Vilgalys R, Gale LR, Cubeta MA: Detecting migrants in populations of Rhizoctonia solani anastomosis group 3 from potato in North Carolina using multilocus genotype probabilities. Phytopathology. 2003, 93: 610-615. 10.1094/PHYTO.2003.93.5.610.View ArticlePubMedGoogle Scholar
- Knowles LL: The burgeoning field of statistical phylogeography. Journal of Evolutionary Biology. 2004, 17 (1): 1-10. 10.1046/j.1420-9101.2003.00644.x.View ArticlePubMedGoogle Scholar
- Carbone I, Kohn LM: Inferring process from pattern in fungal population genetics. Applied Mycology and Biotechnology. Edited by: Arora DK, Khachatourians G. 2004, Amsterdam: Elsevier Science BV, 4: 29-58.Google Scholar
- Hudson RR, Kaplan NL: Statistical properties of the number of recombination events in the history of a sample of DNA sequences. Genetics. 1985, 111: 147-164.PubMed CentralPubMedGoogle Scholar
- Hudson RR: Estimating the recombination parameter of a finite population model without selection. Genetical Research. 1987, 50 (3): 245-250.View ArticlePubMedGoogle Scholar
- Posada D, Crandall KA: Intraspecific gene genealogies: trees grafting into networks. Trends in Ecology & Evolution. 2001, 16 (1): 37-45. 10.1016/S0169-5347(00)02026-7.View ArticleGoogle Scholar
- Fu YX: Statistical tests of neutrality of mutations against population growth, hitchhiking and background selection. Genetics. 1997, 147: 915-925.PubMed CentralPubMedGoogle Scholar
- Kuhner MK: LAMARC 2.0: maximum likelihood and Bayesian estimation of population parameters. Bioinformatics. 2006, 22 (6): 768-770. 10.1093/bioinformatics/btk051.View ArticlePubMedGoogle Scholar
- Kuhner MK, Smith LP: Comparing likelihood and bayesian coalescent estimation of population parameters. Genetics. 2007, 175 (1): 155-165. 10.1534/genetics.106.056457.PubMed CentralView ArticlePubMedGoogle Scholar
- Anderson JB, Kohn LM: Clonality in soilborne, plant-pathogenic fungi. Annual Review of Phytopathology. 1995, 33 (1): 369-391. 10.1146/annurev.py.33.090195.002101.View ArticlePubMedGoogle Scholar
- Milgroom MG: Recombination and the multilocus structure of fungal populations. Annual Review of Phytopathology. 1996, 34: 457-477. 10.1146/annurev.phyto.34.1.457.View ArticlePubMedGoogle Scholar
- Taylor JW, Jacobson DJ, Fisher MC: The evolution of asexual fungi: Reproduction, Speciation and Classification. Annual Review of Phytopathology. 1999, 37 (1): 197-246. 10.1146/annurev.phyto.37.1.197.View ArticlePubMedGoogle Scholar
- Taylor JW, Jacobson DJ, Kroken S, Kasuga T, Geiser DM, Hibbett DS, Fisher MC: Phylogenetic species recognition and species concepts in fungi. Fungal Genetics and Biology. 2000, 31 (1): 21-32. 10.1006/fgbi.2000.1228.View ArticlePubMedGoogle Scholar
- Templeton AR: Nested clade analyses of phylogeographic data: testing hypotheses about gene flow and population history. Molecular Ecology. 1998, 7 (4): 381-397. 10.1046/j.1365-294x.1998.00308.x.View ArticlePubMedGoogle Scholar
- Templeton AR, Routman E, Phillips CA: Separating population structure from population history: A cladistic analysis of the geographical distribution of mitochondrial DNA haplotypes in the tiger salamander, Ambystoma tigrinum. Genetics. 1995, 140 (2): 767-782.PubMed CentralPubMedGoogle Scholar
- Crandall KA, Templeton AR: Empirical tests of some predictions from coalescent theory with applications to intraspecific phylogeny reconstruction. Genetics. 1993, 134 (3): 959-969.PubMed CentralPubMedGoogle Scholar
- Karl SA, Avise JC: Balancing selection at allozyme loci in oysters – Implications from nuclear RFLPs. Science. 1992, 256 (5053): 100-102. 10.1126/science.1348870.View ArticlePubMedGoogle Scholar
- Latta RG, Mitton JB: A comparison of population differentiation across four classes of gene marker in limber pine (Pinus flexilis James). Genetics. 1997, 146 (3): 1153-1163.PubMed CentralPubMedGoogle Scholar
- McDonald JH: Detecting natural selection by comparing geographic variation in protein and DNA polymorphisms. 1994, New York: Chapman & HallView ArticleGoogle Scholar
- Volis S, Yakubov B, Shulgina I, Ward D, Mendlinger S: Distinguishing adaptive from nonadaptive genetic differentiation: comparison of QST and FST at two spatial scales. Heredity. 2005, 95 (6): 466-475. 10.1038/sj.hdy.6800745.View ArticlePubMedGoogle Scholar
- Carling MD, Brumfield RT: Gene sampling strategies for multi-locus population estimates of genetic diversity (θ). PLoS ONE. 2007, 2 (1): e160-10.1371/journal.pone.0000160.PubMed CentralView ArticlePubMedGoogle Scholar
- Beerli P, Felsenstein J: Maximum-likelihood estimation of migration rates and effective population numbers in two populations using a coalescent approach. Genetics. 1999, 152 (2): 763-773.PubMed CentralPubMedGoogle Scholar
- Kohn LM: Mechanisms of fungal speciation. Annual Review of Phytopathology. 2005, 43 (1): 279-308. 10.1146/annurev.phyto.43.040204.135958.View ArticlePubMedGoogle Scholar
- Dobzhansky T: Genetics and the origin of species. 1937, New York: Columbia University PressGoogle Scholar
- Couch BC, Fudal I, Lebrun M-H, Tharreau D, Valent B, van Kim P, Notteghem J-L, Kohn LM: Origins of host-specific populations of the blast pathogen Magnaporthe oryzae in crop domestication with subsequent expansion of pandemic clones on rice and weeds of rice. Genetics. 2005, 170 (2): 613-630. 10.1534/genetics.105.041780.PubMed CentralView ArticlePubMedGoogle Scholar
- Stukenbrock EH, Banke S, Javan-Nikkhah M, McDonald BA: Origin and domestication of the fungal wheat pathogen Mycosphaerella graminicola via sympatric speciation. Molecular Biology and Evolution. 2007, 24 (2): 398-411. 10.1093/molbev/msl169.View ArticlePubMedGoogle Scholar
- Dettman JR, Jacobson DJ, Taylor JW: A multilocus genealogical approach to phylogenetic species recognition in the model eukaryote Neurospora. Evolution. 2003, 57 (12): 2703-2720.View ArticlePubMedGoogle Scholar
- Dettman JR, Jacobson DJ, Turner E, Pringle A, Taylor JW: Reproductive isolation and phylogenetic divergence in Neurospora: comparing methods of species recognition in a model eukaryote. Evolution. 2003, 57 (12): 2721-2741.View ArticlePubMedGoogle Scholar
- Vilgalys RJ, Johnson JL: Extensive genetic divergence associated with speciation in filamentous fungi. PNAS. 1987, 84 (8): 2355-2358. 10.1073/pnas.84.8.2355.PubMed CentralView ArticlePubMedGoogle Scholar
- Horgen PA, Arthur R, Davy O, Moum A, Herr F, Straus N, Anderson J: The nucleotide sequence homologies of unique DNA's of some cultivated and wild mushrooms. Canadian Journal of Microbiology. 1984, 30: 587-593.View ArticleGoogle Scholar
- Bicknell JN, Douglas HC: Nucleic acid homologies among species of Saccharomyces. Journal of Bacteriology. 1970, 101 (2): 505-512.PubMed CentralPubMedGoogle Scholar
- Symon DE: Gondwanan elements of the Solanaceae. Solanaceae III: Taxonomy – Chemistry – Evolution. Edited by: Hawkes JG, Lester RN, Nee M, Estrada N. 1991, Richmond: Royal Botanical Gardens Kew, 139-150.Google Scholar
- Altschul S, Madden T, Schaffer A, Zhang J, Zhang Z, Miller W, Lipman D: Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Research. 1997, 25 (17): 3389-3402. 10.1093/nar/25.17.3389.PubMed CentralView ArticlePubMedGoogle Scholar
- Rozas J, Rozas R: DnaSP version 3: an integrated program for molecular population genetics and molecular evolution analysis. Bioinformatics. 1999, 15 (2): 174-175. 10.1093/bioinformatics/15.2.174.View ArticlePubMedGoogle Scholar
- Nei M: Molecular Evolutionary Genetics. 1987, New York: Columbia Univ. PressGoogle Scholar
- Swofford DL: PAUP* 4.0: Phylogenetic Analysis Using Parsimony, 4.0. 2002, Champaign: Illinois Natural History SurveyGoogle Scholar
- Huelsenbeck JP, Ronquist F: MRBAYES: Bayesian inference of phylogenetic trees. Bioinformatics. 2001, 17 (8): 754-755. 10.1093/bioinformatics/17.8.754.View ArticlePubMedGoogle Scholar
- Posada D, Crandall KA: MODELTEST: testing the model of DNA substitution. Bioinformatics. 1998, 14 (9): 817-818. 10.1093/bioinformatics/14.9.817.View ArticlePubMedGoogle Scholar
- Templeton AR, Crandall KA, Sing CF: A Cladistic analysis of phenotypic associations with haplotypes inferred from restriction endonuclease mapping and DNA sequence data. III. Cladogram estimation. Genetics. 1992, 132 (2): 619-633.PubMed CentralPubMedGoogle Scholar
- Crandall KA, Templeton AR: Statistical methods for detecting recombination. The Evolution of HIV. Edited by: Crandall KA. 1999, Baltimore: The Johns Hopkins Univ. Press, 153-176.Google Scholar
- Clement M, Posada D, Crandall KA: TCS: a computer program to estimate gene genealogies. Molecular Ecology. 2000, 9 (10): 1657-1659. 10.1046/j.1365-294x.2000.01020.x.View ArticlePubMedGoogle Scholar
- Wright S: The genetical structure of populations. Annals of Eugenics. 1951, 15: 323-354.View ArticlePubMedGoogle Scholar
- Excoffier L, Smouse PE, Quattro JM: Analysis of molecular variance inferred from metric distances among DNA haplotypes: Application to human mitochondrial DNA restriction data. Genetics. 1992, 131 (2): 479-491.PubMed CentralPubMedGoogle Scholar
- Schneider S, Roeslli D, Excoffier L: ARLEQUIN: A software for population genetics data analysis, Version 2.000. 2000, Geneva: University of GenevaGoogle Scholar
- Weir BS, Cockerman CC: Estimating F-statistics for the analysis of population structure. Evolution. 1984, 38: 1358-1370. 10.2307/2408641.View ArticleGoogle Scholar
- Weir BS: Genetic Data Analysis. 1996, Sunderland: Sinauer, 2Google Scholar
- Posada D, Crandall KA, Templeton AR: GeoDis: A program for the cladistic nested analysis of the geographical distribution of genetic haplotypes. Molecular Ecology. 2000, 9 (4): 487-488. 10.1046/j.1365-294x.2000.00887.x.View ArticlePubMedGoogle Scholar
- Price EW, Carbone I: SNAP: workbench management tool for evolutionary population genetic analysis. Bioinformatics. 2005, 21 (3): 402-404. 10.1093/bioinformatics/bti003.View ArticlePubMedGoogle Scholar
- Aylor DL, Price EW, Carbone I: SNAP: Combine and Map modules for multilocus population genetic analysis. Bioinformatics. 2006, 22 (11): 1399-1401. 10.1093/bioinformatics/btl136.View ArticlePubMedGoogle Scholar
- Tajima F: Statistical method for testing the neutral mutation hypothesis by DNA polymorphism. Genetics. 1989, 123 (3): 585-595.PubMed CentralPubMedGoogle Scholar
- Fu YX, Li WH: Statistical tests of neutrality of mutations. Genetics. 1993, 693-709. 133Google Scholar
- Statistical Tests of Neutrality of Mutations against Excess of Recent Mutations (Rare Allels). [http://hgc.sph.uth.tmc.edu/fu/genealogy/test2/welcome.html]
- Carbone I, Kohn LM: A microbial population-species interface: nested cladistic and coalescent inference with multilocus data. Molecular Ecology. 2001, 10 (4): 947-964. 10.1046/j.1365-294X.2001.01244.x.View ArticlePubMedGoogle Scholar
- Beerli P, Felsenstein J: Maximum likelihood estimation of a migration matrix and effective population sizes in n subpopulations by using a coalescent approach. PNAS. 2001, 98 (8): 4563-4568. 10.1073/pnas.081068098.PubMed CentralView ArticlePubMedGoogle Scholar
- Griffiths RC, Tavare S: Sampling Theory for Neutral Alleles in a Varying Environment. Philosophical Transactions Royal Society of London: Biological Sciences. 1994, 344 (1310): 403-410. 10.1098/rstb.1994.0079.View ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.