Nanopore-based direct sequencing of RNA transcripts with 10 different modified nucleotides reveals gaps in existing technology.

RNA undergoes complex posttranscriptional processing including chemical modifications of the nucleotides. The resultant-modified nucleotides are an integral part of RNA sequences that must be considered in studying the biology of RNA and in the design of RNA therapeutics. However, the current "RNA-sequencing" methods primarily sequence complementary DNA rather than RNA itself, which means that the modifications present in RNA are not captured in the sequencing results. Emerging direct RNA-sequencing technologies, such as those offered by Oxford Nanopore, aim to address this limitation. In this study, we synthesized and used Nanopore technology to sequence RNA transcripts consisting of canonical nucleotides and 10 different modifications in various concentrations. The results show that direct RNA sequencing still has a baseline error rate of >10%, and although some modifications can be detected, many remain unidentified. Thus, there is a need to develop sequencing technologies and analysis methods that can comprehensively capture the total complexity of RNA. The RNA sequences obtained through this project are made available for benchmarking analysis methods.

A common transcriptional mechanism involving R-loop and RNA abasic site regulates an enhancer RNA of APOE.

RNA is modified by hundreds of chemical reactions and folds into innumerable shapes. However, the regulatory role of RNA sequence and structure and how dysregulation leads to diseases remain largely unknown. Here, we uncovered a mechanism where RNA abasic sites in R-loops regulate transcription by pausing RNA polymerase II. We found an enhancer RNA, AANCR, that regulates the transcription and expression of apolipoprotein E (APOE). In some human cells such as fibroblasts, AANCR is folded into an R-loop and modified by N-glycosidic cleavage; in this form, AANCR is a partially transcribed nonfunctional enhancer and APOE is not expressed. In contrast, in other cell types including hepatocytes and under stress, AANCR does not form a stable R-loop as its sequence is not modified, so it is transcribed into a full-length enhancer that promotes APOE expression. DNA sequence variants in AANCR are associated significantly with APOE expression and Alzheimer's Disease, thus AANCR is a modifier of Alzheimer's Disease. Besides AANCR, thousands of noncoding RNAs are regulated by abasic sites in R-loops. Together our data reveal the essentiality of the folding and modification of RNA in cellular regulation and demonstrate that dysregulation underlies common complex diseases such as Alzheimer's disease.

RNA abasic sites in yeast and human cells.

RNA abasic sites and the mechanisms involved in their regulation are mostly unknown; in contrast, DNA abasic sites are well-studied. We found surprisingly that, in yeast and human cells, RNA abasic sites are prevalent. When a base is lost from RNA, the remaining ribose is found as a closed-ring or an open-ring sugar with a reactive C1' aldehyde group. Using primary amine-based reagents that react with the aldehyde group, we uncovered evidence for abasic sites in nascent RNA, messenger RNA, and ribosomal RNA from yeast and human cells. Mass spectroscopic analysis confirmed the presence of RNA abasic sites. The RNA abasic sites were found to be coupled to R-loops. We show that human methylpurine DNA glycosylase cleaves N-glycosidic bonds on RNA and that human apurinic/apyrimidinic endonuclease 1 incises RNA abasic sites in RNA-DNA hybrids. Our results reveal that, in yeast and human cells, there are RNA abasic sites, and we identify a glycosylase that generates these sites and an AP endonuclease that processes them.

cis Elements that Mediate RNA Polymerase II Pausing Regulate Human Gene Expression.

Aberrant gene expression underlies many human diseases. RNA polymerase II (Pol II) pausing is a key regulatory step in transcription. Here, we mapped the locations of RNA Pol II in normal human cells and found that RNA Pol II pauses in a consistent manner across individuals and cell types. At more than 1,000 genes including MYO1E and SESN2, RNA Pol II pauses at precise nucleotide locations. Characterization of these sites shows that RNA Pol II pauses at GC-rich regions that are marked by a sequence motif. Sixty-five percent of the pause sites are cytosines. By differential allelic gene expression analysis, we showed in our samples and a population dataset from the Genotype-Tissue Expression (GTEx) consortium that genes with more paused polymerase have lower expression levels. Furthermore, mutagenesis of the pause sites led to a significant increase in promoter activities. Thus, our data uncover that RNA Pol II pauses precisely at sites with distinct sequence features that in turn regulate gene expression.

