Human YRNA 4 (HY4) plasma levels are a prognostic indicator of SARS-CoV-2 infection clinical severity.

SARS-CoV-2 infection can result in a range of outcomes from asymptomatic/mild disease to severe COVID-19/fatality. In this study, we investigated the differential expression of small noncoding RNAs (sncRNAs) between patient cohorts defined by disease severity. We collected plasma samples, stratified these based on clinical outcomes, and sequenced their circulating sncRNAs. Excitingly, we found YRNA HY4 displays significant differential expression (p=0.025) between patients experiencing mild and severe disease. In agreement with recent reports identifying plasma YRNAs as indicators of influenza infection severity, our results strongly suggest that circulating HY4 levels represent a powerful prognostic indicator of likely SARS-CoV-2 patient infection outcome.

PFKFB3 Inhibits Fructose Metabolism in Pulmonary Microvascular Endothelial Cells.

Pulmonary microvascular endothelial cells contribute to the integrity of the lung gas exchange interface, and they are highly glycolytic. Although glucose and fructose represent discrete substrates available for glycolysis, pulmonary microvascular endothelial cells prefer glucose over fructose, and the mechanisms involved in this selection are unknown. 6-Phosphofructo-2-kinase/fructose-2, 6-bisphosphatase 3 (PFKFB3) is an important glycolytic enzyme that drives glycolytic flux against negative feedback and links glycolytic and fructolytic pathways. We hypothesized that PFKFB3 inhibits fructose metabolism in pulmonary microvascular endothelial cells. We found that PFKFB3 knockout cells survive better than wild-type cells in fructose-rich medium under hypoxia. Seahorse assays, lactate and glucose measurements, and stable isotope tracing showed that PFKFB3 inhibits fructose-hexokinase-mediated glycolysis and oxidative phosphorylation. Microarray analysis revealed that fructose upregulates PFKFB3, and PFKFB3 knockout cells increase fructose-specific GLUT5 (glucose transporter 5) expression. Using conditional endothelial-specific PFKFB3 knockout mice, we demonstrated that endothelial PFKFB3 knockout increases lung tissue lactate production after fructose gavage. Last, we showed that pneumonia increases fructose in BAL fluid in mechanically ventilated ICU patients. Thus, PFKFB3 knockout increases GLUT5 expression and the hexokinase-mediated fructose use in pulmonary microvascular endothelial cells that promotes their survival. Our findings indicate that PFKFB3 is a molecular switch that controls glucose versus fructose use in glycolysis and help better understand lung endothelial cell metabolism during respiratory failure.

RUNX1 Is Regulated by Androgen Receptor to Promote Cancer Stem Markers and Chemotherapy Resistance in Triple Negative Breast Cancer.

Triple negative breast cancer (TNBC) is an aggressive breast cancer subtype for which no effective targeted therapies are available. Growing evidence suggests that chemotherapy-resistant cancer cells with stem-like properties (CSC) may repopulate the tumor. The androgen receptor (AR) is expressed in up to 50% of TNBCs, and AR inhibition decreases CSC and tumor initiation. Runt-related transcription factor 1 (RUNX1) correlates with poor prognosis in TNBC and is regulated by the AR in prostate cancer. Our group has shown that RUNX1 promotes TNBC cell migration and regulates tumor gene expression. We hypothesized that RUNX1 is regulated by the AR and that both may work together in TNBC CSC to promote disease recurrence following chemotherapy. Chromatin immunoprecipitation sequencing (ChIP-seq) experiments in MDA-MB-453 revealed AR binding to RUNX1 regulatory regions. RUNX1 expression is upregulated by dihydrotestosterone (DHT) in MDA-MB-453 and in an AR+-TNBC HCI-009 patient-derived xenograft (PDX) tumors (p < 0.05). RUNX1 is increased in a CSC-like experimental model in MDA-MB-453 and SUM-159PT cells (p < 0.05). Inhibition of RUNX1 transcriptional activity reduced the expression of CSC markers. Interestingly, RUNX1 inhibition reduced cell viability and enhanced paclitaxel and enzalutamide sensitivity. Targeting RUNX1 may be an attractive strategy to potentiate the anti-tumor effects of AR inhibition, specifically in the slow-growing CSC-like populations that resist chemotherapy which lead to metastatic disease.

A single N6-methyladenosine site regulates lncRNA HOTAIR function in breast cancer cells.

