Systemic CLIP-seq analysis and game theory approach to model microRNA mode of binding.

microRNAs (miRNAs) associate with Ago proteins to post-transcriptionally silence gene expression by targeting mRNAs. To characterize the modes of miRNA-binding, we developed a novel computational framework, called optiCLIP, which considers the reproducibility of the identified peaks among replicates based on the peak overlap. We identified 98 999 binding sites for mouse and human miRNAs, from eleven Ago2 CLIP-seq datasets. Clustering the binding preferences, we found heterogeneity of the mode of binding for different miRNAs. Finally, we set up a quantitative model, named miRgame, based on an adaptation of the game theory. We have developed a new algorithm to translate the miRgame into a score that corresponds to a miRNA degree of occupancy for each Ago2 peak. The degree of occupancy summarizes the number of miRNA-binding sites and miRNAs targeting each binding site, and binding energy of each miRNA::RNA heteroduplex in each peak. Ago peaks were stratified accordingly to the degree of occupancy. Target repression correlates with higher score of degree of occupancy and number of miRNA-binding sites within each Ago peak. We validated the biological performance of our new method on miR-155-5p. In conclusion, our data demonstrate that miRNA-binding sites within each Ago2 CLIP-seq peak synergistically interplay to enhance target repression.

Subcellular Heterogeneity of the microRNA Machinery.

Different methods have recently been developed to understand the subcellular localization and role of microRNAs (miRNAs) as well as small RNAs associated with Argonaute (AGO) proteins. The heterogeneity of the protein complexes associated with miRNAs, along with their subcellular localization, provides clues into their biochemical mechanism of function. Subcellular diversity indicates that miRNAs localized to different cellular regions could have different functions, including transcriptional regulation on chromatin or post-transcriptional control, providing global regulation of gene expression by miRNAs. Herein, I review the current knowledge and most recent discoveries relating to the subcellular function of miRNAs and other AGO-associated small RNAs, revealing the emergence of a multitude of functions of the miRNA pathway to control different steps of the gene expression program(s).

Recent computational developments on CLIP-seq data analysis and microRNA targeting implications.

Cross-Linking Immunoprecipitation associated to high-throughput sequencing (CLIP-seq) is a technique used to identify RNA directly bound to RNA-binding proteins across the entire transcriptome in cell or tissue samples. Recent technological and computational advances permit the analysis of many CLIP-seq samples simultaneously, allowing us to reveal the comprehensive network of RNA-protein interaction and to integrate it to other genome-wide analyses. Therefore, the design and quality management of the CLIP-seq analyses are of critical importance to extract clean and biological meaningful information from CLIP-seq experiments. The application of CLIP-seq technique to Argonaute 2 (Ago2) protein, the main component of the microRNA (miRNA)-induced silencing complex, reveals the direct binding sites of miRNAs, thus providing insightful information about the role played by miRNA(s). In this review, we summarize and discuss the most recent computational methods for CLIP-seq analysis, and discuss their impact on Ago2/miRNA-binding site identification and prediction with a regard toward human pathologies.

Regulation of stimulus-inducible gene expression in myeloid cells.

One of the best-characterized and biologically important gene expression programmes in myeloid cells is their response to pro-inflammatory stimuli. Macrophages and DCs in particular are key mediators of immune responses, and are widely-used as prototypes to understand and define the determinants of specific and inducible gene expression. In this review we summarize advances and concepts which have been made towards the understanding of inducible gene expression, with a particular focus on insights gained using the myeloid system as a model. We discuss the emerging concept of layered control of gene regulation and cell identity by different functional classes of transcription factors; and examine recent progress to understanding the molecular processes involved, including the involvement of nucleosome positioning, chromatin modifications, and nuclear architecture. We also address the exciting but less-well understood role of non-coding RNAs in controlling specific gene expression programmes in myeloid and other cell-types.

From benchmarking HITS-CLIP peak detection programs to a new method for identification of miRNA-binding sites from Ago2-CLIP data.

