Intron retention-dependent gene regulation in Cryptococcus neoformans

The spliceosome impacts morphogenesis in the human fungal pathogen Candida albicans

At human body temperature, the fungal pathogen Candida albicans can transition from yeast to filamentous morphologies in response to host-relevant cues. Additionally, elevated temperatures encountered during febrile episodes can independently induce C. albicans filamentation. However, the underlying genetic pathways governing this developmental transition in response to elevated temperatures remain largely unexplored. Here, we conducted a functional genomic screen to unravel the genetic mechanisms orchestrating C. albicans filamentation specifically in response to elevated temperature, implicating 45% of genes associated with the spliceosome or pre-mRNA splicing in this process. Employing RNA-Seq to elucidate the relationship between mRNA splicing and filamentation, we identified greater levels of intron retention in filaments compared to yeast, which correlated with reduced expression of the affected genes. Intriguingly, homozygous deletion of a gene encoding a spliceosome component important for filamentation (PRP19) caused even greater levels of intron retention compared with wild type and displayed globally dysregulated gene expression. This suggests that intron retention is a mechanism for fine-tuning gene expression during filamentation, with perturbations of the spliceosome exacerbating this process and blocking filamentation. Overall, this study unveils a novel biological process governing C. albicans filamentation, providing new insights into the complex regulation of this key virulence trait.IMPORTANCEFungal pathogens such as Candida albicans can cause serious infections with high mortality rates in immunocompromised individuals. When C. albicans is grown at temperatures encountered during human febrile episodes, yeast cells undergo a transition to filamentous cells, and this process is key to its virulence. Here, we expanded our understanding of how C. albicans undergoes filamentation in response to elevated temperature and identified many genes involved in mRNA splicing that positively regulate filamentation. Through transcriptome analyses, we found that intron retention is a mechanism for fine-tuning gene expression in filaments, and perturbation of the spliceosome exacerbates intron retention and alters gene expression substantially, causing a block in filamentation. This work adds to the growing body of knowledge on the role of introns in fungi and provides new insights into the cellular processes that regulate a key virulence trait in C. albicans.

Alternative TSS use is widespread in Cryptococcus fungi in response to environmental cues and regulated genome-wide by the transcription factor Tur1

Alternative transcription start site (TSS) usage regulation has been identified as a major means of gene expression regulation in metazoans. However, in fungi, its impact remains elusive as its study has thus far been restricted to model yeasts. Here, we first re-analyzed TSS-seq data to define genuine TSS clusters in 2 species of pathogenic Cryptococcus. We identified 2 types of TSS clusters associated with specific DNA sequence motifs. Our analysis also revealed that alternative TSS usage regulation in response to environmental cues is widespread in Cryptococcus, altering gene expression and protein targeting. Importantly, we performed a forward genetic screen to identify a unique transcription factor (TF) named Tur1, which regulates alternative TSS (altTSS) usage genome-wide when cells switch from exponential phase to stationary phase. ChiP-Seq and DamID-Seq analyses suggest that at some loci, the role of Tur1 might be direct. Tur1 has been previously shown to be essential for virulence in C. neoformans. We demonstrated here that a tur1Δ mutant strain is more sensitive to superoxide stress and phagocytosed more efficiently by macrophages than the wild-type (WT) strain.

Cryptococcus neoformans Slu7 ensures nuclear positioning during mitotic progression through RNA splicing

The position of the nucleus before it divides during mitosis is variable in different budding yeasts. Studies in the pathogenic intron-rich fungus Cryptococcus neoformans reveal that the nucleus moves entirely into the daughter bud before its division. Here, we report functions of a zinc finger motif containing spliceosome protein C. neoformans Slu7 (CnSlu7) in cell cycle progression. The budding yeast and fission yeast homologs of Slu7 have predominant roles for intron 3' splice site definition during pre-mRNA splicing. Using a conditional knockdown strategy, we show CnSlu7 is an essential factor for viability and is required for efficient cell cycle progression with major role during mitosis. Aberrant nuclear migration, including improper positioning of the nucleus as well as the spindle, were frequently observed in cells depleted of CnSlu7. However, cell cycle delays observed due to Slu7 depletion did not activate the Mad2-dependent spindle assembly checkpoint (SAC). Mining of the global transcriptome changes in the Slu7 knockdown strain identified downregulation of transcripts encoding several cell cycle regulators and cytoskeletal factors for nuclear migration, and the splicing of specific introns of these genes was CnSlu7 dependent. To test the importance of splicing activity of CnSlu7 on nuclear migration, we complemented Slu7 knockdown cells with an intron less PAC1 minigene and demonstrated that the nuclear migration defects were significantly rescued. These findings show that CnSlu7 regulates the functions of diverse cell cycle regulators and cytoskeletal components, ensuring timely cell cycle transitions and nuclear division during mitosis.

Lessons on fruiting body morphogenesis from genomes and transcriptomes of Agaricomycetes

Fruiting bodies (sporocarps, sporophores or basidiomata) of mushroom-forming fungi (Agaricomycetes) are among the most complex structures produced by fungi. Unlike vegetative hyphae, fruiting bodies grow determinately and follow a genetically encoded developmental program that orchestrates their growth, tissue differentiation and sexual sporulation. In spite of more than a century of research, our understanding of the molecular details of fruiting body morphogenesis is still limited and a general synthesis on the genetics of this complex process is lacking. In this paper, we aim at a comprehensive identification of conserved genes related to fruiting body morphogenesis and distil novel functional hypotheses for functionally poorly characterised ones. As a result of this analysis, we report 921 conserved developmentally expressed gene families, only a few dozens of which have previously been reported to be involved in fruiting body development. Based on literature data, conserved expression patterns and functional annotations, we provide hypotheses on the potential role of these gene families in fruiting body development, yielding the most complete description of molecular processes in fruiting body morphogenesis to date. We discuss genes related to the initiation of fruiting, differentiation, growth, cell surface and cell wall, defence, transcriptional regulation as well as signal transduction. Based on these data we derive a general model of fruiting body development, which includes an early, proliferative phase that is mostly concerned with laying out the mushroom body plan (via cell division and differentiation), and a second phase of growth via cell expansion as well as meiotic events and sporulation. Altogether, our discussions cover 1 480 genes of Coprinopsis cinerea, and their orthologs in Agaricus bisporus, Cyclocybe aegerita, Armillaria ostoyae, Auriculariopsis ampla, Laccaria bicolor, Lentinula edodes, Lentinus tigrinus, Mycena kentingensis, Phanerochaete chrysosporium, Pleurotus ostreatus, and Schizophyllum commune, providing functional hypotheses for ~10 % of genes in the genomes of these species. Although experimental evidence for the role of these genes will need to be established in the future, our data provide a roadmap for guiding functional analyses of fruiting related genes in the Agaricomycetes. We anticipate that the gene compendium presented here, combined with developments in functional genomics approaches will contribute to uncovering the genetic bases of one of the most spectacular multicellular developmental processes in fungi. Citation: Nagy LG, Vonk PJ, Künzler M, Földi C, Virágh M, Ohm RA, Hennicke F, Bálint B, Csernetics Á, Hegedüs B, Hou Z, Liu XB, Nan S, M. Pareek M, Sahu N, Szathmári B, Varga T, Wu W, Yang X, Merényi Z (2023). Lessons on fruiting body morphogenesis from genomes and transcriptomes of Agaricomycetes. Studies in Mycology 104: 1-85. doi: 10.3114/sim.2022.104.01.

