De Novo Transcriptome Analysis of Plant Pathogenic Fungus Myrothecium roridum and Identification of Genes Associated with Trichothecene Mycotoxin Biosynthesis

A reference transcriptome for walnut anthracnose pathogen, Ophiognomonia leptostyla, guides the discovery of candidate virulence genes

Despite the economic losses due to the walnut anthracnose, Ophiognomonia leptostyla is an orphan fungus with respect to genomic resources. In the present study, the transcriptome of O. leptostyla was assembled for the first time. RNA sequencing was conducted for the fungal mycelia grown in a liquid media, and the inoculated leaf samples of walnut with the fungal conidia sampled at 48, 96 and 144 h post inoculation (hpi). The completeness, correctness, and contiguity of the de novo transcriptome assemblies generated with Trinity, Oases, SOAPdenovo-Trans and Bridger were compared to identify a single superior reference assembly. In most of the assessment criteria including N50, Transrate score, number of ORFs with known description in gene bank, the percentage of reads mapped back to the transcript (RMBT), BUSCO score, Swiss-Prot coverage bin and RESM-EVAL score, the Bridger assembly was the superior and thus used as a reference for profiling the O. leptostyla transcriptome in liquid media vs. during walnut infection. The k-means clustering of transcripts resulted in four distinct transcription patterns across the three sampling time points. Most of the detected CAZy transcripts had elevated transcription at 96 hpi that is hypothetically concurrent with the start of intracellular growth. The in-silico analysis revealed 103 candidate effectors of which six were members of Necrosis and Ethylene Inducing Like Protein (NLP) gene family belonging to three distinct k-means clusters. This study provided a complex and temporal pattern of the CAZys and candidate effectors transcription during six days post O. leptostyla inoculation on walnut leaves, introducing a list of candidate virulence genes for validation in future studies.

Shaping effects of rice, wheat, maize, and soybean seedlings on their rhizosphere microbial community

The rhizosphere microbiome plays critical roles in plant growth and is an important interface for resource exchange between plants and the soil environment. Crops at various growing stages, especially the seedling stage, have strong shaping effects on the rhizosphere microbial community, and such community reconstruction will positively feed back to the plant growth. In the present study, we analyzed the variations of bacterial and fungal communities in the rhizosphere of four crop species: rice, soybean, maize, and wheat during successive cultivations (three repeats for the seedling stages) using 16S rRNA gene and internal transcribed spacer (ITS) high-throughput sequencing. We found that the relative abundances of specific microorganisms decreased after different cultivation times, e.g., Sphingomonas, Pseudomonas, Rhodanobacter, and Caulobacter, which have been reported as plant-growth beneficial bacteria. The relative abundances of potential plant pathogenic fungi Myrothecium and Ascochyta increased with the successive cultivation times. The co-occurrence network analysis showed that the bacterial and fungal communities under maize were much more stable than those under rice, soybean, and wheat. The present study explored the characteristics of bacteria and fungi in crop seedling rhizosphere and indicated that the characteristics of indigenous soil flora might determine the plant growth status. Further study will focus on the use of the critical microorganisms to control the growth and yield of specific crops.

The excavation of novel toxin-resistance proteins against trichothecenes toxins in Paramyrothecium roridum

Trichothecene toxins cause serious hazard towards human health and economical crops. However, there are no sufficient molecular strategies to reduce the hazard of trichothecene toxins. Thus it is urgent to exploit novel approaches to control the hazard of trichothecenes. In this study, four trichothecene toxin-resistance genes including mfs1, GNAT1, TRP1 and tri12 in Paramyrothecium roridum were excavated based on genome sequencing results, and then expressed in toxin-sensitive Saccharomyces cerevisiae BJ5464, the toxin resistance genes pdr5, pdr10 and pdr15 of which were firstly knocked out simultaneously by the introduction of TAA stop codon employing CRISPR/Cas9 system. Therefore, three novel hazardous toxin-resistance genes mfs1, GNAT1, TRP1 in P. roridum were firstly excavated by the co-incubation of DON toxin and toxin resistant genes-containing BJ5464 strains. The in vitro function and properties of novel toxin-resistance genes coding proteins including GNAT1, MFS1 and TRP1 were identified by heterologous expression and cellular location analysis as well as in vitro biochemical reaction. The excavation of novel trichothecene toxin-resistance genes provide novel molecular clues for controlling the harm of trichothecenes, meanwhile, this study will also pave a new way for the yield improvement of trichothecenes by heterologous expression to facilitate the development of trichothecenes as anti-tumor lead compounds.

