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.

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.

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.

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.

Diagnostic microRNA profiling in cutaneous T-cell lymphoma (CTCL).

Cutaneous T-cell lymphomas (CTCLs) are the most frequent primary skin lymphomas. Nevertheless, diagnosis of early disease has proven difficult because of a clinical and histologic resemblance to benign inflammatory skin diseases. To address whether microRNA (miRNA) profiling can discriminate CTCL from benign inflammation, we studied miRNA expression levels in 198 patients with CTCL, peripheral T-cell lymphoma (PTL), and benign skin diseases (psoriasis and dermatitis). Using microarrays, we show that the most induced (miR-326, miR-663b, and miR-711) and repressed (miR-203 and miR-205) miRNAs distinguish CTCL from benign skin diseases with > 90% accuracy in a training set of 90 samples and a test set of 58 blinded samples. These miRNAs also distinguish malignant and benign lesions in an independent set of 50 patients with PTL and skin inflammation and in experimental human xenograft mouse models of psoriasis and CTCL. Quantitative (q)RT-PCR analysis of 103 patients with CTCL and benign skin disorders validates differential expression of 4 of the 5 miRNAs and confirms previous reports on miR-155 in CTCL. A qRT-PCR-based classifier consisting of miR-155, miR-203, and miR-205 distinguishes CTCL from benign disorders with high specificity and sensitivity, and with a classification accuracy of 95%, indicating that miRNAs have a high diagnostic potential in CTCL.

Impact of climate change on vector- and rodent-borne infectious diseases.

Background: Endemic and imported vector- and rodent-borne infectious agents can be linked to high morbidity and mortality. Therefore, vector- and rodent-borne human diseases and the effects of climate change are important public health issues. Methods: For this review, the relevant literature was identified and evaluated according to the thematic aspects and supplemented with an analysis of surveillance data for Germany. Results: Factors such as increasing temperatures, changing precipitation patterns, and human behaviour may influence the epidemiology of vector- and rodent-borne infectious diseases in Germany. Conclusions: The effects of climatic changes on the spread of vector- and rodent-borne infectious diseases need to be further studied in detail and considered in the context of climate adaptation measures.

Comparison of microRNA expression using different preservation methods of matched psoriatic skin samples.

MicroRNAs are non-coding RNA molecules modulating gene expression post-transcriptionally. Formalin-fixed, paraffin-embedding (FFPE) is a standard preservation method often used in clinical practices, but induces RNA degradation. Extracting high-quality RNA from human skin can be challenging as skin contains high levels of RNases. As microRNAs are 19-23 nucleotides long and lack a poly-A tail, they may be less prone to RNA degradation than mRNAs. We investigated whether microRNAs in psoriatic (FFPE) samples reliably reflect microRNA expression in samples less prone to RNA degradation such as fresh-frozen (FS) and Tissue-Tek-embedding (OCT). We found a strong correlation of the microRNA expression levels between all preservation methods of matched psoriatic skin samples (r(s) ranging from 0.91 to 0.95 (P < 0.001)). These observations were further confirmed with qRT-PCR. Our results demonstrate that microRNA detection in human skin is robust irrespective of preservation method; thus, microRNAs offer an appropriate and flexible approach in clinical practices and for diagnostic purposes in skin disorders.

Allergic contact dermatitis induces upregulation of identical microRNAs in humans and mice.

