Histone H3K27me3 demethylases regulate human Th17 cell development and effector functions by impacting on metabolism.

T helper (Th) cells are CD4+ effector T cells that play a critical role in immunity by shaping the inflammatory cytokine environment in a variety of physiological and pathological situations. Using a combined chemico-genetic approach, we identify histone H3K27 demethylases KDM6A and KDM6B as central regulators of human Th subsets. The prototypic KDM6 inhibitor GSK-J4 increases genome-wide levels of the repressive H3K27me3 chromatin mark and leads to suppression of the key transcription factor RORγt during Th17 differentiation. In mature Th17 cells, GSK-J4 induces an altered transcriptional program with a profound metabolic reprogramming and concomitant suppression of IL-17 cytokine levels and reduced proliferation. Single-cell analysis reveals a specific shift from highly inflammatory cell subsets toward a resting state upon demethylase inhibition. The root cause of the observed antiinflammatory phenotype in stimulated Th17 cells is reduced expression of key metabolic transcription factors, such as PPRC1. Overall, this leads to reduced mitochondrial biogenesis, resulting in a metabolic switch with concomitant antiinflammatory effects. These data are consistent with an effect of GSK-J4 on Th17 T cell differentiation pathways directly related to proliferation and include regulation of effector cytokine profiles. This suggests that inhibiting KDM6 demethylases may be an effective, even in the short term, therapeutic target for autoimmune diseases, including ankylosing spondylitis.

Identifying and avoiding off-target effects of RNase H-dependent antisense oligonucleotides in mice.

Antisense oligonucleotides that are dependent on RNase H for cleavage and subsequent degradation of complementary RNA are being developed as therapeutics. Besides the intended RNA target, such oligonucleotides may also cause degradation of unintended RNA off-targets by binding to partially complementary target sites. Here, we characterized the global effects on the mouse liver transcriptome of four oligonucleotides designed as gapmers, two targeting Apob and two targeting Pcsk9, all in different regions on their respective intended targets. This study design allowed separation of intended- and off-target effects on the transcriptome for each gapmer. Next, we used sequence analysis to identify possible partially complementary binding sites among the potential off-targets, and validated these by measurements of melting temperature and RNase H-cleavage rates. Generally, our observations were as expected in that fewer mismatches or bulges in the gapmer/transcript duplexes resulted in a higher chance of those duplexes being effective substrates for RNase H. Follow-up experiments in mice and cells show, that off-target effects can be mitigated by ensuring that gapmers have minimal sequence complementarity to any RNA besides the intended target, and that they do not have exaggerated binding affinity to the intended target.

Inhibition of histone H3K27 demethylases selectively modulates inflammatory phenotypes of natural killer cells.

Natural killer (NK) cells are innate lymphocytes, important in immune surveillance and elimination of stressed, transformed, or virus-infected cells. They critically shape the inflammatory cytokine environment to orchestrate interactions of cells of the innate and adaptive immune systems. Some studies have reported that NK cell activation and cytokine secretion are controlled epigenetically but have yielded only limited insight into the mechanisms. Using chemical screening with small-molecule inhibitors of chromatin methylation and acetylation, further validated by knockdown approaches, we here identified Jumonji-type histone H3K27 demethylases as key regulators of cytokine production in human NK cell subsets. The prototypic JMJD3/UTX (Jumonji domain-containing protein 3) H3K27 demethylase inhibitor GSK-J4 increased global levels of the repressive H3K27me3 mark around transcription start sites of effector cytokine genes. Moreover, GSK-J4 reduced IFN-γ, TNFα, granulocyte-macrophage colony-stimulating factor (GM-CSF), and interleukin-10 levels in cytokine-stimulated NK cells while sparing their cytotoxic killing activity against cancer cells. The anti-inflammatory effect of GSK-J4 in NK cell subsets, isolated from peripheral blood or tissue from individuals with rheumatoid arthritis (RA), coupled with an inhibitory effect on formation of bone-resorbing osteoclasts, suggested that histone demethylase inhibition has broad utility for modulating immune and inflammatory responses. Overall, our results indicate that H3K27me3 is a dynamic and important epigenetic modification during NK cell activation and that JMJD3/UTX-driven H3K27 demethylation is critical for NK cell function.