A high-throughput RNA-Seq approach to elucidate the transcriptional response of Piriformospora indica to high salt stress

Piriformospora indica, a root endophytic fungus, augments plant nutrition and productivity as well as protects plants against pathogens and abiotic stresses. High salinity is a major problem faced by plants as well as by microbes. Until now, the precise mechanism of salt stress tolerance in P. indica has remained elusive. In this study, the transcriptomes of control and salt-treated (0.5 M NaCl) P. indica were sequenced via the RNA-seq approach. A total of 30,567 transcripts and 15,410 unigenes for P. indica were obtained from 7.3 Gb clean reads. Overall 661 differentially expressed genes (DEGs) between control and treated samples were retrieved. Gene ontology (GO) and EuKaryotic Orthologous Groups (KOG) enrichments revealed that DEGs were specifically involved in metabolic and molecular processes, such as "response to salt stress", "oxidoreductase activity", "ADP binding", "translation, ribosomal structure and biogenesis", "cytoskeleton", and others. The unigenes involved in "cell wall integrity", "sterol biosynthesis", and "oxidative stress" such as Rho-type GTPase, hydroxymethylglutaryl-CoA synthase, and thioredoxin peroxidase were up-regulated in P. indica subjected to salt stress. The salt-responsive DEGs have shown that they might have a potential role in salt stress regulation. Our study on the salt-responsive DEGs established a foundation for the elucidation of molecular mechanisms related to P. indica stress adaptation and a future reference for comparative functional genomics studies of biotechnologically important fungal species.

Genetic bases for variation in structure and biological activity of trichothecene toxins produced by diverse fungi

Trichothecenes are sesquiterpene toxins produced by diverse but relatively few fungal species in at least three classes of Ascomycetes: Dothideomycetes, Eurotiomycetes, and Sordariomycetes. Approximately 200 structurally distinct trichothecene analogs have been described, but a given fungal species typically produces only a small subset of analogs. All trichothecenes share a core structure consisting of a four-ring nucleus known as 12,13-epoxytrichothec-9-ene. This structure can be substituted at various positions with hydroxyl, acyl, or keto groups to give rise to the diversity of trichothecene structures that has been described. Over the last 30 years, the genetic and biochemical pathways required for trichothecene biosynthesis in several species of the fungi Fusarium and Trichoderma have been elucidated. In addition, phylogenetic and functional analyses of trichothecene biosynthetic (TRI) genes from fungi in multiple genera have provided insights into how acquisition, loss, and changes in functions of TRI genes have given rise to the diversity of trichothecene structures. These analyses also suggest both divergence and convergence of TRI gene function during the evolutionary history of trichothecene biosynthesis. What has driven trichothecene structural diversification remains an unanswered question. However, insight into the role of trichothecenes in plant pathogenesis of Fusarium species and into plant glucosyltransferases that detoxify the toxins by glycosylating them point to a possible driver. Because the glucosyltransferases can have substrate specificity, changes in trichothecene structures produced by a fungus could allow it to evade detoxification by the plant enzymes. Thus, it is possible that advantages conferred by evading detoxification have contributed to trichothecene structural diversification. KEY POINTS : • TRI genes have evolved by diverse processes: loss, acquisition and changes in function. • Some TRI genes have acquired the same function by convergent evolution. • Some other TRI genes have evolved divergently to have different functions. • Some TRI genes were acquired or resulted from diversification in function of other genes. • Substrate specificity of plant glucosyltransferases could drive trichothecene diversity.

Biocontrol of Root Diseases and Growth Promotion of the Tuberous Plant Aconitum carmichaelii Induced by Actinomycetes Are Related to Shifts in the Rhizosphere Microbiota

