Global analysis of mRNA isoform half-lives reveals stabilizing and destabilizing elements in yeast.

We measured half-lives of 21,248 mRNA 3' isoforms in yeast by rapidly depleting RNA polymerase II from the nucleus and performing direct RNA sequencing throughout the decay process. Interestingly, half-lives of mRNA isoforms from the same gene, including nearly identical isoforms, often vary widely. Based on clusters of isoforms with different half-lives, we identify hundreds of sequences conferring stabilization or destabilization upon mRNAs terminating downstream. One class of stabilizing element is a polyU sequence that can interact with poly(A) tails, inhibit the association of poly(A)-binding protein, and confer increased stability upon introduction into ectopic transcripts. More generally, destabilizing and stabilizing elements are linked to the propensity of the poly(A) tail to engage in double-stranded structures. Isoforms engineered to fold into 3' stem-loop structures not involving the poly(A) tail exhibit even longer half-lives. We suggest that double-stranded structures at 3' ends are a major determinant of mRNA stability.

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.

Htt CAG repeat expansion confers pleiotropic gains of mutant huntingtin function in chromatin regulation.

The CAG repeat expansion in the Huntington's disease gene HTT extends a polyglutamine tract in mutant huntingtin that enhances its ability to facilitate polycomb repressive complex 2 (PRC2). To gain insight into this dominant gain of function, we mapped histone modifications genome-wide across an isogenic panel of mouse embryonic stem cell (ESC) and neuronal progenitor cell (NPC) lines, comparing the effects of Htt null and different size Htt CAG mutations. We found that Htt is required in ESC for the proper deposition of histone H3K27me3 at a subset of 'bivalent' loci but in NPC it is needed at 'bivalent' loci for both the proper maintenance and the appropriate removal of this mark. In contrast, Htt CAG size, though changing histone H3K27me3, is prominently associated with altered histone H3K4me3 at 'active' loci. The sets of ESC and NPC genes with altered histone marks delineated by the lack of huntingtin or the presence of mutant huntingtin, though distinct, are enriched in similar pathways with apoptosis specifically highlighted for the CAG mutation. Thus, the manner by which huntingtin function facilitates PRC2 may afford mutant huntingtin with multiple opportunities to impinge upon the broader machinery that orchestrates developmentally appropriate chromatin status.

Noncoding RNAs that associate with YB-1 alter proliferation in prostate cancer cells.