Blue Light-Dependent Pre-mRNA Splicing Controls Pigment Biosynthesis in the Mushroom Terana caerulea

Light induces the production of ink-blue pentacyclic natural products, the corticin pigments, in the cobalt crust mushroom Terana caerulea. Here, we describe the genetic locus for corticin biosynthesis and provide evidence for a light-dependent dual transcriptional/cotranscriptional regulatory mechanism. Light selectively induces the expression of the corA gene encoding the gateway enzyme, the first described mushroom polyporic acid synthetase CorA, while other biosynthetic genes for modifying enzymes necessary to complete corticin assembly are induced only at lower levels. The strongest corA induction was observed following exposure to blue and UV light. A second layer of regulation is provided by the light-dependent splicing of the three introns in the pre-mRNA of corA. Our results provide insight into the fundamental organization of how mushrooms regulate natural product biosynthesis. IMPORTANCE The regulation of natural product biosyntheses in mushrooms in response to environmental cues is poorly understood. We addressed this knowledge gap and chose the cobalt crust mushroom Terana caerulea as our model. Our work discovered a dual-level regulatory mechanism that connects light as an abiotic stimulus with a physiological response, i.e., the production of dark-blue pigments. Exposure to blue light elicits strongly increased transcription of the gene encoding the gateway enzyme, the polyporic acid synthetase CorA, that catalyzes the formation of the pigment core structure. Additionally, light is a prerequisite for the full splicing of corA pre-mRNA and, thus, its proper maturation. Dual transcriptional/cotranscriptional light-dependent control of fungal natural product biosynthesis has previously been unknown. As it allows the tight control of a key metabolic step, it may be a much more prevalent mechanism among these organisms.

Organization and Unconventional Integration of the Mating-Type Loci in Morchella Species

True morels (Morchella spp.) are a group of delicious fungi in high demand worldwide, and some species of morels have been successfully cultivated in recent years. To better understand the sexual reproductive mechanisms of these fungi, we characterized the structure of the mating-type loci from ten morel species, and seven of them were obtained using long-range PCR amplification. Among the studied species, eight were heterothallic, two were homothallic, and four types of composition were observed in the MAT loci. In three of the five black morel species, the MAT1-1-1, MAT1-1-10, and MAT1-1-11 genes were in the MAT1-1 idiomorph, and only the MAT1-2-1 gene was in the MAT1-2 idiomorph, while an integration event occurred in the other two species and resulted in the importation of the MAT1-1-11 gene into the MAT1-2 idiomorph and survival as a truncated fragment in the MAT1-1 idiomorph. However, the MAT1-1-11 gene was not available in the four yellow morels and one blushing morel species. M. rufobrunnea, a representative species of the earliest diverging branch of true morels, along with another yellow morel Mes-15, were confirmed to be homothallic, and the MAT1-1-1, MAT1-1-10, and MAT1-2-1 genes were arranged in a tandem array. Therefore, we hypothesized that homothallism should be the ancestral reproductive state in Morchella. RT-PCR analyses revealed that four mating genes could be constitutively expressed, while the MAT1-1-10 gene underwent alternative splicing to produce different splice variants.

Alternative splicing as a source of phenotypic diversity

A major goal of evolutionary genetics is to understand the genetic processes that give rise to phenotypic diversity in multicellular organisms. Alternative splicing generates multiple transcripts from a single gene, enriching the diversity of proteins and phenotypic traits. It is well established that alternative splicing contributes to key innovations over long evolutionary timescales, such as brain development in bilaterians. However, recent developments in long-read sequencing and the generation of high-quality genome assemblies for diverse organisms has facilitated comparisons of splicing profiles between closely related species, providing insights into how alternative splicing evolves over shorter timescales. Although most splicing variants are probably non-functional, alternative splicing is nonetheless emerging as a dynamic, evolutionarily labile process that can facilitate adaptation and contribute to species divergence.

Landscape and regulation of alternative splicing and alternative polyadenylation in a plant pathogenic fungus

Alternative splicing (AS) and alternative polyadenylation (APA) contribute significantly to the regulation of gene expression in higher eukaryotes. Their biological impact in filamentous fungi, however, is largely unknown. Here we combine PacBio Isoform-Sequencing and strand-specific RNA-sequencing of multiple tissues and mutant characterization to reveal the landscape and regulation of AS and APA in Fusarium graminearum. We generated a transcript annotation comprising 51 617 isoforms from 17 189 genes. In total, 4997 and 11 133 genes are alternatively spliced and polyadenylated, respectively. Majority of the AS events alter coding sequences. Unexpectedly, the AS transcripts containing premature-termination codons are not sensitive to nonsense-mediated messenger RNA decay. Unlike in yeasts and animals, distal APA sites have strong signals, but proximal APA isoforms are highly expressed in F. graminearum. The 3'-end processing factors FgRNA15, FgHRP1, and FgFIP1 play roles in promoting proximal APA site usage and intron splicing. A genome-wide increase in intron inclusion and distal APA site usage and downregulation of the spliceosomal and 3'-end processing factors were observed in older and quiescent tissues, indicating intron inclusion and 3'-untranslated region lengthening as novel mechanisms in regulating aging and dormancy in fungi. This study provides new insights into the complexity and regulation of AS and APA in filamentous fungi.

Native RNA sequencing in fission yeast reveals frequent alternative splicing isoforms

The unicellular yeast Schizosaccharomyces pombe (fission yeast) retains many of the splicing features observed in humans and is thus an excellent model to study the basic mechanisms of splicing. Nearly half the genes contain introns, but the impact of alternative splicing in gene regulation and proteome diversification remains largely unexplored. Here we leverage Oxford Nanopore Technologies native RNA sequencing (dRNA), as well as ribosome profiling data, to uncover the full range of polyadenylated transcripts and translated open reading frames. We identify 332 alternative isoforms affecting the coding sequences of 262 different genes, 97 of which occur at frequencies >20%, indicating that functional alternative splicing in S. pombe is more prevalent than previously suspected. Intron retention events make ∼80% of the cases; these events may be involved in the regulation of gene expression and, in some cases, generate novel protein isoforms, as supported by ribosome profiling data in 18 of the intron retention isoforms. One example is the rpl22 gene, in which intron retention is associated with the translation of a protein of only 13 amino acids. We also find that lowly expressed transcripts tend to have longer poly(A) tails than highly expressed transcripts, highlighting an interdependence between poly(A) tail length and transcript expression level. Finally, we discover 214 novel transcripts that are not annotated, including 158 antisense transcripts, some of which also show translation evidence. The methodologies described in this work open new opportunities to study the regulation of splicing in a simple eukaryotic model.

