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

Hallmarks of Basidiomycete Soft- and White-Rot in Wood-Decay -Omics Data of Two Armillaria Species

Wood-decaying Basidiomycetes are among the most efficient degraders of plant cell walls, making them key players in forest ecosystems, global carbon cycle, and in bio-based industries. Recent insights from -omics data revealed a high functional diversity of wood-decay strategies, especially among the traditional white-rot and brown-rot dichotomy. We examined the mechanistic bases of wood-decay in the conifer-specialists Armillaria ostoyae and Armillaria cepistipes using transcriptomic and proteomic approaches. Armillaria spp. (Fungi, Basidiomycota) include devastating pathogens of temperate forests and saprotrophs that decay wood. They have been discussed as white-rot species, though their response to wood deviates from typical white-rotters. While we observed an upregulation of a diverse suite of plant cell wall degrading enzymes, unlike white-rotters, they possess and express an atypical wood-decay repertoire in which pectinases and expansins are enriched, whereas lignin-decaying enzymes (LDEs) are generally downregulated. This combination of wood decay genes resembles the soft-rot of Ascomycota and appears widespread among Basidiomycota that produce a superficial white rot-like decay. These observations are consistent with ancestral soft-rot decay machinery conserved across asco- and Basidiomycota, a gain of efficient lignin-degrading ability in white-rot fungi and repeated, complete, or partial losses of LDE encoding gene repertoires in brown- and secondarily soft-rot fungi.

Genome sequencing of Rigidoporus microporus provides insights on genes important for wood decay, latex tolerance and interspecific fungal interactions

Fungal plant pathogens remain a serious threat to the sustainable agriculture and forestry, despite the extensive efforts undertaken to control their spread. White root rot disease is threatening rubber tree (Hevea brasiliensis) plantations throughout South and Southeast Asia and Western Africa, causing tree mortality and severe yield losses. Here, we report the complete genome sequence of the basidiomycete fungus Rigidoporus microporus, a causative agent of the disease. Our phylogenetic analysis confirmed the position of R. microporus among the members of Hymenochaetales, an understudied group of basidiomycetes. Our analysis further identified pathogen's genes with a predicted role in the decay of plant cell wall polymers, in the utilization of latex components and in interspecific interactions between the pathogen and other fungi. We also detected putative horizontal gene transfer events in the genome of R. microporus. The reported first genome sequence of a tropical rubber tree pathogen R. microporus should contribute to the better understanding of how the fungus is able to facilitate wood decay and nutrient cycling as well as tolerate latex and utilize resinous extractives.

De novo transcriptomic profiling of the clonal Leymus chinensis response to long-term overgrazing-induced memory

Sheepgrass (Leymus chinensis) is one of the dominant grass species present on typical steppes of the Inner Mongolia Plateau. However, L. chinensis has developed a dwarfing phenotype in response to the stressful habitat in grasslands that are severely degraded due to heavy grazing. The lack of transcriptomic and genomic information has prohibited the understanding of the transgenerational effect on physiological alterations in clonal L. chinensis at the molecular level in response to livestock grazing. To solve this problem, transcriptomic information from the leaves of clonal L. chinensis obtained from overgrazed (GR) and non-grazed (NG) grasslands was studied using a paired-end Illumina HiSeq 2500 sequencing platform. First, despite the influence of grazing being absent during the growth of clonal offspring in our hydroponic experiment, compared with those from the NG group, clonal L. chinensis from the GR group exhibited significant dwarf-type morphological traits. A total of 116,356 unigenes were subsequently generated and assembled de novo, of which 55,541 could be annotated to homologous matches in the NCBI non-redundant (Nr), Swiss-Prot, Clusters of Orthologous Groups (COG), gene ontology (GO), or Kyoto Encyclopedia of Genes and Genomes (KEGG) databases. The expression of 3,341 unigenes significantly differed between the GR group and the NG group with an absolute value of Log₂ ratio ≥ 1. The altered expression of genes involved in defence and immune responses, pathogenic resistance and cell development indicates that livestock grazing induces a transgenerational effect on the growth inhibition of clonal L. chinensis. The results of the present study will provide important large-scale transcriptomic information on L. chinensis. Furthermore, the results facilitated our investigation of grazing-induced transgenerational effects on both the morphological and physiological characteristics of L. chinensis at the molecular levels.

