- Research article
- Open Access
MicroRNAs and essential components of the microRNA processing machinery are not encoded in the genome of the ctenophore Mnemiopsis leidyi
BMC Genomicsvolume 13, Article number: 714 (2012)
MicroRNAs play a vital role in the regulation of gene expression and have been identified in every animal with a sequenced genome examined thus far, except for the placozoan Trichoplax. The genomic repertoires of metazoan microRNAs have become increasingly endorsed as phylogenetic characters and drivers of biological complexity.
In this study, we report the first investigation of microRNAs in a species from the phylum Ctenophora. We use short RNA sequencing and the assembled genome of the lobate ctenophore Mnemiopsis leidyi to show that this species appears to lack any recognizable microRNAs, as well as the nuclear proteins Drosha and Pasha, which are critical to canonical microRNA biogenesis. This finding represents the first reported case of a metazoan lacking a Drosha protein.
Recent phylogenomic analyses suggest that Mnemiopsis may be the earliest branching metazoan lineage. If this is true, then the origins of canonical microRNA biogenesis and microRNA-mediated gene regulation may postdate the last common metazoan ancestor. Alternatively, canonical microRNA functionality may have been lost independently in the lineages leading to both Mnemiopsis and the placozoan Trichoplax, suggesting that microRNA functionality was not critical until much later in metazoan evolution.
MicroRNAs (miRNAs) are a class of small RNA molecules derived from transcribed mRNA hairpin structures and spliced introns [1–3] that play a key role in mRNA targeting, leading to the degradation or translational repression of the target transcript. The regulatory functions of miRNAs are essential to many key biological processes in metazoans, including development, cell growth and death, stem cell maintenance, hematopoiesis, and neurogenesis. Aberrations in miRNA regulation have been linked to blood disorders, oncogenesis, and other malignancies in humans . The hairpin structures in mRNA transcripts that give rise to primary microRNAs (pri-miRNAs) are not unique to miRNAs or metazoans; these hairpins can form much more frequently than functional pri-miRNAs [3, 5] and can arise from inverted duplications, transposable elements, and genomic repeats [3, 6, 7]. Metazoans, however, possess a unique complement of cellular machinery for processing and transporting mature miRNAs to their targets that has not been identified in any non-metazoan species to date [8–11]. It has been observed that once novel miRNAs emerge in a metazoan lineage, they are very rarely lost. Thus, miRNAs are thought to represent strong phylogenetic markers and, through their ability to fine-tune gene expression, appear to be major drivers of biological complexity [8, 12–14].
The canonical miRNA biogenesis pathway in metazoans is part of the larger RNA interference (RNAi) pathway, which includes the closely related siRNA pathway (Figure 1). The miRNA pathway is distinct from the ancestral siRNA pathway in that it is initiated by the cleavage of hairpin structures (i.e., pri-miRNAs) from mRNAs in the nucleus by the Drosha/Pasha complex (also known as the Microprocessor complex), producing precursor-miRNAs (i.e., pre-miRNAs) that can be exported into the cytosol via the Exportin-5—Ran-GTP complex. After being transported into the cytosol, miRNAs and siRNAs undergo the same processing and targeting steps, initiated by Dicer cleavage and loading into the RNA-induced silencing complex (RISC) with Argonaute . The siRNA pathway is an ancient biological defense mechanism used to ward off the integration of foreign nucleic acids, such as double stranded RNAs (dsRNAs) introduced by viruses, and is known to have existed in the oldest eukaryotes [7, 10]. Thus, the emergence of the metazoan canonical miRNA biogenesis pathway most likely coincided with the evolution of the Drosha/Pasha complex found only in metazoans [10, 11]. Functionally, the Drosha/Pasha complex enables cleavage of pri-miRNA hairpins that are subsequently exported out of the nucleus and processed by the pre-existing RNAi pathway.
Given the differences in molecular machinery, processing, and target recognition, miRNAs are thought to have evolved separately and exclusively in animals and plants [3, 7, 9, 16]. However, a number of recent studies have reported identification of miRNAs in unicellular eukaryotes, including several thought to be homologs of miRNAs specific to animal and plant lineages [17–29]. These studies imply that miRNAs evolved once, early in eukaryotic evolution. Nevertheless, a recent report  reexamined these studies and found that, of the cumulative 232 reported miRNAs, none of the putative plant or animal homologs met established criteria for miRNA annotation; they were, instead, likely traces of other small RNAs (e.g., siRNAs, rRNAs, or snoRNAs) that happened to fit the length spectrum of mature miRNA sequences. Additionally, only 28 of the putative novel miRNAs passed the annotation criteria, and those were restricted to green and brown algae. In light of this evidence, it appears most likely that miRNAs evolved independently in multiple eukaryotic lineages, with the metazoan pathway being dependent upon the Drosha/Pasha protein complex.
Here, we describe an in-depth characterization of both the miRNA biogenesis pathway proteins and genomic regions that may correspond to pri-miRNA loci in the recently sequenced genome of Mnemiopsis leidyi (http://research.nhgri.nih.gov/mnemiopsis/). Recent phylogenomic analyses suggest that Ctenophora may be the earliest branching metazoan lineage [31, 32], and genomic studies of a number of gene superclasses [33, 34] and signaling pathways  in Mnemiopsis are consistent with this theory. If ctenophores are, indeed, the earliest metazoan branch, examining the genome of Mnemiopsis provides us a rare opportunity to better understand the origin of miRNA processing in metazoans. Alternatively, if ctenophores branched later in evolution and Porifera is the most basal metazoan lineage , Mnemiopsis still provides a valuable model from which to study the early evolution of this important small RNA processing pathway. Putative miRNAs (and the pathway proteins involved in their canonical biogenesis) have been studied in other non-bilaterian metazoans, including Nematostella vectensis, Hydra magnipapillata, Trichoplax adhaerens, and Amphimedon queenslandica[9, 13, 37]. The complete processing pathway was identified in all cases except Trichoplax, which lacks a Pasha homolog and recognizable miRNAs [6, 9, 38]. However, the presence of Drosha, Pasha, and miRNAs in Amphimedon, a metazoan lineage that branched prior to Trichoplax, suggests that Trichoplax must have lost miRNA functionality .
Results and discussion
In order to understand the increasing complexity observed in the early evolution of animals, we have sequenced, annotated, and performed a preliminary analysis of the Mnemiopsis genome. During this process, we were able to map 99.4% of the 15,752 publicly available Mnemiopsis EST sequences to our genome assembly. These data are available through our Mnemiopsis Genome Project Web site (http://research.nhgri.nih.gov/mnemiopsis/). This Web site provides access to the assembled genome scaffolds, predicted protein models, transcriptome data, and EST data. The Web site also provides access to the Mnemiopsis Genome Browser, a BLAST utility, a gene-centric Wiki, protein domain annotations, and information on gene clusters mapped to human KEGG pathways via an intuitive and easy-to-use interface.
Through our examination of the Mnemiopsis genome and its predicted proteome, we were able to identify multiple RNAi pathway proteins necessary for miRNA and siRNA processing, including Dicer, Argonaute, Ran, and exportin-5, but the miRNA-specific biogenesis pathway proteins Drosha and Pasha are strikingly absent. To our knowledge, this is the first reported case of a metazoan genome lacking a Drosha homolog. Since Dicer and Drosha are both members of the ribonuclease III (RNase III) protein family (Figure 2), we focused our analysis on the RNase III protein domain to better characterize the Mnemiopsis Dicer protein and to yield insight into how, through the evolution of this protein family in the Metazoa, the canonical miRNA biogenesis pathway may have emerged.
Drosha and Dicer belong to subclasses 2 (Drosha) and 3 (Dicer) of the RNase III protein family . Both proteins are characterized by tandem RNase III domains that cleave dsRNA to a specific length, often producing cleavage products with a two-nucleotide 3′ overhang. However, distinct differences have been observed in the dsRNA-binding specificity and cellular localization of these two RNase III subclasses . Class 3 RNase III enzymes have a PAZ domain that recognizes dsRNA ends with the distinctive two-nucleotide 3′ overhang indicative of prior RNase III cleavage. Class 2 RNase III enzymes do not appear to contain a domain with specific affinity for dsRNA and, instead, rely on complex formation in the nucleus with a co-factor (Pasha, or DGCR8 in vertebrates) that recognizes the ssRNA-dsRNA junctions characteristic of pri-miRNA hairpins . RNase III class 3 Dicer-like proteins that lack a PAZ domain (and have a domain structure more similar to Drosha) have been identified in non-metazoans but function as part of an unrelated pathway ; they have also been identified in early branching metazoans, but their function has not been confirmed experimentally . Since deletion of the PAZ domain in a functional Dicer has been shown to produce an RNase III enzyme without target specificity , there are likely functional binding domains other than PAZ within the RNase III class 3 subfamily.
To determine which class(es) of RNase III enzymes the Mnemiopsis Dicer protein is most closely related to, we performed a phylogenetic analysis on the RNase III domains of early-branching metazoan Dicer and Drosha proteins. We used HMMER  to search available non-bilaterian animal protein sequences (i.e., Mnemiopsis, Nematostella, Hydra, Trichoplax, and Amphimedon) to identify all candidate class 2 or class 3 RNase III proteins containing tandem RNase III domains. Our search yielded only one Dicer protein in Mnemiopsis and numbers of proteins consistent with other reports on the early-branching Metazoa [9, 43]. We included a sample of bilaterian Dicer and Drosha sequences in our analysis to ensure each protein class was monophyletic across the Metazoa. We separated the RNase IIIa and RNase IIIb domains of each protein (Figure 2), aligned the domains, trimmed the poorly conserved and flanking regions, and used the resulting alignment as the basis for further phylogenetic analysis (see Additional file 1: Dataset 1a-b).
The tree generated from this alignment (Figure 3a) contains separate clades for each RNase III domain subgroup, confirming the characterization of the Mnemiopsis RNase III protein as a Dicer protein. Importantly, the topology unites the Drosha RNase IIIa and RNase IIIb domains with the respective Dicer RNase III domains. Given that RNase III class 2 (Drosha) proteins are restricted to the Metazoa [10, 11], whereas RNase III class 3 (Dicer) proteins are found in the RNAi pathways of ancestral eukaryotes [7, 10, 43], this topology suggests that Drosha evolved from Dicer via a duplication event early in the evolution of the Metazoa, roughly coinciding with the emergence of miRNA functionality (Figure 3b). This observation contradicts the less parsimonious argument that these double RNase III domain-containing enzymes evolved independently from separate eubacterial RNase III domains  (Additional file 2: Figure S1).
It is possible that Mnemiopsis utilizes alternative methods for producing miRNAs for transcriptional regulation. Therefore, we searched for miRNAs using data from short RNA sequencing runs on two Mnemiopsis samples. We were unable to identify any known metazoan miRNAs that mapped to the Mnemiopsis genome. While we were able to predict several novel miRNA candidates using two methods, no predictions were reproducible across all samples and methods. In addition, even the highest-scoring predictions exhibited atypical read mapping signatures. Thus, we have classified all of these predictions as false positives, as they do not appear to be processed by the canonical miRNA machinery (see Methods).
Some spliced introns can correctly fold into pre-miRNAs, called mirtrons, independent of cleavage by Drosha and Pasha [1, 2, 6] (Figure 1). However, within the Mnemiopsis genome, only a handful of introns have predicted secondary structures suggestive of mirtron-coding potential, and none of these have read mapping signatures to indicate that they are functional mirtrons. The presence of exportin-5 and downstream RNAi pathway proteins Dicer and Argonaute in Mnemiopsis could indicate the existence of an alternative mechanism for miRNA production that predates the canonical miRNA pathway. The lack of recognizable miRNAs in our small RNA sequences, however, suggests that this scenario is unlikely. Recently, cases of functional exogenous miRNAs acquired via ingestion were identified in animals , suggesting a possible dietary mechanism by which Mnemiopsis could utilize miRNA regulatory functions in the absence of a functional endogenous canonical pathway. However, the mechanism for exogenous miRNA activity remains poorly understood.
It has been hypothesized that mirtrons may have predated the Drosha/Pasha-mediated pathway, based on the observation that the mechanistic requirements for their evolution may have been fairly simple [1, 2]. The identification of mirtrons in rice [3, 45] and the presence of the necessary machinery in Mnemiopsis (described above) are consistent with this hypothesis. However, given the absence of functional mirtrons in Mnemiopsis, it appears more likely that miRNA functionality evolved alongside the Drosha/Pasha-mediated pathway, independently of the mirtron pathway. Discerning the point in evolutionary time in which mirtrons became functional will require a thorough analysis of the genomes of additional species beyond nematodes, mammals, and avians [3, 45].
The implications of these results depend upon the phylogenetic position of Ctenophora. If ctenophores are the most basal metazoan clade, the most parsimonious explanation for our observations is that metazoan miRNA functionality originated after ctenophores diverged from the rest of animals (Figure 4a). Alternatively, if poriferans are the most basal metazoan clade, then Drosha, Pasha and canonical miRNA functionality must have been lost in the Mnemiopsis lineage (Figure 4b). If the latter were true, then canonical microRNAs and their machinery would have been independently lost in both Ctenophora and Placozoa. This, along with the large-scale losses of miRNAs described in acoelomorphs  and cnidarians , would contradict the premise that miRNAs are ultraconserved, canalized characters that are continuously added, but rarely lost – and, as such, would challenge their usefulness as phylogenetic markers [12, 13].
Our data supports a scenario in which the role of miRNAs in fine-tuning gene expression was not solidified until more recently in metazoan evolution and thus indicates that miRNA regulatory functions were, perhaps, non-essential during early metazoan diversification. Given this, the lack of recognizable miRNA functionality in Mnemiopsis supports a scenario with Ctenophora branching at the base of the Metazoa, prior to the emergence of miRNA functionality (Figure 4a). It may also indicate that a novel RNA-based regulatory pathway evolved either within the ctenophore lineage or as a precursor to the canonical miRNA pathway recognizable in the rest of the Metazoa. In either case, ctenophores represent an intriguing model for better understanding the early evolution of small RNA-based regulatory functions, shedding light on a point in evolutionary time that may have predated the need for additional plasticity in key molecular systems inherent to animals. We expect that further exploration of the genomes of other ctenophores, early branching metazoans, and closely related non-metazoans will help determine the exact point in evolutionary history at which both canonical and mirtron-based miRNA pathways (and their components) emerged.
Two RNA sources were used for sequencing miRNAs. Sample 1 was collected in Woods Hole, MA from mixed stage late embryos 15–30 hours post-fertilization. Total RNA was prepped with TRI-Reagent. Sample 2 was collected in Miami, FL from mixed stage embryos 0–30 hours post-fertilization. Total RNA was prepped with TRIzol Reagent and resuspended in 50 μl of THE RNA solution spiked with RNAsecure.
Sequencing of short RNAs and genome mapping
Libraries of small RNAs were prepared from 5 μg total RNA using Illumina’s Small RNA Alternative v1.5 Sample Prep Protocol with the following modifications. Adapter ligation times were increased from 1 hour to 6 hours, a total of 15 PCR cycles were used, and a 10% acrylamide gel was used for better resolution of properly ligated sequences from unligated free adapters. Sequencing of adapter libraries was performed on an Illumina GAiix using version 5 chemistry and RTA version 18.104.22.168. Both runs were 36-cycle single read. Raw sequencing data was post-processed using CASAVA 1.7.0 and deposited in the NCBI Short Read Archive (http://www.ncbi.nlm.nih.gov/sra/), accession SRA057204.
The 3′ adapter sequence ATCTCGTATGCCGTCTTCTGCTTGT was trimmed from reads using Novocraft’s Novoalign v2.07.18. After filtering reads of low quality, we mapped the trimmed reads to the Mnemiopsis genome independently with both Novoalign and Bowtie v0.12.  (allowing up to two mismatches). Novoalign successfully mapped 65.9% of reads from sample 1 (out of 14,965,804 reads after removal of an overrepresented, unannotated rRNA transcript) and 58.5% of reads from sample 2 to the genome (out of 30,311,098 reads). Bowtie mapped 68.3% and 66.7% of reads from each sample, respectively. Rough estimates showed that ~94% of read mappings from sample 1 were represented in sample 2 and, conversely, ~91% of read mappings from sample 2 were represented in sample 1. This indicates that differences in samples and sequencing protocols did not significantly affect read sources.
Canonical miRNA prediction
miRDeep2  and miRanalyzer (version 0.2)  were used to predict miRNAs from our short RNA sequence data and the Mnemiopsis genome. Candidate predictions were restricted to those present in both samples in at least one read. Next, candidate miRNAs were ranked by the number of methods predicting them, where identification in both methods was considered most confident and predictions by miRDeep2-only were favored over miRanalyzer-only. This ranking is a result of noise filtering to reduce false positives in miRDeep2, producing fewer predictions (143 in sample 1 and 248 in sample 2 with miRDeep2, versus 4197 in sample 1 and 9056 in sample 2 with miRanalyzer).
For miRDeep2, we used all metazoan mature miRNA sequences in miRBase (http://mirbase.org/ftp.shtml) as the input set of known miRNAs. This is used to identify potentially conserved miRNAs, in addition to providing a template for estimating the false positive rate and signal-to-noise ratio at different score cutoffs . No known metazoan miRNAs, including those of other early branching metazoans studied in this work, were identified in the Mnemiopsis samples based on strict sequence similarity having identical seed sequences (nucleotides 2–7) and a maximum of three mismatches in the remaining mature or mature-star arm . The reported signal-to-noise distributions for each sample were notably dissimilar to those reported in other species with known miRNAs . The signal-to-noise ratio is expected to be roughly monotonically increasing with respect to miRDeep2 scores and, in other species including Nematostella, should provide a true positive score cutoff at which signal-to-noise is 10:1, or in the worst case (sea squirt), at least 3.5:1. In our samples, the signal-to-noise ratio peaks at 1.6:1 and 1.3:1, respectively at a score cutoff of 4, and drops off at higher scores (Additional file 3: Dataset 2e & 2h). Although in those experiments the input set of known miRNAs was specific to a single species, opposed to all metazoans, the distributions of signal-to-noise ratio versus score cutoffs does not appear high enough to make any positive predictions in our experiments. Further, our top predictions were sample-specific.
For miRanalyzer, we used all Rfam sequences, provided automatically by the program, to identify known miRNAs and to filter short RNA sequences from other sources. In both samples, no known miRNA mature or mature-star sequences were identified. We did not use miRanalyzer predictions alone to identify novel miRNAs because of the immense number of predictions made. Manual analysis showed that the most highly expressed predictions corresponded to rRNA sequences. We therefore only used miRanalyzer predictions to support miRDeep2 predictions.
The best predictions over all samples and methods were made by miRDeep2 on sample 2. Thus, in addition to looking at the top predictions using the combinatorial criteria described above, we also looked at miRDeep2 predictions for each sample independently. No predicted miRNA had the ideal combination of read mapping signature and secondary structure to be considered a confident miRNA. Top miRDeep2 predictions for each sample are summarized in Additional file 3: Dataset 2a-b. Raw prediction outputs are provided in Additional file 3: Dataset 2c-h.
Finally, in the absence of confident miRNA predictions by the methods described above, we searched the Mnemiopsis genome specifically for miR-100 and miR-2022, as these miRNAs are the only known miRNAs (to our knowledge) thought to be conserved outside of the Bilateria; miR-100 appears to be conserved between Nematostella and bilaterians, while miR-2022 appears to be conserved between Nematostella and Hydra. Querying the Mnemiopsis genome with BLASTN using the conserved portions of the respective mature sequences (miR-100: ACCCGTAGATCCGAACTTGTG, miR-2022: TTTGCTAGTTGCTTTTGTCCC) yielded partial hits in both cases (14 and 16-nucleotide identity, respectively). However, only one hit (for miR-2022 on scaffold ML1502) covered the expected seed site, and no short RNA sequencing reads from either sample mapped to this region. In all, these results support the absence of miR-100 and miR-2022 in Mnemiopsis in addition to all other canonical miRNAs.
The basis of our mirtron prediction method was the combination of an absolute count of mapping reads from Bowtie  and predicted secondary structures by UNAFold  scored using an SVM approach trained on fly mirtrons . All introns of length 50 to 120 nt in Mnemiopsis were considered candidate mirtrons (3953 total, Additional file 3: Dataset 2k) and scored by the SVM based on secondary structure alone. For every candidate mirtron, we independently counted the number of reads pooled from both samples mapping in the correct orientation to the 3′ or 5′ splice sites, with a three-nucleotide buffer in both directions. Our strict read mapping criteria was meant to identify the most likely candidates; while mirtron reads can be found further from the splice sites in other species, the majority of reads tend to fall in this range. We produced three rankings of candidate mirtrons based on the highest scored secondary structures, most correctly mapping reads, and finally by the intersection of the two. Our results did not uncover any high-confidence mirtron candidates. Scoring of the secondary structures resulted in noticeably fewer and lower quality predictions compared to scores reported on Drosophila melanogaster and Caenorhabditis elegans introns  (Additional file 4: Figure S2).
We analyzed introns up to length 150 nt (7324 additional introns from those length 50–120 nt) in the case that Mnemiopsis mirtrons, like Amphimedon miRNAs , were longer than those of flies. The intron length distribution can be seen in Additional file 5: Figure S3. We produced a ranked list based on read counts and manually analyzed the secondary structures of the most highly expressed. Again, no acceptable mirtron candidates were identified.
The best candidates had very low read counts and generally hit only one of the two splice sites; if they are truly functional mirtrons, they are not expressed at high enough levels to be concluded as functional. In addition, their secondary structure predictions were less than ideal relative to known mirtrons in other species. The best identified mirtron candidate (scaffold ML4098, from 40399–40490 on the ‘+’ strand) contains only seven reads total from a single sample (sample 2), six at the 5′ splice site and one at the 3′ splice site, and does not have a characteristic loop or 5′/3′ overhang structure. See Additional file 6: Figures S4-S8 for a summary of the best manually curated predictions, based on the combination of predicted secondary structure and read mappings.
Annotation of miRNA pathway proteins
RNAi pathway proteins identified in Mnemiopsis throughout the course of this study have been deposited in GenBank (http://www.ncbi.nlm.nih.gov/Genbank/), with accessions JQ437405 (Dicer), JQ437406 (Argonaute), JQ437407 (Exportin-5), and JQ437408 (Ran). Two additional Argonaute family members were annotated: JX483728 and JX483729. Identification and annotation of Mnemiopsis proteins was based on high-scoring reciprocal BLASTP hits to the human RefSeq protein set. TBLASTN was also used but did not identify any better candidates. Human Dicer and Drosha both hit uniquely to the same Mnemiopsis protein, but reciprocal BLASTP results favored Dicer. The protein models of all species represented in Figure 4 were searched with HMMER 3.0  for tandem RNase III domains; no Dicer or Drosha candidates were identified in the closest non-metazoan outgroups (i.e., Monosiga brevicollis, Salpingoeca rosetta, Capsapora owczarzaki and Sphaeroforma arctica). Nematostella, Hydra, Trichoplax, and Amphimedon protein sequence data were downloaded from the Joint Genome Institute (JGI) Web site and protein sequence data for the closest non-metazoan outgroups were downloaded from the Origins of Multicellularity Sequencing Project Web site of the Broad Institute of Harvard and MIT (http://www.broadinstitute.org/) in November 2011. In some of these species, the RNase III domains of Dicer and Drosha proteins were not properly annotated. In these cases, we instead used published, manually curated sequences  or the appropriate RefSeq entries when those were not available. Other RNase III sequences from the bilateria and eubacteria included in our analysis were selected from sequences used in a previous study  or sampled from RefSeq and GenBank. All accession numbers for RNase III enzymes included in our final analysis are reported in Additional file 7: Table S1. The trimmed RNase III domain sequences used to build the phylogenetic tree in Figure 3 were aligned with HMMER 3.0  and manually padded in cases where terminal gaps could be reliably filled. Residues 59–98 were manually trimmed from the alignment based on poor conservation. Both alignments are reported in Additional file 1: Dataset 1a-b.
Figure 3 was generated to better-categorize the Mnemiopsis RNase III enzyme as a Dicer or Drosha and to better-understand the origin of Drosha. This phylogenetic tree was built on the trimmed alignment described above. ProtTest v2.4  was used to pick the best model of evolution and selected the LG model with optimization of substitution rates, gamma model of rate heterogeneity, and empirical amino acid frequencies (PROTGAMMAILGF model). We used RAxML v7.2.8a  to build trees seeded on 24 random starting trees and 24 maximum parsimony trees. We also ran MrBayes v3.1.2  to construct a Bayesian tree, using five million iterations on five chains with a burn-in factor of 25%. MrBayes was run using the second best model selected by ProtTest since the LG model is not available in MrBayes: RtRev with optimized substitution rates, gamma model of rate heterogeneity, and empirical amino acid frequencies. All 49 trees were compared in a maximum likelihood framework, and we reported the tree with the highest likelihood (RAxML with maximum parsimony starting tree, log likelihood = −5895.384778). Support for clades was assessed using 1000 bootstrap replicates and posterior probabilities computed with MrBayes. NEWICK formatted trees are provided in Additional file 1: Dataset 1c-d with bootstraps and Bayesian posterior probabilities.
Ruby JG, Jan CH, Bartel DP: Intronic microRNA precursors that bypass Drosha processing. Nature. 2007, 448: 83-86. 10.1038/nature05983.
Berezikov E, Chung W-J, Willis J, Cuppen E, Lai EC: Mammalian mirtron genes. Mol Cell. 2007, 28: 328-336. 10.1016/j.molcel.2007.09.028.
Axtell MJ, Westholm JO, Lai EC: Vive la différence: biogenesis and evolution of microRNAs in plants and animals. Genome Biol. 2011, 12: 221-10.1186/gb-2011-12-4-221.
Schickel R, Boyerinas B, Park S-M, Peter ME: MicroRNAs: key players in the immune system, differentiation, tumorigenesis and cell death. Oncogene. 2008, 27: 5959-5974. 10.1038/onc.2008.274.
Liu N, Okamura K, Tyler DM, Phillips MD, Chung W-J, Lai EC: The evolution and functional diversification of animal microRNA genes. Cell Res. 2008, 18: 985-996. 10.1038/cr.2008.278.
Berezikov E: Evolution of microRNA diversity and regulation in animals. Nat Rev Genet. 2011, 12: 846-860. 10.1038/nrg3079.
Shabalina SA, Koonin EV: Origins and evolution of eukaryotic RNA interference. Trends Ecol Evol. 2008, 23: 578-587. 10.1016/j.tree.2008.06.005.
Kosik KS: MicroRNAs tell an evo-devo story. Nat Rev Neurosci. 2009, 10: 754-759. 10.1038/nrn2713.
Grimson A, Srivastava M, Fahey B, Woodcroft BJ, Chiang HR, King N, Degnan BM, Rokhsar DS, Bartel DP: Early origins and evolution of microRNAs and piwi-interacting RNAs in animals. Nature. 2008, 455: 1193-1197. 10.1038/nature07415.
Cerutti H, Casas-Mollano JA: On the origin and functions of RNA-mediated silencing: from protists to man. Curr Genet. 2006, 50: 81-99. 10.1007/s00294-006-0078-x.
Kim VN, Han J, Siomi MC: Biogenesis of small RNAs in animals. Nat Rev Mol Cell Bio. 2009, 10: 126-139. 10.1038/nrm2632.
Peterson KJ, Dietrich MR, McPeek MA: MicroRNAs and metazoan macroevolution: insights into canalization, complexity, and the Cambrian explosion. BioEssays. 2009, 31: 736-747. 10.1002/bies.200900033.
Wheeler BM, Heimberg AM, Moy VN, Sperling EA, Holstein TW, Heber S, Peterson KJ: The deep evolution of metazoan microRNAs. EvoDevo. 2009, 11: 50-68.
Niwa R, Slack FJ: The evolution of animal microRNA function. Curr Opin Genet Dev. 2007, 17: 145-150. 10.1016/j.gde.2007.02.004.
Bartel DP: MicroRNAs: genomics, biogenesis, mechanism, and function. Cell. 2004, 116: 281-297. 10.1016/S0092-8674(04)00045-5.
Jones-Rhoades MW, Bartel DP, Bartel B: MicroRNAs and their regulatory roles in plants. Annu Rev Plant Biol. 2006, 57: 19-53. 10.1146/annurev.arplant.57.032905.105218.
Hinas A, Reimegard J, Wagner EGH, Nellen W, Ambros VR, Soderbom F: The small RNA repertoire of Dictyostelium discoideum and its regulation by components of the RNAi pathway. Nucleic Acids Res. 2007, 35: 6714-6726. 10.1093/nar/gkm707.
Cock JM, Sterck L, Rouzé P, Scornet D, Allen AE, Amoutzias G, Anthouard V, Artiguenave F, Aury J-M, Badger JH, et al: The Ectocarpus genome and the independent evolution of multicellularity in brown algae. Nature. 2010, 465: 617-621. 10.1038/nature09016.
Huang A, He L, Wang G: Identification and characterization of microRNAs from Phaeodactylum tricornutum by high-throughput sequencing and bioinformatics analysis. BMC Genomics. 2011, 12: 337-10.1186/1471-2164-12-337.
Lin W-C, Li S-C, Lin W-C, Shin J-W, Hu S-N, Yu X-M, Huang T-Y, Chen S-C, Chen H-C, Chen S-J, et al: Identification of microRNA in the protist Trichomonas vaginalis. Genomics. 2009, 93: 487-493. 10.1016/j.ygeno.2009.01.004.
Chen XS, Collins LJ, Biggs PJ, Penny D: High throughput genome-wide survey of small RNAs from the parasitic protists Giardia intestinalis and Trichomonas vaginalis. Genome Biol Evol. 2009, 1: 165-175.
Huang P-J, Lin W-C, Chen S-C, Lin Y-H, Sun C-H, Lyu P-C, Tang P: Identification of putative miRNAs from the deep-branching unicellular flagellates. Genomics. 2012, 99: 101-107. 10.1016/j.ygeno.2011.11.002.
Lin W-C, Huang K-Y, Chen S-C, Huang T-Y, Chen S-J, Huang P-J, Tang P: Malate dehydrogenase is negatively regulated by miR-1 in Trichomonas vaginalis. Parasitol Res. 2009, 105: 1683-1689. 10.1007/s00436-009-1616-5.
Li W, Saraiya AA, Wang CC: Gene regulation in Giardia lamblia involves a putative microRNA derived from a small nucleolar RNA. PLoS Negl Trop Dis. 2011, 5: e1338-10.1371/journal.pntd.0001338.
Saraiya AA, Li W, Wang CC: A microRNA derived from an apparent canonical biogenesis pathway regulates variant surface protein gene expression in Giardia lamblia. RNA. 2011, 17: 2152-2164. 10.1261/rna.028118.111.
Saraiya AA, Wang CC: SnoRNA, a novel precursor of microRNA in Giardia lamblia. PLoS Pathog. 2008, 4: e1000224-10.1371/journal.ppat.1000224.
Braun L, Cannella D, Ortet P, Barakat M, Sautel CF, Kieffer S, Garin J, Bastien O, Voinnet O, Hakimi MA: A complex small RNA repertoire is generated by a plant/fungal-like machinery and effected by a metazoan-like Argonaute in the single-cell human parasite Toxoplasma gondii. PLoS Pathog. 2010, 6: e1000920-10.1371/journal.ppat.1000920.
Molnár A, Schwach F, Studholme DJ, Thuenemann EC, Baulcombe DC: MiRNAs control gene expression in the single-cell alga Chlamydomonas reinhardtii. Nature. 2007, 447: 1126-1129. 10.1038/nature05903.
Zhao T, Li G, Mi S, Li S, Hannon GJ, Wang XJ, Qi Y: A complex system of small RNAs in the unicellular green alga Chlamydomonas reinhardtii. Genes Dev. 2007, 21: 1190-1203. 10.1101/gad.1543507.
Tarver JE, Donoghue PCJ, Peterson KJ: Do miRNAs have a deep evolutionary history?. BioEssays. 2012, 34: 857-866. 10.1002/bies.201200055.
Dunn CW, Hejnol A, Matus DQ, Pang K, Browne WE, Smith SA, Seaver E, Rouse GW, Obst M, Edgecombe GD, et al: Broad phylogenomic sampling improves resolution of the animal tree of life. Nature. 2008, 452: 745-749. 10.1038/nature06614.
Hejnol A, Obst M, Stamatakis A, Ott M, Rouse GW, Edgecombe GD, Martinez P, Baguñà J, Bailly X, Jondelius U, et al: Assessing the root of bilaterian animals with scalable phylogenomic methods. Proc Biol Sci. 2009, 276: 4261-4270. 10.1098/rspb.2009.0896.
Ryan JF, Pang K, Comparative Sequencing Program N, Mullikin JC, Martindale MQ, Baxevanis AD: The homeodomain complement of the ctenophore Mnemiopsis leidyi suggests that Ctenophora and Porifera diverged prior to the ParaHoxozoa. EvoDevo. 2010, 1: 9-10.1186/2041-9139-1-9.
Reitzel AM, Pang K, Ryan JF, Mullikin JC, Martindale MQ, Baxevanis AD, Tarrant AM: Nuclear receptors from the ctenophore Mnemiopsis leidyi lack a zinc-finger DNA-binding domain: lineage-specific loss or ancestral condition in the emergence of the nuclear receptor superfamily?. EvoDevo. 2011, 2: 3-10.1186/2041-9139-2-3.
Pang K, Ryan JF, Baxevanis AD, Martindale MQ: Evolution of the TGF-β signaling pathway and its potential role in the ctenophore, Mnemiopsis leidyi. PLoS ONE. 2011, 6: e24152-10.1371/journal.pone.0024152.
Pick KS, Philippe H, Schreiber F, Erpenbeck D, Jackson DJ, Wrede P, Wiens M, Alie A, Morgenstern B, Manuel M, et al: Improved phylogenomic taxon sampling noticeably affects nonbilaterian relationships. Mol Biol Evol. 2010, 27: 1983-1987. 10.1093/molbev/msq089.
Chapman JA, Kirkness EF, Simakov O, Hampson SE, Mitros T, Weinmaier T, Rattei T, Balasubramanian PG, Borman J, Busam D, et al: The dynamic genome of Hydra. Nature. 2010, 464: 592-596. 10.1038/nature08830.
Hertel J, de Jong D, Marz M, Rose D, Tafer H, Tanzer A, Schierwater B, Stadler PF: Non-coding RNA annotation of the genome of Trichoplax adhaerens. Nucleic Acids Res. 2009, 37: 1602-1615. 10.1093/nar/gkn1084.
MacRae IJ, Doudna JA: Ribonuclease revisited: structural insights into ribonuclease III family enzymes. Cur Opin Struct Biol. 2007, 17: 138-145. 10.1016/j.sbi.2006.12.002.
Mochizuki K: A Dicer-like protein in Tetrahymena has distinct functions in genome rearrangement, chromosome segregation, and meiotic prophase. Genes Dev. 2005, 19: 77-89. 10.1101/gad.1265105.
MacRae IJ, Zhou K, Doudna JA: Structural determinants of RNA recognition and cleavage by Dicer. Nat Struct Mol Biol. 2007, 14: 934-940. 10.1038/nsmb1293.
Finn RD, Clements J, Eddy SR: HMMER web server: interactive sequence similarity searching. Nucleic Acids Res. 2011, 39: W29-W37. 10.1093/nar/gkr367.
de Jong D, Eitel M, Jakob W, Osigus H-J, Hadrys H, Desalle R, Schierwater B: Multiple Dicer genes in the early-diverging metazoa. Mol Biol Evol. 2009, 26: 1333-1340. 10.1093/molbev/msp042.
Zhang L, Hou D, Chen X, Li D, Zhu L, Zhang Y, Li J, Bian Z, Liang X, Cai X, et al: Exogenous plant MIR168a specifically targets mammalian LDLRAP1: evidence of cross-kingdom regulation by microRNA. Cell Res. 2011, 1-20.
Flynt AS, Greimann JC, Chung W-J, Lima CD, Lai EC: MicroRNA biogenesis via splicing and exosome-mediated trimming in Drosophila. Molecular Cell. 2010, 38: 900-907. 10.1016/j.molcel.2010.06.014.
Philippe H, Brinkmann H, Copley RR, Moroz LL, Nakano H, Poustka AJ, Wallberg A, Peterson KJ, Telford MJ: Acoelomorph flatworms are deuterostomes related to Xenoturbella. Nature. 2011, 470: 255-258. 10.1038/nature09676.
Langmead B, Trapnell C, Pop M, Salzberg SL: Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol. 2009, 10: R25-10.1186/gb-2009-10-3-r25.
Friedländer MR, Mackowiak SD, Li N, Chen W, Rajewsky N: miRDeep2 accurately identifies known and hundreds of novel microRNA genes in seven animal clades. Nucleic Acids Res. 2012, 40: 37-52. 10.1093/nar/gkr688.
Hackenberg M, Sturm M, Langenberger D, Falcon-Perez JM, Aransay AM: miRanalyzer: a microRNA detection and analysis tool for next-generation sequencing experiments. Nucleic Acids Res. 2009, 37: W68-W76. 10.1093/nar/gkp347.
Markham NR, Zuker M: UNAFold: software for nucleic acid folding and hybridization. Methods Mol Biol. 2008, 453: 3-31. 10.1007/978-1-60327-429-6_1.
Chung WJ, Agius P, Westholm JO, Chen M, Okamura K, Robine N, Leslie CS, Lai EC: Computational and experimental identification of mirtrons in Drosophila melanogaster and Caenorhabditis elegans. Genome Res. 2011, 21: 286-300. 10.1101/gr.113050.110.
Abascal F, Zardoya R, Posada D: ProtTest: selection of best-fit models of protein evolution. Bioinformatics. 2005, 21: 2104-2105. 10.1093/bioinformatics/bti263.
Stamatakis A, Hoover P, Rougemont J: A rapid bootstrap algorithm for the RAxML Web servers. Syst Biol. 2008, 57: 758-771. 10.1080/10635150802429642.
Ronquist F, Huelsenbeck JP: MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics. 2003, 19: 1572-1574. 10.1093/bioinformatics/btg180.
The authors would like to thank the NIH Intramural Sequencing Center, particularly A. Young, for performing the small RNA sequencing and describing the protocol, K. Pang for providing samples, M. Martindale for reviewing the manuscript, D. Gildea for assistance with short RNA sequencing analysis, J. Fekecs and D. Leja for assistance in the creation and editing of figures, M. Srivastava for input regarding selection of protein sequences, N. Trivedi for input on figure design, and A. Nguyen for assistance with miRNA predictions. This work was supported by an NIH Graduate Research Fellowship (E.K.M.), by the Intramural Research Program of the National Human Genome Research Institute, National Institutes of Health (E.K.M., J.F.R., C.E.S., A.D.B.), and by the University of Miami, College of Arts and Sciences and Provost Research Award (W.E.B.).
The authors declare no competing interests.
EKM performed the majority of computational analyses and was primary author of the manuscript. JFR, CES, WEB and ADB contributed to performing the miRNA predictions, protein/pathway identification, and phylogenetic analyses. WEB performed experimental analysis. All authors contributed to the design of the study and preparation of the manuscript.