Genome-Wide Annotation and Comparative Analysis of Cytochrome P450 Monooxygenases in Basidiomycete Biotrophic Plant Pathogens

Transcriptional effects of carbon and nitrogen starvation on Ganoderma boninense, an oil palm phytopathogen

BACKGROUND

Ganoderma boninense is a phytopathogen of oil palm, causing basal and upper stem rot diseases.

METHODS

The genome sequence was used as a reference to study gene expression during growth in a starved carbon (C) and nitrogen (N) environment with minimal sugar and sawdust as initial energy sources. This study was conducted to mimic possible limitations of the C-N nutrient sources during the growth of G. boninense in oil palm plantations.

RESULTS

Genome sequencing of an isolate collected from a palm tree in West Malaysia generated an assembly of 67.12 Mb encoding 19,851 predicted genes. Transcriptomic analysis from a time course experiment during growth in this starvation media identified differentially expressed genes (DEGs) that were found to be associated with 29 metabolic pathways. During the active growth phase, 26 DEGs were related to four pathways, including secondary metabolite biosynthesis, carbohydrate metabolism, glycan metabolism and mycotoxin biosynthesis. G. boninense genes involved in the carbohydrate metabolism pathway that contribute to the degradation of plant cell walls were up-regulated. Interestingly, several genes associated with the mycotoxin biosynthesis pathway were identified as playing a possible role in pathogen-host interaction. In addition, metabolomics analysis revealed six metabolites, maltose, xylobiose, glucooligosaccharide, glycylproline, dimethylfumaric acid and arabitol that were up-regulated on Day2 of the time course experiment.

CONCLUSIONS

This study provides information on genes expressed by G. boninense in metabolic pathways that may play a role in the initial infection of the host.

Evolutionary history of the cytochrome P450s from Colletotrichum species and prediction of their putative functional roles during host-pathogen interactions

The genomes of species belonging to the genus Colletotrichum harbor a substantial number of cytochrome P450 monooxygenases (CYPs) encoded by a broad diversity of gene families. However, the biological role of their CYP complement (CYPome) has not been elucidated. Here, we investigated the putative evolutionary scenarios that occurred during the evolution of the CYPome belonging to the Colletotrichum Graminicola species complex (s.c.) and their biological implications. The study revealed that most of the CYPome gene families belonging to the Graminicola s.c. experienced gene contractions. The reductive evolution resulted in species restricted CYPs are predominant in each CYPome of members from the Graminicola s.c., whereas only 18 families are absolutely conserved among these species. However, members of CYP families displayed a notably different phylogenetic relationship at the tertiary structure level, suggesting a putative convergent evolution scenario. Most of the CYP enzymes of the Graminicola s.c. share redundant functions in secondary metabolite biosynthesis and xenobiotic metabolism. Hence, this current work suggests that the presence of a broad CYPome in the genus Colletotrichum plays a critical role in the optimization of the colonization capability and virulence.

First genome assembly and annotation of Sanghuangporus weigelae uncovers its medicinal functions, metabolic pathways, and evolution

Sanghuangporus, also known as "Sanghuang" in China, is a well-known genus of traditional Chinese medicinal macrofungi. To make more effective use of Sanghuangporus resources, we completed the first genome assembly and annotation of a monokaryon strain of S. weigelae in the present study. A 33.96-Mb genome sequence was assembled as 13 contigs, leading to prediction of 9377 protein-coding genes. Phylogenetic and average nucleotide identity analyses indicated that the S. weigelae genome is closely related to those of other Sanghuangporus species in evolutionary tree, which clustered in one clade. Collinearity analysis revealed a high level of collinearity of S. weigelae with S. baumii, S. vaninii, and S. sanghuang. Biosynthesis pathways potentially involved in medicinal properties, including terpenoid and polysaccharide synthesis, were identified in S. weigelae, while polysaccharides were identified as the main medicinal metabolites in S. weigelae, with flavonoids more important in Sanghuangporus than other medicinal mushroom groups. Genes encoding 332 carbohydrate-active enzymes were identified in the S. weigelae genome, including major glycoside hydrolases and glycosyltransferases predicted, revealing the robust lignocellulose degradation capacity of S. weigelae. Further, 130 genes, clustered in seven classes were annotated to encode cytochromes P450 in the S. weigelae genome. Overall, our results reveal the remarkably medicinal capacity of S. weigelae and provide new insights that will inform the study of evolution and medicinal application of S. weigelae. The data are a reference resource for the formulation of scientific and rational ecological protection policies for Sanghuangporus species.

Pezizomycetes Genomes Reveal Diverse P450 Complements Characteristic of Saprotrophic and Ectomycorrhizal Lifestyles

Cytochrome P450 monooxygenases (CYPs/P450s) are heme proteins that play a role in organisms' primary and secondary metabolism. P450s play an important role in organism adaptation since lifestyle influences P450 composition in their genome. This phenomenon is well-documented in bacteria but less so in fungi. This study observed this phenomenon where diverse P450 complements were identified in saprophytic and ectomycorrhizal Pezizomycetes. Genome-wide data mining, annotation, and phylogenetic analysis of P450s in 19 Pezizomycetes revealed 668 P450s that can be grouped into 153 P450 families and 245 P450 subfamilies. Only four P450 families, namely, CYP51, CYP61, CYP5093, and CYP6001, are conserved across 19 Pezizomycetes, indicating their important role in these species. A total of 5 saprophyte Pezizomycetes have 103 P450 families, whereas 14 ectomycorrhizal Pezizomycetes have 89 P450 families. Only 39 P450 families were common, and 50 and 64 P450 families, respectively, were unique to ectomycorrhizal and saprophytic Pezizomycetes. These findings suggest that the switch from a saprophytic to an ectomycorrhizal lifestyle led to both the development of diverse P450 families as well as the loss of P450s, which led to the lowest P450 family diversity, despite the emergence of novel P450 families in ectomycorrhizal Pezizomycetes.

