SnoRNAs from the filamentous fungus Neurospora crassa: structural, functional and evolutionary insights

The fungal snoRNAome

Small nucleolar RNAs (snoRNAs) are essential players in the rRNA biogenesis due to their involvement in the nucleolytic processing of the precursor and the subsequent guidance of nucleoside modifications. Within the kingdom Fungi, merely a few species-specific surveys have explored their snoRNA repertoire. However, the wide range of the snoRNA landscape spanning all major fungal lineages has not been mapped so far, mainly because of missing tools for automatized snoRNA detection and functional analysis. For the first time, we report here a comprehensive inventory of fungal snoRNAs together with a functional analysis and an in-depth investigation of their evolutionary history including innovations, deletions, and target switches. This large-scale analysis, incorporating more than 120 snoRNA families with more than 7700 individual snoRNA sequences, catalogs and clarifies the landscape of fungal snoRNA families, assigns functions to previously orphan snoRNAs, and increases the number of sequences by 450%. We also show that the snoRNAome is subject to ongoing rearrangements and adaptations, e.g., through lineage-specific targets and redundant guiding functions.

Identification of 15 candidate structured noncoding RNA motifs in fungi by comparative genomics

BACKGROUND

With the development of rapid and inexpensive DNA sequencing, the genome sequences of more than 100 fungal species have been made available. This dataset provides an excellent resource for comparative genomics analyses, which can be used to discover genetic elements, including noncoding RNAs (ncRNAs). Bioinformatics tools similar to those used to uncover novel ncRNAs in bacteria, likewise, should be useful for searching fungal genomic sequences, and the relative ease of genetic experiments with some model fungal species could facilitate experimental validation studies.

RESULTS

We have adapted a bioinformatics pipeline for discovering bacterial ncRNAs to systematically analyze many fungal genomes. This comparative genomics pipeline integrates information on conserved RNA sequence and structural features with alternative splicing information to reveal fungal RNA motifs that are candidate regulatory domains, or that might have other possible functions. A total of 15 prominent classes of structured ncRNA candidates were identified, including variant HDV self-cleaving ribozyme representatives, atypical snoRNA candidates, and possible structured antisense RNA motifs. Candidate regulatory motifs were also found associated with genes for ribosomal proteins, S-adenosylmethionine decarboxylase (SDC), amidase, and HexA protein involved in Woronin body formation. We experimentally confirm that the variant HDV ribozymes undergo rapid self-cleavage, and we demonstrate that the SDC RNA motif reduces the expression of SAM decarboxylase by translational repression. Furthermore, we provide evidence that several other motifs discovered in this study are likely to be functional ncRNA elements.

CONCLUSIONS

Systematic screening of fungal genomes using a computational discovery pipeline has revealed the existence of a variety of novel structured ncRNAs. Genome contexts and similarities to known ncRNA motifs provide strong evidence for the biological and biochemical functions of some newly found ncRNA motifs. Although initial examinations of several motifs provide evidence for their likely functions, other motifs will require more in-depth analysis to reveal their functions.

Phylogenetic distribution of plant snoRNA families

Background

Small nucleolar RNAs (snoRNAs) are one of the most ancient families amongst non-protein-coding RNAs. They are ubiquitous in Archaea and Eukarya but absent in bacteria. Their main function is to target chemical modifications of ribosomal RNAs. They fall into two classes, box C/D snoRNAs and box H/ACA snoRNAs, which are clearly distinguished by conserved sequence motifs and the type of chemical modification that they govern. Similarly to microRNAs, snoRNAs appear in distinct families of homologs that affect homologous targets. In animals, snoRNAs and their evolution have been studied in much detail. In plants, however, their evolution has attracted comparably little attention.

Results

