Recombination between Clonal Lineages of the Asexual Fungus Verticillium dahliae Detected by Genotyping by Sequencing

Most asexual species of fungi have either lost sexuality recently, or they experience recombination by cryptic sexual reproduction. Verticillium dahliae is a plant-pathogenic, ascomycete fungus with no known sexual stage, even though related genera have well-described sexual reproduction. V. dahliae reproduces mitotically and its population structure is highly clonal. However, previously described discrepancies in phylogenetic relationships among clonal lineages may be explained more parsimoniously by recombination than mutation; therefore, we looked for evidence of recombination within and between clonal lineages. Genotyping by sequencing was performed on 141 V. dahliae isolates from diverse geographic and host origins, resulting in 26,748 single-nucleotide polymorphisms (SNPs). We found a strongly clonal population structure with the same lineages as described previously by vegetative compatibility groups (VCGs) and molecular markers. We detected 443 recombination events, evenly distributed throughout the genome. Most recombination events detected were between clonal lineages, with relatively few recombinant haplotypes detected within lineages. The only three isolates with mating type MAT1-1 had recombinant SNP haplotypes; all other isolates had mating type MAT1-2. We found homologs of eight meiosis-specific genes in the V. dahliae genome, all with conserved or partially conserved protein domains. The extent of recombination and molecular signs of sex in (mating-type and meiosis-specific genes) suggest that V. dahliae clonal lineages arose by recombination, even though the current population structure is markedly clonal. Moreover, the detection of new lineages may be evidence that sexual reproduction has occurred recently and may potentially occur under some circumstances. We speculate that the current clonal population structure, despite the sexual origin of lineages, has arisen, in part, as a consequence of agriculture and selection for adaptation to agricultural cropping systems.     

Genetic diversification of the chestnut blight fungus Cryphonectria parasitica and its associated hypovirus in Germany

Chestnut blight in south-western Germany was first reported in 1992 and is since expanding in distribution. Here we investigated the invasion history of Cryphonectria parasitica and its associated hypovirus. For this, we characterized 284 isolates collected between 1992 and 2012 for hypovirulence, vegetative compatibility (vc), mating type, and microsatellite haplotype. A total of 27 haplotypes and 15 vc types were observed, although the C. parasitica population analyzed is currently dominated to 50 % by one haplotype and to 64 % by the vc type EU-2. Structure analysis indicated two divergent genetic pools. Over 66 % of the haplotypes belonged to a pool probably originating from northern Italy. Further diversification is expected due to ongoing sexual recombination, but also to new migration and additional introductions. Cryphonectria hypovirus 1 (CHV-1) was found in four of five C. parasitica populations from Baden-Württemberg. Genetic analysis of the 35 CHV-1 isolates obtained revealed that they all belong to the German subtype, although they have clearly diverged from the first German hypovirus isolated in 1992. Our study suggests that C. parasitica has been introduced into Germany several times from two different gene pools, whereas the hypovirus most probably has a single origin.     

Clonal population structure of the chestnut blight fungus in expanding ranges in southeastern Europe

Expanding populations are often less genetically diverse at their margins than at the centre of a species' range. Established, older populations of the chestnut blight fungus, Cryphonectria parasitica, are more variable for vegetative compatibility (vc) types than in expanding populations in southeastern Europe where C. parasitica has colonized relatively recently. To test whether vc types represent clones, we genotyped 373 isolates of C. parasitica from southern Italy, Romania, Bulgaria, Macedonia, Greece and Turkey using 11 sequence-characterized amplified region (SCAR) markers. Ten SCAR loci and six vegetative incompatibility (vic) loci were polymorphic in these samples. These populations are clonal by all criteria tested: (i) among 373 isolates, we found only eight multilocus haplotypes, and the same haplotypes were found in multiple countries, sometimes separated in time by as much as 12 years; (ii) the number of haplotypes observed was significantly less than expected under random mating; (iii) populations are in linkage disequilibrium; (iv) the two sets of independent markers, SCARs and vc types, are highly correlated; and (v) sexual structures of C. parasitica were found only in Bulgaria and Romania. One mating type (MAT-1) was found in 98% of the isolates sampled. In contrast, a population in northern Italy, in the central part of the range in Europe, had 12 multilocus haplotypes among 19 isolates. The spread of a few clones could be the result either of founder effect and restricted migration, or these clones have greater fitness than others and spread because they are better adapted to conditions in southeastern Europe.     