Genomic Comparisons of Two Armillaria Species with Different Ecological Behaviors and Their Associated Soil Microbial Communities

Armillaria species show considerable variation in ecological roles and virulence, from mycorrhizae and saprophytes to important root pathogens of trees and horticultural crops. We studied two Armillaria species that can be found in coniferous forests of northwestern USA and southwestern Canada. Armillaria altimontana not only is considered as a weak, opportunistic pathogen of coniferous trees, but it also appears to exhibit in situ biological control against A. solidipes, formerly North American A. ostoyae, which is considered a virulent pathogen of coniferous trees. Here, we describe their genome assemblies and present a functional annotation of the predicted genes and proteins for the two Armillaria species that exhibit contrasting ecological roles. In addition, the soil microbial communities were examined in association with the two Armillaria species within a 45-year-old plantation of western white pine (Pinus monticola) in northern Idaho, USA, where A. altimontana was associated with improved tree growth and survival, while A. solidipes was associated with reduced growth and survival. The results from this study reveal a high similarity between the genomes of the beneficial/non-pathogenic A. altimontana and pathogenic A. solidipes; however, many relatively small differences in gene content were identified that could contribute to differences in ecological lifestyles and interactions with woody hosts and soil microbial communities.

Contrasting Health Effects of Bacteroidetes and Firmicutes Lies in Their Genomes: Analysis of P450s, Ferredoxins, and Secondary Metabolite Clusters

Species belonging to the bacterial phyla Bacteroidetes and Firmicutes represent over 90% of the gastrointestinal microbiota. Changes in the ratio of these two bacterial groups were found to have contrasting health effects, including obesity and inflammatory diseases. Despite the availability of many bacterial genomes, comparative genomic studies on the gene pools of these two bacterial groups concerning cytochrome P450 monooxygenases (P450s), ferredoxins, and secondary metabolite biosynthetic gene clusters (smBGCs) are not reported. This study is aimed to address this research gap. The study revealed the presence of diverse sets of P450s, ferredoxins, and smBGCs in their genomes. Bacteroidetes species have the highest number of P450 families, ferredoxin cluster-types, and smBGCs compared to Firmicutes species. Only four P450 families, three ferredoxin cluster types, and five smBGCs are commonly shared between these two bacterial groups. Considering the above facts, we propose that the contrasting effects of these two bacterial groups on the host are partly due to the distinct nature of secondary metabolites produced by these organisms. Thus, the cause of the contrasting health effects of these two bacterial groups lies in their gene pools.

Pseudomonas aeruginosa cytochrome P450 CYP168A1 is a fatty acid hydroxylase that metabolizes arachidonic acid to the vasodilator 19-HETE

Pseudomonas aeruginosa is a Gram-negative opportunistic human pathogen that is highly prevalent in individuals with cystic fibrosis (CF). A major problem in treating CF patients infected with P. aeruginosa is the development of antibiotic resistance. Therefore, the identification of novel P. aeruginosa antibiotic drug targets is of the utmost urgency. The genome of P. aeruginosa contains four putative cytochrome P450 enzymes (CYPs) of unknown function that have never before been characterized. Analogous to some of the CYPs from Mycobacterium tuberculosis, these P. aeruginosa CYPs may be important for growth and colonization of CF patients’ lungs. In this study, we cloned, expressed, and characterized CYP168A1 from P. aeruginosa and identified it as a subterminal fatty acid hydroxylase. Spectral binding data and computational modeling of substrates and inhibitors suggest that CYP168A1 has a large, expansive active site and preferentially binds long chain fatty acids and large hydrophobic inhibitors. Furthermore, metabolic experiments confirm that the enzyme is capable of hydroxylating arachidonic acid, an important inflammatory signaling molecule present in abundance in the CF lung, to 19-hydroxyeicosatetraenoic acid (19-HETE; K_m = 41 μM, V_max = 220 pmol/min/nmol P450), a potent vasodilator, which may play a role in the pathogen’s ability to colonize the lung. Additionally, we found that the in vitro metabolism of arachidonic acid is subject to substrate inhibition and is also inhibited by the presence of the antifungal agent ketoconazole. This study identifies a new metabolic pathway in this important human pathogen that may be of utility in treating P. aeruginosa infections.

The First Whole Genome Sequencing of Sanghuangporus sanghuang Provides Insights into Its Medicinal Application and Evolution

Sanghuangporus is a medicinal macrofungal genus typified by S. sanghuang, the very species utilized in traditional Chinese medicines by Chinese ancient people. To facilitate the medicinal application of S. sanghuang, we, for the first time, perform its genome sequencing and analyses from a monokaryon strain. A 33.34 Mb genome sequence was assembled to 26 contigs, which lead to the prediction of 8278 protein-coding genes. From these genes, the potential biosynthesis pathway of sesquiterpenoids was, for the first time, identified from Sanghuangporus, besides that of triterpenoids. While polysaccharides are the main medicinal metabolites in S. sanghuang, flavonoids are especially abundant medicinal metabolites comparing with other medicinal macrofungal groups. From the genomic perspective, S. sanghuang has a tetrapolar heterothallic mating system, and has its special nutritional strategy and advantageous medicinal properties compared with S. baumii and S. vaninii. A phylogenomics analysis indicates that Sanghuangporus emerged 15.39 million years ago and S. sanghuang has a closer phylogenetic relationship with S. baumii than S. vaninii. However, S. sanghuang shares a higher region of synteny and more orthologous genes, including carbohydrate-active enzymes with S. vaninii than S. baumii. A comparative genomics analysis with S. baumii and S. vaninii indicates that species diversification within Sanghuangporus may be driven by the translocation and translocation plus inversion of genome sequences, while the expansion and contraction of gene families may contribute to the host specificity of Sanghuangporus species. In general, the genome sequence of S. sanghuang provides insights into its medicinal application and evolution.

