Chemical modification patterns for microRNA therapeutic mimics: a structure-activity relationship (SAR) case-study on miR-200c.

microRNA (miRNA) mimics are an emerging class of oligonucleotide therapeutics, with a few compounds already in clinical stages. Synthetic miRNAs are able to restore downregulated levels of intrinsic miRNAs, allowing for parallel regulation of multiple genes involved in a particular disease. In this work, we examined the influence of chemical modifications patterns in miR-200c mimics, assessing the regulation of a selection of target messenger RNAs (mRNA) and, subsequently, of the whole transcriptome in A549 cells. We have probed 37 mimics and provided an initial set of instructions for designing miRNA mimics with potency and selectivity similar to an unmodified miRNA duplex. Additionally, we have examined the stability of selected mimics in serum. Finally, the selected two modification patterns were translated to two other miRNAs, miR-34a and miR-155. To differing degrees, these designs acted on target mRNAs in a similar manner to the unmodified mimic. Here, for the first time, we describe a structured overview of 'miRNA mimics modification templates' that are chemically stabilised and optimised for use in an in vitro set up and highlight the need of further sequence specific optimization when mimics are to be used beyond in vitro tool experiments.

Tissue pharmacokinetics of antisense oligonucleotides.

Pharmacokinetics (PK) of antisense oligonucleotides (ASOs) is characterized by rapid distribution from plasma to tissue and slow terminal plasma elimination driven by re-distribution from tissue. Quantitative understanding of tissue PK and RNA knockdown for various ASO chemistries, conjugations, and administration routes is critical for successful drug discovery. Here, we report concentration-time and RNA knockdown profiles for a gapmer ASO with locked nucleic acid ribose chemistry in mouse liver, kidney, heart, and lung after subcutaneous and intratracheal administration. Additionally, the same ASO with liver targeting conjugation (galactosamine-N-acetyl) is evaluated for subcutaneous administration. Data indicate that exposure and knockdown differ between tissues and strongly depend on administration route and conjugation. In a second study, we show that tissue PK is similar between the three different ribose chemistries locked nucleic acid, constrained ethyl and 2'-O-methoxyethyl, both after subcutaneous and intratracheal administration. Further, we show that the half-life in mouse liver may vary with ASO sequence. Finally, we report less than dose-proportional increase in liver concentration in the dose range of 3-30 µmol/kg. Overall, our studies contribute pivotal data to support design and interpretation of ASO in vivo studies, thereby increasing the probability of delivering novel ASO therapies to patients.

A universal model of RNA.DNA:DNA triplex formation accurately predicts genome-wide RNA-DNA interactions.

RNA.DNA:DNA triple helix (triplex) formation is a form of RNA-DNA interaction which regulates gene expression but is difficult to study experimentally in vivo. This makes accurate computational prediction of such interactions highly important in the field of RNA research. Current predictive methods use canonical Hoogsteen base pairing rules, which whilst biophysically valid, may not reflect the plastic nature of cell biology. Here, we present the first optimization approach to learn a probabilistic model describing RNA-DNA interactions directly from motifs derived from triplex sequencing data. We find that there are several stable interaction codes, including Hoogsteen base pairing and novel RNA-DNA base pairings, which agree with in vitro measurements. We implemented these findings in TriplexAligner, a program that uses the determined interaction codes to predict triplex binding. TriplexAligner predicts RNA-DNA interactions identified in all-to-all sequencing data more accurately than all previously published tools in human and mouse and also predicts previously studied triplex interactions with known regulatory functions. We further validated a novel triplex interaction using biophysical experiments. Our work is an important step towards better understanding of triplex formation and allows genome-wide analyses of RNA-DNA interactions.

Deep sequencing of short capped RNAs reveals novel families of noncoding RNAs.

In eukaryotes, capped RNAs include long transcripts such as messenger RNAs and long noncoding RNAs, as well as shorter transcripts such as spliceosomal RNAs, small nucleolar RNAs, and enhancer RNAs. Long capped transcripts can be profiled using cap analysis gene expression (CAGE) sequencing and other methods. Here, we describe a sequencing library preparation protocol for short capped RNAs, apply it to a differentiation time course of the human cell line THP-1, and systematically compare the landscape of short capped RNAs to that of long capped RNAs. Transcription initiation peaks associated with genes in the sense direction have a strong preference to produce either long or short capped RNAs, with one out of six peaks detected in the short capped RNA libraries only. Gene-associated short capped RNAs have highly specific 3' ends, typically overlapping splice sites. Enhancers also preferentially generate either short or long capped RNAs, with 10% of enhancers observed in the short capped RNA libraries only. Enhancers producing either short or long capped RNAs show enrichment for GWAS-associated disease SNPs. We conclude that deep sequencing of short capped RNAs reveals new families of noncoding RNAs and elucidates the diversity of transcripts generated at known and novel promoters and enhancers.