In order to chart the phylogenetic distribution of individual snoRNA families in plants, we applied a sophisticated approach for identifying homologs of known plant snoRNAs across the plant kingdom. In response to the relatively fast evolution of snoRNAs, information on conserved sequence boxes, target sequences, and secondary structure is combined to identify additional snoRNAs. We identified 296 families of snoRNAs in 24 species and traced their evolution throughout the plant kingdom. Many of the plant snoRNA families comprise paralogs. We also found that targets are well-conserved for most snoRNA families.

Conclusions

The sequence conservation of snoRNAs is sufficient to establish homologies between phyla. The degree of this conservation tapers off, however, between land plants and algae. Plant snoRNAs are frequently organized in highly conserved spatial clusters. As a resource for further investigations we provide carefully curated and annotated alignments for each snoRNA family under investigation.

Electronic supplementary material

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

Identification of Non-Coding RNAs in the Candida parapsilosis Species Group

The Candida CTG clade is a monophyletic group of fungal species that translates CTG as serine, and includes the pathogens Candida albicans and Candida parapsilosis. Research has typically focused on identifying protein-coding genes in these species. Here, we use bioinformatic and experimental approaches to annotate known classes of non-coding RNAs in three CTG-clade species, Candida parapsilosis, Candida orthopsilosis and Lodderomyces elongisporus. We also update the annotation of ncRNAs in the C. albicans genome. The majority of ncRNAs identified were snoRNAs. Approximately 50% of snoRNAs (including most of the C/D box class) are encoded in introns. Most are within mono- and polycistronic transcripts with no protein coding potential. Five polycistronic clusters of snoRNAs are highly conserved in fungi. In polycistronic regions, splicing occurs via the classical pathway, as well as by nested and recursive splicing. We identified spliceosomal small nuclear RNAs, the telomerase RNA component, signal recognition particle, RNase P RNA component and the related RNase MRP RNA component in all three genomes. Stem loop IV of the U2 spliceosomal RNA and the associated binding proteins were lost from the ancestor of C. parapsilosis and C. orthopsilosis, following the divergence from L. elongisporus. The RNA component of the MRP is longer in C. parapsilosis, C. orthopsilosis and L. elongisporus than in S. cerevisiae, but is substantially shorter than in C. albicans.

Identification and Characterization of Small Noncoding RNAs in Genome Sequences of the Edible Fungus Pleurotus ostreatus

Noncoding RNAs (ncRNAs) have been identified in many fungi. However, no genome-scale identification of ncRNAs has been inventoried for basidiomycetes. In this research, we detected 254 small noncoding RNAs (sncRNAs) in a genome assembly of an isolate (CCEF00389) of Pleurotus ostreatus, which is a widely cultivated edible basidiomycetous fungus worldwide. The identified sncRNAs include snRNAs, snoRNAs, tRNAs, and miRNAs. SnRNA U1 was not found in CCEF00389 genome assembly and some other basidiomycetous genomes by BLASTn. This implies that if snRNA U1 of basidiomycetes exists, it has a sequence that varies significantly from other organisms. By analyzing the distribution of sncRNA loci, we found that snRNAs and most tRNAs (88.6%) were located in pseudo-UTR regions, while miRNAs are commonly found in introns. To analyze the evolutionary conservation of the sncRNAs in P. ostreatus, we aligned all 254 sncRNAs to the genome assemblies of some other Agaricomycotina fungi. The results suggest that most sncRNAs (77.56%) were highly conserved in P. ostreatus, and 20% were conserved in Agaricomycotina fungi. These findings indicate that most sncRNAs of P. ostreatus were not conserved across Agaricomycotina fungi.

Endogenous Small RNA Mediates Meiotic Silencing of a Novel DNA Transposon

