Statistical Assessment of Depth Normalization for Small RNA Sequencing.

PURPOSE: Methods for depth normalization have been assessed primarily with simulated data or cell-line-mixture data. There is a pressing need for benchmark data enabling a more realistic and objective assessment, especially in the context of small RNA sequencing. METHODS: We collected a unique pair of microRNA sequencing data sets for the same set of tumor samples; one data set was collected with and the other without uniform handling and balanced design. The former provided a benchmark for evaluating evidence of differential expression and the latter served as a test bed for normalization. Next, we developed a data perturbation algorithm to simulate additional data set pairs. Last, we assembled a set of computational tools to visualize and quantify the assessment. RESULTS: We validated the quality of the benchmark data and showed the need for normalization of the test data. For illustration, we applied the data and tools to assess the performance of 9 existing normalization methods. Among them, trimmed mean of M-values was a better scaling method, whereas the median and the upper quartiles were consistently the worst performers; one variation of remove unwanted variation had the best chance of capturing true positives but at the cost of increased false positives. In general, these methods were, at best, moderately helpful when the level of differential expression was extensive and asymmetric. CONCLUSION: Our study (1) provides the much-needed benchmark data and computational tools for assessing depth normalization, (2) shows the dependence of normalization performance on the underlying pattern of differential expression, and (3) calls for continued research efforts to develop more effective normalization methods.

The transcription factor ERG recruits CCR4-NOT to control mRNA decay and mitotic progression.

Control of mRNA levels, a fundamental aspect in the regulation of gene expression, is achieved through a balance between mRNA synthesis and decay. E26-related gene (Erg) proteins are canonical transcription factors whose previously described functions are confined to the control of mRNA synthesis. Here, we report that ERG also regulates gene expression by affecting mRNA stability and identify the molecular mechanisms underlying this function in human cells. ERG is recruited to mRNAs via interaction with the RNA-binding protein RBPMS, and it promotes mRNA decay by binding CNOT2, a component of the CCR4-NOT deadenylation complex. Transcriptome-wide mRNA stability analysis revealed that ERG controls the degradation of a subset of mRNAs highly connected to Aurora signaling, whose decay during S phase is necessary for mitotic progression. Our data indicate that control of gene expression by mammalian transcription factors may follow a more complex scheme than previously anticipated, integrating mRNA synthesis and degradation.

Structural basis underlying CAC RNA recognition by the RRM domain of dimeric RNA-binding protein RBPMS.

RNA-binding protein with multiple splicing (designated RBPMS) is a higher vertebrate mRNA-binding protein containing a single RNA recognition motif (RRM). RBPMS has been shown to be involved in mRNA transport, localization and stability, with key roles in axon guidance, smooth muscle plasticity, as well as regulation of cancer cell proliferation and migration. We report on structure-function studies of the RRM domain of RBPMS bound to a CAC-containing single-stranded RNA. These results provide insights into potential topologies of complexes formed by the RBPMS RRM domain and the tandem CAC repeat binding sites as detected by photoactivatable-ribonucleoside-enhanced crosslinking and immunoprecipitation. These studies establish that the RRM domain of RBPMS forms a symmetrical dimer in the free state, with each monomer binding sequence-specifically to all three nucleotides of a CAC segment in the RNA bound state. Structure-guided mutations within the dimerization and RNA-binding interfaces of RBPMS RRM on RNA complex formation resulted in both disruption of dimerization and a decrease in RNA-binding affinity as observed by size exclusion chromatography and isothermal titration calorimetry. As anticipated from biochemical binding studies, over-expression of dimerization or RNA-binding mutants of Flag-HA-tagged RBPMS were no longer able to track with stress granules in HEK293 cells, thereby documenting the deleterious effects of such mutations in vivo.

Deep MicroRNA sequencing reveals downregulation of miR-29a in neuroblastoma central nervous system metastasis.

