Comparative genomics of the extremophile Cryomyces antarcticus and other psychrophilic Dothideomycetes.

Over a billion years of fungal evolution has enabled representatives of this kingdom to populate almost all parts of planet Earth and to adapt to some of its most uninhabitable environments including extremes of temperature, salinity, pH, water, light, or other sources of radiation. Cryomyces antarcticus is an endolithic fungus that inhabits rock outcrops in Antarctica. It survives extremes of cold, humidity and solar radiation in one of the least habitable environments on Earth. This fungus is unusual because it produces heavily melanized, meristematic growth and is thought to be haploid and asexual. Due to its growth in the most extreme environment, it has been suggested as an organism that could survive on Mars. However, the mechanisms it uses to achieve its extremophilic nature are not known. Comparative genomics can provide clues to the processes underlying biological diversity, evolution, and adaptation. This effort has been greatly facilitated by the 1000 Fungal Genomes project and the JGI MycoCosm portal where sequenced genomes have been assembled into phylogenetic and ecological groups representing different projects, lifestyles, ecologies, and evolutionary histories. Comparative genomics within and between these groups provides insights into fungal adaptations, for example to extreme environmental conditions. Here, we analyze two Cryomyces genomes in the context of additional psychrophilic fungi, as well as non-psychrophilic fungi with diverse lifestyles selected from the MycoCosm database. This analysis identifies families of genes that are expanded and contracted in Cryomyces and other psychrophiles and may explain their extremophilic lifestyle. Higher GC contents of genes and of bases in the third positions of codons may help to stabilize DNA under extreme conditions. Numerous smaller contigs in C. antarcticus suggest the presence of an alternative haplotype that could indicate the sequenced isolate is diploid or dikaryotic. These analyses provide a first step to unraveling the secrets of the extreme lifestyle of C. antarcticus.

Phyllosticta paracitricarpa is synonymous with the EU quarantine fungus P. citricarpa based on phylogenomic analyses.

Phyllosticta citricarpa is an important citrus-pathogen and a quarantine organism in the European Union. Its recently described relative, P. paracitricarpa, is very closely related and not listed as a quarantine organism. P. paracitricarpa is very difficult to distinguish from P. citricarpa, since its morphological features overlap and the barcoding gene sequences that were originally used to delimit them as distinct species have a low number of species-specific polymorphisms that have subsequently been shown to overlap between the two clades. Therefore, we performed extensive genomic analyses to determine whether the genetic variation between P. citricarpa and P. paracitricarpa strains should be considered to represent infraspecific variation within P. citricarpa, or whether it is indicative of distinct species. Using a phylogenomic analysis with 3,000 single copy ortholog genes and whole-genome comparisons, we determined that the variation between P. citricarpa and P. paracitricarpa can be considered as infraspecies variation within P. citricarpa. We also determined the level of variation in mitochondrial assemblies of several Phyllosticta species and concluded there are only minimal differences between the assemblies of P. citricarpa and P. paracitricarpa. Thus, using several orthogonal approaches, we here demonstrate that variation within the nuclear and mitochondrial genomes of other Phyllosticta species is larger than variation between genomes obtained from P. citricarpa and P. paracitricarpa strains. Thus, P. citricarpa and P. paracitricarpa should be considered as conspecific.

Comparative genomic analysis of thermophilic fungi reveals convergent evolutionary adaptations and gene losses.

