Microevolution in the pansecondary metabolome of Aspergillus flavus and its potential macroevolutionary implications for filamentous fungi.

Fungi produce a wealth of pharmacologically bioactive secondary metabolites (SMs) from biosynthetic gene clusters (BGCs). It is common practice for drug discovery efforts to treat species' secondary metabolomes as being well represented by a single or a small number of representative genomes. However, this approach misses the possibility that intraspecific population dynamics, such as adaptation to environmental conditions or local microbiomes, may harbor novel BGCs that contribute to the overall niche breadth of species. Using 94 isolates of Aspergillus flavus, a cosmopolitan model fungus, sampled from seven states in the United States, we dereplicate 7,821 BGCs into 92 unique BGCs. We find that more than 25% of pangenomic BGCs show population-specific patterns of presence/absence or protein divergence. Population-specific BGCs make up most of the accessory-genome BGCs, suggesting that different ecological forces that maintain accessory genomes may be partially mediated by population-specific differences in secondary metabolism. We use ultra-high-performance high-resolution mass spectrometry to confirm that these genetic differences in BGCs also result in chemotypic differences in SM production in different populations, which could mediate ecological interactions and be acted on by selection. Thus, our results suggest a paradigm shift that previously unrealized population-level reservoirs of SM diversity may be of significant evolutionary, ecological, and pharmacological importance. Last, we find that several population-specific BGCs from A. flavus are present in Aspergillus parasiticus and Aspergillus minisclerotigenes and discuss how the microevolutionary patterns we uncover inform macroevolutionary inferences and help to align fungal secondary metabolism with existing evolutionary theory.

Genetic Differentiation of Verticillium dahliae Populations Recovered from Symptomatic and Asymptomatic Hosts.

Verticillium dahliae is a soilborne fungal pathogen affecting many economically important crops that can also infect weeds and rotational crops with no apparent disease symptoms. The main research goal was to test the hypothesis that V. dahliae populations recovered from asymptomatic rotational crops and weed species are evolutionarily and genetically distinct from symptomatic hosts. We collected V. dahliae isolates from symptomatic and asymptomatic hosts growing in fields with histories of Verticillium wilt of potato in Israel and Pennsylvania (United States), and used genotyping-by-sequencing to analyze the evolutionary history and genetic differentiation between populations of different hosts. A phylogeny inferred from 26,934 single-nucleotide polymorphisms (SNPs) in 126 V. dahliae isolates displayed a highly clonal structure correlated with vegetative compatibility groups, and isolates grouped in lineages 2A, 2B824, 4A, and 4B, with 77% of the isolates in lineage 4B. The lineages identified in this study were differentiated by host of origin; we found 2A, 2B824, and 4A only in symptomatic hosts but isolates from asymptomatic hosts (weeds, oat, and sorghum) grouped exclusively within lineage 4B, and were genetically indistinguishable from 4B isolates sampled from symptomatic hosts (potato, eggplant, and avocado). Using coalescent analysis of 158 SNPs of lineage 4B, we inferred a genealogy with clades that correlated with geographic origin. In contrast, isolates from asymptomatic and symptomatic hosts shared some of the same haplotypes and were not differentiated. We conclude that asymptomatic weeds and rotational hosts may be potential reservoirs for V. dahliae populations of lineage 4B, which are pathogenic to many cultivated hosts.

Genome Sequence of the Chestnut Blight Fungus Cryphonectria parasitica EP155: A Fundamental Resource for an Archetypical Invasive Plant Pathogen.

Cryphonectria parasitica is the causal agent of chestnut blight, a fungal disease that almost entirely eliminated mature American chestnut from North America over a 50-year period. Here, we formally report the genome of C. parasitica EP155 using a Sanger shotgun sequencing approach. After finishing and integration with simple-sequence repeat markers, the assembly was 43.8 Mb in 26 scaffolds (L50 = 5; N50 = 4.0Mb). Eight chromosomes are predicted: five scaffolds have two telomeres and six scaffolds have one telomere sequence. In total, 11,609 gene models were predicted, of which 85% show similarities to other proteins. This genome resource has already increased the utility of a fundamental plant pathogen experimental system through new understanding of the fungal vegetative incompatibility system, with significant implications for enhancing mycovirus-based biological control.

