Cytochrome P450 complement (CYPome) of Candida oregonensis, a gut-associated yeast of bark beetle, Dendroctonus rhizophagus

Comparative metabarcoding and biodiversity of gut-associated fungal assemblages of Dendroctonus species (Curculionidae: Scolytinae)

The genus Dendroctonus is a Holarctic taxon composed of 21 nominal species; some of these species are well known in the world as disturbance agents of forest ecosystems. Under the bark of the host tree, these insects are involved in complex and dynamic associations with phoretic ectosymbiotic and endosymbiotic communities. Unlike filamentous fungi and bacteria, the ecological role of yeasts in the bark beetle holobiont is poorly understood, though yeasts were the first group to be recorded as microbial symbionts of these beetles. Our aim was characterize and compare the gut fungal assemblages associated to 14 species of Dendroctonus using the internal transcribed spacer 2 (ITS2) region. A total of 615,542 sequences were recovered yielding 248 fungal amplicon sequence variants (ASVs). The fungal diversity was represented by 4 phyla, 16 classes, 34 orders, 54 families, and 71 genera with different relative abundances among Dendroctonus species. The α-diversity consisted of 32 genera of yeasts and 39 genera of filamentous fungi. An analysis of β-diversity indicated differences in the composition of the gut fungal assemblages among bark beetle species, with differences in species and phylogenetic diversity. A common core mycobiome was recognized at the genus level, integrated mainly by Candida present in all bark beetles, Nakazawaea, Cladosporium, Ogataea, and Yamadazyma. The bipartite networks confirmed that these fungal genera showed a strong association between beetle species and dominant fungi, which are key to maintaining the structure and stability of the fungal community. The functional variation in the trophic structure was identified among libraries and species, with pathotroph-saprotroph-symbiotroph represented at the highest frequency, followed by saprotroph-symbiotroph, and saprotroph only. The overall network suggested that yeast and fungal ASVs in the gut of these beetles showed positive and negative associations among them. This study outlines a mycobiome associated with Dendroctonus nutrition and provides a starting point for future in vitro and omics approaches addressing potential ecological functions and interactions among fungal assemblages and beetle hosts.

Microbial P450 repertoire (P450ome) and its application feasibility in pharmaceutical industry, chemical industry, and environmental protection

Cytochrome P450s (CYPs) are heme-thiolated enzymes that catalyze the oxidation of C-H bonds in a regio- and stereo-selective manner. CYPs are widely present in the biological world. With the completion of more biological genome sequencing, the number and types of P450 enzymes have increased rapidly. P450 in microorganisms is easy to clone and express, rich in catalytic types, and strong in substrate adaptability, which has good application potential. Although the number of P450 enzymes found in microorganisms is huge, the function of most of the microorganism P450s has not been studied, and it contains a large number of excellent biocatalysts to be developed. This review is based on the P450 groups in microorganisms. First, it reviews the distribution of P450 groups in different microbial species, and then studies the application of microbial P450 enzymes in the pharmaceutical industry, chemical industry and environmental pollutant treatment in recent years. And focused on the application fields of P450 enzymes of different families to guide the selection of suitable P450s from the huge P450 library. In view of the current shortcomings of microbial P450 in the application process, the final solution is the most likely to assist the application of P450 enzymes in large-scale, that is, whole cell transformation combined with engineering, fusion P450 combined with immobilization technology.

A Symbiotic Fungus Sistotrema Benefits Blueberry Rejuvenation and Abiotic Stress Tolerance

