Evolutionary analysis of hydrophobin gene family in two wood-degrading basidiomycetes, Phlebia brevispora and Heterobasidion annosum s.l

In Silico Evaluation, Phylogenetic Analysis, and Structural Modeling of the Class II Hydrophobin Family from Different Fungal Phytopathogens

The class II hydrophobin group (HFBII) is an extracellular group of proteins that contain the HFBII domain and eight conserved cysteine residues. These proteins are exclusively secreted by fungi and have multiple functions with a probable role as effectors. In the present study, a total of 45 amino acid sequences of hydrophobin class II proteins from different phytopathogenic fungi were retrieved from the NCBI database. We used the integration of well-designed bioinformatic tools to characterize and predict their physicochemical parameters, novel motifs, 3D structures, multiple sequence alignment (MSA), evolution, and functions as effector proteins through molecular docking. The results revealed new features for these protein members. The ProtParam tool detected the hydrophobicity properties of all proteins except for one hydrophilic protein (KAI3335996.1). Out of 45 proteins, six of them were detected as GPI-anchored proteins by the PredGPI server. Different 3D structure templates with high pTM scores were designed by Multifold v1, AlphaFold2, and trRosetta. Most of the studied proteins were anticipated as apoplastic effectors and matched with the ghyd5 gene of Fusarium graminearum as virulence factors. A protein-protein interaction (PPI) analysis unraveled the molecular function of this group as GTP-binding proteins, while a molecular docking analysis detected a chitin-binding effector role. From the MSA analysis, it was observed that the HFBII sequences shared conserved 2 Pro (P) and 2 Gly (G) amino acids besides the known eight conserved cysteine residues. The evolutionary analysis and phylogenetic tree provided evidence of episodic diversifying selection at the branch level using the aBSREL tool. A detailed in silico analysis of this family and the present findings will provide a better understanding of the HFBII characters and evolutionary relationships, which could be very useful in future studies.

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.

Cysteine-Rich Hydrophobin Gene Family: Genome Wide Analysis, Phylogeny and Transcript Profiling in Cordyceps militaris

Hydrophobins are a family of small secreted proteins found exclusively in fungi, and they play various roles in the life cycle. In the present study, genome wide analysis and transcript profiling of the hydrophobin family in Cordyceps militaris, a well-known edible and medicinal mushroom, were studied. The distribution of hydrophobins in ascomycetes with different lifestyles showed that pathogenic fungi had significantly more hydrophobins than saprotrophic fungi, and class II members accounted for the majority. Phylogenetic analysis of hydrophobin proteins from the species of Cordyceps s.l. indicated that there was more variability among the class II members than class I. Only a few hydrophobin-encoding genes evolved by duplication in Cordyceps s.l., which was inconsistent with the important role of gene duplication in basidiomycetes. Different transcript patterns of four hydrophobin-encoding genes during the life cycle indicated the possible different functions for each. The transcripts of Cmhyd2, 3 and 4 can respond to light and were related with the photoreceptors. CmQHYD, with four hydrophobin II domains, was first found in C. militaris, and multi-domain hydrophobins were only distributed in the species of Cordycipitaceae and Clavicipitaceae. These results could be helpful for further function research of hydrophobins and could provide valuable information for the evolution of hydrophobins.

Expression, purification, and refolding of diverse class IB hydrophobins

Hydrophobins are low molecular weight proteins secreted by fungi that are extremely surface-active and able to self-assemble into larger structures. Due to their unusual biochemical properties, hydrophobins are an attractive target for commercial applications such as drug emulsification and surface modification. When produced in E. coli, hydrophobins are often not soluble and need to be refolded. In this work we use SHuffle T7 Express E. coli coupled with glutathione redox buffers to produce and refold four distinct class IB hydrophobins that originate from Phanerochaete carnosa (PC1), Wallemia ichthyophaga (WI1), Serpula lacrymans (SL1), and Schizophyllum commune (SC16). Proper refolding and function of these purified hydrophobins was confirmed using nuclear magnetic resonance spectroscopy and thioflavin T assays. These results indicate that class IB hydrophobins can be consistently produced and purified from E. coli, aiding future structural and biochemical studies that require highly pure hydrophobins.

The Conifer Root and Stem Rot Pathogen (Heterobasidion parviporum): Effectome Analysis and Roles in Interspecific Fungal Interactions