Population Subdivision and the Frequency of Aflatoxigenic Isolates in Aspergillus flavus in the United States.

Consumption of food contaminated with aflatoxin, from crops infected by Aspergillus flavus, is associated with acute toxicosis, cancer, and stunted growth. Although such contamination is more common in the lower latitudes of the United States, it is unclear whether this pattern is associated with differences in the relative frequencies of aflatoxigenic individuals of A. flavus. To determine whether the frequency of the aflatoxin-producing ability of A. flavus increases as latitude decreases, we sampled 281 isolates from field soils in two north-south transects in the United States and tested them for aflatoxin production. We also genotyped 161 isolates using 10 microsatellite markers to assess population structure. Although the population density of A. flavus was highest at lower latitudes, there was no difference in the frequency of aflatoxigenic A. flavus isolates in relation to latitude. We found that the U.S. population of A. flavus is subdivided into two genetically differentiated subpopulations that are not associated with the chemotype or geographic origin of the isolates. The two populations differ markedly in allelic and genotypic diversity. The less diverse population is more abundant and may represent a clonal lineage derived from the more diverse population. Overall, increased aflatoxin contamination in lower latitudes may be explained partially by differences in the population density of A. flavus, not genetic population structure.

Balancing selection at nonself recognition loci in the chestnut blight fungus, Cryphonectria parasitica, demonstrated by trans-species polymorphisms, positive selection, and even allele frequencies.

Balancing selection has been inferred in diverse organisms for nonself recognition genes, including those involved in immunity, mating compatibility, and vegetative incompatibility. Although selective forces maintaining polymorphisms are known for genes involved in immunity and mating, mechanisms of balancing selection for vegetative incompatibility genes in fungi are being debated. We hypothesized that allorecognition and its consequent inhibition of virus transmission contribute to the maintenance of polymorphisms in vegetative incompatibility loci (vic) in the chestnut blight fungus, Cryphonectria parasitica. Balancing selection was demonstrated at two loci, vic2 and vic6, by trans-species polymorphisms in C. parasitica, C. radicalis, and C. japonica and signatures of positive selection in gene sequences. In addition, more than half (31 of 54) of allele frequency estimates at six vic loci in nine field populations of C. parasitica from Asia and the eastern US were not significantly different from 0.5, as expected at equilibrium for two alleles per locus under balancing selection. At three vic loci, deviations from 0.5 were predicted based on the effects of heteroallelism on virus transmission. Twenty-five of 27 allele frequency estimates were greater than or equal to 0.5 for the allele that confers significantly stronger inhibition of virus transmission at three loci with asymmetric transmission. These results are consistent with the allorecognition hypothesis that vegetative incompatibility genes are under selection because of their role in reducing infection by viruses.

Population Genetics of Verticillium dahliae in Iran Based on Microsatellite and Single Nucleotide Polymorphism Markers.

Verticillium dahliae is a plant pathogenic fungus that reproduces asexually and its population structure is highly clonal. In the present study, 78 V. dahliae isolates from Iran were genotyped for mating type, single nucleotide polymorphisms (SNPs), and microsatellites to assign them to clonal lineages and to determine population genetic structure in Iran. The mating type of all isolates was MAT1-2. Based on neighbor-joining analysis and minimum spanning networks constructed from SNPs and microsatellite genotypes, respectively, all but four isolates were assigned to lineage 2B824; four isolates were assigned to lineage 4B. The inferred coalescent genealogy of isolates in lineage 2B824 showed a clear divergence into two clades that corresponded to geographic origin and host. Haplotypes of cotton and pistachio isolates sampled from central Iran were in one clade, and those of isolates from Prunus spp. sampled from northwestern Iran were in the other. The strong divergence in haplotypes between the two clades suggests that there were at least two separate introductions of lineage 2B824 to different parts of Iran. Given the history of cotton and pistachio cultivation and Verticillium wilt in Iran, these results are consistent with the hypothesis that cotton was historically a likely source inoculum causing Verticillium wilt in pistachio.