Experimental evidence indicates that about 60% of miRNA-binding activity does not follow the canonical rule about the seed matching between miRNA and target mRNAs, but rather a non-canonical miRNA targeting activity outside the seed or with a seed-like motifs. Here, we propose a new unbiased method to identify canonical and non-canonical miRNA-binding sites from peaks identified by Ago2 Cross-Linked ImmunoPrecipitation associated to high-throughput sequencing (CLIP-seq). Since the quality of peaks is of pivotal importance for the final output of the proposed method, we provide a comprehensive benchmarking of four peak detection programs, namely CIMS, PIPE-CLIP, Piranha and Pyicoclip, on four publicly available Ago2-HITS-CLIP datasets and one unpublished in-house Ago2-dataset in stem cells. We measured the sensitivity, the specificity and the position accuracy toward miRNA binding sites identification, and the agreement with TargetScan. Secondly, we developed a new pipeline, called miRBShunter, to identify canonical and non-canonical miRNA-binding sites based on de novo motif identification from Ago2 peaks and prediction of miRNA::RNA heteroduplexes. miRBShunter was tested and experimentally validated on the in-house Ago2-dataset and on an Ago2-PAR-CLIP dataset in human stem cells. Overall, we provide guidelines to choose a suitable peak detection program and a new method for miRNA-target identification.

RNY-derived small RNAs as a signature of coronary artery disease.

BACKGROUND: Data from next generation sequencing technologies uncovered the existence of many classes of small RNAs. Recent studies reported that small RNAs are released by cells and can be detected in the blood. In this report, we aimed to discover the occurrence of novel circulating small RNAs in coronary artery disease (CAD). METHODS: We used high-throughput sequencing of small RNAs from human and mouse apoptotic primary macrophages, and analyzed the data by empirical Bayes moderated t-statistics to assess differential expression and the Benjamini and Hochberg method to control the false discovery rate. Results were then confirmed by Northern blot and RT-qPCR in foam cells and in two animal models for atherosclerosis, namely ApoE(-/-) and Ldlr(-/-) mouse lines. Quantitative RT-PCR to detect identified small RNAs, the RNY-derived small RNAs, was performed using sera of 263 patients with CAD compared to 514 matched healthy controls; the Student t-test was applied to statistically assess differences. Associations of small RNAs with clinical characteristics and biological markers were tested using Spearman's rank correlations, while multivariate logistic regressions were performed to test the statistical association of small RNA levels with CAD. RESULTS: Here, we report that, in macrophages stimulated with pro-apoptotic or pro-atherogenic stimuli, the Ro-associated non-coding RNAs, called RNYs or Y-RNAs, are processed into small RNAs (~24-34 nt) referred to as small-RNYs (s-RNYs), including s-RNY1-5p processed from RNY1. A significant upregulation of s-RNY expression was found in aortic arches and blood plasma from ApoE(-/-) and Ldlr(-/-) mice and in serum from CAD patients (P <0.001). Biostatistical analysis revealed a positive association of s-RNY1-5p with hs-CRP and ApoB levels; however, no statistical interaction was found between either of these two markers and s-RNY1-5p in relation to the CAD status. Levels of s-RNY1-5p were also independent from statin and fibrate therapies. CONCLUSION: Our results position the s-RNY1-5p as a relevant novel independent diagnostic biomarker for atherosclerosis-related diseases. Measurement of circulating s-RNY expression would be a valuable companion diagnostic to monitor foam cell apoptosis during atherosclerosis pathogenesis and to evaluate patient's responsiveness to future therapeutic strategies aiming to attenuate apoptosis in foam cells in advanced atherosclerotic lesions.

The RNA-binding protein KSRP promotes the biogenesis of a subset of microRNAs.

Consistent with the role of microRNAs (miRNAs) in down-regulating gene expression by reducing the translation and/or stability of target messenger RNAs, the levels of specific miRNAs are important for correct embryonic development and have been linked to several forms of cancer. However, the regulatory mechanisms by which primary miRNAs (pri-miRNAs) are processed first to precursor miRNAs (pre-miRNAs) and then to mature miRNAs by the multiprotein Drosha and Dicer complexes, respectively, remain largely unknown. The KH-type splicing regulatory protein (KSRP, also known as KHSRP) interacts with single-strand AU-rich-element-containing mRNAs and is a key mediator of mRNA decay. Here we show in mammalian cells that KSRP also serves as a component of both Drosha and Dicer complexes and regulates the biogenesis of a subset of miRNAs. KSRP binds with high affinity to the terminal loop of the target miRNA precursors and promotes their maturation. This mechanism is required for specific changes in target mRNA expression that affect specific biological programs, including proliferation, apoptosis and differentiation. These findings reveal an unexpected mechanism that links KSRP to the machinery regulating maturation of a cohort of miRNAs that, in addition to its role in promoting mRNA decay, independently serves to integrate specific regulatory programs of protein expression.

