Genes in prokaryotic genomes are often arranged into clusters and co-transcribed into polycistronic RNAs. Isolated examples of polycistronic RNAs were also reported in some higher eukaryotes but their presence was generally considered rare. Here we developed a long-read sequencing strategy to identify polycistronic transcripts in several mushroom forming fungal species including Plicaturopsis crispa, Phanerochaete chrysosporium, Trametes versicolor, and Gloeophyllum trabeum. We found genome-wide prevalence of polycistronic transcription in these Agaricomycetes, involving up to 8% of the transcribed genes. Unlike polycistronic mRNAs in prokaryotes, these co-transcribed genes are also independently transcribed. We show that polycistronic transcription may interfere with expression of the downstream tandem gene. Further comparative genomic analysis indicates that polycistronic transcription is conserved among a wide range of mushroom forming fungi. In summary, our study revealed, for the first time, the genome prevalence of polycistronic transcription in a phylogenetic range of higher fungi. Furthermore, we systematically show that our long-read sequencing approach and combined bioinformatics pipeline is a generic powerful tool for precise characterization of complex transcriptomes that enables identification of mRNA isoforms not recovered via short-read assembly.

As decomposers, fungi are key players in recycling plant material in global carbon cycles. We hypothesized that genomes of early diverging fungi may have inherited pectinases from an ancestral species that had been able to extract nutrients from pectin-containing land plants and their algal allies (Streptophytes). We aimed to infer, based on pectinase gene expansions and on the organismal phylogeny, the geological timing of the plant-fungus association. We analyzed 40 fungal genomes, three of which, including Gonapodya prolifera, were sequenced for this study. In the organismal phylogeny from 136 housekeeping loci, Rozella diverged first from all other fungi. Gonapodya prolifera was included among the flagellated, predominantly aquatic fungal species in Chytridiomycota. Sister to Chytridiomycota were the predominantly terrestrial fungi including zygomycota I and zygomycota II, along with the ascomycetes and basidiomycetes that comprise Dikarya. The Gonapodya genome has 27 genes representing five of the seven classes of pectin-specific enzymes known from fungi. Most of these share a common ancestry with pectinases from Dikarya. Indicating functional and sequence similarity, Gonapodya, like many Dikarya, can use pectin as a carbon source for growth in pure culture. Shared pectinases of Dikarya and Gonapodya provide evidence that even ancient aquatic fungi had adapted to extract nutrients from the plants in the green lineage. This implies that 750 million years, the estimated maximum age of origin of the pectin-containing streptophytes represents a maximum age for the divergence of Chytridiomycota from the lineage including Dikarya.

The draft genome sequence of the smut fungus Tilletiaria anomala UBC 951 (Basidiomycota, Ustilaginomycotina) is presented. The sequenced genome size is 18.7 Mb, consisting of 289 scaffolds and a total of 6,810 predicted genes. This is the first genome sequence published for a fungus in the order Georgefisheriales (Exobasidiomycetes).

Campylobacter jejuni accounts for one of the leading causes of foodborne bacterial enteritis in humans. Despite being considered an obligate microaerobic microorganism, C. jejuni is regularly exposed to oxidative stress. However, its adaptive strategies to survive the atmospheric oxygen level during transmission to humans remain unclear. Recently, the clinical C. jejuni strain Bf was singled out for its unexpected ability to grow under ambient atmosphere. Here, we aimed to understand better the biological mechanisms underlying its atypical aerotolerance trait using two-dimensional protein electrophoresis, gene expression, and enzymatic activities. Forty-seven proteins were identified with a significantly different abundance between cultivation under microaerobic and aerobic conditions. The over-expressed proteins in aerobiosis belonged mainly to the oxidative stress response, enzymes of the tricarboxylic acid cycle, iron uptake, and regulation, and amino acid uptake when compared to microaerobic conditions. The higher abundance of proteins related to oxidative stress was correlated to dramatically higher transcript levels of the corresponding encoding genes in aerobic conditions compared to microaerobic conditions. In addition, a higher catalase-equivalent activity in strain Bf was observed. Despite the restricted catabolic capacities of C. jejuni, this study reveals that strain Bf is equipped to withstand oxidative stress. This ability could contribute to emergence and persistence of particular strains of C. jejuni throughout food processing or macrophage attack during human infection.