Balancing selection for aflatoxin in Aspergillus flavus is maintained through interference competition with, and fungivory by insects.

The role of microbial secondary metabolites in the ecology of the organisms that produce them remains poorly understood. Variation in aflatoxin production by Aspergillus flavus is maintained by balancing selection, but the ecological function and impact on fungal fitness of this compound are unknown. We hypothesize that balancing selection for aflatoxin production in A. flavus is driven by interaction with insects. To test this, we competed naturally occurring aflatoxigenic and non-aflatoxigenic fungal isolates against Drosophila larvae on medium containing 0-1750 ppb aflatoxin, using quantitative PCR to quantify A. flavus DNA as a proxy for fungal fitness. The addition of aflatoxin across this range resulted in a 26-fold increase in fungal fitness. With no added toxin, aflatoxigenic isolates caused higher mortality of Drosophila larvae and had slightly higher fitness than non-aflatoxigenic isolates. Additionally, aflatoxin production increased an average of 1.5-fold in the presence of a single larva and nearly threefold when the fungus was mechanically damaged. We argue that the role of aflatoxin in protection from fungivory is inextricably linked to its role in interference competition. Our results, to our knowledge, provide the first clear evidence of a fitness advantage conferred to A. flavus by aflatoxin when interacting with insects.

Aphid vector population density determines the emergence of necrogenic satellite RNAs in populations of cucumber mosaic virus.

The satellite RNAs of cucumber mosaic virus (CMV) that induce systemic necrosis in tomato plants (N-satRNA) multiply to high levels in the infected host while severely depressing CMV accumulation and, hence, its aphid transmission efficiency. As N-satRNAs are transmitted into CMV particles, the conditions for N-satRNA emergence are not obvious. Model analyses with realistic parameter values have predicted that N-satRNAs would invade CMV populations only when transmission rates are high. Here, we tested this hypothesis experimentally by passaging CMV or CMV+N-satRNAs at low or high aphid densities (2 or 8 aphids/plant). As predicted, high aphid densities were required for N-satRNA emergence. The results showed that at low aphid densities, random effects due to population bottlenecks during transmission dominate the epidemiological dynamics of CMV/CMV+N-satRNA. The results suggest that maintaining aphid populations at low density will prevent the emergence of highly virulent CMV+N-satRNA isolates.

Vertical transmission selects for reduced virulence in a plant virus and for increased resistance in the host.

For the last three decades, evolutionary biologists have sought to understand which factors modulate the evolution of parasite virulence. Although theory has identified several of these modulators, their effect has seldom been analysed experimentally. We investigated the role of two such major factors-the mode of transmission, and host adaptation in response to parasite evolution-in the evolution of virulence of the plant virus Cucumber mosaic virus (CMV) in its natural host Arabidopsis thaliana. To do so, we serially passaged three CMV strains under strict vertical and strict horizontal transmission, alternating both modes of transmission. We quantified seed (vertical) transmission rate, virus accumulation, effect on plant growth and virulence of evolved and non-evolved viruses in the original plants and in plants derived after five passages of vertical transmission. Our results indicated that vertical passaging led to adaptation of the virus to greater vertical transmission, which was associated with reductions of virus accumulation and virulence. On the other hand, horizontal serial passages did not significantly modify virus accumulation and virulence. The observed increases in CMV seed transmission, and reductions in virus accumulation and virulence in vertically passaged viruses were due also to reciprocal host adaptation during vertical passages, which additionally reduced virulence and multiplication of vertically passaged viruses. This result is consistent with plant-virus co-evolution. Host adaptation to vertically passaged viruses was traded-off against reduced resistance to the non-evolved viruses. Thus, we provide evidence of the key role that the interplay between mode of transmission and host-parasite co-evolution has in determining the evolution of virulence.