Intron retention coupled with nonsense-mediated decay is involved in cellulase biosynthesis in cellulolytic fungi

BACKGROUND

Knowledge on regulatory networks associated with cellulase biosynthesis is prerequisite for exploitation of such regulatory systems in enhancing cellulase production with low cost. The biological functions of intron retention (IR) and nonsense-mediated mRNA decay (NMD) in filamentous fungi is lack of study, let alone their roles in cellulase biosynthesis.

RESULTS

We found that major cellulase genes (cel7a, cel7b, and cel3a) exhibited concomitant decrease in IR rates and increase in their gene expression in T. reesei under cellulase-producing condition (cellulose and lactose) that was accompanied with a more active NMD pathway, as compared to cellulase non-producing condition (glucose). In the presence of the NMD pathway inhibitor that successfully repressed the NMD pathway, the mRNA levels of cellulase genes were sharply down-regulated, but the rates of IR in these genes were significantly up-regulated. Consistently, the cellulase activities were severely inhibited. In addition, the NMD pathway inhibitor caused the downregulated mRNA levels of two important genes of the target of rapamycin (TOR) pathway, trfkbp12 and trTOR1. The absence of gene trfkbp12 made the cellulase production in T. reesei more sensitive to the NMD pathway inhibitor.

CONCLUSIONS

All these findings suggest that the IR of cellulase genes regulates their own gene expression by coupling with the NMD pathway, which might involve the TOR pathway. Our results provide better understanding on intron retention, the NMD pathway, and cellulase production mechanism in filamentous fungi.

Alternative splicing diversifies the transcriptome and proteome of the rice blast fungus during host infection

Alternative splicing (AS) contributes to diversifying and regulating cellular responses to environmental conditions and developmental cues by differentially producing multiple mRNA and protein isoforms from a single gene. Previous studies on AS in pathogenic fungi focused on profiling AS isoforms under a limited number of conditions. We analysed AS profiles in the rice blast fungus Magnaporthe oryzae, a global threat to rice production, using high-quality transcriptome data representing its vegetative growth (mycelia) and multiple host infection stages. We identified 4,270 AS isoforms derived from 2,413 genes, including 499 genes presumably regulated by infection-specific AS. AS appears to increase during infection, with 32.7% of the AS isoforms being produced during infection but absent in mycelia. Analysis of the isoforms observed at each infection stage showed that 636 AS isoforms were more abundant than corresponding annotated mRNAs, especially after initial hyphal penetration into host cell. Many such dominant isoforms were predicted to encode regulatory proteins such as transcription factors and phospho-transferases. We also identified the genes encoding distinct proteins via AS and confirmed the translation of some isoforms via a proteomic analysis, suggesting potential AS-mediated neo-functionalization of some genes during infection. Comprehensive profiling of the pattern of genome-wide AS during multiple stages of rice-M. oryzae interaction established a foundational resource that will help investigate the role and regulation of AS during rice infection.

Intron distribution and emerging role of alternative splicing in fungi

Spliceosomal introns are noncoding sequences that are spliced from pre-mRNA. They are ubiquitous in eukaryotic genomes, although the average number of introns per gene varies considerably between different eukaryotic species. Fungi are diverse in terms of intron numbers ranging from 4% to 99% genes with introns. Alternative splicing is one of the most common modes of posttranscriptional regulation in eukaryotes, giving rise to multiple transcripts from a single pre-mRNA and is widespread in metazoans and drives extensive proteome diversity. Earlier, alternative splicing was considered to be rare in fungi, but recently, increasing numbers of studies have revealed that alternative splicing is also widespread in fungi and has been implicated in the regulation of fungal growth and development, protein localization and the improvement of survivability, likely underlying their unique capacity to adapt to changing environmental conditions. However, the role of alternative splicing in pathogenicity and development of drug resistance is only recently gaining attention. In this review, we describe the intronic landscape in fungi. We also present in detail the newly discovered functions of alternative splicing in various cellular processes and outline areas particularly in pathogenesis and clinical drug resistance for future studies that could lead to the development of much needed new therapeutics.

Post-transcriptional and translational control of the morphology and virulence in human fungal pathogens

Host-pathogen interactions at the molecular level are the key to fungal pathogenesis. Fungal pathogens utilize several mechanisms such as adhesion, invasion, phenotype switching and metabolic adaptations, to survive in the host environment and respond. Post-transcriptional and translational regulations have emerged as key regulatory mechanisms ensuring the virulence and survival of fungal pathogens. Through these regulations, fungal pathogens effectively alter their protein pool, respond to various stress, and undergo morphogenesis, leading to efficient and comprehensive changes in fungal physiology. The regulation of virulence through post-transcriptional and translational regulatory mechanisms is mediated through mRNA elements (cis factors) or effector molecules (trans factors). The untranslated regions upstream and downstream of the mRNA, as well as various RNA-binding proteins involved in translation initiation or circularization of the mRNA, play pivotal roles in the regulation of morphology and virulence by influencing protein synthesis, protein isoforms, and mRNA stability. Therefore, post-transcriptional and translational mechanisms regulating the morphology, virulence and drug-resistance processes in fungal pathogens can be the target for new therapeutics. With improved "omics" technologies, these regulatory mechanisms are increasingly coming to the forefront of basic biology and drug discovery. This review aims to discuss various modes of post-transcriptional and translation regulations, and how these mechanisms exert influence in the virulence and morphogenesis of fungal pathogens.

Cryptococcus extracellular vesicles properties and their use as vaccine platforms

Whereas extracellular vesicle (EV) research has become commonplace in different biomedical fields, this field of research is still in its infancy in mycology. Here we provide a robust set of data regarding the structural and compositional aspects of EVs isolated from the fungal pathogenic species Cryptococcus neoformans, C. deneoformans and C. deuterogattii. Using cutting-edge methodological approaches including cryogenic electron microscopy and cryogenic electron tomography, proteomics, and flow cytometry, we revisited cryptococcal EV features and suggest a new EV structural model, in which the vesicular lipid bilayer is covered by mannoprotein-based fibrillar decoration, bearing the capsule polysaccharide as its outer layer. About 10% of the EV population is devoid of fibrillar decoration, adding another aspect to EV diversity. By analysing EV protein cargo from the three species, we characterized the typical Cryptococcus EV proteome. It contains several membrane-bound protein families, including some Tsh proteins bearing a SUR7/PalI motif. The presence of known protective antigens on the surface of Cryptococcus EVs, resembling the morphology of encapsulated virus structures, suggested their potential as a vaccine. Indeed, mice immunized with EVs obtained from an acapsular C. neoformans mutant strain rendered a strong antibody response in mice and significantly prolonged their survival upon C. neoformans infection.