In Silico Structural Modeling and Analysis of Interactions of Tremellomycetes Cytochrome P450 Monooxygenases CYP51s with Substrates and Azoles

Cytochrome P450 monooxygenase CYP51 (sterol 14α-demethylase) is a well-known target of the azole drug fluconazole for treating cryptococcosis, a life-threatening fungal infection in immune-compromised patients in poor countries. Studies indicate that mutations in CYP51 confer fluconazole resistance on cryptococcal species. Despite the importance of CYP51 in these species, few studies on the structural analysis of CYP51 and its interactions with different azole drugs have been reported. We therefore performed in silico structural analysis of 11 CYP51s from cryptococcal species and other Tremellomycetes. Interactions of 11 CYP51s with nine ligands (three substrates and six azoles) performed by Rosetta docking using 10,000 combinations for each of the CYP51-ligand complex (11 CYP51s × 9 ligands = 99 complexes) and hierarchical agglomerative clustering were used for selecting the complexes. A web application for visualization of CYP51s' interactions with ligands was developed (http://bioshell.pl/azoledocking/). The study results indicated that Tremellomycetes CYP51s have a high preference for itraconazole, corroborating the in vitro effectiveness of itraconazole compared to fluconazole. Amino acids interacting with different ligands were found to be conserved across CYP51s, indicating that the procedure employed in this study is accurate and can be automated for studying P450-ligand interactions to cater for the growing number of P450s.

Comparative genome analyses suggest a hemibiotrophic lifestyle and virulence differences for the beech bark disease fungal pathogens Neonectria faginata and Neonectria coccinea

Neonectria faginata and Neonectria coccinea are the causal agents of the insect-fungus disease complex known as beech bark disease (BBD), known to cause mortality in beech forest stands in North America and Europe. These fungal species have been the focus of extensive ecological and disease management studies, yet less progress has been made toward generating genomic resources for both micro- and macro-evolutionary studies. Here, we report a 42.1 and 42.7 mb highly contiguous genome assemblies of N. faginata and N. coccinea, respectively, obtained using Illumina technology. These species share similar gene number counts (12,941 and 12,991) and percentages of predicted genes with assigned functional categories (64 and 65%). Approximately 32% of the predicted proteomes of both species are homologous to proteins involved in pathogenicity, yet N. coccinea shows a higher number of predicted mitogen-activated protein kinase genes, virulence determinants possibly contributing to differences in disease severity between N. faginata and N. coccinea. A wide range of genes encoding for carbohydrate-active enzymes capable of degradation of complex plant polysaccharides and a small number of predicted secretory effector proteins, secondary metabolite biosynthesis clusters and cytochrome oxidase P450 genes were also found. This arsenal of enzymes and effectors correlates with, and reflects, the hemibiotrophic lifestyle of these two fungal pathogens. Phylogenomic analysis and timetree estimations indicated that the N. faginata and N. coccinea species divergence may have occurred at ∼4.1 million years ago. Differences were also observed in the annotated mitochondrial genomes as they were found to be 81.7 kb (N. faginata) and 43.2 kb (N. coccinea) in size. The mitochondrial DNA expansion observed in N. faginata is attributed to the invasion of introns into diverse intra- and intergenic locations. These first draft genomes of N. faginata and N. coccinea serve as valuable tools to increase our understanding of basic genetics, evolutionary mechanisms and molecular physiology of these two nectriaceous plant pathogenic species.

Conserved white-rot enzymatic mechanism for wood decay in the Basidiomycota genus Pycnoporus

White-rot (WR) fungi are pivotal decomposers of dead organic matter in forest ecosystems and typically use a large array of hydrolytic and oxidative enzymes to deconstruct lignocellulose. However, the extent of lignin and cellulose degradation may vary between species and wood type. Here, we combined comparative genomics, transcriptomics and secretome proteomics to identify conserved enzymatic signatures at the onset of wood-decaying activity within the Basidiomycota genus Pycnoporus. We observed a strong conservation in the genome structures and the repertoires of protein-coding genes across the four Pycnoporus species described to date, despite the species having distinct geographic distributions. We further analysed the early response of P. cinnabarinus, P. coccineus and P. sanguineus to diverse (ligno)-cellulosic substrates. We identified a conserved set of enzymes mobilized by the three species for breaking down cellulose, hemicellulose and pectin. The co-occurrence in the exo-proteomes of H₂O₂-producing enzymes with H₂O₂-consuming enzymes was a common feature of the three species, although each enzymatic partner displayed independent transcriptional regulation. Finally, cellobiose dehydrogenase-coding genes were systematically co-regulated with at least one AA9 lytic polysaccharide monooxygenase gene, indicative of enzymatic synergy in vivo. This study highlights a conserved core white-rot fungal enzymatic mechanism behind the wood-decaying process.

In Silico Analysis of P450s and Their Role in Secondary Metabolism in the Bacterial Class Gammaproteobacteria

The impact of lifestyle on shaping the genome content of an organism is a well-known phenomenon and cytochrome P450 enzymes (CYPs/P450s), heme-thiolate proteins that are ubiquitously present in organisms, are no exception. Recent studies focusing on a few bacterial species such as Streptomyces, Mycobacterium, Cyanobacteria and Firmicutes revealed that the impact of lifestyle affected the P450 repertoire in these species. However, this phenomenon needs to be understood in other bacterial species. We therefore performed genome data mining, annotation, phylogenetic analysis of P450s and their role in secondary metabolism in the bacterial class Gammaproteobacteria. Genome-wide data mining for P450s in 1261 Gammaproteobacterial species belonging to 161 genera revealed that only 169 species belonging to 41 genera have P450s. A total of 277 P450s found in 169 species grouped into 84 P450 families and 105 P450 subfamilies, where 38 new P450 families were found. Only 18% of P450s were found to be involved in secondary metabolism in Gammaproteobacterial species, as observed in Firmicutes as well. The pathogenic or commensal lifestyle of Gammaproteobacterial species influences them to such an extent that they have the lowest number of P450s compared to other bacterial species, indicating the impact of lifestyle on shaping the P450 repertoire. This study is the first report on comprehensive analysis of P450s in Gammaproteobacteria.

