Evolution of Fungal Carbohydrate-Active Enzyme Portfolios and Adaptation to Plant Cell-Wall Polymers

Genome Characteristics of the Endophytic Fungus Talaromyces sp. DC2 Isolated from Catharanthus roseus (L.) G. Don

Talaromyces sp. DC2 is an endophytic fungus that was isolated from the stem of Catharanthus roseus (L.) G. Don in Hanoi, Vietnam and is capable of producing vinca alkaloids. This study utilizes the PacBio Sequel technology to completely sequence the whole genome of Talaromyces sp. DC2The genome study revealed that DC2 contains a total of 34.58 Mb spanned by 156 contigs, with a GC content of 46.5%. The identification and prediction of functional protein-coding genes, tRNA, and rRNA were comprehensively predicted and highly annotated using various BLAST databases, including non-redundant (Nr) protein sequence, Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG), Clusters of Orthologous Groups (COG), and Carbohydrate-Active Enzymes (CAZy) databases. The genome of DC2 has a total of 149, 227, 65, 153, 53, and 6 genes responsible for cellulose, hemicellulose, lignin, pectin, chitin, starch, and inulin degradation, respectively. The Antibiotics and Secondary Metabolites Analysis Shell (AntiSMASH) analyses revealed that strain DC2 possesses 20 biosynthetic gene clusters responsible for producing secondary metabolites. The strain DC2 has also been found to harbor the DDC gene encoding aromatic L-amino acid decarboxylase enzyme. Conclusively, this study has provided a comprehensive understanding of the processes involved in secondary metabolites and the ability of the Talaromyces sp. DC2 strain to degrade plant cell walls.

Evaluating the effects of mefenoxam on taxonomic and functional dynamics of nontarget fungal communities during carrot cultivation

Ridomil Gold SL (45.3% a.i. mefenoxam) is a widely used chemical fungicide for the control of oomycetes. However, its impact on fungal communities remains unexplored. Therefore, the goal of this study was to examine the effects of mefenoxam on the temporal dynamics of fungal taxonomic and functional diversities during carrot cultivation under four treatment groups: mefenoxam application with and without Pythium inoculation, and untreated control groups with and without Pythium inoculation. Our in vitro sensitivity assay showed that the maximum recommended concentration of mefenoxam, 0.24 ppm, did not suppress the mycelial growth of P. irregulare. At 100 ppm, mycelial growth was only reduced by 11.4%, indicating that the isolate was resistant to mefenoxam. MiSeq sequencing data revealed transient taxonomic variations among treatments 2 weeks post-treatment. Mortierella dominated the fungal community in the mefenoxam-Pythium combination treatment, as confirmed through PCR using our newly designed Mortierella-specific primers. Conversely, mefenoxam-Pythium combination had adverse effects on Penicillium, Trichoderma, and Fusarium, and decrease the overall alpha diversity. However, these compositional changes gradually reverted to those observed in the control by the 12th week. The predicted ecological functions of fungal communities in all Pythium and mefenoxam treatments shifted, leading to a decrease in symbiotrophs and plant pathogen functional groups. Moreover, the community-level physiological profiling approach, utilizing 96-well Biolog FF microplates, showed discernible variations in the utilization of 95 diverse carbon sources among the treatments. Notably, arbutin, L-arabinose, Tween 80, and succinamic acid demonstrated a strong positive association with Mortierella. Our findings demonstrate that a single application of mefenoxam at its recommended rate triggers substantial taxonomic and functional shifts in the soil fungal community. Considering this impact, the conventional agricultural practice of repeated mefenoxam application is likely to exert considerable shifts on the soil ecosystem that may affect agricultural sustainability.

Structural and functional analysis of the active cow rumen's microbial community provides a catalogue of genes and microbes participating in the deconstruction of cardoon biomass

BACKGROUND

Ruminal microbial communities enriched on lignocellulosic biomass have shown considerable promise for the discovery of microorganisms and enzymes involved in digesting cell wall compounds, a key bottleneck in the development of second-generation biofuels and bioproducts, enabling a circular bioeconomy. Cardoon (Cynara cardunculus) is a promising inedible energy crop for current and future cellulosic biorefineries and the emerging bioenergy and bioproducts industries. The rumen microbiome can be considered an anaerobic "bioreactor", where the resident microbiota carry out the depolymerization and hydrolysis of plant cell wall polysaccharides (PCWPs) through the catalytic action of fibrolytic enzymes. In this context, the rumen microbiota represents a potential source of microbes and fibrolytic enzymes suitable for biofuel production from feedstocks. In this study, metatranscriptomic and 16S rRNA sequencing were used to profile the microbiome and to investigate the genetic features within the microbial community adherent to the fiber fractions of the rumen content and to the residue of cardoon biomass incubated in the rumen of cannulated cows.

RESULTS

The metatranscriptome of the cardoon and rumen fibre-adherent microbial communities were dissected in their functional and taxonomic components. From a functional point of view, transcripts involved in the methanogenesis from CO2 and H2, and from methanol were over-represented in the cardoon-adherent microbial community and were affiliated with the Methanobrevibacter and Methanosphaera of the Euryarchaeota phylum. Transcripts encoding glycoside hydrolases (GHs), carbohydrate-binding modules (CBMs), carbohydrate esterases (CEs), polysaccharide lyases (PLs), and glycoside transferases (GTs) accounted for 1.5% (6,957) of the total RNA coding transcripts and were taxonomically affiliated to major rumen fibrolytic microbes, such as Oscillospiraceae, Fibrobacteraceae, Neocallimastigaceae, Prevotellaceae, Lachnospiraceae, and Treponemataceae. The comparison of the expression profile between cardoon and rumen fiber-adherent microbial communities highlighted that specific fibrolytic enzymes were potentially responsible for the breakdown of cardoon PCWPs, which was driven by specific taxa, mainly Ruminococcus, Treponema, and Neocallimastigaceae.

CONCLUSIONS

Analysis of 16S rRNA and metatranscriptomic sequencing data revealed that the cow rumen microbiome harbors a repertoire of new enzymes capable of degrading PCWPs. Our results demonstrate the feasibility of using metatranscriptomics of enriched microbial RNA as a potential approach for accelerating the discovery of novel cellulolytic enzymes that could be harnessed for biotechnology. This research contributes a relevant perspective towards degrading cellulosic biomass and providing an economical route to the production of advanced biofuels and high-value bioproducts.