Clonal Expansion and Migration of a Highly Virulent, Defoliating Lineage of Verticillium dahliae.

We used a population genomics approach to test the hypothesis of clonal expansion of a highly fit genotype in populations of Verticillium dahliae. This fungal pathogen has a broad host range and can be dispersed in contaminated seed or other plant material. It has a highly clonal population structure, with several lineages having nearly worldwide distributions in agricultural crops. Isolates in lineage 1A are highly virulent and cause defoliation in cotton, okra, and olive (denoted 1A/D), whereas those in other lineages cause wilting but not defoliation (ND). We tested whether the highly virulent lineage 1A/D could have spread from the southwestern United States to the Mediterranean basin, as predicted from historical records. We found 187 single-nucleotide polymorphisms (SNPs), determined by genotyping by sequencing, among 91 isolates of lineage 1A/D and 5 isolates in the closely related lineage 1B/ND. Neighbor-joining and maximum-likelihood analyses on the 187 SNPs showed a clear divergence between 1A/D and 1B/ND haplotypes. Data for only 77 SNPs were obtained for all 96 isolates (no missing data); lineages 1A/D and 1B/ND differed by 27 of these 77 SNPs, confirming a clear divergence between the two lineages. No evidence of recombination was detected within or between these two lineages. Phylogenetic and genealogical analyses resulted in five distinct subclades of 1A/D isolates that correlated closely with geographic origins in the Mediterranean basin, consistent with the hypothesis that the D pathotype was introduced at least five times in independent founder events into this region from a relatively diverse source population. The inferred ancestral haplotype was found in two isolates sampled before 1983 from the southwestern United States, which is consistent with historical records that 1A/D originated in North America. The five subclades coalesce with the ancestral haplotype at the same time, consistent with a hypothesis of rapid population expansion in the source population during the emergence of 1A/D as a severe pathogen of cotton in the United States.

Mechanisms of Resistance to an Azole Fungicide in the Grapevine Powdery Mildew Fungus, Erysiphe necator.

We studied the mechanisms of azole resistance in Erysiphe necator by quantifying the sensitivity to myclobutanil (EC50) in 65 isolates from the eastern United States and 12 from Chile. From each isolate, we sequenced the gene for sterol 14α-demethylase (CYP51), and measured the expression of CYP51 and homologs of four putative efflux transporter genes, which we identified in the E. necator transcriptome. Sequence variation in CYP51 was relatively low, with sequences of 40 U.S. isolates identical to the reference sequence. Nine U.S. isolates and five from Chile carried a previously identified A to T nucleotide substitution in position 495 (A495T), which results in an amino acid substitution in codon 136 (Y136F) and correlates with high levels of azole resistance. We also found a nucleotide substitution in position 1119 (A1119C) in 15 U.S. isolates, whose mean EC50 value was equivalent to that for the Y136F isolates. Isolates carrying mutation A1119C had significantly greater CYP51 expression, even though A1119C does not affect the CYP51 amino acid sequence. Regression analysis showed no significant effects of the expression of efflux transporter genes on EC50. Both the Y136F mutation in CYP51 and increased CYP51 expression appear responsible for azole resistance in eastern U.S. populations of E. necator.

Linkage disequilibrium and spatial aggregation of genotypes in sexually reproducing populations of Erysiphe necator.

Random mating and recombination in heterothallic ascomycetes should result in high genotypic diversity, 1:1 mating-type ratios, and random associations of alleles, or linkage equilibrium, at different loci. To test for random mating in populations of the grape powdery mildew fungus Erysiphe necator, we sampled isolates from vineyards of Vitis vinifera in Burdett, NY (NY09) and Winchester, VA (VA09) at the end of the epidemic in fall 2009. We also sampled isolates from the same Winchester, VA vineyard in spring 2010 at the onset of the next epidemic. Isolates were genotyped for mating type and 11 microsatellite markers. In the spring sample, which originated from ascospore infections, nearly every isolate had a unique genotype. In contrast, fall populations were less diverse. In all, 9 of 45 total genotypes in VA09 were represented by two or more isolates; 3 of 40 total genotypes in NY09 were represented by two or more isolates, with 1 genotype represented by 20 isolates. After clone correction, mating-type ratios in the three populations did not deviate from 1:1. However, even with clone correction, we detected significant linkage disequilibrium (LD) in all populations. Mantel tests detected positive correlations between genetic and physical distances within vineyards. Spatial autocorrelation showed aggregations up to 42 and 3 m in VA09 and NY09, respectively. Spatial autocorrelation most likely results from short dispersal distances. Overall, these results suggest that spatial genetic aggregation and clonal genotypes that arise during the asexual phase of the epidemic contribute to persistent LD even though populations undergo sexual reproduction annually.