Alternaria brassicicola is a successful saprophyte and necrotrophic plant pathogen. Several A. brassicicola genes have been characterized as affecting pathogenesis of Brassica species. To study regulatory mechanisms of pathogenesis, we mined 421 genes in silico encoding putative transcription factors in a machine-annotated, draft genome sequence of A. brassicicola. In this study, targeted gene disruption mutants for 117 of the transcription factor genes were produced and screened. Three of these genes were associated with pathogenesis. Disruption mutants of one gene (AbPacC) were nonpathogenic and another gene (AbVf8) caused lesions less than half the diameter of wild-type lesions. Unexpectedly, mutants of the third gene, Amr1, caused lesions with a two-fold larger diameter than the wild type and complementation mutants. Amr1 is a homolog of Cmr1, a transcription factor that regulates melanin biosynthesis in several fungi. We created gene deletion mutants of Δamr1 and characterized their phenotypes. The Δamr1 mutants used pectin as a carbon source more efficiently than the wild type, were melanin-deficient, and more sensitive to UV light and glucanase digestion. The AMR1 protein was localized in the nuclei of hyphae and in highly melanized conidia during the late stage of plant pathogenesis. RNA-seq analysis revealed that three genes in the melanin biosynthesis pathway, along with the deleted Amr1 gene, were expressed at low levels in the mutants. In contrast, many hydrolytic enzyme-coding genes were expressed at higher levels in the mutants than in the wild type during pathogenesis. The results of this study suggested that a gene important for survival in nature negatively affected virulence, probably by a less efficient use of plant cell-wall materials. We speculate that the functions of the Amr1 gene are important to the success of A. brassicicola as a competitive saprophyte and plant parasite.

Saccharomyces cerevisiae strain M3707 was isolated from a sample of commercial distillers yeast, and its genome sequence together with the genome sequences for the four derived haploid strains M3836, M3837, M3838, and M3839 has been determined. Yeasts have potential for consolidated bioprocessing (CBP) for biofuel production, and access to these genome sequences will facilitate their development.

White-rot basidiomycete fungi are potent degraders of plant biomass, with the ability to mineralize all lignocellulose components. Recent comparative genomics studies showed that these fungi use a wide diversity of enzymes for wood degradation. Deeper functional analyses are however necessary to understand the enzymatic mechanisms leading to lignocellulose breakdown. The Polyporale fungus Pycnoporus coccineus BRFM310 grows well on both coniferous and deciduous wood. In the present study, we analyzed the early response of the fungus to softwood (pine) and hardwood (aspen) feedstocks and tested the effect of the secreted enzymes on lignocellulose deconstruction.

Black Sigatoka or black leaf streak disease, caused by the Dothideomycete fungus Pseudocercospora fijiensis (previously: Mycosphaerella fijiensis), is the most significant foliar disease of banana worldwide. Due to the lack of effective host resistance, management of this disease requires frequent fungicide applications, which greatly increase the economic and environmental costs to produce banana. Weekly applications in most banana plantations lead to rapid evolution of fungicide-resistant strains within populations causing disease-control failures throughout the world. Given its extremely high economic importance, two strains of P. fijiensis were sequenced and assembled with the aid of a new genetic linkage map. The 74-Mb genome of P. fijiensis is massively expanded by LTR retrotransposons, making it the largest genome within the Dothideomycetes. Melting-curve assays suggest that the genomes of two closely related members of the Sigatoka disease complex, P. eumusae and P. musae, also are expanded. Electrophoretic karyotyping and analyses of molecular markers in P. fijiensis field populations showed chromosome-length polymorphisms and high genetic diversity. Genetic differentiation was also detected using neutral markers, suggesting strong selection with limited gene flow at the studied geographic scale. Frequencies of fungicide resistance in fungicide-treated plantations were much higher than those in untreated wild-type P. fijiensis populations. A homologue of the Cladosporium fulvum Avr4 effector, PfAvr4, was identified in the P. fijiensis genome. Infiltration of the purified PfAVR4 protein into leaves of the resistant banana variety Calcutta 4 resulted in a hypersensitive-like response. This result suggests that Calcutta 4 could carry an unknown resistance gene recognizing PfAVR4. Besides adding to our understanding of the overall Dothideomycete genome structures, the P. fijiensis genome will aid in developing fungicide treatment schedules to combat this pathogen and in improving the efficiency of banana breeding programs.

