Abstract:It has been known for some time that there are many weakly expressed and fast-evolving microRNAs (miRNAs)1-6. After analyzing more than 100,000 reads of small RNAs obtained from Drosophila heads, we suggested that most of these miRNAs were born and then died with the evolutionary dynamics of neutrally evolving sequences7. The birth and death rates of miRNAs were estimated to be about 12 and 11.7 genes per Myr, respectively, resulting in a fairly modest net gain of only 0.3 genes per Myr. Berezikov et al. revised the estimate of net gain to be about 1 gene per Myr but did not provide an estimate of the birth and death rates separately8. Instead, they argued that our estimates for both the birth and death rates were too high by claiming that many of the miRNAs in our analysis were not miRNAs at all but were merely degraded products of mRNAs. Their main concern was that many of the newly born miRNAs in our observation were singletons. In their expanded data, many of these singletons were missing and some were accompanied by sequences of similar length in the same hairpin, both of which suggested RNA degradation. Being mindful of the limitation of low coverage in our study7, we have collected 18 million small RNA reads in Drosophila heads by sequencing of oligonucleotide ligation and detection (SOLiD) sequencing. We analyzed only reads that appeared at least 50 times in the arm of a hairpin and for which the accurate processing rate of the 5′ end of the miRNAs was ≥90%. The new dataset shows that the conclusion of Lu et al.7 was not biased by low coverage, as suggested by Berezikov et al.8. What, then, may be the reasons for the discrepancy between the analysis of Berezikov et al.8 and our analysis7? Due to space limitations, we shall only give a brief account and will present a more thorough comparison on our website (http://pondside.uchicago.edu/wulab/microRNA/). First, Berezikov et al.8 defined miRNAs much more narrowly than we believe is reasonable. They did not consider hairpins on exons—a legitimate source of miRNAs, as many known miRNAs share sequences with exons9-11. In addition, they used criteria derived from conserved miRNAs to screen out candidates. These criteria include a much lower rate of evolution in the arms of the hairpin than in the loop and also the narrow size distribution of mature miRNAs. If miRNAs are defined by conservation, then there would be few that are evolutionarily transient. Second, different sequencing platforms are known to yield fairly different results2,12. By sequencing platform, we mean more than just sequencing chemistry; rather, we include the entire protocols, from upstream library preparation to base callings, recommended by the manufacturers of Roche 454 GS-FLX, IlluminaGA or SOLiD-ABI DNA sequencers. Lu et al.7, Berezikov et al.8 and our new study (unpublished) used the 454, GA and SOLiD methods, respectively. Even the most abundant miRNAs were recorded with surprisingly large disparity by different methods, and the rare miRNAs were not detected by all methods. Our SOLiD datasets indeed contain many genes with multiple reads that were absent in the GA data Berezikov et al. used8. In another accompanying paper (Zhou et al., unpublished), we show that GA and SOLiD sequencing also perform very differently in SNP detection. Third, Berezikov et al.8 incorrectly used the argument of proportionality, which further confounded the interpretation of rare miRNAs. When there exists a large number of weakly expressed genes, a tenfold increase in coverage might increase the observed number of all rare genes tenfold, but each new gene discovery is likely to be different from the previously observed ones13. The same rare genes should not be expected to occur ten times as often, as Berezikov et al. asserted8. Fourth, the variation in miRNA expression within the same species should not be ignored. In humans, with GA sequencing, we have observed 300 conserved miRNAs that could be detected in fewer than 6 of the 12 human kidney libraries (Lu et al., unpublished data). Because the level of genetic variation in D. melanogaster is 20 times greater than that in humans, we expect the level of miRNA expression polymorphism to be at least as large. Berezikov et al.8 and Lu et al.7 used different fly strains in their surveys. The account above shows that low-frequency small RNAs (including miRNAs) are often seen in some studies (such as Lu et al.7 and our new SOLiD data) but not others (such as Berezikov et al.8). The hasty conclusion that RNA degradation must be at issue seems inappropriate. Finally, we believe that the debate should shift to the rate of net gene gain (birth rate – death rate). After all, there is no disagreement about the conclusion that the transient genes, miRNAs or not, do not appear to be biologically important. We suspect that our original estimate of 0.3 genes per Myr7 might be an underestimate. Because of the moderate number of reads from a single platform (454), many miRNAs were not observed or were observed in one species only. As a result, birth rate might have been underestimated and death rate overestimated. We will provide a revised estimate on the basis of the new and expanded dataset elsewhere.

Meta-Research: understudied genes are lost in a leaky pipeline between genome-wide assays and reporting of results

Meta-Research: Understudied genes are lost in a leaky pipeline between genome-wide assays and reporting of results

Gene Set linkage analysis: a tool for interpreting the overall functional impacts of observed transcriptomic changes

Functional unknomics: Systematic screening of conserved genes of unknown function

Imputation of truncated p-values for meta-analysis methods and its genomic application

Identifying overrepresented concepts in gene lists from literature: a statistical approach based on Poisson mixture model

The landscape of gene loss and missense variation across the mammalian tree informs on gene essentiality

A Longitudinal Analysis of Function Annotations of the Human Proteome Reveals Consistently High Biases

Identification of human gene research articles with wrongly identified nucleotide sequences

Enhanced recovery of single-cell RNA-sequencing reads for missing gene expression data

The Possibility of Systematic Research Fraud Targeting Under-Studied Human Genes: Causes, Consequences, and Potential Solutions

Why an integrated view of gene expression studies on hematopoiesis in mouse aging is better than the sum of their parts

Evolutionary Flux of Canonical Micrornas and Mirtrons in Drosophila

The mutational constraint spectrum quantified from variation in 141,456 humans

Software tools identify forgotten genes

Gene annotation bias impedes biomedical research

Gene length and detection bias in single cell RNA sequencing protocols

No country for old methods: New tools for studying microproteins

Adapting functional genomic tools to metagenomic analyses: investigating the role of gut bacteria in relation to obesity.

A New Machine Learning-Based Framework for Mapping Uncertainty Analysis in RNA-Seq Read Alignment and Gene Expression Estimation

A Re-Analysis of an Existing Dataset Reveals a New Set of Genes Involved in Post-Mating Response