Genome defense likely evolved to curtail the spread of transposable elements and invading viruses. A combination of effective defense mechanisms has been shown to limit colonization of the Neurospora crassa genome by transposable elements. A novel DNA transposon named Sly1-1 was discovered in the genome of the most widely used laboratory "wild-type" strain FGSC 2489 (OR74A). Meiotic silencing by unpaired DNA, also simply called meiotic silencing, prevents the expression of regions of the genome that are unpaired during karyogamy. This mechanism is posttranscriptional and is proposed to involve the production of small RNA, so-called masiRNAs, by proteins homologous to those involved in RNA interference-silencing pathways in animals, fungi, and plants. Here, we demonstrate production of small RNAs when Sly1-1 was unpaired in a cross between two wild-type strains. These small RNAs are dependent on SAD-1, an RNA-dependent RNA polymerase necessary for meiotic silencing. We present the first case of endogenously produced masiRNA from a novel N. crassa DNA transposable element.

Conservation and divergence of transcriptional coregulations between box C/D snoRNA and ribosomal protein genes in Ascomycota

Coordinated assembly of the ribosome is essential for proper translational activity in eukaryotic cells. It is therefore critical to coordinate the expression of components of ribosomal programs with the cell's nutritional status. However, coordinating expression of these components is poorly understood. Here, by combining experimental and computational approaches, we systematically identified box C/D snoRNAs in four fission yeasts and found that the expression of box C/D snoRNA and ribosomal protein (RP) genes were orchestrated by a common Homol-D box, thereby ensuring a constant balance of these two genetic components. Interestingly, such transcriptional coregulations could be observed in most Ascomycota species and were mediated by different cis-regulatory elements. Via the reservation of cis elements, changes in spatial configuration, the substitution of cis elements, and gain or loss of cis elements, the regulatory networks of box C/D snoRNAs evolved to correspond with those of the RP genes, maintaining transcriptional coregulation between box C/D snoRNAs and RP genes. Our results indicate that coregulation via common cis elements is an important mechanism to coordinate expression of the RP and snoRNA genes, which ensures a constant balance of these two components.

Identification and characterisation of non-coding small RNAs in the pathogenic filamentous fungus Trichophyton rubrum

BACKGROUND

Accumulating evidence demonstrates that non-coding RNAs (ncRNAs) are indispensable components of many organisms and play important roles in cellular events, regulation, and development.

RESULTS

Here, we analysed the small non-coding RNA (ncRNA) transcriptome of Trichophyton rubrum by constructing and sequencing a cDNA library from conidia and mycelia. We identified 352 ncRNAs and their corresponding genomic loci. These ncRNA candidates included 198 entirely novel ncRNAs and 154 known ncRNAs classified as snRNAs, snoRNAs and other known ncRNAs. Further bioinformatic analysis detected 96 snoRNAs, including 56 snoRNAs that had been annotated in other organisms and 40 novel snoRNAs. All snoRNAs belonged to two major classes--C/D box snoRNAs and H/ACA snoRNAs--and their potential target sites in rRNAs and snRNAs were predicted. To analyse the evolutionary conservation of the ncRNAs in T. rubrum, we aligned all 352 ncRNAs to the genomes of six dermatophytes and to the NCBI non-redundant nucleotide database (NT). The results showed that most of the identified snRNAs were conserved in dermatophytes. Of the 352 ncRNAs, 102 also had genomic loci in other dermatophytes, and 27 were dermatophyte-specific.

CONCLUSIONS

Our systematic analysis may provide important clues to the function and evolution of ncRNAs in T. rubrum. These results also provide important information to complement the current annotation of the T. rubrum genome, which primarily comprises protein-coding genes.

Transcriptome analysis of the filamentous fungus Aspergillus nidulans directed to the global identification of promoters

BACKGROUND

The filamentous fungus Aspergillus nidulans has been a tractable model organism for cell biology and genetics for over 60 years. It is among a large number of Aspergilli whose genomes have been sequenced since 2005, including medically and industrially important species. In order to advance our knowledge of its biology and increase its utility as a genetic model by improving gene annotation we sequenced the transcriptome of A. nidulans with a focus on 5' end analysis.

RESULTS