BACKGROUND: MicroRNAs are short, endogenous RNA molecules that can bind to parts of target mRNAs, thus inhibiting their translation and causing accelerated turnover or degradation of transcripts, thereby regulating gene expression. Several microRNAs have been found to be upregulated in atopic dermatitis and psoriasis, indicating a role in inflammatory skin diseases. However, there have been no studies on the expression of microRNAs in allergic contact dermatitis. OBJECTIVES: To investigate expression of microRNAs in allergic contact dermatitis. Methods. Lesional and non-lesional skin biopsies were collected from subjects with allergic responses to diphenylcyclopropenone (DPCP). Additional samples for profiling were collected from an experimental mouse model by use of the strong allergen dinitrofluorobenzene. RNA was purified from all samples, and locked nucleic acid microarray analysis was performed, followed by validation with quantitative polymerase chain reaction (PCR). RESULTS: In humans sensitized with DPCP, we found significant upregulation of miR-21, miR-142-3p, miR-142-5p and miR-223 in challenged skin. The same microRNAs were significantly upregulated in the skin of mice in a mouse model of contact allergy. The upregulation of microRNA was confirmed by quantitative PCR. CONCLUSION: These are the first results indicating that microRNAs may be involved in the pathogenesis of allergic contact dermatitis, and they show that mouse models are valuable tools for further study of the involvement of microRNAs in allergic contact dermatitis.

Acute Neurotoxicity of Antisense Oligonucleotides After Intracerebroventricular Injection Into Mouse Brain Can Be Predicted from Sequence Features.

Antisense oligonucleotides are a relatively new therapeutic modality and safety evaluation is still a developing area of research. We have observed that some oligonucleotides can produce acute, nonhybridization dependent, neurobehavioral side effects after intracerebroventricular (ICV) dosing in mice. In this study, we use a combination of in vitro, in vivo, and bioinformatics approaches to identify a sequence design algorithm, which can reduce the number of acutely toxic molecules synthesized and tested in mice. We find a cellular assay measuring spontaneous calcium oscillations in neuronal cells can predict the behavioral side effects after ICV dosing, and may provide a mechanistic explanation for these observations. We identify sequence features that are overrepresented or underrepresented among oligonucleotides causing these reductions in calcium oscillations. A weighted linear combination of the five most informative sequence features predicts the outcome of ICV dosing with >80% accuracy. From this, we develop a bioinformatics tool that allows oligonucleotide designs with acceptable acute neurotoxic potential to be identified, thereby reducing the number of toxic molecules entering drug discovery pipelines. The informative sequence features we identified also suggest areas in which to focus future medicinal chemistry efforts.

Identification and characterization of a MAPT-targeting locked nucleic acid antisense oligonucleotide therapeutic for tauopathies.

Tau is a microtubule-associated protein (MAPT, tau) implicated in the pathogenesis of tauopathies, a spectrum of neurodegenerative disorders characterized by accumulation of hyperphosphorylated, aggregated tau. Because tau pathology can be distinct across diseases, a pragmatic therapeutic approach may be to intervene at the level of the tau transcript, as it makes no assumptions to mechanisms of tau toxicity. Here we performed a large library screen of locked-nucleic-acid (LNA)-modified antisense oligonucleotides (ASOs), where careful tiling of the MAPT locus resulted in the identification of hot spots for activity in the 3' UTR. Further modifications to the LNA design resulted in the generation of ASO-001933, which selectively and potently reduces tau in primary cultures from hTau mice, monkey, and human neurons. ASO-001933 was well tolerated and produced a robust, long-lasting reduction in tau protein in both mouse and cynomolgus monkey brain. In monkey, tau protein reduction was maintained in brain for 20 weeks post injection and corresponded with tau protein reduction in the cerebrospinal fluid (CSF). Our results demonstrate that LNA-ASOs exhibit excellent drug-like properties and sustained efficacy likely translating to infrequent, intrathecal dosing in patients. These data further support the development of LNA-ASOs against tau for the treatment of tauopathies.

Independent component and pathway-based analysis of miRNA-regulated gene expression in a model of type 1 diabetes.