RNase H sequence preferences influence antisense oligonucleotide efficiency.

RNase H cleaves RNA in RNA-DNA duplexes. It is present in all domains of life as well as in multiple viruses and is essential for mammalian development and for human immunodeficiency virus replication. Here, we developed a sequencing-based method to measure the cleavage of thousands of different RNA-DNA duplexes and thereby comprehensively characterized the sequence preferences of HIV-1, human and Escherichia coli RNase H enzymes. We find that the catalytic domains of E. coli and human RNase H have nearly identical sequence preferences, which correlate with the efficiency of RNase H-recruiting antisense oligonucleotides. The sequences preferred by HIV-1 RNase H are distributed in the HIV genome in a way suggesting selection for efficient RNA cleavage during replication. Our findings can be used to improve the design of RNase H-recruiting antisense oligonucleotides and show that sequence preferences of HIV-1 RNase H may have shaped evolution of the viral genome and contributed to the use of tRNA-Lys3 as primer during viral replication.

Managing the sequence-specificity of antisense oligonucleotides in drug discovery.

All drugs perturb the expression of many genes in the cells that are exposed to them. These gene expression changes can be divided into effects resulting from engaging the intended target and effects resulting from engaging unintended targets. For antisense oligonucleotides, developments in bioinformatics algorithms, and the quality of sequence databases, allow oligonucleotide sequences to be analyzed computationally, in terms of the predictability of their interactions with intended and unintended RNA targets. Applying these tools enables selection of sequence-specific oligonucleotides where no- or only few unintended RNA targets are expected. To evaluate oligonucleotide sequence-specificity experimentally, we recommend a transcriptomics protocol where two or more oligonucleotides targeting the same RNA molecule, but with entirely different sequences, are evaluated together. This helps to clarify which changes in cellular RNA levels result from downstream processes of engaging the intended target, and which are likely to be related to engaging unintended targets. As required for all classes of drugs, the toxic potential of oligonucleotides must be evaluated in cell- and animal models before clinical testing. Since potential adverse effects related to unintended targeting are sequence-dependent and therefore species-specific, in vitro toxicology assays in human cells are especially relevant in oligonucleotide drug discovery.

Dissecting the target specificity of RNase H recruiting oligonucleotides using massively parallel reporter analysis of short RNA motifs.

Processing and post-transcriptional regulation of RNA often depend on binding of regulatory molecules to short motifs in RNA. The effects of such interactions are difficult to study, because most regulatory molecules recognize partially degenerate RNA motifs, embedded in a sequence context specific for each RNA. Here, we describe Library Sequencing (LibSeq), an accurate massively parallel reporter method for completely characterizing the regulatory potential of thousands of short RNA sequences in a specific context. By sequencing cDNA derived from a plasmid library expressing identical reporter genes except for a degenerate 7mer subsequence in the 3'UTR, the regulatory effects of each 7mer can be determined. We show that LibSeq identifies regulatory motifs used by RNA-binding proteins and microRNAs. We furthermore apply the method to cells transfected with RNase H recruiting oligonucleotides to obtain quantitative information for >15000 potential target sequences in parallel. These comprehensive datasets provide insights into the specificity requirements of RNase H and allow a specificity measure to be calculated for each tested oligonucleotide. Moreover, we show that inclusion of chemical modifications in the central part of an RNase H recruiting oligonucleotide can increase its sequence-specificity.

MicroRNA profiling identifies microRNA-155 as an adverse mediator of cardiac injury and dysfunction during acute viral myocarditis.