The highly conserved, multifunctional YB-1 is a powerful breast cancer prognostic indicator. We report on a pervasive role for YB-1 in which it associates with thousands of nonpolyadenylated short RNAs (shyRNAs) that are further processed into small RNAs (smyRNAs). Many of these RNAs have previously been identified as functional noncoding RNAs (http://www.johnlab.org/YB1). We identified a novel, abundant, 3'-modified short RNA antisense to Dicer1 (Shad1) that colocalizes with YB-1 to P-bodies and stress granules. The expression of Shad1 was shown to correlate with that of YB-1 and whose inhibition leads to an increase in cell proliferation. Additionally, Shad1 influences the expression of additional prognostic markers of cancer progression such as DLX2 and IGFBP2. We propose that the examination of these noncoding RNAs could lead to better understanding of prostate cancer progression.

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.

Attomole-level Genomics with Single-molecule Direct DNA, cDNA and RNA Sequencing Technologies.

With the introduction of next-generation sequencing (NGS) technologies in 2005, the domination of microarrays in genomics quickly came to an end due to NGS's superior technical performance and cost advantages. By enabling genetic analysis capabilities that were not possible previously, NGS technologies have started to play an integral role in all areas of biomedical research. This chapter outlines the low-quantity DNA and cDNA sequencing capabilities and applications developed with the Helicos single molecule DNA sequencing technology.

The Lkb1 metabolic sensor maintains haematopoietic stem cell survival.

Haematopoietic stem cells (HSCs) can convert between growth states that have marked differences in bioenergetic needs. Although often quiescent in adults, these cells become proliferative upon physiological demand. Balancing HSC energetics in response to nutrient availability and growth state is poorly understood, yet essential for the dynamism of the haematopoietic system. Here we show that the Lkb1 tumour suppressor is critical for the maintenance of energy homeostasis in haematopoietic cells. Lkb1 inactivation in adult mice causes loss of HSC quiescence followed by rapid depletion of all haematopoietic subpopulations. Lkb1-deficient bone marrow cells exhibit mitochondrial defects, alterations in lipid and nucleotide metabolism, and depletion of cellular ATP. The haematopoietic effects are largely independent of Lkb1 regulation of AMP-activated protein kinase (AMPK) and mammalian target of rapamycin (mTOR) signalling. Instead, these data define a central role for Lkb1 in restricting HSC entry into cell cycle and in broadly maintaining energy homeostasis in haematopoietic cells through a novel metabolic checkpoint.

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.

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.

Oncogenic MITF dysregulation in clear cell sarcoma: defining the MiT family of human cancers.

Clear cell sarcoma (CCS) harbors a pathognomonic chromosomal translocation fusing the Ewing's sarcoma gene (EWS) to the CREB family transcription factor ATF1 and exhibits melanocytic features. We show that EWS-ATF1 occupies the MITF promoter, mimicking melanocyte-stimulating hormone (MSH) signaling to induce expression of MITF, the melanocytic master transcription factor and an amplified oncogene in melanoma. Knockdown/rescue studies revealed that MITF mediates the requirement of EWS-ATF1 for CCS survival in vitro and in vivo as well as for melanocytic differentiation. Moreover, MITF and TFE3 reciprocally rescue one another in lines derived from CCS or pediatric renal carcinoma. Seemingly unrelated tumors thus employ distinct strategies to oncogenically dysregulate the MiT family, collectively broadening the definition of MiT-associated human cancers.

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.

Prenatal AAV9-GFP administration in fetal lambs results in transduction of female germ cells and maternal exposure to virus.

Prenatal somatic cell gene therapy (PSCGT) could potentially treat severe, early-onset genetic disorders such as spinal muscular atrophy (SMA) or muscular dystrophy. Given the approval of adeno-associated virus serotype 9 (AAV9) vectors in infants with SMA by the U.S. Food and Drug Administration, we tested the safety and biodistribution of AAV9-GFP (clinical-grade and dose) in fetal lambs to understand safety and efficacy after umbilical vein or intracranial injection on embryonic day 75 (E75) . Umbilical vein injection led to widespread biodistribution of vector genomes in all examined lamb tissues and in maternal uteruses at harvest (E96 or E140; term = E150). There was robust GFP expression in brain, spinal cord, dorsal root ganglia (DRGs), without DRG toxicity and excellent transduction of diaphragm and quadriceps muscles. However, we found evidence of systemic toxicity (fetal growth restriction) and maternal exposure to the viral vector (transient elevation of total bilirubin and a trend toward elevation in anti-AAV9 antibodies). There were no antibodies against GFP in ewes or lambs. Analysis of fetal gonads demonstrated GFP expression in female (but not male) germ cells, with low levels of integration-specific reads, without integration in select proto-oncogenes. These results suggest potential therapeutic benefit of AAV9 PSCGT for neuromuscular disorders, but warrant caution for exposure of female germ cells.

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.

Gene expression profiles and differential cytoglobin expression in atrophy and adenocarcinoma of the prostate.

BACKGROUND: Proliferative inflammatory atrophy (PIA) has been proposed as a potential precursor for prostate cancer. The precise molecular abnormalities in prostatic atrophy compared to high-grade prostatic intraepithelial neoplasia (HGPIN) and carcinoma have not been fully defined. METHODS: We utilized laser capture microdissection and microarray analysis to characterize cells of PIA, HGPIN, invasive prostatic carcinoma, and non-atrophic benign prostatic epithelium (NABE). Cytoglobin was selected for immunohistochemistry (IHC) validation. IHC stains were evaluated for proportion of positive glands, and intensity of cytoglobin staining. An immunoreactive score (IR score) was determined as the product of the percentage of positive staining and intensity. RESULTS: Microarray analysis revealed probe sets that separated the microdissected cell types. Several genes showed overlapping expression patterns between PIA and PIN, and HGPIN and invasive carcinoma. Cytoglobin protein expression was detected in 57/93 (61%) of NABE and BPH cases, 92/93 atrophy (99%), 3/34 (9%) of PIN, and 23/61 carcinoma (37%) samples. The highest IHC scores were calculated for atrophy foci. A subset (33%) of atrophy cases showed the same low-cytoglobin expression level as PIN and carcinoma. CONCLUSIONS: Prostatic epithelium can be stratified into normal, atrophic, PIN, and invasive carcinoma categories based on differential genetic signatures. Cytoglobin, a protein that can be induced in response to oxidative stress, was elevated in most atrophy foci, suggesting hypoxic, and/or oxidative damage. The lower level of cytoglobin seen in neoplastic cells and 33% of atrophy foci may indicate a shared susceptibility to oxidative damage for this subset of atrophy cases and prostatic neoplasia.

The majority of total nuclear-encoded non-ribosomal RNA in a human cell is 'dark matter' un-annotated RNA.

BACKGROUND: Discovery that the transcriptional output of the human genome is far more complex than predicted by the current set of protein-coding annotations and that most RNAs produced do not appear to encode proteins has transformed our understanding of genome complexity and suggests new paradigms of genome regulation. However, the fraction of all cellular RNA whose function we do not understand and the fraction of the genome that is utilized to produce that RNA remain controversial. This is not simply a bookkeeping issue because the degree to which this un-annotated transcription is present has important implications with respect to its biologic function and to the general architecture of genome regulation. For example, efforts to elucidate how non-coding RNAs (ncRNAs) regulate genome function will be compromised if that class of RNAs is dismissed as simply 'transcriptional noise'. RESULTS: We show that the relative mass of RNA whose function and/or structure we do not understand (the so called 'dark matter' RNAs), as a proportion of all non-ribosomal, non-mitochondrial human RNA (mt-RNA), can be greater than that of protein-encoding transcripts. This observation is obscured in studies that focus only on polyA-selected RNA, a method that enriches for protein coding RNAs and at the same time discards the vast majority of RNA prior to analysis. We further show the presence of a large number of very long, abundantly-transcribed regions (100's of kb) in intergenic space and further show that expression of these regions is associated with neoplastic transformation. These overlap some regions found previously in normal human embryonic tissues and raises an interesting hypothesis as to the function of these ncRNAs in both early development and neoplastic transformation. CONCLUSIONS: We conclude that 'dark matter' RNA can constitute the majority of non-ribosomal, non-mitochondrial-RNA and a significant fraction arises from numerous very long, intergenic transcribed regions that could be involved in neoplastic transformation.

High-throughput mapping of the chromatin structure of human promoters.

Our understanding of how chromatin structure influences cellular processes such as transcription and replication has been limited by a lack of nucleosome-positioning data in human cells. We describe a high-resolution microarray approach combined with an analysis algorithm to examine nucleosome positioning in 3,692 promoters within seven human cell lines. Unlike unexpressed genes without transcription-preinitiation complexes at their promoters, expressed genes or genes containing preinitiation complexes exhibit characteristic nucleosome-free regions at their transcription start sites. The combination of these nucleosome data with chromatin immunoprecipitation-chip analyses reveals that the melanocyte master regulator microphthalmia-associated transcription factor (MITF) predominantly binds nucleosome-free regions, supporting the model that nucleosomes limit sequence accessibility. This study presents a global view of human nucleosome positioning and provides a high-throughput tool for analyzing chromatin structure in development and disease.

Quantitative polyadenylation site mapping with single-molecule direct RNA sequencing.

The known regulatory role of 3' untranslated regions (3'UTRs) and poly(A) tails in RNA localization, stability, and translation, and polyadenylation regulation defects leading to human diseases such as oculopharyngeal muscular dystrophy, thalassemias, thrombophilia, and IPEX syndrome underline the need to fully characterize genome-wide polyadenylation states and mechanisms across normal physiological and disease states. This chapter outlines the quantitative polyadenylation site mapping and analysis strategies developed with the single-molecule direct RNA sequencing technology.