Developmental Dynamics of Long Noncoding RNA Expression during Sexual Fruiting Body Formation in Fusarium graminearum

Phylogenetic and functional analyses of N6-methyladenosine RNA methylation factors in the wheat scab fungus Fusarium graminearum

In eukaryotes, N6-methyladenosine (m6A) RNA modification plays a crucial role in governing the fate of RNA molecules and has been linked to various developmental processes. However, the phyletic distribution and functions of genetic factors responsible for m6A modification remain largely unexplored in fungi. To get insights into the evolution of m6A machineries, we reconstructed global phylogenies of potential m6A writers, readers, and erasers in fungi. Substantial copy number variations were observed, ranging from up to five m6A writers in early-diverging fungi to a single copy in the subphylum Pezizomycotina, which primarily comprises filamentous fungi. To characterize m6A factors in a phytopathogenic fungus Fusarium graminearum, we generated knockout mutants lacking potential m6A factors including the sole m6A writer MTA1. However, the resulting knockouts did not exhibit any noticeable phenotypic changes during vegetative and sexual growth stages. As obtaining a homozygous knockout lacking MTA1 was likely hindered by its essential role, we generated MTA1-overexpressing strains (MTA1-OE). The MTA1-OE5 strain showed delayed conidial germination and reduced hyphal branching, suggesting its involvement during vegetative growth. Consistent with these findings, the expression levels of MTA1 and a potential m6A reader YTH1 were dramatically induced in germinating conidia, followed by the expression of potential m6A erasers at later vegetative stages. Several genes including transcription factors, transporters, and various enzymes were found to be significantly upregulated and downregulated in the MTA1-OE5 strain. Overall, our study highlights the functional importance of the m6A methylation during conidial germination in F. graminearum and provides a foundation for future investigations into m6A modification sites in filamentous fungi.IMPORTANCEN6-methyladenosine (m6A) RNA methylation is a reversible posttranscriptional modification that regulates RNA function and plays a crucial role in diverse developmental processes. This study addresses the knowledge gap regarding phyletic distribution and functions of m6A factors in fungi. The identification of copy number variations among fungal groups enriches our knowledge regarding the evolution of m6A machinery in fungi. Functional characterization of m6A factors in a phytopathogenic filamentous fungus Fusarium graminearum provides insights into the essential role of the m6A writer MTA1 in conidial germination and hyphal branching. The observed effects of overexpressing MTA1 on fungal growth and gene expression patterns of m6A factors throughout the life cycle of F. graminearum further underscore the importance of m6A modification in conidial germination. Overall, this study significantly advances our understanding of m6A modification in fungi, paving the way for future research into its roles in filamentous growth and potential applications in disease control.

A non-pheromone GPCR is essential for meiosis and ascosporogenesis in the wheat scab fungus

Meiosis is essential for generating genetic diversity and sexual spores, but the regulation of meiosis and ascosporogenesis is not clear in filamentous fungi, in which dikaryotic and diploid cells formed inside fruiting bodies are not free living and independent of pheromones or pheromone receptors. In this study, Gia1, a non-pheromone GPCR (G protein-coupled receptor) with sexual-specific expression in Fusarium graminearum, is found to be essential for ascosporogenesis. The gia1 mutant was normal in perithecium development, crozier formation, and karyogamy but failed to undergo meiosis, which could be partially rescued by a dominant active mutation in GPA1 and activation of the Gpmk1 pathway. GIA1 orthologs have conserved functions in regulating meiosis and ascosporogenesis in Sordariomycetes. GIA1 has a paralog, GIP1, in F. graminearum and other Hypocreales species which is essential for perithecium formation. GIP1 differed from GIA1 in expression profiles and downstream signaling during sexual reproduction. Whereas the C-terminal tail and IR3 were important for intracellular signaling, the N-terminal region and EL3 of Gia1 were responsible for recognizing its ligand, which is likely a protein enriched in developing perithecia, particularly in the gia1 mutant. Taken together, these results showed that GIA1 encodes a non-pheromone GPCR that regulates the entry into meiosis and ascosporogenesis via the downstream Gpmk1 MAP kinase pathway in F. graminearum and other filamentous ascomycetes.

The Sordariomycetes: an expanding resource with Big Data for mining in evolutionary genomics and transcriptomics

Advances in genomics and transcriptomics accompanying the rapid accumulation of omics data have provided new tools that have transformed and expanded the traditional concepts of model fungi. Evolutionary genomics and transcriptomics have flourished with the use of classical and newer fungal models that facilitate the study of diverse topics encompassing fungal biology and development. Technological advances have also created the opportunity to obtain and mine large datasets. One such continuously growing dataset is that of the Sordariomycetes, which exhibit a richness of species, ecological diversity, economic importance, and a profound research history on amenable models. Currently, 3,574 species of this class have been sequenced, comprising nearly one-third of the available ascomycete genomes. Among these genomes, multiple representatives of the model genera Fusarium, Neurospora, and Trichoderma are present. In this review, we examine recently published studies and data on the Sordariomycetes that have contributed novel insights to the field of fungal evolution via integrative analyses of the genetic, pathogenic, and other biological characteristics of the fungi. Some of these studies applied ancestral state analysis of gene expression among divergent lineages to infer regulatory network models, identify key genetic elements in fungal sexual development, and investigate the regulation of conidial germination and secondary metabolism. Such multispecies investigations address challenges in the study of fungal evolutionary genomics derived from studies that are often based on limited model genomes and that primarily focus on the aspects of biology driven by knowledge drawn from a few model species. Rapidly accumulating information and expanding capabilities for systems biological analysis of Big Data are setting the stage for the expansion of the concept of model systems from unitary taxonomic species/genera to inclusive clusters of well-studied models that can facilitate both the in-depth study of specific lineages and also investigation of trait diversity across lineages. The Sordariomycetes class, in particular, offers abundant omics data and a large and active global research community. As such, the Sordariomycetes can form a core omics clade, providing a blueprint for the expansion of our knowledge of evolution at the genomic scale in the exciting era of Big Data and artificial intelligence, and serving as a reference for the future analysis of different taxonomic levels within the fungal kingdom.