RATIONALE: Viral myocarditis results from an adverse immune response to cardiotropic viruses, which causes irreversible myocyte destruction and heart failure in previously healthy people. The involvement of microRNAs and their usefulness as therapeutic targets in this process are unknown. OBJECTIVE: To identify microRNAs involved in viral myocarditis pathogenesis and susceptibility. METHODS AND RESULTS: Cardiac microRNAs were profiled in both human myocarditis and in Coxsackievirus B3-injected mice, comparing myocarditis-susceptible with nonsusceptible mouse strains longitudinally. MicroRNA responses diverged depending on the susceptibility to myocarditis after viral infection in mice. MicroRNA-155, -146b, and -21 were consistently and strongly upregulated during acute myocarditis in both humans and susceptible mice. We found that microRNA-155 expression during myocarditis was localized primarily in infiltrating macrophages and T lymphocytes. Inhibition of microRNA-155 by a systemically delivered LNA-anti-miR attenuated cardiac infiltration by monocyte-macrophages, decreased T lymphocyte activation, and reduced myocardial damage during acute myocarditis in mice. These changes were accompanied by the derepression of the direct microRNA-155 target PU.1 in cardiac inflammatory cells. Beyond the acute phase, microRNA-155 inhibition reduced mortality and improved cardiac function during 7 weeks of follow-up. CONCLUSIONS: Our data show that cardiac microRNA dysregulation is a characteristic of both human and mouse viral myocarditis. The inflammatory microRNA-155 is upregulated during acute myocarditis, contributes to the adverse inflammatory response to viral infection of the heart, and is a potential therapeutic target for viral myocarditis.

LNA-mediated microRNA silencing in non-human primates.

microRNAs (miRNAs) are small regulatory RNAs that are important in development and disease and therefore represent a potential new class of targets for therapeutic intervention. Despite recent progress in silencing of miRNAs in rodents, the development of effective and safe approaches for sequence-specific antagonism of miRNAs in vivo remains a significant scientific and therapeutic challenge. Moreover, there are no reports of miRNA antagonism in primates. Here we show that the simple systemic delivery of a unconjugated, PBS-formulated locked-nucleic-acid-modified oligonucleotide (LNA-antimiR) effectively antagonizes the liver-expressed miR-122 in non-human primates. Acute administration by intravenous injections of 3 or 10 mg kg(-1) LNA-antimiR to African green monkeys resulted in uptake of the LNA-antimiR in the cytoplasm of primate hepatocytes and formation of stable heteroduplexes between the LNA-antimiR and miR-122. This was accompanied by depletion of mature miR-122 and dose-dependent lowering of plasma cholesterol. Efficient silencing of miR-122 was achieved in primates by three doses of 10 mg kg(-1) LNA-antimiR, leading to a long-lasting and reversible decrease in total plasma cholesterol without any evidence for LNA-associated toxicities or histopathological changes in the study animals. Our findings demonstrate the utility of systemically administered LNA-antimiRs in exploring miRNA function in rodents and primates, and support the potential of these compounds as a new class of therapeutics for disease-associated miRNAs.

mRNA expression profiling reveals conserved and non-conserved miR-140 targets.

microRNAs are non-coding RNAs that regulate gene expression. A significant proportion of microRNAs is perfectly conserved across the vertebrate clade, including miR-140, which is specifically expressed in cartilage. Although it has been computationally predicted that a large majority of microRNA targets are conserved, experimental evidence for this hypothesis remains scarce. In this work we use mRNA expression profiles obtained after manipulation of miR-140 activity levels in human and chicken primary chondrocytes to explore the extent of miR-140 target conservation. Our data suggest that miR-140 has a large number of targets conserved between human and chicken and we validate one of these, BMP2. However, we also found a significant number of non-conserved targets in the two species. In addition, we found that a commercially available scrambled siRNA, which is regularly used as a negative control, regulate the accumulation of many genes.

A large-scale genome-wide gene expression analysis in peripheral blood identifies very few differentially expressed genes related to antidepressant treatment and response in patients with major depressive disorder.