Heterobasidion parviporum Niemelä & Korhonen is an economically important basidiomycete, causing root and stem rot disease of Norway spruce (Picea abies (L.) Karst) in Northern Europe. The H. parviporum genome encodes numerous small secreted proteins, which might be of importance for interacting with mycorrhiza symbionts, endophytes, and other saprotrophs. We hypothesized that small secreted proteins from H. parviporum (HpSSPs) are involved in interspecific fungal interaction. To identify HpSSP-coding genes potentially involved, we screened the H. parviporum effectome and compared their transcriptomic profiles during fungal development and in planta tree infection. We further conducted phylogenetic analysis, and identified a subset of hypothetical proteins with nonpredicted domain or unknown function as HpSSPs candidates for further characterization. The HpSSPs candidates were selected based on high-quality sequence, cysteine residue frequency, protein size, and in planta expression. We subsequently explored their roles during in vitro interaction in paired cultures of H. parviporum with ectomycorrhizal Cortinarius gentilis, endophytic Phialocephala sphaeroides, saprotrophs (Mycena sp., Phlebiopsis gigantea, and Phanerochaete chrysosporium), respectively. The transcriptomic profile revealed that a large proportion of effector candidates was either barely expressed or highly expressed under all growth conditions. In vitro dual-culture test showed that P. sphaeroides and C. gentilis were overgrown by H. parviporum. The barrage zone formation or no physical contact observed in paired cultures with the saprotrophs suggest they had either combative interaction or antibiosis effect with H. parviporum. Several HpSSPs individuals were up- or downregulated during the nonself interactions. The results of HpSSPs gene expression patterns provide additional insights into the diverse roles of SSPs in tree infection and interspecific fungal interactions.

Genome description of Phlebia radiata 79 with comparative genomics analysis on lignocellulose decomposition machinery of phlebioid fungi

BACKGROUND

The white rot fungus Phlebia radiata, a type species of the genus Phlebia, is an efficient decomposer of plant cell wall polysaccharides, modifier of softwood and hardwood lignin, and is able to produce ethanol from various waste lignocellulose substrates. Thus, P. radiata is a promising organism for biotechnological applications aiming at sustainable utilization of plant biomass. Here we report the genome sequence of P. radiata isolate 79 originally isolated from decayed alder wood in South Finland. To better understand the evolution of wood decay mechanisms in this fungus and the Polyporales phlebioid clade, gene content and clustering of genes encoding specific carbohydrate-active enzymes (CAZymes) in seven closely related fungal species was investigated. In addition, other genes encoding proteins reflecting the fungal lifestyle including peptidases, transporters, small secreted proteins and genes involved in secondary metabolism were identified in the genome assembly of P. radiata.

RESULTS

The PACBio sequenced nuclear genome of P. radiata was assembled to 93 contigs with 72X sequencing coverage and annotated, revealing a dense genome of 40.4 Mbp with approximately 14 082 predicted protein-coding genes. According to functional annotation, the genome harbors 209 glycoside hydrolase, 27 carbohydrate esterase, 8 polysaccharide lyase, and over 70 auxiliary redox enzyme-encoding genes. Comparisons with the genomes of other phlebioid fungi revealed shared and specific properties among the species with seemingly similar saprobic wood-decay lifestyles. Clustering of especially GH10 and AA9 enzyme-encoding genes according to genomic localization was discovered to be conserved among the phlebioid species. In P. radiata genome, a rich repertoire of genes involved in the production of secondary metabolites was recognized. In addition, 49 genes encoding predicted ABC proteins were identified in P. radiata genome together with 336 genes encoding peptidases, and 430 genes encoding small secreted proteins.

CONCLUSIONS

The genome assembly of P. radiata contains wide array of carbohydrate polymer attacking CAZyme and oxidoreductase genes in a composition identifiable for phlebioid white rot lifestyle in wood decomposition, and may thus serve as reference for further studies. Comparative genomics also contributed to enlightening fungal decay mechanisms in conversion and cycling of recalcitrant organic carbon in the forest ecosystems.

Comparative genomics and expression levels of hydrophobins from eight mycorrhizal genomes

