Search tips
Search criteria 


Logo of aemPermissionsJournals.ASM.orgJournalAEM ArticleJournal InfoAuthorsReviewers
Appl Environ Microbiol. 2016 November 1; 82(21): 6317–6325.
Published online 2016 October 14. Prepublished online 2016 August 19. doi:  10.1128/AEM.01708-16
PMCID: PMC5066348

Population Genetics of Hirsutella rhossiliensis, a Dominant Parasite of Cyst Nematode Juveniles on a Continental Scale

A. A. Brakhage, Editor
HKI and University of Jena


Hirsutella rhossiliensis is a parasite of juvenile nematodes, effective against a diversity of plant-parasitic nematodes. Its global distribution on various nematode hosts and its genetic variation for several geographic regions have been reported, while the global population genetic structure and factors underlying patterns of genetic variation of H. rhossiliensis are unclear. In this study, 87 H. rhossiliensis strains from five nematode species (Globodera sp., Criconemella xenoplax, Rotylenchus robustus, Heterodera schachtii, and Heterodera glycines) in Europe, the United States, and China were investigated by multilocus sequence analyses. A total of 280 variable sites (frequency, 0.6%) at eight loci and six clustering in high accordance with geographic populations or host nematode-associated populations were identified. Although H. rhossiliensis is currently recognized as an asexual fungus, recombination events were frequently detected. In addition, significant genetic isolation by geography and nematode hosts was revealed. Overall, our analyses showed that recombination, geographic isolation, and nematode host adaptation have played significant roles in the evolutionary history of H. rhossiliensis.

IMPORTANCE H. rhossiliensis has great potential for use as a biocontrol agent to control nematodes in a sustainable manner as an endoparasitic fungus. Therefore, this study has important implications for the use of H. rhossiliensis as a biocontrol agent and provides interesting insights into the biology of this species.


It is well known that fungi possess diverse strategies to obtain nutrients, as saprobes, pathogens (of plant and animal hosts), and symbionts (with other microbes, plants, or animals) (1). One of the important living strategies is parasitoidism, where parasitoids sterilize and/or kill their hosts before the hosts reach reproductive age. Plant-parasitic nematodes are worldwide agricultural pests and are responsible for global agricultural losses amounting to an estimated $157 billion each year (2). Endoparasitic fungi are predominantly soil-dwelling organisms that have the ability to sterilize plant-parasitic nematodes (3). They produce characteristic spores that adhere to, penetrate, and consume nematode bodies (4, 5).

Hirsutella rhossiliensis (Ophiocordycipitaceae, Hypocreales, Ascomycota) (6) lives in the soil, and it has been mainly isolated from five nematodes including Criconemella xenoplax (7), Heterodera schachtii (8), Heterodera glycines (9), Rotylenchus robustus (10), and Globodera (11). Though a dominant parasitoid of H. glycines J2s (second-stage juveniles) in the United States, mainly from Minnesota, H. rhossiliensis occurs less frequently in northeast China (12). In the agricultural fields, the plants, nematodes, and endoparasitic fungi constitute a tight tripartite system with the plants providing food to nematodes and the nematodes serving as nitrogen sources to endoparasitic fungi. This existing natural interaction suggests that endoparasitic fungi could have great potential for use as biocontrol agents to control nematodes in a sustainable manner (13,15).

To date, H. rhossiliensis has been detected in many localities in Europe, the United States, and China (11, 12). Genetic variation among different isolates of H. rhossiliensis has been detected using randomly amplified polymorphic DNA (RAPD) markers (16) and analyses of two housekeeping genes, internal transcribed spacer (ITS) of ribosomal DNA (rDNA) and the mitogen-activated protein kinase (MAPK) gene (13). Because of its global distribution and wide host range, H. rhossiliensis provides an ideal model to study the effects of geography and host nematodes on fungal population genetic divergence and factors driving genetic differentiation of fungi.

Investigating the evolutionary patterns and processes of nematode endoparasitic fungi over their geographical and ecological distribution provides valuable information on the radiation between nematode endoparasitic fungi and hosts (17). While phenotypic plasticity and genomic heterogeneity can lead to broad geographic and ecological distribution, low levels of gene flow often result in microbial speciation (18, 19). Population genetic analyses based on genotyping and genealogical analyses can help to infer past evolutionary events and evaluate the potential processes that generated the patterns of genetic variation (20, 21). An analysis of the population genetics of Hirsutella minnesotensis indicated that this dominant species in China had a clonal population structure and likely experienced a founder effect after colonization of China but without a significant bottleneck (22). However, the population genetic structure, possible geographic origins, and host nematode adaptation of H. rhossiliensis are largely unknown.

In recent years, single-nucleotide polymorphism (SNP) analyses have been widely used to detect intraspecific variation in microorganisms, including endoparasitic fungi (23, 24). When analyzing the genetic diversities of human fungal pathogens, for instance, Histoplasma capsulatum (25), Candida albicans (26), Aspergillus fumigatus (23), and Cryptococcus neoformans (27), it has been demonstrated that SNP analyses were highly discriminative in these species. Based on multilocus sequence analyses (MLSAs), genetic diversity in plant-pathogenic fungi was also reported for Magnaporthe grisea and Ustilaginoidea virens (28, 29).

Field experiments demonstrated that H. rhossiliensis was effective for control of H. glycines in greenhouses (30). Current biological control applications using H. rhossiliensis, however, have shown relatively limited success in agricultural fields. We think that an effective application of H. rhossiliensis or other fungi for biological control requires a thorough understanding of their population structure and their modes of reproduction. In this study, we used a multilocus sequence analysis (MLSA) scheme based on SNP markers to analyze population genetics of H. rhossiliensis by sampling from three geographic populations and five host nematode-associated populations. We aimed to understand (i) the level and partitioning of genetic variation within H. rhossiliensis, (ii) genetic structure and modes of reproduction in H. rhossiliensis, and (iii) whether geographic isolation and host nematode adaptation lead to rapid differentiation among populations of H. rhossiliensis. This study not only provides essential insights into the population differentiation and host adaptation of H. rhossiliensis but also helps scientists to better understand the ecological mechanism underlying the natural control of H. rhossiliensis in its role as a biocontrol agent against nematodes.