N6-methyladenosine (m6A) modification of RNA regulates normal and cancer biology, but knowledge of its function on long noncoding RNAs (lncRNAs) remains limited. Here, we reveal that m6A regulates the breast cancer-associated human lncRNA HOTAIR. Mapping m6A in breast cancer cell lines, we identify multiple m6A sites on HOTAIR, with 1 single consistently methylated site (A783) that is critical for HOTAIR-driven proliferation and invasion of triple-negative breast cancer (TNBC) cells. Methylated A783 interacts with the m6A "reader" YTHDC1, promoting chromatin association of HOTAIR, proliferation and invasion of TNBC cells, and gene repression. A783U mutant HOTAIR induces a unique antitumor gene expression profile and displays loss-of-function and antimorph behaviors by impairing and, in some cases, causing opposite gene expression changes induced by wild-type (WT) HOTAIR. Our work demonstrates how modification of 1 base in an lncRNA can elicit a distinct gene regulation mechanism and drive cancer-associated phenotypes.

MicroRNA-like snoRNA-Derived RNAs (sdRNAs) Promote Castration-Resistant Prostate Cancer.

We have identified 38 specifically excised, differentially expressed snoRNA fragments (sdRNAs) in TCGA prostate cancer (PCa) patient samples as compared to normal prostate controls. SnoRNA-derived fragments sdRNA-D19b and -A24 emerged among the most differentially expressed and were selected for further experimentation. We found that the overexpression of either sdRNA significantly increased PC3 (a well-established model of castration-resistant prostate cancer (CRPC)) cell proliferation, and that sdRNA-D19b overexpression also markedly increased the rate of PC3 cell migration. In addition, both sdRNAs provided drug-specific resistances with sdRNA-D19b levels correlating with paclitaxel resistance and sdRNA-24A conferring dasatinib resistance. In silico and in vitro analyses revealed that two established PCa tumor suppressor genes, CD44 and CDK12, represent targets for sdRNA-D19b and sdRNA-A24, respectively. This outlines a biologically coherent mechanism by which sdRNAs downregulate tumor suppressors in AR-PCa to enhance proliferative and metastatic capabilities and to encourage chemotherapeutic resistance. Aggressive proliferation, rampant metastasis, and recalcitrance to chemotherapy are core characteristics of CRPC that synergize to produce a pathology that ranks second in cancer-related deaths for men. This study defines sdRNA-D19b and -A24 as contributors to AR-PCa, potentially providing novel biomarkers and therapeutic targets of use in PCa clinical intervention.

Establishing RNA-RNA interactions remodels lncRNA structure and promotes PRC2 activity.

Human Polycomb Repressive Complex 2 (PRC2) catalysis of histone H3 lysine 27 methylation at certain loci depends on long noncoding RNAs (lncRNAs). Yet, in apparent contradiction, RNA is a potent catalytic inhibitor of PRC2. Here, we show that intermolecular RNA-RNA interactions between the lncRNA HOTAIR and its targets can relieve RNA inhibition of PRC2. RNA bridging is promoted by heterogeneous nuclear ribonucleoprotein B1, which uses multiple protein domains to bind HOTAIR regions via multivalent protein-RNA interactions. Chemical probing demonstrates that establishing RNA-RNA interactions changes HOTAIR structure. Genome-wide HOTAIR/PRC2 activity occurs at genes whose transcripts can make favorable RNA-RNA interactions with HOTAIR. We demonstrate that RNA-RNA matches of HOTAIR with target gene RNAs can relieve the inhibitory effect of a single lncRNA for PRC2 activity after B1 dissociation. Our work highlights an intrinsic switch that allows PRC2 activity in specific RNA contexts, which could explain how many lncRNAs work with PRC2.

Identification of m6A residues at single-nucleotide resolution using eCLIP and an accessible custom analysis pipeline.