Strand-specific whole transcriptome sequencing showed that 80-95% of annotated genes appear to be expressed across the conditions tested. We estimate that the total gene number should be increased by approximately 1000, to 11,800. With respect to splicing 8.3% of genes had multiple alternative transcripts, but alternative splicing by exon-skipping was very rare. 75% of annotated genes showed some level of antisense transcription and for one gene, meaB, we demonstrated the antisense transcript has a regulatory role. Specific sequencing of the 5' ends of transcripts was used for genome wide mapping of transcription start sites, allowing us to interrogate over 7000 promoters and 5' untranslated regions.

CONCLUSIONS

Our data has revealed the complexity of the A. nidulans transcriptome and contributed to improved genome annotation. The data can be viewed on the AspGD genome browser.

Massively parallel sequencing technology in pathogenic microbes

Next-Generation Sequencing (NGS) methods have revolutionized various aspects of genomics including transcriptome analysis. Digital expression analysis is all set to replace analog expression analysis that uses microarray chips through their cost-effectiveness, reproducibility, accuracy, and speed. The last 2 years have seen a surge in the development of statistical methods and software tools for analysis and visualization of NGS data. Large amounts of NGS data are available for pathogenic fungi and oomycetes. As the analysis results start pouring in, it brings about a paradigm shift in the understanding of host pathogen interactions with discovery of new transcripts, splice variants, mutations, regulatory elements, and epigenetic controls. Here we describe the core technology of the new sequencing platforms, the methodology of data analysis, and different aspects of applications.

Box C/D snoRNP catalysed methylation is aided by additional pre-rRNA base-pairing

2'-O-methylation of eukaryotic ribosomal RNA (r)RNA, essential for ribosome function, is catalysed by box C/D small nucleolar (sno)RNPs. The RNA components of these complexes (snoRNAs) contain one or two guide sequences, which, through base-pairing, select the rRNA modification site. Adjacent to the guide sequences are protein-binding sites (the C/D or C'/D' motifs). Analysis of >2000 yeast box C/D snoRNAs identified additional conserved sequences in many snoRNAs that are complementary to regions adjacent to the rRNA methylation site. This 'extra base-pairing' was also found in many human box C/D snoRNAs and can stimulate methylation by up to five-fold. Sequence analysis, combined with RNA-protein crosslinking in Saccharomyces cerevisiae, identified highly divergent box C'/D' motifs that are bound by snoRNP proteins. In vivo rRNA methylation assays showed these to be active. Our data suggest roles for non-catalytic subunits (Nop56 and Nop58) in rRNA binding and support an asymmetric model for box C/D snoRNP organization. The study provides novel insights into the extent of the snoRNA-rRNA interactions required for efficient methylation and the structural organization of the snoRNPs.

Evolution of yeast noncoding RNAs reveals an alternative mechanism for widespread intron loss

The evolutionary forces responsible for intron loss are unresolved. Whereas research has focused on protein-coding genes, here we analyze noncoding small nucleolar RNA (snoRNA) genes in which introns, rather than exons, are typically the functional elements. Within the yeast lineage exemplified by the human pathogen Candida albicans, we find--through deep RNA sequencing and genome-wide annotation of splice junctions--extreme compaction and loss of associated exons, but retention of snoRNAs within introns. In the Saccharomyces yeast lineage, however, we find it is the introns that have been lost through widespread degeneration of splicing signals. This intron loss, perhaps facilitated by innovations in snoRNA processing, is distinct from that observed in protein-coding genes with respect to both mechanism and evolutionary timing.

SnoPatrol: how many snoRNA genes are there?

Small nucleolar RNAs (snoRNAs) are among the most evolutionarily ancient classes of small RNA. Two experimental screens published in BMC Genomics expand the eukaryotic snoRNA catalog, but many more snoRNAs remain to be found.

See research articles http://www.biomedcentral.com/1471-2164/10/515 and http://www.biomedcentral.com/1471-2164/11/61.