Mycena s.s. is a ubiquitous mushroom genus whose members degrade multiple dead plant substrates and opportunistically invade living plant roots. Having sequenced the nuclear genomes of 24 Mycena species, we find them to defy the expected patterns for fungi based on both their traditionally perceived saprotrophic ecology and substrate specializations. Mycena displayed massive genome expansions overall affecting all gene families, driven by novel gene family emergence, gene duplications, enlarged secretomes encoding polysaccharide degradation enzymes, transposable element (TE) proliferation, and horizontal gene transfers. Mainly due to TE proliferation, Arctic Mycena species display genomes of up to 502 Mbp (2-8× the temperate Mycena), the largest among mushroom-forming Agaricomycetes, indicating a possible evolutionary convergence to genomic expansions sometimes seen in Arctic plants. Overall, Mycena show highly unusual, varied mosaic-like genomic structures adaptable to multiple lifestyles, providing genomic illustration for the growing realization that fungal niche adaptations can be far more fluid than traditionally believed.

BACKGROUND: Colletotrichum fungi infect a wide diversity of monocot and dicot hosts, causing diseases on almost all economically important plants worldwide. Colletotrichum is also a suitable model for studying gene family evolution on a fine scale to uncover events in the genome associated with biological changes.
RESULTS: Here we present the genome sequences of 30 Colletotrichum species covering the diversity within the genus. Evolutionary analyses revealed that the Colletotrichum ancestor diverged in the late Cretaceous in parallel with the diversification of flowering plants. We provide evidence of independent host jumps from dicots to monocots during the evolution of Colletotrichum, coinciding with a progressive shrinking of the plant cell wall degradative arsenal and expansions in lineage-specific gene families. Comparative transcriptomics of 4 species adapted to different hosts revealed similarity in gene content but high diversity in the modulation of their transcription profiles on different plant substrates. Combining genomics and transcriptomics, we identified a set of core genes such as specific transcription factors, putatively involved in plant cell wall degradation.
CONCLUSIONS: These results indicate that the ancestral Colletotrichum were associated with dicot plants and certain branches progressively adapted to different monocot hosts, reshaping the gene content and its regulation.

A. fumigatus is a deadly fungal pathogen, responsible for >400,000 infections/year and high mortality rates. A. fumigatus strains exhibit variation in infection-relevant traits, including in their virulence. However, most A. fumigatus protein-coding genes, including those that modulate its virulence, are shared between A. fumigatus strains and closely related non-pathogenic relatives. We hypothesized that A. fumigatus genes exhibit substantial genetic variation in the non-coding regions immediately upstream to the start codons of genes, which could reflect differences in gene regulation between strains. To begin testing this hypothesis, we identified 5,812 single-copy orthologs across the genomes of 263 A. fumigatus strains. In general, A. fumigatus non-coding regions showed higher levels of sequence variation compared to their corresponding protein-coding regions. Focusing on 2,482 genes whose protein-coding sequence identity scores ranged between 75% and 99%, we identified 478 total genes with signatures of positive selection only in their non-coding regions and 65 total genes with signatures only in their protein-coding regions. 28 of the 478 non-coding regions and 5 of the 65 protein-coding regions under selection are associated with genes known to modulate A. fumigatus virulence. Non-coding region variation between A. fumigatus strains included single nucleotide polymorphisms and insertions or deletions of at least a few nucleotides. These results show that non-coding regions of A. fumigatus genes harbor greater sequence variation than protein-coding regions, raising the hypothesis that this variation may contribute to A. fumigatus phenotypic heterogeneity.

We sequenced and annotated the genomes of the ascomycete fungi Trichoderma harzianum, Trichoderma aggressivum f. europaeum, and Purpureocillium lilacinum. Moreover, we developed a website to allow users to interactively analyze the assemblies, gene predictions, and functional annotations of these species and 70+ previously sequenced fungi.

