Cap Analysis Of Gene Expression Research Articles

Application of a CAGE Method to an Avian Development Study.

Cap analysis of gene expression (CAGE) is a convenient approach for genome-wide identification of promoter regions at single base-pair resolution level and accurate expression estimation of the corresponding transcripts. Depending on the initial biomaterial amount and sequencing technology, different computational pipelines for data processing are available, as well as variations of the CAGE protocol that improve sensitivity and accuracy. Therefore, this chapter elucidates the key steps of sample preparation, sequencing, and data analysis via an example of a promoter expression estimation study in chicken development. We also describe the applicability of this approach for studying other avian and reptilian species.

The FANTOM5 Computation Ecosystem: Genomic Information Hub for Promoters and Active Enhancers.

The Functional Annotation of the Mammalian Genome 5 (FANTOM5) project conducted transcriptome analysis of various mammalian cell types and provided a comprehensive resource to understand transcriptome and transcriptional regulation in individual cellular states encoded in the genome.FANTOM5 used cap analysis of gene expression (CAGE) with single-molecule sequencing to map transcription start sites (TSS) and measured their expression in a diverse range of samples. The main results from FANTOM5 were published as a promoter-level mammalian expression atlas and an atlas of active enhancers across human cell types. The FANTOM5 dataset is composed of raw experimental data and the results of bioinformatics analyses. In this chapter, we give a detailed description of the content of the FANTOM5 dataset and elaborate on different computing applications developed to publish the data and enable reproducibility and discovery of new findings. We present use cases in which the FANTOM5 dataset has been reused, leading to new findings.

NanoCAGE: A Method for the Analysis of Coding and Noncoding 5'-Capped Transcriptomes.

Transcripts in all eukaryotes are characterized by the 5'-end specific cap structure in mRNAs. Cap Analysis Gene Expression or CAGE makes use of these caps to specifically obtain cDNA fragments from the 5'-end of RNA and sequences those at high throughput for transcript identification and genome-wide mapping of transcription start sites for coding and noncoding genes. Here, we provide an improved version of our nanoCAGE protocol that has been developed for preparing CAGE libraries from as little as 50 ng of total RNA within three standard working days. Key steps in library preparation have been improved over our previously published protocol to obtain libraries having a good 5'-end selection and a more equal size distribution for higher sequencing efficiency on Illumina MiSeq and HiSeq sequencers. We recommend nanoCAGE as the method of choice for transcriptome profiling projects even from limited amounts of RNA, and as the best approach for genome-wide mapping of transcription start sites within promoter regions.

Cap Analysis of Gene Expression (CAGE) Sequencing Reveals Alterations of the Transcriptional Signatures of FLT3-ITD with Secondary D835 TKD Mutations in Acute Myeloid Leukemia