BACKGROUND: Several approaches have been developed for miRNA target prediction, including methods that incorporate expression profiling. However the methods are still in need of improvements due to a high false discovery rate. So far, none of the methods have used independent component analysis (ICA). Here, we developed a novel target prediction method based on ICA that incorporates both seed matching and expression profiling of miRNA and mRNA expressions. The method was applied on a cellular model of type 1 diabetes. RESULTS: Microarray profiling identified eight miRNAs (miR-124/128/192/194/204/375/672/708) with differential expression. Applying ICA on the mRNA profiling data revealed five significant independent components (ICs) correlating to the experimental conditions. The five ICs also captured the miRNA expressions by explaining > 97% of their variance. By using ICA, seven of the eight miRNAs showed significant enrichment of sequence predicted targets, compared to only four miRNAs when using simple negative correlation. The ICs were enriched for miRNA targets that function in diabetes-relevant pathways e.g. type 1 and type 2 diabetes and maturity onset diabetes of the young (MODY). CONCLUSIONS: In this study, ICA was applied as an attempt to separate the various factors that influence the mRNA expression in order to identify miRNA targets. The results suggest that ICA is better at identifying miRNA targets than negative correlation. Additionally, combining ICA and pathway analysis constitutes a means for prioritizing between the predicted miRNA targets. Applying the method on a model of type 1 diabetes resulted in identification of eight miRNAs that appear to affect pathways of relevance to disease mechanisms in diabetes.

Integrative analysis correlates donor transcripts to recipient autoantibodies in primary graft dysfunction after lung transplantation.

Up to one in four lung-transplanted patients develop pulmonary infiltrates and impaired oxygenation within the first days after lung transplantation. Known as primary graft dysfunction (PGD), this condition increases mortality significantly. Complex interactions between donor lung and recipient immune system are the suspected cause. We took an integrative, systems-level approach by first exploring whether the recipient's immune response to PGD includes the development of long-lasting autoreactivity. We next explored whether proteins displaying such differential autoreactivity also display differential gene expression in donor lungs that later develop PGD compared with those that did not. We evaluated 39 patients from whom autoantibody profiles were already available for PGD based on chest radiographs and oxygenation data. An additional nine patients were evaluated for PGD based on their medical records and set aside for validation. From two recent donor lung gene expression studies, we reanalysed and paired gene profiles with autoantibody profiles. Primary graft dysfunction can be distinguished by a profile of differentially reactive autoantibodies binding to 17 proteins. Functional analysis showed that 12 of these proteins are part of a protein-protein interaction network (P=3 x 10⁻6) involved in proliferative processes. A nearest centroid classifier assigned correct PGD grades to eight out of the nine patients in the validation cohort (P=0·048). We observed significant positive correlation (r=0·63, P=0·011) between differences in IgM reactivity and differences in gene expression levels. This connection between donor lung gene expression and long-lasting recipient IgM autoantibodies towards a specific set of proteins suggests a mechanism for the development of autoimmunity in PGD.

Factors affecting patterns of tick parasitism on forest rodents in tick-borne encephalitis risk areas, Germany.

Identifying factors affecting individual vector burdens is essential for understanding infectious disease systems. Drawing upon data of a rodent monitoring programme conducted in nine different forest patches in southern Hesse, Germany, we developed models which predict tick (Ixodes spp. and Dermacentor spp.) burdens on two rodent species Apodemus flavicollis and Myodes glareolus. Models for the two rodent species were broadly similar but differed in some aspects. Patterns of Ixodes spp. burdens were influenced by extrinsic factors such as season, unexplained spatial variation (both species), relative humidity and vegetation cover (A. flavicollis). We found support for the 'body mass' (tick burdens increase with body mass/age) and for the 'dilution' hypothesis (tick burdens decline with increasing rodent densities) and little support for the 'sex-bias' hypothesis (both species). Surprisingly, roe deer densities were not correlated with larvae counts on rodents. Factors influencing the mean burden did not significantly explain the observed dispersion of tick counts. Co-feeding aggregations, which are essential for tick-borne disease transmission, were mainly found in A. flavicollis of high body mass trapped in areas with fast increase in spring temperatures. Locally, Dermacentor spp. appears to be an important parasite on A. flavicollis and M. glareolus. Dermacentor spp. was rather confined to areas with higher average temperatures during the vegetation period. Nymphs of Dermacentor spp. mainly fed on M. glareolus and were seldom found on A. flavicollis. Whereas Ixodes spp. is the dominant tick genus in woodlands of our study area, the distribution and epidemiological role of Dermacentor spp. should be monitored closely.

