Comprehensive polyadenylation site maps in yeast and human reveal pervasive alternative polyadenylation.

The emerging discoveries on the link between polyadenylation and disease states underline the need to fully characterize genome-wide polyadenylation states. Here, we report comprehensive maps of global polyadenylation events in human and yeast generated using refinements to the Direct RNA Sequencing technology. This direct approach provides a quantitative view of genome-wide polyadenylation states in a strand-specific manner and requires only attomole RNA quantities. The polyadenylation profiles revealed an abundance of unannotated polyadenylation sites, alternative polyadenylation patterns, and regulatory element-associated poly(A)(+) RNAs. We observed differences in sequence composition surrounding canonical and noncanonical human polyadenylation sites, suggesting novel noncoding RNA-specific polyadenylation mechanisms in humans. Furthermore, we observed the correlation level between sense and antisense transcripts to depend on gene expression levels, supporting the view that overlapping transcription from opposite strands may play a regulatory role. Our data provide a comprehensive view of the polyadenylation state and overlapping transcription.

RNA sequencing of pancreatic circulating tumour cells implicates WNT signalling in metastasis.

Circulating tumour cells (CTCs) shed into blood from primary cancers include putative precursors that initiate distal metastases. Although these cells are extraordinarily rare, they may identify cellular pathways contributing to the blood-borne dissemination of cancer. Here, we adapted a microfluidic device for efficient capture of CTCs from an endogenous mouse pancreatic cancer model and subjected CTCs to single-molecule RNA sequencing, identifying Wnt2 as a candidate gene enriched in CTCs. Expression of WNT2 in pancreatic cancer cells suppresses anoikis, enhances anchorage-independent sphere formation, and increases metastatic propensity in vivo. This effect is correlated with fibronectin upregulation and suppressed by inhibition of MAP3K7 (also known as TAK1) kinase. In humans, formation of non-adherent tumour spheres by pancreatic cancer cells is associated with upregulation of multiple WNT genes, and pancreatic CTCs revealed enrichment for WNT signalling in 5 out of 11 cases. Thus, molecular analysis of CTCs may identify candidate therapeutic targets to prevent the distal spread of cancer.

Deadenylase depletion protects inherited mRNAs in primordial germ cells.

A crucial event in animal development is the specification of primordial germ cells (PGCs), which become the stem cells that create sperm and eggs. How PGCs are created provides a valuable paradigm for understanding stem cells in general. We find that the PGCs of the sea urchin Strongylocentrotus purpuratus exhibit broad transcriptional repression, yet enrichment for a set of inherited mRNAs. Enrichment of several germline determinants in the PGCs requires the RNA-binding protein Nanos to target the transcript that encodes CNOT6, a deadenylase, for degradation in the PGCs, thereby creating a stable environment for RNA. Misexpression of CNOT6 in the PGCs results in their failure to retain Seawi transcripts and Vasa protein. Conversely, broad knockdown of CNOT6 expands the domain of Seawi RNA as well as exogenous reporters. Thus, Nanos-dependent spatially restricted CNOT6 differential expression is used to selectively localize germline RNAs to the PGCs. Our findings support a 'time capsule' model of germline determination, whereby the PGCs are insulated from differentiation by retaining the molecular characteristics of the totipotent egg and early embryo.

RNA sequencing: advances, challenges and opportunities.

In the few years since its initial application, massively parallel cDNA sequencing, or RNA-seq, has allowed many advances in the characterization and quantification of transcriptomes. Recently, several developments in RNA-seq methods have provided an even more complete characterization of RNA transcripts. These developments include improvements in transcription start site mapping, strand-specific measurements, gene fusion detection, small RNA characterization and detection of alternative splicing events. Ongoing developments promise further advances in the application of RNA-seq, particularly direct RNA sequencing and approaches that allow RNA quantification from very small amounts of cellular materials.

Genome-wide mapping of 5-hydroxymethylcytosine in embryonic stem cells.