Methylation at the N6 position of adenosine (m6A) is one of the most abundant RNA modifications found in eukaryotes; however, accurate detection of specific m6A nucleotides within transcripts has been historically challenging due to m6A and unmodified adenosine having virtually indistinguishable chemical properties. While previous strategies such as methyl-RNA immunoprecipitation and sequencing (MeRIP-seq) have relied on m6A-specific antibodies to isolate RNA fragments containing the modification, these methods do not allow for precise identification of individual m6A residues. More recently, modified cross-linking and immunoprecipitation (CLIP)-based approaches that rely on inducing specific mutations during reverse transcription via UV cross-linking of the anti-m6A antibody to methylated RNA have been used to overcome this limitation. However, the most utilized version of this approach, miCLIP, can be technically challenging to use for achieving high-complexity libraries. Here we present an improved methodology that yields high library complexity and allows for the straightforward identification of individual m6A residues with reliable confidence metrics. Based on enhanced CLIP (eCLIP), our m6A-eCLIP (meCLIP) approach couples the improvements of eCLIP with the inclusion of an input sample and an easy-to-use computational pipeline to allow for precise calling of m6A sites at true single-nucleotide resolution. As the effort to accurately identify m6As in an efficient and straightforward way intensifies, this method is a valuable tool for investigators interested in unraveling the m6A epitranscriptome.

A Positive Feedback Loop Between TGFß and Androgen Receptor Supports Triple-negative Breast Cancer Anoikis Resistance.

Triple-negative breast cancer (TNBC) is an aggressive subtype with peak recurrence as metastatic disease within the first few years of diagnosis. Androgen receptor (AR) expression is increased in anchorage-independent cells in TNBC preclinical models. Both AR knockdown and inhibition lead to reduced TNBC invasion in vitro, reduced tumorigenicity, and less recurrence in vivo in preclinical models. Transforming growth factor ß (TGFß) pathway gene signatures also increased during anchorage-independent survival both in vitro and in vivo in preclinical models and in circulating tumor cells (CTCs) from patients during emergence of chemo resistant disease. We hypothesized that a positive loop between AR and TGFß signaling facilitates TNBC anchorage-independent survival. We find that multiple components of the TGFß pathway, including TGFß1 and 3, as well as pathway activity measured by nuclear localization and transcriptional activity of phosphorylated Smad3, are enhanced in anchorage-independent conditions. Further, exogenous TGFß increased AR protein while TGFß inhibition decreased AR and TNBC viability, particularly under anchorage-independent culture conditions. ChIP-seq experiments revealed AR binding to TGFB1 and SMAD3 regulatory regions in MDA-MB-453 cells. In clinical datasets, TGFB3 and AR positively correlate and high expression of both genes together corresponded to significantly worse recurrence-free and overall survival in both ER-negative and basal-like breast cancer. Finally, inhibiting both AR and TGFß decreased cell survival, particularly under anchorage-independent conditions. These findings warrant further investigations into whether combined inhibition of AR and TGFß pathways might decrease metastatic recurrence rates and mortality from TNBC.

RNA matchmaking in chromatin regulation.

Beyond being the product of gene expression, RNA can also influence the regulation of chromatin. The majority of the human genome has the capacity to be transcribed and the majority of the non-protein-coding transcripts made by RNA Polymerase II are enriched in the nucleus. Many chromatin regulators can bind to these ncRNAs in the nucleus; in some cases, there are clear examples of direct RNA-mediated chromatin regulation mechanisms stemming from these interactions, while others have yet to be determined. Recent studies have highlighted examples of chromatin regulation via RNA matchmaking, a term we use broadly here to describe intermolecular base-pairing interactions between one RNA molecule and an RNA or DNA match. This review provides examples of RNA matchmaking that regulates chromatin processes and summarizes the technical approaches used to capture these events.

Characterization of novel small RNAs (sRNAs) contributing to the desiccation response of Salmonella enterica serovar Typhimurium.

Noncoding RNA (ncRNA) modulation of gene expression has now been ubiquitously observed across all domains of life. An increasingly apparent role of ncRNAs is to coordinate changes in gene expressions in response to environmental stress. Salmonella enterica, a common food-born pathogen, is known for its striking ability to survive, adapt, and thrive in various unfavourable environments which makes it a particularly difficult pathogen to eliminate as well as an interesting model in which to study ncRNA contributions to cellular stress response. Mounting evidence now suggests that small RNAs (sRNAs) represent key regulators of Salmonella stress adaptation. Approximately 50-500 nucleotides in length, sRNAs regulate gene expression through complementary base pairing with molecular targets and have recently been suggested to outnumber protein-coding genes in bacteria. In this work, we employ small RNA transcriptome sequencing to characterize changes in the sRNA profiles of Salmonella in response to desiccation. In all, we identify 102 previously annotated sRNAs significantly differentially expressed during desiccation; and excitingly, 71 novel sRNAs likewise differentially expressed. Small transcript northern blotting and qRT-PCRs confirm the identities and expressions of several of our novel sRNAs, and computational analyses indicate the majority are highly conserved and structurally related to characterized sRNAs. Predicted sRNA targets include several proteins necessary for desiccation survival and this, in part, suggests a role for desiccation-regulated sRNAs in this stress response. Furthermore, we find individual knock-outs of two of the novel sRNAs identified herein, either sRNA1320429 or sRNA3981754, significantly impairs the ability of Salmonella to survive desiccation, confirming their involvements (and suggesting the potential involvements of other sRNAs we identify in this work) in the Salmonella response to desiccation.