Methionine oxidation of carbohydrate-active enzymes during white-rot wood decay

White-rot fungi employ secreted carbohydrate-active enzymes (CAZymes) along with reactive oxygen species (ROS), like hydrogen peroxide (H2O2), to degrade lignocellulose in wood. H2O2 serves as a co-substrate for key oxidoreductases during the initial decay phase. While the degradation of lignocellulose by CAZymes is well documented, the impact of ROS on the oxidation of the secreted proteins remains unclear, and the identity of the oxidized proteins is unknown. Methionine (Met) can be oxidized to Met sulfoxide (MetO) or Met sulfone (MetO2) with potential deleterious, antioxidant, or regulatory effects. Other residues, like proline (Pro), can undergo carbonylation. Using the white-rot Pycnoporus cinnabarinus grown on aspen wood, we analyzed the Met content of the secreted proteins and their susceptibility to oxidation combining H218O2 with deep shotgun proteomics. Strikingly, their overall Met content was significantly lower (1.4%) compared to intracellular proteins (2.1%), a feature conserved in fungi but not in metazoans or plants. We evidenced that a catalase, widespread in white-rot fungi, protects the secreted proteins from oxidation. Our redox proteomics approach allowed the identification of 49 oxidizable Met and 40 oxidizable Pro residues within few secreted proteins, mostly CAZymes. Interestingly, many of them had several oxidized residues localized in hotspots. Some Met, including those in GH7 cellobiohydrolases, were oxidized up to 47%, with a substantial percentage of sulfone (13%). These Met are conserved in fungal homologs, suggesting important functional roles. Our findings reveal that white-rot fungi safeguard their secreted proteins by minimizing their Met content and by scavenging ROS and pinpoint redox-active residues in CAZymes.IMPORTANCEThe study of lignocellulose degradation by fungi is critical for understanding the ecological and industrial implications of wood decay. While carbohydrate-active enzymes (CAZymes) play a well-established role in lignocellulose degradation, the impact of hydrogen peroxide (H2O2) on secreted proteins remains unclear. This study aims at evaluating the effect of H2O2 on secreted proteins, focusing on the oxidation of methionine (Met). Using the model white-rot fungi Pycnoporus cinnabarinus grown on aspen wood, we showed that fungi protect their secreted proteins from oxidation by reducing their Met content and utilizing a secreted catalase to scavenge exogenous H2O2. The research identified key oxidizable Met within secreted CAZymes. Importantly, some Met, like those of GH7 cellobiohydrolases, undergone substantial oxidation levels suggesting important roles in lignocellulose degradation. These findings highlight the adaptive mechanisms employed by white-rot fungi to safeguard their secreted proteins during wood decay and emphasize the importance of these processes in lignocellulose breakdown.

Analysis of five near-complete genome assemblies of the tomato pathogen Cladosporium fulvum uncovers additional accessory chromosomes and structural variations induced by transposable elements effecting the loss of avirulence genes

BACKGROUND

Fungal plant pathogens have dynamic genomes that allow them to rapidly adapt to adverse conditions and overcome host resistance. One way by which this dynamic genome plasticity is expressed is through effector gene loss, which enables plant pathogens to overcome recognition by cognate resistance genes in the host. However, the exact nature of these loses remains elusive in many fungi. This includes the tomato pathogen Cladosporium fulvum, which is the first fungal plant pathogen from which avirulence (Avr) genes were ever cloned and in which loss of Avr genes is often reported as a means of overcoming recognition by cognate tomato Cf resistance genes. A recent near-complete reference genome assembly of C. fulvum isolate Race 5 revealed a compartmentalized genome architecture and the presence of an accessory chromosome, thereby creating a basis for studying genome plasticity in fungal plant pathogens and its impact on avirulence genes.

RESULTS

Here, we obtained near-complete genome assemblies of four additional C. fulvum isolates. The genome assemblies had similar sizes (66.96 to 67.78 Mb), number of predicted genes (14,895 to 14,981), and estimated completeness (98.8 to 98.9%). Comparative analysis that included the genome of isolate Race 5 revealed high levels of synteny and colinearity, which extended to the density and distribution of repetitive elements and of repeat-induced point (RIP) mutations across homologous chromosomes. Nonetheless, structural variations, likely mediated by transposable elements and effecting the deletion of the avirulence genes Avr4E, Avr5, and Avr9, were also identified. The isolates further shared a core set of 13 chromosomes, but two accessory chromosomes were identified as well. Accessory chromosomes were significantly smaller in size, and one carried pseudogenized copies of two effector genes. Whole-genome alignments further revealed genomic islands of near-zero nucleotide diversity interspersed with islands of high nucleotide diversity that co-localized with repeat-rich regions. These regions were likely generated by RIP, which generally asymmetrically affected the genome of C. fulvum.

CONCLUSIONS

Our results reveal new evolutionary aspects of the C. fulvum genome and provide new insights on the importance of genomic structural variations in overcoming host resistance in fungal plant pathogens.

Enzymatic machinery of wood-inhabiting fungi that degrade temperate tree species

Deadwood provides habitat for fungi and serves diverse ecological functions in forests. We already have profound knowledge of fungal assembly processes, physiological and enzymatic activities, and resulting physico-chemical changes during deadwood decay. However, in situ detection and identification methods, fungal origins, and a mechanistic understanding of the main lignocellulolytic enzymes are lacking. This study used metaproteomics to detect the main extracellular lignocellulolytic enzymes in 12 tree species in a temperate forest that have decomposed for 8 ½ years. Mainly white-rot (and few brown-rot) Basidiomycota were identified as the main wood decomposers, with Armillaria as the dominant genus; additionally, several soft-rot xylariaceous Ascomycota were identified. The key enzymes involved in lignocellulolysis included manganese peroxidase, peroxide-producing alcohol oxidases, laccase, diverse glycoside hydrolases (cellulase, glucosidase, xylanase), esterases, and lytic polysaccharide monooxygenases. The fungal community and enzyme composition differed among the 12 tree species. Ascomycota species were more prevalent in angiosperm logs than in gymnosperm logs. Regarding lignocellulolysis as a function, the extracellular enzyme toolbox acted simultaneously and was interrelated (e.g. peroxidases and peroxide-producing enzymes were strongly correlated), highly functionally redundant, and present in all logs. In summary, our in situ study provides comprehensive and detailed insight into the enzymatic machinery of wood-inhabiting fungi in temperate tree species. These findings will allow us to relate changes in environmental factors to lignocellulolysis as an ecosystem function in the future.