5-hydroxymethylcytosine (5hmC) is a modified base present at low levels in diverse cell types in mammals. 5hmC is generated by the TET family of Fe(II) and 2-oxoglutarate-dependent enzymes through oxidation of 5-methylcytosine (5mC). 5hmC and TET proteins have been implicated in stem cell biology and cancer, but information on the genome-wide distribution of 5hmC is limited. Here we describe two novel and specific approaches to profile the genomic localization of 5hmC. The first approach, termed GLIB (glucosylation, periodate oxidation, biotinylation) uses a combination of enzymatic and chemical steps to isolate DNA fragments containing as few as a single 5hmC. The second approach involves conversion of 5hmC to cytosine 5-methylenesulphonate (CMS) by treatment of genomic DNA with sodium bisulphite, followed by immunoprecipitation of CMS-containing DNA with a specific antiserum to CMS. High-throughput sequencing of 5hmC-containing DNA from mouse embryonic stem (ES) cells showed strong enrichment within exons and near transcriptional start sites. 5hmC was especially enriched at the start sites of genes whose promoters bear dual histone 3 lysine 27 trimethylation (H3K27me3) and histone 3 lysine 4 trimethylation (H3K4me3) marks. Our results indicate that 5hmC has a probable role in transcriptional regulation, and suggest a model in which 5hmC contributes to the 'poised' chromatin signature found at developmentally-regulated genes in ES cells.

New class of gene-termini-associated human RNAs suggests a novel RNA copying mechanism.

Small (<200 nucleotide) RNA (sRNA) profiling of human cells using various technologies demonstrates unexpected complexity of sRNAs with hundreds of thousands of sRNA species present. Genetic and in vitro studies show that these RNAs are not merely degradation products of longer transcripts but could indeed have a function. Furthermore, profiling of RNAs, including the sRNAs, can reveal not only novel transcripts, but also make clear predictions about the existence and properties of novel biochemical pathways operating in a cell. For example, sRNA profiling in human cells indicated the existence of an unknown capping mechanism operating on cleaved RNA, a biochemical component of which was later identified. Here we show that human cells contain a novel type of sRNA that has non-genomically encoded 5' poly(U) tails. The presence of these RNAs at the termini of genes, specifically at the very 3' ends of known mRNAs, strongly argues for the presence of a yet uncharacterized endogenous biochemical pathway in cells that can copy RNA. We show that this pathway can operate on multiple genes, with specific enrichment towards transcript-encoding components of the translational machinery. Finally, we show that genes are also flanked by sense, 3' polyadenylated sRNAs that are likely to be capped.

New models of collaboration in genome-wide association studies: the Genetic Association Information Network.

The Genetic Association Information Network (GAIN) is a public-private partnership established to investigate the genetic basis of common diseases through a series of collaborative genome-wide association studies. GAIN has used new approaches for project selection, data deposition and distribution, collaborative analysis, publication and protection from premature intellectual property claims. These demonstrate a new commitment to shared scientific knowledge that should facilitate rapid advances in understanding the genetics of complex diseases.

Single-step capture and sequencing of natural DNA for detection of BRCA1 mutations.

Genetic testing for disease risk is an increasingly important component of medical care. However, testing can be expensive, which can lead to patients and physicians having limited access to the genetic information needed for medical decisions. To simplify DNA sample preparation and lower costs, we have developed a system in which any gene can be captured and sequenced directly from human genomic DNA without amplification, using no proteins or enzymes prior to sequencing. Extracted whole-genome DNA is acoustically sheared and loaded in a flow cell channel for single-molecule sequencing. Gene isolation, amplification, or ligation is not necessary. Accurate and low-cost detection of DNA sequence variants is demonstrated for the BRCA1 gene. Disease-causing mutations as well as common variants from well-characterized samples are identified. Single-molecule sequencing generates very reproducible coverage patterns, and these can be used to detect any size insertion or deletion directly, unlike PCR-based methods, which require additional assays. Because no gene isolation or amplification is required for sequencing, the exceptionally low costs of sample preparation and analysis could make genetic tests more accessible to those who wish to know their own disease susceptibility. Additionally, this approach has applications for sequencing integration sites for gene therapy vectors, transposons, retroviruses, and other mobile DNA elements in a more facile manner than possible with other methods.

Direct RNA sequencing.