ADAR Mediated RNA Editing Modulates MicroRNA Targeting in Human Breast Cancer.

RNA editing by RNA specific adenosine deaminase acting on RNA (ADAR) is increasingly being found to alter microRNA (miRNA) regulation. Editing of miRNA transcripts can affect their processing, as well as which messenger RNAs (mRNAs) they target. Further, editing of target mRNAs can also affect their complementarity to miRNAs. Notably, ADAR editing is often increased in malignancy with the effect of these RNA changes being largely unclear. In addition, numerous reports have now identified an array of miRNAs that directly contribute to various malignancies although the majority of their targets remain largely undefined. Here we propose that modulating the targets of miRNAs via mRNA editing is a frequent occurrence in cancer and an underappreciated participant in pathology. In order to more accurately characterize the relationship between these two regulatory processes, this study examined RNA editing events within mRNA sequences of two breast cancer cell lines (MCF-7 and MDA-MB-231) and determined whether or not these edits could modulate miRNA associations. Computational analyses of RNA-Seq data from these two cell lines identified over 50,000 recurrent editing sites within human mRNAs, and many of these were located in 3' untranslated regions (UTRs). When these locations were screened against the list of currently-annotated miRNAs we discovered that editing caused a subset (~9%) to have significant alterations to mRNA complementarity. One miRNA in particular, miR-140-3p, is known to be misexpressed in many breast cancers, and we found that mRNA editing allowed this miRNA to directly target the apoptosis inducing gene DFFA in MCF-7, but not in MDA-MB-231 cells. As these two cell lines are known to have distinct characteristics in terms of morphology, invasiveness and physiological responses, we hypothesized that the differential RNA editing of DFFA in these two cell lines could contribute to their phenotypic differences. Indeed, we confirmed through western blotting that inhibiting miR-140-3p increases expression of the DFFA protein product in MCF-7, but not MDA-MB-231, and further that inhibition of miR-140-3p also increases cellular growth in MCF-7, but not MDA-MB-231. Broadly, these results suggest that the creation of miRNA targets may be an underappreciated function of ADAR and may help further elucidate the role of RNA editing in tumor pathogenicity.

Global profiling of hnRNP A2/B1-RNA binding on chromatin highlights LncRNA interactions.

Long noncoding RNAs (lncRNAs) often carry out their functions through associations with adaptor proteins. We recently identified heterogeneous ribonucleoprotein (hnRNP) A2/B1 as an adaptor of the human HOTAIR lncRNA. hnRNP A2 and B1 are splice isoforms of the same gene. The spliced version of HOTAIR preferentially associates with the B1 isoform, which we hypothesize contributes to RNA-RNA matching between HOTAIR and transcripts of target genes in breast cancer. Here we used enhanced cross-linking immunoprecipitation (eCLIP) to map the direct interactions between A2/B1 and RNA in breast cancer cells. Despite differing by only twelve amino acids, the A2 and B1 splice isoforms associate preferentially with distinct populations of RNA in vivo. Through cellular fractionation experiments we characterize the pattern of RNA association in chromatin, nucleoplasm, and cytoplasm. We find that a majority of interactions occur on chromatin, even those that do not contribute to co-transcriptional splicing. A2/B1 binding site locations on multiple RNAs hint at a contribution to the regulation and function of lncRNAs. Surprisingly, the strongest A2/B1 binding site occurs in a retained intron of HOTAIR, which interrupts an RNA-RNA interaction hotspot. In vitro eCLIP experiments highlight additional exonic B1 binding sites in HOTAIR which also surround the RNA-RNA interaction hotspot. Interestingly, a version of HOTAIR with the intron retained is still capable of making RNA-RNA interactions in vitro through the hotspot region. Our data further characterize the multiple functions of a repurposed splicing factor with isoform-biased interactions, and highlight that the majority of these functions occur on chromatin-associated RNA.

Human snoRNA-93 is processed into a microRNA-like RNA that promotes breast cancer cell invasion.