Recombination of repeat elements generates somatic complexity in human genomes.

Non-allelic recombination between homologous repetitive elements contributes to evolution and human genetic disorders. Here, we combine short- and long-DNA read sequencing of repeat elements with a new bioinformatics pipeline to show that somatic recombination of Alu and L1 elements is widespread in the human genome. Our analysis uncovers tissue-specific non-allelic homologous recombination hallmarks; moreover, we find that centromeres and cancer-associated genes are enriched for retroelements that may act as recombination hotspots. We compare recombination profiles in human-induced pluripotent stem cells and differentiated neurons and find that the neuron-specific recombination of repeat elements accompanies chromatin changes during cell-fate determination. Finally, we report that somatic recombination profiles are altered in Parkinson's and Alzheimer's disease, suggesting a link between retroelement recombination and genomic instability in neurodegeneration. This work highlights a significant contribution of the somatic recombination of repeat elements to genomic diversity in health and disease.

Analysis of Enhancer-Promoter Interactions using CAGE and RADICL-Seq Technologies.

Regulation of gene expression is a key feature for higher eukaryotes and how chromatin topology relates to gene activation is an intense area of research. Enhancer-promoter interactions are believed to mediate activation of target genes. Bidirectional transcription represents one hallmark of active enhancers that can be measured using transcriptome technologies such as Cap analysis of gene expression (CAGE). Recently, we have developed RNA and DNA interacting complexes ligated and sequenced (RADICL-Seq) a novel methodology to map genome-wide RNA-chromatin interactions in intact nuclei. Here, we describe how CAGE and RADICL-Seq data can be used to characterize enhancer elements and identify their target genes.

FANTOM enters 20th year: expansion of transcriptomic atlases and functional annotation of non-coding RNAs.

The Functional ANnoTation Of the Mammalian genome (FANTOM) Consortium has continued to provide extensive resources in the pursuit of understanding the transcriptome, and transcriptional regulation, of mammalian genomes for the last 20 years. To share these resources with the research community, the FANTOM web-interfaces and databases are being regularly updated, enhanced and expanded with new data types. In recent years, the FANTOM Consortium's efforts have been mainly focused on creating new non-coding RNA datasets and resources. The existing FANTOM5 human and mouse miRNA atlas was supplemented with rat, dog, and chicken datasets. The sixth (latest) edition of the FANTOM project was launched to assess the function of human long non-coding RNAs (lncRNAs). From its creation until 2020, FANTOM6 has contributed to the research community a large dataset generated from the knock-down of 285 lncRNAs in human dermal fibroblasts; this is followed with extensive expression profiling and cellular phenotyping. Other updates to the FANTOM resource includes the reprocessing of the miRNA and promoter atlases of human, mouse and chicken with the latest reference genome assemblies. To facilitate the use and accessibility of all above resources we further enhanced FANTOM data viewers and web interfaces. The updated FANTOM web resource is publicly available at https://fantom.gsc.riken.jp/.

Functional annotation of human long noncoding RNAs via molecular phenotyping.

Long noncoding RNAs (lncRNAs) constitute the majority of transcripts in the mammalian genomes, and yet, their functions remain largely unknown. As part of the FANTOM6 project, we systematically knocked down the expression of 285 lncRNAs in human dermal fibroblasts and quantified cellular growth, morphological changes, and transcriptomic responses using Capped Analysis of Gene Expression (CAGE). Antisense oligonucleotides targeting the same lncRNAs exhibited global concordance, and the molecular phenotype, measured by CAGE, recapitulated the observed cellular phenotypes while providing additional insights on the affected genes and pathways. Here, we disseminate the largest-to-date lncRNA knockdown data set with molecular phenotyping (over 1000 CAGE deep-sequencing libraries) for further exploration and highlight functional roles for ZNF213-AS1 and lnc-KHDC3L-2.

The New Frontier of Functional Genomics: From Chromatin Architecture and Noncoding RNAs to Therapeutic Targets.