Molecular phylogenetics of Candida albicans

We analyzed data on multilocus sequence typing (MLST), ABC typing, mating type-like locus (MAT) status, and antifungal susceptibility for a panel of 1,391 Candida albicans isolates. Almost all (96.7%) of the isolates could be assigned by MLST to one of 17 clades. eBURST analysis revealed 53 clonal clusters. Diploid sequence type 69 was the most common MLST strain type and the founder of the largest clonal cluster, and examples were found among isolates from all parts of the world. ABC types and geographical origins showed statistically significant variations among clades by univariate analysis of variance, but anatomical source and antifungal susceptibility data were not significantly associated. A separate analysis limited to European isolates, thereby minimizing geographical effects, showed significant differences in the proportions of isolates from blood, commensal carriage, and superficial infections among the five most populous clades. The proportion of isolates with low antifungal susceptibility was highest for MAT homozygous a/a types and then alpha/alpha types and was lowest for heterozygous a/alpha types. The tree of clades defined by MLST was not congruent with trees generated from the individual gene fragments sequenced, implying a separate evolutionary history for each fragment. Analysis of nucleic acid variation among loci and within loci supported recombination. Computational haplotype analysis showed a high frequency of recombination events, suggesting that isolates had mixed evolutionary histories resembling those of a sexually reproducing species.     

Evaluation of the potential for sexual reproduction in field populations of Cercospora beticola from USA

Cercospora leaf spot, caused by the hemibiotrophic fungal pathogen Cercospora beticola, is the most economically damaging foliar disease of sugarbeet worldwide. Although most C. beticola populations display characteristics reminiscent of sexual recombination, no teleomorph has been described. To assess whether populations in northern United States have characteristics consistent with sexual reproduction, 1024 isolates collected over a 3-y period were analyzed for frequency and distribution of mating type genes. After clone correction, an approximately equal distribution of mating types was found for each sampling year. Mating type frequency was also assessed in individual lesions. Lesions always consisted of isolates with a single mating type and microsatellite haplotype, but both mating types and up to five microsatellite haplotypes could be found on an individual leaf. The MAT1-1-1 and MAT1-2-1 genes were sequenced from 28 MAT1-1 and 28 MAT1-2 isolates, respectively. Three MAT1-1-1 nucleotide haplotypes were identified that encoded a single amino acid sequence. For MAT1-2-1, five nucleotide haplotypes were identified that encoded four protein variants. MAT1-1-1 and MAT1-2-1 gene expression analyses were conducted on plants inoculated with either or both mating types. MAT1-1-1 expression remained low, but MAT1-2-1 spiked during late stages of colonization. A segment of the MAT1-2-1 coding sequence was also found in MAT1-1 isolates. Taken together, these results suggest that C. beticola has the potential for sexual reproduction.     

Diversity in and evidence for selection on the mitochondrial genome of Phytophthora infestans

Two extant nomenclature systems were reconciled to relate six mitochondrial DNA (mtDNA) haplotypes of Phytophthora infestans, the oomycete pathogen causing late blight disease on potato and tomato. Carter's haplotypes I-a and I-b were included in Goodwin's haplotype A, while Carter's haplotypes II-a and II-b were included in Goodwin's haplotype B. In addition, haplotypes E and F were included in Carter's haplotype I-b. The mutational differences separating the various haplotypes were determined, and we propose that either haplotype I-b(A) or haplotype I-a(A) is the putative ancestral mtDNA of P. infestans, because either can center all the other haplotypes in a logical stepwise network of mutational changes. The occurrence of the six haplotypes in 548 isolates worldwide was determined. Haplotypes I-a and II-a were associated with diverse genotypes worldwide. As previously suggested, haplotype I-b was found only in the US-1 clonal lineage and its variants (n = 99 isolates from 16 countries on 5 continents), and haplotype II-b was limited to the US-6 clonal lineage and its derivatives (n = 36). In a confirmation of a previous suggestion, the randomly mating population in the Toluca Valley of central Mexico (n = 78) was monomorphic for mtDNA haplotype I-a(A). We hypothesize that selection there may be driving the dominance of that single mtDNA haplotype.     