Sample collection and DNA extraction.

We employed all available H. rhossiliensis isolates (87 in total) around the world for this study. They were classified into three geographic populations (Europe, the United States, and China) and five host nematode-associated populations (Globodera, H. schachtii, C. xenoplax, R. robusta, and H. glycines). Except for 8 isolates from the CBS-KNAW Fungal Biodiversity Centre and 17 from the USDA-ARS Collection of Entomopathogenic Fungal Cultures (ARSEF), the other 62 isolates were all recovered by us from parasitized H. glycines in Heilongjiang Province, China, or northern Minnesota, USA (Table 1 and Fig. 1; see also Table S1 in the supplemental material).

Summary information for the Hirsutella rhossiliensis populations analyzed in this study
Map of the geographic distribution of H. rhossiliensis: United States, Europe, and China.

To isolate H. rhossiliensis from parasitized H. glycines, we collected rhizosphere soil samples from soybean fields in the main soybean-producing regions across northeast China from 2007 to 2008 and across Minnesota (USA) from 2010 to 2011. The sucrose flotation and centrifugation method (31) was used to extract the second-stage juveniles (J2s) that naturally occur in the soil, and these J2s were then transferred to wells of a 6-well tissue culture plate. The extracted J2s were examined with an inverted microscope (Olympus, Japan) at ×40 magnification, washed twice with sterilized water, and then placed on Difco potato dextrose agar (PDA) plates at 25°C for 2 weeks. All of the filamentous fungi isolated from these J2s were cultured on PDA plates covered with a piece of cellophane paper, and 0.05 to 0.1 g of mycelia was harvested and used to extract genomic DNA with the cetyltrimethylammonium bromide (CTAB) method (32). Fungal cultures were identified based on their morphology (30) and nonribosomal DNA (nrDNA) ITS sequences (13).

Molecular marker selection.