Our understanding of human biology and disease is ultimately dependent on a complete understanding of the genome and its functions. The recent application of microarray and sequencing technologies to transcriptomics has changed the simplistic view of transcriptomes to a more complicated view of genome-wide transcription where a large fraction of transcripts emanates from unannotated parts of genomes, and underlined our limited knowledge of the dynamic state of transcription. Most of this broad body of knowledge was obtained indirectly because current transcriptome analysis methods typically require RNA to be converted to complementary DNA (cDNA) before measurements, even though the cDNA synthesis step introduces multiple biases and artefacts that interfere with both the proper characterization and quantification of transcripts. Furthermore, cDNA synthesis is not particularly suitable for the analysis of short, degraded and/or small quantity RNA samples. Here we report direct single molecule RNA sequencing without prior conversion of RNA to cDNA. We applied this technology to sequence femtomole quantities of poly(A)(+) Saccharomyces cerevisiae RNA using a surface coated with poly(dT) oligonucleotides to capture the RNAs at their natural poly(A) tails and initiate sequencing by synthesis. We observed transcript 3' end heterogeneity and polyadenylated small nucleolar RNAs. This study provides a path to high-throughput and low-cost direct RNA sequencing and achieving the ultimate goal of a comprehensive and bias-free understanding of transcriptomes.

On the importance of small changes in RNA expression.

The analysis of the differential expression of genes has been the key goal of many molecular biology methods for decades and will remain with us for decades to come. It constitutes a fundamental resource at our disposal for determining the relationship between products of transcription, biology and disease. The completed genome sequencing of many common species allowed microarrays and RNA sequencing (RNAseq) to become major tools in Systems Biology. However, we estimate that at least half of all experiments ignore transcripts that change less than some subjectively chosen threshold, typically around 2-3 fold. Here we show that a majority of the informative RNAs and differentially expressed transcripts can exhibit fold changes less than 2. We use highly quantitative single-molecule sequencing of total cellular RNA derived from a time course of inflammatory response, a process critical to a large number of diseases. Furthermore, we show that enrichment of biologically-relevant functions occurs even at very low fold changes in RNA levels. In addition, we show that most of the common statistical methods can reliably detect transcripts with low fold change when as few as 3 biological replicates are sequenced using single-molecule based RNAseq. In conclusion, given the prevalence of expression profiling in current research, the loss of data in half of all expression studies results in a significant, yet needless drain on the discovery process.

An in-depth map of polyadenylation sites in cancer.

We present a comprehensive map of over 1 million polyadenylation sites and quantify their usage in major cancers and tumor cell lines using direct RNA sequencing. We built the Expression and Polyadenylation Database to enable the visualization of the polyadenylation maps in various cancers and to facilitate the discovery of novel genes and gene isoforms that are potentially important to tumorigenesis. Analyses of polyadenylation sites indicate that a large fraction (∼30%) of mRNAs contain alternative polyadenylation sites in their 3' untranslated regions, independent of the cell type. The shortest 3' untranslated region isoforms are preferentially upregulated in cancer tissues, genome-wide. Candidate targets of alternative polyadenylation-mediated upregulation of short isoforms include POLR2K, and signaling cascades of cell-cell and cell-extracellular matrix contact, particularly involving regulators of Rho GTPases. Polyadenylation maps also helped to improve 3' untranslated region annotations and identify candidate regulatory marks such as sequence motifs, H3K36Me3 and Pabpc1 that are isoform dependent and occur in a position-specific manner. In summary, these results highlight the need to go beyond monitoring only the cumulative transcript levels for a gene, to separately analysing the expression of its RNA isoforms.

Digital transcriptome profiling from attomole-level RNA samples.

Accurate profiling of minute quantities of RNA in a global manner can enable key advances in many scientific and clinical disciplines. Here, we present low-quantity RNA sequencing (LQ-RNAseq), a high-throughput sequencing-based technique allowing whole transcriptome surveys from subnanogram RNA quantities in an amplification/ligation-free manner. LQ-RNAseq involves first-strand cDNA synthesis from RNA templates, followed by 3' polyA tailing of the single-stranded cDNA products and direct single molecule sequencing. We applied LQ-RNAseq to profile S. cerevisiae polyA+ transcripts, demonstrate the reproducibility of the approach across different sample preparations and independent instrument runs, and establish the absolute quantitative power of this method through comparisons with other reported transcript profiling techniques and through utilization of RNA spike-in experiments. We demonstrate the practical application of this approach to define the transcriptional landscape of mouse embryonic and induced pluripotent stem cells, observing transcriptional differences, including over 100 genes exhibiting differential expression between these otherwise very similar stem cell populations. This amplification-independent technology, which utilizes small quantities of nucleic acid and provides quantitative measurements of cellular transcripts, enables global gene expression measurements from minute amounts of materials and offers broad utility in both basic research and translational biology for characterization of rare cells.