Mating-type heterokaryosis and selfing in Cryphonectria parasitica

Selfing in the chestnut blight fungus, Cryphonectria parasitica, occurs by two different genetic mechanisms. Most self-fertile isolates of C. parasitica are heterokaryotic for mating type, and the progeny from selfing segregate for mating type. Further, we resolved mating-type (MAT) heterokaryons into homokaryons of both mating types by isolating uninucleate asexual spores (conidia). However, because ascospore progeny, with rare exceptions, are not MAT heterokaryons, C. parasitica must lack a regular mechanism to maintain heterokaryosis by selfing. We hypothesize that heterokaryon formation may occur either because of recurrent biparental inbreeding, or by mating-type switching, possibly one involving some kind of parasexual process. The second mechanism found for selfing in C. parasitica occurred less frequently. Three single-conidial isolates (MAT-1 and MAT-2) selfed and produced progeny that did not segregate for mating type. It is currently not known if meiosis occurs during ascospore formation by this mechanism.     

Heterokaryons and parasexual recombinants of Cryphonectria parasitica in two clonal populations in southeastern Europe

Evidence for parasexuality in natural populations of haploid fungi requires the demonstration of diploids or heterokaryons and recombinant genotypes in the absence of sex. We studied clonal populations of the chestnut blight fungus, Cryphonectria parasitica, in southeastern Europe and found evidence of parasexuality in two locations. In Osoj, Macedonia, we found one isolate (Os05-66) that had two alleles at six codominant loci, giving a haplotype that was a composite of two clones in this population. Six single-conidial isolates from Os05-66 had two alleles at some loci, suggesting partial diploidy or aneuploidy, and we found four recombinant haplotypes among single-conidial isolates from hyphal-tip isolates of the same isolate. In Teano, Italy, we found two heterokaryon isolates that were partial composites of two dominant clones. Single-conidial isolates from hyphal-tip isolates had recombinant haplotypes. These results provide evidence that is consistent with the hypothesis of parasexuality in C. parasitica in Europe, similar to an earlier report in a natural population in the USA.     

Significant difference in pathogenicity between MAT1-1 and MAT1-2 isolates in the wheat pathogen Mycosphaerella graminicola

Five Mycosphaerella graminicola populations from four geographic regions (Australia, Israel, Switzerland, and the USA) were assayed for neutral RFLP markers and mating type idiomorphs. On average, 25-30 genetically distinct isolates were selected from each population and their pathogenicity was measured on two wheat cultivars in a common garden experiment conducted in a greenhouse. A significant difference in pathogenicity was found between MAT1-1 and MAT1-2 isolates. On average, MAT1-1 isolates had 14-22% greater pathogenicity than MAT1-2 isolates. The pattern of higher pathogenicity in MAT1-1 isolates was consistent across four geographical populations and on two wheat cultivars. A uniform and continuous variation in pathogenicity was found among isolates within each mating type, but no genetic differentiation in selectively neutral RFLP loci was found between mating types, consistent with the hypothesis that differences in pathogenicity were not due to the effects of specific pathogenicity genes or non-random genetic backgrounds.     

Chickpea Ascochyta Blight: Disease Status and Pathogen Mating Type Distribution in Syria

Chickpea fields were surveyed in nine major chickpea‐growing provinces of Syria in 2008 and 2009 to determine the prevalence and severity of Ascochyta blight, and the distribution of Didymella rabiei mating types (MATs) in the country. A total of 133 Ascochyta rabiei isolates were assayed for mating type, including isolates from older collections that date back to 1982. Multiplex MAT‐specific PCR with three primers was used for MAT analysis. Out of the 133 tested isolates, 64% were MAT1‐1 and 36% were MAT1‐2. Both MATs were found in six provinces but MAT1‐1 alone was found in three provinces. Chi‐squared analysis was used to test for a 1 : 1 ratio of MAT frequencies in all samples. The MAT ratios in the six provinces were not significantly different from 1 : 1, suggesting that there is random mating of the pathogen population under natural conditions. The presence of the two MATs is expected to play a role in the evolution of novel virulence genes that could threaten currently resistant chickpea varieties. Overall analysis of the 133 isolates showed a significant deviation from the 1 : 1 ratio with almost twice as many MAT1‐1 isolates than MAT1‐2 isolates, which indicates a competitive advantage associated with MAT1‐1 in Syria. However, the overall picture of an unequal frequency in MATs indicates that there may be limited sexual recombination occurring in the Syrian population.     