De novo transcriptomic assembly and mRNA expression patterns of Botryosphaeria dothidea infection with mycoviruses chrysovirus 1 (BdCV1) and partitivirus 1 (BdPV1)

BACKGROUND

Pear ring rot, caused by Botryosphaeria species, is responsible for substantial economic losses by causing severe recession of pear tree growth in China. Mycovirus-mediated hypovirulence in plant pathogenic fungi is a crucial biological method to control fungal diseases.

METHODS

We conducted a large-scale and comprehensive transcriptome analysis to identify mRNA in B. dothidea in response to mycovirus. De novo sequencing technology from four constructed libraries of LW-C (Botryosphaeria dothidea chrysovirus 1, BdCV1), LW-P (Botryosphaeria dothidea partitivirus 1, BdPV1), LW-CP (LW-1 strain infection with BdCV1 and BdPV1), and Mock (free virus) was used to investigate and compare gene expression changes in B.dothidea strains infected with mycovirus.

RESULTS

In total, 30,058 Unigenes with an average length of 2128 bp were obtained from 4 libraries of B. dothidea strains. These were annotated to specify their classified function. We demonstrate that mRNAs of B. dothidea strains in response to mycovirus are differentially expressed. In total, 5598 genes were up-regulated and 3298 were down-regulated in the LW-CP group, 4468 were up-regulated and 4291 down-regulated in the LW-C group, and 2590 were up-regulated and 2325 down-regulated in the LW-P group. RT-qPCR was used to validate the expression of 9 selected genes. The B. dothidea transcriptome was more affected by BdCV1 infection than BdPV1. We conducted GO enrichment analysis to characterize gene functions regulated by B. dothidea with mycovirus infection. These involved metabolic process, cellular process, catalytic activity, transporter activity, signaling, and other biological pathways. KEGG function analysis demonstrated that the enriched differentially expressed genes are involved in metabolism, transcription, signal transduction, and ABC transport. mRNA is therefore involved in the interaction between fungi and mycovirus. In addition, changes in differential accumulation levels of cp and RdRp of BdCV1 and BdPV1 in B. dothidea strains were evaluated, revealing that the accumulation of BdCV1 and BdPV1 is related to the phenotype and virulence of B. dothidea strain LW-1.

CONCLUSIONS

The identification and analysis of mRNAs from B. dothidea was first reported at the transcriptome level. Our analysis provides further insight into the interaction of B. dothidea strains infection with chrysovirus 1 (BdCV1) and partitivirus 1 (BdPV1) at the transcriptome level.

Transcriptomic analysis of basidiocarp development in Ustilago maydis (DC) Cda

Previously, we demonstrated that when Ustilago maydis (DC) Cda., a phytopathogenic basidiomycete and the causal agent of corn smut, is grown in the vicinity of maize embryogenic calli in a medium supplemented with the herbicide Dicamba, it developed gastroid-like basidiocarps. To elucidate the molecular mechanisms involved in the basidiocarp development by the fungus, we proceeded to analyze the transcriptome of the process, identifying a total of 2002 and 1064 differentially expressed genes at two developmental stages, young and mature basidiocarps, respectively. Function of these genes was analyzed with the use of different databases. MIPS analysis revealed that in the stage of young basidiocarp, among the ca. two thousand differentially expressed genes, there were some previously described for basidiocarp development in other fungal species. Additional elements that operated at this stage included, among others, genes encoding the transcription factors FOXO3, MIG3, PRO1, TEC1, copper and MFS transporters, and cytochromes P450. During mature basidiocarp development, important up-regulated genes included those encoding hydrophobins, laccases, and ferric reductase (FRE/NOX). The demonstration that a mapkk mutant was unable to form basidiocarps, indicated the importance of the MAPK signaling pathway in this developmental process.