Common diseases are complex, multifactorial disorders whose pathogenesis is influenced by the interplay of genetic predisposition and environmental factors. Genome-wide association studies have interrogated genetic polymorphisms across genomes of individuals to test associations between genotype and susceptibility to specific disorders, providing insights into the genetic architecture of several complex disorders. However, genetic variants associated with the susceptibility to common diseases are often located in noncoding regions of the genome, such as tissue-specific enhancers or long noncoding RNAs, suggesting that regulatory elements might play a relevant role in human diseases. Enhancers are cis-regulatory genomic sequences that act in concert with promoters to regulate gene expression in a precise spatiotemporal manner. They can be located at a considerable distance from their cognate target promoters, increasing the difficulty of their identification. Genomes are organized in domains of chromatin folding, namely topologically associating domains (TADs). Identification of enhancer-promoter interactions within TADs has revealed principles of cell-type specificity across several organisms and tissues. The vast majority of mammalian genomes are pervasively transcribed, accounting for a previously unappreciated complexity of the noncoding RNA fraction. Particularly, long noncoding RNAs have emerged as key players for the establishment of chromatin architecture and regulation of gene expression. In this perspective, we describe the new advances in the fields of transcriptomics and genome organization, focusing on the role of noncoding genomic variants in the predisposition of common diseases. Finally, we propose a new framework for the identification of the next generation of pharmacological targets for common human diseases.

RADICL-seq identifies general and cell type-specific principles of genome-wide RNA-chromatin interactions.

Mammalian genomes encode tens of thousands of noncoding RNAs. Most noncoding transcripts exhibit nuclear localization and several have been shown to play a role in the regulation of gene expression and chromatin remodeling. To investigate the function of such RNAs, methods to massively map the genomic interacting sites of multiple transcripts have been developed; however, these methods have some limitations. Here, we introduce RNA And DNA Interacting Complexes Ligated and sequenced (RADICL-seq), a technology that maps genome-wide RNA-chromatin interactions in intact nuclei. RADICL-seq is a proximity ligation-based methodology that reduces the bias for nascent transcription, while increasing genomic coverage and unique mapping rate efficiency compared with existing methods. RADICL-seq identifies distinct patterns of genome occupancy for different classes of transcripts as well as cell type-specific RNA-chromatin interactions, and highlights the role of transcription in the establishment of chromatin structure.

Birth, coming of age and death: The intriguing life of long noncoding RNAs.

Mammalian genomes are pervasively transcribed, with long noncoding RNAs being the most abundant fraction. Recent studies have highlighted the central role played by these transcripts in several physiological and pathological processes. Despite several metabolic features shared between coding and noncoding transcripts, these two classes of RNAs exhibit multiple differences regarding their biogenesis and processing. Here we review such distinctions, focusing on the unique features of specific long noncoding RNAs.

FANTOM5 CAGE profiles of human and mouse samples.

In the FANTOM5 project, transcription initiation events across the human and mouse genomes were mapped at a single base-pair resolution and their frequencies were monitored by CAGE (Cap Analysis of Gene Expression) coupled with single-molecule sequencing. Approximately three thousands of samples, consisting of a variety of primary cells, tissues, cell lines, and time series samples during cell activation and development, were subjected to a uniform pipeline of CAGE data production. The analysis pipeline started by measuring RNA extracts to assess their quality, and continued to CAGE library production by using a robotic or a manual workflow, single molecule sequencing, and computational processing to generate frequencies of transcription initiation. Resulting data represents the consequence of transcriptional regulation in each analyzed state of mammalian cells. Non-overlapping peaks over the CAGE profiles, approximately 200,000 and 150,000 peaks for the human and mouse genomes, were identified and annotated to provide precise location of known promoters as well as novel ones, and to quantify their activities.

Discovery and functional analysis of lncRNAs: Methodologies to investigate an uncharacterized transcriptome.

It is known that more than 70% of mammalian genomes are transcribed, yet the vast majority of transcripts do not code for proteins. Are these noncoding transcripts merely transcriptional noise, or do they serve a biological purpose? Recent developments in genomic analysis technologies, especially sequencing methods, have allowed researchers to create a large atlas of transcriptomes, study subcellular localization, and investigate potential interactions with proteins for a growing number of transcripts. Here, we review the current methodologies available for discovering and investigating functions of long noncoding RNAs (lncRNAs), which require a wide variety of applications to study their potential biological roles. This article is part of a Special Issue entitled: Clues to long noncoding RNA taxonomy1, edited by Dr. Tetsuro Hirose and Dr. Shinichi Nakagawa.

CAGE profiling of ncRNAs in hepatocellular carcinoma reveals widespread activation of retroviral LTR promoters in virus-induced tumors.