A better understanding of the biological factors underlying antidepressant treatment in patients with major depressive disorder (MDD) is needed. We perform gene expression analyses and explore sources of variability in peripheral blood related to antidepressant treatment and treatment response in patients suffering from recurrent MDD at baseline and after 8 weeks of treatment. The study includes 281 patients, which were randomized to 8 weeks of treatment with vortioxetine (N = 184) or placebo (N = 97). To our knowledge, this is the largest dataset including both gene expression in blood and placebo-controlled treatment response measured by a clinical scale in a randomized clinical trial. We identified three novel genes whose RNA expression levels at baseline and week 8 are significantly (FDR < 0.05) associated with treatment response after 8 weeks of treatment. Among these genes were SOCS3 (FDR = 0.0039) and PROK2 (FDR = 0.0028), which have previously both been linked to depression. Downregulation of these genes was associated with poorer treatment response. We did not identify any genes that were differentially expressed between placebo and vortioxetine groups at week 8 or between baseline and week 8 of treatment. Nor did we replicate any genes identified in previous peripheral blood gene expression studies examining treatment response. Analysis of genome-wide expression variability showed that type of treatment and treatment response explains very little of the variance, a median of <0.0001% and 0.05% in gene expression across all genes, respectively. Given the relatively large size of the study, the limited findings suggest that peripheral blood gene expression might not be the best approach to explore the biological factors underlying antidepressant treatment.

miRMaid: a unified programming interface for microRNA data resources.

BACKGROUND: MicroRNAs (miRNAs) are endogenous small RNAs that play a key role in post-transcriptional regulation of gene expression in animals and plants. The number of known miRNAs has increased rapidly over the years. The current release (version 14.0) of miRBase, the central online repository for miRNA annotation, comprises over 10.000 miRNA precursors from 115 different species. Furthermore, a large number of decentralized online resources are now available, each contributing with important miRNA annotation and information. RESULTS: We have developed a software framework, designated here as miRMaid, with the goal of integrating miRNA data resources in a uniform web service interface that can be accessed and queried by researchers and, most importantly, by computers. miRMaid is built around data from miRBase and is designed to follow the official miRBase data releases. It exposes miRBase data as inter-connected web services. Third-party miRNA data resources can be modularly integrated as miRMaid plugins or they can loosely couple with miRMaid as individual entities in the World Wide Web. miRMaid is available as a public web service but is also easily installed as a local application. The software framework is freely available under the LGPL open source license for academic and commercial use. CONCLUSION: miRMaid is an intuitive and modular software platform designed to unify miRBase and independent miRNA data resources. It enables miRNA researchers to computationally address complex questions involving the multitude of miRNA data resources. Furthermore, miRMaid constitutes a basic framework for further programming in which microRNA-interested bioinformaticians can readily develop their own tools and data sources.

Experimental identification of microRNA-140 targets by silencing and overexpressing miR-140.

MicroRNAs (miRNAs) are short noncoding RNA molecules regulating the expression of mRNAs. Target identification of miRNAs is computationally difficult due to the relatively low homology between miRNAs and their targets. We present here an experimental approach to target identification where the cartilage-specific miR-140 was overexpressed and silenced in cells it is normally expressed in separate experiments. Expression of mRNAs was profiled in both experiments and the intersection of mRNAs repressed by miR-140 overexpression and derepressed by silencing of miR-140 was identified. The intersection contained only 49 genes, although both treatments affected the accumulation of hundreds of mRNAs. These 49 genes showed a very strong enrichment for the miR-140 seed sequence implying that the approach is efficient and specific. Twenty-one of these 49 genes were predicted to be direct targets based on the presence of the seed sequence. Interestingly, none of these were predicted by the published target prediction methods we used. One of the potential target mRNAs, Cxcl12, was experimentally validated by Northern blot analysis and a luciferase reporter assay.