Amplification-free digital gene expression profiling from minute cell quantities.

Generating reliable expression profiles from minute cell quantities is critical for scientific discovery and potential clinical applications. Here we present low-quantity digital gene expression (LQ-DGE), an amplification-free approach involving capture of poly(A)(+) RNAs from cellular lysates onto poly(dT)-coated sequencing surfaces, followed by on-surface reverse transcription and sequencing. We applied LQ-DGE to profile malignant and nonmalignant mouse and human cells, demonstrating its quantitative power and potential applicability to archival specimens.

Chromatin profiling by directly sequencing small quantities of immunoprecipitated DNA.

Chromatin structure and transcription factor localization can be assayed genome-wide by sequencing genomic DNA fractionated by protein occupancy or other properties, but current technologies involve multiple steps that introduce bias and inefficiency. Here we apply a single-molecule approach to directly sequence chromatin immunoprecipitated DNA with minimal sample manipulation. This method is compatible with just 50 pg of DNA and should thus facilitate charting chromatin maps from limited cell populations.

Association of torsades de pointes with novel and known single nucleotide polymorphisms in long QT syndrome genes.

BACKGROUND: Reduction of drug-induced adverse events may be achievable through a better understanding of the underlying causes of such events. Identifying phenotypes and genotypes that allow event prediction would provide greater safety margins for new therapeutics. Torsades de pointes (TdP) is one such life-threatening adverse event and can arise from excessive lengthening of the QT interval. This study was designed to better understand the role of genetics in the development of TdP and to determine whether genotypes can be used to predict susceptibility and thus reduce adverse events. METHODS: Seven known familial long QT syndrome genes were scanned for sequence variations in 34 patients with TdP. This group of patients is the largest such cohort ever assembled for this type of analysis. The allele frequencies for novel and known polymorphisms in these patients were compared with those in healthy control subjects. RESULTS: Six novel mutations--4 in ANK2, 1 in KCNQ1, and 1 in SCN5A--were found in the patients with TdP. Two mutations were also found in 595 healthy control subjects, whereas the others were unique to patients with TdP. Two common single nucleotide polymorphisms may be associated with the risk of TdP. The entire ANK2 gene had not been screened in a population this large previously. CONCLUSIONS: Genotypes alone could not be used to completely predict susceptibility to TdP, even when used with phenotypes. The best model using genotypic and phenotypic variables was unable to predict all events. It is unclear what other risk genes or environmental effects might be necessary to predict such cases.

Orthogonal NGS for High Throughput Clinical Diagnostics.

Next generation sequencing is a transformative technology for discovering and diagnosing genetic disorders. However, high-throughput sequencing remains error-prone, necessitating variant confirmation in order to meet the exacting demands of clinical diagnostic sequencing. To address this, we devised an orthogonal, dual platform approach employing complementary target capture and sequencing chemistries to improve speed and accuracy of variant calls at a genomic scale. We combined DNA selection by bait-based hybridization followed by Illumina NextSeq reversible terminator sequencing with DNA selection by amplification followed by Ion Proton semiconductor sequencing. This approach yields genomic scale orthogonal confirmation of ~95% of exome variants. Overall variant sensitivity improves as each method covers thousands of coding exons missed by the other. We conclude that orthogonal NGS offers improvements in variant calling sensitivity when two platforms are used, better specificity for variants identified on both platforms, and greatly reduces the time and expense of Sanger follow-up, thus enabling physicians to act on genomic results more quickly.

Association between single-nucleotide polymorphisms in the endothelial lipase (LIPG) gene and high-density lipoprotein cholesterol levels.