Comparative Transcriptomics of Fusarium graminearum and Magnaporthe oryzae Spore Germination Leading up To Infection

For fungal plant pathogens, the germinating spore provides the first interaction with the host. Spore germlings move across the plant surface and use diverse penetration strategies for ingress into plant surfaces. Penetration strategies include pressurized melanized appressoria, which facilitate physically punching through the plant cuticle, and nonmelanized appressoria, which penetrate with the help of enzymes or cuticular damage to breach the plant surface. Two well-studied plant pathogens, Fusarium graminearum and Magnaporthe oryzae, are typical of these two modes of penetration. We applied comparative transcriptomics to Fusarium graminearum and Magnaporthe oryzae to characterize the genetic programming of the early host-pathogen interface. Four sequential stages of development following spore localization on the plant surface, from spore swelling to appressorium formation, were sampled for each species on culture medium and on barley sheaths, and transcriptomic analyses were performed. Gene expression in the prepenetration stages in both species and under both conditions was similar. In contrast, gene expression in the final stage was strongly influenced by the environment. Appressorium formation involved the greatest number of differentially expressed genes. Laser-dissection microscopy was used to perform detailed transcriptomics of initial infection points by F. graminearum. These analyses revealed new and important aspects of early fungal ingress in this species. Expression of the trichothecene genes involved in biosynthesis of deoxynivalenol by F. graminearum implies that toxisomes are not fully functional until after penetration and indicates that deoxynivalenol is not essential for penetration under our conditions. The use of comparative gene expression of divergent fungi promises to advance highly effective targets for antifungal strategies. IMPORTANCE Fusarium graminearum and Magnaporthe oryzae are two of the most important pathogens of cereal grains worldwide. Despite years of research, strong host resistance has not been identified for F. graminearum, so other methods of control are essential. The pathogen takes advantage of multiple entry points to infect the host, including breaches in the florets due to senescence of flower parts and penetration of the weakened trichome bases to breach the epidermis. In contrast, M. oryzae directly punctures leaves that it infects, and resistant cultivars have been characterized. The threat of either pathogen causing a major disease outbreak is ever present. Comparative transcriptomics demonstrated its potential to reveal novel and effective disease prevention strategies that affect the initial stages of disease. Shedding light on the basis of this diversity of infection strategies will result in development of increasingly specific control strategies.

Identification of Long Non-Coding RNAs and Their Target Genes from Mycelium and Primordium in Model Mushroom Schizophyllum commune

Schizophyllum commune has emerged as the most promising model mushroom to study developmental stages (mycelium, primordium), which are two primary processes of fruit body development. Long non-coding RNA (lncRNA) has been proved to participate in fruit development and sex differentiation in fungi. However, potential lncRNAs have not been identified in S. commune from mycelium to primordium developmental stages. In this study, lncRNA-seq was performed in S. commune and 61.56 Gb clean data were generated from mycelium and primordium developmental stages. Furthermore, 191 lncRNAs had been obtained and a total of 49 lncRNAs were classified as differently expressed lncRNAs. Additionally, 26 up-regulated differently expressed lncRNAs and 23 down-regulated between mycelium and primordia libraries were detected. Further, Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis showed that differentially expressed lncRNAs target genes from the MAPK pathway, phosphatidylinositol signal, ubiquitin-mediated proteolysis, autophagy, and cell cycle. This study provides a new resource for further research on the relationship between lncRNA and two developmental stages (mycelium, primordium) in S. commune.

Transcriptional Divergence Underpinning Sexual Development in the Fungal Class Sordariomycetes

Gene expression divergence through evolutionary processes is thought to be important for achieving programmed development in multicellular organisms. To test this premise in filamentous fungi, we investigated transcriptional profiles of 3,942 single-copy orthologous genes (SCOGs) in five related sordariomycete species that have morphologically diverged in the formation of their flask-shaped perithecia. We compared expression of the SCOGs to inferred gene expression levels of the most recent common ancestor of the five species, ranking genes from their largest increases to smallest increases in expression during perithecial development in each of the five species. We found that a large proportion of the genes that exhibited evolved increases in gene expression were important for normal perithecial development in Fusarium graminearum. Many of these genes were previously uncharacterized, encoding hypothetical proteins without any known functional protein domains. Interestingly, the developmental stages during which aberrant knockout phenotypes appeared largely coincided with the elevated expression of the deleted genes. In addition, we identified novel genes that affected normal perithecial development in Magnaporthe oryzae and Neurospora crassa, which were functionally and transcriptionally diverged from the orthologous counterparts in F. graminearum. Furthermore, comparative analysis of developmental transcriptomes and phylostratigraphic analysis suggested that genes encoding hypothetical proteins are generally young and transcriptionally divergent between related species. This study provides tangible evidence of shifts in gene expression that led to acquisition of novel function of orthologous genes in each lineage and demonstrates that several genes with hypothetical function are crucial for shaping multicellular fruiting bodies. IMPORTANCE The fungal class Sordariomycetes includes numerous important plant and animal pathogens. It also provides model systems for studying fungal fruiting body development, as its members develop fruiting bodies with a few well-characterized tissue types on common growth media and have rich genomic resources that enable comparative and functional analyses. To understand transcriptional divergence of key developmental genes between five related sordariomycete fungi, we performed targeted knockouts of genes inferred to have evolved significant upward shifts in expression. We found that many previously uncharacterized genes play indispensable roles at different stages of fruiting body development, which have undergone transcriptional activation in specific lineages. These novel genes are predicted to be phylogenetically young and tend to be involved in lineage- or species-specific function. Transcriptional activation of genes with unknown function seems to be more frequent than ever thought, which may be crucial for rapid adaption to changing environments for successful sexual reproduction.

Efficient and specific DNA oligonucleotide rRNA probe-based rRNA removal in Talaromyces marneffei