Mutations in the gene encoding Fms-like tyrosine kinase-3 (FLT3) are found in approximately 30% of acute myeloid leukemia (AML) patients, with two distinct groups of mutations; the internal tandem duplication mutations in juxtamembrane domain (FLT3-ITD) and the point mutations within the tyrosine kinase domains (FLT3-TKD).The second-generation tyrosine kinase inhibitor (TKI) quizartinib which primarily targets an inactive (type II) kinase conformation, has demonstrated responses in relapsed/refractory FLT3-ITD AML patients. However, duration of response to quizartinib is limited largely by secondary mutations in FLT3-TKD, most commonly at D835 activation loop, emerging in at least 20% of patients.To elucidate the alterations of transcript signatures of quizartinib-resistant FLT3-ITD/D835 mutants in AML, we performed cap analysis of gene expression (CAGE) sequencing utilizing 11 FLT3-ITD and 6 FLT3-ITD/D835 primary AML cells with wild type NPM1. CAGE identifies and quantifies the 5’ ends of capped mRNA transcripts, which enables the identification of transcription start sites (TSS) and allows investigating promoter structures necessary for gene expression. The TSS of genes altered in AML cells with FLT3-ITD/D835 compared with FLT3-ITD were mapped, and CAGE detected upregulation of 124 genes and downregulation of 25 genes in FLT3-ITD/D835 cells (FDR < 0.05). Specifically, several genes closely linked with activation of MAPK signaling (MAPK6, TCF4), AKT/mTOR signaling (NBN) and both (TNFAIP6) were upregulated in the D835 acquired cells. Expression of inflammatory response cytokine genes IL1B and IL8 were also increased in FLT3-ITD/D835 cells.The FLT related functional pathway analysis by KEGG ontology indicated that Ras pathway with its downstream Raf/MAPK and PI3K/AKT/mTOR signaling including transcription factor NFKB are the most significantly activated pathways in FLT3-ITD/D835 cells compared to FLT3-ITD cells. No additional activation of STAT5 pathway, a direct downstream signaling of mutated FLT3, was found in FLT3-ITD/D835 cells.To assess the underlying molecular basis of transcription signature, we performed an unbiased search for enriched sequence motifs using HOMER software. HOMER revealed upregulation of 80 motifs and downregulation of 51 motifs in AML cells with FLT3-ITD/D835 compared with FLT3-ITD (q < 0.05). Motifs of AP-1, JUN, CEBP, ATF and NFKB were substantially enriched in FLT3-ITD/D835 compared to FLT3-ITD cells. ETS, SPI-1 and RUNX1 were the most common binding motifs altered (both increased and decreased) by D835 mutation acquisition. Of note, AP-1 modulates inflammatory response and is tightly regulated by MAPK.To validate the transcriptional changes and the highlighted gene ontology results of AML primary samples, we utilized paired isogenic FLT3-ITD or FLT3-ITD/D835Y transfected Ba/F3 cell lines. As predicted, type II kinase inhibitor quizartinib was active against FLT3-ITD Ba/F3 cells (IC50; 0.6nM) but not FLT3-ITD/D835Y Ba/F3 cells (IC50; > 5000nM). Type I kinase inhibitor crenolanib was highly active against both FLT3-ITD (IC50; 2.1nM) and FLT3-ITD/D835Y cells (IC50; 16.4nM). Focusing on FLT related Ras pathway, we first investigated MAPK activation in these cells. Immunoblot analysis demonstrated substantial reduction of phospho- (p-) ERK after quizartinib (5nM) or crenolanib (50nM) treatment both in FLT3-ITD and FLT3-ITD/D835Y cells, irrespective of their susceptibility. We next examined p-S6, the target of PI3K/AKT/mTOR and MAPK signaling, and observed that crenolanib but not quizartinib successfully diminished p-S6 expression in FLT3-ITD/D835Y cells. FLT3-ITD/D835 cells showed higher baseline expression of NFkB regulator p-IkB compared to FLT3-ITD cells, which was not repressed by crenolanib.In conclusion, our study showed that the acquired D835 mutation in FLT3-ITD AML cells facilitates transcriptional network activation to induce proinflammatory genes through Ras-MAPK-NF-κB/AP-1 pathway and Ras-PI3K-AKT-mTOR signaling. This transcriptional modulation promotes oncogenic pro-survival and proliferative signaling in AML cells and defines differential sensitivity against Type I and Type II kinase inhibitors. DisclosuresShah:ARIAD: Research Funding; BMS: Research Funding; Daiichi-Sankyo: Research Funding; Pfizer: Research Funding; Plexxikon: Research Funding. Konopleva:Reata Pharmaceuticals: Equity Ownership; Abbvie: Consultancy, Research Funding; Genentech: Consultancy, Research Funding; Stemline: Consultancy, Research Funding; Eli Lilly: Research Funding; Cellectis: Research Funding; Calithera: Research Funding.

A Transcriptional Switch Point During Hematopoietic Stem and Progenitor Cell Ontogeny.

During mammalian embryogenesis, hematopoietic stem and progenitor cells (HSPCs) originate from mesoderm-derived endothelial cells in the aorta-gonad-mesonephros (AGM) region and placenta (PL). Later, HSPCs expand in fetal liver (FL) and migrate to bone marrow (BM) shortly before birth. Understanding global transcriptional regulation governing HSPC emergence from embryonic stem/induced pluripotent stem cells is necessary to devise clinical applications, such as novel transplantation approaches. In this study, to assess transcriptional dynamics during development, we performed cap analysis of gene expression on 10 developmental murine HSPC populations isolated from the AGM region, PL, FL, and BM and identified 15,681 transcripts across HSPC ontogeny. We performed microarray analysis of AGM-derived HSPCs at 9.5 and 10.5 days postcoitum (dpc) and identified 40 differentially expressed genes, 23 confirmed as significantly changed by real-time polymerase chain reaction. We conclude that a transcriptional switch point occurs in HSPC ontogeny between 9.5 and 10.5 dpc in the AGM region.