Hydrophobins are small secreted proteins that are present as several gene copies in most fungal genomes. Their properties are now well understood: they are amphiphilic and assemble at hydrophilic/hydrophobic interfaces. However, their physiological functions remain largely unexplored, especially within mycorrhizal fungi. In this study, we identified hydrophobin genes and analysed their distribution in eight mycorrhizal genomes. We then measured their expression levels in three different biological conditions (mycorrhizal tissue vs. free-living mycelium, organic vs. mineral growth medium and aerial vs. submerged growth). Results confirmed that the size of the hydrophobin repertoire increased in the terminal orders of the fungal evolutionary tree. Reconciliation analysis predicted that in 41% of the cases, hydrophobins evolved from duplication events. Whatever the treatment and the fungal species, the pattern of expression of hydrophobins followed a reciprocal function, with one gene much more expressed than others from the same repertoire. These most-expressed hydrophobin genes were also among the most expressed of the whole genome, which suggests that they play a role as structural proteins. The fine-tuning of the expression of hydrophobin genes in each condition appeared complex because it differed considerably between species, in a way that could not be explained by simple ecological traits. Hydrophobin gene regulation in mycorrhizal tissue as compared with free-living mycelium, however, was significantly associated with a calculated high exposure of hydrophilic residues.

Characterization of a Basidiomycota hydrophobin reveals the structural basis for a high-similarity Class I subdivision

Class I hydrophobins are functional amyloids secreted by fungi. They self-assemble into organized films at interfaces producing structures that include cellular adhesion points and hydrophobic coatings. Here, we present the first structure and solution properties of a unique Class I protein sequence of Basidiomycota origin: the Schizophyllum commune hydrophobin SC16 (hyd1). While the core β-barrel structure and disulphide bridging characteristic of the hydrophobin family are conserved, its surface properties and secondary structure elements are reminiscent of both Class I and II hydrophobins. Sequence analyses of hydrophobins from 215 fungal species suggest this structure is largely applicable to a high-identity Basidiomycota Class I subdivision (IB). To validate this prediction, structural analysis of a comparatively distinct Class IB sequence from a different fungal order, namely the Phanerochaete carnosa PcaHyd1, indicates secondary structure properties similar to that of SC16. Together, these results form an experimental basis for a high-identity Class I subdivision and contribute to our understanding of functional amyloid formation.

CYPome of the conifer pathogen Heterobasidion irregulare: Inventory, phylogeny, and transcriptional analysis of the response to biocontrol

The molecular mechanisms underlying the interaction of the pathogen, Heterobasidion annosum s.l., the conifer tree and the biocontrol fungus, Phlebiopsis gigantea have not been fully elucidated. Members of the cytochrome P450 (CYP) protein family may contribute to the detoxification of components of chemical defence of conifer trees by H. annosum during infection. Additionally, they may also be involved in the interaction between H. annosum and P. gigantea. A genome-wide analysis of CYPs in Heterobasidion irregulare was carried out alongside gene expression studies. According to the Standardized CYP Nomenclature criteria, the H. irregulare genome has 121 CYP genes and 17 CYP pseudogenes classified into 11 clans, 35 families, and 64 subfamilies. Tandem CYP arrays originating from gene duplications and belonging to the same family and subfamily were found. Phylogenetic analysis showed that all the families of H. irregulare CYPs were monophyletic groups except for the family CYP5144. Microarray analysis revealed the transcriptional pattern for 130 transcripts of CYP-encoding genes during growth on culture filtrate produced by P. gigantea. The high level of P450 gene diversity identified in this study could result from extensive gene duplications presumably caused by the high metabolic demands of H. irregulare in its ecological niches.

De novo transcriptomic assembly and profiling of Rigidoporus microporus during saprotrophic growth on rubber wood

BACKGROUND

The basidiomycete Rigidoporus microporus is a fungus that causes the white rot disease of the tropical rubber tree, Hevea brasiliensis, the major source of commercial natural rubber. Besides its lifestyle as a pathogen, the fungus is known to switch to saprotrophic growth on wood with the ability to degrade both lignin and cellulose. There is almost no genomic or transcriptomic information on the saprotrophic abilities of this fungus. In this study, we present the fungal transcriptomic profiles during saprotrophic growth on rubber wood.

RESULTS