Emerging evidence showed that lncRNAs play important roles in a wide range of biological processes of fungi such as Saccharomyces cerevisiae. However, systemic identification of lncRNAs in non-model fungi is a challenging task as the efficiency of rRNA removal has been proved to be affected by mismatches of universal rRNA-targeting probes of commercial kits, which forces deeper sequencing depth and increases costs. Here, we developed a low-cost and simple rRNA depletion method (rProbe) that could efficiently remove more than 99% rRNA in both yeast and mycelium samples of Talaromyces marneffei. The efficiency and robustness of rProbe were demonstrated to outperform the Illumina Ribo-Zero kit. Using rProbe RNA-seq, we identified 115 differentially expressed lncRNAs and constructed lncRNA-mRNA co-expression network related to dimorphic switch of T. marneffei. Our rRNA removal method has the potential to be a useful tool to explore non-coding transcriptomes of non-model fungi by adjusting rRNA probe sequences species specifically.

Long Non-Coding RNAs in Cryptococcus neoformans: Insights Into Fungal Pathogenesis

Long non-coding RNAs (lncRNAs) are highly expressed and can modulate multiple cellular processes including transcription, splicing, translation, and many diverse signaling events. LncRNAs can act as sponges for miRNAs, RNA and DNA binding proteins, functioning as competitive endogenous RNAs. The contribution of lncRNAs to microbial pathogenesis is largely neglected in eukaryotic pathogens despite the abundance of RNA sequencing datasets encompassing conditions of stress, gene deletions and conditions that mimic the host environment. The human fungal pathogen Cryptococcus neoformans encodes 6975 (84%) protein-coding and 1359 (16%) non-protein-coding RNAs, of which 1182 (14.2%) are lncRNAs defined by a threshold of greater than 200 nucleotides in length. Here, we discuss the current state of knowledge in C. neoformans lncRNA biology. Utilizing existing RNA seq datasets, we examine trends in lncRNA expression and discuss potential implications for pathogenesis.

Sexual Crossing, Chromosome-Level Genome Sequences, and Comparative Genomic Analyses for the Medicinal Mushroom Taiwanofungus Camphoratus (Syn. Antrodia Cinnamomea, Antrodia Camphorata)

Taiwanofungus camphoratus mushrooms are a complementary and alternative medicine for hangovers, cancer, hypertension, obesity, diabetes, and inflammation. Though Taiwanofungus camphoratus has attracted considerable biotechnological and pharmacological attention, neither classical genetic nor genomic approaches have been properly established for it. We isolated four sexually competent monokaryons from two T. camphoratus dikaryons used for the commercial cultivation of orange-red (HC1) and milky-white (SN1) mushrooms, respectively. We also sequenced, annotated, and comparatively analyzed high-quality and chromosome-level genome sequences of these four monokaryons. These genomic resources represent a valuable basis for understanding the biology, evolution, and secondary metabolite biosynthesis of this economically important mushrooms. We demonstrate that T. camphoratus has a tetrapolar mating system and that HC1 and SN1 represent two intraspecies isolates displaying karyotypic variation. Compared with several edible mushroom model organisms, T. camphoratus underwent a significant contraction in the gene family and individual gene numbers, most notably for plant, fungal, and bacterial cell-wall-degrading enzymes, explaining why T. camphoratus mushrooms are rare in natural environments, are difficult and time-consuming to artificially cultivate, and are susceptible to fungal and bacterial infections. Our results lay the foundation for an in-depth T. camphoratus study, including precise genetic manipulation, improvements to mushroom fruiting, and synthetic biology applications for producing natural medicinal products. IMPORTANCETaiwanofungus camphoratus (Tc) is a basidiomycete fungus that causes brown heart rot of the aromatic tree Cinnamomum kanehirae. The Tc fruiting bodies have been used to treat hangovers, abdominal pain, diarrhea, hypertension, and other diseases first by aboriginal Taiwanese and later by people in many countries. To establish classical genetic and genomic approaches for this economically important medicinal mushroom, we first isolated and characterized four sexually competent monokaryons from two dikaryons wildly used for commercial production of Tc mushrooms. We applied PacBio single molecule, real-time sequencing technology to determine the near-completed genome sequences of four monokaryons. These telomere-to-telomere and gapless haploid genome sequences reveal all genomic variants needed to be studied and discovered, including centromeres, telomeres, retrotransposons, mating type loci, biosynthetic, and metabolic gene clusters. Substantial interspecies diversities are also discovered between Tc and several other mushroom model organisms, including Agrocybe aegerita, Coprinopsis cinerea, and Schizophyllum commune, and Ganoderma lucidum.

Integrative Analysis of lncRNA-mRNA Co-expression Provides Novel Insights Into the Regulation of Developmental Transitions in Female Varroa destructor

Varroa destructor is a major pathogenic driver of the Western honeybee colony losses globally. Understanding the developmental regulation of V. destructor is critical to develop effective control measures. Development is a complex biological process regulated by numerous genes and long non-coding RNAs (lncRNAs); however, the underlying regulation of lncRNAs in the development of V. destructor remains unknown. In this study, we analyzed the RNA sequencing (RNA-Seq) data derived from the four stages of female V. destructor in the reproductive phase (i.e., egg, protonymph, deutonymph, and adult). The identified differentially expressed mRNAs and lncRNAs exhibited a stage-specific pattern during developmental transitions. Further functional enrichment established that fat digestion and absorption, ATP-binding cassette (ABC) transporters, mitogen-activated protein kinase (MAPK) signaling pathway, and ubiquitin-proteasome pathway play key roles in the maturation of female V. destructor. Moreover, the lncRNAs and mRNAs of some pivotal genes were significantly upregulated at the deutonymph stage, such as cuticle protein 65/6.4/63/38 and mucin 5AC, suggesting that deutonymph is the key stage of metamorphosis development and pathogen resistance acquisition for female V. destructor. Our study provides novel insights into a foundational understanding of V. destructor biology.

Genome-wide profiling of long non-coding RNA of the rice blast fungus Magnaporthe oryzae during infection

Background

Long non-coding RNAs (lncRNAs) play essential roles in developmental processes and disease development at the transcriptional and post-transcriptional levels across diverse taxa. However, only few studies have profiled fungal lncRNAs in a genome-wide manner during host infection.

Results