Endothelial lipase (LIPG) is the latest addition to the triglyceride lipase family of genes that includes pancreatic lipase (PL), hepatic lipase (HL), and lipoprotein lipase (LPL). These lipolytic enzymes demonstrate both triglyceride lipase as well as phospholipase activities and are integrally involved in lipid absorption, transport, and metabolism. Several studies have demonstrated that LIPG is important for affecting lipid levels in mice but the data in humans is less complete. To more thoroughly characterize the LIPG gene, we resequenced it from an ethnically diverse population. Thirteen novel single-nucleotide polymorphisms (SNPs) were identified and seven others confirmed. High linkage disequilibrium was found among these SNPs spanning the length of the transcript, allowing interrogation of the entire gene for functional variation. Subjects with either high or low HDL cholesterol were used to investigate its association with LIPG gene variation. Associations were found with the most significant being the intronic variants C+42T/In5 and T+2864C/In8 (P=0.007 and 0.004, respectively). A trend for an association of the same SNPs with fewer myocardial infarctions (P=0.03) was also observed but was not significant after correction for multiple testing. The results of this study provide data linking variation in the human LIPG gene with HDL cholesterol levels as well as further evidence in support of LIPG as a potential target for therapeutic intervention.

Highly polymorphic repeat region in the CETP promoter induces unusual DNA structure.

Genetic variation in the human cholesteryl ester transfer protein (CETP) promoter is associated with HDL cholesterol levels and cardiovascular disease with much of the genetic variation in CETP attributed to the promoter region. In this region, there are several single nucleotide polymorphisms as well as a variable length tandem repeat located 1946 base pairs upstream of the CETP transcription start that is highly polymorphic with respect to both length and sequence. There are more than 10 different long alleles and these vary in their repeat structure. We find that the short allele of this repeat is associated with high HDL cholesterol levels in vivo (P<0.0001). In males, this association is independent of the nearby -629 polymorphism. In addition, the variable length GAAA repeat can stimulate an adjacent GGGGA repeat to form a structure that hinders DNA amplification and sequencing. This structure also has an effect in vivo as shown by orientation effects and cloning efficiency in Escherichia coli.

CETP polymorphisms associated with HDL cholesterol may differ from those associated with cardiovascular disease.

To better understand the role of cholesteryl ester transfer protein (CETP) in cardiovascular disease, nine polymorphisms spanning the gene from the upstream promoter region to beyond the 3'UTR were genotyped in 2553 individuals from multiple ethnic groups and with different cardiovascular disease profiles. The frequency of four of these SNPs varied by 40-300% between Caucasians and African Americans. SNPs in each ethnic group fell into two haploblocks with significant linkage disequilibrium within each block. SNPs in the 5' haploblock were significantly associated with HDL cholesterol while SNPs in the 3' haploblock were, at best, only weakly associated with HDL-C. One SNP in the 3' haploblock (rs1800774 in intron 12) was highly associated with history of myocardial infarction even though it was not associated with HDL-C. This association was driven by the effect in Caucasian women where 11.9% of the women with no history of MI are homozygous for the less common allele while 23.7% of those with a history of MI share this genotype. In addition, this SNP was highly associated with BMI among Caucasians (p < 0.0001). The association of HDL-C with CETP genotype was found to be independent of smoking or alcohol consumption. These results replicate some earlier findings and also help to explain some of the apparent contradictions in the literature surrounding the role of CETP in modulating HDL-C and cardiovascular disease.

Polymorphisms in the CETP gene and association with CETP mass and HDL levels.

The cholesteryl ester transfer protein (CETP) gene has been implicated in the variation of HDL levels but most studies have focused on only one or a few genetic variations. In order to properly understand the role of CETP in determining phenotype, it is necessary to examine the entire gene and all its common polymorphisms. The coding regions, adjacent introns, and proximal 5' and 3' regions were resequenced from an ethnically diverse population. Novel and previously known polymorphisms were then characterized and associations with HDL and CETP mass levels determined. The polymorphism most highly associated with CETP was 629 bp upstream of the transcription start site while the polymorphism most highly associated with HDL was a VNTR 1946 bp upstream of the transcription start site. Genetic variation in the CETP gene is associated with protective HDL levels. The ethnic diversity of some SNPs and complex interplay among them dictate careful analysis of the whole gene prior to conclusions about the role of individual polymorphisms.