Phenotypic and Genotypic Changes in the Phytophthora infestans Population in Taiwan – 1991 to 2006

Late blight, caused by Phytophthora infestans , is one of the most destructive diseases of tomato in Taiwan. A total of 655 isolates of P. infestans , including 29 isolates from potato, was collected from major tomato and potato production areas of Taiwan during 1991 to 2006. Isolates were characterized for their pathogenicity, mating type, in vitro metalaxyl sensitivity and molecular genotype (including allozyme pattern, mitochondrial genomic haplotype and DNA fingerprint) to monitor population changes in P. infestans . The population of P. infestans in Taiwan underwent a dramatic genetic shift in the 1997–1998 cool growing season. Isolates collected from tomato before 1997 were aggressive to tomato but not potato; most isolates obtained after 1998, were aggressive to both hosts. Metalaxyl sensitivity of isolates changed from sensitive/intermediate before 1997 to resistant since 1998. Similarly, the isolates obtained before 1997 were all US-1 clonal lineage (including variants US-1.1, US-1.2, US-1.3 and US-1.4). During the 1997–1998 cool growing season, the US-11 clonal lineage and the TW-1 genotype appeared, possibly introduced on imported table potatoes. The US-11 lineage spread rapidly and since 1999 has almost completely displaced the old population in Taiwan. Mating type determined by pairing with A1 and A2 reference isolages of P. parasitica , showed all isolates were of the A1 mating type, suggesting that the A2 mating type had not become established in Taiwan. The increasing percentage (up to 42.3% in 2006) of the US-11 variants (including US-11.l, US-11.2, US-11.3 and US-11.4) implied that genomic diversity of the pathogen is changing quickly. Therefore, it is important to continuously monitor the population changes of P. infestans and develop an integrated management strategy for this disease.     

Both mating types of the wheat pathogen Mycosphaerella graminicola are present in Morocco

Septoria tritici blotch caused by the heterothallic ascomycete Mycosphaerella graminicola is one of the most currently damaging diseases on wheat crops worldwide. So far, no information was reported about the status of sexual reproduction of this pathogen under Moroccan conditions. We investigated here for the first time the occurrence of the two mating types (MAT1-1 and MAT1-2) of M. graminicola in Morocco by sampling 141 single-conidial isolates from 4 important wheat producing regions (Gharb, Sa?s, Chaouia and Tadla). The mating type of each isolate was determined by amplification with multiplex PCR of a partial sequence from the corresponding idiomorph. Overall, 43% out of the assessed isolates were MAT1-1 and 57 % were MAT1-2. Both mating types were identified within the 3 sampled regions Gharb, Sa?s and Chaouia, but not in Tadla, where only MAT1-2 isolates were found. The presence of the two mating types highlighted here offers a suitable genetic condition for M. graminicola to occur sexual reproduction in Morocco. The potential of sexual recombination will be examined by the study of mating type frequencies using a large sample size as well as by searching and quantification of pseudothecia in the field.     

Distribution of Mating Types and Genetic Diversity of Ascochyta rabiei Populations in Tunisia Revealed by Mating-type-specific PCR and Random Amplified Polymorphic DNA Markers