Infection-associated lncRNAs were identified using lncRNA profiling over six stages of host infection (e.g., vegetative growth, pre-penetration, biotrophic, and necrotrophic stages) in the model pathogenic fungus, Magnaporthe oryzae. We identified 2,601 novel lncRNAs, including 1,286 antisense lncRNAs and 980 intergenic lncRNAs. Among the identified lncRNAs, 755 were expressed in a stage-specific manner and 560 were infection-specifically expressed lncRNAs (ISELs). To decipher the potential roles of lncRNAs during infection, we identified 365 protein-coding genes that were associated with 214 ISELs. Analysis of the predicted functions of these associated genes suggested that lncRNAs regulate pathogenesis-related genes, including xylanases and effectors.

Conclusions

The ISELs and their associated genes provide a comprehensive view of lncRNAs during fungal pathogen-plant interactions. This study expands new insights into the role of lncRNAs in the rice blast fungus, as well as other plant pathogenic fungi.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08380-4.

Gene age shapes the transcriptional landscape of sexual morphogenesis in mushroom forming fungi (Agaricomycetes)

Multicellularity has been one of the most important innovations in the history of life. The role of gene regulatory changes in driving transitions to multicellularity is being increasingly recognized; however, factors influencing gene expression patterns are poorly known in many clades. Here, we compared the developmental transcriptomes of complex multicellular fruiting bodies of eight Agaricomycetes and Cryptococcus neoformans, a closely related human pathogen with a simple morphology. In-depth analysis in Pleurotus ostreatus revealed that allele-specific expression, natural antisense transcripts, and developmental gene expression, but not RNA editing or a ‘developmental hourglass,’ act in concert to shape its transcriptome during fruiting body development. We found that transcriptional patterns of genes strongly depend on their evolutionary ages. Young genes showed more developmental and allele-specific expression variation, possibly because of weaker evolutionary constraint, suggestive of nonadaptive expression variance in fruiting bodies. These results prompted us to define a set of conserved genes specifically regulated only during complex morphogenesis by excluding young genes and accounting for deeply conserved ones shared with species showing simple sexual development. Analysis of the resulting gene set revealed evolutionary and functional associations with complex multicellularity, which allowed us to speculate they are involved in complex multicellular morphogenesis of mushroom fruiting bodies.

lncRsp1, a long noncoding RNA, influences Fgsp1 expression and sexual reproduction in Fusarium graminearum

Long noncoding RNAs (lncRNAs) are crucial regulators of gene expression in many biological processes, but their biological functions remain largely unknown, especially in fungi. Fusarium graminearum is an important pathogen that causes the destructive disease Fusarium head blight (FHB) or head scab disease on wheat and barley. In our previous RNA sequencing (RNA-Seq) study, we discovered that lncRsp1 is an lncRNA that is located +99 bp upstream of a putative sugar transporter gene, Fgsp1, with the same transcription direction. Functional studies revealed that ΔlncRsp1 and ΔFgsp1 were normal in growth and conidiation but had defects in ascospore discharge and virulence on wheat coleoptiles. Moreover, lncRsp1 and Fgsp1 were shown to negatively regulate the expression of several deoxynivalenol (DON) biosynthesis genes, TRI4, TRI5, TRI6, and TRI13, as well as DON production. Further analysis showed that the overexpression of lncRsp1 enhanced the ability of ascospore release and increased the mRNA expression level of the Fgsp1 gene, while lncRsp1-silenced strains reduced ascospore discharge and inhibited Fgsp1 expression during the sexual reproduction stage. In addition, the lncRsp1 complementary strains lncRsp1-LC-1 and lncRsp1-LC-2 restored ascospore discharge to the level of the wild-type strain PH-1. Taken together, our results reveal the distinct and specific functions of lncRsp1 and Fgsp1 in F. graminearum and principally demonstrate that lncRsp1 can affect the release of ascospores by regulating the expression of Fgsp1.

The long non-coding RNA landscape of Candida yeast pathogens

Long non-coding RNAs (lncRNAs) constitute a poorly studied class of transcripts with emerging roles in key cellular processes. Despite efforts to characterize lncRNAs across a wide range of species, these molecules remain largely unexplored in most eukaryotic microbes, including yeast pathogens of the Candida clade. Here, we analyze thousands of publicly available sequencing datasets to infer and characterize the lncRNA repertoires of five major Candida pathogens: Candida albicans, Candida tropicalis, Candida parapsilosis, Candida auris and Candida glabrata. Our results indicate that genomes of these species encode hundreds of lncRNAs that show levels of evolutionary constraint intermediate between those of intergenic genomic regions and protein-coding genes. Despite their low sequence conservation across the studied species, some lncRNAs are syntenic and are enriched in shared sequence motifs. We find co-expression of lncRNAs with certain protein-coding transcripts, hinting at potential functional associations. Finally, we identify lncRNAs that are differentially expressed during infection of human epithelial cells for four of the studied species. Our comprehensive bioinformatic analyses of Candida lncRNAs pave the way for future functional characterization of these transcripts.

Long non-coding RNAs (lncRNAs) play roles in key cellular processes, but remain largely unexplored in fungal pathogens such as Candida. Here, Hovhannisyan and Gabaldón analyze thousands of sequencing datasets to infer and characterize the lncRNA repertoires of five Candida species, paving the way for their future functional characterization.

Transcriptional Landscapes of Long Non-coding RNAs and Alternative Splicing in Pyricularia oryzae Revealed by RNA-Seq

Pyricularia oryzae causes the rice blast, which is one of the most devastating crop diseases worldwide, and is a model fungal pathogen widely used for dissecting the molecular mechanisms underlying fungal virulence/pathogenicity. Although the whole genome sequence of P. oryzae is publicly available, its current transcriptomes remain incomplete, lacking the information on non-protein coding genes and alternative splicing. Here, we performed and analyzed RNA-Seq of conidia and hyphae, resulting in the identification of 3,374 novel genes. Interestingly, the vast majority of these novel genes likely transcribed long non-coding RNAs (lncRNAs), and most of them were localized in the intergenic regions. Notably, their expressions were concomitant with the transcription of neighboring genes thereof in conidia and hyphae. In addition, 2,358 genes were found to undergo alternative splicing events. Furthermore, we exemplified that a lncRNA was important for hyphal growth likely by regulating the neighboring protein-coding gene and that alternative splicing of the transcription factor gene CON7 was required for appressorium formation. In summary, results from this study indicate that lncRNA transcripts and alternative splicing events are two important mechanisms for regulating the expression of genes important for conidiation, hyphal growth, and pathogenesis, and provide new insights into transcriptomes and gene regulation in the rice blast fungus.