Application of an optimized annotation pipeline to the Cryptococcus deuterogattii genome reveals dynamic primary metabolic gene clusters and genomic impact of RNAi loss

Evaluating the quality of a de novo annotation of a complex fungal genome based on RNA-seq data remains a challenge. In this study, we sequentially optimized a Cufflinks-CodingQuary-based bioinformatics pipeline fed with RNA-seq data using the manually annotated model pathogenic yeasts Cryptococcus neoformans and Cryptococcus deneoformans as test cases. Our results show that the quality of the annotation is sensitive to the quantity of RNA-seq data used and that the best quality is obtained with 5–10 million reads per RNA-seq replicate. We also showed that the number of introns predicted is an excellent a priori indicator of the quality of the final de novo annotation. We then used this pipeline to annotate the genome of the RNAi-deficient species Cryptococcus deuterogattii strain R265 using RNA-seq data. Dynamic transcriptome analysis revealed that intron retention is more prominent in C. deuterogattii than in the other RNAi-proficient species C. neoformans and C. deneoformans. In contrast, we observed that antisense transcription was not higher in C. deuterogattii than in the two other Cryptococcus species. Comparative gene content analysis identified 21 clusters enriched in transcription factors and transporters that have been lost. Interestingly, analysis of the subtelomeric regions in these three annotated species identified a similar gene enrichment, reminiscent of the structure of primary metabolic clusters. Our data suggest that there is active exchange between subtelomeric regions, and that other chromosomal regions might participate in adaptive diversification of Cryptococcus metabolite assimilation potential.

Population genomic analysis of Cryptococcus Brazilian isolates reveals an African type subclade distribution

The genomes of a large number of Cryptococcus neoformans isolates have been sequenced and analyzed in recent years. These genomes have been used to understand the global population structure of this opportunistic pathogen. However, only a small number of South American isolates have been considered in these studies, and the population structure of C. neoformans in this part of the world remains elusive. Here, we analyzed the genomic sequences of 53 Brazilian Cryptococcus isolates and deciphered the C. neoformans population structure in this country. Our data reveal an African-like structure that suggested repeated intercontinental transports from Africa to South America. We also identified a mutator phenotype in one VNBII Brazilian isolate, exemplifying how fast-evolving isolates can shape the Cryptococcus population structure. Finally, phenotypic analyses revealed wide diversity but not lineage specificity in the expression of classical virulence traits within the set of isolates.

Alternative splicing regulates the α-glucosidase synthesis in Aspergillus neoniger NCIM 1400

Aspergillus neoniger NCIM 1400 whose cell-free fraction was earlier established for transglycosylation activity conferred by α-glucosidase gene (agdA), was subjected to sequence analysis. Preliminary results revealed certain dynamics in the intron splicing mechanism, and to ascertain these molecular events, a detailed study was carried. The electrophoresis results from the cDNA portion (B-fragment) of agdA showed multiple bands, indicating the amplification of one or more fragments. The sequence results of cDNA cloned vector revealed the retention type of alternative splicing in the agdA. The splicing mechanism of agdA in NCIM 1400 was compared to different A. niger strains, which harbours agdA orthologues, using PCR. It was observed that effective intron splicing leads to higher α-glucosidase activity from these selected Aspergillus spp. To explore the dynamics of intron retention in A. neoniger NCIM 1400, time-course analysis of intron retention, enzyme activity, and sugar consumption were carried over a period of 168 h of fungal growth. RT-qPCR results revealed that introns retention was not detected during the initial growth phase when the maltose and its hydrolysed product, glucose were consumed. Here we demonstrate that exhaustion of maltose causes increase in retention of introns in the mRNA transcripts of agdA gene, and this could be the possible mode of regulating this gene.

Genome-wide profiling of alternative splicing genes in hybrid poplar (P.alba×P.glandulosa cv.84K) leaves

Alternative splicing (AS) is a post-transcriptional process common in plants and essential for regulation of environmental fitness of plants. In the present study, we focus on the AS events in poplar leaves to understand their effects on plant growth and development. The hybrid poplar (P.alba×P.glandulosa cv.84K) leaves were collected for RNA extraction. The extracted RNA was sequenced using on an Illumina HiSeq™ 2000 platform. Using the Populus trichocarpa genome as the reference, a total of 3810 AS genes were identified (9225 AS events), which accounted for 13.51% of all the expressed genes. Intron retention was the most common AS event, accounting for 43.86% of all the AS events, followed by alternative 3′ splice sites (23.75%), alternative 5′ splice sites (23.71%), and exon skipping (8.68%). Chromosomes 10 had the most condensed AS events (33.67 events/Mb) and chromosome 19 had the least (12.42 events/Mb). Association analysis showed that AS in the poplar leaves was positively correlated with intron length, exon number, exon length, and gene expression level, and was negatively correlated with GC content. AS genes in the poplar leaves were associated mainly with inositol phosphate metabolism and phosphatidylinositol signaling system pathways that would be significant on wooden plant production.

Identification of Genomewide Alternative Splicing Events in Sequential, Isogenic Clinical Isolates of Candida albicans Reveals a Novel Mechanism of Drug Resistance and Tolerance to Cellular Stresses

The emergence of resistance in Candida albicans, an opportunistic pathogen, against the commonly used antifungals is becoming a major obstacle in its treatment. The necessity to identify new drug targets demands fundamental insights into the mechanisms used by this organism to develop drug resistance. C. albicans has introns in 4 to 6% of its genes, the functions of which remain largely unknown. Using the RNA-sequencing data from isogenic pairs of azole-sensitive and -resistant isolates of C. albicans, here, we show how C. albicans uses modulations in mRNA splicing to overcome antifungal drug stress.

Alternative splicing (AS)—a process by which a single gene gives rise to different protein isoforms in eukaryotes—has been implicated in many basic cellular processes, but little is known about its role in drug resistance and fungal pathogenesis. The most common human fungal pathogen, Candida albicans, has introns in 4 to 6% of its genes, the functions of which remain largely unknown. Here, we report AS regulating drug resistance in C. albicans. Comparative RNA-sequencing of two different sets of sequential, isogenic azole-sensitive and -resistant isolates of C. albicans revealed differential expression of splice isoforms of 14 genes. One of these was the superoxide dismutase gene SOD3, which contains a single intron. The sod3Δ/Δ mutant was susceptible to the antifungals amphotericin B (AMB) and menadione (MND). While AMB susceptibility was rescued by overexpression of both the spliced and unspliced SOD3 isoforms, only the spliced isoform could overcome MND susceptibility, demonstrating the functional relevance of this splicing in developing drug resistance. Furthermore, unlike AMB, MND inhibits SOD3 splicing and acts as a splicing inhibitor. Consistent with these observations, MND exposure resulted in increased levels of unspliced SOD3 isoform that are unable to scavenge reactive oxygen species (ROS), resulting in increased drug susceptibility. Collectively, these observations suggest that AS is a novel mechanism for stress adaptation and overcoming drug susceptibility in C. albicans.