Draft genome sequencing of Tilletia caries inciting common bunt of wheat provides pathogenicity-related genes

Common bunt of wheat caused by Tilletia caries is an important disease worldwide. The T. caries TC1_MSG genome was sequenced using the Illumina HiSeq 2500 and Nanopore ONT platforms. The Nanopore library was prepared using the ligation sequencing kit SQK-LSK110 to generate approximately 24 GB for sequencing. The assembly size of 38.18 Mb was generated with a GC content of 56.10%. The whole genome shotgun project was deposited at DDBJ/ENA/GenBank under the accession number JALUTQ000000000. Forty-six contigs were obtained with N50 of 1,798,756 bp. In total, 10,698 genes were predicted in the assembled genome. Out of 10,698 genes, 10,255 genes were predicted significantly in the genome. The repeat sequences made up approximately 1.57% of the genome. Molecular function, cellular components, and biological processes for predicted genes were mapped into the genome. In addition, repeat elements in the genome were assessed. In all, 0.89% of retroelements were observed, followed by long terminal repeat elements (0.86%) in the genome. In simple sequence repeat (SSR) analysis, 8,582 SSRs were found in the genome assembly. The trinucleotide SSR type (3,703) was the most abundant. Few putative secretory signal peptides and pathogenicity-related genes were predicted. The genomic information of T. caries will be valuable in understanding the pathogenesis mechanism as well as developing new methods for the management of the common bunt disease of wheat.

Comparative genome analysis of the freshwater fungus Filosporella fistucella indicates potential for plant-litter degradation at cold temperatures

Freshwater fungi play an important role in the decomposition of organic matter of leaf litter in rivers and streams. They also possess the necessary mechanisms to endure lower temperatures caused by habitat and weather variations. This includes the production of cold-active enzymes and antifreeze proteins. To better understand the physiological activities of freshwater fungi in their natural environment, different methods are being applied, and genome sequencing is one in the spotlight. In our study, we sequenced the first genome of the freshwater fungus Filosporella fistucella (45.7 Mb) and compared the genome with the evolutionary close-related species Tricladium varicosporioides (48.2 Mb). The genomes were annotated using the carbohydrate-active enzyme database where we then filtered for leaf-litter degradation-related enzymes (cellulase, hemicellulase, laccase, pectinase, cutinase, amylase, xylanase, and xyloglucanase). Those enzymes were analyzed for antifreeze properties using a machine-learning approach. We discovered that F. fistucella has more enzymes to participate in the breakdown of sugar, leaf, and wood than T. varicosporioides (855 and 719, respectively). Filosporella fistucella shows a larger set of enzymes capable of resisting cold temperatures than T. varicosporioides (75 and 66, respectively). Our findings indicate that in comparison with T. varicosporioides, F. fistucella has a greater capacity for aquatic growth, adaptability to freshwater environments, and resistance to low temperatures.

Genomes of fungi and relatives reveal delayed loss of ancestral gene families and evolution of key fungal traits

Fungi are ecologically important heterotrophs that have radiated into most niches on Earth and fulfil key ecological services. Despite intense interest in their origins, major genomic trends of their evolutionary route from a unicellular opisthokont ancestor to derived multicellular fungi remain poorly known. Here we provide a highly resolved genome-wide catalogue of gene family changes across fungal evolution inferred from the genomes of 123 fungi and relatives. We show that a dominant trend in early fungal evolution has been the gradual shedding of protist genes and the punctuated emergence of innovation by two main gene duplication events. We find that the gene content of non-Dikarya fungi resembles that of unicellular opisthokonts in many respects, owing to the conservation of protist genes in their genomes. The most rapidly duplicating gene groups included extracellular proteins and transcription factors, as well as ones linked to the coordination of nutrient uptake with growth, highlighting the transition to a sessile osmotrophic feeding strategy and subsequent lifestyle evolution as important elements of early fungal history. These results suggest that the genomes of pre-fungal ancestors evolved into the typical filamentous fungal genome by a combination of gradual gene loss, turnover and several large duplication events rather than by abrupt changes. Consequently, the taxonomically defined Fungi represents a genomically non-uniform assemblage of species.

Insights into the mechanisms involved in the fungal degradation of plastics

Fungi are considered among the most efficient microbial degraders of plastics, as they produce salient enzymes and can survive on recalcitrant compounds with limited nutrients. In recent years, studies have reported numerous species of fungi that can degrade different types of plastics, yet there remain many gaps in our understanding of the processes involved in biodegradation. In addition, many unknowns need to be resolved regarding the fungal enzymes responsible for plastic fragmentation and the regulatory mechanisms which fungi use to hydrolyse, assimilate and mineralize synthetic plastics. This review aims to detail the main methods used in plastic hydrolysis by fungi, key enzymatic and molecular mechanisms, chemical agents that enhance the enzymatic breakdown of plastics, and viable industrial applications. Considering that polymers such as lignin, bioplastics, phenolics, and other petroleum-based compounds exhibit closely related characteristics in terms of hydrophobicity and structure, and are degraded by similar fungal enzymes as plastics, we have reasoned that genes that have been reported to regulate the biodegradation of these compounds or their homologs could equally be involved in the regulation of plastic degrading enzymes in fungi. Thus, this review highlights and provides insight into some of the most likely regulatory mechanisms by which fungi degrade plastics, target enzymes, genes, and transcription factors involved in the process, as well as key limitations to industrial upscaling of plastic biodegradation and biological approaches that can be employed to overcome these challenges.

Whole genome sequencing and analysis of Armillaria gallica Jzi34 symbiotic with Gastrodia elata

BACKGROUND

Armillaria species are plant pathogens, but a few Armillaria species can establish a symbiotic relationship with Gastrodia elata, a rootless and leafless orchid, that is used as a Chinese herbal medicine. Armillaria is a nutrient source for the growth of G. elata. However, there are few reports on the molecular mechanism of symbiosis between Armillaria species and G. elata. The genome sequencing and analysis of Armillaria symbiotic with G. elata would provide genomic information for further studying the molecular mechanism of symbiosis.