Extracting novel hypotheses and findings from RNA-seq data

Over the past decade, improvements in technology and methods have enabled rapid and relatively inexpensive generation of high-quality RNA-seq datasets. These datasets have been used to characterize gene expression for several yeast species and have provided systems-level insights for basic biology, biotechnology and medicine. Herein, we discuss new techniques that have emerged and existing techniques that enable analysts to extract information from multifactorial yeast RNA-seq datasets. Ultimately, this minireview seeks to inspire readers to query datasets, whether previously published or freshly obtained, with creative and diverse methods to discover and support novel hypotheses.

Long noncoding RNAs are potentially involved in the degeneration of virulence in an aphid-obligate pathogen, Conidiobolus obscurus (Entomophthoromycotina)

Virulence attenuation frequently occurs in in vitro culturing of pathogenic microbes. In this study, we investigated the total putative long noncoding RNAs (lncRNAs) in an aphid-obligate pathogen, Conidiobolus obscurus, and screened the differentially expressed (DE) lncRNAs and protein-coding genes involved in the virulence decline. The virulence was significantly attenuated after eight subculturing events, in which the median lethal concentration of the conidia ejected from mycelial mats relative to the bamboo aphid, Takecallis taiwanus, increased from 36.1 to 126.1 conidia mm^–2, four days after inoculation. In total, 1,252 lncRNAs were identified based on the genome-wide transcriptional analysis. By characterizing their molecular structures and expression patterns, we found that the lncRNAs possessed shorter transcripts, lower expression, and fewer exons than did protein-coding genes in C. obscurus. A total of 410 DE genes of 329 protein-coding genes and 81 lncRNAs were identified. The functional enrichment analysis showed the DE genes were enriched in peptidase activity, protein folding, autophagy, and metabolism. Moreover, target prediction analysis of the 81 lncRNAs revealed 3,111 cis-regulated and 23 trans-regulated mRNAs, while 121 DE lncRNA-mRNA pairs were possibly involved in virulence decline. Moreover, the DE lncRNA-regulated target genes mainly encoded small heat shock proteins, secretory proteins, transporters, autophagy proteins, and other stress response-related proteins. This implies that the decline in virulence regulated by lncRNAs was likely associated with the environmental stress response of C. obscurus. Hence, these findings can provide insights into the lncRNA molecules of Entomophthoromycotina, with regards to virulence regulators of entomopathogens.

Stage-specific regulation of purine metabolism during infectious growth and sexual reproduction in Fusarium graminearum

Ascospores generated during sexual reproduction are the primary inoculum for the wheat scab fungus Fusarium graminearum. Purine metabolism is known to play important roles in fungal pathogens but its lifecycle stage-specific regulation is unclear. By characterizing the genes involved in purine de novo and salvage biosynthesis pathways, we showed that de novo syntheses of inosine, adenosine and guanosine monophosphates (IMP, AMP and GMP) are important for vegetative growth, sexual/asexual reproduction, and infectious growth, whereas purine salvage synthesis is dispensable for these stages in F. graminearum. Addition of GMP rescued the defects of the Fgimd1 mutant in vegetative growth and conidiation but not sexual reproduction, whereas addition of AMP rescued all of these defects of the Fgade12 mutant, suggesting that the function of de novo synthesis of GMP rather than AMP is distinct in sexual stages. Moreover, Acd1, an ortholog of AMP deaminase, is dispensable for growth but essential for ascosporogenesis and pathogenesis, suggesting that AMP catabolism has stage-specific functions during sexual reproduction and infectious growth. The expression of almost all the genes involved in de novo purine synthesis is downregulated during sexual reproduction and infectious growth relative to vegetative growth. This study revealed that F. graminearum has stage-specific regulation of purine metabolism during infectious growth and sexual reproduction.

In Silico Analysis of Common Long Noncoding RNAs in Schistosoma mansoni and Schistosoma haematobium

Background

Schistosomiasis caused by Schistosoma parasites is one of the most common parasitic infections worldwide. Genetic regulation of the genus Schistosoma, which has different developmental stages throughout its life, is quite complex. In these parasites, thousands of long noncoding RNAs (lncRNAs) estimated to be functional were identified. Identifying the transcripts expressed in common and detecting their functions for better understanding of the role of these lncRNAs require a comparative study.

Methods

Assembled RNA-seq datasets belonging to S. mansoni and S. haematobium were obtained from the National Center for Biotechnology. A basic local alignment search tool (BLASTN) analysis was conducted against previously constructed lncRNA library to identify the common lncRNAs between two species. LncRNAs target genes and their gene ontology annotation was performed.

Results

In S. mansoni and S. haematobium, 5132 and 3589 lncRNA transcripts were detected, respectively. These two species had 694 lncRNAs in common. A significant number of lncRNAs was determined to be transcribed from sex chromosomes. The frequently expressed lncRNAs appear to be involved in metabolic and biological regulation processes.

Conclusions

These two species share similar lncRNAs; thus, this finding is a clue that they might have similar functions. In sexual development, they especially might play important roles. Our results will provide important clues to further studies about interactions between human hosts and parasites and the infection mechanisms of Schistosoma parasites.

Upregulation of Long Non-Coding RNA GCC2-AS1 Facilitates Malignant Phenotypes and Correlated With Unfavorable Prognosis for Lung Adenocarcinoma

Background

The role played by long noncoding RNA GCC2-AS1 in primary malignant tumors remains poorly understood. This study aimed to determine the expression levels and evaluate the clinical significance and biological effects of GCC2-AS1 in lung adenocarcinoma (LUAD).

Methods

We used data obtained from tissue samples and the TCGA database to determine the levels of GCC2-AS1 expression LUAD patients, and the prognostic value of those levels. Functional experiments were performed to investigate the effect of GCC2-AS1 on LUAD cells. Genes that were differentially expressed in GCC2-AS1-low and -high groups were analyzed by an enrichment analysis. Additionally, a nomogram model was created and subgroup analyses were performed to further determine the prognostic value of GCC2-AS1 in LUAD.

