Identification of 15 candidate structured noncoding RNA motifs in fungi by comparative genomics
BMC Genomics volume 18, Article number: 785 (2017)
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.
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.
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.
The fungal kingdom is an abundantly represented division of the eukaryotic domain of life. It has been estimated that over 1.5 million fungal species could exist, although only about 74,000 have been documented [1, 2]. Numerous species of fungi are major contributors to the composition of many ecosystems where they play important roles in the environment ranging from the simple promotion of organic material decomposition to the formation of essential symbiotic relationships with bacteria, plants and animals . Fungi also serve as important sources for certain human foods, medicines and industrial agents. However, fungi can also be pathogens of many important crop plants, agricultural animals, and of humans. Given the wide-ranging importance of fungal species, deeper explorations of their conserved genetic elements are warranted.
The emergence of inexpensive and efficient DNA sequencing methods has facilitated the determination of the genomic DNA sequences of many fungal species. The yeast Saccharomyces cerevisiae provided the first genome of a eukaryote to be completely sequenced . This was followed quickly by the sequencing of additional fungal genomes, including Schizosaccharomyces pombe , and Neurospora crassa . N. crassa was the first filamentous species to have its genome sequenced, and it is estimated to possess about double the number of genes predicted to be present in the 6,000-gene S. cerevisiae. So far more than 100 fungal species have been fully sequenced . These genome sequences serve as excellent resources for conducting comparative genomics analyses to find common features among diverse fungal species.
Of particular interest to us are structured ncRNAs, such as riboswitches  and ribozymes , which perform gene regulation, RNA processing, or other important biochemical functions. Such ncRNAs are likely to remain hidden in the genomes of organisms because there is no simple strategy known that can be implemented to comprehensively predict their existence. By contrast, many protein-coding genes contain long open reading frames (ORFs) that are easily revealed by computer-assisted search algorithms. Initially, simple clues regarding the existence of ncRNAs were used to identify novel RNA motifs in bacteria. For example, large gaps between protein-coding regions along with high sequence conservation among closely related bacterial species typically are signatures of ncRNAs [10, 11]. Similarly, orphan promoter or terminator sequences not immediately associated with an ORF [10, 12], GC-rich regions in an organism with a high AT genomic sequence composition , and conserved RNA secondary structures [14,15,16] also have been used to identify structured ncRNAs.
A variety of complex approaches and algorithms that exploit these and other characteristics have since been developed to search for and annotate ncRNAs in species from all three domains of life [17, 18]. Some of these algorithms have been successfully applied to identify additional representatives of known ncRNAs and to discover novel classes present in fungal genomes [12, 19,20,21,22,23]. Even from the earliest computer-assisted searches for ncRNAs in fungi conducted in the 1990s, evidence emerged that many such RNAs remained to be discovered. For example, by detecting RNA polymerase III transcripts and analyzing large noncoding gaps present in the S. cerevisiae genome, 16 candidate ncRNAs were identified . Soon thereafter, a total of 22 novel small nucleolar RNAs (snoRNAs) that guide RNA methylation were identified from this same genome by using comparative sequence analysis methods . More recently, RNA transcriptomics methods have been used to identify dozens of additional ncRNAs from fungal species such as Aspergillus fumigatus  and Trichophyton rubrum . These findings suggest that additional ncRNA classes remain to be discovered among the many diverse species of fungi.
Each of the previous computational searches noted above largely focused on a single fungal genome, or a few very closely related genomes. However, given the growth in fungal DNA sequence databases, we sought to identify structured ncRNA candidates that are present in a large diversity of fungal species. Most widely-distributed ncRNAs tend to carry only short stretches of conserved sequences interspersed among structural elements that are only poorly conserved at the nucleotide level . As a result, some successful computational search strategies for structured ncRNA candidates use algorithms such as QRNA , RNAz , or CMfinder , which search for nucleotide covariation or other evidence of structure formation and conservation, rather than just searching for conserved sequence. A computational search across numerous fungal species increases the opportunity to observe nucleotide covariations indicative of structure conservation.
In the current study, we used a computational pipeline similar to that used previously to identify numerous structured ncRNAs in bacteria [29,30,31,32]. This pipeline relies on CMfinder , which uses an expectation maximization algorithm analogous to MEME  to identify probable RNA secondary structures from unaligned sequences via covariance models [26, 34]. To improve the accuracy of consensus secondary structure predictions, a thermodynamic model for base-pairing prediction is also included in CMfinder. Our searches were limited to noncoding regions of 2811 fungal genomes (see the databases for species names), including both intergenic regions and intronic regions. We believe this to be the ideal place to search because eukaryotic TPP riboswitches (the only riboswitch class found in fungi and plants) are commonly located in introns where they regulate alternative splicing upon binding this essential coenzyme .
Our pipeline revealed the existence of 15 classes of candidate structured ncRNAs (for a complete list of motifs and species, see Additional file 1: Table S1). These motifs include novel HDV ribozyme variants, non-typical snoRNAs, two structured antisense RNA motifs associated with the chromatin remodeling complex, and various other RNA motifs associated with genes coding for SAM decarboxylase, amidase, and Woronin body proteins. We experimentally confirmed that HDV variants are self-cleaving ribozymes that exhibit catalytic characteristics similar to known HDV ribozymes in other organisms. We also confirmed that the SDC motif is a translational repressor that coordinates the expression of a small upstream open reading frame (uORF) to regulate SAM decarboxylase biosynthesis. Reverse transcription and polymerase chain reaction (RT-PCR) analyses of intronic examples reveal that they likely are involved in the regulation of alternative splicing. These findings demonstrate that numerous well-conserved structured RNA elements are used by a diverse collection of fungi to serve important biochemical roles. Furthermore, this ncRNA discovery pipeline could be used to search through the genomes of other eukaryotic species to discover additional novel structured ncRNAs.
Results and discussion
Identification of candidate structured ncRNAs
Promising ncRNA motifs were identified by initially applying our computational pipeline to fungal genome sequences present in NCBI RefSeq Release 29 , and later data from RefSeq Release 62 were incorporated into this study. Refinement of the list of representatives, sequence alignments and conserved sequence/structural models were completed with the sequence database as updated in RefSeq Release 75.
Briefly, the workflow for discovering structured ncRNAs (Fig. 1) involved extracting noncoding fungal DNA sequences as informed by pre-existing genome annotations, clustering of similar sequence regions by using BLAST , and filtering to remove clusters matching known ncRNAs present in Rfam  or to remove those with extensive protein-coding potential by using RNACode . This process yielded many “pre-candidate” motifs that required further analysis to assess their relative likelihoods of functioning as common fungal ncRNAs. Each pre-candidate was subjected to analysis by CMfinder, which in part uncovers evidence for nucleotide sequence covariation that is indicative of secondary structure formation. If nucleotides in a predicted base-paired stem frequently co-vary in a manner that suggests conservation of the stem, then the cluster of RNAs was considered a strong candidate for having ncRNA function. Iterative analysis of conserved sequences and substructures, augmented by the discovery of additional representatives, was conducted using Infernal  to yield a refined sequence and structure model for the ncRNA candidate.
Analysis of candidate structured ncRNAs
Each ncRNA candidate was evaluated based on the complexity of its conserved sequence and structure model, its common genomic locations, and its phylogenetic distribution, among other features. Below are described the most promising candidate motifs that likely serve as functional ncRNAs. The alignments of these high-ranking motifs in Stockholm format are presented as supplemental material (Additional files 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15 and 16).
HDV self-cleaving ribozyme variants
The first member of the HDV ribozyme class [40, 41] to be reported was identified in the antigenomic sequence of the Hepatitis Delta Virus . Since this initial HDV ribozyme discovery, numerous self-cleaving RNAs that conform to this same general consensus sequence and secondary structure model have been discovered [43,44,45]. Despite the relatively large size of this ribozyme, there are only a few highly conserved nucleotides interspersed among the four extended base-paired regions that form its unique nested pseudoknot architecture  (Fig. 2a). The vast majority of HDV ribozyme representatives have been discovered by using a computational strategy that only seeks RNAs conforming to this unique secondary structure, while largely ignoring the few conserved nucleotides [44, 45]. Thus, bioinformatics search methods that also take into consideration conserved nucleotides have the potential to uncover additional representatives.
Our computational pipeline has uncovered a total of 230 representatives of a motif with considerable similarity in sequence and secondary structure to the known members of the HDV self-cleaving ribozyme class (Fig. 2b). These RNAs are present among 26 fungal species. The most noteworthy differences between these new RNAs and the previously published consensus (Fig. 2a) are (i) a C-G base-pair typically represents the otherwise more general R-Y base-pair at the ribozyme cleavage site, (ii) the P4 stem is partially replaced by an E loop [47,48,49] RNA motif, and (iii) an additional nucleotide (most commonly a C residue) is inserted between the P2 and P3 stems. The latter two differences violate the criteria used previously [41, 44, 45] to find additional members of the HDV ribozyme class, which in part explains why these distinctive fungal examples remained undiscovered. Only 11 fungal HDV ribozyme examples that conform to the published consensus had been discovered previously, suggesting that members of the unusual variant type revealed by our search method predominate in fungi.
Given the distinctive features of the newly-found fungal representatives of this motif, and given their widely variable gene associations, we chose to determine whether members can cleave RNA. Bimolecular substrate-ribozyme complexes were constructed for two examples derived from different organisms. For example, a bimolecular construct based on a representative from Penicillium chrysogenum (Fig. 2c) promotes cleavage of the RNA substrate strand at the base of the P1 stem (Fig. 2d). This cleavage site precisely matches that expected for HDV ribozymes, based on the similarity between the predicted structure of the fungal RNAs and the previously-established architecture of self-cleaving ribozymes belonging to the HDV class .
Similar results were observed for another bimolecular RNA construct derived from the fungal species Aspergillus niger (Fig. 3a). Again, the cleavage site for this fungal HDV ribozyme variant corresponds to that expected for more conventional HDV ribozymes as determined by gel mobility (Fig. 3b) and analysis by mass spectrometry (Fig. 3c) of the substrate cleavage products. Moreover, the masses of the products are consistent with a ribozyme mechanism wherein a 2ˊ-oxygen atom serves as a nucleophile to attack the adjacent phosphorus center to yield cleavage products with 2ˊ, 3ˊ-cyclic phosphate and 5ˊ hydroxyl termini on the 5ˊ and 3ˊ cleavage fragments, respectively. This is the same mechanism used by all other small self-cleaving ribozymes discovered to date .
The SDC motif, represented by 34 distinct examples from 26 fungal species, typically forms a small hairpin with an 11-base-pair stem (Fig. 4a). However, SDC motifs can exhibit some variation in stem integrity and length (e.g. see Fig. 4b for the representative from N. crassa). Each representative is located immediately upstream of a SAM (S-adenosylmethionine) decarboxylase gene, suggesting that hairpin formation and the associated conserved nucleotides are important for regulation of the SDC gene. In N. crassa, there appear to be two uORFs (Fig. 4c), which are short peptide-coding regions located in some mRNAs a short distance upstream of a main open reading frame. Such uORF regions are commonly involved in controlling translation initiation of the adjoining gene [51, 52]. For example, in N. crassa, high arginine concentrations cause ribosomes to increase stalling within the arg-2 uORF, which reduces translation initiation at the main ORF located immediately downstream [53, 54]. Both arginine and newly synthesized uORF-derived peptides are required for ribosomal stalling.
The SDC enzyme catalyzes the synthesis of S-adenosylmethionine to provide propylamine for the production of polyamines, such as spermidine and spermine [55, 56]. The production of polyamines is known to be highly regulated by cells at multiple levels including transcription, translation, enzyme activation and protein degradation. At the translational level, mammalian and plant SDC expression is primarily regulated by the expression of uORFs located in their 5´ UTRs [55, 56], and to a lesser extent by secondary structures formed by the 5ˊ UTR . Indeed, polyamine-responsive gene control very commonly involves uORF elements .
The two uORFs associated with the SDC motif of N. crassa are located nearby two consensus 5ˊ splice sites and one 3ˊ splice site (Fig. 4c). Given that some TPP riboswitches in fungi control both alternative splicing and uORF expression , we initially speculated that the SDC motif might be part of a complex set of regulatory elements that similarly controls SDC gene expression. To examine the biological function of the SDC motif, the 5′ UTR of the N. crassa SDC gene including a portion of the SDC ORF was fused in-frame with a luciferase reporter gene. This construct (Pre, Fig. 4c) was transformed into N. crassa cells and the presence of various RNA transcripts was examined by RT-PCR (Fig. 4d). The mRNA precursor was found to be efficiently alternatively spliced into two predominant forms called Sp-I, which contains two predicted uORFs, and Sp-II, which lacks these two uORFs (though the 3ˊ half of uORF 2 is still present) (Fig 4c, d). Therefore, any effects on expression of the SDC ORF possibly caused by the uORFs cannot occur with the Sp-II form of the processed mRNA.
We then assessed the importance of the SDC motif hairpin to gene expression by examining luciferase activity in fungal preparations carrying the WT and M1 through M4 reporter-fusion constructs. Disruption of the P1 stem by changing one (M1), two (M3), or three (M5) base-pairs to mismatches causes a substantial increase in reporter gene expression relative to that observed with the WT construct (Fig. 4e). By introducing additional mutations that restore base-pairing (M2, M4), these mutant constructs exhibit gene expression levels that approach that of the WT construct. These findings indicate that the SDC motif behaves as a negative regulatory element to inhibit SDC gene expression.
RT-PCR analysis on the variant constructs revealed that the disruptive SDC motif mutations M1 and M3 do not influence the levels of Sp-I and Sp-II mRNA products (data not shown), suggesting that the SDC motif regulates gene expression in a manner that does not directly involve alternative splicing. Moreover, we did not observe evidence for specific binding of spermine or spermidine to a representative of the RNA motif by using in-line probing [59, 60] assays (data not shown). Therefore, the mechanism of regulation by the SDC motif appears to differ from that observed for some fungal TPP riboswitches that regulate alternative splicing and the translation of uORFs . Additional studies will be necessary to determine (i) if alternative splicing of the SDC precursor RNA is regulated, (ii) if uORF-mediated regulation of SDC gene expression occurs, (iii) how the SDC motif hairpin suppresses gene expression, and (iv) how SDC motif structure can be naturally manipulated to affect expression.
A total of 23 unique examples of the amd motif (Fig. 5) have been identified among 20 fungal species. The consensus sequence and secondary structure model based on these sequences reveal that the motif likely adopts an elongated two-stem junction wherein the most-highly conserved nucleotides reside in the internal loop and in sections of P1 and P2 nearest to this loop. Moreover, the nucleotides immediately upstream of P1, along with several others extending into the internal loop, appear to code for a short uORF. This indicates that the ORF and overlapping RNA structure of the amd motif might collaboratively regulate the downstream gene. Of the 23 unique amd motif examples, nine members carry a uORF that codes for eight amino acids, and 14 members carry a uORF that codes for only six amino acids (Fig. 5). The putative peptides are highly conserved, and the last five amino acids of both groups carry the sequence A(V/F)(A/V)EL.
Each amd motif example resides upstream of a protein-coding gene of unknown function, making it difficult to formulate hypotheses regarding the function and mechanism of this putative ncRNA element and associated uORF. In more than nine species, the associated gene is annotated as coding for an amidase enzyme, whose function is to hydrolyze amide functional groups to yield ammonia and a carboxylic acid group. Given the presence of a glutamate codon in the penultimate position of the uORF, the amidase activity might be related to the production of this amino acid. For example, a shortage of glutamine might cause the uORF system to trigger expression of a glutaminase ORF to convert glutamine to glutamate. A previous attempt to create an amd gene (NCU05182) knockout in N. crassa was unsuccessful . The heterokaryotic strain, but not the homokaryotic strain, survived, indicating that the gene associated with the amd motif is critical for cell survival. Given this experimental complication, we did not pursue additional validation studies.
The ies6 motif (Fig. 6a) includes 34 examples from 20 species of fungi. This RNA motif forms an extended hairpin structure with one small and one large internal loop. The hairpin loop conforms to a GNRA tetraloop sequence , which is frequently found in structured RNAs. GNRA tetraloops can form tertiary interactions with tetraloop receptor structures [64, 65]. This fact, coupled with the extensive conservation on each side of the hairpin structure suggests that the consensus motif as presented might represent only a portion of a more complex RNA architecture.
These predicted RNA structures are almost always located adjacent to a gene coding for a protein similar to the chromatin remodeling complex subunit Ies6 (or Ino80 Subunit 6), which is involved in chromatin modification, chromosome segregation and the regulation of telomere length [66,67,68]. However, the polarity of the RNA motif as depicted (Fig. 6a) is opposite to that of the ies6 mRNA, because this configuration retains a GNRA tetraloop commonly found in structured RNAs. Therefore we speculated that this motif is likely present in an antisense RNA produced from the same genomic location as the ies6 gene. To determine if both sense and antisense RNAs corresponding to the ies6 motif were produced by N. crassa cells, we conducted RT-PCR assays. As expected, RT-PCR product bands corresponding to both the sense and antisense RNA transcripts were observed by gel electrophoresis (Fig. 6b).
Antisense transcripts are found in many species where they often participate in the gene regulation. They are able to base-pair to their targets to regulate gene expression through various mechanisms including transcription attenuation, translation inhibition, primer maturation inhibition, splicing regulation, nuclear retention, and mRNA degradation and stabilization [69,70,71,72]. However, additional experiments will be required to determine if the ies6 motif is indeed part of an antisense RNA transcript, and what role this motif might play in regulation of antisense production and action.
The hexA motif (Fig. 7a) includes 19 examples from 10 fungal species. The predicted secondary structure includes a long hairpin interrupted by two small internal loops. In addition, a long region with sequence conservation lacking evidence of structure extends downstream. The motif is found in the putative 5´ UTRs or introns of genes annotated as encoding either hypothetical proteins (40%) or Woronin body [73, 74] protein Hex subunits (60%). Woronin bodies are fungal-specific organelles that plug the septal pores quickly to prevent a cell from losing its contents during physical damage. Hex subunit proteins are derived from differently spliced mRNA forms and are one of the major and essential components of the Woronin body [74, 75].
Intriguingly, a portion of the sequence forming the right shoulder of the hexA motif hairpin closely approximates the consensus for a 3ˊ splice site (Fig. 7a). Moreover, nucleotides within the long 3ˊ tail exhibit evidence that they code for protein. Specifically, the least conserved nucleotides in this region occur at every third position, which might correspond to the wobble position of codons. These characteristics suggest that the hexA motif RNA structure might regulate splicing to create an altered ORF sequence for translation. However, at least seven of the representatives of the hexA motif are in the opposite direction as the genes encoding Woronin body proteins, whereas at least ten are in the same orientation. This variability in orientation weakens the hypothesis that the motif might be involved in gene regulation by controlling alternative splicing of the adjacent coding region.
The SART-1 motif consensus model (Fig. 7b) is based on 12 unique examples from 11 fungal species. The secondary structure can potentially include at least two hairpins, although only P1 is supported by extensive evidence of covariation and the frequent presence of a UNCG tetraloop element. RNA hairpin loops that conform to this consensus are known to be structurally stabilizing to adjoining base-paired regions [76, 77]. Nearly all representatives are located in the same location as the 5´ UTRs of genes that produce a protein of unknown function similar to the mammalian SART-1 (squamous cell carcinoma antigen recognized by T cells) protein . SART-1 is similar to the yeast Snu66 spliceosomal protein. One SART-1 RNA representative resides in the location of an annotated intron for the SART-1 gene. In mammals, the SART-1 gene encodes split ORFs that are possibly translated by a mechanism of -1 frameshifting.
Importantly, the SART-1 RNA motif is predicted to be formed by antisense transcripts of the associated gene. This hypothesis is supported by the fact that numerous G-U wobble pairs in the consensus model would otherwise become A-C mismatches if the RNA structure were formed by the sense transcript, which is unlikely. Furthermore, the commonly-occurring UNCG tetraloop structure in P1 would become a CGNA tetraloop in the sense direction, which is not known to confer the same structural benefit as a UNCG tetraloop.
AU-rich hairpin motif
Only 12 examples of the AU-rich hairpin motif (Fig. 7c) were identified from three fungal species, and eight of these representatives were from a single species, Rhizoctonia solani. Interestingly, two additional examples with considerable sequence and structural similarity were found in the bacterial genome of Orientia tsutsugamushi. These two bacterial examples were originally annotated as members of the mraW class of putative ncRNAs as reported previously . However, only a small portion of the consensus sequence for mraW motif RNAs is similar to the consensus sequence derived from the fungal examples of the AU-rich hairpin motif. Therefore it is not certain whether these similarities are biologically relevant or coincidental.
The function of the bacterial mraW motif RNAs remains unknown, and unfortunately the gene associations for the fungal examples are highly variable. As a result there are no compelling clues regarding the possible biological roles of the fungal motif. If additional examples are found in the future, perhaps genomic location data might provide the insight necessary to better formulate testable hypotheses.
Atypical snoRNA motif
The atypical snoRNA motif (Fig. 7d) is represented by 81 examples from 52 fungal species. Most of the representatives are located in introns, with the exception of seven that appear to reside apart from introns. For example, in the yeast Schizosaccharomyces pombe (NC_003424.3), the motif resides in the intron of the cpc2 gene from nucleotides 2440331 to 2440419. In most instances, this RNA motif is associated with genes encoding guanine nucleotide-binding proteins (G proteins), and more rarely are located near several genes for proteins of unknown function or for sterol 24-C-methyltransferase. G proteins are important signal transduction components whose functions are well established in a diversity of signaling pathways [79, 80], although the functions of the specific G proteins associated with this ncRNA motif are unknown.
These fungal ncRNA candidates exhibit considerable similarity to certain snoRNAs, which guide chemical modifications on other RNAs, including ribosomal RNAs (rRNAs), transfer RNA (tRNAs), and small nuclear RNAs (snRNAs) . One type, called box C/D snoRNAs, have a conserved C box (RUGAUGA) and D box (CUGA). Additional conserved regions called the Cˊ box and Dˊ box mimic the sequence of the C box and D box, respectively. The atypical snoRNA motif examples we identified carry two regions that closely approximate the C box (AUGAUGY) and D box (CUGA), although the apparent Cˊ box (AUGAGAC) and Dˊ box (CAGA) consensus sequences correspond to the consensus snoRNA sequences more poorly.
RT-PCR was used to evaluate the production of the A. nidulans representative of this RNA class. Only the RT-PCR product corresponding to the spliced RNA was observed (data not shown), indicating that the intron carrying the atypical snoRNA is efficiently removed from the original transcript. This result might indicate that the intron is always removed, rather than undergoing regulation by the structured element in the intron. If true, then this atypical snoRNA might have a function similar to other snoRNAs. Consistent with this hypothesis is the fact that the highly-conserved nucleotides in the guide sequence region of the atypical snoRNA are complementary to regions within 5.8S rRNA, 18S rRNA, and 28S rRNA, suggesting that members of this ncRNA class direct modifications to these rRNAs.
Group I ribozymes
Bioinformatics searches that are guided by consensus sequences and structure models for known RNA classes are likely to miss many representatives, particularly those that can vary considerably from the consensus model. Searches that can uncover conserved sequences and structures without relying on pre-existing consensus models can reveal distal variants of known ncRNA classes. As noted above, our bioinformatics search strategy already has revealed numerous additional representatives of the HDV class of self-cleaving ribozymes (Fig. 2) and members of an atypical snoRNA (Fig. 7d). Similarly, we have identified 208 examples of what appear to be previously unannotated group I self-splicing ribozymes  that are present in 114 fungal species (Fig. 7e).
Numerous additional examples of this motif were uncovered in bacteria by searching for sequences conforming to the resulting consensus model based on these fungal representatives (unpublished observations). Most of these newly-found representatives carry readily recognizable structural elements of group I ribozymes, including stems P3 through P7, and the conserved guanosine binding site. Since there is considerable sequence and structural variability at both the 5ˊ and 3ˊ termini of group I ribozymes, we did not further examine each representative to determine if they also carry stems P1, P2, P9, and P10 that are typical of this ribozyme class. However, these RNAs are most likely group I ribozymes that have previously escaped annotation in the genomes we analyzed.
We identified many distinct types of candidate structured RNA domains in close association with fungal genes coding for ribosomal proteins. The rps0 motif is one such representative (Fig. 8a) that included 41 representatives from 25 fungal species. Ribosomal protein genes commonly use a feedback auto-regulation mechanism to regulate their expression in bacteria [83,84,85]. These systems involve the binding of the ribosomal protein to a special RNA structure commonly located in the 5ˊ UTR of its corresponding mRNA. For example, E. coli uses at least 12 distinct RNA structures to regulate the expression of numerous ribosomal proteins . This general mechanism is also used by some eukaryotic species for regulating pre-mRNA splicing  or translation . In S. cerevisiae, ribosomal protein L32 (RPL32) binds to a structured RNA formed by intronic and exonic sequences within its own mRNA to cause alternative splicing. Moreover, the spliced product can also fold into a very similar structure that still binds the RPL32 protein and inhibits translation [89, 90].
The rps0 motif identified in the current study consists of two hairpins, one of which carries two conserved 5´-GGGGAAAG sequence elements partly located on each side of an internal loop. All the representatives of this motif are located in the 3´ UTRs of genes encoding 40S ribosomal protein subunit S0 (RPS0). RPS0 proteins are required for the processing of the precursor of 18S rRNAs and the formation of active 40S ribosomal subunits . Given the apparent symmetry of the rps0 RNA sequence and its location adjacent to the rpo0 gene, it seems possible that the motif might bind two or more RPS0 proteins to regulate expression of this ribosomal protein factor.
The rps2 motif (Fig. 8b) is represented by 13 examples from nine fungal species belonging to four subclasses including Pichia, Candida, Lodderomyces and Clavispora. All examples are located in the 5′ UTRs of the genes encoding 40S ribosomal protein S2. Although there are few examples, there are four predicted base-pairs distributed between P1 and P2 that covary in a manner consistent with the predicted secondary structure.
There are 20 representatives of the rps20 motif (Fig. 8c) derived from 13 species distributed among the Saccharomyces, Candida, Lachancea, Kazachstania and Tetrapisispora genera. All the representatives are located in the 3′ UTRs of the genes encoding 40S ribosomal protein S20. Again, given the limited number and distribution of members of this motif, we observe only two base-pairs with evidence of covariation in the proposed P3 stem. The other two predicted base-paired regions are formed by highly-conserved sequences and therefore further evidence of their formation is currently lacking.
The rpl7-l8-s3 motif (Fig. 8d) forms two long hairpins with internal loops. A total of 29 representatives from 25 species are located mostly in the introns of genes encoding 60S ribosomal proteins L7 and L8, as well as the 40S ribosomal protein S3. L7 and L8 proteins are required for the processing of 60S pre-rRNA, which includes the peptide bond formation center of 28S rRNA. The motif includes a well-conserved purine-rich region located between two internal loops of P2, which is similar to a previously reported ribosomal protein binding site . However, this previously reported motif in S. cerevisiae, which is located in the pre-mRNA for ribosomal protein L32, is proposed to be bound by the L32 protein as an unpaired bulge. In contrast, the conserved purine-rich region of the rpl7-l8-s3 motif is predicted to be part of an extended base-paired structure that includes evidence for some nucleotide covariation, implying that this region might not exist as unpaired RNA.
The rpl7 motif (Fig. 8e) is represented by 22 examples from 17 species. It is a small motif of approximately 30 nucleotides containing some conserved nucleotides in and around a central loop. All examples of this RNA motif are located in the introns of 60S ribosomal protein L7 mRNAs. In the absence of L7, other ribosomal proteins such as L6, L14, L20, and L33 are greatly diminished , and so this motif might be important for the coordination of ribosomal protein production.
We also used RT-PCR to test whether splicing occurs. If so, the motif might participate in regulating splicing events, as has been observed for an RNA structure associated with S. cerevisiae ribosomal protein L32 (RPL32) . RT-PCR analysis (Fig. 8e) reveals that the L7 pre-mRNA that contains the rpl7 motif produces at least one major splicing product. This result, along with the location of the rpl7 motif in introns, is consistent with the hypothesis that the rpl7 motif might regulate pre-mRNA splicing, perhaps by directly binding to the L7 ribosomal protein.
The rpl30 motif (Fig. 8f) is represented by 47 examples from 26 species. This motif is comprised of a long 5ˊ region with little evidence for structure formation, followed by a region that appears to form a large hairpin structure with a well-conserved purine-rich internal loop. Most representatives are located in the 5´ UTRs of genes encoding 60S ribosomal protein subunit L30. Notably, expression of S. cerevisiae ribosomal protein L32 is autogenously regulated by a mechanism wherein the protein binds to a purine-rich internal loop . This precedence strengthens the working hypothesis that the rpl30 motif is a regulatory RNA structure that serves as a binding partner for ribosomal protein L30 for expression autoregulation.
Comparative genomics is a powerful approach to evaluate the degree of conservation between even distantly-related genomic DNA sequences and has been widely applied to compare genes that encode proteins . We have used such approaches to discover numerous novel RNA motifs in bacteria [29,30,31,32]. However, using comparative genomics to discover structured ncRNAs in eukaryotes has not been commonly pursued, particularly among evolutionarily distant species [20, 95, 96]. Technical advances in analyzing numerous large genomes, such as improvements in computer processing speeds and more effective computational pipelines [28, 39], have created new opportunities to efficiently discover novel eukaryotic ncRNAs. Another advance is the availability of many sequenced eukaryotic genomes, which overcomes the problems otherwise caused by examining only a few phylogenetically similar eukaryotes.
To identify novel fungal ncRNAs, we first extracted IGRs of fungal genomes using available genome annotations. The elimination of coding regions greatly diminishes the computational challenge, and yet enriches the sequence data searched for candidates because known structured RNAs such as ribozymes and riboswitches almost never overlap with large ORFs. Comparative sequence analyses on this select sequence dataset using our pipeline reveal both conserved sequence and secondary structure features that define each of the novel RNA motifs. The list of representatives of promising ncRNA candidates was further expanded by additional homology searches, which iteratively refine the consensus models and allow for additional representatives to be discovered .
Through these efforts, we have identified 15 classes of structured ncRNA motifs in fungal genomes, including variants of HDV self-cleaving ribozymes, non-typical snoRNA candidates, the SDC motif, two antisense structured RNAs (ies6 and SART-1), a motif including a uORF that might regulate a gene that is critical for cell survival, a motif involved in Woronin body formation, and six motifs likely involved in regulating ribosomal protein expression. Although the functions of some motifs have been analyzed by using biochemical or genetic methods, many additional experiments will be needed to precisely define the roles and mechanisms for these ncRNA candidates.
Already, it appears that the functional roles of these ncRNAs are diverse. For example, representatives of the variant HDV ribozymes were found to efficiently self-cleave (Fig. 2) by using the same mechanism for phosphoester transfer as other self-cleaving ribozymes. The SDC motif appears to function as a translation repressor that likely coordinates with a uORF in the same 5´ UTR to regulate SAM decarboxylase synthesis (Fig. 4). It appears that the amd motif also encodes a conserved oligopeptide wherein this uORF overlaps with a domain of extensive RNA sequence and structure conservation. Possibly, uORF expression regulates the formation of amd motif RNA structure bring about the desired biological effect. Finally, RT-PCR results indicate that the function of the ies6 motif might be performed as part of an antisense transcript (Fig. 5), and that the rpl7 motif is present in a pre-mRNA that undergoes splicing (Fig. 8).
These findings demonstrate how the analysis of multiple fungal genomes can be used to reveal the existence of numerous structured ncRNA candidates. Similar approaches could be used to reveal the existence of additional structured ncRNAs in other eukaryotic systems whose biochemical and biological functions are highly diverse.
Identification of candidate RNA motifs
Fungal genome sequences including annotations were initially downloaded from the RefSeq Release 29 obtained from NCBI  and The Fungal Genomics group at the Broad Institute (http://www.broadinstitute.org/science/projects/fungal-genome-initiative/current-fgi-sequence-projects). The fungal sequences available as of November 2013 (RefSeq Release 62) were incorporated into this study (see also Additional file 17: Table S2).
The procedure to discover structured ncRNAs is outlined in Fig. 1. Briefly, IGRs from all genomes and introns from each annotated gene were extracted. The sequences also were selected by considering alternative splicing information according to the fungal genome annotations. IGRs and introns were compared among species of fungi by NCBI BLAST (version 2.2.23)  using parameters –m8 –W7 –e 1e–5. Similar sequences were automatically grouped by BLAST and a PERL script to form different clusters. Each cluster was filtered using BLAST against the Rfam database to remove any known RNAs with parameters –m8 –W10 –e 1e–10. The resulting new candidate clusters were further filtered by RNAcode  to remove clusters whose score is greater than 15, indicating that they are potentially protein coding. The remaining clusters were fed into CMfinder  to predict RNA structures based on covarying mutations among RNA sequences. Clusters with little covariation or with very few examples (less than three sequences in a cluster) were discarded. Finally, Infernal  was used to search for more homologs according to both sequence and structure similarity and to refine the structure of the RNA if more examples were discovered. During the Infernal search, the information about the RNA motif location and associated upstream and downstream genes were also inspected.
Assessing the novelty of motifs
To determine whether the predicted structured RNA motifs were reported previously, we took out every sequence of each alignment and compared them to sequences in the Rfam database using default parameters . Novel RNA candidates are those not found in Rfam, not highly homologous to RNAs in Rfam, or highly homologous to any known RNAs in other known databases such as NCBI gene bank databases. Notably, 11 sequences of putative HDV-like ribozymes were reported previously  and some of them overlapped with our collection, but the rest of the 223 HDV variant representatives were newly identified.
Ribozyme self-cleavage assays
Self-cleaving ribozyme assays were conducted using reaction conditions similar to those described previously . Synthetic RNA substrates were purchased (Sigma-Aldrich), whereas the enzyme RNA strands were prepared by in vitro transcription . Substrate RNAs were dephosphorylated, 5´-radiolabeled by using γ-32P ATP and T4 polynucleotide kinase (New England Biolabs), and purified by denaturing (8 M urea) 20% polyacrylamide gel electrophoresis (PAGE). 5 nM radiolabeled substrate and 100 nM enzyme RNA strands were combined in a reaction mixture containing 30 mM Tris-HCl (pH 7.5 at 23°C), 100 mM KCl, and 20 mM MgCl2, and incubated at 23°C for 30 min. The reaction was stopped by adding an equal volume of stop solution containing 90% formamide, 50 mM EDTA, 0.05% xylene cyanol and 0.05% bromophenol blue. Partial alkaline and RNase T1 digests of RNA samples were prepared for marker lanes T1 and –OH as described previously . The reaction products were separated by denaturing 20% PAGE and imaged/quantitated by a PhosphorImager (GE Healthcare).
Mass spectrometry analysis of cleavage products
Twenty pmol each of the HDV substrate and enzyme RNAs were incubated as described above in a 20 μL reaction for one hour. The reaction products were assayed by monoisotopic (exact mass) spectrometry (Novatia LLC).
Plasmids and strains
Plasmid pLL07-2-1 was constructed as described previously . To create an in-frame fusion of the luciferase reporter gene to the NCU01083 ORF, a DNA fragment of 967 base-pairs (including the first 62 base-pairs of the main ORF) was amplified by primers sdc-F and sdc-R (Additional file 18: Table S3), which include EcoRI and XbaI restriction sites, respectively. The fragment was sub-cloned into a pCR2.1-TOPO vector (Life Technologies). After confirmation by sequencing (Keck Foundation Biotechnology Resource Center at Yale University), the fragment was digested, purified by agarose gel electrophoresis, and inserted into the plasmid pLL07-2-1 at the EcoRI and XbaI sites. For mutant constructs, site-directed mutagenesis (Stratagene) and two-step PCR were used to make mutations at the appropriate sites.
N. crassa transformation and luciferase assay
Reverse transcriptase polymerase chain reaction (RT-PCR)
RT-PCR was carried out as described previously . Primers that are complementary to exon sequences located immediately upstream or downstream of the intronic motifs were used to detect alternative splicing products. For antisense motif ies6 two pairs of primers were designed to detect antisense and sense transcripts independently. Sequences of primers and predicted sizes of RT-PCR products are listed in Additional file 18: Table S3. During the RT-PCR process, a negative control, in which the RT reaction lacked reverse transcriptase, was included to ensure there was no DNA contamination.
Open reading frame
Polyacrylamide gel electrophoresis
Reverse transcription – polymerase chain reaction
small nucleolar RNAs
small nuclear RNAs
Hawksworth DL. The fungal dimension of biodiversity - magnitude, significance, and conservation. Mycol. Res. 1991;95(6):641–55.
Hawksworth DL. The magnitude of fungal diversity: the 1.5 million species estimate revisited. Mycol Res. 2001;105(12):1422–32.
Galagan JE, Henn MR, Ma LJ, Cuomo CA, Birren B. Genomics of the fungal kingdom: insights into eukaryotic biology. Genome Res. 2005;15(12):1620–31.
Goffeau A, Barrell BG, Bussey H, Davis RW, Dujon B, Feldmann H, et al. Life with 6000 genes. Science. 1996;274(5287):546–67.
Wood V, Gwilliam R, Rajandream MA, Lyne M, Lyne R, Stewart A, et al. The genome sequence of Schizosaccharomyces pombe. Nature. 2002;415(6874):871–80.
Galagan JE, Calvo SE, Borkovich KA, Selker EU, Read ND, Jaffe D, et al. The genome sequence of the filamentous fungus Neurospora crassa. Nature. 2003;422(6934):859–68.
Hedeler C, Wong HM, Cornell MJ, Alam I, Soanes DM, Rattray M, et al. e-Fungi: a data resource for comparative analysis of fungal genomes. BMC Genomics. 2007;8:426.
Breaker RR. Riboswitches and the RNA world. Cold Spring Harb. Perspect. Biol. 2012;4(2):pii:a003566.
Lilley DM. Structure, folding and mechanisms of ribozymes. Curr. Opin. Struct. Biol. 2005;15(3):313–23.
Argaman L, Hershberg R, Vogel J, Bejerano G, Wagner EG, Margalit H, et al. Novel small RNA-encoding genes in the intergenic regions of Escherichia coli. Curr. Biol. 2001;11(12):941–50.
Wassarman KM, Repoila F, Rosenow C, Storz G, Gottesman S. Identification of novel small RNAs using comparative genomics and microarrays. Genes Dev. 2001;15(13):1637–51.
Olivas WM, Muhlrad D, Parker R. Analysis of the yeast genome: identification of new non-coding and small ORF-containing RNAs. Nucleic Acids Res. 1997;25(22):4619–25.
Klein RJ, Misulovin Z, Eddy SR. Noncoding RNA genes identified in AT-rich hyperthermophiles. Proc. Natl. Acad. Sci. USA. 2002;99(11):7542–7.
Carter RJ, Dubchak I, Holbrook SRA. computational approach to identify genes for functional RNAs in genomic sequences. Nucleic Acids Res. 2001;29(19):3928–38.
Rivas E, Eddy SR. Noncoding RNA gene detection using comparative sequence analysis. BMC Bioinformatics. 2001;2:8.
Rivas E, Klein RJ, Jones TA, Eddy SR. Computational identification of noncoding RNAs in E. coli by comparative genomics. Curr. Biol. 2001;11(17):1369–73.
Achar A, Sætrom PRNA. motif discovery: a computational overview. Biol. Direct. 2015;10:61.
Daub J, Eberhardt RY, Tate JG, Burge SW. Rfam: annotating families of non-coding RNA sequences. Methods Mol. Biol. 2015;1269:349–63.
Lowe TM, Eddy SRA. computational screen for methylation guide snoRNAs in yeast. Science. 1999;283(5405):1168–71.
McCutcheon JP, Eddy SR. Computational identification of non-coding RNAs in Saccharomyces cerevisiae by comparative genomics. Nucleic Acids Res. 2003;31(14):4119–28.
Steigele S, Huber W, Stocsits C, Stadler PF, Nieselt K. Comparative analysis of structured RNAs in S. cerevisiae indicates a multitude of different functions. BMC Biol. 2007;5:25.
Schattner P, Decatur WA, Davis CA, Ares M, Jr, Fournier MJ, Lowe TM. Genome-wide searching for pseudouridylation guide snoRNAs: analysis of the Saccharomyces cerevisiae genome. Nucleic Acids Res. 2004;32(14):4281-4296.
Liu N, Xiao ZD, CH Y, Shao P, Liang YT, Guan DG, et al. SnoRNAs from the filamentous fungus Neurospora crassa: structural, functional and evolutionary insights. BMC Genomics. 2009;10:15.
Jochl C, Rederstorff M, Hertel J, Stadler PF, Hofacker IL, Schrettl M, et al. Small ncRNA transcriptome analysis from Aspergillus fumigatus suggests a novel mechanism for regulation of protein synthesis. Nucleic Acids Res. 2008;36(8):2677–89.
Liu T, Ren X, Xiao T, Yang J, Xu X, Dong J, et al. Identification and characterisation of non-coding small RNAs in the pathogenic filamentous fungus Trichophyton rubrum. BMC Genomics. 2013;14:931.
Eddy SR, Durbin RRNA. sequence analysis using covariance models. Nucleic Acids Res. 1994;22(11):2079–88.
Washietl S, Hofacker IL, Stadler PF. Fast and reliable prediction of noncoding RNAs. Proc. Natl. Acad. Sci. USA. 2005;102(7):2454–9.
Yao Z, Weinberg Z, Ruzzo WL. CMfinder-a covariance model based RNA motif finding algorithm. Bioinformatics. 2006;22(4):445–52.
Weinberg Z, Barrick JE, Yao Z, Roth A, Kim JN, Gore J, et al. Identification of 22 candidate structured RNAs in bacteria using the CMfinder comparative genomics pipeline. Nucleic Acids Res. 2007;35(14):4809–19.
Ames TD, Rodionov DA, Weinberg Z, Breaker RR. A eubacterial riboswitch class that senses the coenzyme tetrahydrofolate. Chem. Biol. 2010;17(7):681–5.
Weinberg Z, Wang JX, Bogue J, Yang J, Corbino K, Moy RH, et al. Comparative genomics reveals 104 candidate structured RNAs from bacteria, archaea, and their metagenomes. Genome Biol. 2010;11:R31.
Roth A, Weinberg Z, Chen AG, Kim PB, Ames TD, Breaker RRA. widespread self-cleaving ribozyme class is revealed by bioinformatics. Nat. Chem. Biol. 2014;10(1):56–60.
Bailey TL. Discovering novel sequence motifs with MEME. Curr. Protoc. Bioinformatics. Chapter 2:Unit. 2002;2:4.
Storz G. An expanding universe of noncoding RNAs. Science. 2002;296(5571):1260–3.
Pruitt KD, Tatusova T, Maglott DR. NCBI reference sequences (RefSeq): a curated non-redundant sequence database of genomes, transcripts and proteins. Nucleic Acids Res. 2007;35(database issue):D61–5.
Altschul SF, Madden TL, Schaffer AA, Zhang JH, Zhang Z, Miller W, et al. Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 1997;25(17):3389–402.
Gardner PP, Daub J, Tate JG, Nawrocki EP, Kolbe DL, Lindgreen S, et al. Rfam: updates to the RNA families database. Nucleic Acids Res. 2009;37(Database issue):D136–40.
Washietl S, Findeiss S, Muller SA, Kalkhof S, von Bergen M, Hofacker IL, et al. RNAcode: Robust discrimination of coding and noncoding regions in comparative sequence data. RNA. 2011;17(4):578–94.
Nawrocki EP, Eddy SR. Infernal 1.1: 100-fold faster RNA homology searches. Bioinformatics. 2013;29(22):2933–5.
Been MD. HDV ribozymes. Curr. Top. Microbiol. Immunol. 2006;307:47–65.
Riccitelli N, Lupták AHDV. family of self-cleaving ribozymes. Prog. Mol. Biol. Transl. Sci. 2013;120:123–71.
Sharmeen L, Kuo MYP, Dinter-Gottieb G, Taylor J, Antigenomic RNA. of human hepatitis delta virus can undergo self-cleavage. J. Virol. 1988;62(8):2674–9.
Salehi-Ashtiani K, Lupták A, Litovchick A, Szostak JWA. genomewide search for ribozymes reveals an HDV-like sequence in the human CPEB3 gene. Science. 2006;313(5794):1788–92.
Webb C-HT, Riccitelli NJ, Ruminski DJ, Lupták A. Widespread occurrence of self-cleaving ribozymes. Science. 2009;326(5955):953.
Webb C-HT, Lupták A. HDV-like self-cleaving ribozymes. RNA Biol. 2011;8:719–27.
Wadkins TS, Perrotta AT, Ferre-D'Amare AR, Doudna JA, Been MD. A nested double pseudoknot is required for self-cleavage activity of both the genomic and antigenomic hepatitis delta virus ribozymes. RNA 1999;5(6):720-727.
Fox GE, Woese CR. 5S rRNA secondary structure. Nature. 1975;256(5517):505–7.
Correll CC, Freeborn B, Moore PB, Steitz TA. Metals, motifs, and recognition in the crystal structure of a 5S rRNA domain. Cell. 1997;91(5):705–12.
Dallas A, Moore PB. The loop E-loop D region of Escherichia coli 5S rRNA: the solution structure reveals an unusual loop that may be important for binding ribosomal proteins. Curr. Biol. 1997;5(12):1639–53.
Weinberg Z, Kim PB, Chen TH, Li S, Harris KA, Lünse CE, et al. New classes of self-cleaving ribozymes revealed by comparative genomics analysis. Nat. Chem. Biol. 2015;11(8):606–10.
Lovett PS, Rogers EJ. Ribosome regulation by the nascent peptide. Microbiol. Rev. 1996;60(2):366–85.
Vilela C, McCarthy JE. Regulation of fungal gene expression via short open reading frames in the mRNA 5ˊ untranslated region. Mol. Microbiol. 2003;49(4):859–67.
Luo ZL, Sachs MS. Role of an upstream open reading frame in mediating arginine-specific translational control in Neurospora crassa. J. Bacteriol. 1996;178(8):2172–7.
Wang Z, Sachs MS. Ribosome stalling is responsible for arginine-specific translational attenuation in Neurospora crassa. Mol. Cell Biol. 1997;17(9):4904–13.
WW H, Gong H, Pua EC. The pivotal roles of the plant S-adenosylmethionine decarboxylase 5ˊ untranslated leader sequence in regulation of gene expression at the transcriptional and posttranscriptional levels. Plant Physiol. 2005;138(1):276–86.
Pegg AE, Xiong H, Feith DJ, Shantz LM. S-adenosylmethionine decarboxylase: structure, function and regulation by polyamines. Biochem. Soc. Trans. 1998;26:580–6.
Shantz LM, Viswanath R, Pegg AE. Role of the 5'-untranslated region of mRNA in the synthesis of S-adenosylmethionine decarboxylase and its regulation by spermine. Biochem. J. 1994;302(Pt 3):765–72.
Ivanov IP, Atkins JF, Michael AJA. profusion of upstream open reading frame mechanisms in polyamine-responsive translational regulation. Nucleic Acids Res. 2010;38(2):353–9.
Soukup GA, Breaker RR. Relationship between internucleotide linkage geometry and the stability of RNA. RNA 1999;5(10):1308-1325.
Regulski EE, Breaker RR. In-line probing analysis of riboswitches. Methods Mol Biol. 2008;419:53–67.
Cheah MT, Wachter A, Sudarsan N, Breaker RR. Control of alternative RNA splicing and gene expression by eukaryotic riboswitches. Nature. 2007;447(7143):497–500.
Colot HV, Park G, Turner GE, Ringelberg C, Crew CM, Litvinkova L, et al. A high-throughput gene knockout procedure for Neurospora reveals functions for multiple transcription factors. Proc. Natl. Acad. Sci. USA. 2006;103(27):10352–7.
Woese CR, Winkler S, Gutell RR. Architecture of ribosomal RNA: constraints on the sequence of “tetra-loops”. Proc. Natl. Acad. Sci. USA. 1990;87(21):8467–71.
Costa M, Michel F. Rules for RNA recognition of GNRA tetraloops deduced by in vitro selection: comparison with in vivo evolution. EMBO J. 1997;16(11):3289–302.
Fiore JL, Nesbitt DJ, An RNA. folding motif: GNRA tetraloop-receptor interactions. Q. Rev. Biophys. 2013;46(3):223–64.
Chambers AL, Ormerod G, Durley SC, Sing TL, Brown GW, Kent NA, et al. The INO80 chromatin remodeling complex prevents polyploidy and maintains normal chromatin structure at centromeres. Genes Dev. 2012;26(23):2590–603.
Morrison AJ, Highland J, Krogan NJ, Arbel-Eden A, Greenblatt JF, Haber JE, et al. INO80 and gamma-H2AX interaction links ATP-dependent chromatin remodeling to DNA damage repair. Cell. 2004;119(6):767–75.
EY Y, Steinberg-Neifach O, Dandjinou AT, Kang F, Morrison AJ, Shen X, et al. Regulation of telomere structure and functions by subunits of the INO80 chromatin remodeling complex. Mol Cell Biol. 2007;27(16):5639–49.
Faghihi MA, Wahlestedt C. Regulatory roles of natural antisense transcripts. Nat. Rev. Mol. Cell Biol. 2009;10(9):637–43.
Katayama S, Tomaru Y, Kasukawa T, Waki K, Nakanishi M, Nakamura M, et al. Antisense transcription in the mammalian transcriptome. Science. 2005;309(5740):1564–6.
Donaldson ME, Saville BJ. Natural antisense transcripts in fungi. Mol. Microbiol. 2012;85(3):405–17.
Brantl S. Regulatory mechanisms employed by cis-encoded antisense RNAs. Curr. Opin. Microbiol. 2007;10(2):102–9.
Markham P, Collinge AJ. Woronin bodies of filamentous fungi. FEMS Microbiol. Lett. 1987;46(1):1–11.
Beck J, Ebel F. Characterization of the major Woronin body protein HexA of the human pathogenic mold Aspergillus fumigatus. Int. J. Med. Microbiol. 2013;303(2):90–7.
Yuan P, Jedd G, Kumaran D, Swaminathan S, Shio H, Hewitt D, et al. A HEX-1 crystal lattice required for Woronin body function in Neurospora crassa. Nat. Struct. Biol. 2003;10(4):264–70.
Varani G, Cheong C, Tinoco I Jr. Structure of an unusually stable RNA hairpin. Biochemistry. 1991;30(13):3280–9.
Molinaro M, Tinoco I Jr. Use of ultra stable UNCG tetraloop hairpins to fold RNA structures: thermodynamic and spectroscopic applications. Nucleic Acids Res. 1995;23(15):3056–63.
Shichijo S, Nakao M, Imai Y, Takasu H, Kawamoto M, et al. A gene encoding peptides of human squamous cell carcinoma recognized by cytotoxic T lymphocytes. J. Exp. Med. 1998;187(3):277–88.
Neer EJ, Clapham DE. Roles of G protein subunits in transmembrane signalling. Nature. 1988;333(6169):129–34.
Neves SR, Ram PT, Iyengar RG. protein pathways. Science. 2002;296(5573):1636–9.
Kiss T. Small nucleolar RNAs: an abundant group of noncoding RNAs with diverse cellular functions. Cell. 2002;109(2):145–8.
Hausner G, Hafez M, Edgell DR. Bacterial group I introns: mobile RNA catalysts. Mobile DNA. 2014;5(1):8.
Nomura M, Yates JL, Dean D, Post LE. Feedback regulation of ribosomal protein gene expression in Escherichia coli: structural homology of ribosomal RNA and ribosomal protein mRNA. Proc. Nat. Acad. Sci. USA. 1980;77(12):7084–8.
Zengel JM, Lindahl L. Diverse mechanisms for regulating ribosomal protein synthesis in Escherichia coli. Prog. Nucleic Acid Res. Mol. Biol. 1994;47:331–70.
Deiorio-Haggar K, Anthony J, Meyer MMRNA. structures regulating ribosomal protein biosynthesis in bacilli. RNA Biol. 2013;10(7):1180–4.
Fu Y, Deiorio-Haggar K, Anthony J, Meyer MM. Most RNAs regulating ribosomal protein biosynthesis in Escherichia coli are narrowly distributed to Gammaproteobacteria. Nucleic Acids Res. 2013;41(6):3491–503.
Ivanov AV, Malygin AA, Karpova GG. Eukaryotic ribosomal proteins: interactions with their own pre-mRNAs and their involvement in splicing regulation. Mol. Biol. 2006;40(4):570–8.
Amaldi F, Bozzoni I, Beccari E, Pierandrei-Amaldi P. Expression of ribosomal protein genes and regulation of ribosome biosynthesis in Xenopus development. Trends Biochem. Sci. 1989;14(5):175–8.
Vilardell J, Warner JR. Regulation of splicing at an intermediate step in the formation of the spliceosome. Genes Dev. 1994;8(2):211–20.
Vilardell J, Warner JR. Ribosomal protein L32 of Saccharomyces cerevisiae influences both the splicing of its own transcript and the processing of rRNA. Mol. Cell Biol. 1997;17(4):1959–65.
Ford CL, Randal-Whitis L, Ellis SR. Yeast proteins related to the p40/laminin receptor precursor are required for 20S ribosomal RNA processing and the maturation of 40S ribosomal subunits. Cancer Res. 1999;59(3):704–10.
Li H, Dalal S, Kohler J, Vilardell J, White SA. Characterization of the pre-mRNA binding site for yeast ribosomal protein L32: the importance of a purine-rich internal loop. J. Mol. Biol. 1995;250(4):447–59.
Jakovljevic J, Ohmayer U, Gamalinda M, Talkish J, Alexander L, Linnemann J, Milkereit P, Woolford JL. Ribosomal proteins L7 and L8 function in concert with six A3 assembly factors to propagate assembly of domains I and II of 25S rRNA in yeast 60S ribosomal subunits. RNA. 2012;18(10):1805–22.
Hardison RC. Comparative genomics. PLoS Biol. 2003;1(2):e58.
Rubin GM, Yandell MD, Wortman JR, Miklos GLG, Nelson CR, Hariharan IK, Fortini ME, Li PW, Apweiler R, Fleischmann W, et al. Comparative genomics of the eukaryotes. Science. 2000;287(5461):2204–15.
Will S, Yu M, Berger B. Structure-based whole-genome realignment reveals many novel noncoding RNAs. Genome Res. 2013;23(6):1018–27.
Roth A, Nahvi A, Lee M, Jona I, Breaker RR. Characteristics of the glmS ribozyme suggest only structural roles for divalent metal ions. RNA. 2006;12(4):607–19.
Li S, Breaker RR, Eukaryotic TPP. riboswitch regulation of alternative splicing involving long-distance base pairing. Nucleic Acids Res. 2013;41(5):3022–31.
Froehlich AC, Loros JJ, Dunlap JC. Rhythmic binding of a WHITE COLLAR-containing complex to the frequency promoter is inhibited by FREQUENCY. Proc. Natl. Acad. Sci. U.S.A. 2003;100(10):5914–9.
We thank Z. Weinberg, A. Roth, N. Sudarsan, K. Harris and other members of the Breaker laboratory for helpful discussions and recommended edits to the manuscript, Z. Weinberg for advice on computational searches, and N. Carriero and R. Bjornson for assistance with the Yale Life Sciences High Performance Computing Center. All authors have read and approved of the final version of this manuscript.
Funding was supplied by NIH grant RR19895 to the Yale Life Sciences High Performance Computing Center and NIH grant (GM022778 to R.R.B. R.R.B. is an Investigator with the Howard Hughes Medical Institute.
Availability of data and materials
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Discovered ncRNA Information. Properties of and species for each ncRNA discovered. (DOCX 24 kb)
HDV ribozyme sequences. sto file for HDV ribozyme sequences and their alignment. (STO 40 kb)
SDC motif sequences. sto file for SDC motif RNAs and their alignment. (STO 3 kb)
amd motif sequences. sto file for amd motif sequences and their alignment. (STO 4 kb)
ies6 motif sequences. sto file for ies6 motif sequences and their alignment. (STO 5 kb)
hexA motif sequences. sto file for hexA motif sequences and their alignment. (STO 2 kb)
SART-1 motif sequences. sto file for SART-1 motif sequences and their alignment. (STO 2 kb)
AU-rich motif sequences. sto file for AU-rich motif sequences and their alignment. (STO 1 kb)
Atypical snoRNA motif sequences. sto file for atypical snoRNA motif sequences and their alignment. (STO 15 kb)
Group I ribozyme sequences. sto file for group I ribozyme sequences and their alignment. (STO 59 kb)
rps0 motif sequences. sto file for rps0 motif sequences and their alignment. (STO 9 kb)
rps2 motif sequences. sto file for rps2 motif sequences and their alignment. (STO 1 kb)
rps20 motif sequences. sto file for rps20 motif sequences and their alignment. (STO 2 kb)
rps7-l8-s3 motif sequences. sto file for rps7-l8-s3 motif sequences and their alignment. (STO 5 kb)
rpl7 motif sequences. sto file for rpl7 motif sequences and their alignment. (STO 1 kb)
rpl30 motif sequences. sto file for rpl30 motif sequences and their alignment. (STO 12 kb)
Genome database lists. A listing of the genomic database files used for computational searches. (DOCX 13 kb)
PCR primers. DNA sequences of PCR primers used for the study. (DOCX 17 kb)
About this article
Cite this article
Li, S., Breaker, R.R. Identification of 15 candidate structured noncoding RNA motifs in fungi by comparative genomics. BMC Genomics 18, 785 (2017). https://doi.org/10.1186/s12864-017-4171-y