Soils in boreal forests contain large stocks of carbon. Plants are the main source of this carbon through tissue residues and root exudates. A major part of the exudates are allocated to symbiotic ectomycorrhizal fungi. In return, the plant receives nutrients, in particular nitrogen from the mycorrhizal fungi. To capture the nitrogen, the fungi must at least partly disrupt the recalcitrant organic matter-protein complexes within which the nitrogen is embedded. This disruption process is poorly characterized. We used spectroscopic analyses and transcriptome profiling to examine the mechanism by which the ectomycorrhizal fungus Paxillus involutus degrades organic matter when acquiring nitrogen from plant litter. The fungus partially degraded polysaccharides and modified the structure of polyphenols. The observed chemical changes were consistent with a hydroxyl radical attack, involving Fenton chemistry similar to that of brown-rot fungi. The set of enzymes expressed by Pa. involutus during the degradation of the organic matter was similar to the set of enzymes involved in the oxidative degradation of wood by brown-rot fungi. However, Pa. involutus lacked transcripts encoding extracellular enzymes needed for metabolizing the released carbon. The saprotrophic activity has been reduced to a radical-based biodegradation system that can efficiently disrupt the organic matter-protein complexes and thereby mobilize the entrapped nutrients. We suggest that the released carbon then becomes available for further degradation and assimilation by commensal microbes, and that these activities have been lost in ectomycorrhizal fungi as an adaptation to symbiotic growth on host photosynthate. The interdependence of ectomycorrhizal symbionts and saprophytic microbes would provide a key link in the turnover of nutrients and carbon in forest ecosystems.

Current genomic perspectives on animal diversity neglect two prominent phyla, the molluscs and annelids, that together account for nearly one-third of known marine species and are important both ecologically and as experimental systems in classical embryology. Here we describe the draft genomes of the owl limpet (Lottia gigantea), a marine polychaete (Capitella teleta) and a freshwater leech (Helobdella robusta), and compare them with other animal genomes to investigate the origin and diversification of bilaterians from a genomic perspective. We find that the genome organization, gene structure and functional content of these species are more similar to those of some invertebrate deuterostome genomes (for example, amphioxus and sea urchin) than those of other protostomes that have been sequenced to date (flies, nematodes and flatworms). The conservation of these genomic features enables us to expand the inventory of genes present in the last common bilaterian ancestor, establish the tripartite diversification of bilaterians using multiple genomic characteristics and identify ancient conserved long- and short-range genetic linkages across metazoans. Superimposed on this broadly conserved pan-bilaterian background we find examples of lineage-specific genome evolution, including varying rates of rearrangement, intron gain and loss, expansions and contractions of gene families, and the evolution of clade-specific genes that produce the unique content of each genome.

The class Dothideomycetes is one of the largest groups of fungi with a high level of ecological diversity including many plant pathogens infecting a broad range of hosts. Here, we compare genome features of 18 members of this class, including 6 necrotrophs, 9 (hemi)biotrophs and 3 saprotrophs, to analyze genome structure, evolution, and the diverse strategies of pathogenesis. The Dothideomycetes most likely evolved from a common ancestor more than 280 million years ago. The 18 genome sequences differ dramatically in size due to variation in repetitive content, but show much less variation in number of (core) genes. Gene order appears to have been rearranged mostly within chromosomal boundaries by multiple inversions, in extant genomes frequently demarcated by adjacent simple repeats. Several Dothideomycetes contain one or more gene-poor, transposable element (TE)-rich putatively dispensable chromosomes of unknown function. The 18 Dothideomycetes offer an extensive catalogue of genes involved in cellulose degradation, proteolysis, secondary metabolism, and cysteine-rich small secreted proteins. Ancestors of the two major orders of plant pathogens in the Dothideomycetes, the Capnodiales and Pleosporales, may have had different modes of pathogenesis, with the former having fewer of these genes than the latter. Many of these genes are enriched in proximity to transposable elements, suggesting faster evolution because of the effects of repeat induced point (RIP) mutations. A syntenic block of genes, including oxidoreductases, is conserved in most Dothideomycetes and upregulated during infection in L. maculans, suggesting a possible function in response to oxidative stress.

Plant biomass conversion for green chemistry and bio-energy is a current challenge for a modern sustainable bioeconomy. The complex polyaromatic lignin polymers in raw biomass feedstocks (i.e., agriculture and forestry by-products) are major obstacles for biomass conversions. White-rot fungi are wood decayers able to degrade all polymers from lignocellulosic biomass including cellulose, hemicelluloses, and lignin. The white-rot fungus Polyporus brumalis efficiently breaks down lignin and is regarded as having a high potential for the initial treatment of plant biomass in its conversion to bio-energy. Here, we describe the extraordinary ability of P. brumalis for lignin degradation using its enzymatic arsenal to break down wheat straw, a lignocellulosic substrate that is considered as a biomass feedstock worldwide.