Results

GCC2-AS1 expression was significantly upregulated in lung adenocarcinoma tissues as compared with normal tissues. Depletion of GCC2-AS1 inhibited the proliferation and invasion of LUAD cells in vitro. An elevated level of GCC2-AS1 was strongly correlated with shorter overall survival time and was identified as an independent prognostic marker for LUAD patients. Enrichment analyses conducted using GO, KEGG, and GSEA databases were performed to identify biological pathways that might involve GCC2-AS1. Several subgroups were found to have a significant prognostic value for patients in the GCC2-AS1-low and -high groups.

Conclusions

Our findings suggest that GCC2-AS1 can serve as a diagnostic and prognostic biomarker for LUAD patients.

A novel antisense long non-coding RNA participates in asexual and sexual reproduction by regulating the expression of GzmetE in Fusarium graminearum

Fusarium graminearum is an important worldwide pathogen that causes Fusarium head blight in wheat, barley, maize and other grains. LncRNAs play important roles in many biological processes, but little is known about their functions and mechanisms in filamentous fungi. Here, we report that a natural antisense RNA, GzmetE-AS, is transcribed from the opposite strand of GzmetE. GzmetE encodes a homoserine O-acetyltransferase, which is important for sexual development and plant infection. The expression of GzmetE-AS was increased significantly during the conidiation stage, while GzmetE was upregulated in the late stage of sexual reproduction. Overexpression of GzmetE-AS inhibited the transcription of GzmetE. In contrast, the expression of GzmetE was significantly increased in GzmetE-AS transcription termination strain GzmetE-AS-T. Furthermore, GzmetE-AS-T produced more perithecia and facilitated the ascospore discharge, resembling the phenotype of GzmetE overexpressing strains. However, overexpression of GzmetE-AS in ∆dcl1/2 strain cannot inhibit the expression of GzmetE, and the GzmetE nat-siRNA is also significantly reduced in ∆dcl1/2 mutant. Taken together, we have identified a novel antisense lncRNA GzmetE-AS, which is involved in asexual and sexual reproduction by regulating its antisense gene GzmetE through RNAi pathway. Our findings reveal that the lncRNA plays critical roles in the development of F. graminearum.

Advances in understanding the evolution of fungal genome architecture

Diversity within the fungal kingdom is evident from the wide range of morphologies fungi display as well as the various ecological roles and industrial purposes they serve. Technological advances, particularly in long-read sequencing, coupled with the increasing efficiency and decreasing costs across sequencing platforms have enabled robust characterization of fungal genomes. These sequencing efforts continue to reveal the rampant diversity in fungi at the genome level. Here, we discuss studies that have furthered our understanding of fungal genetic diversity and genomic evolution. These studies revealed the presence of both small-scale and large-scale genomic changes. In fungi, research has recently focused on many small-scale changes, such as how hypermutation and allelic transmission impact genome evolution as well as how and why a few specific genomic regions are more susceptible to rapid evolution than others. High-throughput sequencing of a diverse set of fungal genomes has also illuminated the frequency, mechanisms, and impacts of large-scale changes, which include chromosome structural variation and changes in chromosome number, such as aneuploidy, polyploidy, and the presence of supernumerary chromosomes. The studies discussed herein have provided great insight into how the architecture of the fungal genome varies within species and across the kingdom and how modern fungi may have evolved from the last common fungal ancestor and might also pave the way for understanding how genomic diversity has evolved in all domains of life.

An upstream open reading frame regulates vasculogenic mimicry of glioma via ZNRD1-AS1/miR-499a-5p/ELF1/EMI1 pathway

Increasing evidence has suggested that gliomas can supply blood through vasculogenic mimicry. In this study, the expression and function of ZNRD1‐AS1‐144aa‐uORF (144aa‐uORF) and some non‐coding RNAs in gliomas were assessed. Real‐time quantitative PCR or Western blot was used to discover the expression of 144aa‐uORF, ZNRD1‐AS1, miR‐499a‐5p, ELF1 and EMI1 in gliomas. In addition, RIP and RNA pull‐down assays were applied to explore the interrelationship between 144aa‐uORF and ZNRD1‐AS1. The role of the 144aa‐uORF\ZNRD1‐AS1\miR‐499a‐5p\ELF1\EMI1 axis in vasculogenic mimicry formation of gliomas was analysed. This study illustrates the reduced expression of the 144aa‐uORF in glioma tissues and cells. Up‐regulation of 144aa‐uORF inhibits proliferation, migration, invasion and vasculogenic mimicry formation within glioma cells. The up‐regulated 144aa‐uORF can increase the degradation of ZNRD1‐AS1 through the nonsense‐mediated RNA decay (NMD) pathway. Knockdown of ZNRD1‐AS1 inhibits vasculogenic mimicry in glioma cells by modulating miR‐499a‐5p. At the same time, miR‐499a‐5p is down‐regulated and has a tumour‐suppressive effect in gliomas. In addition, ZNRD1‐AS1 serves as a competitive endogenous RNA (ceRNA) and regulates the expression of ELF1 by binding to miR‐499a‐5p. Notably, ELF1 binds to the promoter region of EMI1 and up‐regulates EMI1 expression, while simultaneously promoting vasculogenic mimicry in glioma cells. This study suggests that the 144aa‐uORF\ZNRD1‐AS1\miR‐499a‐5p\ELF1\EMI1 axis takes key part in regulating the formation of vasculogenic mimicry in gliomas and may provide a potential target for glioma treatment.

LncRNAs in molluscan and mammalian stages of parasitic schistosomes are developmentally-regulated and coordinately expressed with protein-coding genes