A total of 266.6 million RNA-Seq reads were generated from six libraries of the fungus growing either on rubber wood or without wood. De novo assembly produced 34, 518 unigenes with an average length of 2179 bp. Annotation of unigenes using public databases; GenBank, Swiss-Prot, Kyoto Encyclopedia of Genes and Genomes (KEGG), Cluster of Orthologous Groups (COG) and Gene Ontology (GO) produced 25, 880 annotated unigenes. Transcriptomic profiling analysis revealed that the fungus expressed over 300 genes encoding lignocellulolytic enzymes. Among these, 175 genes were up-regulated in rubber wood. These include three members of the glycoside hydrolase family 43, as well as various glycosyl transferases, carbohydrate esterases and polysaccharide lyases. A large number of oxidoreductases which includes nine manganese peroxidases were also significantly up-regulated in rubber wood. Several genes involved in fatty acid metabolism and degradation as well as natural rubber degradation were expressed in the transcriptome. Four genes (acyl-CoA synthetase, enoyl-CoA hydratase, 3-hydroxyacyl-CoA dehydrogenase and acyl-CoA acetyltransferase) potentially involved in rubber latex degradation pathway were also induced. A number of ATP binding cassette (ABC) transporters and hydrophobin genes were significantly expressed in the transcriptome during saprotrophic growth. Some genes related to energy metabolism were also induced.

CONCLUSIONS

The analysed data gives an insight into the activation of lignocellulose breakdown machinery of R. microporus. This study also revealed genes with relevance in antibiotic metabolism (e.g. cephalosporin esterase) as well as those with potential applications in fatty acid degradation. This is the first study on the transcriptomic analysis of R. microporus on rubber wood and should serve as a pioneering resource for future studies of the fungus at the genomic or transcriptomic level.

Transcriptomic responses of Phanerochaete chrysosporium to oak acetonic extracts: focus on a new glutathione transferase

The first steps of wood degradation by fungi lead to the release of toxic compounds known as extractives. To better understand how lignolytic fungi cope with the toxicity of these molecules, a transcriptomic analysis of Phanerochaete chrysosporium genes was performed in the presence of oak acetonic extracts. It reveals that in complement to the extracellular machinery of degradation, intracellular antioxidant and detoxification systems contribute to the lignolytic capabilities of fungi, presumably by preventing cellular damages and maintaining fungal health. Focusing on these systems, a glutathione transferase (P. chrysosporium GTT2.1 [PcGTT2.1]) has been selected for functional characterization. This enzyme, not characterized so far in basidiomycetes, has been classified first as a GTT2 compared to the Saccharomyces cerevisiae isoform. However, a deeper analysis shows that the GTT2.1 isoform has evolved functionally to reduce lipid peroxidation by recognizing high-molecular-weight peroxides as substrates. Moreover, the GTT2.1 gene has been lost in some non-wood-decay fungi. This example suggests that the intracellular detoxification system evolved concomitantly with the extracellular ligninolytic machinery in relation to the capacity of fungi to degrade wood.

Gene expression profiling of candidate virulence factors in the laminated root rot pathogen Phellinus sulphurascens

BACKGROUND

Phellinus sulphurascens is a fungal pathogen that causes laminar root rot in conifers, one of the most damaging root diseases in western North America. Despite its importance as a forest pathogen, this fungus is still poorly studied at the genomic level. An understanding of the molecular events involved in establishment of the disease should help to develop new methods for control of this disease.

RESULTS

We generated over 4600 expressed sequence tags from two cDNA libraries constructed using either mycelia grown on cellophane sheets and exposed to Douglas-fir roots or tissues from P. sulphurascens-infected Douglas-fir roots. A total of 890 unique genes were identified from the two libraries, and functional classification of 636 of these genes was possible using the Functional Catalogue (FunCat) annotation scheme. cDNAs were identified that encoded 79 potential virulence factors, including numerous genes implicated in virulence in a variety of phytopathogenic fungi. Many of these putative virulence factors were also among 82 genes identified as encoding putatively secreted proteins. The expression patterns of 86 selected fungal genes over 7 days of infection of Douglas-fir were examined using real-time PCR, and those significantly up-regulated included rhamnogalacturonan acetylesterase, 1,4-benzoquinone reductase, a cyclophilin, a glucoamylase, 3 hydrophobins, a lipase, a serine carboxypeptidase, a putative Ran-binding protein, and two unknown putatively secreted proteins called 1 J04 and 2 J12. Significantly down-regulated genes included a manganese-superoxide dismutase, two metalloproteases, and an unknown putatively secreted protein called Ps0058.

CONCLUSIONS

This first collection of Phellinus sulphurascens EST sequences and its annotation provide an important resource for future research aimed at understanding key virulence factors of this forest pathogen. We examined the expression patterns of numerous fungal genes with potential roles in virulence, and found a collection of functionally diverse genes that are significantly up- or down-regulated during infection of Douglas-fir seedling roots by P. sulphurascens.