Thermophily is a trait scattered across the fungal tree of life, with its highest prevalence within three fungal families (Chaetomiaceae, Thermoascaceae, and Trichocomaceae), as well as some members of the phylum Mucoromycota. We examined 37 thermophilic and thermotolerant species and 42 mesophilic species for this study and identified thermophily as the ancestral state of all three prominent families of thermophilic fungi. Thermophilic fungal genomes were found to encode various thermostable enzymes, including carbohydrate-active enzymes such as endoxylanases, which are useful for many industrial applications. At the same time, the overall gene counts, especially in gene families responsible for microbial defense such as secondary metabolism, are reduced in thermophiles compared to mesophiles. We also found a reduction in the core genome size of thermophiles in both the Chaetomiaceae family and the Eurotiomycetes class. The Gene Ontology terms lost in thermophilic fungi include primary metabolism, transporters, UV response, and O-methyltransferases. Comparative genomics analysis also revealed higher GC content in the third base of codons (GC3) and a lower effective number of codons in fungal thermophiles than in both thermotolerant and mesophilic fungi. Furthermore, using the Support Vector Machine classifier, we identified several Pfam domains capable of discriminating between genomes of thermophiles and mesophiles with 94% accuracy. Using AlphaFold2 to predict protein structures of endoxylanases (GH10), we built a similarity network based on the structures. We found that the number of disulfide bonds appears important for protein structure, and the network clusters based on protein structures correlate with the optimal activity temperature. Thus, comparative genomics offers new insights into the biology, adaptation, and evolutionary history of thermophilic fungi while providing a parts list for bioengineering applications.

Stress responses in an Arctic microalga (Pelagophyceae) following sudden salinity change revealed by gene expression analysis.

Marine microbes that have for eons been adapted to stable salinity regimes are confronted with sudden decreases in salinity in the Arctic Ocean. The episodic freshening is increasing due to climate change with melting multi-year sea-ice and glaciers, greater inflows from rivers, and increased precipitation. To investigate algal responses to lowered salinity, we analyzed the responses and acclimatation over 24 h in a non-model Arctic marine alga (pelagophyte CCMP2097) following transfer to realistic lower salinities. Using RNA-seq transcriptomics, here we show rapid differentially expressed genes related to stress oxidative responses, proteins involved in the photosystem and circadian clock, and those affecting lipids and inorganic ions. After 24 h the pelagophyte adjusted to the lower salinity seen in the overexpression of genes associated with freezing resistance, cold adaptation, and salt tolerance. Overall, a suite of ancient widespread pathways is recruited enabling the species to adjust to the stress of rapid salinity change.

Neurospora intermedia from a traditional fermented food enables waste-to-food conversion.

Fungal fermentation of food and agricultural by-products holds promise for improving food sustainability and security. However, the molecular basis of fungal waste-to-food upcycling remains poorly understood. Here we use a multi-omics approach to characterize oncom, a fermented food traditionally produced from soymilk by-products in Java, Indonesia. Metagenomic sequencing of samples from small-scale producers in Western Java indicated that the fungus Neurospora intermedia dominates oncom. Further transcriptomic, metabolomic and phylogenomic analysis revealed that oncom-derived N. intermedia utilizes pectin and cellulose degradation during fermentation and belongs to a genetically distinct subpopulation associated with human-generated by-products. Finally, we found that N. intermedia grew on diverse by-products such as fruit and vegetable pomace and plant-based milk waste, did not encode mycotoxins, and could create foods that were positively perceived by consumers outside Indonesia. These results showcase the traditional significance and future potential of fungal fermentation for creating delicious and nutritious foods from readily available by-products.

Symmetric and asymmetric DNA N6-adenine methylation regulates different biological responses in Mucorales.

DNA N6-adenine methylation (6mA) has recently gained importance as an epigenetic modification in eukaryotes. Its function in lineages with high levels, such as early-diverging fungi (EDF), is of particular interest. Here, we investigated the biological significance and evolutionary implications of 6mA in EDF, which exhibit divergent evolutionary patterns in 6mA usage. The analysis of two Mucorales species displaying extreme 6mA usage reveals that species with high 6mA levels show symmetric methylation enriched in highly expressed genes. In contrast, species with low 6mA levels show mostly asymmetric 6mA. Interestingly, transcriptomic regulation throughout development and in response to environmental cues is associated with changes in the 6mA landscape. Furthermore, we identify an EDF-specific methyltransferase, likely originated from endosymbiotic bacteria, as responsible for asymmetric methylation, while an MTA-70 methylation complex performs symmetric methylation. The distinct phenotypes observed in the corresponding mutants reinforced the critical role of both types of 6mA in EDF.