Let-7b/c enhance the stability of a tissue-specific mRNA during mammalian organogenesis as part of a feedback loop involving KSRP.

Gene silencing mediated by either microRNAs (miRNAs) or Adenylate/uridylate-rich elements Mediated mRNA Degradation (AMD) is a powerful way to post-transcriptionally modulate gene expression. We and others have reported that the RNA-binding protein KSRP favors the biogenesis of select miRNAs (including let-7 family) and activates AMD promoting the decay of inherently labile mRNAs. Different layers of interplay between miRNA- and AMD-mediated gene silencing have been proposed in cultured cells, but the relationship between the two pathways in living organisms is still elusive. We conditionally deleted Dicer in mouse pituitary from embryonic day (E) 9.5 through Cre-mediated recombination. In situ hybridization, immunohistochemistry, and quantitative reverse transcriptase-PCR revealed that Dicer is essential for pituitary morphogenesis and correct expression of hormones. Strikingly, αGSU (alpha glycoprotein subunit, common to three pituitary hormones) was absent in Dicer-deleted pituitaries. αGSU mRNA is unstable and its half-life increases during pituitary development. A transcriptome-wide analysis of microdissected E12.5 pituitaries revealed a significant increment of KSRP expression in conditional Dicer-deleted mice. We found that KSRP directly binds to αGSU mRNA, promoting its rapid decay; and, during pituitary development, αGSU expression displays an inverse temporal relationship to KSRP. Further, let-7b/c downregulated KSRP expression, promoting the degradation of its mRNA by directly binding to the 3'UTR. Therefore, we propose a model in which let-7b/c and KSRP operate within a negative feedback loop. Starting from E12.5, KSRP induces the maturation of let-7b/c that, in turn, post-transcriptionally downregulates the expression of KSRP itself. This event leads to stabilization of αGSU mRNA, which ultimately enhances the steady-state expression levels. We have identified a post-transcriptional regulatory network active during mouse pituitary development in which the expression of the hormone αGSU is increased by let7b/c through downregulation of KSRP. Our study unveils a functional crosstalk between miRNA- and AMD-dependent gene regulation during mammalian organogenesis events.

Impact of testicular cancer on sperm small non-coding RNA signature: a pilot study.

Testicular germ cell tumours (TGCTs) are the most common tumours in young adults of European ancestry. The high heritability and the constantly increased incidence, which has doubled over the last 20 years, strongly suggest that both genetic and environmental factors are likely to shape the TGCT susceptibility. While genome-wide association studies have identified loci associated with TGCT susceptibility, the role played by environmental molecular vectors in TGCT susceptibility remains unclear. Evidence shows that sperm non-coding RNAs provide a good vision of the environmental stresses experienced by men. Here, to determine whether TGCT impacts the abundance of specific non-coding RNAs in sperm, small RNA deep sequencing analysis of sperm of 25 men aged between 19 and 42 years, diagnosed with (n = 16) or without (n = 9) TGCT was performed. The primary analysis showed no statistical significance in the sncRNA population between the TGCT and non-TGCT groups. However, when sperm physiological parameters were considered to look for differentially expressed sncRNA, we evidenced 11 differentially expressed sncRNA between patients and control which allow a clear discrimination between control and TGCT samples after Hierarchical Clustering analysis. Together, these findings indicate that sperm small non-coding RNAs abundance may have the potential for diagnosing men with TGCT. However, specific care should be taken regarding sperm physiological parameters of the TGCT patients. Hence, larger studies are needed to confirm our findings and to determine whether such a signature associates with the risks to develop TGCT.

KSRP Promotes the Maturation of a Group of miRNA Precuresors.

microrNNA (mirNAs) are small noncoding RNAs that down-regulate gene expression by reducing stability and/or translation of target mRNAs. In animals, miRNAs arise from sequential processing of hairpin primary transcripts by two rNase III domain-containing enzymes, namely Drosha and Dicer, to generate a mature form of about 22 nucleotides. In this chapter we discuss our latest fndings indicating that KSRP is an integral component of both Drosha and Dicer complexes. KSRP binds to the terminal loop sequence of a subset of miRNA precursors promoting their maturation. our data indicate that the terminal loop is a pivotal structure where activators of miRNA processing as well as repressors of miRNA processing act in a coordinated way to convert cellular signals into changes in miRNA expression processing. This uncovers a new level of complexity of miRNA mechanisms for gene expression regulation.