A highly contiguous reference genome assembly for Colletotrichum falcatum pathotype Cf08 causing red rot disease in sugarcane

Among the biotic factors, which affect the productivity and quality of sugarcane, red rot disease caused by the fungal pathogen, Colletotrichum falcatum is the most devastating that cause enormous loss to millers as well as cane growers. We present a highly contiguous genome assembly of C. falcatum pathotype Cf08 which is virulent to popular sugarcane varieties grown in more than 3 million hectares in sub-tropical India. By performing long read sequencing on PacBio RSII system, 56.06 Mb assemblies with 238 contigs having N50 of 0.51 Mb and L50 of 34 was produced. A BUSCO completeness score of 97.24% (including 4.1% fragmented) of the entire C. falcatum Cf08 nuclear genome, greatly improved contiguity compared to an existing highly fragmented draft of C. falcatum Cf671 genome (48.13 Mb) was obtained. This Cf08 assembly had 54.14% GC content and possessed < 1% repetitive elements. A total of 18,635 protein-coding genes were predicted compared with 12,270 for Cf671. Among 617 CAZymes predicted, glycoside hydrolases were the predominant (298), and among 7264 genes associated with pathogenicity/virulence, 77 genes having effector functions were identified. The assembled genome showed its similarity with the genome of C. graminicola and C. higginsianum, the causal organisms of anthracnose in maize and in members of Brassicaceae, respectively. A total of 94 large sequences (> 100 kb) of Cf08 were mapped over C. higginsianum 10 of 12 chromosomes with 106 synteny blocks. Results discussed here would provide an important tool for future studies of evolutionary and functional genomics in C. falcatum.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s13205-021-02695-x.

Gene family expansions and transcriptome signatures uncover fungal adaptations to wood decay

Because they comprise some of the most efficient wood‐decayers, Polyporales fungi impact carbon cycling in forest environment. Despite continuous discoveries on the enzymatic machinery involved in wood decomposition, the vision on their evolutionary adaptation to wood decay and genome diversity remains incomplete. We combined the genome sequence information from 50 Polyporales species, including 26 newly sequenced genomes and sought for genomic and functional adaptations to wood decay through the analysis of genome composition and transcriptome responses to different carbon sources. The genomes of Polyporales from different phylogenetic clades showed poor conservation in macrosynteny, indicative of genome rearrangements. We observed different gene family expansion/contraction histories for plant cell wall degrading enzymes in core polyporoids and phlebioids and captured expansions for genes involved in signalling and regulation in the lineages of white rotters. Furthermore, we identified conserved cupredoxins, thaumatin‐like proteins and lytic polysaccharide monooxygenases with a yet uncharacterized appended module as new candidate players in wood decomposition. Given the current need for enzymatic toolkits dedicated to the transformation of renewable carbon sources, the observed genomic diversity among Polyporales strengthens the relevance of mining Polyporales biodiversity to understand the molecular mechanisms of wood decay.

Phosphate availability and ectomycorrhizal symbiosis with Pinus sylvestris have independent effects on the Paxillus involutus transcriptome

Many plant species form symbioses with ectomycorrhizal fungi, which help them forage for limiting nutrients in the soil such as inorganic phosphate (Pi). The transcriptional responses to symbiosis and nutrient-limiting conditions in ectomycorrhizal fungal hyphae, however, are largely unknown. An artificial system was developed to study ectomycorrhizal basidiomycete Paxillus involutus growth in symbiosis with its host tree Pinus sylvestris at different Pi concentrations. RNA-seq analysis was performed on P. involutus hyphae growing under Pi-limiting conditions, either in symbiosis or alone. We show that Pi starvation and ectomycorrhizal symbiosis have an independent effect on the P. involutus transcriptome. Notably, low Pi availability induces expression of newly identified putative high-affinity Pi transporter genes, while reducing the expression of putative organic acid transporters. Additionally, low Pi availability induces a close transcriptional interplay between P and N metabolism. GTP-related signalling was found to have a positive effect in the maintenance of ectomycorrhizal symbiosis, whereas multiple putative cytochrome P450 genes were found to be downregulated, unlike arbuscular mycorrhizal fungi. We provide the first evidence of global transcriptional changes induced by low Pi availability and ectomycorrhizal symbiosis in the hyphae of P. involutus, revealing both similarities and differences with better-characterized arbuscular mycorrhizal fungi.

Impact of lifestyle on cytochrome P450 monooxygenase repertoire is clearly evident in the bacterial phylum Firmicutes

Cytochrome P450 monooxygenases (CYPs/P450s), heme thiolate proteins, are well known for their role in organisms' primary and secondary metabolism. Research on eukaryotes such as animals, plants, oomycetes and fungi has shown that P450s profiles in these organisms are affected by their lifestyle. However, the impact of lifestyle on P450 profiling in bacteria is scarcely reported. This study is such an example where the impact of lifestyle seems to profoundly affect the P450 profiles in the bacterial species belonging to the phylum Firmicutes. Genome-wide analysis of P450s in 972 Firmicutes species belonging to 158 genera revealed that only 229 species belonging to 37 genera have P450s; 38% of Bacilli species, followed by 14% of Clostridia and 2.7% of other Firmicutes species, have P450s. The pathogenic or commensal lifestyle influences P450 content to such an extent that species belonging to the genera Streptococcus, Listeria, Staphylococcus, Lactobacillus, Lactococcus and Leuconostoc do not have P450s, with the exception of a handful of Staphylococcus species that have a single P450. Only 18% of P450s are found to be involved in secondary metabolism and 89 P450s that function in the synthesis of specific secondary metabolites are predicted. This study is the first report on comprehensive analysis of P450s in Firmicutes.