Lavender (Lavandula angustifolia) and lavandin (a sterile hybrid of L. angustifolia × L. latifolia) essential oils are among those most commonly used in the world for various industrial purposes, including perfumes, pharmaceuticals and cosmetics. The solid residues from aromatic plant distillation such as lavender- and lavandin-distilled straws are generally considered as wastes, and consequently either left in the fields or burnt. However, lavender- and lavandin-distilled straws are a potentially renewable plant biomass as they are cheap, non-food materials that can be used as raw feedstocks for green chemistry industry. The objective of this work was to assess different pathways of valorization of these straws as bio-based platform chemicals and fungal enzymes of interest in biorefinery.

The fungal genus Aspergillus is of critical importance to humankind. Species include those with industrial applications, important pathogens of humans, animals and crops, a source of potent carcinogenic contaminants of food, and an important genetic model. The genome sequences of eight aspergilli have already been explored to investigate aspects of fungal biology, raising questions about evolution and specialization within this genus.

Corynespora cassiicola is an Ascomycetes fungus with a broad host range and diverse life styles. Mostly known as a necrotrophic plant pathogen, it has also been associated with rare cases of human infection. In the rubber tree, this fungus causes the Corynespora leaf fall (CLF) disease, which increasingly affects natural rubber production in Asia and Africa. It has also been found as an endophyte in South American rubber plantations where no CLF outbreak has yet occurred. The C. cassiicola species is genetically highly diverse, but no clear relationship has been evidenced between phylogenetic lineage and pathogenicity. Cassiicolin, a small glycosylated secreted protein effector, is thought to be involved in the necrotrophic interaction with the rubber tree but some virulent C. cassiicola isolates do not have a cassiicolin gene. This study set out to identify other putative effectors involved in CLF. The genome of a highly virulent C. cassiicola isolate from the rubber tree (CCP) was sequenced and assembled. In silico prediction revealed 2870 putative effectors, comprising CAZymes, lipases, peptidases, secreted proteins and enzymes associated with secondary metabolism. Comparison with the genomes of 44 other fungal species, focusing on effector content, revealed a striking proximity with phylogenetically unrelated species (Colletotrichum acutatum, Colletotrichum gloesporioides, Fusarium oxysporum, nectria hematococca, and Botrosphaeria dothidea) sharing life style plasticity and broad host range. Candidate effectors involved in the compatible interaction with the rubber tree were identified by transcriptomic analysis. Differentially expressed genes included 92 putative effectors, among which cassiicolin and two other secreted singleton proteins. Finally, the genomes of 35 C. cassiicola isolates representing the genetic diversity of the species were sequenced and assembled, and putative effectors identified. At the intraspecific level, effector-based classification was found to be highly consistent with the phylogenomic trees. Identification of lineage-specific effectors is a key step toward understanding C. cassiicola virulence and host specialization mechanisms.

Microbial eukaryotes are integral components of natural microbial communities, and their inclusion is critical for many ecosystem studies, yet the majority of published metagenome analyses ignore eukaryotes. In order to include eukaryotes in environmental studies, we propose a method to recover eukaryotic genomes from complex metagenomic samples. A key step for genome recovery is separation of eukaryotic and prokaryotic fragments. We developed a k-mer-based strategy, EukRep, for eukaryotic sequence identification and applied it to environmental samples to show that it enables genome recovery, genome completeness evaluation, and prediction of metabolic potential. We used this approach to test the effect of addition of organic carbon on a geyser-associated microbial community and detected a substantial change of the community metabolism, with selection against almost all candidate phyla bacteria and archaea and for eukaryotes. Near complete genomes were reconstructed for three fungi placed within the Eurotiomycetes and an arthropod. While carbon fixation and sulfur oxidation were important functions in the geyser community prior to carbon addition, the organic carbon-impacted community showed enrichment for secreted proteases, secreted lipases, cellulose targeting CAZymes, and methanol oxidation. We demonstrate the broader utility of EukRep by reconstructing and evaluating relatively high-quality fungal, protist, and rotifer genomes from complex environmental samples. This approach opens the way for cultivation-independent analyses of whole microbial communities.