Identification of race-specific resistance in North American Vitis spp. limiting Erysiphe necator hyphal growth.

Race-specific resistance against powdery mildews is well documented in small grains but, in other crops such as grapevine, controlled analysis of host-pathogen interactions on resistant plants is uncommon. In the current study, we attempted to confirm powdery mildew resistance phenotypes through vineyard, greenhouse, and in vitro inoculations for test cross-mapping populations for two resistance sources: (i) a complex hybrid breeding line, 'Bloodworth 81-107-11', of at least Vitis rotundifolia, V. vinifera, V. berlandieri, V. rupestris, V. labrusca, and V. aestivalis background; and (ii) Vitis hybrid 'Tamiami' of V. aestivalis and V. vinifera origin. Statistical analysis of vineyard resistance data suggested the segregation of two and three race-specific resistance genes from the two sources, respectively. However, in each population, some resistant progeny were susceptible in greenhouse or in vitro screens, which suggested the presence of Erysiphe necator isolates virulent on progeny segregating for one or more resistance genes. Controlled inoculation of resistant and susceptible progeny with a diverse set of E. necator isolates clearly demonstrated the presence of fungal races differentially interacting with race-specific resistance genes, providing proof of race specificity in the grape powdery mildew pathosystem. Consistent with known race-specific resistance mechanisms, both resistance sources were characterized by programmed cell death of host epidermal cells under appressoria, which arrested or slowed hyphal growth; this response was also accompanied by collapse of conidia, germ tubes, appressoria, and secondary hyphae. The observation of prevalent isolates virulent on progeny with multiple race-specific resistance genes before resistance gene deployment has implications for grape breeding strategies. We suggest that grape breeders should characterize the mechanisms of resistance and pyramid multiple resistance genes with different mechanisms for improved durability.

Identification and structure of the mating-type locus and development of PCR-based markers for mating type in powdery mildew fungi.

In ascomycetes, mating compatibility is regulated by the mating-type locus, MAT1. The objectives of this study were to identify and sequence genes at the MAT1 locus in the grape powdery mildew fungus, Erysiphe necator, to develop a PCR-based marker for determining mating type in E. necator, and to develop degenerate primers for amplification by PCR of conserved regions of mating-type idiomorphs in other powdery mildew fungi. We identified MAT1-2-1 of the MAT1-2 idiomorph in E. necator based on the homologous sequence in the genome of Blumeria graminis f. sp. hordei and we found MAT1-1-1 and MAT1-1-3 of the MAT1-1 idiomorph from transcriptome sequences of E. necator. We developed and applied a reliable PCR-based multiplex marker to confirm that genotype correlated with mating phenotype, which was determined by pairing with mating-type tester isolates. Additionally, we used the marker to genotype populations of E. necator from different Vitis spp. from throughout the USA. We found both mating types were present in all populations and mating-type ratios did not deviate from 1:1. The mating-type genes in E. necator are similar to those of other Leotiomycetes; however, the structure of the MAT1 locus in E. necator, like the MAT1-2 idiomorph of B. graminis, is markedly different from other ascomycetes in that it is greatly expanded and may contain a large amount of repetitive DNA. As a result, we were unable to amplify and sequence either idiomorph in its entirety. We designed degenerate primers that amplify conserved regions of MAT1-1 and MAT1-2 in E. necator, Podosphaera xanthii, Microsphaera syringae, and B. graminis, representing the major clades of the Erysiphales. These degenerate primers or sequences obtained in this study from these species can be used to identify and sequence MAT1 genes or design mating-type markers in other powdery mildew fungi as well.