Central nervous system (CNS) is an increasingly common site of isolated metastasis for patients with Stage 4 neuroblastoma. To explore the microRNA (miRNA) profile of this metastatic process, miRNA sequencing was performed to identify miRNA sequence families with differential expression between tumor pairs (pre-CNS primary and CNS metastasis) from 13 patients with Stage 4 neuroblastoma. Seven miRNA sequence families had distinct expression in CNS metastases when compared with their corresponding pre-CNS primaries. MiR-7 was upregulated (3.75-fold), and miR-21, miR-22, miR-29a, miR-143, miR-199a-1-3p, and miR-199a-1-5p were downregulated (3.5-6.1-fold), all confirmed by quantitative reverse transcription-PCR. MiR-29a, previously shown to be downregulated in a broad spectrum of solid tumors including neuroblastoma, had the most significant decrease in all 13 CNS metastases (P = 0.001). Its known onco-targets CDC6, CDK6, and DNMT3A, as well as B7-H3, an inhibitory ligand for T cells, and natural killer cells, were found to have higher differential expression in these 13 CNS metastases when compared with their paired primaries. Additionally, miR-29a expression in primary tumors was significantly lower among patients who eventually relapsed in the CNS. Irrespective of the amplification status of MYCN, which is known to be associated with metastasis, pre-CNS primaries, and CNS metastases had significantly lower miR-29a expression than non-CNS primary tumors. Among MYCN amplified cell lines, those from CNS relapse also had lower miR-29a expression than non-CNS relapse. These findings raised the hypothesis that miR-29a could be a biomarker for neuroblastoma CNS metastasis, and its downregulation may play a pivotal role in CNS progression.

Identification of the RNA recognition element of the RBPMS family of RNA-binding proteins and their transcriptome-wide mRNA targets.

Recent studies implicated the RNA-binding protein with multiple splicing (RBPMS) family of proteins in oocyte, retinal ganglion cell, heart, and gastrointestinal smooth muscle development. These RNA-binding proteins contain a single RNA recognition motif (RRM), and their targets and molecular function have not yet been identified. We defined transcriptome-wide RNA targets using photoactivatable-ribonucleoside-enhanced crosslinking and immunoprecipitation (PAR-CLIP) in HEK293 cells, revealing exonic mature and intronic pre-mRNA binding sites, in agreement with the nuclear and cytoplasmic localization of the proteins. Computational and biochemical approaches defined the RNA recognition element (RRE) as a tandem CAC trinucleotide motif separated by a variable spacer region. Similar to other mRNA-binding proteins, RBPMS family of proteins relocalized to cytoplasmic stress granules under oxidative stress conditions suggestive of a support function for mRNA localization in large and/or multinucleated cells where it is preferentially expressed.

Human CLP1 mutations alter tRNA biogenesis, affecting both peripheral and central nervous system function.

CLP1 is a RNA kinase involved in tRNA splicing. Recently, CLP1 kinase-dead mice were shown to display a neuromuscular disorder with loss of motor neurons and muscle paralysis. Human genome analyses now identified a CLP1 homozygous missense mutation (p.R140H) in five unrelated families, leading to a loss of CLP1 interaction with the tRNA splicing endonuclease (TSEN) complex, largely reduced pre-tRNA cleavage activity, and accumulation of linear tRNA introns. The affected individuals develop severe motor-sensory defects, cortical dysgenesis, and microcephaly. Mice carrying kinase-dead CLP1 also displayed microcephaly and reduced cortical brain volume due to the enhanced cell death of neuronal progenitors that is associated with reduced numbers of cortical neurons. Our data elucidate a neurological syndrome defined by CLP1 mutations that impair tRNA splicing. Reduction of a founder mutation to homozygosity illustrates the importance of rare variations in disease and supports the clan genomics hypothesis.

Identification of distinct miRNA target regulation between breast cancer molecular subtypes using AGO2-PAR-CLIP and patient datasets.