A.fumigatus is a deadly fungal pathogen, responsible for >400,000 infections/year and high mortality rates. A. fumigatus strains exhibit variation in infection-relevant traits, including in their virulence. However, most A. fumigatus protein-coding genes, including those that modulate its virulence, are shared between A. fumigatus strains and closely related non-pathogenic relatives. We hypothesized that A. fumigatus genes exhibit substantial genetic variation in the non-coding regions immediately upstream to the start codons of genes, which could reflect differences in gene regulation between strains. To begin testing this hypothesis, we identified 5,812 single-copy orthologs across the genomes of 263 A. fumigatus strains. A. fumigatus non-coding regions showed higher levels of sequence variation compared to their corresponding protein-coding regions. Specifically, we found that 1,274 non-coding regions exhibited <75% nucleotide sequence similarity (compared to 928 protein-coding regions) and 3,721 non-coding regions exhibited between 75% and 99% similarity (compared to 2,482 protein-coding regions) across strains. Only 817 non-coding regions exhibited ≥99% sequence similarity compared to 2,402 protein-coding regions. By examining 2,482 genes whose protein-coding sequence identity scores ranged between 75% and 99%, we identified 478 total genes with signatures of positive selection only in their non-coding regions and 65 total genes with signatures only in their protein-coding regions. 28 of the 478 non-coding regions and 5 of the 65 protein-coding regions under selection are associated with genes known to modulate A. fumigatus virulence. Non-coding region variation between A. fumigatus strains included single nucleotide polymorphisms and insertions or deletions of at least a few nucleotides. These results show that non-coding regions of A. fumigatus genes harbor greater sequence variation than protein-coding regions, raising the hypothesis that this variation may contribute to A. fumigatus phenotypic heterogeneity.

BACKGROUND: Filamentous fungi are prolific producers of bioactive molecules and enzymes with important applications in industry. Yet, the vast majority of fungal species remain undiscovered or uncharacterized. Here we focus our attention to a wild fungal isolate that we identified as Anthostomella pinea. The fungus belongs to a complex polyphyletic genus in the family of Xylariaceae, which is known to comprise endophytic and pathogenic fungi that produce a plethora of interesting secondary metabolites. Despite that, Anthostomella is largely understudied and only two species have been fully sequenced and characterized at a genomic level.
RESULTS: In this work, we used long-read sequencing to obtain the complete 53.7 Mb genome sequence including the full mitochondrial DNA. We performed extensive structural and functional annotation of coding sequences, including genes encoding enzymes with potential applications in biotechnology. Among others, we found that the genome of A. pinea encodes 91 biosynthetic gene clusters, more than 600 CAZymes, and 164 P450s. Furthermore, untargeted metabolomics and molecular networking analysis of the cultivation extracts revealed a rich secondary metabolism, and in particular an abundance of sesquiterpenoids and sesquiterpene lactones. We also identified the polyketide antibiotic xanthoepocin, to which we attribute the anti-Gram-positive effect of the extracts that we observed in antibacterial plate assays.
CONCLUSIONS: Taken together, our results provide a first glimpse into the potential of Anthstomella pinea to provide new bioactive molecules and biocatalysts and will facilitate future research into these valuable metabolites.

Fungi that are edible or fermentative were domesticated through selective cultivation of their desired traits. Domestication is often associated with inbreeding or selfing, which may fix traits other than those under selection, and causes an overall decrease in heterozygosity. A hallucinogenic mushroom, Psilocybe cubensis, was domesticated from its niche in livestock dung for production of psilocybin. It has caused accidental poisonings since the 1940s in Australia, which is a population hypothesized to be introduced from an unknown center of origin. We sequenced genomes of 38 isolates from Australia and compared them with 86 genomes of commercially available cultivars to determine (1) whether P. cubensis was introduced to Australia, and (2) how domestication has impacted commercial cultivars. Our analyses of genome-wide SNPs and single-copy orthologs showed that the Australian population is naturalized, having recovered its effective population size after a bottleneck when it was introduced, and it has maintained relatively high genetic diversity based on measures of nucleotide and allelic diversity. In contrast, domesticated cultivars generally have low effective population sizes and hallmarks of selfing and clonal propagation, including low genetic diversity, low heterozygosity, high linkage disequilibrium, and low allelic diversity of mating-compatibility genes. Analyses of kinship show that most cultivars are founded from related populations. Alleles in the psilocybin gene cluster are identical across most cultivars of P. cubensis with low diversity across coding sequence; however, unique allelic diversity in Australia and some cultivars may translate to differences in biosynthesis of psilocybin and its analogs.