Single-Nucleotide Resolution Mapping of Hepatitis B Virus Promoters in Infected Human Livers and Hepatocellular Carcinoma.

ABSTRACTHepatitis 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.

Discriminative identification of transcriptional responses of promoters and enhancers after stimulus.

Promoters and enhancers regulate the initiation of gene expression and maintenance of expression levels in spatial and temporal manner. Recent findings stemming from the Cap Analysis of Gene Expression (CAGE) demonstrate that promoters and enhancers, based on their expression profiles after stimulus, belong to different transcription response subclasses. One of the most promising biological features that might explain the difference in transcriptional response between subclasses is the local chromatin environment. We introduce a novel computational framework, PEDAL, for distinguishing effectively transcriptional profiles of promoters and enhancers using solely histone modification marks, chromatin accessibility and binding sites of transcription factors and co-activators. A case study on data from MCF-7 cell-line reveals that PEDAL can identify successfully the transcription response subclasses of promoters and enhancers from two different stimulations. Moreover, we report subsets of input markers that discriminate with minimized classification error MCF-7 promoter and enhancer transcription response subclasses. Our work provides a general computational approach for identifying effectively cell-specific and stimulation-specific promoter and enhancer transcriptional profiles, and thus, contributes to improve our understanding of transcriptional activation in human.

FARNA: knowledgebase of inferred functions of non-coding RNA transcripts.

Non-coding RNA (ncRNA) genes play a major role in control of heterogeneous cellular behavior. Yet, their functions are largely uncharacterized. Current available databases lack in-depth information of ncRNA functions across spectrum of various cells/tissues. Here, we present FARNA, a knowledgebase of inferred functions of 10,289 human ncRNA transcripts (2,734 microRNA and 7,555 long ncRNA) in 119 tissues and 177 primary cells of human. Since transcription factors (TFs) and TF co-factors (TcoFs) are crucial components of regulatory machinery for activation of gene transcription, cellular processes and diseases in which TFs and TcoFs are involved suggest functions of the transcripts they regulate. In FARNA, functions of a transcript are inferred from TFs and TcoFs whose genes co-express with the transcript controlled by these TFs and TcoFs in a considered cell/tissue. Transcripts were annotated using statistically enriched GO terms, pathways and diseases across cells/tissues based on guilt-by-association principle. Expression profiles across cells/tissues based on Cap Analysis of Gene Expression (CAGE) are provided. FARNA, having the most comprehensive function annotation of considered ncRNAs across widest spectrum of human cells/tissues, has a potential to greatly contribute to our understanding of ncRNA roles and their regulatory mechanisms in human. FARNA can be accessed at: http://cbrc.kaust.edu.sa/farna

Promoter Architecture and Sex-Specific Gene Expression in Daphnia pulex.

Large-scale transcription start site (TSS) profiling produces a high-resolution, quantitative picture of transcription initiation and core promoter locations within a genome. However, application of TSS profiling to date has largely been restricted to a small set of prominent model systems. We sought to characterize the cis-regulatory landscape of the water flea Daphnia pulex, an emerging model arthropod that reproduces both asexually (via parthenogenesis) and sexually (via meiosis). We performed Cap Analysis of Gene Expression (CAGE) with RNA isolated from D. pulex within three developmental states: sexual females, asexual females, and males. Identified TSSs were utilized to generate a "Daphnia Promoter Atlas," i.e., a catalog of active promoters across the surveyed states. Analysis of the distribution of promoters revealed evidence for widespread alternative promoter usage in D. pulex, in addition to a prominent fraction of compactly-arranged promoters in divergent orientations. We carried out de novo motif discovery using CAGE-defined TSSs and identified eight candidate core promoter motifs; this collection includes canonical promoter elements (e.g., TATA and Initiator) in addition to others lacking obvious orthologs. A comparison of promoter activities found evidence for considerable state-specific differential gene expression between states. Our work represents the first global definition of transcription initiation and promoter architecture in crustaceans. The Daphnia Promoter Atlas presented here provides a valuable resource for comparative study of cis-regulatory regions in metazoans, as well as for investigations into the circuitries that underpin meiosis and parthenogenesis.