IMPORTANCE The emergence of resistance in Candida albicans, an opportunistic pathogen, against the commonly used antifungals is becoming a major obstacle in its treatment. The necessity to identify new drug targets demands fundamental insights into the mechanisms used by this organism to develop drug resistance. C. albicans has introns in 4 to 6% of its genes, the functions of which remain largely unknown. Using the RNA-sequencing data from isogenic pairs of azole-sensitive and -resistant isolates of C. albicans, here, we show how C. albicans uses modulations in mRNA splicing to overcome antifungal drug stress.

Quantitative global studies reveal differential translational control by start codon context across the fungal kingdom

Eukaryotic protein synthesis generally initiates at a start codon defined by an AUG and its surrounding Kozak sequence context, but the quantitative importance of this context in different species is unclear. We tested this concept in two pathogenic Cryptococcus yeast species by genome-wide mapping of translation and of mRNA 5' and 3' ends. We observed thousands of AUG-initiated upstream open reading frames (uORFs) that are a major contributor to translation repression. uORF use depends on the Kozak sequence context of its start codon, and uORFs with strong contexts promote nonsense-mediated mRNA decay. Transcript leaders in Cryptococcus and other fungi are substantially longer and more AUG-dense than in Saccharomyces. Numerous Cryptococcus mRNAs encode predicted dual-localized proteins, including many aminoacyl-tRNA synthetases, in which a leaky AUG start codon is followed by a strong Kozak context in-frame AUG, separated by mitochondrial-targeting sequence. Analysis of other fungal species shows that such dual-localization is also predicted to be common in the ascomycete mould, Neurospora crassa. Kozak-controlled regulation is correlated with insertions in translational initiation factors in fidelity-determining regions that contact the initiator tRNA. Thus, start codon context is a signal that quantitatively programs both the expression and the structures of proteins in diverse fungi.

Improved Homologous Expression of the Acidic Lipase from Aspergillus niger

In this study, the acidic lipase from Aspergillus niger (ANL) was homologously expressed in A. niger. The expression of ANL was significantly improved by the expression of the native ANL with the introns, the addition of the Kozak sequence and the optimization of the signal sequences. When the cDNA sequence of ANL fused with the glaA signal was expressed under the gpdA promoter in A. niger, no lipase activity could be detected. We then tried to improve the expression by using the full-length ANL gene containing three introns, and the lipase activity in the supernatant reached 75.80 U/ml, probably as a result of a more stable mRNA structure. The expression was further improved to 100.60 U/ml by introducing a Kozak sequence around the start codon due to a higher translation efficiency. Finally, the effects of three signal sequences including the cbhI signal, the ANL signal and the glaA signal on the lipase expression were evaluated. The transformant with the cbhI signal showed the highest lipase activity (314.67 U/ml), which was 1.90-fold and 3.13-fold higher than those with the ANL signal and the glaA signal, respectively. The acidic lipase was characterized and its highest activity was detected at pH 3.0 and a temperature of 45°C. These results provided promising strategies for the production of the acidic lipase from A. niger.

Spliceosomal Prp8 intein at the crossroads of protein and RNA splicing

The spliceosome is a large ribonucleoprotein complex that removes introns from pre-mRNAs. At its functional core lies the essential pre-mRNA processing factor 8 (Prp8) protein. Across diverse eukaryotes, this protein cofactor of RNA catalysis harbors a self-splicing element called an intein. Inteins in Prp8 are extremely pervasive and are found at 7 different sites in various species. Here, we focus on the Prp8 intein from Cryptococcus neoformans (Cne), a human fungal pathogen. We solved the crystal structure of this intein, revealing structural homology among protein splicing sequences in eukaryotes, including the Hedgehog C terminus. Working with the Cne Prp8 intein in a reporter assay, we find that the biologically relevant divalent metals copper and zinc inhibit intein splicing, albeit by 2 different mechanisms. Copper likely stimulates reversible modifications on a catalytically important cysteine, whereas zinc binds at the terminal asparagine and the same critical cysteine. Importantly, we also show that copper treatment inhibits Prp8 protein splicing in Cne. Lastly, an intein-containing Prp8 precursor model is presented, suggesting that metal-induced protein splicing inhibition would disturb function of both Prp8 and the spliceosome. These results indicate that Prp8 protein splicing can be modulated, with potential functional implications for the spliceosome.

Studying fungal pathogens of humans and fungal infections: fungal diversity and diversity of approaches

Seminal work by Louis Pasteur revealed the contribution of fungi - yeasts and microsporidia to agroindustry and disease in animals, respectively. More than 150 years later, the impact of fungi on human health and beyond is an ever-increasing issue, although often underestimated. Recent studies estimate that fungal infections, especially those caused by Candida, Cryptococcus and Aspergillus species, kill more than one million people annually. Indeed, these neglected infections are in general very difficult to cure and the associated mortality remains very high even when antifungal treatments exist. The development of new antifungals and diagnostic tools that are both necessary to fight fungal diseases efficiently, requires greater insights in the biology of the fungal pathogens of humans in the context of the infection, on their epidemiology, and on their role in the human mycobiota. We also need a better understanding of the host immune responses to fungal pathogens as well as the genetic basis for the increased sensitivity of some individuals to fungal infections. Here, we highlight some recent progress made in these different areas of research, in particular based on work conducted in our own laboratories. These progresses should lay the ground for better management of fungal infections, as they provide opportunities for better diagnostic, vaccination, the development of classical antifungals but also strategies for targeting virulence factors or the host.

The Rhizophagus irregularis Genome Encodes Two CTR Copper Transporters That Mediate Cu Import Into the Cytosol and a CTR-Like Protein Likely Involved in Copper Tolerance

Arbuscular mycorrhizal fungi increase fitness of their host plants under Cu deficient and toxic conditions. In this study, we have characterized two Cu transporters of the CTR family (RiCTR1 and RiCTR2) and a CTR-like protein (RiCTR3A) of Rhizophagus irregularis. Functional analyses in yeast revealed that RiCTR1 encodes a plasma membrane Cu transporter, RiCTR2 a vacuolar Cu transporter and RiCTR3A a plasma membrane protein involved in Cu tolerance. RiCTR1 was more highly expressed in the extraradical mycelia (ERM) and RiCTR2 in the intraradical mycelia (IRM). In the ERM, RiCTR1 expression was up-regulated by Cu deficiency and down-regulated by Cu toxicity. RiCTR2 expression increased only in the ERM grown under severe Cu-deficient conditions. These data suggest that RiCTR1 is involved in Cu uptake by the ERM and RiCTR2 in mobilization of vacuolar Cu stores. Cu deficiency decreased mycorrhizal colonization and arbuscule frequency, but increased RiCTR1 and RiCTR2 expression in the IRM, which suggest that the IRM has a high Cu demand. The two alternatively spliced products of RiCTR3, RiCTR3A and RiCTR3B, were more highly expressed in the ERM. Up-regulation of RiCTR3A by Cu toxicity and the yeast complementation assays suggest that RiCTR3A might function as a Cu receptor involved in Cu tolerance.

