Small RNA populations revealed by blocking rRNA fragments in Drosophila melanogaster reproductive tissues

RNA interference (RNAi) is a complex and highly conserved regulatory mechanism mediated via small RNAs (sRNAs). Recent technical advances in high throughput sequencing have enabled an increasingly detailed analysis of sRNA abundances and profiles in specific body parts and tissues. This enables investigations of the localized roles of microRNAs (miRNAs) and small interfering RNAs (siRNAs). However, variation in the proportions of non-coding RNAs in the samples being compared can hinder these analyses. Specific tissues may vary significantly in the proportions of fragments of longer non-coding RNAs (such as ribosomal RNA or transfer RNA) present, potentially reflecting tissue-specific differences in biological functions. For example, in Drosophila, some tissues contain a highly abundant 30nt rRNA fragment (the 2S rRNA) as well as abundant 5’ and 3’ terminal rRNA fragments. These can pose difficulties for the construction of sRNA libraries as they can swamp the sequencing space and obscure sRNA abundances. Here we addressed this problem and present a modified “rRNA blocking” protocol for the construction of high-definition (HD) adapter sRNA libraries, in D. melanogaster reproductive tissues. The results showed that 2S rRNAs targeted by blocking oligos were reduced from >80% to < 0.01% total reads. In addition, the use of multiple rRNA blocking oligos to bind the most abundant rRNA fragments allowed us to reveal the underlying sRNA populations at increased resolution. Side-by-side comparisons of sequencing libraries of blocked and non-blocked samples revealed that rRNA blocking did not change the miRNA populations present, but instead enhanced their abundances. We suggest that this rRNA blocking procedure offers the potential to improve the in-depth analysis of differentially expressed sRNAs within and across different tissues

MysiRNA-designer: a workflow for efficient siRNA design

The design of small interfering RNA (siRNA) is a multi factorial problem that has gained the attention of many researchers in the area of therapeutic and functional genomics. MysiRNA score was previously introduced that improves the correlation of siRNA activity prediction considering state of the art algorithms. In this paper, a new program, MysiRNA-Designer, is described which integrates several factors in an automated work-flow considering mRNA transcripts variations, siRNA and mRNA target accessibility, and both near-perfect and partial off-target matches. It also features the MysiRNA score, a highly ranked correlated siRNA efficacy prediction score for ranking the designed siRNAs, in addition to top scoring models Biopredsi, DISR, Thermocomposition21 and i-Score, and integrates them in a unique siRNA score-filtration technique. This multi-score filtration layer filters siRNA that passes the 90% thresholds calculated from experimental dataset features. MysiRNA-Designer takes an accession, finds conserved regions among its transcript space, finds accessible regions within the mRNA, designs all possible siRNAs for these regions, filters them based on multi-scores thresholds, and then performs SNP and off-target filtration. These strict selection criteria were tested against human genes in which at least one active siRNA was designed from 95.7% of total genes. In addition, when tested against an experimental dataset, MysiRNA-Designer was found capable of rejecting 98% of the false positive siRNAs, showing superiority over three state of the art siRNA design programs. MysiRNA is a freely accessible (Microsoft Windows based) desktop application that can be used to design siRNA with a high accuracy and specificity. We believe that MysiRNA-Designer has the potential to play an important role in this area

Quantitative Prediction of miRNA-mRNA Interaction Based on Equilibrium Concentrations

MicroRNAs (miRNAs) suppress gene expression by forming a duplex with a target messenger RNA (mRNA), blocking translation or initiating cleavage. Computational approaches have proven valuable for predicting which mRNAs can be targeted by a given miRNA, but currently available prediction methods do not address the extent of duplex formation under physiological conditions. Some miRNAs can at low concentrations bind to target mRNAs, whereas others are unlikely to bind within a physiologically relevant concentration range. Here we present a novel approach in which we find potential target sites on mRNA that minimize the calculated free energy of duplex formation, compute the free energy change involved in unfolding these sites, and use these energies to estimate the extent of duplex formation at specified initial concentrations of both species. We compare our predictions to experimentally confirmed miRNA-mRNA interactions (and non-interactions) in Drosophila melanogaster and in human. Although our method does not predict whether the targeted mRNA is degraded and/or its translation to protein inhibited, our quantitative estimates generally track experimentally supported results, indicating that this approach can be used to predict whether an interaction occurs at specified concentrations. Our approach offers a more-quantitative understanding of post-translational regulation in different cell types, tissues, and developmental conditions

The Caenorhabditis elegans HEN1 Ortholog, HENN-1, Methylates and Stabilizes Select Subclasses of Germline Small RNAs