Dark septate endophytes (DSE) are a form-group of root endophytic fungi with elusive functions. Here, the genomes of two common DSE of semiarid areas, Cadophora sp. and Periconia macrospinosa were sequenced and analyzed with another 32 ascomycetes of different lifestyles. Cadophora sp. (Helotiales) and P. macrospinosa (Pleosporales) have genomes of 70.46 Mb and 54.99 Mb with 22,766 and 18,750 gene models, respectively. The majority of DSE-specific protein clusters lack functional annotation with no similarity to characterized proteins, implying that they have evolved unique genetic innovations. Both DSE possess an expanded number of carbohydrate active enzymes (CAZymes), including plant cell wall degrading enzymes (PCWDEs). Those were similar in three other DSE, and contributed a signal for the separation of root endophytes in principal component analyses of CAZymes, indicating shared genomic traits of DSE fungi. Number of secreted proteases and lipases, aquaporins, and genes linked to melanin synthesis were also relatively high in our fungi. In spite of certain similarities between our two DSE, we observed low levels of convergence in their gene family evolution. This suggests that, despite originating from the same habitat, these two fungi evolved along different evolutionary trajectories and display considerable functional differences within the endophytic lifestyle.

The rumen harbors a complex microbial mixture of archaea, bacteria, protozoa, and fungi that efficiently breakdown plant biomass and its complex dietary carbohydrates into soluble sugars that can be fermented and subsequently converted into metabolites and nutrients utilized by the host animal. While rumen bacterial populations have been well documented, only a fraction of the rumen eukarya are taxonomically and functionally characterized, despite the recognition that they contribute to the cellulolytic phenotype of the rumen microbiota. To investigate how anaerobic fungi actively engage in digestion of recalcitrant fiber that is resistant to degradation, we resolved genome-centric metaproteome and metatranscriptome datasets generated from switchgrass samples incubated for 48 h in nylon bags within the rumen of cannulated dairy cows. Across a gene catalog covering anaerobic rumen bacteria, fungi and viruses, a significant portion of the detected proteins originated from fungal populations. Intriguingly, the carbohydrate-active enzyme (CAZyme) profile suggested a domain-specific functional specialization, with bacterial populations primarily engaged in the degradation of hemicelluloses, whereas fungi were inferred to target recalcitrant cellulose structures via the detection of a number of endo- and exo-acting enzymes belonging to the glycoside hydrolase (GH) family 5, 6, 8, and 48. Notably, members of the GH48 family were amongst the highest abundant CAZymes and detected representatives from this family also included dockerin domains that are associated with fungal cellulosomes. A eukaryote-selected metatranscriptome further reinforced the contribution of uncultured fungi in the ruminal degradation of recalcitrant fibers. These findings elucidate the intricate networks of in situ recalcitrant fiber deconstruction, and importantly, suggest that the anaerobic rumen fungi contribute a specific set of CAZymes that complement the enzyme repertoire provided by the specialized plant cell wall degrading rumen bacteria.

From bacteria to humans, individual cells within isogenic populations can show significant variation in stress tolerance, but the nature of this heterogeneity is not clear. To investigate this, we used single-cell RNA sequencing to quantify transcript heterogeneity in single Saccharomyces cerevisiae cells treated with and without salt stress to explore population variation and identify cellular covariates that influence the stress-responsive transcriptome. Leveraging the extensive knowledge of yeast transcriptional regulation, we uncovered significant regulatory variation in individual yeast cells, both before and after stress. We also discovered that a subset of cells appears to decouple expression of ribosomal protein genes from the environmental stress response in a manner partly correlated with the cell cycle but unrelated to the yeast ultradian metabolic cycle. Live-cell imaging of cells expressing pairs of fluorescent regulators, including the transcription factor Msn2 with Dot6, Sfp1, or MAP kinase Hog1, revealed both coordinated and decoupled nucleocytoplasmic shuttling. Together with transcriptomic analysis, our results suggest that cells maintain a cellular filter against decoupled bursts of transcription factor activation but mount a stress response upon coordinated regulation, even in a subset of unstressed cells.

Many heritable mutualisms, in which beneficial symbionts are transmitted vertically between host generations, originate as antagonisms with parasite dispersal constrained by the host. Only after the parasite gains control over its transmission is the symbiosis expected to transition from antagonism to mutualism. Here, we explore this prediction in the mutualism between the fungus Rhizopus microsporus (Rm, Mucoromycotina) and a beta-proteobacterium Burkholderia, which controls host asexual reproduction. We show that reproductive addiction of Rm to endobacteria extends to mating, and is mediated by the symbiont gaining transcriptional control of the fungal ras2 gene, which encodes a GTPase central to fungal reproductive development. We also discover candidate G-protein-coupled receptors for the perception of trisporic acids, mating pheromones unique to Mucoromycotina. Our results demonstrate that regulating host asexual proliferation and modifying its sexual reproduction are sufficient for the symbiont's control of its own transmission, needed for antagonism-to-mutualism transition in heritable symbioses. These properties establish the Rm-Burkholderia symbiosis as a powerful system for identifying reproductive genes in Mucoromycotina.