Despite the low level expression of some long noncoding RNAs (lncRNAs), the differential expression of specific lncRNAs plays important roles during the development of many organisms. Schistosomes, parasitic flatworms that are responsible for schistosomiasis, infects over 200 million people resulting in chronic disease and hundreds of thousands of deaths. Schistosomes have a complex life cycle that transitions between molluscan and mammalian hosts. In a molluscan snail host, the sporocyst stage develops over 5 weeks undergoing asexual reproduction to give rise to free-swimming and infectious cercariae that penetrate human skin and eventually mature into egg producing worms in mammals. The tight integration of the sporocyst to the snail host hepatopancreas hinders the -omics study in the molluscan stage, so the sporocyst transcriptome has only been examined for lncRNAs in immature in vitro samples. Here we analyzed the in vivo mature sporocyst transcriptome to identify 4,930 total lncRNAs between the molluscan and mammalian stages of the parasite. We further demonstrate that the lncRNAs are differentially expressed in a development-dependent manner. In addition, we constructed a co-expression correlation network between lncRNAs and protein-coding (PC) genes that was used to identify clusters of lncRNA transcripts with potential functional relevance. We also describe lncRNA-lncRNA and lncRNA-kinome correlations that identify lncRNAs with prospective roles in gene regulation. Finally, our results show clear differential expression patterns of lncRNAs in host-dependent development stages of S. mansoni and ascribe potential functional roles in development based on predicted intracellular interaction.

The novel bZIP transcription factor Fpo1 negatively regulates perithecial development by modulating carbon metabolism in the ascomycete fungus Fusarium graminearum

Fungal sexual reproduction requires complex cellular differentiation processes of hyphal cells. The plant pathogenic fungus Fusarium graminearum produces fruiting bodies called perithecia via sexual reproduction, and perithecia forcibly discharge ascospores into the air for disease initiation and propagation. Lipid metabolism and accumulation are closely related to perithecium formation, yet the molecular mechanisms that regulate these processes are largely unknown. Here, we report that a novel fungal specific bZIP transcription factor, F. graminearum perithecium overproducing 1 (Fpo1), plays a role as a global transcriptional repressor during perithecium production and maturation in F. graminearum. Deletion of FPO1 resulted in reduced vegetative growth, asexual sporulation and virulence and overproduced perithecium, which reached maturity earlier, compared with the wild type. Intriguingly, the hyphae of the fpo1 mutant accumulated excess lipids during perithecium production. Using a combination of molecular biological, transcriptomic and biochemical approaches, we demonstrate that repression of FPO1 after sexual induction leads to reprogramming of carbon metabolism, particularly fatty acid production, which affects sexual reproduction of this fungus. This is the first report of a perithecium-overproducing F. graminearum mutant, and the findings provide comprehensive insight into the role of modulation of carbon metabolism in the sexual reproduction of fungi.

A novel lncRNA as a positive regulator of carotenoid biosynthesis in Fusarium

The fungi Fusarium oxysporum and Fusarium fujikuroi produce carotenoids, lipophilic terpenoid pigments of biotechnological interest, with xanthophyll neurosporaxanthin as the main end product. Their carotenoid biosynthesis is activated by light and negatively regulated by the RING-finger protein CarS. Global transcriptomic analysis identified in both species a putative 1-kb lncRNA that we call carP, referred to as Fo-carP and Ff-carP in each species, upstream to the gene carS and transcribed from the same DNA strand. Fo-carP and Ff-carP are poorly transcribed, but their RNA levels increase in carS mutants. The deletion of Fo-carP or Ff-carP in the respective species results in albino phenotypes, with strong reductions in mRNA levels of structural genes for carotenoid biosynthesis and higher mRNA content of the carS gene, which could explain the low accumulation of carotenoids. Upon alignment, Fo-carP and Ff-carP show 75-80% identity, with short insertions or deletions resulting in a lack of coincident ORFs. Moreover, none of the ORFs found in their sequences have indications of possible coding functions. We conclude that Fo-carP and Ff-carP are regulatory lncRNAs necessary for the active expression of the carotenoid genes in Fusarium through an unknown molecular mechanism, probably related to the control of carS function or expression.

Integrative Activity of Mating Loci, Environmentally Responsive Genes, and Secondary Metabolism Pathways during Sexual Development of Chaetomium globosum

Fungal diversity has amazed evolutionary biologists for decades. One societally important aspect of this diversity manifests in traits that enable pathogenicity. The opportunistic pathogen Chaetomium globosum is well adapted to a high-humidity environment and produces numerous secondary metabolites that defend it from predation. Many of these chemicals can threaten human health. Understanding the phases of the C. globosum life cycle in which these products are made enables better control and even utilization of this fungus. Among its intriguing traits is that it both is self-fertile and lacks any means of propagule-based asexual reproduction. By profiling genome-wide gene expression across the process of sexual reproduction in C. globosum and comparing it to genome-wide gene expression in the model filamentous fungus N. crassa and other closely related fungi, we revealed associations among mating-type genes, sexual developmental genes, sexual incompatibility regulators, environmentally responsive genes, and secondary metabolic pathways.

The origins and maintenance of the rich fungal diversity have been longstanding issues in evolutionary biology. To investigate how differences in expression regulation contribute to divergences in development and ecology among closely related species, transcriptomes were compared between Chaetomium globosum, a homothallic pathogenic fungus thriving in highly humid ecologies, and Neurospora crassa, a heterothallic postfire saprotroph. Gene expression was quantified in perithecia at nine distinct morphological stages during nearly synchronous sexual development. Unlike N. crassa, expression of all mating loci in C. globosum was highly correlated. Key regulators of the initiation of sexual development in response to light stimuli—including orthologs of N. crassasub-1, sub-1-dependent gene NCU00309, and asl-1—showed regulatory dynamics matching between C. globosum and N. crassa. Among 24 secondary metabolism gene clusters in C. globosum, 11—including the cochliodones biosynthesis cluster—exhibited highly coordinated expression across perithecial development. C. globosum exhibited coordinately upregulated expression of histidine kinases in hyperosmotic response pathways—consistent with gene expression responses to high humidity we identified in fellow pathogen Fusarium graminearum. Bayesian networks indicated that gene interactions during sexual development have diverged in concert with the capacities both to reproduce asexually and to live a self-compatible versus self-incompatible life cycle, shifting the hierarchical roles of genes associated with conidiation and heterokaryon incompatibility in N. crassa and C. globosum. This divergence supports an evolutionary history of loss of conidiation due to unfavorable combinations of heterokaryon incompatibility in homothallic species.

XRN1-associated long non-coding RNAs may contribute to fungal virulence and sexual development in entomopathogenic fungus Cordyceps militaris

BACKGROUND