Candida albicans is an opportunistic fungal pathogen of humans, yet the within-host dynamics of C. albicans infection are not clear. While C. albicans is commonly diploid, it exhibits a range of ploidies, including tetraploidy. Previous work found that tetraploid C. albicans populations exhibited rapid adaptation and significant genome instability when evolved in vitro. Host immune function alters the rate and magnitude of C. albicans virulence evolution, but the effects of the host immunity on tetraploid C. albicans populations are unclear. Here, we tested the effects of the host immunity on genome stability and virulence evolution of tetraploid C. albicans using experimental evolution. We selected for C. albicans increased virulence within either immunocompetent or immunocompromised Caenorhabditis elegans hosts. After nine passages we observed a response to selection for increased virulence. Both populations exposed to either immunocompetent or immunocompromised hosts increased virulence after passage through C. elegans hosts. However, the C. albicans populations passaged through immunocompetent hosts under selection exhibited unique temporal dynamics, a rapid increase in virulence and then subsequent loss of virulence. Most C. albicans populations exhibited genome size reduction within six passages, however populations exposed to immunocompetent hosts exhibited the most rapid transition to ~diploid. Therefore, we found that tetraploids rapidly increase in virulence and decrease genome size within host environments. Further, the combination of selection for greater virulence in the presence of immunocompetent hosts results in major virulence fluctuations and genome size changes. Thus, host immunity significantly impacts the evolutionary trajectories of tetraploid C. albicans.

Understanding the origins of variation in agricultural pathogens is of fundamental interest and practical importance, especially for diseases that threaten food security. Fusarium oxysporum is among the most important of soil-borne pathogens, with a global distribution and an extensive host range. The pathogen is considered to be asexual, with horizontal transfer of chromosomes providing an analog of assortment by meiotic recombination. Here, we challenge those assumptions based on the results of population genomic analyses, describing the pathogen\'s diversity and inferring its origins and functional consequences in the context of a single, long-standing agricultural system. We identify simultaneously low nucleotide distance among strains, and unexpectedly high levels of genetic and genomic variability. We determine that these features arise from a combination of genome-scale recombination, best explained by widespread sexual reproduction, and presence-absence variation consistent with chromosomal rearrangement. Pangenome analyses document an accessory genome more than twice the size of the core genome, with contrasting evolutionary dynamics. The core genome is stable, with low diversity and high genetic differentiation across geographic space, while the accessory genome is paradoxically more diverse and unstable but with lower genetic differentiation and hallmarks of contemporary gene flow at local scales. We suggest a model in which episodic sexual reproduction generates haplotypes that are selected and then maintained through clone-like dynamics, followed by contemporary genomic rearrangements that reassort the accessory genome among sympatric strains. Taken together, these processes contribute unique genome content, including reassortment of virulence determinants that may explain observed variation in pathogenic potential.

UNASSIGNED: Nosema is a diverse genus of unicellular microsporidian parasites of insects and other arthropods. Nosema muscidifuracis infects parasitoid wasp species of Muscidifurax zaraptor and M. raptor (Hymenoptera: Pteromalidae), causing ~50% reduction in longevity and ~90% reduction in fecundity.
UNASSIGNED: Here, we report the first assembly of the N. muscidifuracis genome (14,397,169 bp in 28 contigs) of high continuity (contig N50 544.3 Kb) and completeness (BUSCO score 97.0%). A total of 2,782 protein-coding genes were annotated, with 66.2% of the genes having two copies and 24.0% of genes having three copies. These duplicated genes are highly similar, with a sequence identity of 99.3%. The complex pattern suggests extensive gene duplications and rearrangements across the genome. We annotated 57 rDNA loci, which are highly GC-rich (37%) in a GC-poor genome (25% genome average). Nosema-specific qPCR primer sets were designed based on 18S rDNA annotation as a diagnostic tool to determine its titer in host samples. We discovered high Nosema titers in Nosema-cured M. raptor and M. zaraptor using heat treatment in 2017 and 2019, suggesting that the remedy did not completely eliminate the Nosema infection. Cytogenetic analyses revealed heavy infections of N. muscidifuracis within the ovaries of M. raptor and M. zaraptor, consistent with the titer determined by qPCR and suggesting a heritable component of infection and per ovum vertical transmission.
UNASSIGNED: The parasitoids-Nosema system is laboratory tractable and, therefore, can serve as a model to inform future genome manipulations of Nosema-host system for investigations of Nosemosis.