Blueberry (Vaccinium spp.) rhizosphere microorganisms can significantly increase the absorption area and improve the efficiency of rhizospheric nutrient uptake. However, there has been little research on blueberry rhizosphere microorganisms, especially those that can complement root function deficiency. In this study, we analyzed the rhizosphere fungi of 'O'Neal,' 'Sharpblue,' and 'Premier' blueberry cultivars and found that 'Premier' blueberries showed strong growth potential and relatively high root regulation ability. The dominant symbiotic fungus Sistotrema was correlated with the strong growth of 'Premier' and was directionally screened and isolated based on conserved gene structures and COG function analysis. This fungus was reinoculated onto the roots of 'Gulfcoast' and 'Star' blueberry cultivars. Sistotrema promoted the growth of blueberries and improved their ability to resist stress and grow under adverse conditions, as indicated by maintained or increased chlorophyll content under such conditions. Further analyses showed that Sistotrema has certain functional characteristics such as the ability to dissolve iron in its insoluble form and then release it, to fix nitrogen, and to inhibit nitrification in soil. Thus, it effectively doubled the soil nitrogen content and increased the soluble iron content in soil by 50%. This investigation indicates sistotrema inoculation as an approach to increase blueberry stress tolerance and complete their root nutrition deficiency.

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.

Predicting interactions of the frass-associated yeast Hyphopichia heimii with Olea europaea subsp. cuspidata and twig-boring bark beetles

Bark beetles are destructive insect pests known to form symbioses with different fungal taxa, including yeasts. The aim of this study was to (1) determine the prevalence of the rare yeast Hyphopichia heimii in bark beetle frass from wild olive trees in South Africa and to (2) predict the potential interaction of this yeast with trees and bark beetles. Twenty-eight culturable yeast species were isolated from frass in 35 bark beetle galleries, including representatives of H. heimii from nine samples. Physiological characterization of H. heimii isolates revealed that none was able to degrade complex polymers present in hemicellulose; however, all were able to assimilate sucrose and cellobiose, sugars associated with an arboreal habitat. All isolates were able to produce the auxin indole acetic acid, indicative of a potential symbiosis with the tree. Sterol analysis revealed that the isolates possessed ergosterol quantities ranging from 3.644 ± 0.119 to 13.920 ± 1.230 mg/g dry cell weight, which suggested that H. heimii could serve as a source of sterols in bark beetle diets, as is known for other bark beetle-associated fungi. In addition, gas chromatography-mass spectrometry demonstrated that at least one of the isolates, Hyphopichia heimii CAB 1614, was able to convert the insect pheromone cis-verbenol to the anti-aggregation pheromone verbenone. This indicated that H. heimii could potentially influence beetle behaviour. These results support the contention of a tripartite symbiosis between H. heimii, olive trees, and bark beetles.

Gene excavation and expression analysis of CYP and UGT related to the post modifying stage of gypenoside biosynthesis in Gynostemma pentaphyllum (Thunb.) Makino by comprehensive analysis of RNA and proteome sequencing

Previous studies have revealed that gypenosides produced from Gynostemma pentaphyllum (Thunb.) Makino are mainly dammarane-type triterpenoid saponins with diverse structures and important biological activities, but the mechanism of diversity for gypenoside biosynthesis is still unclear. In this study, a combination of isobaric tags for relative and absolute quantification (iTRAQ) proteome analysis and RNA sequencing transcriptome analysis was performed to identify the proteins and genes related to gypenoside biosynthesis. A total of 3925 proteins were identified by proteomic sequencing, of which 2537 were quantified. Seventeen cytochrome P450 (CYP) and 11 uridine 5’-diphospho-glucuronosyltransferase (UDP-glucuronosyltransferase, UGT) candidate genes involved in the side chain synthesis and modification of gypenosides were found. Seven putative CYPs (CYP71B19, CYP77A3, CYP86A7, CYP86A8, CYP89A2, CYP90A1, CYP94A1) and five putative UGTs (UGT73B4, UGT76B1, UGT74F2, UGT91C1 and UGT91A1) were selected as candidate structural modifiers of triterpenoid saponins, which were cloned for gene expression analysis. Comprehensive analysis of RNA sequencing and proteome sequencing showed that some CYPs and UGTs were found at both the transcription and translation levels. In this study, an expression analysis of 7 CYPs and 5 UGTs that contributed to gypenoside biosynthesis and distribution in G. pentaphyllum was performed, providing consistent results that will inspire more future research on vital genes/proteins involved in gypenoside biosynthesis.