Small RNAs regulate diverse biological processes by directing effector proteins called Argonautes to silence complementary mRNAs. Maturation of some classes of small RNAs involves terminal 2′-O-methylation to prevent degradation. This modification is catalyzed by members of the conserved HEN1 RNA methyltransferase family. In animals, Piwi-interacting RNAs (piRNAs) and some endogenous and exogenous small interfering RNAs (siRNAs) are methylated, whereas microRNAs are not. However, the mechanisms that determine animal HEN1 substrate specificity have yet to be fully resolved. In Caenorhabditis elegans, a HEN1 ortholog has not been studied, but there is evidence for methylation of piRNAs and some endogenous siRNAs. Here, we report that the worm HEN1 ortholog, HENN-1 (HEN of Nematode), is required for methylation of C. elegans small RNAs. Our results indicate that piRNAs are universally methylated by HENN-1. In contrast, 26G RNAs, a class of primary endogenous siRNAs, are methylated in female germline and embryo, but not in male germline. Intriguingly, the methylation pattern of 26G RNAs correlates with the expression of distinct male and female germline Argonautes. Moreover, loss of the female germline Argonaute results in loss of 26G RNA methylation altogether. These findings support a model wherein methylation status of a metazoan small RNA is dictated by the Argonaute to which it binds. Loss of henn-1 results in phenotypes that reflect destabilization of substrate small RNAs: dysregulation of target mRNAs, impaired fertility, and enhanced somatic RNAi. Additionally, the henn-1 mutant shows a weakened response to RNAi knockdown of germline genes, suggesting that HENN-1 may also function in canonical RNAi. Together, our results indicate a broad role for HENN-1 in both endogenous and exogenous gene silencing pathways and provide further insight into the mechanisms of HEN1 substrate discrimination and the diversity within the Argonaute family

MicroRNA interactome analysis predicts post-transcriptional regulation of ADRB2 and PPP3R1 in the hypercholesterolemic myocardium

Little is known about the molecular mechanism including microRNAs (miRNA) in hypercholesterolemia-induced cardiac dysfunction. We aimed to explore novel hypercholesterolemia-induced pathway alterations in the heart by an unbiased approach based on miRNA omics, target prediction and validation. With miRNA microarray we identified forty-seven upregulated and ten downregulated miRNAs in hypercholesterolemic rat hearts compared to the normocholesterolemic group. Eleven mRNAs with at least 4 interacting upregulated miRNAs were selected by a network theoretical approach, out of which 3 mRNAs (beta-2 adrenergic receptor [Adrb2], calcineurin B type 1 [Ppp3r1] and calcium/calmodulin-dependent serine protein kinase [Cask]) were validated with qRT-PCR and Western blot. In hypercholesterolemic hearts, the expression of Adrb2 mRNA was significantly decreased. ADRB2 and PPP3R1 protein were significantly downregulated in hypercholesterolemic hearts. The direct interaction of Adrb2 with upregulated miRNAs was demonstrated by luciferase reporter assay. Gene ontology analysis revealed that the majority of the predicted mRNA changes may contribute to the hypercholesterolemia-induced cardiac dysfunction. In summary, the present unbiased target prediction approach based on global cardiac miRNA expression profiling revealed for the first time in the literature that both the mRNA and protein product of Adrb2 and PPP3R1 protein are decreased in the hypercholesterolemic heart

Characterization and mitigation of gene expression burden in mammalian cells

Despite recent advances in circuit engineering, the design of genetic networks in mammalian cells is still painstakingly slow and fraught with inexplicable failures. Here, we demonstrate that transiently expressed genes in mammalian cells compete for limited transcriptional and translational resources. This competition results in the coupling of otherwise independent exogenous and endogenous genes, creating a divergence between intended and actual function. Guided by a resource-aware mathematical model, we identify and engineer natural and synthetic miRNA-based incoherent feedforward loop (iFFL) circuits that mitigate gene expression burden. The implementation of these circuits features the use of endogenous miRNAs as elementary components of the engineered iFFL device, a versatile hybrid design that allows burden mitigation to be achieved across different cell-lines with minimal resource requirements. This study establishes the foundations for context-aware prediction and improvement of in vivo synthetic circuit performance, paving the way towards more rational synthetic construct design in mammalian cells

Quantitative analysis of MicroRNAs in vaccinia virus infection reveals diversity in their susceptibility to modification and suppression

Vaccinia virus (VACV) is a large cytoplasmic DNA virus that causes dramatic alterations to many cellular pathways including microRNA biogenesis. The virus encodes a poly(A) polymerase which was previously shown to add poly(A) tails to the 3' end of cellular miRNAs, resulting in their degradation by 24 hours post infection (hpi). Here we used small RNA sequencing to quantify the impact of VACV infection on cellular miRNAs in human cells at both early (6 h) and late (24 h) times post infection. A detailed quantitative analysis of individual miRNAs revealed marked diversity in the extent of their modification and relative change in abundance during infection. Some miRNAs became highly modified (e.g. miR-29a-3p, miR-27b-3p) whereas others appeared resistant (e.g. miR-16-5p). Furthermore, miRNAs that were highly tailed at 6 hpi were not necessarily among the most reduced at 24 hpi. These results suggest that intrinsic features of human cellular miRNAs cause them to be differentially polyadenylated and altered in abundance during VACV infection. We also demonstrate that intermediate and late VACV gene expression are required for optimal repression of some miRNAs including miR-27-3p. Overall this work reveals complex and varied consequences of VACV infection on host miRNAs and identifies miRNAs which are largely resistant to VACV-induced polyadenylation and are therefore present at functional levels during the initial stages of infection and replication