Differential gene expression during conidiation in the grape powdery mildew pathogen, Erysiphe necator.

Asexual sporulation (conidiation) is coordinately regulated in the grape powdery mildew pathogen Erysiphe necator but nothing is known about its genetic regulation. We hypothesized that genes required for conidiation in other fungi would be upregulated at conidiophore initiation or full conidiation (relative to preconidiation vegetative growth and development of mature ascocarps), and that the obligate biotrophic lifestyle of E. necator would necessitate some novel gene regulation. cDNA amplified fragment length polymorphism analysis with 45 selective primer combinations produced ≈1,600 transcript-derived fragments (TDFs), of which 620 (39%) showed differential expression. TDF sequences were annotated using BLAST analysis of GenBank and of a reference transcriptome for E. necator developed by 454-FLX pyrosequencing of a normalized cDNA library. One-fourth of the differentially expressed, annotated sequences had similarity to fungal genes of unknown function. The remaining genes had annotated function in metabolism, signaling, transcription, transport, and protein fate. As expected, a portion of orthologs known in other fungi to be involved in developmental regulation was upregulated immediately prior to or during conidiation; particularly noteworthy were several genes associated with the light-dependent VeA regulatory system, G-protein signaling (Pth11 and a kelch repeat), and nuclear transport (importin-ß and Ran). This work represents the first investigation into differential gene expression during morphogenesis in E. necator and identifies candidate genes and hypotheses for characterization in powdery mildews. Our results indicate that, although control of conidiation in powdery mildews may share some basic elements with established systems, there are significant points of divergence as well, perhaps due, in part, to the obligate biotrophic lifestyle of powdery mildews.

Phylogeography and population structure of the grape powdery mildew fungus, Erysiphe necator, from diverse Vitis species.

BACKGROUND: The grape powdery mildew fungus, Erysiphe necator, was introduced into Europe more than 160 years ago and is now distributed everywhere that grapes are grown. To understand the invasion history of this pathogen we investigated the evolutionary relationships between introduced populations of Europe, Australia and the western United States (US) and populations in the eastern US, where E. necator is thought to be native. Additionally, we tested the hypothesis that populations of E. necator in the eastern US are structured based on geography and Vitis host species. RESULTS: We sequenced three nuclear gene regions covering 1803 nucleotides from 146 isolates of E. necator collected from the eastern US, Europe, Australia, and the western US. Phylogeographic analyses show that the two genetic groups in Europe represent two separate introductions and that the genetic groups may be derived from eastern US ancestors. Populations from the western US and Europe share haplotypes, suggesting that the western US population was introduced from Europe. Populations in Australia are derived from European populations. Haplotype richness and nucleotide diversity were significantly greater in the eastern US populations than in the introduced populations. Populations within the eastern US are geographically differentiated; however, no structure was detected with respect to host habitat (i.e., wild or cultivated). Populations from muscadine grapes, V. rotundifolia, are genetically distinct from populations from other Vitis host species, yet no differentiation was detected among populations from other Vitis species. CONCLUSIONS: Multilocus sequencing analysis of the grape powdery mildew fungus is consistent with the hypothesis that populations in Europe, Australia and the western US are derived from two separate introductions and their ancestors were likely from native populations in the eastern US. The invasion history of E. necator follows a pattern consistent with plant-mediated dispersal, however, more exhaustive sampling is required to make more precise conclusions as to origin. E. necator shows no genetic structure across Vitis host species, except with respect to V. rotundifolia.

Variation in pathogenicity and aggressiveness of Erysiphe necator from different Vitis spp. and geographic origins in the eastern United States.