Comparative Analysis, Structural Insights, and Substrate/Drug Interaction of CYP128A1 in Mycobacterium tuberculosis

Cytochrome P450 monooxygenases (CYPs/P450s) are well known for their role in organisms' primary and secondary metabolism. Among 20 P450s of the tuberculosis-causing Mycobacterium tuberculosis H37Rv, CYP128A1 is particularly important owing to its involvement in synthesizing electron transport molecules such as menaquinone-9 (MK9). This study employs different in silico approaches to understand CYP128 P450 family's distribution and structural aspects. Genome data-mining of 4250 mycobacterial species has revealed the presence of 2674 CYP128 P450s in 2646 mycobacterial species belonging to six different categories. Contrast features were observed in the CYP128 gene distribution, subfamily patterns, and characteristics of the secondary metabolite biosynthetic gene cluster (BGCs) between M. tuberculosis complex (MTBC) and other mycobacterial category species. In all MTBC species (except one) CYP128 P450s belong to subfamily A, whereas subfamily B is predominant in another four mycobacterial category species. Of CYP128 P450s, 78% was a part of BGCs with CYP124A1, or together with CYP124A1 and CYP121A1. The CYP128 family ranked fifth in the conservation ranking. Unique amino acid patterns are present at the EXXR and CXG motifs. Molecular dynamic simulation studies indicate that the CYP128A1 bind to MK9 with the highest affinity compared to the azole drugs analyzed. This study provides comprehensive comparative analysis and structural insights of CYP128A1 in M. tuberculosis.

More P450s Are Involved in Secondary Metabolite Biosynthesis in Streptomyces Compared to Bacillus, Cyanobacteria, and Mycobacterium

Unraveling the role of cytochrome P450 monooxygenases (CYPs/P450s), heme-thiolate proteins present in living and non-living entities, in secondary metabolite synthesis is gaining momentum. In this direction, in this study, we analyzed the genomes of 203 Streptomyces species for P450s and unraveled their association with secondary metabolism. Our analyses revealed the presence of 5460 P450s, grouped into 253 families and 698 subfamilies. The CYP107 family was found to be conserved and highly populated in Streptomyces and Bacillus species, indicating its key role in the synthesis of secondary metabolites. Streptomyces species had a higher number of P450s than Bacillus and cyanobacterial species. The average number of secondary metabolite biosynthetic gene clusters (BGCs) and the number of P450s located in BGCs were higher in Streptomyces species than in Bacillus, mycobacterial, and cyanobacterial species, corroborating the superior capacity of Streptomyces species for generating diverse secondary metabolites. Functional analysis via data mining confirmed that many Streptomyces P450s are involved in the biosynthesis of secondary metabolites. This study was the first of its kind to conduct a comparative analysis of P450s in such a large number (203) of Streptomyces species, revealing the P450s’ association with secondary metabolite synthesis in Streptomyces species. Future studies should include the selection of Streptomyces species with a higher number of P450s and BGCs and explore the biotechnological value of secondary metabolites they produce.

Puccinia triticina pathotypes THTT and THTS display complex transcript profiles on wheat cultivar Thatcher

BACKGROUND

Wheat leaf rust is an important disease worldwide. Understanding the pathogenic molecular mechanism of Puccinia triticina Eriks. (Pt) and the inconstant toxic region is critical for managing the disease. The present study aimed to analyze the pathogenic divergence between Pt isolates.

RESULTS

Total RNA was extracted from the wheat cultivar Thatcher infected by two Pt isolates, Tc361_1 (THTT) and Tc284_2 (THTS), at 144 h post inoculation (hpi). The mRNA was then sequenced, and a total of 2784 differentially expressed genes (DEGs) were detected. Forty-five genes were specifically expressed in THTT; these genes included transcription initiation factors and genes with transmembrane transporter activity and other genes. Twenty-six genes were specifically expressed in THTS, including genes with GTPase activity, ABC transporters and other genes. Fifty-four differentially expressed candidate effectors were screened from the two isolates. Two candidate effectors were chosen and validated on tobacco, and the results showed that they could inhibit necrosis induced by BAX. qRT-PCR of 12 significant DEGs was carried out to validate that the results are similar to those of RNA-seq at 144 hpi, to show the expression levels of these DEGs in the early stage and to elucidate the differences in expression between the two Pt pathotypes.

CONCLUSION

The results obtained in this study showed that although the two pathotypes of THTT and THTS contribute similar virulence to wheat, there are a large number of genes participate in the interaction with the susceptible wheat cultivar Thatcher, and revealed the pathogenicity of rust is very complicated.

Genome comparison and transcriptome analysis of the invasive brown root rot pathogen, Phellinus noxius, from different geographic regions reveals potential enzymes associated with degradation of different wood substrates

Phellinus noxius is a root-decay pathogen with a pan-tropical/subtropical distribution that attacks a wide range of tree hosts. For this study, genomic sequencing was conducted on P. noxius isolate P919-02W.7 from Federated States of Micronesia (Pohnpei), and its gene expression profile was analyzed using different host wood (Acer, Pinus, Prunus, and Salix) substrates. The assembled genome was 33.92 Mbp with 2954 contigs and 9389 predicted genes. Only small differences were observed in size and gene content in comparison with two other P. noxius genome assemblies (isolates OVT-YTM/97 from Hong Kong, China and FFPRI411160 from Japan, respectively). Genome analysis of P. noxius isolate P919-02W.7 revealed 488 genes encoding proteins related to carbohydrate and lignin metabolism, many of these enzymes are associated with degradation of plant cell wall components. Most of the transcripts expressed by P. noxius isolate P919-02W.7 were similar regardless of wood substrates. This study highlights the vast suite of decomposing enzymes produced by P. noxius, which suggests potential for degrading diverse wood substrates, even from temperate host trees. This information contributes to our understanding of pathogen ecology, mechanisms of wood decomposition, and pathogenic/saprophytic lifestyle.