Abstract 2897: Recurrent transcriptome alterations across multiple cancer types

Abstract Genes that are frequently deregulated in cancer are clinically attractive as both potential pan-cancer diagnostic markers and therapeutic targets. Here we compared Cap Analysis of Gene Expression (CAGE) profiles from 225 cancer cell lines and 339 corresponding primary cell samples to identify transcripts that are recurrently deregulated in a broad range of cancer types. CAGE is a 5’ sequence tag technology that globally determines transcription start sites (TSS) in the genome and their expression levels. This allowed us to assess novel aspects of the cancer transcriptome. First, we identified (at promoter resolution) hundreds of protein coding and long non-coding transcripts that are commonly de-regulated in cancer cell lines. Next, we showed that promoters that overlap repetitive elements (especially SINE/Alu and LTR/ERV1 elements) are often upregulated in cancer. In particular, a specific repeat family, REP522 (largely palindromic, unclassified interspersed repeat of ∼1.8Kb in size), was strongly enriched for the most up-regulated promoters. To our knowledge this is the first report implicating REP522 activation in cancer. Then, taking advantage of the fact that CAGE data can be used to estimate the activity of enhancers from balanced bidirectional transcription we identified 90 enhancer RNA producing regions that are recurrently activated in cancer cell lines. With ENCODE ChIA-PET data, we linked 16 of the cancer-activated enhancers to promoters of known cancer related genes. Finally, to confirm that our results are relevant to clinical tumors we performed complementary analysis in RNA-seq data from 4,055 tumors and 563 normal tissues profiled by The Cancer Genome Atlas (TCGA) and we identified a core set of pan-cancer biomarkers (of both coding and non-coding transcripts) that are recurrently perturbed in both the FANTOM5 and TCGA datasets. In summary, our extensive transcriptome analysis identified a comprehensive set of candidate biomarkers with pan-cancer potential, and extended the perspective of enhancers and repetitive elements that are recurrently activated during carcinogenesis. Citation Format: Bogumil Kaczkowski, Yuji Tanaka, Hideya Kawaji, Albin Sandelin, Robin Andersson, Masayoshi Itoh, Timo Lassmann, Yoshihide Hayashizaki, Piero Carninci, Alistair R.R. Forrest, FANTOM5 Consortium. Recurrent transcriptome alterations across multiple cancer types. [abstract]. In: Proceedings of the 107th Annual Meeting of the American Association for Cancer Research; 2016 Apr 16-20; New Orleans, LA. Philadelphia (PA): AACR; Cancer Res 2016;76(14 Suppl):Abstract nr 2897.

Global transcript structure resolution of high gene density genomes through multi-platform data integration.

Annotation of herpesvirus genomes has traditionally been undertaken through the detection of open reading frames and other genomic motifs, supplemented with sequencing of individual cDNAs. Second generation sequencing and high-density microarray studies have revealed vastly greater herpesvirus transcriptome complexity than is captured by existing annotation. The pervasive nature of overlapping transcription throughout herpesvirus genomes, however, poses substantial problems in resolving transcript structures using these methods alone. We present an approach that combines the unique attributes of Pacific Biosciences Iso-Seq long-read, Illumina short-read and deepCAGE (Cap Analysis of Gene Expression) sequencing to globally resolve polyadenylated isoform structures in replicating Epstein-Barr virus (EBV). Our method, Transcriptome Resolution through Integration of Multi-platform Data (TRIMD), identifies nearly 300 novel EBV transcripts, quadrupling the size of the annotated viral transcriptome. These findings illustrate an array of mechanisms through which EBV achieves functional diversity in its relatively small, compact genome including programmed alternative splicing (e.g. across the IR1 repeats), alternative promoter usage by LMP2 and other latency-associated transcripts, intergenic splicing at the BZLF2 locus, and antisense transcription and pervasive readthrough transcription throughout the genome.

Heat*seq: an interactive web tool for high-throughput sequencing experiment comparison with public data.