Studying fungal pathogens of humans and fungal infections: fungal diversity and diversity of approaches

Seminal work by Louis Pasteur revealed the contribution of fungi-yeasts and microsporidia to agroindustry and disease in animals, respectively. More than 150 years later, the impact of fungi on human health and beyond is an ever-increasing issue, although often underestimated. Recent studies estimate that fungal infections, especially those caused by Candida, Cryptococcus and Aspergillus species, kill more than one million people annually. Indeed, these neglected infections are in general very difficult to cure and the associated mortality remains very high even when antifungal treatments exist. The development of new antifungals and diagnostic tools that are both necessary to fight fungal diseases efficiently, requires greater insights in the biology of the fungal pathogens of humans in the context of the infection, on their epidemiology, and on their role in the human mycobiota. We also need a better understanding of the host immune responses to fungal pathogens as well as the genetic basis for the increased sensitivity of some individuals to fungal infections. Here, we highlight some recent progress made in these different areas of research, in particular based on work conducted in our own laboratories. These progress should lay the ground for better management of fungal infections, as they provide opportunities for better diagnostic, vaccination, the development of classical antifungals but also strategies for targeting virulence factors or the host.

Spliceosome Profiling Visualizes Operations of a Dynamic RNP at Nucleotide Resolution

Tools to understand how the spliceosome functions in vivo have lagged behind advances in the structural biology of the spliceosome. Here, methods are described to globally profile spliceosome-bound pre-mRNA, intermediates, and spliced mRNA at nucleotide resolution. These tools are applied to three yeast species that span 600 million years of evolution. The sensitivity of the approach enables the detection of canonical and non-canonical events, including interrupted, recursive, and nested splicing. This application of statistical modeling uncovers independent roles for the size and position of the intron and the number of introns per transcript in substrate progression through the two catalytic stages. These include species-specific inputs suggestive of spliceosome-transcriptome coevolution. Further investigations reveal the ATP-dependent discard of numerous endogenous substrates after spliceosome assembly in vivo and connect this discard to intron retention, a form of splicing regulation. Spliceosome profiling is a quantitative, generalizable global technology used to investigate an RNP central to eukaryotic gene expression.

Pan-genome analyses of model fungal species

The concept of the species ‘pan-genome’, the union of ‘core’ conserved genes and all ‘accessory’ non-conserved genes across all strains of a species, was first proposed in prokaryotes to account for intraspecific variability. Species pan-genomes have been extensively studied in prokaryotes, but evidence of species pan-genomes has also been demonstrated in eukaryotes such as plants and fungi. Using a previously published methodology based on sequence homology and conserved microsynteny, in addition to bespoke pipelines, we have investigated the pan-genomes of four model fungal species: Saccharomyces cerevisiae, Candida albicans, Cryptococcus neoformans var. grubii and Aspergillus fumigatus. Between 80 and 90 % of gene models per strain in each of these species are core genes that are highly conserved across all strains of that species, many of which are involved in housekeeping and conserved survival processes. In many of these species, the remaining ‘accessory’ gene models are clustered within subterminal regions and may be involved in pathogenesis and antimicrobial resistance. Analysis of the ancestry of species core and accessory genomes suggests that fungal pan-genomes evolve by strain-level innovations such as gene duplication as opposed to wide-scale horizontal gene transfer. Our findings lend further supporting evidence to the existence of species pan-genomes in eukaryote taxa.

Transcriptome Sequencing Approaches to Elucidate Host-Microbe Interactions in Opportunistic Human Fungal Pathogens

Infections caused by opportunistic human fungal pathogens are a source of increasing medical concern, due to their growing incidence, the emergence of novel pathogenic species, and the lack of effective diagnostics tools. Fungal pathogens are phylogenetically diverse, and their virulence mechanisms can differ widely across species. Despite extensive efforts, the molecular bases of virulence in pathogenic fungi and their interactions with the human host remain poorly understood for most species. In this context, next-generation sequencing approaches hold the promise of helping to close this knowledge gap. In particular, high-throughput transcriptome sequencing (RNA-Seq) enables monitoring the transcriptional profile of both host and microbes to elucidate their interactions and discover molecular mechanisms of virulence and host defense. Here, we provide an overview of transcriptome sequencing techniques and approaches, and survey their application in studying the interplay between humans and fungal pathogens. Finally, we discuss novel RNA-Seq approaches in studying host-pathogen interactions and their potential role in advancing the clinical diagnostics of fungal infections.

Golgi Reassembly and Stacking Protein (GRASP) Participates in Vesicle-Mediated RNA Export in Cryptococcus Neoformans

Golgi reassembly and stacking protein (GRASP) is required for polysaccharide secretion and virulence in Cryptococcus neoformans. In fungal species, extracellular vesicles (EVs) participate in the export of polysaccharides, proteins and RNA. In the present work, we investigated if EV-mediated RNA export is functionally connected with GRASP in C. neoformans using a graspΔ mutant. Since GRASP-mediated unconventional secretion involves autophagosome formation in yeast, we included the atg7Δ mutant with defective autophagic mechanisms in our analysis. All fungal strains exported EVs but deletion of GRASP or ATG7 profoundly affected vesicular dimensions. The mRNA content of the graspΔ EVs differed substantially from that of the other two strains. The transcripts associated to the endoplasmic reticulum were highly abundant transcripts in graspΔ EVs. Among non-coding RNAs (ncRNAs), tRNA fragments were the most abundant in both mutant EVs but graspΔ EVs alone concentrated 22 exclusive sequences. In general, our results showed that the EV RNA content from atg7Δ and WT were more related than the RNA content of graspΔ, suggesting that GRASP, but not the autophagy regulator Atg7, is involved in the EV export of RNA. This is a previously unknown function for a key regulator of unconventional secretion in eukaryotic cells.

Conservation of mRNA quality control factor Ski7 and its diversification through changes in alternative splicing and gene duplication