Fungal P450 monooxygenases - the diversity in catalysis and their promising roles in biocontrol activity

The fungal P_450s catalyze vital monooxygenation reactions in primary and secondary metabolism, which may lead to the production of diverse secondary metabolites. Many of these, such as from the family of trichothecenes, involve in biocontrol activities. The diversified nature of fungal P₄₅₀ monooxygenases makes their host organisms adoptable to various ecological niches. The available genome data analysis provided an insight into the activity and mechanisms of the fungal P_450s. However, still more structural and functional studies are needed to elucidate the details of its catalytic mechanism, and the advance studies are also required to decipher further about their dynamic role in various aspects of trichothecene oxygenations. This mini review will provide updated information on different fungal P₄₅₀ monooxygenases, their genetic diversity, and their role in catalyzing various biochemical reactions leading to the production of plant growth promoting secondary metabolites.

Heterologous coexpression of the benzoate-para-hydroxylase CYP53B1 with different cytochrome P450 reductases in various yeasts

Cytochrome P450 monooxygenases (P450) are enzymes with high potential as biocatalysts for industrial applications. Their large-scale applications are, however, limited by instability and requirement for coproteins and/or expensive cofactors. These problems are largely overcome when whole cells are used as biocatalysts. We previously screened various yeast species heterologously expressing self-sufficient P450s for their potential as whole-cell biocatalysts. Most P450s are, however, not self-sufficient and consist of two or three protein component systems. Therefore, in the present study, we screened different yeast species for coexpression of P450 and P450-reductase (CPR) partners, using CYP53B1 from Rhodotorula minuta as an exemplary P450. The abilities of three different coexpressed CPR partners to support P450 activity were investigated, two from basidiomycetous origin and one from an ascomycete. The various P450-CPR combinations were cloned into strains of Saccharomyces cerevisiae, Kluyveromyces marxianus, Hansenula polymorpha, Yarrowia lipolytica and Arxula adeninivorans, using a broad-range yeast expression vector. The results obtained supported the previous finding that recombinant A. adeninivorans strains perform excellently as whole-cell biocatalysts. This study also demonstrated for the first time the P450 reductase activity of the CPRs from R. minuta and U. maydis. A very interesting observation was the variation in the supportive activity provided by the different reductase partners tested and demonstrated better P450 activity enhancement by a heterologous CPR compared to its natural partner CPR. This study highlights reductase selection as a critical variable for consideration in the pursuit of optimal P450-based catalytic systems. The usefulness of A. adeninivorans as both a host for recombinant P450s and whole-cell biocatalyst was emphasized, supporting earlier findings.

Fungal evolution: major ecological adaptations and evolutionary transitions

Fungi are a highly diverse group of heterotrophic eukaryotes characterized by the absence of phagotrophy and the presence of a chitinous cell wall. While unicellular fungi are far from rare, part of the evolutionary success of the group resides in their ability to grow indefinitely as a cylindrical multinucleated cell (hypha). Armed with these morphological traits and with an extremely high metabolical diversity, fungi have conquered numerous ecological niches and have shaped a whole world of interactions with other living organisms. Herein we survey the main evolutionary and ecological processes that have guided fungal diversity. We will first review the ecology and evolution of the zoosporic lineages and the process of terrestrialization, as one of the major evolutionary transitions in this kingdom. Several plausible scenarios have been proposed for fungal terrestralization and we here propose a new scenario, which considers icy environments as a transitory niche between water and emerged land. We then focus on exploring the main ecological relationships of Fungi with other organisms (other fungi, protozoans, animals and plants), as well as the origin of adaptations to certain specialized ecological niches within the group (lichens, black fungi and yeasts). Throughout this review we use an evolutionary and comparative-genomics perspective to understand fungal ecological diversity. Finally, we highlight the importance of genome-enabled inferences to envision plausible narratives and scenarios for important transitions.

Distribution and Diversity of Cytochrome P450 Monooxygenases in the Fungal Class Tremellomycetes

Tremellomycetes, a fungal class in the subphylum Agaricomycotina, contain well-known opportunistic and emerging human pathogens. The azole drug fluconazole, used in the treatment of diseases caused by some species of Tremellomycetes, inhibits cytochrome P450 monooxygenase CYP51, an enzyme that converts lanosterol into an essential component of the fungal cell membrane ergosterol. Studies indicate that mutations and over-expression of CYP51 in species of Tremellomycetes are one of the reasons for fluconazole resistance. Moreover, the novel drug, VT-1129, that is in the pipeline is reported to exert its effect by binding and inhibiting CYP51. Despite the importance of CYPs, the CYP repertoire in species of Tremellomycetes has not been reported to date. This study intends to address this research gap. Comprehensive genome-wide CYP analysis revealed the presence of 203 CYPs (excluding 16 pseudo-CYPs) in 23 species of Tremellomycetes that can be grouped into 38 CYP families and 72 CYP subfamilies. Twenty-three CYP families are new and three CYP families (CYP5139, CYP51 and CYP61) were conserved across 23 species of Tremellomycetes. Pathogenic cryptococcal species have 50% fewer CYP genes than non-pathogenic species. The results of this study will serve as reference for future annotation and characterization of CYPs in species of Tremellomycetes.

Similarities, variations, and evolution of cytochrome P450s in Streptomyces versus Mycobacterium