Comparative Sequence Analysis of TRI1 of Fusarium

Trichothecene mycotoxins are a class of secondary metabolites produced by multiple genera of fungi, including certain plant pathogenic Fusarium species. Functional variation in the TRI1 gene produces a novel Type A trichothecene called NX-2 in strains of F. graminearum. Using a bioinformatics approach, a systematic analysis of 52 translated TRI1 sequences of Fusarium species, including five F. graminearum NX-2 producers and four F. graminearum non-NX-2 producers, was conducted to explain the functional difference of TRI1p of FGNX-2. An assessment of several signature motifs of fungal P450s revealed amino acid substitutions in addition to the post-translational N-X-S/T sequons motif, which is indicative of N-linked glycosylation of this TRI1-encoded protein characteristic of NX-2 producers. There was evidence of selection bias, where TRI1 gene sequences were found to be under positive selection and, therefore, under functional constraints. The cumulative amino acid changes in the TRI1p sequences were reflected in the phylogenetic analyses which revealed species-specific clustering with a distinct separation of FGNX-2 from FG-non-NX-2 producers with high bootstrap support. Together, our findings provide insight into the amino acid sequence features responsible for the functional diversification of this TRI1p.

An Overview of Genes From Cyberlindnera americana, a Symbiont Yeast Isolated From the Gut of the Bark Beetle Dendroctonus rhizophagus (Curculionidae: Scolytinae), Involved in the Detoxification Process Using Genome and Transcriptome Data

Bark beetles from Dendroctonus genus promote ecological succession and nutrient cycling in coniferous forests. However, they can trigger outbreaks leading to important economic losses in the forest industry. Conifers have evolved resistance mechanisms that can be toxic to insects but at the same time, bark beetles are capable of overcoming tree barriers and colonize these habitats. In this sense, symbiont yeasts present in the gut of bark beetles have been suggested to play a role in the detoxification process of tree defensive chemicals. In the present study, genes related to this process were identified and their response to a terpene highly toxic to bark beetles and their symbionts was analyzed in the Cyberlindnera americana yeast. The genome and transcriptome of C. americana (ChDrAdgY46) isolated from the gut of Dendroctonus rhizophagus were presented. Genome analysis identified 5752 protein-coding genes and diverse gene families associated with the detoxification process. The most abundant belonged to the Aldo-Keto Reductase Superfamily, ATP-binding cassette Superfamily, and the Major Facilitator Superfamily transporters. The transcriptome analysis of non-α-pinene stimulated and α-pinene stimulated yeasts showed a significant expression of genes belonging to these families. The activities demonstrated by the genes identified as Aryl-alcohol dehydrogenase and ABC transporter under (+)-α-pinene suggest that they are responsible, that C. americana is a dominant symbiont that resists high amounts of monoterpenes inside the gut of bark beetles.

The CYPome of the model xenobiotic-biotransforming fungus Cunninghamella elegans

The fungus Cunninghamella elegans is recognised as a microbial model of mammalian drug metabolism owing to its ability to catabolise xenobiotic compounds in an analogous fashion to animals. Its ability to produce phase I (oxidative) metabolites of drugs is associated with cytochrome P450 (CYP) activity; however, almost nothing is known about these enzymes in the fungus. In this paper we report the in silico analysis of the genome sequence of C. elegans B9769, which contains 32 genes putatively coding for CYPs. Based on their predicted amino acid sequences these were classified as belonging to CYP509, 5203, 5208, 5313, 5210, 61 and 51 families. Reverse transcription-quantitative PCR revealed that the gene coding for CYP5313D1 was significantly upregulated when C. elegans DSM1908 was cultivated in sabouraud dextrose in contrast to its expression in cells grown in Roswell Park Memorial Institute medium. This corresponded to the fungus' xenobiotic biotransformation ability when grown in the two media. Heterologous expression of cyp5313D1 in Pichia pastoris resulted in a recombinant strain that biotransformed flurbiprofen to 4'-hydroxyflurbiprofen, the same metabolite generated by C. elegans cultures. This is the first report of a xenobiotic-biotransforming CYP from this biotechnologically important fungus.

