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.

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.

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.

Integration of cytogenetic landmarks into the draft sequence of the human genome

We have placed 7,600 cytogenetically defined landmarks on the draft sequence of the human genome to help with the characterization of genes altered by gross chromosomal aberrations that cause human disease. The landmarks are large-insert clones mapped to chromosome bands by fluorescence in situ hybridization. Each clone contains a sequence tag that is positioned on the genomic sequence. This genome-wide set of sequence-anchored clones allows structural and functional analyses of the genome. This resource represents the first comprehensive integration of cytogenetic, radiation hybrid, linkage and sequence maps of the human genome; provides an independent validation of the sequence map and framework for contig order and orientation; surveys the genome for large-scale duplications, which are likely to require special attention during sequence assembly; and allows a stringent assessment of sequence differences between the dark and light bands of chromosomes. It also provides insight into large-scale chromatin structure and the evolution of chromosomes and gene families and will accelerate our understanding of the molecular bases of human disease and cancer.

A resource of mapped human bacterial artificial chromosome clones

To date, despite the increasing number of genomic tools, there is no repository of ordered human BAC clones that covers entire chromosomes. This project presents a resource of mapped large DNA fragments that span eight human chromosomes at approximately 1-Mb resolution. These DNA fragments are bacterial artificial chromosome (BAC) clones anchored to sequence tagged site (STS) markers. This clone collection, which currently contains 759 mapped clones, is useful in a wide range of applications from microarray-based gene mapping to identification of chromosomal mutations. In addition to the clones themselves, we describe a database, GenMapDB (http://genomics.med.upenn.edu/genmapdb), that contains information about each clone in our collection.

Comprehensive allelotyping of human renal cell carcinomas using microsatellite DNA probes

The von Hippel-Lindau locus on chromosome 3p is a tumor suppressor gene known to be involved in nonpapillary renal cell carcinoma. A previous loss of heterozygosity (LOH) study aimed at determining the allelotype of kidney tumors has indicated that in addition to 3p, chromosome arms 5q, 6q, 10q, 11q, 17p, and 19p may also harbor tumor suppressor genes. However, cytogenetic studies reveal that chromosomes 3p, 6q, 8p, 9pq, and 14q most frequently undergo karyotypic changes in renal tumors. To resolve these differences, a collection of microsatellite DNA probes has been used to scan for LOH so that 90% of individual tumor genomes were rendered informative for allele loss. The assay is capable of detecting quantitative genomic alterations in tumor cells as well. We find that LOH is most frequent for chromosome arm 3p. However, in no tumor is 3p exclusively affected. LOH for 6q, 8p, 9pq, and 14q is also distinctly elevated for both nonpapillary as well as papillary tumors and suggest that many of the tumor suppressor loci involved may be common to the etiology of both forms of kidney cancer.

Universal mapping probes and the origin of human chromosome 3

Universal mapping probes (UMPs) are defined as short segments of human DNA that are useful for physical and genetic mapping in a wide variety of mammals. The most useful UMPs contain a conserved DNA sequence immediately adjoined to a highly polymorphic CA repeat. The conserved region determines physical gene location, whereas the CA repeat facilitates genetic mapping. Both the CA repeat and its neighboring sequence are highly conserved in evolution. This permits molecular, cytogenetic, and genetic mapping of UMPs throughout mammalia. UMPs are significant because they make genetic information cumulative among well-studied species and because they transfer such information from "map rich" organisms to those that are "map poor." As a demonstration of the utility of UMPs, comparative maps between human chromosome 3 (HSA3) and the rat genome have been constructed. HSA3 is defined by at least 12 syntenic clusters located on seven different rat chromosomes. These data, together with previous comparative mapping information between human, mouse, and bovine genomes, allow us to propose a distinct evolutionary pathway that connects HSA3 with the chromosomes of rodents, artiodactyls, and primates. The model predicts a parsimonious phylogenetic tree, is readily testable, and will be of considerable use for determining the pathways of mammalian evolution.

Isolation and characterization of the gene coding for cytosolic phosphoenolpyruvate carboxykinase (GTP) from the rat

The gene for cytosolic phosphoenolpyruvate carboxykinase (GTP) [GTP:oxaloacetate carboxy-lyase (transphosphorylating), EC 4.1.1.32] from the rat was isolated from a recombinant library containing the rat genome in phage lambda Charon 4A. The isolated clone, lambda PCK1, contains the complete gene for phosphoenolpyruvate carboxykinase and approximately equal to 7 kilobases (kb) of flanking sequence at the 5' end and 1 kb at the 3' terminus. Restriction endonuclease mapping, R-loop mapping, and partial DNA sequence assay indicate that the gene is approximately equal to 6.0 kb in length (coding for a mRNA of 2.8 kb) and contains eight introns. Southern blotting of rat DNA digested with various restriction enzymes shows a pattern predicted from the restriction map of lambda PCK1. A control region at the 5' end of the gene contained in a 1.2-kb restriction fragment was isolated and subcloned into pBR322. This segment of the gene contains the usual transcription start sequences and a 24-base sequence virtually identical to the sequence found in the 5'-flanking region of the human proopiomelonocortin gene, which is known to be regulated by glucocorticoids. The 1.2-kb fragment of the phosphoenolpyruvate carboxykinase gene can be transcribed into a unique RNA fragment of predicted size by an in vitro transcription assay.