We initially chose three housekeeping genes as molecular markers: the nrDNA ITS region, the beta-tubulin gene (tub), and the translation elongation factor 1α gene (tef). However, the number of phylogenetically informative sites provided by the three fragments was insufficient. Therefore, novel DNA markers were determined by searching the genome data (15). Specifically, 27 fragments of genes putatively involved in the infection processes of H. rhossiliensis into nematodes and 15 fragments of interest as a result of a genomic comparison between H. minnesotensis and H. rhossiliensis were chosen as candidate markers. PCR primers were designed for each of the 42 fragments using the online software Primer3 ( From 12 randomly chosen isolates of different origins (CBS104.94, CBS105.94, ARSEF2006, ARSEF3746, ARSEF2788, ARSEF2789, ARSEF3755, ARSEF3756, CN-13-1, CN-30-1, USA-87-2, and USA-92-1) (see Table S1 in the supplemental material), we amplified and sequenced each of the 42 DNA fragments. Fragments with a high frequency of nucleotide variation were finally chosen as markers used in this study. A neutrality test for the selected markers was performed by DnaSP 5.10 (33).

PCR amplification and sequence determination.

Nine molecular markers, including ITS, tub, tef, and six newly identified markers, were amplified from each of the 87 H. rhossiliensis isolates. Each PCR mixture (50 μl) was composed of 5 μl of 10× EasyTaq buffer, 5 μl of deoxynucleoside triphosphate (dNTP) mixture (2.5 mM), 2 μl of each primer (10 μM), 1 μl of genomic DNA, 0.5 μl of EasyTaq DNA polymerase (TransGen Biotech, Beijing, China), and 36.5 μl of double-distilled water (ddH2O). The PCRs were performed in a TGradient thermocycler (Biometra, Germany) under the following conditions: denaturation at 94°C for 3 min; 35 cycles of denaturation at 94°C for 40 s, annealing at variable temperatures for 50 s, and elongation at 72°C for 1 min; and a final 10-min elongation step at 72°C. After confirmation of the PCR products by agarose gel electrophoresis, the amplicons were cleaned using the 3S Spin PCR product purification kit (Biocolor Bioscience & Technology Company, China). The purified PCR products were then sequenced using an ABI 3730 XL DNA sequencer (Applied Biosystems, USA) with Sanger dideoxy terminator sequencing. In the case of poor sequencing quality resulting from the direct sequencing of PCR products, the representative PCR fragments were cloned, and individually cloned fragments were then sequenced again.

Nucleotide variations.

Sequences of the fragments were aligned with the Muscle algorithm of MEGA 6.0 (34) and trimmed to remove ambiguously aligned regions. The nucleotide diversity (π), haplotype diversity, number of polymorphic sites, and allele numbers were estimated for both individual fragment and combined sequences using the program DnaSP 5.10 (33). By employing DnaSP 5.10, alleles at each locus were assigned and combined into an allelic profile designated a haplotype (HA) for each strain.

Phylogenetic and network analyses.

We performed a partition homogeneity test (PHT) to detect potential phylogenetic conflicts between different loci using PAUP version 4.0b10 (Sinauer Associates, Sunderland, MA, USA; D. Swofford, 2003). PHT used informative characters and simple stepwise-addition heuristic searches with 1,000 replicates. A P value lower than 0.001 indicated statistically significant differences (35). After identification of loci with no significant conflict, the phylogenic relationship among 87 strains was constructed using two methods: Bayesian inference (BI) and maximum likelihood (ML). To conduct phylogenetic analyses, we used the program PartitionFinder 1.1.1 (36) to evaluate the best partitioning scheme and determine suitable substitution models under the Bayesian information criterion (see Table S5 in the supplemental material). For partitioned data sets, model parameters across different partitions were unlinked, and the overall rate of evolution among partitions was allowed to vary. The Bayesian analysis, which was performed in MrBayes 3.2.2 (37), was initiated with random starting trees, run for 100 million generations with four incrementally heated chains (Metropolis-coupled Markov chain Monte Carlo [38]), and sampled at intervals of 10,000 generations. To avoid entrapment in the local optima, two independent Bayesian analyses were run, and the log-likelihood scores were compared for convergence (38). An ML analysis was performed in RAxML version 8.0.0 (39) and conducted under the GTRGAMMA model, and the ML support was assessed by 1,000 bootstrap replicates. Two H. minnesotensis strains, CN3608 and CBS113353, were used as the outgroup.

A phylogenetic network was constructed using the median-joining method implemented in the program SplitsTree 4.12 (40) to further assess the relationships among worldwide H. rhossiliensis haplotypes. This model-free method uses a parsimony approach based on pairwise differences to connect each sequence to its closest neighbor and allows for the creation of internal nodes (median vectors), which are interpreted as unsampled or extinct ancestral genotypes to link the existing genotypes in the most parsimonious way (41).

Population genetic analyses.

To determine whether the number of loci was sufficient to represent the genotypic diversity of the populations, we plotted the relationship between the number of loci and the genotype diversity using MultiLocus1.3b (42). A population genetic analysis was also performed within and between the three geographic populations and within and between the five host nematode-associated populations. For the within-population genetic analysis, we analyzed the genetic variation patterns with DnaSP 5.10, including the number of haplotypes, nucleotide diversity, and haplotype diversity. For the cross-population genetic analysis, pairwise differentiation coefficients (Fst) (43) were calculated and tested for significance against 1,000 bootstrap replicates using Arlequin 3.1 (44), and the average gene flow (Nm) (45) was calculated using DnaSP 5.10.

Tests for recombination.

Sexual reproduction generates random associations between alleles at different loci because genes from different individuals are mixed in a given population by sexual reproduction. Therefore, when recombination occurs among populations, random associations between alleles may be observed at different loci, even if sexual reproduction has not been observed in the natural populations (46). To identify evidence for recombination in H. rhossiliensis, two methods were used: the proportion of compatible pairs of loci (PrCP) (47) and the Φw test (or pairwise homoplasy index [PHI]) (48). Tests for recombination based on the principle of compatibility have been proven to be the most powerful (49, 50). Two loci are compatible (PrCP = 1) if it is possible to account for all the observed genotypes by mutations without inferring homoplasy (reversals, parallelisms, or convergences) or recombination; otherwise, the loci are incompatible (PrCP < 1). The Φw test means that phylogenetic incompatibility is an indicator of recombination at the population level; the lack of phylogenetic incompatibility, in contrast, implies asexual reproduction. The tests for PrCP were calculated using MultiLocus1.3b (42) with 1,000 randomizations; the Φw test was calculated by SplitsTree 4.12.

Tests for genetic isolation according to geography and host nematodes.

An analysis of molecular variance (AMOVA) was performed to determine the relative contribution of genetic variation among and within populations using GenAlEx 6.5 (51). The AMOVA was conducted for both the three geographic populations and the five host nematode-associated populations.

Mantel tests were performed in GenAlEx 6.5 for three groups of correlations: between the fungal genetic distance and longitudinal-latitudinal coordinates of geographic distance, between the fungal genetic distance and longitudinal distance, and between the fungal genetic distance and latitudinal distance. For the genetic distances, we used the pairwise population PhiPT values, which represent genetic distances among all pairs of populations as a tri-matrix, and these distances were then compared with the geographic distances between populations.

Accession number(s).

All of the nucleotide sequences obtained in this study have been submitted to the GenBank database, and the accession numbers are KP885720 to KP886014 (see also Table S2 in the supplemental material).


Nucleotide variations.

To determine appropriate markers, we successfully amplified the 42 candidate fragments from each of the 12 test isolates. Six fragments (HIR_00714, HIR_06686, HIR_08121, HIR_08679, HIR_01569, and HIR_09423) showed a high frequency of nucleotide variation and existed as single copies as indicated by BLAST searches against the H. rhossiliensis genome. The six fragments, together with three commonly used markers (nrDNA ITS, tub, and tef), were finally chosen to amplify from each of the 87 isolates (Table 2). Direct sequencing of PCR products from any of the nine fragments indicated no heterogeneity in sequencing chromatograms. Three to 11 alleles were identified in our samples depending on the locus. Moreover, the result of a neutrality test showed that all nine markers were neutral (Table 2).

Patterns of genetic variation at the individual loci and combined-locus data sets

A significant conflict was detected by PHT between HIR_00714 and the other loci (see Table S4 in the supplemental material). Therefore, HIR_00714 was excluded from the concatenated data set in the following analyses. From the concatenation of the remaining eight loci, a total of 280 variable sites within the combined 4,658-bp sequence and 44 multilocus haplotypes were identified from the 87 isolates (Table 2).

Phylogenetic and population structure analyses of H. rhossiliensis.

Bayesian (BI) and maximum likelihood (ML) phylogenetic trees were constructed based on the 8-locus concatenated data set. H. rhossiliensis strains tend to cluster by geography and nematode hosts, and each of the two phylogenetic approaches identified six clusters: European strains from Globodera sp. (NL-GL), American strains from C. xenoplax (USA-CX), American strains from R. robusta (USA-RR), American strains from H. schachtii (USA-HS), American strains from H. glycines (USA-HG), and Chinese (northeastern) strains from H. glycines (CN-HG). However, there were subtle topological differences between the two approaches. For example, two American strains from H. schachtii (CBS567.92 and ARSEF3761) did not group into the USA-HS subclade, and the positions of two Chinese strains from H. glycines (CN-30-2-JX and CN-3-10-JL) and one American strain from C. xenoplax (ARSEF3746) had inconsistent placements between the BI and ML trees (Fig. 2; see also Fig. S2 in the supplemental material).

Bayesian tree of the H. rhossiliensis strains based on the concatenated 8-locus data set. The strict consensus tree from phylogenetic analysis was constructed using a Bayesian inference method. The numbers above the branches indicate the posterior probability ...

The network diagram showed a linear structure among the 87 strains (Fig. 3), with Chinese strains from H. glycines, American strains from H. glycines, European strains from Globodera sp., American strains from H. schachtii, American strains from R. robusta, and American strains from C. xenoplax each forming an individual cluster.

Median joining of phylogenetic network inferred with H. rhossiliensis strain concatenated 8-locus data set. This network includes all of the most parsimonious trees linking the isolates. Each unique haplotype is represented by a block, with the block's ...

Genetic differentiation and gene flow among geographic or host nematode-associated populations.

A high value for the genetic differentiation coefficient Fst (P < 0.01) revealed that significant genetic differentiation occurred among all of the geographic or host nematode-associated populations except between NL-P and USA-P, the value for which was not significant (P > 0.01) (Tables 3 and and4).4). Consistently, the estimated gene flow was low among the three geographic populations (Nm = 0.14) and among the five host nematode-associated populations (Nm = 0.28).