The distribution of mating types of Ascochyta rabiei (teleomorph: Didymella rabiei) was determined in Tunisia using a MAT‐specific PCR assay. Among 123 isolates tested, 80% were MAT1‐1 and 20%MAT1‐2. Only MAT1‐1 isolates were present in the Beja and Bizerte regions of Tunisia, whereas both mating types were present in Nabeul, Kef and Jendouba. In the latter three regions, the hypothesis of random mating could not be rejected based on chi‐squared tests of mating‐type ratios (P > 0.05). The lower frequency of the MAT1‐2 coupled with the restricted distribution of this mating type in Tunisia may indicate a recent introduction of MAT1‐2 in Tunisia. This speculation is consistent with the recent (2001) observation of D. rabiei pseudothecia on chickpea debris in Tunisia. Forty isolates representative of the five regions were genetically analysed using 10 random amplified polymorphic DNA (RAPD) primers to provide a preliminary estimate of genetic diversity of the pathogen in Tunisia. Among 129 putative RAPD loci amplified, 81% were polymorphic and 32 unique RAPD fingerprints were detected. A high level of genetic differentiation was detected among subpopulations (G_ST = 0.33). Cluster analyses revealed that isolates from Bizerte, Beja and Jendouba were genetically similar and distinct from isolates sampled in Nabeul and Kef. MAT1‐1 isolates were clustered separately from MAT1‐2 isolates in Jendouba and Nabeul suggesting that recombination may not yet be occurring in these regions despite the occurrence of both mating types in equal frequency in these regions. This lack of recombination between MAT1‐1 and MAT1‐2 also supports the hypothesis of a recent introduction of MAT1‐2 into Tunisia.     

High diversity of vegetative compatibility types in Cryphonectria parasitica in Japan and China

We found high diversity of vegetative compatibility (vc) types in native populations of the chestnut blight fungus, Cryphonectria parasitica, in Japan and China; almost every isolate was in a unique vc type. In Japan we found 71 vc types in a sample of 79 isolates pooled from six populations. Within two populations in China, all isolates (n=28 and 11) had unique vc types; we found 15 vc types among 25 isolates in a third Chinese population where multiple isolates were collected from some trees. None of the isolates from China and only three isolates in the 71 vc types from Japan were compatible with any of 64 vc type testers from Europe, which have known vegetative incompatibility genotypes. To our knowledge this is the first report of vc type diversity for C. parasitica in Japan or of any comparisons of vc types between Asia and Europe. The most significant result of this survey is the identification fungal isolates for expanding knowledge of the genetics of vegetative incompatibility.     

Both mating types in the heterothallic fungus Ophiostoma quercus contain MAT1-1 and MAT1-2 genes

In heterothallic Ascomycota, two opposite but distinct mating types control all sexual processes. Using mating crosses, mating types were assigned to ten isolates of the heterothallic fungal species Ophiostoma quercus. Primers were subsequently designed to target the MAT1-1-1, MAT1-1-3 (of the mating type 1 idiomorph), and MAT1-2-1 (of the mating type 2 idiomorph) genes in these isolates. Results showed that all isolates contained the full gene sequence for the MAT1-2-1 gene. In addition, fragments of the MAT1-1-1 and MAT1-1-3 genes were sequenced from all isolates. These results were unexpected, as each isolate from a heterothallic species would typically contain only one of the two possible MAT idiomorphs.     

Founder events influence structures of Aspergillus flavus populations

In warm regions, agricultural fields are occupied by complex Aspergillus flavus communities composed of isolates in many vegetative compatibility groups (VCGs) with varying abilities to produce highly toxic, carcinogenic aflatoxins. Aflatoxin contamination is reduced with biocontrol products that enable atoxigenic isolates from atoxigenic VCGs to dominate the population. Shifts in VCG frequencies similar to those caused by the introduction of biocontrol isolates were detected in Sonora, Mexico, where biocontrol is not currently practiced. The shifts were attributed to founder events. Although VCGs reproduce clonally, significant diversity exists within VCGs. Simple sequence repeat (SSR) fingerprinting revealed that increased frequencies of VCG YV150 involved a single haplotype. This is consistent with a founder event. Additionally, great diversity was detected among 82 YV150 isolates collected over 20 years across Mexico and the United States. Thirty-six YV150 haplotypes were separated into two populations by Structure and SplitsTree analyses. Sixty-five percent of isolates had MAT1-1 and belonged to one population. The remaining had MAT1-2 and belonged to the second population. SSR alleles varied within populations, but recombination between populations was not detected despite co-occurrence at some locations. Results suggest that YV150 isolates with opposite mating-type have either strongly restrained or lost sexual reproduction among themselves.     

Phylogenetic Analyses on the Diversity of Aspergillus fumigatus Sensu Lato Based on Five Orthologous Loci