RESULTS

The de novo genome assembly was performed with the PacBio Sequel platform and Illumina NovaSeq PE150 for the A. gallica Jzi34 strain, which was symbiotic with G. elata. Its genome assembly contained ~ 79.9 Mbp and consisted of 60 contigs with an N50 of 2,535,910 bp. There were only 4.1% repetitive sequences in the genome assembly. Functional annotation analysis revealed a total of 16,280 protein coding genes. Compared with the other five genomes of Armillaria, the carbohydrate enzyme gene family of the genome was significantly contracted, while it had the largest set of glycosyl transferase (GT) genes. It also had an expansion of auxiliary activity enzymes AA3-2 gene subfamily and cytochrome P450 genes. The synteny analysis result of P450 genes reveals that the evolutionary relationship of P450 proteins between A. gallica Jzi34 and other four Armillaria was complex.

CONCLUSIONS

These characteristics may be beneficial for establishing a symbiotic relationship with G. elata. These results explore the characteristics of A. gallica Jzi34 from a genomic perspective and provide an important genomic resource for further detailed study of Armillaria. This will help to further study the symbiotic mechanism between A. gallica and G. elata.

Carbohydrate-active enzymes in animal feed

Considering an ever-growing global population, which hit 8 billion people in the fall of 2022, it is essential to find solutions to avoid the competition between human food and animal feed for croplands. Agricultural co-products have become important components of the circular economy with their use in animal feed. Their implementation was made possible by the addition of exogenous enzymes in the diet, especially carbohydrate-active enzymes (CAZymes). In this review, we describe the diversity and versatility of microbial CAZymes targeting non-starch polysaccharides to improve the nutritional potential of diets containing cereals and protein meals. We focused our attention on cellulases, hemicellulases, pectinases which were often found to be crucial in vivo. We also highlight the performance and health benefits brought by the exogenous addition of enzymatic cocktails containing CAZymes in the diets of monogastric animals. Taking the example of the well-studied commercial cocktail Rovabio™, we discuss the evolution, constraints and future challenges faced by feed enzymes suppliers. We hope that this review will promote the use and development of enzyme solutions for industries to sustainably feed humans in the future.

Applying molecular and genetic methods to trees and their fungal communities

Forests provide invaluable economic, ecological, and social services. At the same time, they are exposed to several threats, such as fragmentation, changing climatic conditions, or increasingly destructive pests and pathogens. Trees, the inherent species of forests, cannot be viewed as isolated organisms. Manifold (micro)organisms are associated with trees playing a pivotal role in forest ecosystems. Of these organisms, fungi may have the greatest impact on the life of trees. A multitude of molecular and genetic methods are now available to investigate tree species and their associated organisms. Due to their smaller genome sizes compared to tree species, whole genomes of different fungi are routinely compared. Such studies have only recently started in forest tree species. Here, we summarize the application of molecular and genetic methods in forest conservation genetics, tree breeding, and association genetics as well as for the investigation of fungal communities and their interrelated ecological functions. These techniques provide valuable insights into the molecular basis of adaptive traits, the impacts of forest management, and changing environmental conditions on tree species and fungal communities and can enhance tree-breeding cycles due to reduced time for field testing. It becomes clear that there are multifaceted interactions among microbial species as well as between these organisms and trees. We demonstrate the versatility of the different approaches based on case studies on trees and fungi. KEY POINTS: • Current knowledge of genetic methods applied to forest trees and associated fungi. • Genomic methods are essential in conservation, breeding, management, and research. • Important role of phytobiomes for trees and their ecosystems.

Hybrid De Novo Whole-Genome Assembly, Annotation, and Identification of Secondary Metabolite Gene Clusters in the Ex-Type Strain of Chrysosporium keratinophilum

Chrysosporium is a polyphyletic genus belonging (mostly) to different families of the order Onygenales (Eurotiomycetes, Ascomycota). Certain species, such as Chrysosporium keratinophilum, are pathogenic for animals, including humans, but are also a source of proteolytic enzymes (mainly keratinases) potentially useful in bioremediation. However, only a few studies have been published regarding bioactive compounds, of which the production is mostly unpredictable due to the absence of high-quality genomic sequences. During the development of our study, the genome of the ex-type strain of Chrysosporium keratinophilum, CBS 104.66, was sequenced and assembled using a hybrid method. The results showed a high-quality genome of 25.4 Mbp in size spread across 25 contigs, with an N50 of 2.0 Mb, 34,824 coding sequences, 8002 protein sequences, 166 tRNAs, and 24 rRNAs. The functional annotation of the predicted proteins was performed using InterProScan, and the KEGG pathway mapping using BlastKOALA. The results identified a total of 3529 protein families and 856 superfamilies, which were classified into six levels and 23 KEGG categories. Subsequently, using DIAMOND, we identified 83 pathogen-host interactions (PHI) and 421 carbohydrate-active enzymes (CAZymes). Finally, the analysis using AntiSMASH showed that this strain has a total of 27 biosynthesis gene clusters (BGCs), suggesting that it has a great potential to produce a wide variety of secondary metabolites. This genomic information provides new knowledge that allows for a deeper understanding of the biology of C. keratinophilum, and offers valuable new information for further investigations of the Chrysosporium species and the order Onygenales.

Botryosphaeriaceae gene machinery: Correlation between diversity and virulence