Genomic Analysis of Aspergillus Section Terrei Reveals a High Potential in Secondary Metabolite Production and Plant Biomass Degradation.

Aspergillus terreus has attracted interest due to its application in industrial biotechnology, particularly for the production of itaconic acid and bioactive secondary metabolites. As related species also seem to possess a prosperous secondary metabolism, they are of high interest for genome mining and exploitation. Here, we present draft genome sequences for six species from Aspergillus section Terrei and one species from Aspergillus section Nidulantes. Whole-genome phylogeny confirmed that section Terrei is monophyletic. Genome analyses identified between 70 and 108 key secondary metabolism genes in each of the genomes of section Terrei, the highest rate found in the genus Aspergillus so far. The respective enzymes fall into 167 distinct families with most of them corresponding to potentially unique compounds or compound families. Moreover, 53% of the families were only found in a single species, which supports the suitability of species from section Terrei for further genome mining. Intriguingly, this analysis, combined with heterologous gene expression and metabolite identification, suggested that species from section Terrei use a strategy for UV protection different to other species from the genus Aspergillus. Section Terrei contains a complete plant polysaccharide degrading potential and an even higher cellulolytic potential than other Aspergilli, possibly facilitating additional applications for these species in biotechnology.

Comparative transcriptomics provides insights into molecular mechanisms of zinc tolerance in the ectomycorrhizal fungus Suillus luteus.

Zinc (Zn) is a major soil contaminant and high Zn levels can disrupt growth, survival, and reproduction of fungi. Some fungal species evolved Zn tolerance through cell processes mitigating Zn toxicity, although the genes and detailed mechanisms underlying mycorrhizal fungal Zn tolerance remain unexplored. To fill this gap in knowledge, we investigated the gene expression of Zn tolerance in the ectomycorrhizal fungus Suillus luteus. We found that Zn tolerance in this species is mainly a constitutive trait that can also be environmentally dependent. Zinc tolerance in S. luteus is associated with differences in the expression of genes involved in metal exclusion and immobilization, as well as recognition and mitigation of metal-induced oxidative stress. Differentially expressed genes were predicted to be involved in transmembrane transport, metal chelation, oxidoreductase activity, and signal transduction. Some of these genes were previously reported as candidates for S. luteus Zn tolerance, while others are reported here for the first time. Our results contribute to understanding the mechanisms of fungal metal tolerance and pave the way for further research on the role of fungal metal tolerance in mycorrhizal associations.

Draft genome sequence of the Tremellomycetes yeast Papiliotrema laurentii 5307AH, isolated from aircraft.

Papiliotrema laurentii 5307AH was isolated from an aircraft polymer-coated surface. The genome size is 19,510,785 bp with a G + C content of 56%. The genome harbors genes encoding oxygenases, cutinases, lipases, and enzymes for styrene degradation, all of which could play a critical role in survival on xenobiotic surfaces.

Genome evolution and transcriptome plasticity is associated with adaptation to monocot and dicot plants in Colletotrichum fungi.

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.

Extreme overall mushroom genome expansion in Mycena s.s. irrespective of plant hosts or substrate specializations.

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.

Expression of dehydroshikimate dehydratase in poplar induces transcriptional and metabolic changes in the phenylpropanoid pathway.