Characterization of Escherichia coli RNase H Discrimination of DNA Phosphorothioate Stereoisomers.

Phosphorothioate (PS) modification of antisense oligonucleotides (ASOs) is a critical factor enabling their therapeutic use. Standard chemical synthesis incorporates this group in a stereorandom manner; however, significant effort was made over the years to establish and characterize the impact of chiral control. In this work, we present our in-depth characterization of interactions between Escherichia coli RNase H and RNA-DNA heteroduplexes carrying chirally defined PS groups. First, using a massive parallel assay, we showed that at least a single Rp-PS group is necessary for efficient RNase H-mediated cleavage. We followed by demonstrating that this group needs to be aligned to the phosphate-binding pocket of RNase H, and that chiral status of other PS groups in close proximity to RNase H does not affect cleavage efficiency. We have shown that RNase H's PS chiral preference can be utilized to guide cleavage to a specific chemical bond. Finally, we present a strategy for ASO optimization by mapping preferred RNase H cleavage sites of a non-thioated compound, followed by introduction of Rp-PS in a strategic position. This results in a cleaner cleavage profile and higher knockdown activity compared with a compound carrying an Sp-PS at the same location.

Anatomical work-up of an individual with multiple muscular variants on both forearms.

Muscular variants of the forearm are common and frequently cause neurovascular compression syndromes, especially when interfering with the compact topography of the carpal tunnel or the Canalis ulnaris. Here, we report on a male body donor with multiple muscular normal variations on both forearms. The two main findings are (1) an accessory variant muscle (AVM) on the right forearm originating from the M. brachioradialis, the distal radius, and the M. flexor pollicis longus. It spanned the wrist beneath the Fascia antebrachia and inserted at the proximal phalanx of the digitus minimus. (2) Moreover, we found a three-headed palmaris longus variant on the left arm with proximal origin tendon and a distal, trifurcated muscle belly, with separated insertions at the palmar aponeurosis, the flexor retinaculum, and, in analogy to the accessory muscle on the contralateral arm, at the base of the proximal phalanx of the digitus minimus. We found a considerable thickening of the left-hand median nerve right before entering the carpal tunnel indicative of a possible chronic compression syndrome adding clinical relevance to this anatomical case. We also discuss the notion that both, the AVM and the contralateral three-headed palmaris variant are developmental descendants of the M. palmaris longus. Additionally, we found a previously not recorded variant of the M. palmaris brevis on the left hand.

Chronic rejection of a lung transplant is characterized by a profile of specific autoantibodies.

SUMMARY: Obliterative bronchiolitis (OB) continues to be the major limitation to long-term survival after lung transplantation. The specific aetiology and pathogenesis of OB are not well understood. To explore the role of autoreactivity in OB, we spotted 751 different self molecules onto glass slides, and used these antigen microarrays to profile 48 human serum samples for immunoglobulin G (IgG) and IgM autoantibodies; 27 patients showed no or mild bronchiolitis obliterans syndrome (BOS; a clinical correlate of OB) and 15 patients showed medium to severe BOS. We now report that these BOS grades could be differentiated by a profile of autoantibodies binding to 28 proteins or their peptides. The informative autoantibody profile included down-regulation as well as up-regulation of both IgM and IgG specific reactivities. This profile was evaluated for robustness using a panel of six independent test patients. Analysis of the functions of the 28 informative self antigens showed that eight of them are connected in an interaction network involved in apoptosis and protein metabolism. Thus, a profile of autoantibodies may reflect pathological processes in the lung allograft, suggesting a role for autoimmunity in chronic rejection leading to OB.

Determinants of tick-borne encephalitis in counties of southern Germany, 2001-2008.