Eastern North America is considered the center of diversity for many Vitis spp. and for the grape powdery mildew pathogen, Erysiphe necator. However, little is known about populations of E. necator from wild Vitis spp. We determined the phenotypic variation in pathogenicity and aggressiveness of E. necator among isolates from wild and domesticated Vitis spp. from diverse geographic regions in the eastern United States. To test pathogenicity, we inoculated 38 E. necator isolates on three wild Vitis spp., two commercially grown hybrids and the European wine grape, Vitis vinifera. V. rotundifolia (muscadine grape) was the only host species on which complete host specialization was evident; it was only susceptible to isolates collected from V. rotundifolia. All isolates, regardless of source host, were pathogenic on the other Vitis spp. We found no differences in components of aggressiveness latent period and lesion size among isolates from different source hosts when inoculated on V. vinifera, which is highly susceptible to powdery mildew. However significant variation was evident among isolates on the more resistant V. labruscana 'Niagara'. Isolates from the wild species V. aestivalis were the most aggressive, whereas isolates from V. vinifera were not more aggressive than isolates from other source hosts. Greater aggressiveness was also detected among isolates from the southeastern United States compared with isolates from the northeastern United States.

The Frequency of Sex: Population Genomics Reveals Differences in Recombination and Population Structure of the Aflatoxin-Producing Fungus Aspergillus flavus.

The apparent rarity of sex in many fungal species has raised questions about how much sex is needed to purge deleterious mutations and how differences in frequency of sex impact fungal evolution. We sought to determine how differences in the extent of recombination between populations of Aspergillus flavus impact the evolution of genes associated with the synthesis of aflatoxin, a notoriously potent carcinogen. We sequenced the genomes of, and quantified aflatoxin production in, 94 isolates of A. flavus sampled from seven states in eastern and central latitudinal transects of the United States. The overall population is subdivided into three genetically differentiated populations (A, B, and C) that differ greatly in their extent of recombination, diversity, and aflatoxin-producing ability. Estimates of the number of recombination events and linkage disequilibrium decay suggest relatively frequent sex only in population A. Population B is sympatric with population A but produces significantly less aflatoxin and is the only population where the inability of nonaflatoxigenic isolates to produce aflatoxin was explained by multiple gene deletions. Population expansion evident in population B suggests a recent introduction or range expansion. Population C is largely nonaflatoxigenic and restricted mainly to northern sampling locations through restricted migration and/or selection. Despite differences in the number and type of mutations in the aflatoxin gene cluster, codon optimization and site frequency differences in synonymous and nonsynonymous mutations suggest that low levels of recombination in some A. flavus populations are sufficient to purge deleterious mutations.IMPORTANCE Differences in the relative frequencies of sexual and asexual reproduction have profound implications for the accumulation of deleterious mutations (Muller's ratchet), but little is known about how these differences impact the evolution of ecologically important phenotypes. Aspergillus flavus is the main producer of aflatoxin, a notoriously potent carcinogen that often contaminates food. We investigated if differences in the levels of production of aflatoxin by A. flavus could be explained by the accumulation of deleterious mutations due to a lack of recombination. Despite differences in the extent of recombination, variation in aflatoxin production is better explained by the demography and history of specific populations and may suggest important differences in the ecological roles of aflatoxin among populations. Furthermore, the association of aflatoxin production and populations provides a means of predicting the risk of aflatoxin contamination by determining the frequencies of isolates from low- and high-production populations.

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.

Speciation of a tropical fungal species pair following transoceanic dispersal.

This study focuses on a pair of fungal species, Moelleriella libera and M. raciborskii (Ascomycota: Clavicipitaceae) from the neotropics and paleotropics, respectively, that are phenotypically nearly indistinguishable. Molecular analyses based on DNA sequences from RNA polymerase II subunit 2 (RPB2), translation elongation factor 1-alpha (EF1-alpha) and beta-tubulin genes confirm that they are recently derived sister species. Speciation appears to have followed an historical transoceanic dispersal event. Models of population structure and migration from TCS, IM, and coalescent-based analyses suggest there is little gene flow between the two species. The direction of dispersal, investigated using the progression rule and coalescent-based gene genealogies, was likely from the New World to the Old World.