Ana Maria Suzuki

Landscape of transcription in human cells

Eukaryotic cells make many types of primary and processed RNAs that are found either in specific subcellular compartments or throughout the cells. A complete catalogue of these RNAs is not yet available and their characteristic subcellular localizations are also poorly understood. Because RNA represents the direct output of the genetic information encoded by genomes and a significant proportion of a cell’s regulatory capabilities are focused on its synthesis, processing, transport, modification and translation, the generation of such a catalogue is crucial for understanding genome function. Here we report evidence that three-quarters of the human genome is capable of being transcribed, as well as observations about the range and levels of expression, localization, processing fates, regulatory regions and modifications of almost all currently annotated and thousands of previously unannotated RNAs. These observations, taken together, prompt a redefinition of the concept of a gene.

Diversity and dynamics of the Drosophila transcriptome

Animal transcriptomes are dynamic, with each cell type, tissue and organ system expressing an ensemble of transcript isoforms that give rise to substantial diversity. Here we have identified new genes, transcripts and proteins using poly(A)+ RNA sequencing from Drosophila melanogaster in cultured cell lines, dissected organ systems and under environmental perturbations. We found that a small set of mostly neural-specific genes has the potential to encode thousands of transcripts each through extensive alternative promoter usage and RNA splicing. The magnitudes of splicing changes are larger between tissues than between developmental stages, and most sex-specific splicing is gonad-specific. Gonads express hundreds of previously unknown coding and long non-coding RNAs (lncRNAs), some of which are antisense to protein-coding genes and produce short regulatory RNAs. Furthermore, previously identified pervasive intergenic transcription occurs primarily within newly identified introns. The fly transcriptome is substantially more complex than previously recognized, with this complexity arising from combinatorial usage of promoters, splice sites and polyadenylation sites.

Two independent transcription initiation codes overlap on vertebrate core promoters

A core promoter is a stretch of DNA surrounding the transcription start site (TSS) that integrates regulatory inputs and recruits general transcription factors to initiate transcription. The nature and causative relationship of the DNA sequence and chromatin signals that govern the selection of most TSSs by RNA polymerase II remain unresolved. Maternal to zygotic transition represents the most marked change of the transcriptome repertoire in the vertebrate life cycle. Early embryonic development in zebrafish is characterized by a series of transcriptionally silent cell cycles regulated by inherited maternal gene products: zygotic genome activation commences at the tenth cell cycle, marking the mid-blastula transition. This transition provides a unique opportunity to study the rules of TSS selection and the hierarchy of events linking transcription initiation with key chromatin modifications. We analysed TSS usage during zebrafish early embryonic development at high resolution using cap analysis of gene expression, and determined the positions of H3K4me3-marked promoter-associated nucleosomes. Here we show that the transition from the maternal to zygotic transcriptome is characterized by a switch between two fundamentally different modes of defining transcription initiation, which drive the dynamic change of TSS usage and promoter shape. A maternal-specific TSS selection, which requires an A/T-rich (W-box) motif, is replaced with a zygotic TSS selection grammar characterized by broader patterns of dinucleotide enrichments, precisely aligned with the first downstream (+1) nucleosome. The developmental dynamics of the H3K4me3-marked nucleosomes reveal their DNA-sequence-associated positioning at promoters before zygotic transcription and subsequent transcription-independent adjustment to the final position downstream of the zygotic TSS. The two TSS-defining grammars coexist, often physically overlapping, in core promoters of constitutively expressed genes to enable their expression in the two regulatory environments. The dissection of overlapping core promoter determinants represents a framework for future studies of promoter structure and function across different regulatory contexts.

Comparative validation of the D. melanogaster modENCODE transcriptome annotation

Accurate gene model annotation of reference genomes is critical for making them useful. The modENCODE project has improved the D. melanogaster genome annotation by using deep and diverse high-throughput data. Since transcriptional activity that has been evolutionarily conserved is likely to have an advantageous function, we have performed large-scale interspecific comparisons to increase confidence in predicted annotations. To support comparative genomics, we filled in divergence gaps in the Drosophila phylogeny by generating draft genomes for eight new species. For comparative transcriptome analysis, we generated mRNA expression profiles on 81 samples from multiple tissues and developmental stages of 15 Drosophila species, and we performed cap analysis of gene expression in D. melanogaster and D. pseudoobscura. We also describe conservation of four distinct core promoter structures composed of combinations of elements at three positions. Overall, each type of genomic feature shows a characteristic divergence rate relative to neutral models, highlighting the value of multispecies alignment in annotating a target genome that should prove useful in the annotation of other high priority genomes, especially human and other mammalian genomes that are rich in noncoding sequences. We report that the vast majority of elements in the annotation are evolutionarily conserved, indicating that the annotation will be an important springboard for functional genetic testing by the Drosophila community.

Dynamic regulation of the transcription initiation landscape at single nucleotide resolution during vertebrate embryogenesis