Population genetic parameters for each of three geographic populations of H. rhossiliensis
Population genetic parameters for each of five host nematode-associated populations of H. rhossiliensis

Evidence for recombination.

The PrCP index provided evidence for recombination except for two populations, American strains from C. xenoplax (USA-CX) and American strains from R. robusta (USA-RR). When the 87 strains were analyzed together, recombination signals were found in 3/5 locus pairs (Table 3; see also Table S6 in the supplemental material).

Similarly, the Φw test found statistically significant evidence for recombination (P < 0.001). However, there was a slight difference in the Φw results for recombination among the three geographic populations and five host nematode-associated populations compared with the PrCP results. For example, the Φw test did not find statistically significant evidence for recombination in the USA-P, USA-CX, USA-HS, or USA-RR population (Table 3).

Genetic isolation by geography or host nematodes.

The AMOVA results showed that 44% of the genetic variation could be attributed to variation among the three geographic populations and 56% could be attributed to variation among the individual isolates within populations (Table 5). For the five host nematode-associated populations, the AMOVA results showed that 61% and 39% of the genetic variation could be attributed to variation among the populations and among individual isolates within populations, respectively. Both of these sources of variation were significant (P < 0.01) (Table 6).

Summary results of the analysis of molecular variance (AMOVA) within and among three geographic populations of H. rhossiliensisa
Summary results of the analysis of molecular variance (AMOVA) within and among five host nematode-associated populations of H. rhossiliensisa

The Mantel tests revealed significant genetic isolation by distance irrespective of the horizontal geographical distances (based on longitudinal-latitudinal coordinates, along a longitudinal gradient, or along a latitudinal gradient) or population divisions (three geographic populations) employed (Table 7).

Mantel test results between genetic distances and geographic distances among three geographic populations


We collected 87 H. rhossiliensis strains from five host nematodes in Europe, the United States, and China. Although a large number of strains were obtained from H. glycines, at least three strains were included from each of the other nematode hosts (Table 1; see also Table S1 in the supplemental material). We would explain that our samples from four nematodes (Globodera sp., H. schachtii, C. xenoplax, and R. robusta) were isolated in the 1980s and 1990s, while samples from H. glycines were isolated during 2007 to 2011. Because some nematodes did not occur in China or we did not find this fungus from them (11, 52,54), we failed to isolate H. rhossiliensis from other host nematodes except from H. glycines. It was also difficult to obtain strains from the above four host nematodes (except H. glycines) in the United States. The cultures from culture collection centers (CBS and ARSEF) were preserved in lyophilized forms, which largely decreased strain mutations over the past 20 to 30 years (see Table S1 in the supplemental material). Although there might be some genetic variations of fungus in nature, we believe that the variation during 20 to 30 years could not affect the main results obtained from this study.

In this study, recombination was found in most of the H. rhossiliensis populations and when the whole samples were considered. It can be inferred and supported by two lines of evidence, the PrCP values and the Φw test (Tables 3 and and4).4). The PrCP test measures whether two loci are compatible, and the Φw test examines phylogenetic incompatibility. So, it is reasonable that there is conflict between the two results (49). Indeed, the sexual cycles of most fungi are difficult to observe in nature (55), and inferences of the potential for a sexual cycle to occur have largely relied on analyses of gene and genotype frequencies in natural populations. The sexual stage of H. rhossiliensis has not been discovered in nature or in the laboratory, although we detected mating type genes in its genome (unpublished data). In plant fungal pathogens, a high level of recombination provides an advantage in rapidly generating many new combinations of virulence genes to counterbalance corresponding resistance genes in the host (56, 57). In arbuscular mycorrhizal fungi, recombination or recombination-like events in addition to clonality have greatly contributed to genetic diversity (58). This study has demonstrated evidence of recombination among H. rhossiliensis strains, and recombination appears to be one of the factors shaping the evolution of H. rhossiliensis. Similar results were also observed in Candida albicans, which showed both clonality and recombination, even though a complete sexuality stage is not known to exist in this fungus (59). However, sexual recombination may not be the sole factor because we cannot exclude the possibility that parasexual recombination has contributed to the observed linkage equilibrium and phylogenetic incompatibility (60).