The Botryosphaeriaceae family comprises numerous fungal pathogens capable of causing economically meaningful diseases in a wide range of crops. Many of its members can live as endophytes and turn into aggressive pathogens following the onset of environmental stress events. Their ability to cause disease may rely on the production of a broad set of effectors, such as cell wall-degrading enzymes, secondary metabolites, and peptidases. Here, we conducted comparative analyses of 41 genomes representing six Botryosphaeriaceae genera to provide insights into the genetic features linked to pathogenicity and virulence. We show that these Botryosphaeriaceae genomes possess a large diversity of carbohydrate-active enzymes (CAZymes; 128 families) and peptidases (45 families). Botryosphaeria, Neofusicoccum, and Lasiodiplodia presented the highest number of genes encoding CAZymes involved in the degradation of the plant cell wall components. The genus Botryosphaeria also exhibited the highest abundance of secreted CAZymes and peptidases. Generally, the secondary metabolites gene cluster profile was consistent in the Botryosphaeriaceae family, except for Diplodia and Neoscytalidium. At the strain level, Neofusicoccum parvum NpBt67 stood out among all the Botryosphaeriaceae genomes, presenting a higher number of secretome constituents. In contrast, the Diplodia strains showed the lowest richness of the pathogenicity- and virulence-related genes, which may correlate with their low virulence reported in previous studies. Overall, these results contribute to a better understanding of the mechanisms underlying pathogenicity and virulence in remarkable Botryosphaeriaceae species. Our results also support that Botryosphaeriaceae species could be used as an interesting biotechnological tool for lignocellulose fractionation and bioeconomy.

Unlocking the magic in mycelium: Using synthetic biology to optimize filamentous fungi for biomanufacturing and sustainability

Filamentous fungi drive carbon and nutrient cycling across our global ecosystems, through its interactions with growing and decaying flora and their constituent microbiomes. The remarkable metabolic diversity, secretion ability, and fiber-like mycelial structure that have evolved in filamentous fungi have been increasingly exploited in commercial operations. The industrial potential of mycelial fermentation ranges from the discovery and bioproduction of enzymes and bioactive compounds, the decarbonization of food and material production, to environmental remediation and enhanced agricultural production. Despite its fundamental impact in ecology and biotechnology, molds and mushrooms have not, to-date, significantly intersected with synthetic biology in ways comparable to other industrial cell factories (e.g. Escherichia coli,Saccharomyces cerevisiae, and Komagataella phaffii). In this review, we summarize a suite of synthetic biology and computational tools for the mining, engineering and optimization of filamentous fungi as a bioproduction chassis. A combination of methods across genetic engineering, mutagenesis, experimental evolution, and computational modeling can be used to address strain development bottlenecks in established and emerging industries. These include slow mycelium growth rate, low production yields, non-optimal growth in alternative feedstocks, and difficulties in downstream purification. In the scope of biomanufacturing, we then detail previous efforts in improving key bottlenecks by targeting protein processing and secretion pathways, hyphae morphogenesis, and transcriptional control. Bringing synthetic biology practices into the hidden world of molds and mushrooms will serve to expand the limited panel of host organisms that allow for commercially-feasible and environmentally-sustainable bioproduction of enzymes, chemicals, therapeutics, foods, and materials of the future.

Mycoparasites, Gut Dwellers, and Saprotrophs: Phylogenomic Reconstructions and Comparative Analyses of Kickxellomycotina Fungi

Improved sequencing technologies have profoundly altered global views of fungal diversity and evolution. High-throughput sequencing methods are critical for studying fungi due to the cryptic, symbiotic nature of many species, particularly those that are difficult to culture. However, the low coverage genome sequencing (LCGS) approach to phylogenomic inference has not been widely applied to fungi. Here we analyzed 171 Kickxellomycotina fungi using LCGS methods to obtain hundreds of marker genes for robust phylogenomic reconstruction. Additionally, we mined our LCGS data for a set of nine rDNA and protein coding genes to enable analyses across species for which no LCGS data were obtained. The main goals of this study were to: 1) evaluate the quality and utility of LCGS data for both phylogenetic reconstruction and functional annotation, 2) test relationships among clades of Kickxellomycotina, and 3) perform comparative functional analyses between clades to gain insight into putative trophic modes. In opposition to previous studies, our nine-gene analyses support two clades of arthropod gut dwelling species and suggest a possible single evolutionary event leading to this symbiotic lifestyle. Furthermore, we resolve the mycoparasitic Dimargaritales as the earliest diverging clade in the subphylum and find four major clades of Coemansia species. Finally, functional analyses illustrate clear variation in predicted carbohydrate active enzymes and secondary metabolites (SM) based on ecology, that is biotroph versus saprotroph. Saprotrophic Kickxellales broadly lack many known pectinase families compared with saprotrophic Mucoromycota and are depauperate for SM but have similar numbers of predicted chitinases as mycoparasitic.

Genome analysis and genomic comparison of a fungal cultivar of the nonsocial weevil Euops chinensis reveals its plant decomposition and protective roles in fungus-farming mutualism

Fungus-farming mutualisms are models for studying co-evolutionary among species. Compared to well-documented fungus-farming in social insects, the molecular aspects of fungus-farming mutualisms in nonsocial insects have been poorly explored. Euops chinensis is a solitary leaf-rolling weevil feeding on Japanese knotweed (Fallopia japonica). This pest has evolved a special proto-farming bipartite mutualism with the fungus Penicillium herquei, which provide nutrition and defensive protection for the E. chinensis larvae. Here, the genome of P. herquei was sequenced, and the structure and specific gene categories in the P. herquei genome were then comprehensively compared with the other two well-studied Penicillium species (P. decumbens and P. chrysogenum). The assembled P. herquei genome had a 40.25  Mb genome size with 46.7% GC content. A diverse set of genes associating with carbohydrate-active enzymes, cellulose and hemicellulose degradation, transporter, and terpenoid biosynthesis were detected in the P. herquei genome. Comparative genomics demonstrate that the three Penicillium species show similar metabolic and enzymatic potential, however, P. herquei has more genes associated with plant biomass degradation and defense but less genes associating with virulence pathogenicity. Our results provide molecular evidence for plant substrate breakdown and protective roles of P. herquei in E. chinensis mutualistic system. Large metabolic potential shared by Penicillium species at the genus level may explain why some Penicillium species are recruited by the Euops weevils as crop fungi.

Role of carbohydrate-active enzymes in mycorrhizal symbioses

Mycorrhizal fungi form mutually beneficial interactions with a wide range of terrestrial plants. During this symbiosis, the associated fungus provides mineral nutrients, such as phosphorus and nitrogen, to its host plant in exchange of photosynthesis-derived carbohydrates. Genome sequencing of mycorrhizal fungi has shown that arbuscular mycorrhizal fungi and ectomycorrhizal fungi have a restricted set of plant-cell wall degrading enzymes (PCWDE) genes, while orchid and ericoid mycorrhizal fungi have an extended PCWDE repertoire similar to soil decomposers and wood-decay fungi. On the other hand, mycorrhizal fungi have retained a substantial set of carbohydrate active enzymes (CAZymes) acting on microbial polysaccharides. Functional analysis has shown that several of the remaining PCWDEs are involved in the fungal root colonization and establishment of the symbiotic interface. In this review, we highlight the current knowledge on the evolution and function of PCWDEs in mycorrhizal fungi.