Summary: Better protocols and decreasing costs have made high-throughput sequencing experiments now accessible even to small experimental laboratories. However, comparing one or few experiments generated by an individual lab to the vast amount of relevant data freely available in the public domain might be limited due to lack of bioinformatics expertise. Though several tools, including genome browsers, allow such comparison at a single gene level, they do not provide a genome-wide view. We developed Heat*seq, a web-tool that allows genome scale comparison of high throughput experiments chromatin immuno-precipitation followed by sequencing, RNA-sequencing and Cap Analysis of Gene Expression) provided by a user, to the data in the public domain. Heat*seq currently contains over 12 000 experiments across diverse tissues and cell types in human, mouse and drosophila. Heat*seq displays interactive correlation heatmaps, with an ability to dynamically subset datasets to contextualize user experiments. High quality figures and tables are produced and can be downloaded in multiple formats.Availability and Implementation: Web application: http://www.heatstarseq.roslin.ed.ac.uk/. Source code: https://github.com/gdevailly.Contact: Guillaume.Devailly@roslin.ed.ac.uk or Anagha.Joshi@roslin.ed.ac.ukSupplementary information: Supplementary data are available at Bioinformatics online.

The rubber tree genome shows expansion of gene family associated with rubber biosynthesis.

Hevea brasiliensis Muell. Arg, a member of the family Euphorbiaceae, is the sole natural resource exploited for commercial production of high-quality natural rubber. The properties of natural rubber latex are almost irreplaceable by synthetic counterparts for many industrial applications. A paucity of knowledge on the molecular mechanisms of rubber biosynthesis in high yield traits still persists. Here we report the comprehensive genome-wide analysis of the widely planted H. brasiliensis clone, RRIM 600. The genome was assembled based on ~155-fold combined coverage with Illumina and PacBio sequence data and has a total length of 1.55 Gb with 72.5% comprising repetitive DNA sequences. A total of 84,440 high-confidence protein-coding genes were predicted. Comparative genomic analysis revealed strong synteny between H. brasiliensis and other Euphorbiaceae genomes. Our data suggest that H. brasiliensis’s capacity to produce high levels of latex can be attributed to the expansion of rubber biosynthesis-related genes in its genome and the high expression of these genes in latex. Using cap analysis gene expression data, we illustrate the tissue-specific transcription profiles of rubber biosynthesis-related genes, revealing alternative means of transcriptional regulation. Our study adds to the understanding of H. brasiliensis biology and provides valuable genomic resources for future agronomic-related improvement of the rubber tree.

Comprehensive promoter level expression quantitative trait loci analysis of the human frontal lobe.

BackgroundExpression quantitative trait loci (eQTL) analysis is a powerful method to detect correlations between gene expression and genomic variants and is widely used to interpret the biological mechanism underlying identified genome wide association studies (GWAS) risk loci. Numerous eQTL studies have been performed on different cell types and tissues of which the majority has been based on microarray technology.MethodsWe present here an eQTL analysis based on cap analysis gene expression sequencing (CAGEseq) data created from human postmortem frontal lobe tissue combined with genotypes obtained through genotyping arrays, exome sequencing, and CAGEseq. Using CAGEseq as an expression profiling technique combined with these different genotyping techniques allows measurement of the molecular effect of variants on individual transcription start sites and increases the resolution of eQTL analysis by also including the non-annotated parts of the genome.ResultsWe identified 2410 eQTLs and show that non-coding transcripts are more likely to contain an eQTL than coding transcripts, in particular antisense transcripts. We provide evidence for how previously identified GWAS loci for schizophrenia (NRGN), Parkinson’s disease, and Alzheimer’s disease (PARK16 and MAPT loci) could increase the risk for disease at a molecular level. Furthermore, we demonstrate that CAGEseq improves eQTL analysis because variants obtained from CAGEseq are highly enriched for having a functional effect and thus are an efficient method towards the identification of causal variants.ConclusionOur data contain both coding and non-coding transcripts and has the added value that we have identified eQTLs for variants directly adjacent to TSS. Future eQTL studies would benefit from combining CAGEseq with RNA sequencing for a more complete interpretation of the transcriptome and increased understanding of eQTL signals.Electronic supplementary materialThe online version of this article (doi:10.1186/s13073-016-0320-1) contains supplementary material, which is available to authorized users.

CAGEd-oPOSSUM: motif enrichment analysis from CAGE-derived TSSs