Modification of lignin in feedstocks via genetic engineering aims to reduce biomass recalcitrance to facilitate efficient conversion processes. These improvements can be achieved by expressing exogenous enzymes that interfere with native biosynthetic pathways responsible for the production of the lignin precursors. In planta expression of a bacterial 3-dehydroshikimate dehydratase in poplar trees reduced lignin content and altered the monomer composition, which enabled higher yields of sugars after cell wall polysaccharide hydrolysis. Understanding how plants respond to such genetic modifications at the transcriptional and metabolic levels is needed to facilitate further improvement and field deployment. In this work, we acquired fundamental knowledge on lignin-modified poplar expressing 3-dehydroshikimate dehydratase using RNA-seq and metabolomics. The data clearly demonstrate that changes in gene expression and metabolite abundance can occur in a strict spatiotemporal fashion, revealing tissue-specific responses in the xylem, phloem, or periderm. In the poplar line that exhibited the strongest reduction in lignin, we found that 3% of the transcripts had altered expression levels and ~19% of the detected metabolites had differential abundance in the xylem from older stems. The changes affected predominantly the shikimate and phenylpropanoid pathways as well as secondary cell wall metabolism, and resulted in significant accumulation of hydroxybenzoates derived from protocatechuate and salicylate.

Diverse fates of ancient horizontal gene transfers in extremophilic red algae.

Horizontal genetic transfer (HGT) is a common phenomenon in eukaryotic genomes. However, the mechanisms by which HGT-derived genes persist and integrate into other pathways remain unclear. This topic is of significant interest because, over time, the stressors that initially favoured the fixation of HGT may diminish or disappear. Despite this, the foreign genes may continue to exist if they become part of a broader stress response or other pathways. The conventional model suggests that the acquisition of HGT equates to adaptation. However, this model may evolve into more complex interactions between gene products, a concept we refer to as the 'Integrated HGT Model' (IHM). To explore this concept further, we studied specialized HGT-derived genes that encode heavy metal detoxification functions. The recruitment of these genes into other pathways could provide clear examples of IHM. In our study, we exposed two anciently diverged species of polyextremophilic red algae from the Galdieria genus to arsenic and mercury stress in laboratory cultures. We then analysed the transcriptome data using differential and coexpression analysis. Our findings revealed that mercury detoxification follows a 'one gene-one function' model, resulting in an indivisible response. In contrast, the arsH gene in the arsenite response pathway demonstrated a complex pattern of duplication, divergence and potential neofunctionalization, consistent with the IHM. Our research sheds light on the fate and integration of ancient HGTs, providing a novel perspective on the ecology of extremophiles.

XA21-mediated resistance to Xanthomonas oryzae pv. oryzae is dose dependent.

The rice receptor kinase XA21 confers broad-spectrum resistance to Xanthomonas oryzae pv. oryzae (Xoo), the causal agent of rice bacterial blight disease. To investigate the relationship between the expression level of XA21 and resulting resistance, we generated independent HA-XA21 transgenic rice lines accumulating the XA21 immune receptor fused with an HA epitope tag. Whole-genome sequence analysis identified the T-DNA insertion sites in sixteen independent T0 events. Through quantification of the HA-XA21 protein and assessment of the resistance to Xoo strain PXO99 in six independent transgenic lines, we observed that XA21-mediated resistance is dose dependent. In contrast, based on the four agronomic traits quantified in these experiments, yield is unlikely to be affected by the expression level of HA-XA21. These findings extend our knowledge of XA21-mediated defense and contribute to the growing number of well-defined genomic landing pads in the rice genome that can be targeted for gene insertion without compromising yield.

Genome-scale model development and genomic sequencing of the oleaginous clade Lipomyces.