Beneficial bacterial-Auricularia cornea interactions fostering growth enhancement identified from microbiota present in spent mushroom substrate

Complex dynamic bacterial-fungal interactions play key roles during mushroom growth, ranging from mutualism to antagonism. These interactions convey a large influence on mushroom’s mycelial and fruiting body formation during mushroom cultivation. In this study, high-throughput amplicon sequencing was conducted to investigate the structure of bacterial communities in spent mushroom substrates obtained from cultivation of two different groups of Auricularia cornea with (A) high yield and (B) low yield of fruiting body production. It was found that species richness and diversity of microbiota in group (A) samples were significantly higher than in group (B) samples. Among the identified 765 bacterial OTUs, 5 bacterial species found to exhibit high differential abundance between group (A) and group (B) were Pseudonocardia mangrovi, Luteimonas composti, Paracoccus pantotrophus, Sphingobium jiangsuense, and Microvirga massiliensis. The co-cultivation with selected bacterial strains showed that A. cornea TBRC 12900 co-cultivated with P. mangrovi TBRC-BCC 42794 promoted a high level of mycelial growth. Proteomics analysis was performed to elucidate the biological activities involved in the mutualistic association between A. cornea TBRC 12900 and P. mangrovi TBRC-BCC 42794. After co-cultivation of A. cornea TBRC 12900 and P. mangrovi TBRC-BCC 42794, 1,616 proteins were detected including 578 proteins of A. cornea origin and 1,038 proteins of P. mangrovi origin. Functional analysis and PPI network construction revealed that the high level of mycelial growth in the co-culture condition most likely resulted from concerted actions of (a) carbohydrate-active enzymes including hydrolases, glycosyltransferases, and carbohydrate esterases important for carbohydrate metabolism and cell wall generation/remodeling, (b) peptidases including cysteine-, metallo-, and serine-peptidases, (c) transporters including the ABC-type transporter superfamily, the FAT transporter family, and the VGP family, and (d) proteins with proposed roles in formation of metabolites that can act as growth-promoting molecules or those normally contain antimicrobial activity (e.g., indoles, terpenes, β-lactones, lanthipeptides, iturins, and ectoines). The findings will provide novel insights into bacterial-fungal interactions during mycelial growth and fruiting body formation. Our results can be utilized for the selection of growth-promoting bacteria to improve the cultivation process of A. cornea with a high production yield, thus conveying potentially high socio-economic impact to mushroom agriculture.

Comparative Genomics of Three Aspergillus Strains Reveals Insights into Endophytic Lifestyle and Endophyte-Induced Plant Growth Promotion

Aspergillus includes both plant pathogenic and beneficial fungi. Although endophytes beneficial to plants have high potential for plant growth promotion and improving stress tolerance, studies on endophytic lifestyles and endophyte-plant interactions are still limited. Here, three endophytes belonging to Aspergillus, AS31, AS33, and AS42, were isolated. They could successfully colonize rice roots and significantly improved rice growth. The genomes of strains AS31, AS33, and AS42 were sequenced and compared with other Aspergillus species covering both pathogens and endophytes. The genomes of AS31, AS33, and AS42 were 36.8, 34.8, and 35.3 Mb, respectively. The endophytic genomes had more genes encoding carbohydrate-active enzymes (CAZymes) and small secreted proteins (SSPs) and secondary metabolism gene clusters involved in indole metabolism than the pathogens. In addition, these endophytes were able to improve Pi (phosphorus) accumulation and transport in rice by inducing the expression of Pi transport genes in rice. Specifically, inoculation with endophytes significantly increased Pi contents in roots at the early stage, while the Pi contents in inoculated shoots were significantly increased at the late stage. Our results not only provide important insights into endophyte-plant interactions but also provide strain and genome resources, paving the way for the agricultural application of Aspergillus endophytes.

Systems biology-guided understanding of white-rot fungi for biotechnological applications: a review

Plant-derived biomass is the most abundant biogenic carbon source on Earth. Despite this, only a small clade of organisms known as white-rot fungi (WRF) can efficiently break down both the polysaccharide and lignin components of plant cell walls. This unique ability imparts a key role for WRF in global carbon cycling and highlights their potential utilization in diverse biotechnological applications. To date, research on WRF has primarily focused on their extracellular ‘digestive enzymes’ whereas knowledge of their intracellular metabolism remains underexplored. Systems biology is a powerful approach to elucidate biological processes in numerous organisms, including WRF. Thus, here we review systems biology methods applied to WRF to date, highlight observations related to their intracellular metabolism, and conduct comparative extracellular proteomic analyses to establish further correlations between WRF species, enzymes, and cultivation conditions. Lastly, we discuss biotechnological opportunities of WRF as well as challenges and future research directions.

Detecting and identifying glycoside hydrolases using cyclophellitol-derived activity-based probes

The ability to detect active enzymes in a complex mixture of folded proteins (e.g., secretome, cell lysate) generally relies on observations of catalytic ability, necessitating the development of an activity assay that is compatible with the sample and selective for the enzyme(s) of interest. Deconvolution of the contributions of different enzymes to an observed catalytic ability further necessitates an often-challenging protein separation. The advent of broadly reactive activity-based probes (ABPs) for retaining glycoside hydrolases (GHs) has enabled an alternative, often complementary, assay for active GHs. Using activity-based protein profiling (ABPP) techniques, many retaining glycoside hydrolases can be separated, detected, and identified with high sensitivity and selectivity. This chapter outlines ABPP methods for the detection and identification of retaining glycoside hydrolases from microbial sources, including protein sample preparation from bacterial lysates and fungal secretomes, enzyme labeling and detection via fluorescence, and enzyme identification using affinity-based enrichment coupled to peptide sequencing following isobaric labeling.

Phylogenomics and Comparative Genomics Highlight Specific Genetic Features in Ganoderma Species