BACKGROUND: Tick-borne encephalitis (TBE) virus can cause severe symptoms in humans. The incidence of this vector-borne pathogen in humans is characterised by spatial and temporal heterogeneity. To explain the variation in reported human TBE cases per county in southern Germany, we designed a time-lagged, spatially-explicit model that incorporates ecological, environmental, and climatic factors. RESULTS: We fitted a logistic regression model to the annual counts of reported human TBE cases in each of 140 counties over an eight year period. The model controlled for spatial autocorrelation and unexplained temporal variation. The occurrence of human TBE was found to be positively correlated with the proportions of broad-leafed, mixed and coniferous forest cover. An index of forest fragmentation was negatively correlated with TBE incidence, suggesting that infection risk is higher in fragmented landscapes. The results contradict previous evidence regarding the relevance of a specific spring-time temperature regime for TBE epidemiology. Hunting bag data of roe deer (Capreolus capreolus) in the previous year was positively correlated with human TBE incidence, and hunting bag density of red fox (Vulpes vulpes) and red deer (Cervus elaphus) in the previous year were negatively correlated with human TBE incidence. CONCLUSIONS: Our approach suggests that a combination of landscape and climatic variables as well as host-species dynamics influence TBE infection risk in humans. The model was unable to explain some of the temporal variation, specifically the high counts in 2005 and 2006. Factors such as the exposure of humans to infected ticks and forest rodent population dynamics, for which we have no data, are likely to be explanatory factors. Such information is required to identify the determinants of TBE more reliably. Having records of TBE infection sites at a finer scale would also be necessary.

Tick burden on European roe deer (Capreolus capreolus).

In our study we assessed the tick burden on roe deer (Capreolus capreolus L.) in relation to age, physical condition, sex, deer density and season. The main objective was to find predictive parameters for tick burden. In September 2007, May, July, and September 2008, and in May and July 2009 we collected ticks on 142 culled roe deer from nine forest departments in Southern Hesse, Germany. To correlate tick burden and deer density we estimated deer density using line transect sampling that accounts for different detectability in March 2008 and 2009, respectively. We collected more than 8,600 ticks from roe deer heads and necks, 92.6% of which were Ixodes spp., 7.4% Dermacentor spp. Among Ixodes, 3.3% were larvae, 50.5% nymphs, 34.8% females and 11.4% males, with significant seasonal deviation. Total tick infestation was high, with considerable individual variation (from 0 to 270 ticks/deer). Adult tick burden was positively correlated with roe deer body indices (body mass, age, hind foot length). Significantly more nymphs were found on deer from forest departments with high roe deer density indices, indicating a positive correlation with deer abundance. Overall, tick burden was highly variable. Seasonality and large scale spatial characteristics appeared to be the most important factors affecting tick burden on roe deer.

Chemical Diversity of Locked Nucleic Acid-Modified Antisense Oligonucleotides Allows Optimization of Pharmaceutical Properties.

The identification of molecules that can modulate RNA or protein function and the subsequent chemical and structural optimization to refine such molecules into drugs is a key activity in drug discovery. Here, we explored the extent to which chemical and structural differences in antisense oligonucleotides, designed as gapmers and capable of recruiting RNase H for target RNA cleavage, can affect their functional properties. To facilitate structure-activity learning, we analyzed two sets of iso-sequential locked nucleic acid (LNA)-modified gapmers, where we systematically varied the number and positions of LNA modifications in the flanks. In total, we evaluated 768 different and architecturally diverse gapmers in HeLa cells for target knockdown activity and cytotoxic potential and found widespread differences in both of these properties. Binding affinity between gapmer and RNA target, as well as the presence of certain short sequence motifs in the gap region, can explain these differences, and we propose statistical and machine-learning models that can be used to predict region-specific, optimal LNA-modification architectures. Once accessible regions in the target of interest have been identified, our results show how to refine and optimize LNA gapmers with improved pharmacological profiles targeting such regions.