Prediction of coronary heart disease incidence in a general male population by circulating non-coding small RNA sRNY1-5p in a nested case-control study.

During the development of atherosclerotic lesion, s-RNYs (small RNAs of about 24/34 nucleotides) are derived by the processing of long Ro-associated non-coding RNAs (RNYs) in macrophages. The levels of serum s-RNYs have been found significantly upregulated in patients with coronary heart disease (CHD) compared to age-matched CHD-free individuals. The present study aimed to examine the predictive value of serum s-RNYs for CHD events in the general male population. Within the frame of nested-case-control study, the GENES study, we measured the absolute expression of a RNY-derived small RNA, the s-RNY1-5p, in the serum of individuals (without CHD at baseline) who encountered a CHD event within 12 years of follow-up (n = 30) (Cases) and compared them to individuals who remained event-free (Controls) (n = 30). The expression of s-RNY1-5p in serum was significantly upregulated in Cases compared to Controls (p = 0.027). The proportion of CHD event-free was significantly higher among individuals with serum s-RNY1-5p below the median value (631 molecules/mL). In a multivariable model adjusted for age, smoking, hypertension, diabetes and dyslipidemia, the risk of CHD events increased more than fourfold in individuals with serum s-RNY1-5p above the median value (HR, 4.36; 95% CI 1.22-15.60). A positive association with CHD events was also observed when considering s-RNY1-5p as a continuous variable (p = 0.022). Based on our results, we conclude that serum s-RNY1-5p is an independent predictor of CHD events in a general male population and might be a relevant biomarker for early detection of cardiovascular diseases.

Paternal multigenerational exposure to an obesogenic diet drives epigenetic predisposition to metabolic diseases in mice.

Obesity is a growing societal scourge. Recent studies have uncovered that paternal excessive weight induced by an unbalanced diet affects the metabolic health of offspring. These reports mainly employed single-generation male exposure. However, the consequences of multigenerational unbalanced diet feeding on the metabolic health of progeny remain largely unknown. Here, we show that maintaining paternal Western diet feeding for five consecutive generations in mice induces an enhancement in fat mass and related metabolic diseases over generations. Strikingly, chow-diet-fed progenies from these multigenerational Western-diet-fed males develop a 'healthy' overweight phenotype characterized by normal glucose metabolism and without fatty liver that persists for four subsequent generations. Mechanistically, sperm RNA microinjection experiments into zygotes suggest that sperm RNAs are sufficient for establishment but not for long-term maintenance of epigenetic inheritance of metabolic pathologies. Progressive and permanent metabolic deregulation induced by successive paternal Western-diet-fed generations may contribute to the worldwide epidemic of metabolic diseases.

LPS induces KH-type splicing regulatory protein-dependent processing of microRNA-155 precursors in macrophages.

The importance of post-transcriptional mechanisms for the regulation of the homoeostasis of the immune system and the response to challenge by microorganisms is becoming increasingly appreciated. We investigated the contribution of microRNAs (miRNAs) to macrophage activation induced by lipopolysaccharide (LPS). We first observed that Dicer knockout in bone marrow-derived macrophages (BMDMs) increases the LPS-induced expression of some inflammation mediators. miRNA microarray analysis in BMDMs revealed that LPS significantly induces the expression of a single miRNA, miR-155, and this induction depends on enhanced miR-155 maturation from its precursors. The single-strand RNA-binding protein KH-type splicing regulatory protein (KSRP) binds to the terminal loop of miR-155 precursors and promotes their maturation. Both inhibition of miR-155 and KSRP knockdown enhance the LPS-induced expression of select inflammation mediators, and the effect of KSRP knockdown is reverted by mature miR-155. Our studies unveil the existence of an LPS-dependent post-transcriptional regulation of miR-155 biogenesis. Once induced, miR-155 finely tunes the expression of select inflammation mediators in response to LPS.

KSRP promotes the maturation of a group of miRNA precursors.

microRNAs (miRNAs) are small noncodingRNAs that down-regulate gene expression by reducing stability and/or translation of target mRNAs. In animals, miRNAs arise from sequential processing of hairpin primary transcripts by two RNAse III domain-containing enzymes, namely Drosha and Dicer, to generate a mature form of about 22 nucleotides. In this chapter we discuss our latest findings indicating that KSRP is an integral component of both Drosha and Dicer complexes. KSRP binds to the terminal loop sequence of a subset of miRNA precursors promoting their maturation. Our data indicate that the terminal loop is a pivotal structure where activators of miRNA processing as well as repressors of miRNA processing act in a coordinated way to convert cellular signals into changes in miRNA expression processing. This uncovers a new level of complexity of miRNA mechanisms for gene expression regulation.