Genetic searches for tumor suppressors have recently linked small nucleolar RNA misregulations with tumorigenesis. In addition to their classically defined functions, several small nucleolar RNAs are now known to be processed into short microRNA-like fragments called small nucleolar RNA-derived RNAs. To determine if any small nucleolar RNA-derived RNAs contribute to breast malignancy, we recently performed a RNA-seq-based comparison of the small nucleolar RNA-derived RNAs of two breast cancer cell lines (MCF-7 and MDA-MB-231) and identified small nucleolar RNA-derived RNAs derived from 13 small nucleolar RNAs overexpressed in MDA-MB-231s. Importantly, we find that inhibiting the most differentially expressed of these small nucleolar RNA-derived RNAs (sdRNA-93) in MDA-MB-231 cells results primarily in a loss of invasiveness, whereas increased sdRNA-93 expression in either cell line conversely results in strikingly enhanced invasion. Excitingly, we recently determined sdRNA-93 expressions in small RNA-seq data corresponding to 116 patient tumors and normal breast controls, and while we find little sdRNA-93 expression in any of the controls and only sporadic expression in most subtypes, we find robust expression of sdRNA-93 in 92.8% of Luminal B Her2+tumors. Of note, our analyses also indicate that at least one of sdRNA-93's endogenous roles is to regulate the expression of Pipox, a sarcosine metabolism-related protein whose expression significantly correlates with distinct molecular subtypes of breast cancer. We find sdRNA-93 can regulate the Pipox 3'UTR via standard reporter assays and that manipulating endogenous sdRNA-93 levels inversely correlates with altered Pipox expression. In summary, our results strongly indicate that sdRNA-93 expression actively contributes to the malignant phenotype of breast cancer through participating in microRNA-like regulation.

Computational Prediction of MicroRNA Target Genes, Target Prediction Databases, and Web Resources.

MicroRNA (miRNA) mediated silencing and repression of mRNA molecules requires complementary base pairing between the "seed" region of the miRNA and the "seed match" region of target mRNAs. While this mechanism is fairly well understood, accurate prediction of valid miRNA targets remains challenging due to factors such as imperfect sequence specificity, target site availability, and the thermodynamic stability of the mRNA structure itself. As knowledge of what genes are being targeted by each miRNA is arguably the most important facet of miRNA biology, many approaches have been developed to address the need for reliable prediction and ranking of putative targets, with most using a combination of various strategies such as evolutionary conservation, statistical inference, and distinct features of the target sequences themselves. This chapter reviews the pros and cons of a number of different prediction algorithms, showcases some databases that store experimentally validated miRNA targets, and also provides a case study that profiles some of the potential microRNA-mRNA interactions predicted by each methodology for various human genes.

Novel small RNA (sRNA) landscape of the starvation-stress response transcriptome of Salmonella enterica serovar typhimurium.

Small RNAs (sRNAs) are short (∼50-200 nucleotides) noncoding RNAs that regulate cellular activities across bacteria. Salmonella enterica starved of a carbon-energy (C) source experience a host of genetic and physiological changes broadly referred to as the starvation-stress response (SSR). In an attempt to identify novel sRNAs contributing to SSR control, we grew log-phase, 5-h C-starved and 24-h C-starved cultures of the virulent Salmonella enterica subspecies enterica serovar Typhimurium strain SL1344 and comprehensively sequenced their small RNA transcriptomes. Strikingly, after employing a novel strategy for sRNA discovery based on identifying dynamic transcripts arising from "gene-empty" regions, we identify 58 wholly undescribed Salmonella sRNA genes potentially regulating SSR averaging an ∼1,000-fold change in expression between log-phase and C-starved cells. Importantly, the expressions of individual sRNA loci were confirmed by both comprehensive transcriptome analyses and northern blotting of select candidates. Of note, we find 43 candidate sRNAs share significant sequence identity to characterized sRNAs in other bacteria, and ∼70% of our sRNAs likely assume characteristic sRNA structural conformations. In addition, we find 53 of our 58 candidate sRNAs either overlap neighboring mRNA loci or share significant sequence complementarity to mRNAs transcribed elsewhere in the SL1344 genome strongly suggesting they regulate the expression of transcripts via antisense base-pairing. Finally, in addition to this work resulting in the identification of 58 entirely novel Salmonella enterica genes likely participating in the SSR, we also find evidence suggesting that sRNAs are significantly more prevalent than currently appreciated and that Salmonella sRNAs may actually number in the thousands.