Cytochrome P450 monooxygenases (P450s) found in all domains of life are known for their catalytic versatility and stereo- and regio-specific activity. While the impact of lifestyle on P450 evolution was reported in many eukaryotes, this remains to be addressed in bacteria. In this report, Streptomyces and Mycobacterium, belonging to the phylum Actinobacteria, were studied owing to their contrasting lifestyles and impacts on human. Analyses of all P450s and those predicted to be associated with secondary metabolism have revealed that different lifestyles have affected the evolution of P450s in these bacterial genera. We have found that while species in both genera have essentially the same number of P450s in the genome, Streptomyces P450s are much more diverse than those of Mycobacterium. Moreover, despite both belonging to Actinobacteria, only 21 P450 families were common, and 123 and 56 families were found to be unique to Streptomyces and Mycobacterium, respectively. The presence of a large and diverse number of P450s in Streptomyces secondary metabolism contributes to antibiotic diversity, helping to secure the niche. Conversely, based on the currently available functional data, types of secondary metabolic pathways and associated P450s, mycobacterial P450s seem to play a role in utilization or synthesis of lipids.

Comparative Analyses of Cytochrome P450s and Those Associated with Secondary Metabolism in Bacillus Species

Cytochrome P450 monooxygenases (CYPs/P450s) are among the most catalytically-diverse enzymes, capable of performing enzymatic reactions with chemo-, regio-, and stereo-selectivity. Our understanding of P450s' role in secondary metabolite biosynthesis is becoming broader. Among bacteria, Bacillus species are known to produce secondary metabolites, and recent studies have revealed the presence of secondary metabolite biosynthetic gene clusters (BGCs) in these species. However, a comprehensive comparative analysis of P450s and P450s involved in the synthesis of secondary metabolites in Bacillus species has not been reported. This study intends to address these two research gaps. In silico analysis of P450s in 128 Bacillus species revealed the presence of 507 P450s that can be grouped into 13 P450 families and 28 subfamilies. No P450 family was found to be conserved in Bacillus species. Bacillus species were found to have lower numbers of P450s, P450 families and subfamilies, and a lower P450 diversity percentage compared to mycobacterial species. This study revealed that a large number of P450s (112 P450s) are part of different secondary metabolite BGCs, and also identified an association between a specific P450 family and secondary metabolite BGCs in Bacillus species. This study opened new vistas for further characterization of secondary metabolite BGCs, especially P450s in Bacillus species.

Rapid Divergence of Genome Architectures Following the Origin of an Ectomycorrhizal Symbiosis in the Genus Amanita

Fungi are evolutionary shape shifters and adapt quickly to new environments. Ectomycorrhizal (EM) symbioses are mutualistic associations between fungi and plants and have evolved repeatedly and independently across the fungal tree of life, suggesting lineages frequently reconfigure genome content to take advantage of open ecological niches. To date analyses of genomic mechanisms facilitating EM symbioses have involved comparisons of distantly related species, but here, we use the genomes of three EM and two asymbiotic (AS) fungi from the genus Amanita as well as an AS outgroup to study genome evolution following a single origin of symbiosis. Our aim was to identify the defining features of EM genomes, but our analyses suggest no clear differentiation of genome size, gene repertoire size, or transposable element content between EM and AS species. Phylogenetic inference of gene gains and losses suggests the transition to symbiosis was dominated by the loss of plant cell wall decomposition genes, a confirmation of previous findings. However, the same dynamic defines the AS species A. inopinata, suggesting loss is not strictly associated with origin of symbiosis. Gene expansions in the common ancestor of EM Amanita were modest, but lineage specific and large gene family expansions are found in two of the three EM extant species. Even closely related EM genomes appear to share few common features. The genetic toolkit required for symbiosis appears already encoded in the genomes of saprotrophic species, and this dynamic may explain the pervasive, recurrent evolution of ectomycorrhizal associations.

Integrative visual omics of the white-rot fungus Polyporus brumalis exposes the biotechnological potential of its oxidative enzymes for delignifying raw plant biomass

BACKGROUND

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

RESULTS

We performed integrative multi-omics analyses by combining data from the fungal genome, transcriptomes, and secretomes. We found that the fungus possessed an unexpectedly large set of genes coding for Class II peroxidases involved in lignin degradation (19 genes) and GMC oxidoreductases/dehydrogenases involved in generating the hydrogen peroxide required for lignin peroxidase activity and promoting redox cycling of the fungal enzymes involved in oxidative cleavage of lignocellulose polymers (36 genes). The examination of interrelated multi-omics patterns revealed that eleven Class II Peroxidases were secreted by the fungus during fermentation and eight of them where tightly co-regulated with redox cycling enzymatic partners.

CONCLUSION

As a peculiar feature of P. brumalis, we observed gene family extension, up-regulation and secretion of an abundant set of versatile peroxidases and manganese peroxidases, compared with other Polyporales species. The orchestrated secretion of an abundant set of these delignifying enzymes and redox cycling enzymatic partners could contribute to the delignification capabilities of the fungus. Our findings highlight the diversity of wood decay mechanisms present in Polyporales and the potentiality of further exploring this taxonomic order for enzymatic functions of biotechnological interest.

Blooming of Unusual Cytochrome P450s by Tandem Duplication in the Pathogenic Fungus Conidiobolus coronatus

While the Zygomycete fungus Conidiobolus coronatus primarily infects insects, it can be pathogenic to mammals as well, including humans. High variability in the treatment of this fungal infection with currently available drugs, including azole drugs is a very common phenomenon. Azoles bind to the cytochrome P450 monooxygenases (P450s/CYP) including CYP51, a sterol 14-α-demethylase, inhibiting the synthesis of cell membrane ergosterol and thus leading to the elimination of infecting fungi. Despite P450’s role as a drug target, to date, no information on C. coronatus P450s has been reported. Genome-wide data mining has revealed the presence of 142 P450s grouped into 12 families and 21 subfamilies in C. coronatus. Except for CYP51, the remaining 11 P450 families are new (CYP5854-CYP5864). Despite having a large number of P450s among entomopathogenic fungi, C. coronatus has the lowest number of P450 families, which suggests blooming P450s. Further analysis has revealed that 79% of the same family P450s is tandemly positioned, suggesting that P450 tandem duplication led to the blooming of P450s. The results of this study; i.e., unravelling the C. coronatus P450 content, will certainly help in designing experiments to understand P450s’ role in C. coronatus physiology, including a highly variable response to azole drugs with respect to P450s.