Soil Actinomycetes have been used as biocontrol agents against soil-borne plant diseases, yet little is known about their effects on the structure of the rhizosphere microbiota and the long-term effects on crop yield and disease intensity after the application of Actinomycetes is stopped. Here, we conducted 3-year plot experiments to investigate the roles of two Actinomycetes strains (Streptomyces pactum Act12 and Streptomyces rochei D74) in the biocontrol of soil-borne root diseases and growth promotion of monkhood (Aconitum carmichaelii). We also examined their long-term effects after soil application of a mixed Actinomycetes preparation (spore powder) was completed. High-throughput sequencing was used to analyze shifts in the rhizosphere microbiota. The antifungal activity and root colonization ability of the two Actinomycetes were also tested. Disease severity of southern blight and root rot decreased following application of the Actinomycetes preparation, whereas biomass yield of tubers increased compared with the control group. Significant effects of disease control and plant growth promotion were also observed after application was stopped. The Actinomycetes preparation induced marked increases in the abundance of beneficial microbes and decreases in the abundance of harmful microbes in rhizosphere soil. Adding cell-free culture filtrates of both strains Act12 and D74 inhibited the growth of fungal pathogens capable of causing southern blight (Sclerotium rolfsii) and root rot (Fusarium oxysporum) in A. carmichaelii. A GFP-labeled strain was used to show that D74 can colonize roots of A. carmichaelii. In conclusion, a preparation of two Actinomycetes plays a role in the biocontrol of root diseases and growth promotion of A. carmichaelii by inhibiting pathogen growth and shaping the rhizosphere microbiota.

Transcriptome analysis of the fungal pathogen Rosellinia necatrix during infection of a susceptible avocado rootstock identifies potential mechanisms of pathogenesis

Background

White root rot disease caused by Rosellinia necatrix is one of the most important threats affecting avocado productivity in tropical and subtropical climates. Control of this disease is complex and nowadays, lies in the use of physical and chemical methods, although none have proven to be fully effective. Detailed understanding of the molecular mechanisms underlying white root rot disease has the potential of aiding future developments in disease resistance and management. In this regard, this study used RNA-Seq technology to compare the transcriptomic profiles of R. necatrix during infection of susceptible avocado ‘Dusa’ roots with that obtained from the fungus cultured in rich medium.

Results

The transcriptomes from three biological replicates of R. necatrix colonizing avocado roots (RGA) and R. necatrix growing on potato dextrose agar media (RGPDA) were analyzed using Illumina sequencing. A total of 12,104 transcripts were obtained, among which 1937 were differentially expressed genes (DEG), 137 exclusively expressed in RGA and 160 in RGPDA. During the root infection process, genes involved in the production of fungal toxins, detoxification and transport of toxic compounds, hormone biosynthesis, gene silencing and plant cell wall degradation were overexpressed. Interestingly, 24 out of the 137 contigs expressed only during R. necatrix growth on avocado roots, were predicted as candidate effector proteins (CEP) with a probability above 60%. The PHI (Pathogen Host Interaction) database revealed that three of the R. necatrix CEP showed homology with previously annotated effectors, already proven experimentally via pathogen-host interaction.

Conclusions

The analysis of the full-length transcriptome of R. necatrix during the infection process is suggesting that the success of this fungus to infect roots of diverse crops might be attributed to the production of different compounds which, singly or in combination, interfere with defense or signaling mechanisms shared among distinct plant families. The transcriptome analysis of R. necatrix during the infection process provides useful information and facilitates further research to a more in -depth understanding of the biology and virulence of this emergent pathogen. In turn, this will make possible to evolve novel strategies for white root rot management in avocado.

Heterologous Biosynthesis of the Fungal Sesquiterpene Trichodermol in Saccharomyces cerevisiae

Trichodermol, a fungal sesquiterpene derived from the farnesyl diphosphate pathway, is the biosynthetic precursor for trichodermin, a member of the trichothecene class of fungal toxins produced mainly by the genera of Trichoderma and Fusarium. Trichodermin is a promising candidate for the development of fungicides and antitumor agents due to its significant antifungal and cytotoxic effects. It can also serve as a scaffold to generate new congeners for structure-activity relationship (SAR) study. We reconstructed the biosynthetic pathway of trichodermol in Saccharomyces cerevisiae BY4741, and investigated the effect of produced trichodermol on the host by de novo RNA sequencing (RNA-Seq) and quantitative Real-time PCR analyses. Co-expression of pESC::FgTRI5 using plasmid pLLeu-tHMGR-UPC2.1 led to trichodiene production of 683 μg L^-1, while integration of only the codon-optimized FgTRI5 into the chromosome of yeast improved the production to 6,535 μg L^-1. Subsequent expression of the codon-optimized cytochrome P450 monooxygenase encoding genes, TaTRI4 and TaTRI11, resulted in trichodermol, with an estimated titer of 252 μg L^-1 at shake flask level. RNA-Seq and qPCR analyses revealed that the produced trichodermol downregulated the expression of the genes involved in ergosterol biosynthesis, but significantly upregulated the expression of PDR5 related to membrane transport pathway in S. cerevisiae. Collectively, we achieved the first heterologous biosynthesis of trichodermol by reconstructing its biosynthetic pathway in yeast, and the reconstructed pathway will serve as a platform to generate trichodermin analogs as potential candidates for agrochemicals and anticancer agents through further optimizations.