Antagonism of microRNA-122 in mice by systemically administered LNA-antimiR leads to up-regulation of a large set of predicted target mRNAs in the liver.

MicroRNA-122 (miR-122) is an abundant liver-specific miRNA, implicated in fatty acid and cholesterol metabolism as well as hepatitis C viral replication. Here, we report that a systemically administered 16-nt, unconjugated LNA (locked nucleic acid)-antimiR oligonucleotide complementary to the 5' end of miR-122 leads to specific, dose-dependent silencing of miR-122 and shows no hepatotoxicity in mice. Antagonism of miR-122 is due to formation of stable heteroduplexes between the LNA-antimiR and miR-122 as detected by northern analysis. Fluorescence in situ hybridization demonstrated uptake of the LNA-antimiR in mouse liver cells, which was accompanied by markedly reduced hybridization signals for mature miR-122 in treated mice. Functional antagonism of miR-122 was inferred from a low cholesterol phenotype and de-repression within 24 h of 199 liver mRNAs showing significant enrichment for miR-122 seed matches in their 3' UTRs. Expression profiling extended to 3 weeks after the last LNA-antimiR dose revealed that most of the changes in liver gene expression were normalized to saline control levels coinciding with normalized miR-122 and plasma cholesterol levels. Combined, these data suggest that miRNA antagonists comprised of LNA are valuable tools for identifying miRNA targets in vivo and for studying the biological role of miRNAs and miRNA-associated gene-regulatory networks in a physiological context.

The utility of LNA in microRNA-based cancer diagnostics and therapeutics.

MicroRNAs (miRNAs) are a novel class of small endogenous non-coding RNAs that regulate gene expression post-transcriptionally by binding to their cognate target mRNAs. Emerging evidence implies that miRNAs play important roles in cancer and thus, miRNAs have rapidly emerged as valuable markers for cancer diagnostics and promising targets for therapeutics. Locked nucleic acid (LNA) is a conformational RNA analoque that binds complementary RNA with unprecedented affinity and specificity. These properties make LNA well suited for miRNA detection and analysis for cancer diagnostics. Furthermore, recent studies on LNA-mediated silencing of miRNA function in vitro and in vivo support the potential of LNA in therapeutic intervention of cancer-associated miRNAs.

Targeting Repeated Regions Unique to a Gene Is an Effective Strategy for Discovering Potent and Efficacious Antisense Oligonucleotides.

Well-validated strategies for discovering potent and efficacious antisense oligonucleotides are central to realize the full therapeutic potential of RNA therapy. In this study, we focus on RNA targets where the same sequence of 16-20 nt is found in several regions across the RNA, and not in any other RNA. Targeting such unique repeated regions with oligonucleotides designed as gapmers and capable of recruiting RNase H has previously been proposed as a strategy for identifying potent gapmers. By sequence analysis of the human and monkey transcriptomes, we find that such unique repeated regions in RNA are often conserved between humans and monkeys, which allow pharmacodynamic effects to be evaluated in non-human primates before testing in humans. For eight potential RNA targets chosen in an unbiased fashion, we targeted their unique repeated regions with locked nucleic acid (LNA)-modified gapmers, and for six of them we identified gapmers that were significantly more potent and efficacious in vitro than non-repeat-targeting gapmer controls. We suggest a stochastic model for repeat-targeting gapmers that explains all effects observed so far and can help guide future work. Our results support the targeting of repeated regions as an effective strategy for discovering gapmer antisense oligonucleotides suitable for therapeutic development.

Targeting of microRNAs for therapeutics.