Eukaryotes maintain fidelity of gene expression by preferential degradation of aberrant mRNAs that arise by errors in RNA processing reactions. In Saccharomyces cerevisiae, Ski7 plays an important role in this mRNA quality control by mediating mRNA degradation by the RNA exosome. Ski7 was initially thought to be restricted to Saccharomyces cerevisiae and close relatives because the SKI7 gene and its paralog HBS1 arose by whole genome duplication (WGD) in a recent ancestor. We have recently shown that the preduplication gene was alternatively spliced and that Ski7 function predates WGD. Here, we use transcriptome analysis of diverse eukaryotes to show that diverse eukaryotes use alternative splicing of SKI7/HBS1 to encode two proteins. Although alternative splicing affects the same intrinsically disordered region of the protein, the pattern of splice site usage varies. This alternative splicing event arose in an early eukaryote that is a common ancestor of plants, animals, and fungi. Remarkably, through changes in alternative splicing and gene duplication, the Ski7 protein has diversified such that different species express one of four distinct Ski7-like proteins. We also show experimentally that the Saccharomyces cerevisiae SKI7 gene has undergone multiple changes that are incompatible with the Hbs1 function and may also have undergone additional changes to optimize mRNA quality control. The combination of transcriptome analysis in diverse eukaryotes and genetic analysis in yeast clarifies the mechanism by which a Ski7-like protein is expressed across eukaryotes and provides a unique view of changes in alternative splicing patterns of one gene over long evolutionary time.

Advances in Cryptococcus genomics: insights into the evolution of pathogenesis

Cryptococcus species are the causative agents of cryptococcal meningitis, a significant source of mortality in immunocompromised individuals. Initial work on the molecular epidemiology of this fungal pathogen utilized genotyping approaches to describe the genetic diversity and biogeography of two species, Cryptococcus neoformans and Cryptococcus gattii. Whole genome sequencing of representatives of both species resulted in reference assemblies enabling a wide array of downstream studies and genomic resources. With the increasing availability of whole genome sequencing, both species have now had hundreds of individual isolates sequenced, providing fine-scale insight into the evolution and diversification of Cryptococcus and allowing for the first genome-wide association studies to identify genetic variants associated with human virulence. Sequencing has also begun to examine the microevolution of isolates during prolonged infection and to identify variants specific to outbreak lineages, highlighting the potential role of hyper-mutation in evolving within short time scales. We can anticipate that further advances in sequencing technology and sequencing microbial genomes at scale, including metagenomics approaches, will continue to refine our view of how the evolution of Cryptococcus drives its success as a pathogen.

Introns in Cryptococcus

In Cryptococcus neoformans, nearly all genes are interrupted by small introns. In recent years, genome annotation and genetic analysis have illuminated the major roles these introns play in the biology of this pathogenic yeast. Introns are necessary for gene expression and alternative splicing can regulate gene expression in response to environmental cues. In addition, recent studies have revealed that C. neoformans introns help to prevent transposon dissemination and protect genome integrity. These characteristics of cryptococcal introns are probably not unique to Cryptococcus, and this yeast likely can be considered as a model for intron-related studies in fungi.

Environmental Triazole Induces Cross-Resistance to Clinical Drugs and Affects Morphophysiology and Virulence of Cryptococcus gattii and C. neoformans

Cryptococcus gattii and Cryptococcus neoformans are environmental fungi that cause cryptococcosis, which is usually treated with amphotericin B and fluconazole. However, therapeutic failure is increasing because of the emergence of resistant strains. Because these species are constantly isolated from vegetal materials and the usage of agrochemicals is growing, we postulate that pesticides could be responsible for the altered susceptibility of these fungi to clinical drugs. Therefore, we evaluated the influence of the pesticide tebuconazole on the susceptibility to clinical drugs, morphophysiology, and virulence of C. gattii and C. neoformans strains. The results showed that tebuconazole exposure caused in vitro cross-resistance (CR) between the agrochemical and clinical azoles (fluconazole, itraconazole, and ravuconazole) but not with amphotericin B. In some strains, CR was observed even after the exposure ceased. Further, tebuconazole exposure changed the morphology, including formation of pseudohyphae in C. neoformans H99, and the surface charge of the cells. Although the virulence of both species previously exposed to tebuconazole was decreased in mice, the tebuconazole-exposed colonies recovered from the lungs were more resistant to azole drugs than the nonexposed cells. This in vivo CR was confirmed when fluconazole was not able to reduce the fungal burden in the lungs of mice. The tolerance to azoles could be due to increased expression of the ERG11 gene in both species and of efflux pump genes (AFR1 and MDR1) in C. neoformans Our study data support the idea that agrochemical usage can significantly affect human pathogens present in the environment by affecting their resistance to clinical drugs.

Functions for fission yeast splicing factors SpSlu7 and SpPrp18 in alternative splice-site choice and stress-specific regulated splicing

Budding yeast spliceosomal factors ScSlu7 and ScPrp18 interact and mediate intron 3'ss choice during second step pre-mRNA splicing. The fission yeast genome with abundant multi-intronic transcripts, degenerate splice signals and SR proteins is an apt unicellular fungal model to deduce roles for core spliceosomal factors in alternative splice-site choice, intron retention and to study the cellular implications of regulated splicing. From our custom microarray data we deduce a stringent reproducible subset of S. pombe alternative events. We examined the role of factors SpSlu7 or SpPrp18 for these splice events and investigated the relationship to growth phase and stress. Wild-type log and stationary phase cells showed ats1+ exon 3 skipped and intron 3 retained transcripts. Interestingly the non-consensus 5'ss in ats1+ intron 3 caused SpSlu7 and SpPrp18 dependent intron retention. We validated the use of an alternative 5'ss in dtd1+ intron 1 and of an upstream alternative 3'ss in DUF3074 intron 1. The dtd1+ intron 1 non-canonical 5'ss yielded an alternative mRNA whose levels increased in stationary phase. Utilization of dtd1+ intron 1 sub-optimal 5' ss required functional SpPrp18 and SpSlu7 while compromise in SpSlu7 function alone hampered the selection of the DUF3074 intron 1 non canonical 3'ss. We analysed the relative abundance of these splice isoforms during mild thermal, oxidative and heavy metal stress and found stress-specific splice patterns for ats1+ and DUF3074 intron 1 some of which were SpSlu7 and SpPrp18 dependent. By studying ats1+ splice isoforms during compromised transcription elongation rates in wild-type, spslu7-2 and spprp18-5 mutant cells we found dynamic and intron context-specific effects in splice-site choice. Our work thus shows the combinatorial effects of splice site strength, core splicing factor functions and transcription elongation kinetics to dictate alternative splice patterns which in turn serve as an additional recourse of gene regulation in fission yeast.

Fungal genome and mating system transitions facilitated by chromosomal translocations involving intercentromeric recombination