An increasing number of noncoding RNAs (ncRNAs) have been implicated in various human diseases including cancer; however, the ncRNA transcriptome of hepatocellular carcinoma (HCC) is largely unexplored. We used CAGE to map transcription start sites across various types of human and mouse HCCs with emphasis on ncRNAs distant from protein-coding genes. Here, we report that retroviral LTR promoters, expressed in healthy tissues such as testis and placenta but not liver, are widely activated in liver tumors. Despite HCC heterogeneity, a subset of LTR-derived ncRNAs were more than 10-fold up-regulated in the vast majority of samples. HCCs with a high LTR activity mostly had a viral etiology, were less differentiated, and showed higher risk of recurrence. ChIP-seq data show that MYC and MAX are associated with ncRNA deregulation. Globally, CAGE enabled us to build a mammalian promoter map for HCC, which uncovers a new layer of complexity in HCC genomics.

Characterization of novel transcripts of human papillomavirus type 16 using cap analysis gene expression technology.

We have performed cap-analysis gene expression (CAGE) sequencing to identify the regulatory networks that orchestrate genome-wide transcription in human papillomavirus type 16 (HPV16)-positive cervical cell lines of different grades: W12E, SiHa, and CaSki. Additionally, a cervical intraepithelial neoplasia grade 1 (CIN1) lesion was assessed for identifying the transcriptome expression profile. Here we have precisely identified a novel antisense noncoding viral transcript in HPV16. In conclusion, CAGE sequencing should pave the way for understanding a diversity of viral transcript expression.

Deep transcriptome profiling of mammalian stem cells supports a regulatory role for retrotransposons in pluripotency maintenance.

The importance of microRNAs and long noncoding RNAs in the regulation of pluripotency has been documented; however, the noncoding components of stem cell gene networks remain largely unknown. Here we investigate the role of noncoding RNAs in the pluripotent state, with particular emphasis on nuclear and retrotransposon-derived transcripts. We have performed deep profiling of the nuclear and cytoplasmic transcriptomes of human and mouse stem cells, identifying a class of previously undetected stem cell-specific transcripts. We show that long terminal repeat (LTR)-derived transcripts contribute extensively to the complexity of the stem cell nuclear transcriptome. Some LTR-derived transcripts are associated with enhancer regions and are likely to be involved in the maintenance of pluripotency.

A promoter-level mammalian expression atlas.

Regulated transcription controls the diversity, developmental pathways and spatial organization of the hundreds of cell types that make up a mammal. Using single-molecule cDNA sequencing, we mapped transcription start sites (TSSs) and their usage in human and mouse primary cells, cell lines and tissues to produce a comprehensive overview of mammalian gene expression across the human body. We find that few genes are truly 'housekeeping', whereas many mammalian promoters are composite entities composed of several closely separated TSSs, with independent cell-type-specific expression profiles. TSSs specific to different cell types evolve at different rates, whereas promoters of broadly expressed genes are the most conserved. Promoter-based expression analysis reveals key transcription factors defining cell states and links them to binding-site motifs. The functions of identified novel transcripts can be predicted by coexpression and sample ontology enrichment analyses. The functional annotation of the mammalian genome 5 (FANTOM5) project provides comprehensive expression profiles and functional annotation of mammalian cell-type-specific transcriptomes with wide applications in biomedical research.

Structure by sequencing: discovery of hidden messages in the noncoding RNA fraction.

A lack of tools has hampered the search for the potential functions of the abundant noncoding RNA (ncRNA) portion of the genome. Wan et al. (2012) report a new sequencing-based method for determining ncRNA secondary structure termed PARTE, promising to accelerate ncRNA studies.

Shrinking lung syndrome successfully treated with rituximab and cyclophosphamide.

Shrinking lung syndrome (SLS) is an uncommon feature of systemic lupus erythematosus (SLE) characterized by dyspnea, pleuritic chest pain, diaphragmatic elevation, restrictive ventilatory defect and reduced respiratory muscle strength as measured by volitional tests. We report the case of a 28-year-old woman with overlapping features of SLE and Sjögren syndrome who developed severe SLS while receiving corticosteroids and azathioprine for severe polyarthritis. She was treated with a combination of rituximab and cyclophosphamide, which led to a dramatic improvement in her clinical condition and respiratory function tests. The increase in vital capacity was one of the highest among 35 published cases of SLS. Thus, restoring a near-normal lung function is an achievable goal in SLS, and the use of rituximab, with or without concomitant cyclophosphamide, certainly deserves further study in this setting.