Numerous long non-coding RNAs (lncRNAs) identified and characterized in mammals, plants, and fungi have been found to play critical regulatory roles in biological processes. However, little is known about the role of lncRNAs in insect pathogenic fungi.

RESULTS

By profiling the transcriptomes of sexual and asexual development in the insect-pathogenic fungus Cordyceps militaris, 4140 lncRNAs were identified and found to be dynamically expressed during fungal development. The lncRNAs had shorter transcript lengths and lower numbers of exons compared to protein-coding genes. The expressed target genes (neighboring and cis-regulated) of various expressed lncRNAs were predicted, and these genes showed significant enrichment in energy metabolism and signaling pathways, such as 'Glycolysis/Gluconeogenesis' and "MAPK signaling pathway". To better understand how lncRNAs function in the fungus, xrn1, the final gene of the NMD pathway, which determines the fate of lncRNAs, was disrupted. The Δxrn1 deletion mutant displayed significant (P < 0.05) attenuation of virulence and a lower growth rate in C. militaris. Quantitative RT-PCR results revealed 10 lncRNAs with significantly higher expression, while 8 of these 10 lncRNA target genes (virulence- and sexual development-related) showed significantly lower expression in Δxrn1 compared to in the wild-type, suggesting that lncRNA expression regulates fungal virulence and sexual development by affecting gene expression.

CONCLUSION

These findings suggest that lncRNAs in C. militaris play important roles in the fungal infection progress and fruiting body production, providing a broad repertoire and resource for further studies of lncRNAs. © 2019 Society of Chemical Industry.

Genome-Wide Identification and Functional Prediction of Long Non-coding RNAs Involved in the Heat Stress Response in Metarhizium robertsii

Long non-coding RNAs (lncRNAs) play a significant role in stress responses. To date, only a few studies have reported the role of lncRNAs in insect-pathogenic fungi. Here, we report a genome-wide transcriptional analysis of lncRNAs produced in response to heat stress in Metarhizium robertsii, a model insect-pathogenic fungus, using strand-specific RNA sequencing. A total of 1655 lncRNAs with 1742 isoforms were identified, of which 1081 differentially expressed (DE) lncRNAs were characterized as being heat responsive. By characterizing their genomic structures and expression patterns, we found that the lncRNAs possessed shorter transcripts, fewer exons, and lower expression levels than the protein-coding genes in M. robertsii. Furthermore, target prediction analysis of the lncRNAs revealed thousands of potential DE lncRNA–messenger RNA (mRNA) pairs, among which 5381 pairs function in the cis-regulatory mode. Further pathway enrichment analysis of the corresponding cis-regulated target genes showed that the targets were significantly enriched in the following biological pathways: the Hippo signaling pathway and cell cycle. This finding suggested that these DE lncRNAs control the expression of their corresponding neighboring genes primarily through environmental information processing and cellular processes. Moreover, only 26 trans-regulated lncRNA–mRNA pairs were determined. In addition, among the targets of heat-responsive lncRNAs, two classic genes that may be involved in the response to heat stress were also identified, including hsp70 (XM_007821830 and XM_007825705). These findings expand our knowledge of lncRNAs as important regulators of the response to heat stress in filamentous fungi, including M. robertsii.

Comparative Genomics and Transcriptomics During Sexual Development Gives Insight Into the Life History of the Cosmopolitan Fungus Fusarium neocosmosporiellum

Fusarium neocosmosporiellum (formerly Neocosmospora vasinfecta) is a cosmopolitan fungus that has been reported from soil, herbivore dung, and as a fruit- and root-rot pathogen of numerous field crops, although it is not known to cause significant losses on any crop. Taking advantage of the fact that this species produces prolific numbers of perithecia in culture, the genome of F. neocosmosporiellum was sequenced and transcriptomic analysis across five stages of perithecium development was performed to better understand the metabolic potential for sexual development and gain insight into its life history. Perithecium morphology together with the genome and transcriptome were compared with those of the plant pathogen F. graminearum, a model for studying perithecium development. Larger ascospores of F. neocosmosporiellum and their tendency to discharge as a cluster demonstrated a duality of dispersal: the majority are passively dispersed through the formation of cirrhi, while a minority of spores are shot longer distances than those of F. graminearum. The predicted gene number in the F. neocosmosporiellum genome was similar to that in F. graminearum, but F. neocosmosporiellum had more carbohydrate metabolism-related and transmembrane transport genes. Many transporter genes were differentially expressed during perithecium development in F. neocosmosporiellum, which may account for its larger perithecia. Comparative analysis of the secondary metabolite gene clusters identified several polyketide synthase genes that were induced during later stages of perithecium development. Deletion of a polyketide synthase gene in F. neocosmosporiellum resulted in a defective perithecium phenotype, suggesting an important role of the corresponding metabolite, which has yet to be identified, in perithecium development. Results of this study have provided novel insights into the genomic underpinning of development in F. neocosmosporiellum, which may help elucidate its ability to occupy diverse ecological niches.

Distribution, Characteristics, and Regulatory Potential of Long Noncoding RNAs in Brown-Rot Fungi

Long noncoding RNAs have been thoroughly studied in plants, animals, and yeasts, where they play important roles as regulators of transcription. Nevertheless, almost nothing is known about their presence and characteristics in filamentous fungi, especially in basidiomycetes. In the present study, we have carried out an exhaustive annotation and characterization of lncRNAs in two lignin degrader basidiomycetes, Coniophora puteana and Serpula lacrymans. We identified 2,712 putative lncRNAs in the former and 2,242 in the latter, mainly originating from intergenic locations of transposon-sparse genomic regions. The lncRNA length, GC content, expression levels, and stability of the secondary structure differ from coding transcripts but are similar in these two species and resemble that of other eukaryotes. Nevertheless, they lack sequence conservation. Also, we found that lncRNAs are transcriptionally regulated in the same proportion as genes when the fungus actively decomposes soil organic matter. Finally, up to 7% of the upstream gene regions of Coniophora puteana and Serpula lacrymans are transcribed and produce lncRNAs. The study of expression trends in these gene-lncRNA pairs uncovered groups with similar and opposite transcriptional profiles which may be the result of cis-transcriptional regulation.