Species within the human pathogenic Cryptococcus species complex are major threats to public health, causing approximately 1 million annual infections globally. Cryptococcus amylolentus is the most closely known related species of the pathogenic Cryptococcus species complex, and it is non-pathogenic. Additionally, while pathogenic Cryptococcus species have bipolar mating systems with a single large mating type (MAT) locus that represents a derived state in Basidiomycetes, C. amylolentus has a tetrapolar mating system with 2 MAT loci (P/R and HD) located on different chromosomes. Thus, studying C. amylolentus will shed light on the transition from tetrapolar to bipolar mating systems in the pathogenic Cryptococcus species, as well as its possible link with the origin and evolution of pathogenesis. In this study, we sequenced, assembled, and annotated the genomes of 2 C. amylolentus isolates, CBS6039 and CBS6273, which are sexual and interfertile. Genome comparison between the 2 C. amylolentus isolates identified the boundaries and the complete gene contents of the P/R and HD MAT loci. Bioinformatic and chromatin immunoprecipitation sequencing (ChIP-seq) analyses revealed that, similar to those of the pathogenic Cryptococcus species, C. amylolentus has regional centromeres (CENs) that are enriched with species-specific transposable and repetitive DNA elements. Additionally, we found that while neither the P/R nor the HD locus is physically closely linked to its centromere in C. amylolentus, and the regions between the MAT loci and their respective centromeres show overall synteny between the 2 genomes, both MAT loci exhibit genetic linkage to their respective centromere during meiosis, suggesting the presence of recombinational suppressors and/or epistatic gene interactions in the MAT-CEN intervening regions. Furthermore, genomic comparisons between C. amylolentus and related pathogenic Cryptococcus species provide evidence that multiple chromosomal rearrangements mediated by intercentromeric recombination have occurred during descent of the 2 lineages from their common ancestor. Taken together, our findings support a model in which the evolution of the bipolar mating system was initiated by an ectopic recombination event mediated by similar repetitive centromeric DNA elements shared between chromosomes. This translocation brought the P/R and HD loci onto the same chromosome, and further chromosomal rearrangements then resulted in the 2 MAT loci becoming physically linked and eventually fusing to form the single contiguous MAT locus that is now extant in the pathogenic Cryptococcus species.

This manuscript explores the evolution of the genomic regions encoding the mating type loci of basidiomycetous fungi. Typically, the mating system is tetrapolar, meaning that it is composed of 2 unlinked mating type (MAT) loci (P/R and HD) that are located on different chromosomes. However, species with bipolar mating systems, in which the P/R and HD loci are located on the same chromosome, have also been identified. Tetrapolar and bipolar species are often closely related, suggesting the transition between these 2 mating systems might occur frequently. For example, the species within the human fungal pathogenic Cryptococcus species complex have bipolar mating systems, with 1 large MAT locus that appears to be a fusion product of the P/R and HD loci. On the other hand, the species that is the closest outgroup to these pathogenic species, Cryptococcus amylolentus, appears to have a classic tetrapolar mating system. Interestingly, the 2 MAT loci of C. amylolentus exhibit centromeric linkage during meiosis, and as a consequence, their resulting meiotic segregation pattern differs from other regions of the genome. Additionally, both pathogenic and non-pathogenic species are found to have large regional centromeres enriched with transposable and repetitive elements. Our genome comparison analyses indicated that these regional centromeres underwent ectopic recombination during the evolution of these 2 lineages. Based on these observations, we propose a model for the transition from the tetrapolar mating system in non-pathogenic C. amylolentus to the bipolar mating system in its related pathogenic species that is initiated by intercentromeric ectopic recombination, followed by chromosomal rearrangements. These events moved the 2 MAT loci closer to each other and eventually fused them to form a single MAT locus. This model is also consistent with recent findings on the organization of MAT loci in other basidiomycetous species.

Transcriptome analysis of smut fungi reveals widespread intergenic transcription and conserved antisense transcript expression

Background

Biotrophic fungal plant pathogens cause billions of dollars in losses to North American crops annually. The model for functional investigation of these fungi is Ustilago maydis. Its 20.5 Mb annotated genome sequence has been an excellent resource for investigating biotrophic plant pathogenesis. Expressed-sequence tag libraries and microarray hybridizations have provided insight regarding the type of transcripts produced by U. maydis but these analyses were not comprehensive and there were insufficient data for transcriptome comparison to other smut fungi. To improve transcriptome annotation and enable comparative analyses, comprehensive strand-specific RNA-seq was performed on cell-types of three related smut species: U. maydis (common smut of corn), Ustilago hordei (covered smut of barley), and Sporisorium reilianum (head smut of corn).

Results

In total, >1 billion paired-end sequence reads were obtained from haploid cell, dikaryon and teliospore RNA of U. maydis, haploid cell RNA of U. hordei, and haploid and dikaryon cell RNA of S. reilianum. The sequences were assembled into transfrags using Trinity, and updated gene models were created using PASA and categorized with Cufflinks Cuffcompare. Representative genes that were predicted for the first time with these RNA-seq analyses and genes with novel annotation features were independently assessed by reverse transcriptase PCR. The analyses indicate hundreds more predicted proteins, relative to the previous genome annotation, could be produced by U. maydis from altered transcript forms, and that the number of non-coding RNAs produced, including transcribed intergenic sequences and natural antisense transcripts, approximately equals the number of mRNAs. This high representation of non-coding RNAs appears to be a conserved feature of the smut fungi regardless of whether they have RNA interference machinery. Approximately 50% of the identified NATs were conserved among the smut fungi.

Conclusions

Overall, these analyses revealed: 1) smut genomes encode a number of transcriptional units that is twice the number of annotated protein-coding genes, 2) a small number of intergenic transcripts may encode proteins with characteristics of fungal effectors, 3) the vast majority of intergenic and antisense transcripts do not contain ORFs, 4) a large proportion of the identified antisense transcripts were detected at orthologous loci among the smut fungi, and 5) there is an enrichment of functional categories among orthologous loci that suggests antisense RNAs could have a genome-wide, non-RNAi-mediated, influence on gene expression in smut fungi.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3720-8) contains supplementary material, which is available to authorized users.

Transcriptome analysis reveals the complexity of alternative splicing regulation in the fungus Verticillium dahliae

Background

Alternative splicing (AS) regulation is extensive and shapes the functional complexity of higher organisms. However, the contribution of alternative splicing to fungal biology is not well studied.

Results

This study provides sequences of the transcriptomes of the plant wilt pathogen Verticillium dahliae, using two different strains and multiple methods for cDNA library preparations. We identified alternatively spliced mRNA isoforms in over a half of the multi-exonic fungal genes. Over one-thousand isoforms involve TopHat novel splice junction; multiple types of combinatory alternative splicing patterns were identified. We showed that one Verticillium gene could use four different 5′ splice sites and two different 3′ donor sites to produce up to five mature mRNAs, representing one of the most sophisticated alternative splicing model in eukaryotes other than animals. Hundreds of novel intron types involving a pair of new splice sites were identified in the V. dahliae genome. All the types of AS events were validated by using RT-PCR. Functional enrichment analysis showed that AS genes are involved in most known biological functions and enriched in ATP biosynthesis, sexual/asexual reproduction, morphogenesis, signal transduction etc., predicting that the AS regulation modulates mRNA isoform output and shapes the V. dahliae proteome plasticity of the pathogen in response to the environmental and developmental changes.

Conclusions

These findings demonstrate the comprehensive alternative splicing mechanisms in a fungal plant pathogen, which argues the importance of this fungus in developing complicate genome regulation strategies in eukaryotes.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3507-y) contains supplementary material, which is available to authorized users.