With the emergence of large-scale Cap Analysis of Gene Expression (CAGE) datasets from individual labs and the FANTOM consortium, one can now analyze the cis-regulatory regions associated with gene transcription at an unprecedented level of refinement. By coupling transcription factor binding site (TFBS) enrichment analysis with CAGE-derived genomic regions, CAGEd-oPOSSUM can identify TFs that act as key regulators of genes involved in specific mammalian cell and tissue types. The webtool allows for the analysis of CAGE-derived transcription start sites (TSSs) either provided by the user or selected from ∼1300 mammalian samples from the FANTOM5 project with pre-computed TFBS predicted with JASPAR TF binding profiles. The tool helps power insights into the regulation of genes through the study of the specific usage of TSSs within specific cell types and/or under specific conditions.Availability and Implementation: The CAGEd-oPOSUM web tool is implemented in Perl, MySQL and Apache and is available at http://cagedop.cmmt.ubc.ca/CAGEd_oPOSSUM.Contacts: anthony.mathelier@ncmm.uio.no or wyeth@cmmt.ubc.caSupplementary information: Supplementary data are available at Bioinformatics online.

Transcriptional, epigenetic and retroviral signatures identify regulatory regions involved in hematopoietic lineage commitment.

Genome-wide approaches allow investigating the molecular circuitry wiring the genetic and epigenetic programs of human somatic stem cells. Hematopoietic stem/progenitor cells (HSPC) give rise to the different blood cell types; however, the molecular basis of human hematopoietic lineage commitment is poorly characterized. Here, we define the transcriptional and epigenetic profile of human HSPC and early myeloid and erythroid progenitors by a combination of Cap Analysis of Gene Expression (CAGE), ChIP-seq and Moloney leukemia virus (MLV) integration site mapping. Most promoters and transcripts were shared by HSPC and committed progenitors, while enhancers and super-enhancers consistently changed upon differentiation, indicating that lineage commitment is essentially regulated by enhancer elements. A significant fraction of CAGE promoters differentially expressed upon commitment were novel, harbored a chromatin enhancer signature, and may identify promoters and transcribed enhancers driving cell commitment. MLV-targeted genomic regions co-mapped with cell-specific active enhancers and super-enhancers. Expression analyses, together with an enhancer functional assay, indicate that MLV integration can be used to identify bona fide developmentally regulated enhancers. Overall, this study provides an overview of transcriptional and epigenetic changes associated to HSPC lineage commitment, and a novel signature for regulatory elements involved in cell identity.

Identification and analysis of mouse non-coding RNA using transcriptome data.

Transcripts are expressed spatially and temporally and they are very complicated, precise and specific; however, most studies are focused on protein-coding related genes. Recently, massively parallel cDNA sequencing (RNA-seq) has emerged to be a new and promising tool for transcriptome research, and numbers of non-coding RNAs, especially lincRNAs, have been widely identified and well characterized as important regulators of diverse biological processes. In this study, we used ultra-deep RNA-seq data from 15 mouse tissues to study the diversity and dynamic of non-coding RNAs in mouse. Using our own criteria, we identified totally 16,249 non-coding genes (21,569 non-coding RNAs) in mouse. We annotated these non-coding RNAs by diverse properties and found non-coding RNAs are generally shorter, have fewer exons, express in lower level and are more strikingly tissue-specific compared with protein-coding genes. Moreover, these non-coding RNAs show significant enrichment with transcriptional initiation and elongation signals including histone modifications (H3K4me3, H3K27me3 and H3K36me3), RNAPII binding sites and CAGE tags. The gene set enrichment analysis (GSEA) result revealed several sets of lincRNAs associated with diverse biological processes such as immune effector process, muscle development and sexual reproduction. Taken together, this study provides a more comprehensive annotation of mouse non-coding RNAs and gives an opportunity for future functional and evolutionary study of mouse non-coding RNAs.

Large-scale profiling of signalling pathways reveals an asthma specific signature in bronchial smooth muscle cells.