One hundred isolates of Aspergillus fumigatus sensu lato mainly from China, as well as from Australia, France, India, Indonesia, Ireland, UK, and USA were analyzed to infer their sequence types (STs) and population diversity based on partial calmodulin, calcineurin regulatory subunit B, beta-tubulin, cytochrome C and calcineurin catalytic subunit A genes as well as their mating types, using ClonalFrame, Structure and MEGA software. Our results inferred 48 STs and showed that most of the STs or lineages evolved independently and without clear population structure among them. Whereas one lineage was recognized that could be a true population and in which one clade might diverge into another distinct lineage, namely, a cryptic species, A. neoellipticus. In addition, we found that mutation, parasexual, and sexual recombination could, respectively, play specific roles in the evolution of these fungi. Our results also showed that MAT1-1/MAT1-2 mating type ratios of A. fumigatus sensu lato was biased to nearly 1:1.4 (20/28) when clone-corrected, but when not clone-corrected, the ratio of MAT1-1/MAT1-2 was so biased as near 1:2 (35/65), which might mean that isolates with MAT1-2 are in the process of losing sexual ability preceding those with MAT1-1.     

Evidence for recombination in Candida glabrata

Despite its clinical importance, little is known of the epidemiology and population structure of Candida glabrata. C. glabrata possesses a mating type system similar to that in Saccharomyces cerevisiae, however mating, meiosis and recombination have not been demonstrated. We performed multilocus sequence typing on a collection of 165 isolates to test for evidence of genetic recombination. A total of 3345 bp from six loci (FKS, LEU2, NMT1, TRP1, UGP1, and URA3) were sequenced for each isolate. The polymorphisms at these loci defined 34 sequence types. Significant evidence for a clonal population was revealed by the index of association and the number of phylogenetically compatible pairs of loci. However, 14 examples of phylogenetic incompatibility were also found. Thus we conclude that although C. glabrata has a predominantly clonal population structure, the multiple phylogenetic incompatibilities found strongly suggest that recombination occurred during the evolution of C. glabrata, and may infrequently still occur.     

Construction of a genetic linkage map of shiitake mushroom Lentinula edodes strain L-54

From fruiting bodies of L. edodes strain L-54, single-spore isolates (SSIs) were collected. Two parental types of L-54 were regenerated via monokaryotization. By means of random-amplified polymorphic DNA (RAPD), DNA samples from L-54, its two parental types, and 32 SSIs were amplified with arbitrary primers. Dedikaryotization was demonstrated, and 91 RAPD-based molecular markers were generated. RAPD markers that were segregated at a 1:1 ratio were used to construct a linkage map of L. edodes. This RAPD-linkage map greatly enhanced the mapping of other inheritable and stable markers [such as those that are linked to a phenotype (the mating type), a known gene (priA) and a sequenced DNA fragment (MAT)] with the aid of mating tests, bulked-segregant analysis, and PCR-single-strand conformational polymorphism. These markers comprised a genetic map of L. edodes with 14 linkage groups and a total length of 622.4 cM.     

Historical and contemporary multilocus population structure of Ascochyta rabiei (teleomorph: Didymella rabiei) in the Pacific Northwest of the United States

The historical and contemporary population genetic structure of the chickpea Ascochyta blight pathogen, Ascochyta rabiei (teleomorph: Didymella rabiei), was determined in the US Pacific Northwest (PNW) using 17 putative AFLP loci, four genetically characterized, sequence-tagged microsatellite loci (STMS) and the mating type locus (MAT). A single multilocus genotype of A. rabiei (MAT1-1) was detected in 1983, which represented the first recorded appearance of Ascochyta blight of chickpea in the PNW. During the following year many additional alleles, including the other mating type allele (MAT1-2), were detected. By 1987, all alleles currently found in the PNW had been introduced. Highly significant genetic differentiation was detected among contemporary subpopulations from different hosts and geographical locations indicating restricted gene flow and/or genetic drift occurring within and among subpopulations and possible selection by host cultivar. Two distinct populations were inferred with high posterior probability which correlated to host of origin and date of sample using Bayesian model-based population structure analyses of multilocus genotypes. Allele frequencies, genotype distributions and population assignment probabilities were significantly different between the historical and contemporary samples of isolates and between isolates sampled from a resistance screening nursery and those sampled from commercial chickpea fields. A random mating model could not be rejected in any subpopulation, indicating the importance of the sexual stage of the fungus both as a source of primary inoculum for Ascochyta blight epidemics and potentially adaptive genotypic diversity.