The Lipomyces clade contains oleaginous yeast species with advantageous metabolic features for biochemical and biofuel production. Limited knowledge about the metabolic networks of the species and limited tools for genetic engineering have led to a relatively small amount of research on the microbes. Here, a genome-scale metabolic model (GSM) of Lipomyces starkeyi NRRL Y-11557 was built using orthologous protein mappings to model yeast species. Phenotypic growth assays were used to validate the GSM (66% accuracy) and indicated that NRRL Y-11557 utilized diverse carbohydrates but had more limited catabolism of organic acids. The final GSM contained 2,193 reactions, 1,909 metabolites, and 996 genes and was thus named iLst996. The model contained 96 of the annotated carbohydrate-active enzymes. iLst996 predicted a flux distribution in line with oleaginous yeast measurements and was utilized to predict theoretical lipid yields. Twenty-five other yeasts in the Lipomyces clade were then genome sequenced and annotated. Sixteen of the Lipomyces species had orthologs for more than 97% of the iLst996 genes, demonstrating the usefulness of iLst996 as a broad GSM for Lipomyces metabolism. Pathways that diverged from iLst996 mainly revolved around alternate carbon metabolism, with ortholog groups excluding NRRL Y-11557 annotated to be involved in transport, glycerolipid, and starch metabolism, among others. Overall, this study provides a useful modeling tool and data for analyzing and understanding Lipomyces species metabolism and will assist further engineering efforts in Lipomyces.

Identification and recombinant expression of a cutinase from Papiliotrema laurentii that hydrolyzes natural and synthetic polyesters.

Given the multitude of extracellular enzymes at their disposal, many of which are designed to degrade nature's polymers (lignin, cutin, cellulose, etc.), fungi are adept at targeting synthetic polyesters with similar chemical composition. Microbial-influenced deterioration of xenobiotic polymeric surfaces is an area of interest for material scientists as these are important for the conservation of the underlying structural materials. Here, we describe the isolation and characterization of the Papiliotrema laurentii 5307AH (P. laurentii) cutinase, Plcut1. P. laurentii is basidiomycete yeast with the ability to disperse Impranil-DLN (Impranil), a colloidal polyester polyurethane, in agar plates. To test whether the fungal factor involved in this clearing was a secreted enzyme, we screened the ability of P. laurentii culture supernatants to disperse Impranil. Using size exclusion chromatography (SEC), we isolated fractions that contained Impranil-clearing activity. These fractions harbored a single ~22 kD band, which was excised and subjected to peptide sequencing. Homology searches using the peptide sequences identified, revealed that the protein Papla1 543643 (Plcut1) displays similarities to serine esterase and cutinase family of proteins. Biochemical assays using recombinant Plcut1 confirmed that this enzyme has the capability to hydrolyze Impranil, soluble esterase substrates, and apple cutin. Finally, we confirmed the presence of the Plcut1 in culture supernatants using a custom antibody that specifically recognizes this protein. The work shown here supports a major role for the Plcut1 in the fungal degradation of natural polyesters and xenobiotic polymer surfaces.IMPORTANCEFungi play a vital role in the execution of a broad range of biological processes that drive ecosystem function through production of a diverse arsenal of enzymes. However, the universal reactivity of these enzymes is a current problem for the built environment and the undesired degradation of polymeric materials in protective coatings. Here, we report the identification and characterization of a hydrolase from Papiliotrema laurentii 5307AH, an aircraft-derived fungal isolate found colonizing a biodeteriorated polymer-coated surface. We show that P. laurentii secretes a cutinase capable of hydrolyzing soluble esters as well as ester-based compounds forming solid surface coatings. These findings indicate that this fungus plays a significant role in biodeterioration through the production of a cutinase adept at degrading ester-based polymers, some of which form the backbone of protective surface coatings. The work shown here provides insights into the mechanisms employed by fungi to degrade xenobiotic polymers.

Draft genome sequencing and assembly of Favolaschia claudopus CIRM-BRFM 2984 isolated from oak limbs.

Favolaschia claudopus, a wood-inhabiting basidiomycete of the Mycenaceae family, is considered an invasive species that has recently spread from Oceania to Europe. The CIRM-BRFM 2984 strain of this fungus was originally isolated from a basidiome collected from the fallen limb of a decayed oak tree in Southwest France. The genome sequence of this strain shared characteristics with other Mycenaceae species, including a large genome size and enriched content of protein-coding genes. The genome sequence provided here will facilitate further investigation on the factors that contribute to the successful global dissemination of F. claudopus.