An oxidative DNA "damage" and repair mechanism localized in the VEGF promoter is important for hypoxia-induced VEGF mRNA expression.

In hypoxia, mitochondria-generated reactive oxygen species not only stimulate accumulation of the transcriptional regulator of hypoxic gene expression, hypoxia inducible factor-1 (Hif-1), but also cause oxidative base modifications in hypoxic response elements (HREs) of hypoxia-inducible genes. When the hypoxia-induced base modifications are suppressed, Hif-1 fails to associate with the HRE of the VEGF promoter, and VEGF mRNA accumulation is blunted. The mechanism linking base modifications to transcription is unknown. Here we determined whether recruitment of base excision DNA repair (BER) enzymes in response to hypoxia-induced promoter modifications was required for transcription complex assembly and VEGF mRNA expression. Using chromatin immunoprecipitation analyses in pulmonary artery endothelial cells, we found that hypoxia-mediated formation of the base oxidation product 8-oxoguanine (8-oxoG) in VEGF HREs was temporally associated with binding of Hif-1α and the BER enzymes 8-oxoguanine glycosylase 1 (Ogg1) and redox effector factor-1 (Ref-1)/apurinic/apyrimidinic endonuclease 1 (Ape1) and introduction of DNA strand breaks. Hif-1α colocalized with HRE sequences harboring Ref-1/Ape1, but not Ogg1. Inhibition of BER by small interfering RNA-mediated reduction in Ogg1 augmented hypoxia-induced 8-oxoG accumulation and attenuated Hif-1α and Ref-1/Ape1 binding to VEGF HRE sequences and blunted VEGF mRNA expression. Chromatin immunoprecipitation-sequence analysis of 8-oxoG distribution in hypoxic pulmonary artery endothelial cells showed that most of the oxidized base was localized to promoters with virtually no overlap between normoxic and hypoxic data sets. Transcription of genes whose promoters lost 8-oxoG during hypoxia was reduced, while those gaining 8-oxoG was elevated. Collectively, these findings suggest that the BER pathway links hypoxia-induced introduction of oxidative DNA modifications in promoters of hypoxia-inducible genes to transcriptional activation.

Burgeoning evidence indicates that microRNAs were initially formed from transposable element sequences.

MicroRNAs (miRNAs) constitute a recently discovered class of noncoding RNAs that play key roles in the regulation of gene expression. Despite being only ~20 nucleotides in length, these highly versatile molecules have been shown to play pivotal roles in development, basic cellular metabolism, apoptosis, and disease. While over 24,000 miRNAs have been characterized since they were first isolated in mammals in 2001, the functions of the majority of these miRNAs remain largely undescribed. That said, many now suggest that characterization of the relationships between miRNAs and transposable elements (TEs) can help elucidate miRNA functionality. Strikingly, over 20 publications have now reported the initial formation of thousands of miRNA loci from TE sequences. In this review we chronicle the findings of these reports, discuss the evolution of the field along with future directions, and examine how this information can be used to ascertain insights into miRNA transcriptional regulation and how it can be exploited to facilitate miRNA target prediction.

Continuing analysis of microRNA origins: Formation from transposable element insertions and noncoding RNA mutations.

MicroRNAs (miRs) are small noncoding RNAs that typically act as regulators of gene expression by base pairing with the 3' UTR of messenger RNAs (mRNAs) and either repressing their translation or initiating degradation. As of this writing over 24,500 distinct miRs have been identified, but the functions of the vast majority of these remain undescribed. This paper represents a summary of our in depth analysis of the genomic origins of miR loci, detailing the formation of 1,213 of the 7,321 recently identified miRs and thereby bringing the total number of miR loci with defined molecular origin to 3,605. Interestingly, our analyses also identify evidence for a second, novel mechanism of miR locus generation through describing the formation of 273 miR loci from mutations to other forms of noncoding RNAs. Importantly, several independent investigations of the genomic origins of miR loci have now supported the hypothesis that miR hairpins are formed by the adjacent genomic insertion of two complementary transposable elements (TEs) into opposing strands. While our results agree that subsequent transcription over such TE interfaces leads to the formation of the majority of functional miR loci, we now also find evidence suggesting that a subset of miR loci were actually formed by an alternative mechanism-point mutations in other structurally complex, noncoding RNAs (e.g., tRNAs and snoRNAs).