Characterization of the cytochrome P450 monooxygenase genes (P450ome) from the carotenogenic yeast Xanthophyllomyces dendrorhous

BACKGROUND

The cytochromes P450 (P450s) are a large superfamily of heme-containing monooxygenases involved in the oxidative metabolism of an enormous diversity of substrates. These enzymes require electrons for their activity, and the electrons are supplied by NAD(P)H through a P450 electron donor system, which is generally a cytochrome P450 reductase (CPR). The yeast Xanthophyllomyces dendrorhous has evolved an exclusive P450-CPR system that specializes in the synthesis of astaxanthin, a carotenoid with commercial potential. For this reason, the aim of this work was to identify and characterize other potential P450 genes in the genome of this yeast using a bioinformatic approach.

RESULTS

Thirteen potential P450-encoding genes were identified, and the analysis of their deduced proteins allowed them to be classified in ten different families: CYP51, CYP61, CYP5139 (with three members), CYP549A, CYP5491, CYP5492 (with two members), CYP5493, CYP53, CYP5494 and CYP5495. Structural analyses of the X. dendrorhous P450 proteins showed that all of them have a predicted transmembrane region at their N-terminus and have the conserved domains characteristic of the P450s, including the heme-binding region (FxxGxRxCxG); the PER domain, with the characteristic signature for fungi (PxRW); the ExxR motif in the K-helix region and the oxygen-binding domain (OBD) (AGxDTT); also, the characteristic secondary structure elements of all the P450 proteins were identified. The possible functions of these P450s include primary, secondary and xenobiotic metabolism reactions such as sterol biosynthesis, carotenoid synthesis and aromatic compound degradation.

CONCLUSIONS

The carotenogenic yeast X. dendrorhous has thirteen P450-encoding genes having potential functions in primary, secondary and xenobiotic metabolism reactions, including some genes of great interest for fatty acid hydroxylation and aromatic compound degradation. These findings established a basis for future studies about the role of P450s in the carotenogenic yeast X. dendrorhous and their potential biotechnological applications.

Discovery of the molecular mechanisms of the novel chalcone-based Magnaporthe oryzae inhibitor C1 using transcriptomic profiling and co-expression network analysis

Background

In our previous studies, we discovered a series of chalcone-based phytopathogenic fungus inhibitors. However, knowledge of their effects, detailed targets and molecular mechanisms in Magnaporthe oryzae (M. oryzae) remained limited.

Methods

To explore the expression and function of differentially expressed genes in M. oryzae after treatment with compound C1, we analyzed the expression profile of mRNAs using a microarray analysis and GO, KEGG and WGCNA analysis, followed by qRT-PCR and Western blots to validate our findings.

Results

A total of 1013 up-regulated and 995 down-regulated mRNAs were differentially expressed after M. oryzae was treated with C1 compared to those of the control samples. Among these, cytochrome P450, glycylpeptide N-myristoyltransferase (NMT) and peroxisomal membrane protein 4 were identified as the most significant DEGs and were validated by experiments.

Conclusion

In conclusion, our study suggests that the combination of transcriptomic microarray, bioinformatics analysis and weighted gene co-expression networks can be used to predict potential therapeutic targets and to map the pathways regulated by small molecular natural product-like drugs.

Electronic supplementary material

The online version of this article (doi:10.1186/s40064-016-3385-9) contains supplementary material, which is available to authorized users.