Our population genetics analyses of H. rhossiliensis were based on two types of population divisions: three geographic populations (NL-P, USA-P, and CN-P) and five host nematode-associated populations (NL-GL, USA-CX, USA-RR, USA-HS, and USA/CN-HG) (Table 1). The significance of geography and the host nematodes in shaping the structure of H. rhossiliensis populations on a large scale was supported by the strong divergent signals in the phylogenetic and haplotype network analyses (Fig. 2 and and3).3). In this study, a significant positive correlation was observed between H. rhossiliensis genetic distance and latitudinal and longitudinal geographic distance, thus indicating that geography is an important barrier to gene flow and has promoted the genetic divergence within fungal species (Table 7). Moreover, differentiation among the three geographic populations was also observed from high Fst values and low Nm values. We suggest that the North Pacific Ocean may be a boundary separating the Chinese population and the American population of H. rhossiliensis because of limited opportunities to disperse via wind or soil. A similar geography-based separation has been reported in Colletotrichum truncatum from chili peppers in China that were genetically differentiated into southern and northern populations (61). Previously, we found significant differences among local populations of H. minnesotensis from China and the United States, where isolation by geography plays a key role in genetic differentiation (22). Based on these results, we conclude that geographic distance contributes to driving population genetic differentiation of H. rhossiliensis.

Among the five host nematode-associated populations studied here, we identified positive correlations between H. rhossiliensis genetic distance at different latitudinal and longitudinal distances. We observed 61% of the molecular genetic variance among those nematode-associated populations. This might result from host-related selection driving population genetic differentiation at neutral and selected loci (62). Another explanation was that host shifts might occur frequently between neighboring nematodes, which increased the mating probability of H. rhossiliensis across the nearby hosts and also increased the genetic variation within and among populations.

Our analyses revealed that both nematode hosts and geography contributed to genetic differentiation among H. rhossiliensis populations (Tables 5 to to7).7). We, however, could not perform a two-way AMOVA to further look at their relative effects because population assignments of our samples did not meet the requirement for performing such an analysis. In this study, we confirmed that difference in host nematode was a key factor in regulating genetic divergence of nematode parasitic fungus from the phylogenetic approaches and the haplotype network (Fig. 2 and and3).3). Similarly, the strong association of genetic groups with different geographic locations suggested the important role of geographic isolation. Previous studies on microorganisms such as Bacillus simplex concluded that genetic differentiation was strongly attributed to geographic locations (63). Numerous studies found that geographic region and host species contribute equally to the population genetic differentiation, for instance, Grosmannia clavigera (64), Claviceps purpurea (65), Venturia inaequalis (66), Arthrobotrys oligospora (67), Colletotrichum gloeosporioides (68), and Metarhizium anisopliae (69).

In conclusion, our study identified the population genetic structure of the nematode endoparasitic fungus H. rhossiliensis and found a high level of genetic variation among its populations worldwide. Using multiple sources of evidence, we demonstrated that recombination among strains, isolation by geography, and isolation based on the host nematode were factors driving the genetic differentiation of H. rhossiliensis. Therefore, the findings of this study on the population structure of H. rhossiliensis and the associated multiple driving factors affecting its genetic diversity provide an essential understanding of the mechanisms of evolution and differentiation of the fungus on a global scale. These findings are expected to provide important insights into nematode control by this fungus.

Supplementary Material

Supplemental material:


This study was supported by the National Natural Science Foundation of China (grant no. 31430071) and the National Basic Research Program of China (973 program, grant no. 2013CB127506).

We thank Richard A. Humber, who kindly provided 15 strains from the USDA-ARS Collection of Entomopathogenic Fungal Cultures (ARSEF) database at the USDA-ARS Biological Integrated Pest Management Research Unit, Cornell University, New York, USA.


Supplemental material for this article may be found at