BackgroundBronchial smooth muscle (BSM) cells from asthmatic patients maintain in vitro a distinct hyper-reactive (“primed”) phenotype, characterized by increased release of pro-inflammatory factors and mediators, as well as hyperplasia and/or hypertrophy. This “primed” phenotype helps to understand pathogenesis of asthma, as changes in BSM function are essential for manifestation of allergic and inflammatory responses and airway wall remodelling.ObjectiveTo identify signalling pathways in cultured primary BSMs of asthma patients and non-asthmatic subjects by genome wide profiling of differentially expressed mRNAs and activated intracellular signalling pathways (ISPs).MethodsTranscriptome profiling by cap-analysis-of-gene-expression (CAGE), which permits selection of preferentially capped mRNAs most likely to be translated into proteins, was performed in human BSM cells from asthmatic (n=8) and non-asthmatic (n=6) subjects and OncoFinder tool were then exploited for identification of ISP deregulations.ResultsCAGE revealed >600 RNAs differentially expressed in asthma vs control cells (p≤0.005), with asthma samples showing a high degree of similarity among them. Comprehensive ISP activation analysis revealed that among 269 pathways analysed, 145 (p<0.05) or 103 (p<0.01) are differentially active in asthma, with profiles that clearly characterize BSM cells of asthmatic individuals. Notably, we identified 7 clusters of coherently acting pathways functionally related to the disease, with ISPs down-regulated in asthma mostly targeting cell death-promoting pathways and up-regulated ones affecting cell growth and proliferation, inflammatory response, control of smooth muscle contraction and hypoxia-related signalization.ConclusionsThese first-time results can now be exploited toward development of novel therapeutic strategies targeting ISP signatures linked to asthma pathophysiology.

The cell wall component lipoteichoic acid of Staphylococcus aureus induces chemokine gene expression in bovine mammary epithelial cells.

Staphylococcus aureus (SA) is a major cause of bovine mastitis, but its pathogenic mechanism remains poorly understood. To evaluate the role of lipoteichoic acid (LTA) in the immune or inflammatory response of SA mastitis, we investigated the gene expression profile in bovine mammary epithelial cells stimulated with LTA alone or with formalin-killed SA (FKSA) using cap analysis of gene expression. Seven common differentially expressed genes related to immune or inflammatory mediators were up-regulated under both LTA and FKSA stimulations. Three of these genes encode chemokines (IL-8, CXCL6 and CCL2) functioning as chemoattractant molecules for neutrophils and macrophages. These results suggest that the initial inflammatory response of SA infection in mammary gland may be related with LTA induced chemokine genes.

Tor Signaling Regulates Transcription of Amino Acid Permeases through a GATA Transcription Factor Gaf1 in Fission Yeast.

In the fission yeast, two Tor isoforms, Tor1 and Tor2, oppositely regulate gene expression of amino acid permeases. To elucidate the transcriptional machinery for these regulations, here we have employed the cap analysis of gene expression (CAGE), a method of analyzing expression profiles and identifying transcriptional start sites (TSSs). The loss of Tor1 decreased, and Tor2 inhibition by its temperature sensitive mutation increased, mRNA expression of isp5 +, per1 +, put4 + and SPBPB2B2.01. In contrast, the loss of Tor1 increased, and Tor2 inhibition decreased, the expression of cat1 +. These changes were confirmed by semi-quantitative RT-PCR. These opposite effects by the loss of Tor1 and Tor2 inhibition appeared to occur evenly across multiple TSSs for the respective genes. The motif discovery analysis based on the CAGE results identified the GATA motifs as a potential cis-regulatory element for Tor-mediated regulation. In the luciferase reporter assay, the loss of Tor1 reduced, and Tor2 inhibition and nitrogen depletion increased, the activity of isp5 + promoter as well as that of a GATAAG reporter. One of the GATAAG motifs in isp5 + promoter was critical for its transcriptional activity, and a GATA transcription factor Gaf1 was critical for the activities of isp5 + promoter and the GATAAG reporter. Furthermore, Tor2 inhibition and nitrogen depletion induced nuclear localization of Gaf1 from the cytosol and its dephosphorylation. These results suggest that Tor2 inhibition, which is known to be induced by nitrogen depletion, promotes nuclear localization of Gaf1, thereby inducing isp5 + transcription through Gaf1 binding to the GATAAG motif in its promoter. Since Gaf1 was also critical for transcription of per1 + and put4 +, Tor-Gaf1 signaling may coordinate transcription of multiple amino acid permeases according to nutrient availability.