BACKGROUND: Various microRNAs (miRNAs) are up- or downregulated in tumors. However, the repression of cognate miRNA targets responsible for the phenotypic effects of this dysregulation in patients remains largely unexplored. To define miRNA targets and associated pathways, together with their relationship to outcome in breast cancer, we integrated patient-paired miRNA-mRNA expression data with a set of validated miRNA targets and pathway inference. RESULTS: To generate a biochemically-validated set of miRNA-binding sites, we performed argonaute-2 photoactivatable-ribonucleoside-enhanced crosslinking and immunoprecipitation (AGO2-PAR-CLIP) in MCF7 cells. We then defined putative miRNA-target interactions using a computational model, which ranked and selected additional TargetScan-predicted interactions based on features of our AGO2-PAR-CLIP binding-site data. We subselected modeled interactions according to the abundance of their constituent miRNA and mRNA transcripts in tumors, and we took advantage of the variability of miRNA expression within molecular subtypes to detect miRNA repression. Interestingly, our data suggest that miRNA families control subtype-specific pathways; for example, miR-17, miR-19a, miR-25, and miR-200b show high miRNA regulatory activity in the triple-negative, basal-like subtype, whereas miR-22 and miR-24 do so in the HER2 subtype. An independent dataset validated our findings for miR-17 and miR-25, and showed a correlation between the expression levels of miR-182 targets and overall patient survival. Pathway analysis associated miR-17, miR-19a, and miR-200b with leukocyte transendothelial migration. CONCLUSIONS: We combined PAR-CLIP data with patient expression data to predict regulatory miRNAs, revealing potential therapeutic targets and prognostic markers in breast cancer.

Mammalian miRNA curation through next-generation sequencing.

Characteristic small RNA biogenesis processing patterns are used for the discovery of novel microRNAs (miRNAs) from next-generation sequencing data. Here, we highlight and discuss key criteria for mammalian - specifically human - miRNA database curation based on small RNA sequencing data. Sequence reads obtained from small RNA cDNA libraries are aligned to reference genomic regions, and miRNA genes are revealed by their distinct read length and bimodal read frequency distribution, the predicted secondary structure of the deduced miRNA stem-loop precursor molecule, and, to a lesser degree, based on evolutionary conservation of small RNAs from other vertebrates. Properly curated miRNA databases are an important resource for investigators interested in miRNA biology, diagnostics, and therapeutics.

MicroRNAs in human cancer.

Mature microRNAs (miRNAs) are single-stranded RNA molecules of 20-23-nucleotide (nt) length that control gene expression in many cellular processes. These molecules typically reduce the translation and stability of mRNAs, including those of genes that mediate processes in tumorigenesis, such as inflammation, cell cycle regulation, stress response, differentiation, apoptosis, and invasion. miRNA targeting is initiated through specific base-pairing interactions between the 5' end ("seed" region) of the miRNA and sites within coding and untranslated regions (UTRs) of mRNAs; target sites in the 3' UTR lead to more effective mRNA destabilization. Since miRNAs frequently target hundreds of mRNAs, miRNA regulatory pathways are complex. To provide a critical overview of miRNA dysregulation in cancer, we first discuss the methods currently available for studying the role of miRNAs in cancer and then review miRNA genomic organization, biogenesis, and mechanism of target recognition, examining how these processes are altered in tumorigenesis. Given the critical role miRNAs play in tumorigenesis processes and their disease specific expression, they hold potential as therapeutic targets and novel biomarkers.

Barcoded cDNA library preparation for small RNA profiling by next-generation sequencing.

The characterization of post-transcriptional gene regulation by small regulatory (20-30 nt) RNAs, particularly miRNAs and piRNAs, has become a major focus of research in recent years. A prerequisite for characterizing small RNAs is their identification and quantification across different developmental stages, and in normal and disease tissues, as well as model cell lines. Here we present a step-by-step protocol for generating barcoded small RNA cDNA libraries compatible with Illumina HiSeq sequencing, thereby facilitating miRNA and other small RNA profiling of large sample collections.

Bioinformatic analysis of barcoded cDNA libraries for small RNA profiling by next-generation sequencing.

The characterization of post-transcriptional gene regulation by small regulatory RNAs of 20-30 nt length, particularly miRNAs and piRNAs, has become a major focus of research in recent years. A prerequisite for the characterization of small RNAs is their identification and quantification across different developmental stages, normal and diseased tissues, as well as model cell lines. Here we present a step-by-step protocol for the bioinformatic analysis of barcoded cDNA libraries for small RNA profiling generated by Illumina sequencing, thereby facilitating miRNA and other small RNA profiling of large sample collections.