De Novo Transcriptome Sequencing of the Deep-Sea-Derived Fungus Dichotomomyces cejpii and Analysis of Gliotoxin Biosynthesis Genes

Gliotoxin, produced by fungi, is an epipolythiodioxopiperazine (ETP) toxin with bioactivities such as anti-liver fibrosis, antitumor, antifungus, antivirus, antioxidation, and immunoregulation. Recently, cytotoxic gliotoxins were isolated from a deep-sea-derived fungus, Dichotomomyces cejpii. However, the biosynthetic pathway for gliotoxins in D. cejpii remains unclear. In this study, the transcriptome of D. cejpii was sequenced using an Illumina Hiseq 2000. A total of 19,125 unigenes for D. cejpii were obtained from 9.73 GB of clean reads. Ten genes related to gliotoxin biosynthesis were annotated. The expression levels of gliotoxin-related genes were detected through quantitative real-time polymerase chain reaction (qRT-PCR). The GliG gene, encoding a glutathione S-transferase (DC-GST); GliI, encoding an aminotransferase (DC-AI); and GliO, encoding an aldehyde reductase (DC-AR), were cloned and expressed, purified, and characterized. The results suggested the important roles of DC-GST, DC-AT, and DC-AR in the biosynthesis of gliotoxins. Our study on the genes related to gliotoxin biosynthesis establishes a molecular foundation for the wider application of gliotoxins from D. cejpii in the biomedical industry in the future.

Evolution of structural diversity of trichothecenes, a family of toxins produced by plant pathogenic and entomopathogenic fungi

Trichothecenes are a family of terpenoid toxins produced by multiple genera of fungi, including plant and insect pathogens. Some trichothecenes produced by the fungus Fusarium are among the mycotoxins of greatest concern to food and feed safety because of their toxicity and frequent occurrence in cereal crops, and trichothecene production contributes to pathogenesis of some Fusarium species on plants. Collectively, fungi produce over 150 trichothecene analogs: i.e., molecules that share the same core structure but differ in patterns of substituents attached to the core structure. Here, we carried out genomic, phylogenetic, gene-function, and analytical chemistry studies of strains from nine fungal genera to identify genetic variation responsible for trichothecene structural diversity and to gain insight into evolutionary processes that have contributed to the variation. The results indicate that structural diversity has resulted from gain, loss, and functional changes of trichothecene biosynthetic (TRI) genes. The results also indicate that the presence of some substituents has arisen independently in different fungi by gain of different genes with the same function. Variation in TRI gene duplication and number of TRI loci was also observed among the fungi examined, but there was no evidence that such genetic differences have contributed to trichothecene structural variation. We also inferred ancestral states of the TRI cluster and trichothecene biosynthetic pathway, and proposed scenarios for changes in trichothecene structures during divergence of TRI cluster homologs. Together, our findings provide insight into evolutionary processes responsible for structural diversification of toxins produced by pathogenic fungi.

Toxins produced by pathogens can contribute to infection and/or colonization of hosts. Some toxins consist of a family of metabolites with similar but distinct chemical structures. This structural variation can affect biological activity, which in turn likely contributes to adaptation to different environments, including to different hosts. Trichothecene toxins consist of over 150 structurally distinct molecules produced by certain fungi, including some plant and insect pathogens. In multiple systems that have been examined, trichothecenes contribute to pathogenesis on plants. To elucidate the evolutionary processes that have given rise to trichothecene structural variation, we conducted comparative analyses of nine fungal genera, most of which produce different trichothecene structures. Using genomic, molecular biology, phylogenetic, and analytical chemistry approaches, we obtained evidence that trichothecene structural variation has arisen primarily from gain, loss, and functional changes of trichothecene biosynthetic genes. Our results also indicate that some structural changes have arisen independently in different fungi. Our findings provide insight into genetic and biochemical changes that can occur in toxin biosynthetic pathways as fungi with the pathways adapt to different environmental conditions.