Oxygenation influences xylose fermentation and gene expression in the yeast genera Spathaspora and Scheffersomyces.

BACKGROUND: Cost-effective production of biofuels from lignocellulose requires the fermentation of D-xylose. Many yeast species within and closely related to the genera Spathaspora and Scheffersomyces (both of the order Serinales) natively assimilate and ferment xylose. Other species consume xylose inefficiently, leading to extracellular accumulation of xylitol. Xylitol excretion is thought to be due to the different cofactor requirements of the first two steps of xylose metabolism. Xylose reductase (XR) generally uses NADPH to reduce xylose to xylitol, while xylitol dehydrogenase (XDH) generally uses NAD+ to oxidize xylitol to xylulose, creating an imbalanced redox pathway. This imbalance is thought to be particularly consequential in hypoxic or anoxic environments. RESULTS: We screened the growth of xylose-fermenting yeast species in high and moderate aeration and identified both ethanol producers and xylitol producers. Selected species were further characterized for their XR and XDH cofactor preferences by enzyme assays and gene expression patterns by RNA-Seq. Our data revealed that xylose metabolism is more redox balanced in some species, but it is strongly affected by oxygen levels. Under high aeration, most species switched from ethanol production to xylitol accumulation, despite the availability of ample oxygen to accept electrons from NADH. This switch was followed by decreases in enzyme activity and the expression of genes related to xylose metabolism, suggesting that bottlenecks in xylose fermentation are not always due to cofactor preferences. Finally, we expressed XYL genes from multiple Scheffersomyces species in a strain of Saccharomyces cerevisiae. Recombinant S. cerevisiae expressing XYL1 from Scheffersomyces xylosifermentans, which encodes an XR without a cofactor preference, showed improved anaerobic growth on xylose as the primary carbon source compared to S. cerevisiae strain expressing XYL genes from Scheffersomyces stipitis. CONCLUSION: Collectively, our data do not support the hypothesis that xylitol accumulation occurs primarily due to differences in cofactor preferences between xylose reductase and xylitol dehydrogenase; instead, gene expression plays a major role in response to oxygen levels. We have also identified the yeast Sc. xylosifermentans as a potential source for genes that can be engineered into S. cerevisiae to improve xylose fermentation and biofuel production.

On the origin of bird's nest fungi: Phylogenomic analyses of fungi in the Nidulariaceae (Agaricales, Basidiomycota).

Nidulariaceae, also known as bird's nest fungi, is an understudied group of mushroom-forming fungi. The common name is derived from their nest-like morphology. Bird's nest fungi are ubiquitous wood decomposers or saprobes on dung. Recent studies showed that species in the Nidulariaceae form a monophyletic group with five sub-clades. However, phylogenetic relationships among genera and placement of Nidulariaceae are still unclear. We present phylogenomic analyses of bird's nest fungi and related Agaricales fungi to gain insight into the evolution of Nidulariaceae. A species tree with 17 newly generated genomes of bird's nest fungi and representatives from all major clades of Agaricales was constructed using 1044 single-copy genes to explore the intergeneric relationships and pinpoint the placement of Nidulariaceae within Agaricales. We corroborated the hypothesis that bird's nest fungi are sister to Squamanitaceae, which includes mushroom-shaped fungi with a stipe and pileus that are saprobes and mycoparasites. Lastly, stochastic character mapping of discrete traits on phylogenies (SIMMAP) suggests that the ancestor of bird's nest fungi likely possessed an evanescent, globose peridium without strings attaching to the spore packets (funiculi). This analysis suggests that the funiculus was gained twice and that the persistent, cupulate peridium form was gained at least four times and lost once. However, alternative coding schemes and datasets with a wider array of Agaricales produced conflicting results during ancestral state reconstruction, indicating that there is some uncertainty in the number of peridium transitions and that taxon sampling may significantly alter ancestral state reconstructions. Overall, our results suggest that several key morphological characters of Nidulariaceae have been subject to homoplasy.