Senataxin Mutation Reveals How R-Loops Promote Transcription by Blocking DNA Methylation at Gene Promoters.

R-loops are three-stranded nucleic acid structures found abundantly and yet often viewed as by-products of transcription. Studying cells from patients with a motor neuron disease (amyotrophic lateral sclerosis 4 [ALS4]) caused by a mutation in senataxin, we uncovered how R-loops promote transcription. In ALS4 patients, the senataxin mutation depletes R-loops with a consequent effect on gene expression. With fewer R-loops in ALS4 cells, the expression of BAMBI, a negative regulator of transforming growth factor ß (TGF-ß), is reduced; that then leads to the activation of the TGF-ß pathway. We uncovered that genome-wide R-loops influence promoter methylation of over 1,200 human genes. DNA methyl-transferase 1 favors binding to double-stranded DNA over R-loops. Thus, in forming R-loops, nascent RNA blocks DNA methylation and promotes further transcription. Hence, our results show that nucleic acid structures, in addition to sequences, influence the binding and activity of regulatory proteins.

RNA-DNA differences are generated in human cells within seconds after RNA exits polymerase II.

RNA sequences are expected to be identical to their corresponding DNA sequences. Here, we found all 12 types of RNA-DNA sequence differences (RDDs) in nascent RNA. Our results show that RDDs begin to occur in RNA chains ~55 nt from the RNA polymerase II (Pol II) active site. These RDDs occur so soon after transcription that they are incompatible with known deaminase-mediated RNA-editing mechanisms. Moreover, the 55 nt delay in appearance indicates that they do not arise during RNA synthesis by Pol II or as a direct consequence of modified base incorporation. Preliminary data suggest that RDD and R-loop formations may be coupled. These findings identify sequence substitution as an early step in cotranscriptional RNA processing.

Widespread RNA and DNA sequence differences in the human transcriptome.

The transmission of information from DNA to RNA is a critical process. We compared RNA sequences from human B cells of 27 individuals to the corresponding DNA sequences from the same individuals and uncovered more than 10,000 exonic sites where the RNA sequences do not match that of the DNA. All 12 possible categories of discordances were observed. These differences were nonrandom as many sites were found in multiple individuals and in different cell types, including primary skin cells and brain tissues. Using mass spectrometry, we detected peptides that are translated from the discordant RNA sequences and thus do not correspond exactly to the DNA sequences. These widespread RNA-DNA differences in the human transcriptome provide a yet unexplored aspect of genome variation.

Monozygotic twins reveal germline contribution to allelic expression differences.

Variation in the level of gene expression is a major determinant of a cell's function and characteristics. Common allelic variants of genes can be expressed at different levels and thus contribute to phenotypic diversity. We have measured allelic expression differences at heterozygous loci in monozygotic twins and in unrelated individuals. We show that the extent of differential allelic expression is highly similar within monozygotic twin pairs for many loci, implying that allelic differences in gene expression are under genetic control. We also show that even subtle departures from equal allelic expression are often genetically determined.

DNA reassociation using oscillating phenol emulsions.

Reassociating double-stranded DNA from single-stranded components is necessary for many molecular genetics experiments. The choice of a DNA reassociation method is dictated by the complexity of the starting material. Reassociation of simple oligomers needs only slow cooling in an aqueous environment, whereas reannealing the many single-stranded DNAs of complex genomic mixtures requires both a phenol emulsion to accelerate DNA reassociation and dedicated equipment to maintain the emulsion. We present a method that is equally suitable for reassociating either simple or complex DNA mixtures. The Oscillating Phenol Emulsion Reassociation Technique (OsPERT) was primarily developed to prepare heteroduplex DNA from alkali-denatured high molecular weight human genomic DNA samples in which hundreds of thousands of fragments need to be reannealed, but the simplicity of the technique makes it practical for less demanding DNA reassociation applications.

Direct IBD mapping: identical-by-descent mapping without genotyping.

Direct identical-by-descent (IBD) mapping is a technique, that combines genomic mismatch scanning (GMS) and DNA microarray technology, for mapping regions shared IBD between two individuals without locus-by-locus genotyping or sequencing. The lack of reagents has limited its widespread application. In particular, two key reagents have been limiting, 1). mismatch repair proteins MutS, L and H, and 2). genomic microarrays for identifying the genomic locations of the GMS-selected IBD fragments. Here, we describe steps that optimized the procedure and resources that will facilitate the development of direct IBD mapping.