1. Cantrell SA, Claribel BF 2010. Fungal molecular diversity of a Puerto Rican subtropical hypersaline microbial mat. Fungal Ecol 3:402–405. doi:.10.1016/j.funeco.2010.04.001 [Cross Ref]
2. Abad P, Gouzy J, Aury J, Castagnone-Sereno P, Danchin E, Deleury E, Perfus-Barbeoch L, Anthouard V, Artiguenave F, Blok VC, Caillaud MC, Coutinho PM, Dasilva C 2008. Genome sequence of the metazoan plant-parasitic nematode Meloidogyne incognita. Nat Biotechnol 26:909–915. doi:.10.1038/nbt.1482 [PubMed] [Cross Ref]
3. Li TF, Zhang KQ, Liu XZ 2000. Taxonomy of nematophagous fungi. China Science and Technology Press, Beijing, China.
4. Liu XZ, Chen SY 2001. Screening isolates of Hirsutella species for biocontrol of Heterodera glycines. Biocontrol Sci Technol 11:151–160. doi:.10.1080/09583150020029826 [Cross Ref]
5. Liu XZ, Xiang MC, Che YS 2009. The living strategy of nematophagous fungi. Mycoscience 50:20–25. doi:.10.1007/S10267-008-0451-3 [Cross Ref]
6. Minter DW, Brady BL 1980. Mononematous species of Hirsutella. Trans Br Mycol Soc 74:271–282. doi:.10.1016/S0007-1536(80)80157-4 [Cross Ref]
7. Jaffee BA, Zehr EI 1983. Effects of certain solutions, osmotic potential, and soil solutions on parasitism of Criconemella xenoplax by Hirsutella rhossiliensis. Phytopathology 73:544–546. doi:.10.1094/Phyto-73-544 [Cross Ref]
8. Müller J. 1984. The influence of two pesticides on fungal parasites of Heterodera schachtii. Colloq INRA 31:225–231.
9. Chen SY, Reese CD 1999. Parasitism of the nematode Heterodera glycines by the fungus Hirsutella rhossiliensis influenced by crop sequence. J Nematol 31:437–444. [PMC free article] [PubMed]
10. Tedford EC, Jaffee BA 1995. In vitro parasitism of Rotylenchus robustus by isolates of Hirsutella rhossiliensis. J Nematol 22:486–489. [PMC free article] [PubMed]
11. Dobosz R, Obrepalska-Steplowska A, Kornobis S 2006. Globodera artemisiae (Eroshenko et Kazachenko, 1972) (Nematoda: Heteroderidae) from Poland. J Plant Prot Res 46:403–407.
12. Ma R, Liu XZ, Jian H, Li SD 2005. Detection of Hirsutella spp. and Pasteuria sp. parasitizing second-stage juveniles of Heterodera glycines in soybean fields in China. Biol Control 33:223–229. doi:.10.1016/j.biocontrol.2005.03.004 [Cross Ref]
13. Xiang MC, Yang XH, Wang ZX, Liu XZ, Chen SY, Xiao QM 2007. Variability of morphology, parasitism, and nucleotide sequences among isolates and species of nematophagous Hirsutella. Biol Control 41:110–119. doi:.10.1016/j.biocontrol.2006.12.016 [Cross Ref]
14. Timper P. 2011. Uitilization of biological control for managing plant-parasitic nematodes, p 259–289. In Davies K, Spiegel Y (ed), Biological control of plant-parasitic nematodes: building coherence between microbial ecology and molecular mechanisms, vol 11. Springer, Dordrecht, Netherlands.
15. Lai YL, Liu KK, Zhang XY, Zhang XL, Li K, Wang NN, Shu C, Wu YP, Wang CS, Bushley KE, Xiang MC, Liu XZ 2014. Comparative genomics and transcriptomics analyses reveal divergent lifestyle features of nematode endoparasitic fungus Hirsutella minnesotensis. Genome Biol Evol 6:3077–3093. doi:.10.1093/gbe/evu241 [PMC free article] [PubMed] [Cross Ref]
16. Tedford EC, Jaffee BA, Muldoon AE 1994. Variability among isolates of the nematophagous fungus Hirsutella rhossiliensis. Mycol Res 98:1127–1136.
17. McDonald BA, Linde C 2002. Pathogen population genetics, evolutionary potential and durable resistance. Annu Rev Phytopathol 40:349–379. doi:.10.1146/annurev.phyto.40.120501.101443 [PubMed] [Cross Ref]
18. Stukenbrock E. 2013. Evolution, selection and isolation: a genomic view of speciation in fungal plant pathogens. New Phytol 199:895–907. doi:.10.1111/nph.12374 [PubMed] [Cross Ref]
19. Cheng K, Rong XY, Pinto-Tomás AA, Fernández-Villalobos M, Murillo-Cruz C, Huang Y 2015. Population genetic analysis of Streptomyces albidoflavus reveals habitat barriers to homologous recombination in the diversification of streptomycetes. Appl Environ Microbiol 81:966–3975. doi:.10.1128/AEM.02925-14 [PMC free article] [PubMed] [Cross Ref]
20. Giraud T, Refrégier G, Le Gac M, de Vienne DM, Hood ME 2008. Speciation in fungi. Fungal Genet Biol 45:791–802. doi:.10.1016/j.fgb.2008.02.001 [PubMed] [Cross Ref]
21. Tittes S, Kane SC 2014. The genomics of adaptation, divergence and speciation: a congealing theory. Mol Ecol 23:3938–3940. doi:.10.1111/mec.12855 [PubMed] [Cross Ref]
22. Shu C, Jiang XZ, Cheng XL, Wang NN, Chen SY, Xiang MC, Liu XZ 2015. Genetic structure and parasitization-related ability divergence of a nematode fungal pathogen Hirsutella minnesotensis following founder effect in China. Fungal Genet Biol 81:212–220. doi:.10.1016/j.fgb.2015.02.005 [PubMed] [Cross Ref]
23. Bain JM, Tavanti A, Davidson AD, Jacobsen MD, Shaw D, Gow NA, Odds FC 2007. Multilocus sequence typing of the pathogenic fungus Aspergillus fumigatus. J Clin Microbiol 45:1469–1477. doi:.10.1128/JCM.00064-07 [PMC free article] [PubMed] [Cross Ref]
24. Nunney L, Vickerman DB, Bromley RE, Russell SA, Hartman JR, Morano LD, Stouthamer R 2013. Recent evolutionary radiation and host plant specialization in the Xylella fastidiosa subspecies native to the United States. Appl Environ Microbiol 79:2189–2200. doi:.10.1128/AEM.03208-12 [PMC free article] [PubMed] [Cross Ref]
25. Kasuga T, White TJ, Koenig G, McEwen J, Restrepo A, Castañeda E, Da Silva Lacaz C, Heins-Vaccari EM, De Freitas RS, Zancopé-Oliveira RM, Qin Z, Negroni R, Carter DA, Mikami Y, Tamura M, Taylor ML, Miller GF, Poonwan N, Taylor JW 2003. Phylogeography of the fungal pathogen Histoplasma capsulatum. Mol Ecol 12:3383–3401. doi:.10.1046/j.1365-294X.2003.01995.x [PubMed] [Cross Ref]
26. Li J, Fan SR, Liu XP, Li DM, Nie ZH, Li F, Lin H, Huang WM, Zong LL, Lei H, Bai FY 2008. Biased genotype distributions of Candida albicans strains associated with vulvovaginal candidosis and candidal balanoposthitis in China. Clin Infect Dis 47:1119–1125. doi:.10.1086/592249 [PubMed] [Cross Ref]
27. Meyer W, Aanensen DM, Boekhout T, Cogliati M, Diaz MR, Esposto MC, Fisher M, Gilgado F, Hagen F, Kaocharoen S, Litvintseva AP, Mitchell TG, Simwami SP, Trilles L, Viviani MA, Kwon-Chung J 2009. Consensus multi-locus sequence typing scheme for Cryptococcus neoformans and Cryptococcus gattii. Med Mycol 47:561–570. doi:.10.1080/13693780902953886 [PMC free article] [PubMed] [Cross Ref]
28. Choi J, Park SY, Kim BR, Roh JH, Oh IS, Han SS, Lee YH 2013. Comparative analysis of pathogenicity and phylogenetic relationship in Magnaporthe grisea species complex. PLoS One 8:e57196. doi:.10.1371/journal.pone.0057196 [PMC free article] [PubMed] [Cross Ref]
29. Sun XY, Kang S, Zhang YJ, Tan XQ, Yu YF, He HY, Zhang XY, Liu YF, Wang S, Sun WX, Cai L, Li SJ 2013. Genetic diversity and population structure of rice pathogen Ustilaginoidea virens in China. PLoS One 8:e76879. doi:.10.1371/journal.pone.0076879 [PMC free article] [PubMed] [Cross Ref]
30. Chen SY, Liu XZ 2005. Control of the soybean cyst nematode by the fungi Hirsutella rhossiliensis and Hirsutella minnesotensis in greenhouse studies. Biol Control 32:208–219. doi:.10.1016/j.biocontrol.2004.09.013 [Cross Ref]
31. Jenkin WR. 1964. A rapid centrifugal-flotation technique for separating nematodes from soil. Plant Dis Rep 48:692.
32. Murray MG, Thompson WF 1980. Rapid isolation of high molecular weight plant DNA. Nucleic Acids Res 8:4321–4326. [PMC free article] [PubMed]
33. Librado P, Rozas J 2009. DnaSP v5: a software for comprehensive analysis of DNA polymorphism data. Bioinformatics 25:1451–1452. doi:.10.1093/bioinformatics/btp187 [PubMed] [Cross Ref]
34. Tamura K, Stecher G, Peterson D, Filipski A, Kumar S 2013. MEGA6: molecular evolutionary genetics analysis version 6.0. Mol Biol Evol 30:2725–2729. doi:.10.1093/molbev/mst197 [PMC free article] [PubMed] [Cross Ref]
35. Cunningham CW. 1997. Can three incongruence tests predict when data should be combined? Mol Biol Evol 14:733–740. [PubMed]
36. Lanfear R, Calcott B, Ho SY, Guindon S 2012. Partitionfinder: combined selection of partitioning schemes and substitution models for phylogenetic analyses. Mol Biol Evol 29:1695–1701. doi:.10.1093/molbev/mss020 [PubMed] [Cross Ref]
37. Ronquist F, Teslenko M, van der Mark P, Ayres DL, Darling A, Höhna S, Larget B, Liu L, Suchard MA, Huelsenbeck JP 2012. MrBayes 3.2: efficient Bayesian phylogenetic inference and model choice across a large model space. Syst Biol 61:539–542. doi:.10.1093/sysbio/sys029 [PMC free article] [PubMed] [Cross Ref]
38. Huelsenbeck JP, Bollback JP 2001. Empirical and hierarchical Bayesian estimation of ancestral states. Syst Biol 50:351–366. doi:.10.1080/106351501300317978 [PubMed] [Cross Ref]
39. Stamatakis A. 2014. RAxML version 8: a tool for phylogenetic analysis and post-analysis of large phylogenies. Bioinformatics 30:1312–1313. doi:.10.1093/bioinformatics/btu033 [PMC free article] [PubMed] [Cross Ref]
40. Huson DH, Bryant D 2006. Application of phylogenetic networks in evolutionary studies. Mol Biol Evol 23:254–267. [PubMed]
41. Bandelt HJ, Forster P, Rohl A 1999. Median-joining networks for inferring intraspecific phylogenies. Mol Biol Evol 16:37–48. [PubMed]
42. Agapow PM, Burt A 2001. Indices of multilocus linkage disequilibrium. Mol Ecol Notes 1:101–102. doi:.10.1046/j.1471-8278.2000.00014.x [Cross Ref]
43. Weir BS, Cockerham CC 1984. Estimating F-statistics for the analysis of population structure. Evolution 38:1358–1370.
44. Excoffier L, Laval G, Schneider S 2005. Arlequin version 3.0: an integrated software package for population genetics data analysis. Evol Bioinform 1:47–50. [PMC free article] [PubMed]
45. Slatkin M. 1987. Gene flow and the geographical structure of natural populations. Science 236:787–792. [PubMed]
46. Xu JP. 2006. Fundamentals of fungal molecular population genetic analyses. Curr Issues Mol Biol 8:75–89. [PubMed]
47. Liang JF, Xu J, Yang ZL 2009. Divergence, dispersal and recombination in Lepiota cristata from China. Fungal Divers 38:105.
48. Bruen TC, Philippe H, Bryant D 2006. A simple and robust statistical test for detecting the presence of recombination. Genetics 172:2665–2681. [PubMed]
49. Posada D, Crandall KA 2001. Evaluation of methods for detecting recombination from DNA sequences: computer simulations. Proc Natl Acad Sci U S A 98:13757–13762. doi:.10.1073/pnas.241370698 [PubMed] [Cross Ref]
50. Wiuf C, Christensen T, Hein J 2001. A simulation study of the reliability of recombination detection methods. Mol Biol Evol 18:1929–1939. doi:.10.1093/oxfordjournals.molbev.a003733 [PubMed] [Cross Ref]
51. Peakall R, Smouse PE 2012. GenAlEx 6.5: genetic analysis in Excel. Population genetic software for teaching and research—an update. Bioinformatics 19:2537–2539. [PMC free article] [PubMed]
52. Nyczepir AP, Zehr EI, Lewis SA, Harshman DC 1983. Short life of peach trees induced by Criconemella xenoplax. Plant Dis 67:507–508.
53. Atighi MR, Pourjam E, Pedram M, Cantalapiedra-Navarrete C, Palomares-Rius JE, Castillo P 2011. Molecular and morphological characterisations of two new species of Rotylenchus (Nematoda: Hoplolaimidae) from Iran. Nematology 13:951–964. doi:.10.1163/138855411X571795 [Cross Ref]
54. Amiri S, Subbotin SA, Moens M 2002. Identification of the beet cyst nematode Heterodera schachtii by PCR. Eur J Plant Pathol 108:497–506. doi:.10.1023/A:1019974101225 [Cross Ref]
55. Calo S, Billmyre RB, Heitman J 2013. Generators of phenotypic diversity in the evolution of pathogenic microorganisms. PLoS Pathog 9:e1003181. doi:.10.1371/journal.ppat.1003181 [PMC free article] [PubMed] [Cross Ref]
56. Zhan J, Pettway RE, McDonald BA 2003. The global genetic structure of the wheat pathogen Mycosphaerella graminicola is characterized by high nuclear diversity, low mitochondrial diversity, regular recombination, and gene flow. Fungal Genet Biol 38:286–297. doi:.10.1016/S1087-1845(02)00538-8 [PubMed] [Cross Ref]
57. Souza EA, Camargo OA Jr, Pinto JM 2010. Sexual recombination in Colletotrichum lindemuthianum occurs on a fine scale. Genet Mol Res 9:1759–1769. doi:.10.4238/vol9-3gmr863 [PubMed] [Cross Ref]
58. Vandenkoornhuyse P, Leyval C, Bonnin I 2001. High genetic diversity in arbuscular mycorrhizal fungi: evidence for recombination events. Heredity 87:243–253. doi:.10.1046/j.1365-2540.2001.00941.x [PubMed] [Cross Ref]
59. Graser Y, Volovsek M, Arrington O, Schonian G, Presber W, Mitchell TG, Vilgalys AR 1996. Molecular markers reveal that population structure of the human pathogen Candida albicans exhibits both clonality and recombination. Proc Natl Acad Sci U S A 93:12473–12477. [PubMed]
60. Noguchi MT, Yasuda N, Fujita Y 2006. Evidence of genetic exchange by parasexual recombination and genetic analysis of pathogenicity and mating type of parasexual recombinants in rice blast fungus, Magnaporthe oryzae. Phytopathology 96:746–750. doi:.10.1094/PHYTO-96-0746 [PubMed] [Cross Ref]
61. Diao YZ, Zhang C, Xu JP, Lin D, Liu L, Mtung'e OG, Liu XL 2015. Genetic differentiation and recombination among geographic populations of the fungal pathogen Colletotrichum truncatum from chili peppers in China. Evol Appl 8:108–118. doi:.10.1111/eva.12233 [PMC free article] [PubMed] [Cross Ref]
62. Orsini L, Vanoverbeke J, Swillen I, Mergeay J, Meester LD 2013. Drivers of population genetic differentiation in the wild: isolation by dispersal limitation, isolation by adaptation and isolation by colonization. Mol Ecol 22:5983–5999. doi:.10.1111/mec.12561 [PubMed] [Cross Ref]
63. Sikorski J, Nevo E 2005. Adaptation and incipient sympatric speciation of Bacillus simplex under microclimatic contrast at “Evolution Canyons” I and II, Israel. Proc Natl Acad Sci U S A 102:15924–15929. doi:.10.1073/pnas.0507944102 [PubMed] [Cross Ref]
64. Alamouti SM, Wang V, DiGuistini S, Six DL, Bohlmann J, Hamelin RC 2011. Gene genealogies reveal cryptic species and host preferences for the pine fungal pathogen Grosmannia clavigera. Mol Ecol 20:2581–2602. doi:.10.1111/j.1365-294X.2011.05109.x [PubMed] [Cross Ref]
65. Douhan GW, Smith ME, Huyrn KL, Westbrook A, Beerli P, Fisher AJ 2008. Multigene analysis suggests ecological speciation in the fungal pathogen Claviceps purpurea. Mol Ecol 17:2276–2286. doi:.10.1111/j.1365-294X.2008.03753.x [PMC free article] [PubMed] [Cross Ref]
66. Gladieux P, Guérin F, Giraud T, Caffier V, Lemaire C, Parisi L, Didelot F, Cam LB 2011. Emergence of novel fungal pathogens by ecological speciation: importance of the reduced viability of immigrants. Mol Ecol 20:4521–4532. doi:.10.1111/j.1365-294X.2011.05288.x [PubMed] [Cross Ref]
67. Zhang Y, Qiao M, Xu J, Cao Y, Zhang KQ, Yu ZF 2013. Genetic diversity and recombination in natural populations of the nematode-trapping fungus Arthrobotrys oligospora from China. Ecol Evol 3:312–325. doi:.10.1002/ece3.450 [PMC free article] [PubMed] [Cross Ref]
68. Gazis R, Rehner S, Chaverri P 2011. Species delimitation in fungal endophyte diversity studies and its implications in ecological and biogeographic inferences. Mol Ecol 20:3001–3013. doi:.10.1111/j.1365-294X.2011.05110.x [PubMed] [Cross Ref]
69. Bidochka MJ, Kamp AM, Lavender TM, Dekoning J, De Croos JA 2001. Habitat association in two genetic groups of the insect-pathogenic fungus Metarhizium anisopliae: uncovering cryptic species? Appl Environ Microbiol 67:1335–1342. doi:.10.1128/AEM.67.3.1335-1342.2001 [PMC free article] [PubMed] [Cross Ref]

Articles from Applied and Environmental Microbiology are provided here courtesy of American Society for Microbiology (ASM)