The Ganoderma species in Polyporales are ecologically and economically relevant wood decayers used in traditional medicine, but their genomic traits are still poorly documented. In the present study, we carried out a phylogenomic and comparative genomic analyses to better understand the genetic blueprint of this fungal lineage. We investigated seven Ganoderma genomes, including three new genomes, G. australe, G. leucocontextum, and G. lingzhi. The size of the newly sequenced genomes ranged from 60.34 to 84.27 Mb and they encoded 15,007 to 20,460 genes. A total of 58 species, including 40 white-rot fungi, 11 brown-rot fungi, four ectomycorrhizal fungi, one endophyte fungus, and two pathogens in Basidiomycota, were used for phylogenomic analyses based on 143 single-copy genes. It confirmed that Ganoderma species belong to the core polyporoid clade. Comparing to the other selected species, the genomes of the Ganoderma species encoded a larger set of genes involved in terpene metabolism and coding for secreted proteins (CAZymes, lipases, proteases and SSPs). Of note, G. australe has the largest genome size with no obvious genome wide duplication, but showed transposable elements (TEs) expansion and the largest set of terpene gene clusters, suggesting a high ability to produce terpenoids for medicinal treatment. G. australe also encoded the largest set of proteins containing domains for cytochrome P450s, heterokaryon incompatibility and major facilitator families. Besides, the size of G. australe secretome is the largest, including CAZymes (AA9, GH18, A01A), proteases G01, and lipases GGGX, which may enhance the catabolism of cell wall carbohydrates, proteins, and fats during hosts colonization. The current genomic resource will be used to develop further biotechnology and medicinal applications, together with ecological studies of the Ganoderma species.

In Silico Predictions of Ecological Plasticity Mediated by Protein Family Expansions in Early-Diverging Fungi

Early-diverging fungi (EDF) are ubiquitous and versatile. Their diversity is reflected in their genome sizes and complexity. For instance, multiple protein families have been reported to expand or disappear either in particular genomes or even whole lineages. The most commonly mentioned are CAZymes (carbohydrate-active enzymes), peptidases and transporters that serve multiple biological roles connected to, e.g., metabolism and nutrients intake. In order to study the link between ecology and its genomic underpinnings in a more comprehensive manner, we carried out a systematic in silico survey of protein family expansions and losses among EDF with diverse lifestyles. We found that 86 protein families are represented differently according to EDF ecological features (assessed by median count differences). Among these there are 19 families of proteases, 43 CAZymes and 24 transporters. Some of these protein families have been recognized before as serine and metallopeptidases, cellulases and other nutrition-related enzymes. Other clearly pronounced differences refer to cell wall remodelling and glycosylation. We hypothesize that these protein families altogether define the preliminary fungal adaptasome. However, our findings need experimental validation. Many of the protein families have never been characterized in fungi and are discussed in the light of fungal ecology for the first time.

Wood degradation by Fomitiporia mediterranea M. Fischer: Physiologic, metabolomic and proteomic approaches

Fomitiporia mediterranea (Fmed) is one of the main fungal species found in grapevine wood rot, also called "amadou," one of the most typical symptoms of grapevine trunk disease Esca. This fungus is functionally classified as a white-rot, able to degrade all wood structure polymers, i.e., hemicelluloses, cellulose, and the most recalcitrant component, lignin. Specific enzymes are secreted by the fungus to degrade those components, namely carbohydrate active enzymes for hemicelluloses and cellulose, which can be highly specific for given polysaccharide, and peroxidases, which enable white-rot to degrade lignin, with specificities relating to lignin composition as well. Furthermore, besides polymers, a highly diverse set of metabolites often associated with antifungal activities is found in wood, this set differing among the various wood species. Wood decayers possess the ability to detoxify these specific extractives and this ability could reflect the adaptation of these fungi to their specific environment. The aim of this study is to better understand the molecular mechanisms used by Fmed to degrade wood structure, and in particular its potential adaptation to grapevine wood. To do so, Fmed was cultivated on sawdust from different origins: grapevine, beech, and spruce. Carbon mineralization rate, mass loss, wood structure polymers contents, targeted metabolites (extractives) and secreted proteins were measured. We used the well-known white-rot model Trametes versicolor for comparison. Whereas no significant degradation was observed with spruce, a higher mass loss was measured on Fmed grapevine culture compared to beech culture. Moreover, on both substrates, a simultaneous degradation pattern was demonstrated, and proteomic analysis identified a relative overproduction of oxidoreductases involved in lignin and extractive degradation on grapevine cultures, and only few differences in carbohydrate active enzymes. These results could explain at least partially the adaptation of Fmed to grapevine wood structural composition compared to other wood species, and suggest that other biotic and abiotic factors should be considered to fully understand the potential adaptation of Fmed to its ecological niche. Proteomics data are available via ProteomeXchange with identifier PXD036889.

Identification of Copper-Containing Oxidoreductases in the Secretomes of Three Colletotrichum Species with a Focus on Copper Radical Oxidases for the Biocatalytic Production of Fatty Aldehydes

Copper radical alcohol oxidases (CRO-AlcOx), which have been recently discovered among fungal phytopathogens, are attractive for the production of fragrant fatty aldehydes. With the initial objective to investigate the secretion of CRO-AlcOx by natural fungal strains, we undertook time course analyses of the secretomes of three Colletotrichum species (C. graminicola, C. tabacum, and C. destructivum) using proteomics. The addition of a copper-manganese-ethanol mixture in the absence of any plant-biomass mimicking compounds to Colletotrichum cultures unexpectedly induced the secretion of up to 400 proteins, 29 to 52% of which were carbohydrate-active enzymes (CAZymes), including a wide diversity of copper-containing oxidoreductases from the auxiliary activities (AA) class (AA1, AA3, AA5, AA7, AA9, AA11, AA12, AA13, and AA16). Under these specific conditions, while a CRO-glyoxal oxidase from the AA5_1 subfamily was among the most abundantly secreted proteins, the targeted AA5_2 CRO-AlcOx were secreted at lower levels, suggesting heterologous expression as a more promising strategy for CRO-AlcOx production and utilization. C. tabacum and C. destructivum CRO-AlcOx were thus expressed in Pichia pastoris, and their preference toward both aromatic and aliphatic primary alcohols was assessed. The CRO-AlcOx from C. destructivum was further investigated in applied settings, revealing a full conversion of C₆ and C₈ alcohols into their corresponding fragrant aldehydes. IMPORTANCE In the context of the industrial shift toward greener processes, the biocatalytic production of aldehydes is of utmost interest owing to their importance for their use as flavor and fragrance ingredients. Copper radical alcohol oxidases (CRO-AlcOx) have the potential to become platform enzymes for the oxidation of alcohols to aldehydes. However, the secretion of CRO-AlcOx by natural fungal strains has never been explored, while the use of crude fungal secretomes is an appealing approach for industrial applications to alleviate various costs pertaining to biocatalyst production. While investigating this primary objective, the secretomics studies revealed unexpected results showing that under the oxidative stress conditions we probed, Colletotrichum species can secrete a broad diversity of copper-containing enzymes (laccases, sugar oxidoreductases, and lytic polysaccharide monooxygenases [LPMOs]) usually assigned to "plant cell wall degradation," despite the absence of any plant-biomass mimicking compound. However, in these conditions, only small amounts of CRO-AlcOx were secreted, pointing out recombinant expression as the most promising path for their biocatalytic application.