miRNAs (microRNAs) comprise a class of small endogenous non-coding RNAs that post-transcriptionally repress gene expression by base-pairing with their target mRNAs. Recent evidence has shown that miRNAs play important roles in a wide variety of human diseases, such as viral infections, cancer and cardiovascular diseases, and thus miRNAs have rapidly emerged as potential targets for therapeutics. LNAs (locked nucleic acids) comprise a class of bicyclic conformational analogues of RNA, which exhibit high binding affinity to complementary RNA molecules and high stability in blood and tissues in vivo. Recent reports on LNA-mediated miRNA silencing in rodents and primates support the potential of LNA-modified oligonucleotides in studying miRNA functions in vivo and in the future development of miRNA-based therapeutics.

Intragenomic matching reveals a huge potential for miRNA-mediated regulation in plants.

microRNAs (miRNAs) are important post-transcriptional regulators, but the extent of this regulation is uncertain, both with regard to the number of miRNA genes and their targets. Using an algorithm based on intragenomic matching of potential miRNAs and their targets coupled with support vector machine classification of miRNA precursors, we explore the potential for regulation by miRNAs in three plant genomes: Arabidopsis thaliana, Populus trichocarpa, and Oryza sativa. We find that the intragenomic matching in conjunction with a supervised learning approach contains enough information to allow reliable computational prediction of miRNA candidates without requiring conservation across species. Using this method, we identify approximately 1,200, approximately 2,500, and approximately 2,100 miRNA candidate genes capable of extensive base-pairing to potential target mRNAs in A. thaliana, P. trichocarpa, and O. sativa, respectively. This is more than five times the number of currently annotated miRNAs in the plants. Many of these candidates are derived from repeat regions, yet they seem to contain the features necessary for correct processing by the miRNA machinery. Conservation analysis indicates that only a few of the candidates are conserved between the species. We conclude that there is a large potential for miRNA-mediated regulatory interactions encoded in the genomes of the investigated plants. We hypothesize that some of these interactions may be realized under special environmental conditions, while others can readily be recruited when organisms diverge and adapt to new niches.

Principles and limitations of computational microRNA gene and target finding.

In 2001 there were four PubMed entries matching the word "microRNA" (miRNA). Interestingly, this number has now far exceeded 1300 and is still rapidly increasing. This more than anything demonstrates the extreme attention this field has had within a short period of time. With the large amounts of sequence data being generated, the need for analysis by computational approaches is obvious. Here, we review the general principles used in computational gene and target finding, and discuss the strengths and weaknesses of the methods. Several methods rely on detection of evolutionary conserved candidates, but recent methods have challenged this paradigm by simultaneously searching for the gene and the corresponding target(s). Whereas the early methods made predictions based on sets of hand-derived rules from precursor-miRNA structure or observed target-miRNA interactions, recent methods apply machine learning techniques. Even though these methods are already powerful, the amount of data they rely on is still limited. Since it is evident that data are continuously being generated, it must be anticipated that these methods will further improve their performance.

Computational evidence for hundreds of non-conserved plant microRNAs.

BACKGROUND: MicroRNAs (miRNA) are small (20-25 nt) non-coding RNA molecules that regulate gene expression through interaction with mRNA in plants and metazoans. A few hundred miRNAs are known or predicted, and most of those are evolutionarily conserved. In general plant miRNA are different from their animal counterpart: most plant miRNAs show near perfect complementarity to their targets. Exploiting this complementarity we have developed a method for identification plant miRNAs that does not rely on phylogenetic conservation. RESULTS: Using the presumed targets for the known miRNA as positive controls, we list and filter all segments of the genome of length approximately 20 that are complementary to a target mRNA-transcript. From the positive control we recover 41 (of 92 possible) of the already known miRNA-genes (representing 14 of 16 families) with only four false positives. Applying the procedure to find possible new miRNAs targeting any annotated mRNA, we predict of 592 new miRNA genes, many of which are not conserved in other plant genomes. A subset of our predicted miRNAs is additionally supported by having more than one target that are not homologues. CONCLUSION: These results indicate that it is possible to reliably predict miRNA-genes without using genome comparisons. Furthermore it suggests that the number of plant miRNAs have been underestimated and points to the existence of recently evolved miRNAs in Arabidopsis.