RNA targets of wild-type and mutant FET family proteins.

FUS, EWSR1 and TAF15, constituting the FET protein family, are abundant, highly conserved RNA-binding proteins with important roles in oncogenesis and neuronal disease, yet their RNA targets and recognition elements are unknown. Using PAR-CLIP, we defined global RNA targets for all human FET proteins and two ALS-causing human FUS mutants. FET members showed similar binding profiles, whereas FUS mutants showed a drastically altered binding pattern, consistent with changes in subcellular localization.

MicroRNA sequence and expression analysis in breast tumors by deep sequencing.

MicroRNAs (miRNA) regulate many genes critical for tumorigenesis. We profiled miRNAs from 11 normal breast tissues, 17 noninvasive, 151 invasive breast carcinomas, and 6 cell lines by in-house-developed barcoded Solexa sequencing. miRNAs were organized in genomic clusters representing promoter-controlled miRNA expression and sequence families representing seed sequence-dependent miRNA target regulation. Unsupervised clustering of samples by miRNA sequence families best reflected the clustering based on mRNA expression available for this sample set. Clustering and comparative analysis of miRNA read frequencies showed that normal breast samples were separated from most noninvasive ductal carcinoma in situ and invasive carcinomas by increased miR-21 (the most abundant miRNA in carcinomas) and multiple decreased miRNA families (including miR-98/let-7), with most miRNA changes apparent already in the noninvasive carcinomas. In addition, patients that went on to develop metastasis showed increased expression of mir-423, and triple-negative breast carcinomas were most distinct from other tumor subtypes due to upregulation of the mir~17-92 cluster. However, absolute miRNA levels between normal breast and carcinomas did not reveal any significant differences. We also discovered two polymorphic nucleotide variations among the more abundant miRNAs miR-181a (T19G) and miR-185 (T16G), but we did not identify nucleotide variations expected for classical tumor suppressor function associated with miRNAs. The differentiation of tumor subtypes and prediction of metastasis based on miRNA levels is statistically possible but is not driven by deregulation of abundant miRNAs, implicating far fewer miRNAs in tumorigenic processes than previously suggested.

miRNAs in human cancer.

Mature microRNAs (miRNAs) are single-stranded RNA molecules of 20-23 nucleotide (nt) length that control gene expression in many cellular processes. These molecules typically reduce the stability of mRNAs, including those of genes that mediate processes in tumorigenesis, such as inflammation, cell cycle regulation, stress response, differentiation, apoptosis and invasion. miRNA targeting is mostly achieved through specific base-pairing interactions between the 5' end ('seed' region) of the miRNA and sites within coding and untranslated regions (UTRs) of mRNAs; target sites in the 3' UTR lead to more effective mRNA destabilization. Since miRNAs frequently target hundreds of mRNAs, miRNA regulatory pathways are complex. To provide a critical overview of miRNA dysregulation in cancer, we first discuss the methods currently available for studying the role of miRNAs in cancer and then review miRNA genomic organization, biogenesis and mechanism of target recognition, examining how these processes are altered in tumorigenesis. Given the critical role miRNAs play in tumorigenesis processes and their disease-specific expression, they hold potential as therapeutic targets and novel biomarkers.

The growing catalog of small RNAs and their association with distinct Argonaute/Piwi family members.

Several distinct classes of small RNAs, some newly identified, have been discovered to play important regulatory roles in diverse cellular processes. These classes include siRNAs, miRNAs, rasiRNAs and piRNAs. Each class binds to distinct members of the Argonaute/Piwi protein family to form ribonucleoprotein complexes that recognize partially, or nearly perfect, complementary nucleic acid targets, and that mediate a variety of regulatory processes, including transcriptional and post-transcriptional gene silencing. Based on the known relationship of Argonaute/Piwi proteins with distinct classes of small RNAs, we can now predict how many new classes of small RNAs or silencing processes remain to be discovered.