From fungal secretomes to enzymes cocktails: The path forward to bioeconomy

Bioeconomy is seen as a way to mitigate the carbon footprint of human activities by reducing at least part of the fossil resources-based economy. In this new paradigm of sustainable development, the use of enzymes as biocatalysts will play an increasing role to provide services and goods. In industry, most of multicomponent enzyme cocktails are of fungal origin. Filamentous fungi secrete complex enzyme sets called "secretomes" that can be utilized as enzyme cocktails to valorize different types of bioresources. In this review, we highlight recent advances in the study of fungal secretomes using improved computational and experimental secretomics methods, the progress in the understanding of industrially important fungi, and the discovery of new enzymatic mechanisms and interplays to degrade renewable resources rich in polysaccharides (e.g. cellulose). We review current biotechnological applications focusing on the benefits and challenges of fungal secretomes for industrial applications with some examples of commercial cocktails of fungal origin containing carbohydrate-active enzymes (CAZymes) and we discuss future trends.

A Putative Lignin Copper Oxidase from Trichoderma reesei

The ability of Trichoderma reesei, a fungus widely used for the commercial production of hemicellulases and cellulases, to grow and modify technical soda lignin was investigated. By quantifying fungal genomic DNA, T. reesei showed growth and sporulation in solid and liquid cultures containing lignin alone. The analysis of released soluble lignin and residual insoluble lignin was indicative of enzymatic oxidative conversion of phenolic lignin side chains and the modification of lignin structure by cleaving the β-O-4 linkages. The results also showed that polymerization reactions were taking place. A proteomic analysis conducted to investigate secreted proteins at days 3, 7, and 14 of growth revealed the presence of five auxiliary activity (AA) enzymes in the secretome: AA6, AA9, two AA3 enzymes), and the only copper radical oxidase encoded in the genome of T. reesei. This enzyme was heterologously produced and characterized, and its activity on lignin-derived molecules was investigated. Phylogenetic characterization demonstrated that this enzyme belonged to the AA5_1 family, which includes characterized glyoxal oxidases. However, the enzyme displayed overlapping physicochemical and catalytic properties across the AA5 family. The enzyme was remarkably stable at high pH and oxidized both, alcohols and aldehydes with preference to the alcohol group. It was also active on lignin-derived phenolic molecules as well as simple carbohydrates. HPSEC and LC-MS analyses on the reactions of the produced protein on lignin dimers (SS ββ, SS βO4 and GG β5) uncovered the polymerizing activity of this enzyme, which was accordingly named lignin copper oxidase (TrLOx). Polymers of up 10 units were formed by hydroxy group oxidation and radical formation. The activations of lignin molecules by TrLOx along with the co-secretion of this enzyme with reductases and FAD flavoproteins oxidoreductases during growth on lignin suggest a synergistic mechanism for lignin breakdown.

Multi-omics analysis provides insights into lignocellulosic biomass degradation by Laetiporus sulphureus ATCC 52600

BACKGROUND

Wood-decay basidiomycetes are effective for the degradation of highly lignified and recalcitrant plant substrates. The degradation of lignocellulosic materials by brown-rot strains is carried out by carbohydrate-active enzymes and non-enzymatic Fenton mechanism. Differences in the lignocellulose catabolism among closely related brown rots are not completely understood. Here, a multi-omics approach provided a global understanding of the strategies employed by L. sulphureus ATCC 52600 for lignocellulose degradation.

RESULTS

The genome of Laetiporus sulphureus ATCC 52600 was sequenced and phylogenomic analysis supported monophyletic clades for the Order Polyporales and classification of this species within the family Laetiporaceae. Additionally, the plasticity of its metabolism was revealed in growth analysis on mono- and disaccharides, and polysaccharides such as cellulose, hemicelluloses, and polygalacturonic acid. The response of this fungus to the presence of lignocellulosic substrates was analyzed by transcriptomics and proteomics and evidenced the occurrence of an integrated oxidative-hydrolytic metabolism. The transcriptomic profile in response to a short cultivation period on sugarcane bagasse revealed 125 upregulated transcripts, which included CAZymes (redox enzymes and hemicellulases) as well as non-CAZy redox enzymes and genes related to the synthesis of low-molecular-weight compounds. The exoproteome produced in response to extended cultivation time on Avicel, and steam-exploded sugarcane bagasse, sugarcane straw, and Eucalyptus revealed 112 proteins. Contrasting with the mainly oxidative profile observed in the transcriptome, the secretomes showed a diverse hydrolytic repertoire including constitutive cellulases and hemicellulases, in addition to 19 upregulated CAZymes. The secretome induced for 7 days on sugarcane bagasse, representative of the late response, was applied in the saccharification of hydrothermally pretreated grass (sugarcane straw) and softwood (pine) by supplementing a commercial cocktail.

CONCLUSION

This study shows the singularity of L. sulphureus ATCC 52600 compared to other Polyporales brown rots, regarding the presence of cellobiohydrolase and peroxidase class II. The multi-omics analysis reinforces the oxidative-hydrolytic metabolism involved in lignocellulose deconstruction, providing insights into the overall mechanisms as well as specific proteins of each step.