The RNA-binding protein KSRP promotes decay of beta-catenin mRNA and is inactivated by PI3K-AKT signaling.

Beta-catenin plays an essential role in several biological events including cell fate determination, cell proliferation, and transformation. Here we report that beta-catenin is encoded by a labile transcript whose half-life is prolonged by Wnt and phosphatidylinositol 3-kinase-AKT signaling. AKT phosphorylates the mRNA decay-promoting factor KSRP at a unique serine residue, induces its association with the multifunctional protein 14-3-3, and prevents KSRP interaction with the exoribonucleolytic complex exosome. This impairs KSRP's ability to promote rapid mRNA decay. Our results uncover an unanticipated level of control of beta-catenin expression pointing to KSRP as a required factor to ensure rapid degradation of beta-catenin in unstimulated cells. We propose KSRP phosphorylation as a link between phosphatidylinositol 3-kinase-AKT signaling and beta-catenin accumulation.

How to control miRNA maturation?

In this point of view we discuss the role of co-activators and co-repressors of miRNA precursors maturation, the possibility that their functions are post translationally regulated by different signaling pathways, and their potential role in the miRNA-dependent control of cell proliferation and differentiation.

Molecular characterization and comparative localization of the mRNAs encoding two glutamic acid decarboxylases (GAD65 and GAD67) in the brain of the African lungfish, Protopterus annectens.

The existence of two distinct genes encoding two isoforms of glutamic acid decarboxylase (GAD65 and GAD67) has been demonstrated in most vertebrate classes, yet little is known about their differential distributions and functions in the central nervous system in nonmammalian vertebrates. In the present study, we have partially sequenced the cDNAs encoding GAD65 and GAD67 in the lungfish Protopterus annectens and determined their relative distributions in the adult brain by in situ hybridization histochemistry. The expression patterns of the GAD65 and GAD67 mRNAs were globally similar; the highest expression levels being observed in the granular layer of the olfactory bulb, the pallium, the subpallium, the anterior preoptic area, the thalamus, the hindbrain central gray, and the rhombencephalic visceral areas. However, striking differential expression was noticed in several structures. Very high to high concentrations of GAD67 mRNA were seen in the dorsal and ventral aspects of the anterior olfactory nucleus, which is in marked contrast to the very low expression of GAD65 in this region. Similarly, high levels of GAD67 mRNA were observed in the intermediate and ventral parts of the medial pallium that were virtually devoid of GAD65 mRNA. In contrast, GAD65 mRNA was found in the periaqueductal gray that did not express GAD67 mRNA. The differential expression of GAD65 and GAD67 mRNAs in these regions of the lungfish CNS indicates that the two GAD isoforms can be differentially regulated and that they may have distinct physiological roles.

Identification of a set of KSRP target transcripts upregulated by PI3K-AKT signaling.

BACKGROUND: KSRP is a AU-rich element (ARE) binding protein that causes decay of select sets of transcripts in different cell types. We have recently described that phosphatidylinositol 3-kinase/AKT (PI3K-AKT) activation induces stabilization and accumulation of the labile beta-catenin mRNA through an impairment of KSRP function. RESULTS: Aim of this study was to identify additional KSRP targets whose stability and steady-state levels are enhanced by PI3K-AKT activation. First, through microarray analyses of the AU-rich transcriptome in pituitary alphaT3-1 cells, we identified 34 ARE-containing transcripts upregulated in cells expressing a constitutively active form of AKT1. In parallel, by an affinity chromatography-based technique followed by microarray analyses, 12 mRNAs target of KSRP, additional to beta-catenin, were identified. Among them, seven mRNAs were upregulated in cells expressing activated AKT1. Both steady-state levels and stability of these new KSRP targets were consistently increased by either KSRP knock-down or PI3K-AKT activation. CONCLUSION: Our study identified a set of transcripts that are targets of KSRP and whose expression is increased by PI3K-AKT activation. These mRNAs encode RNA binding proteins, signaling molecules and a replication-independent histone. The increased expression of these gene products upon PI3K-AKT activation could play a role in the cellular events initiated by this signaling pathway.