PIWI Associated siRNAs and piRNAs Specifically Require the Caenorhabditis elegans HEN1 Ortholog henn-1

Small RNAs—including piRNAs, miRNAs, and endogenous siRNAs—bind Argonaute proteins to form RNA silencing complexes that target coding genes, transposons, and aberrant RNAs. To assess the requirements for endogenous siRNA formation and activity in Caenorhabditis elegans, we developed a GFP-based sensor for the endogenous siRNA 22G siR-1, one of a set of abundant siRNAs processed from a precursor RNA mapping to the X chromosome, the X-cluster. Silencing of the sensor is also dependent on the partially complementary, unlinked 26G siR-O7 siRNA. We show that 26G siR-O7 acts in trans to initiate 22G siRNA formation from the X-cluster. The presence of several mispairs between 26G siR-O7 and the X-cluster mRNA, as well as mutagenesis of the siRNA sensor, indicates that siRNA target recognition is permissive to a degree of mispairing. From a candidate reverse genetic screen, we identified several factors required for 22G siR-1 activity, including the chromatin factors mes-4 and gfl-1, the Argonaute ergo-1, and the 3′ methyltransferase henn-1. Quantitative RT–PCR of small RNAs in a henn-1 mutant and deep sequencing of methylated small RNAs indicate that siRNAs and piRNAs that associate with PIWI clade Argonautes are methylated by HENN-1, while siRNAs and miRNAs that associate with non-PIWI clade Argonautes are not. Thus, PIWI-class Argonaute proteins are specifically adapted to associate with methylated small RNAs in C. elegans

Mechanism of MicroRNA-Target Interaction: Molecular Dynamics Simulations and Thermodynamics Analysis

MicroRNAs (miRNAs) are endogenously produced ∼21-nt riboregulators that associate with Argonaute (Ago) proteins to direct mRNA cleavage or repress the translation of complementary RNAs. Capturing the molecular mechanisms of miRNA interacting with its target will not only reinforce the understanding of underlying RNA interference but also fuel the design of more effective small-interfering RNA strands. To address this, in the present work the RNA-bound (Ago-miRNA, Ago-miRNA-target) and RNA-free Ago forms were analyzed by performing both molecular dynamics simulations and thermodynamic analysis. Based on the principal component analysis results of the simulation trajectories as well as the correlation analysis in fluctuations of residues, we discover that: 1) three important (PAZ, Mid and PIWI) domains exist in Argonaute which define the global dynamics of the protein; 2) the interdomain correlated movements are so crucial for the interaction of Ago-RNAs that they not only facilitate the relaxation of the interactions between residues surrounding the RNA binding channel but also induce certain conformational changes; and 3) it is just these conformational changes that expand the cavity of the active site and open putative pathways for both the substrate uptake and product release. In addition, by thermodynamic analysis we also discover that for both the guide RNA 5′-end recognition and the facilitated site-specific cleavage of the target, the presence of two metal ions (of Mg2+) plays a predominant role, and this conclusion is consistent with the observed enzyme catalytic cleavage activity in the ternary complex (Ago-miRNA-mRNA). Our results find that it is the set of arginine amino acids concentrated in the nucleotide-binding channel in Ago, instead of the conventionally-deemed seed base-paring, that makes greater contributions in stabilizing the binding of the nucleic acids to Ago

Identification and HLA-Tetramer-Validation of Human CD4(+) and CD8(+) T Cell Responses against HCMV Proteins IE1 and IE2

Human cytomegalovirus (HCMV) is an important human pathogen. It is a leading cause of congenital infection and a leading infectious threat to recipients of solid organ transplants as well as of allogeneic hematopoietic cell transplants. Moreover, it has recently been suggested that HCMV may promote tumor development. Both CD4+ and CD8+ T cell responses are important for long-term control of the virus, and adoptive transfer of HCMV-specific T cells has led to protection from reactivation and HCMV disease. Identification of HCMV-specific T cell epitopes has primarily focused on CD8+ T cell responses against the pp65 phosphoprotein. In this study, we have focused on CD4+ and CD8+ T cell responses against the immediate early 1 and 2 proteins (IE1 and IE2). Using overlapping peptides spanning the entire IE1 and IE2 sequences, peripheral blood mononuclear cells from 16 healthy, HLA-typed, donors were screened by ex vivo IFN-γ ELISpot and in vitro intracellular cytokine secretion assays. The specificities of CD4+ and CD8+ T cell responses were identified and validated by HLA class II and I tetramers, respectively. Eighty-one CD4+ and 44 CD8+ T cell responses were identified representing at least seven different CD4 epitopes and 14 CD8 epitopes restricted by seven and 11 different HLA class II and I molecules, respectively, in total covering 91 and 98% of the Caucasian population, respectively. Presented in the context of several different HLA class II molecules, two epitope areas in IE1 and IE2 were recognized in about half of the analyzed donors. These data may be used to design a versatile anti-HCMV vaccine and/or immunotherapy strategy