Fungal Genomes and Insights into the Evolution of the Kingdom

The kingdom Fungi comprises species that inhabit nearly all ecosystems. Fungi exist as both free-living and symbiotic unicellular and multicellular organisms with diverse morphologies. The genomes of fungi encode genes that enable them to thrive in diverse environments, invade plant and animal cells, and participate in nutrient cycling in terrestrial and aquatic ecosystems. The continuously expanding databases of fungal genome sequences have been generated by individual and large-scale efforts such as Génolevures, Broad Institute's Fungal Genome Initiative, and the 1000 Fungal Genomes Project (http://1000.fungalgenomes.org). These efforts have produced a catalog of fungal genes and genomic organization. The genomic datasets can be utilized to better understand how fungi have adapted to their lifestyles and ecological niches. Large datasets of fungal genomic and transcriptomic data have enabled the use of novel methodologies and improved the study of fungal evolution from a molecular sequence perspective. Combined with microscopes, petri dishes, and woodland forays, genome sequencing supports bioinformatics and comparative genomics approaches as important tools in the study of the biology and evolution of fungi.

De Novo Sequencing of a Sparassis latifolia Genome and Its Associated Comparative Analyses

Known to be rich in β-glucan, Sparassis latifolia (S. latifolia) is a valuable edible fungus cultivated in East Asia. A few studies have suggested that S. latifolia is effective on antidiabetic, antihypertension, antitumor, and antiallergen medications. However, it is still unclear genetically why the fungus has these medical effects, which has become a key bottleneck for its further applications. To provide a better understanding of this fungus, we sequenced its whole genome, which has a total size of 48.13 megabases (Mb) and contains 12,471 predicted gene models. We then performed comparative and phylogenetic analyses, which indicate that S. latifolia is closely related to a few species in the antrodia clade including Fomitopsis pinicola, Wolfiporia cocos, Postia placenta, and Antrodia sinuosa. Finally, we annotated the predicted genes. Interestingly, the S. latifolia genome encodes most enzymes involved in carbohydrate and glycoconjugate metabolism and is also enriched in genes encoding enzymes critical to secondary metabolite biosynthesis and involved in indole, terpene, and type I polyketide pathways. As a conclusion, the genome content of S. latifolia sheds light on its genetic basis of the reported medicinal properties and could also be used as a reference genome for comparative studies on fungi.

In silico analysis of cytochrome P450 monooxygenases in chronic granulomatous infectious fungus Sporothrix schenckii: Special focus on CYP51

Sporotrichosis is an emerging chronic, granulomatous, subcutaneous, mycotic infection caused by Sporothrix species. Sporotrichosis is treated with the azole drug itraconazole as ketoconazole is ineffective. It is a well-known fact that azole drugs act by inhibiting cytochrome P450 monooxygenases (P450s), heme-thiolate proteins. To date, nothing is known about P450s in Sporothrix schenckii and the molecular basis of its resistance to ketoconazole. Here we present genome-wide identification, annotation, phylogenetic analysis and comprehensive P450 family-level comparative analysis of S. schenckii P450s with pathogenic fungi P450s, along with a rationale for ketoconazole resistance by S. schenckii based on in silico structural analysis of CYP51. Genome data-mining of S. schenckii revealed 40 P450s in its genome that can be grouped into 32 P450 families and 39 P450 subfamilies. Comprehensive comparative analysis of P450s revealed that S. schenckii shares 11 P450 families with plant pathogenic fungi and has three unique P450 families: CYP5077, CYP5386 and CYP5696 (novel family). Among P450s, CYP51, the main target of azole drugs was also found in S. schenckii. 3D modeling of S. schenckii CYP51 revealed the presence of characteristic P450 motifs with exceptionally large reductase interaction site 2. In silico analysis revealed number of mutations that can be associated with ketoconazole resistance, especially at the channel entrance to the active site. One of possible reason for better stabilization of itraconazole, compared to ketoconazole, is that the more extended molecule of itraconazole may form a hydrogen bond with ASN-230. This in turn may explain its effectiveness against S. schenckii vis-a-vis resistant to ketoconazole. This article is part of a Special Issue entitled: Cytochrome P450 biodiversity and biotechnology, edited by Erika Plettner, Gianfranco Gilardi, Luet Wong, Vlada Urlacher, Jared Goldstone.

Molecular evolutionary dynamics of cytochrome P450 monooxygenases across kingdoms: Special focus on mycobacterial P450s

Since the initial identification of cytochrome P450 monooxygenases (CYPs/P450s), great progress has been made in understanding their structure-function relationship, diversity and application in producing compounds beneficial to humans. However, the molecular evolution of P450s in terms of their dynamics both at protein and DNA levels and functional conservation across kingdoms still needs investigation. In this study, we analyzed 17 598 P450s belonging to 113 P450 families (bacteria -42; fungi -19; plant -28; animal -22; plant and animal -1 and common P450 family -1) and found highly conserved and rapidly evolving P450 families. Results suggested that bacterial P450s, particularly P450s belonging to mycobacteria, are highly conserved both at protein and DNA levels. Mycobacteria possess the highest P450 diversity percentage compared to other microbes and have a high coverage of P450s (≥1%) in their genomes, as found in fungi and plants. Phylogenetic and functional analyses revealed the functional conservation of P450s despite belonging to different biological kingdoms, suggesting the adherence of P450s to their innate function such as their involvement in either generation or oxidation of steroids and structurally related molecules, fatty acids and terpenoids. This study's results offer new understanding of the dynamic structural nature of P450s.