Spatiotemporal control of gene expression is central to animal development. Core promoters represent a previously unanticipated regulatory level by interacting with cis-regulatory elements and transcription initiation in different physiological and developmental contexts. Here, we provide a first and comprehensive description of the core promoter repertoire and its dynamic use during the development of a vertebrate embryo. By using cap analysis of gene expression (CAGE), we mapped transcription initiation events at single nucleotide resolution across 12 stages of zebrafish development. These CAGE-based transcriptome maps reveal genome-wide rules of core promoter usage, structure, and dynamics, key to understanding the control of gene regulation during vertebrate ontogeny. They revealed the existence of multiple classes of pervasive intra- and intergenic post-transcriptionally processed RNA products and their developmental dynamics. Among these RNAs, we report splice donor site-associated intronic RNA (sRNA) to be specific to genes of the splicing machinery. For the identification of conserved features, we compared the zebrafish data sets to the first CAGE promoter map of Tetraodon and the existing human CAGE data. We show that a number of features, such as promoter type, newly discovered promoter properties such as a specialized purine-rich initiator motif, as well as sRNAs and the genes in which they are detected, are conserved in mammalian and Tetraodon CAGE-defined promoter maps. The zebrafish developmental promoterome represents a powerful resource for studying developmental gene regulation and revealing promoter features shared across vertebrates.

Suppression of artifacts and barcode bias in high-throughput transcriptome analyses utilizing template switching

Template switching (TS) has been an inherent mechanism of reverse transcriptase, which has been exploited in several transcriptome analysis methods, such as CAGE, RNA-Seq and short RNA sequencing. TS is an attractive option, given the simplicity of the protocol, which does not require an adaptor mediated step and thus minimizes sample loss. As such, it has been used in several studies that deal with limited amounts of RNA, such as in single cell studies. Additionally, TS has also been used to introduce DNA barcodes or indexes into different samples, cells or molecules. This labeling allows one to pool several samples into one sequencing flow cell, increasing the data throughput of sequencing and takes advantage of the increasing throughput of current sequences. Here, we report TS artifacts that form owing to a process called strand invasion. Due to the way in which barcodes/indexes are introduced by TS, strand invasion becomes more problematic by introducing unsystematic biases. We describe a strategy that eliminates these artifacts in silico and propose an experimental solution that suppresses biases from TS.

Digital expression profiling of the compartmentalized translatome of Purkinje neurons.

Underlying the complexity of the mammalian brain is its network of neuronal connections, but also the molecular networks of signaling pathways, protein interactions, and regulated gene expression within each individual neuron. The diversity and complexity of the spatially intermingled neurons pose a serious challenge to the identification and quantification of single neuron components. To address this challenge, we present a novel approach for the study of the ribosome-associated transcriptome-the translatome-from selected subcellular domains of specific neurons, and apply it to the Purkinje cells (PCs) in the rat cerebellum. We combined microdissection, translating ribosome affinity purification (TRAP) in nontransgenic animals, and quantitative nanoCAGE sequencing to obtain a snapshot of RNAs bound to cytoplasmic or rough endoplasmic reticulum (rER)-associated ribosomes in the PC and its dendrites. This allowed us to discover novel markers of PCs, to determine structural aspects of genes, to find hitherto uncharacterized transcripts, and to quantify biophysically relevant genes of membrane proteins controlling ion homeostasis and neuronal electrical activities.

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.

Deficiency of multidrug resistance 2 contributes to cell transformation through oxidative stress

Summary The transporter of phosphatidylcholine Mdr2/MDR3 not only plays an essential role for bile formation but also is involved in the maintenance of lipid homeostasis. Deficiency of Mdr2 leads to accumulation of ROS, cell transformation and susceptibility to intestinal carcinogenesis.

Single-nucleotide-resolution mapping of HBV promoters in infected human livers and hepatocellular carcinoma

ABSTRACT Hepatitis B virus (HBV) is a major cause of liver diseases, including hepatocellular carcinoma (HCC), and more than 650,000 people die annually due to HBV-associated liver failure. Extensive studies of individual promoters have revealed that heterogeneous RNA 5′ ends contribute to the complexity of HBV transcriptome and proteome. Here, we provide a comprehensive map of HBV transcription start sites (TSSs) in human liver, HCC, and blood, as well as several experimental replication systems, at a single-nucleotide resolution. Using CAGE (cap analysis of gene expression) analysis of 16 HCC/nontumor liver pairs, we identify 17 robust TSSs, including a novel promoter for the X gene located in the middle of the gene body, which potentially produces a shorter X protein translated from the conserved second start codon, and two minor antisense transcripts that might represent viral noncoding RNAs. Interestingly, transcription profiles were similar in HCC and nontumor livers, although quantitative analysis revealed highly variable patterns of TSS usage among clinical samples, reflecting precise regulation of HBV transcription initiation at each promoter. Unlike the variety of TSSs found in liver and HCC, the vast majority of transcripts detected in HBV-positive blood samples are pregenomic RNA, most likely generated and released from liver. Our quantitative TSS mapping using the CAGE technology will allow better understanding of HBV transcriptional responses in further studies aimed at eradicating HBV in chronic carriers. IMPORTANCE Despite the availability of a safe and effective vaccine, HBV infection remains a global health problem, and current antiviral protocols are not able to eliminate the virus in chronic carriers. Previous studies of the regulation of HBV transcription have described four major promoters and two enhancers, but little is known about their activity in human livers and HCC. We deeply sequenced the HBV RNA 5′ ends in clinical human samples and experimental models by using a new, sensitive and quantitative method termed cap analysis of gene expression (CAGE). Our data provide the first comprehensive map of global TSS distribution over the entire HBV genome in the human liver, validating already known promoters and identifying novel locations. Better knowledge of HBV transcriptional activity in the clinical setting has critical implications in the evaluation of therapeutic approaches that target HBV replication.