Cap Analysis Of Gene Expression Research Articles

Cap Analysis of Gene Expression (CAGE) Sequencing Reveals Alterations of the Transcript Signatures in Acute Monocytic Leukemia Cells By Fatty Acid Oxidation Inhibition

Adipocytes, differentiated from bone marrow (BM) stroma-resident mesenchymal stem cells (MSC), are the prevalent stromal cell type in adult BM that increase with aging and cause leukemia cell resistance to chemotherapy (Ehsanipour Cancer Res. 2013).We previously reported that MSC-derived adipocytes prominently inhibited the spontaneous apoptosis of co-cultured acute monocytic leukemia (AMoL) cells, a poor-prognosis subtype of AML, and that a pharmacologic fatty acid oxidation (FAO) inhibitor Etomoxir (EX) reversed the prosurvival effects of adipocytes, indicating the importance of FAO dependent metabolic alterations in AMoL survival (Tabe ASH. 2014).To assess the transcription factors responsible for AMoL cell survival in adipocyte co-culture condition and for apoptosis induction by FAO inhibition, the alterations of transcript signatures were examined by the cap analysis of gene expression (CAGE) sequencing utilizing second-generation sequencing platform (Illumina Genome Analyzer). 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 (Carninci et al. 1996). The TSS of genes altered in U937 and THP1 cells co-cultured with adipocytes in the presence or absence of EX were mapped, and the common alterations observed in both cell types were analyzed. CAGE detected upregulation of 366 genes and downregulation of 219 genes after co-culture with adipocytes (false discovery rate, < 0.05). Ingenuity Pathway Analysis (IPA) revealed that adipocyte co-culture activated the cancer associated transcription factors Myc and FOXM1, and inhibited the p53 transcription regulator IFI16 and FLT1 kinase, an upstream positive regulator of MAPK/ERK and PI3K/AKT signaling. After EX treatment, CAGE-IPA analysis implicated inhibition of the FAO initiation enzyme ACOX1, and activation of the transcription factor ATF4 (Activating Transcription Factor 4), the master coordinator of the integrated stress response (ISR), and of the nuclear receptor PPARG which controls the FAO pathway.To narrow down the specific transcription factors responsible for EX induced apoptosis in AMoL cells co-cultured with adipocytes, CAGE-mapped TTS signature was integrated with the gene expression patterns detected by RNA-seq.CAGE and RNA-seq detected 3 genes consistently upregulated by adipocyte co-culture (KLF9; a transcription factor that activates PPARg2 promoter, FKBP5; HSP90 interacting co-chaperone immunophilin protein, ATP13A2; a member of the P5 subfamily of ATPases) and downregulation of 2 genes (ANPEP; known as CD13 or Alanyl Aminopeptidase, SLC39A10; Zinc transporter which involves in lipid metabolism).EX treatment under adipocyte co-culture condition specifically upregulated 12 genes including ISR mediator ATF4 and its target gene TRIB3. The upregulation of asparagine syntheses gene ASNS, known to be induced by ATF4, was also detected by RNA-seq. Concordant with CAGE-IPA results, EX treatment upregulated lipid accumulation marker PLIN2 and PPARG target of fatty acid binding protein FABP4 likely reflecting the direct feedback of FAO inhibition. Two genes were downregulated by EX treatment (SREBP1; the lipogenic transcription factor, P2RY2; G-protein coupled receptor activated by ATP).Finally, capillary electrophoresis mass spectrometry (CE-MS) detected the upregulation of citric acid by adipocyte co-culture, which was significantly depleted after EX treatment. The EX treatment, however, increased lactic acid along with fructose 6-phosphate and glucose 6-phosphate upregulation, indicating that cellular metabolism shifts from oxidative phosphorylation and FAO to anaerobic glycolysis which is known to activate ISR.Taken together, this study demonstrates that FAO inhibition by EX activates the pro-apoptotic transcriptional program of ISR through the up-regulation of ATF4 in adipocyte co-cultured AMoL cells. The strategies targeting FAO warrant further exploration in AMoL that survives in adipocyte abundant aged adult BM. DisclosuresKonopleva:Novartis: Research Funding; AbbVie: Research Funding; Stemline: Research Funding; Calithera: Research Funding; Threshold: Research Funding.

Remodeling of retrotransposon elements during epigenetic induction of adult visual cortical plasticityby HDAC inhibitors.

BackgroundThe capacity for plasticity in the adult brain is limited by the anatomical traces laid down during early postnatal life. Removing certain molecular brakes, such as histone deacetylases (HDACs), has proven to be effective in recapitulating juvenile plasticity in the mature visual cortex (V1). We investigated the chromatin structure and transcriptional control by genome-wide sequencing of DNase I hypersensitive sites (DHSS) and cap analysis of gene expression (CAGE) libraries after HDAC inhibition by valproic acid (VPA) in adult V1.ResultsWe found that VPA reliably reactivates the critical period plasticity and induces a dramatic change of chromatin organization in V1 yielding significantly greater accessibility distant from promoters, including at enhancer regions. VPA also induces nucleosome eviction specifically from retrotransposon (in particular SINE) elements. The transiently accessible SINE elements overlap with transcription factor-binding sites of the Fox family. Mapping of transcription start site activity using CAGE revealed transcription of epigenetic and neural plasticity-regulating genes following VPA treatment, which may help to re-program the genomic landscape and reactivate plasticity in the adult cortex.ConclusionsTreatment with HDAC inhibitors increases accessibility to enhancers and repetitive elements underlying brain-specific gene expression and reactivation of visual cortical plasticity.Electronic supplementary materialThe online version of this article (doi:10.1186/s13072-015-0043-3) contains supplementary material, which is available to authorized users.

Core Promoter Plasticity Between Maize Tissues and Genotypes Contrasts with Predominance of Sharp Transcription Initiation Sites.

Core promoters are crucial for gene regulation, providing blueprints for the assembly of transcriptional machinery at transcription start sites (TSSs). Empirically, TSSs define the coordinates of core promoters and other regulatory sequences. Thus, experimental TSS identification provides an essential step in the characterization of promoters and their features. Here, we describe the application of CAGE (cap analysis of gene expression) to identify genome-wide TSSs used in root and shoot tissues of two maize (Zea mays) inbred lines (B73 and Mo17). Our studies indicate that most TSS clusters are sharp in maize, similar to mice, but distinct from Arabidopsis thaliana, Drosophila melanogaster, or zebra fish, in which a majority of genes have broad-shaped TSS clusters. We established that ∼38% of maize promoters are characterized by a broader TATA-motif consensus, and this motif is significantly enriched in genes with sharp TSSs. A noteworthy plasticity in TSS usage between tissues and inbreds was uncovered, with ∼1500 genes showing significantly different dominant TSSs, sometimes affecting protein sequence by providing alternate translation initiation codons. We experimentally characterized instances in which this differential TSS utilization results in protein isoforms with additional domains or targeted to distinct subcellular compartments. These results provide important insights into TSS selection and gene expression in an agronomically important crop.

Identification and annotation of conserved promoters and macrophage-expressed genes in the pig genome.

BackgroundThe FANTOM5 consortium used Cap Analysis of Gene Expression (CAGE) tag sequencing to produce a comprehensive atlas of promoters and enhancers within the human and mouse genomes. We reasoned that the mapping of these regulatory elements to the pig genome could provide useful annotation and evidence to support assignment of orthology.ResultsFor human transcription start sites (TSS) associated with annotated human-mouse orthologs, 17% mapped to the pig genome but not to the mouse, 10% mapped only to the mouse, and 55% mapped to both pig and mouse. Around 17% did not map to either species. The mapping percentages were lower where there was not clear orthology relationship, but in every case, mapping to pig was greater than to mouse, and the degree of homology was also greater. Combined mapping of mouse and human CAGE-defined promoters identified at least one putative conserved TSS for >16,000 protein-coding genes. About 54% of the predicted locations of regulatory elements in the pig genome were supported by CAGE and/or RNA-Seq analysis from pig macrophages.ConclusionsComparative mapping of promoters and enhancers from humans and mice can provide useful preliminary annotation of other animal genomes. The data also confirm extensive gain and loss of regulatory elements between species, and the likelihood that pigs provide a better model than mice for human gene regulation and function.Electronic supplementary materialThe online version of this article (doi:10.1186/s12864-015-2111-2) contains supplementary material, which is available to authorized users.

Mapping Mammalian Cell-type-specific Transcriptional Regulatory Networks Using KD-CAGE and ChIP-seq Data in the TC-YIK Cell Line.

Mammals are composed of hundreds of different cell types with specialized functions. Each of these cellular phenotypes are controlled by different combinations of transcription factors. Using a human non islet cell insulinoma cell line (TC-YIK) which expresses insulin and the majority of known pancreatic beta cell specific genes as an example, we describe a general approach to identify key cell-type-specific transcription factors (TFs) and their direct and indirect targets. By ranking all human TFs by their level of enriched expression in TC-YIK relative to a broad collection of samples (FANTOM5), we confirmed known key regulators of pancreatic function and development. Systematic siRNA mediated perturbation of these TFs followed by qRT-PCR revealed their interconnections with NEUROD1 at the top of the regulation hierarchy and its depletion drastically reducing insulin levels. For 15 of the TF knock-downs (KD), we then used Cap Analysis of Gene Expression (CAGE) to identify thousands of their targets genome-wide (KD-CAGE). The data confirm NEUROD1 as a key positive regulator in the transcriptional regulatory network (TRN), and ISL1, and PROX1 as antagonists. As a complimentary approach we used ChIP-seq on four of these factors to identify NEUROD1, LMX1A, PAX6, and RFX6 binding sites in the human genome. Examining the overlap between genes perturbed in the KD-CAGE experiments and genes with a ChIP-seq peak within 50 kb of their promoter, we identified direct transcriptional targets of these TFs. Integration of KD-CAGE and ChIP-seq data shows that both NEUROD1 and LMX1A work as the main transcriptional activators. In the core TRN (i.e., TF-TF only), NEUROD1 directly transcriptionally activates the pancreatic TFs HSF4, INSM1, MLXIPL, MYT1, NKX6-3, ONECUT2, PAX4, PROX1, RFX6, ST18, DACH1, and SHOX2, while LMX1A directly transcriptionally activates DACH1, SHOX2, PAX6, and PDX1. Analysis of these complementary datasets suggests the need for caution in interpreting ChIP-seq datasets. (1) A large fraction of binding sites are at distal enhancer sites and cannot be directly associated to their targets, without chromatin conformation data. (2) Many peaks may be non-functional: even when there is a peak at a promoter, the expression of the gene may not be affected in the matching perturbation experiment.

Abstract 18159: A Functional Variant for Atrial Fibrillation Affects PITX2c Expression by Interacting With TFAP2a

Background: To date, 17 loci associated with atrial fibrillation (AF) have been identified, and there are 4 regions at 4q25 around PITX2 gene strongly associated with AF. However, the functional variants at all known AF loci remain elusive. Methods and Results: Among the 4 AF loci at 4q25, we focused on one associated region of ~84 kb extending from Chr4: 111,520,926 to 111,604,727 that includes the PITX2 gene itself (defined as all SNPs with r 2 &gt; 0.6 from the sentinel SNP rs1448818). To identify potential functional variants , we analyzed DNase hypersensitivity, histone modification and mammalian conservation from ENCODE and Roadmap Epigenomics Project. Additionally, cap analysis of gene expression (CAGE) data was also investigated. We identified 13 candidate regions, 5 of which showed enhancer activity in zebrafish skeletal muscle and heart. Using luciferase assays in a mouse atrial cell line (HL-1), we then identified one SNP (rs2595104) in which the T risk allele was associated with 31% decreased activity compared to the G non-risk allele (p = 0.007). rs2595104 lies ~10 kb upstream of PITX2c transcriptional start site and is in high linkage disequilibrium (LD) (r 2 = 0.67) with rs1448818. In silico investigation of rs2595104 using UniPROBE and HaploReg database revealed differential binding of each allele to activating enhancer binding protein 2 alpha (TFAP2a). Electrophoretic mobility shift assays showed that TFAP2a bound robustly to the non-risk G allele at rs2595104 but not to the risk T allele. Chromatin immunopricipitation followed by allele-specific qPCR in human embryonic stem cell-derived cardiomyocytes (hESC-CM) further suggested that TFAP2a was preferentially recruited to the non-risk G allele at rs2595104 in an AF-relevant cell type. hESC-CMs with homozygous deletion of rs2595104-containing enhancer using CRSPR-Cas9 expressed less PITX2c (35 ± 6%, p = 0.036), compared to the wild-type hESC-CMs. Finally, mutating the G allele at rs2595104 to T allele by CRSPR-Cas9 led to 21% reduction in PITX2c level in hESC-CMs (p = 0.027). Conclusion: We have found that the AF-associated SNP rs2595104 alters PITX2c expression via interaction with TFAP2a. Such a pathway may ultimately contribute to AF susceptibility at the PITX2 locus.

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

Multidrug resistance 2 (Mdr2), also called adenosine triphosphate-binding cassette B4 (ABCB4), is the transporter of phosphatidylcholine (PC) at the canalicular membrane of mouse hepatocytes, which plays an essential role for bile formation. Mutations in human homologue MDR3 are associated with several liver diseases. Knockout of Mdr2 results in hepatic inflammation, liver fibrosis and hepatocellular carcinoma (HCC). Whereas the pathogenesis in Mdr2 (-/-) mice has been largely attributed to the toxicity of bile acids due to the absence of PC in the bile, the question of whether Mdr2 deficiency per se perturbs biological functions in the cell has been poorly addressed. As Mdr2 is expressed in many cell types, we used mouse embryonic fibroblasts (MEF) derived from Mdr2 (-/-) embryos to show that deficiency of Mdr2 increases reactive oxygen species accumulation, lipid peroxidation and DNA damage. We found that Mdr2 (-/-) MEFs undergo spontaneous transformation and that Mdr2 (-/-) mice are more susceptible to chemical carcinogen-induced intestinal tumorigenesis. Microarray analysis in Mdr2-/- MEFs and cap analysis of gene expression in Mdr2 (-/-) HCCs revealed extensively deregulated genes involved in oxidation reduction, fatty acid metabolism and lipid biosynthesis. Our findings imply a close link between Mdr2 (-/-) -associated tumorigenesis and perturbation of these biological processes and suggest potential extrahepatic functions of Mdr2/MDR3.

DeepCAGE Transcriptomics Reveal an Important Role of the Transcription Factor MAFB in the Lymphatic Endothelium.

VEGF-C/VEGFR-3 signaling plays a central role in lymphatic development, regulating the budding oflymphatic progenitor cells from embryonic veins and maintaining the expression of PROX1 during later developmental stages. However, how VEGFR-3 activation translates into target gene expression is still not completely understood. We used cap analysis of gene expression (CAGE) RNAsequencing to characterize the transcriptional changes invoked by VEGF-C in LECs and to identifythe transcription factors (TFs) involved. We found that MAFB, a TF involved in differentiation of various cell types, is rapidly induced and activated by VEGF-C. MAFB induced expression of PROX1 as well as other TFs and markers of differentiated LECs, indicating a role in the maintenance ofthemature LEC phenotype. Correspondingly, E14.5 Mafb(-/-) embryos showed impaired lymphatic patterning in the skin. This suggests that MAFB is an important TF involved in lymphangiogenesis.

Targeting Batf2 for infectious diseases and cancer.

The family members Batf, Batf2 and Batf3 belong to a class of transcription factors containing basic leucine zipper domains that regulate various immunological functions and control the development and differentiation of immune cells. Functional studies by others demonstrated a predominant role for Batf in controlling Th2 cell functions and lineage development of T lymphocytes as well as a critical role of Batf, Batf2 and Batf3 in CD8α+dendritic cell development. Moreover, Batf family member expression was measured in a vast collection of mouse and human cell types by cap analysis gene expression (CAGE), a recent developed sequencing technology, showing reasonable expression spectrum in immune cells consistent with previously published expression profiles. Batf and Batf3 were highly expressed in lymphocytes and the earlier moderately expressed in myeloid lineages. Batf2 was predominantly expressed in monocytes/macrophages. Functional studies in mice demonstrated that Batf2 has a central role in macrophage activation by regulating inflammatory responses during lipopolysaccharides stimulation and mycobacterial infection. Hence, Batf2 could be used as a biomarker and a potential host directed drug target in tuberculosis. Moreover, Batf2 act as a tumor suppressor gene and augmenting Batf2 in malignant cells might be an encouraging therapeutic treatment against cancer.

Evolution of the unspliced transcriptome.

BackgroundDespite their abundance, unspliced EST data have received little attention as a source of information on non-coding RNAs. Very little is know, therefore, about the genomic distribution of unspliced non-coding transcripts and their relationship with the much better studied regularly spliced products. In particular, their evolution has remained virtually unstudied.ResultsWe systematically study the evidence on unspliced transcripts available in EST annotation tracks for human and mouse, comprising 104,980 and 66,109 unspliced EST clusters, respectively. Roughly one third of these are located totally inside introns of known genes (TINs) and another third overlaps exonic regions (PINs). Eleven percent are “intergenic”, far away from any annotated gene. Direct evidence for the independent transcription of many PINs and TINs is obtained from CAGE tag and chromatin data. We predict more than 2000 3’UTR-associated RNA candidates for each human and mouse. Fifteen to twenty percent of the unspliced EST cluster are conserved between human and mouse. With the exception of TINs, the sequences of unspliced EST clusters evolve significantly slower than genomic background. Furthermore, like spliced lincRNAs, they show highly tissue-specific expression patterns.ConclusionsUnspliced long non-coding RNAs are an important, rapidly evolving, component of mammalian transcriptomes. Their analysis is complicated by their preferential association with complex transcribed loci that usually also harbor a plethora of spliced transcripts. Unspliced EST data, although typically disregarded in transcriptome analysis, can be used to gain insights into this rarely investigated transcriptome component. The frequently postulated connection between lack of splicing and nuclear retention and the surprising overlap of chromatin-associated transcripts suggests that this class of transcripts might be involved in chromatin organization and possibly other mechanisms of epigenetic control.Electronic supplementary materialThe online version of this article (doi:10.1186/s12862-015-0437-7) contains supplementary material, which is available to authorized users.

NanoCAGE-XL and CapFilter: an approach to genome wide identification of high confidence transcription start sites.

BackgroundIdentifying the transcription start sites (TSS) of genes is essential for characterizing promoter regions. Several protocols have been developed to capture the 5′ end of transcripts via Cap Analysis of Gene Expression (CAGE) or linker-ligation strategies such as Paired-End Analysis of Transcription Start Sites (PEAT), but often require large amounts of tissue. More recently, nanoCAGE was developed for sequencing on the Illumina GAIIx to overcome these difficulties.ResultsHere we present the first publicly available adaptation of nanoCAGE for sequencing on recent ultra-high throughput platforms such as Illumina HiSeq-2000, and CapFilter, a computational pipeline that greatly increases confidence in TSS identification. We report excellent gene coverage, reproducibility, and precision in transcription start site discovery for samples from Arabidopsis thaliana roots.ConclusionnanoCAGE-XL together with CapFilter allows for genome wide identification of high confidence transcription start sites in large eukaryotic genomes.Electronic supplementary materialThe online version of this article (doi:10.1186/s12864-015-1670-6) contains supplementary material, which is available to authorized users.

Paradigm shifts in genomics through the FANTOM projects.

Big leaps in science happen when scientists from different backgrounds interact. In the past 15 years, the FANTOM Consortium has brought together scientists from different fields to analyze and interpret genomic data produced with novel technologies, including mouse full-length cDNAs and, more recently, expression profiling at single-nucleotide resolution by cap-analysis gene expression. The FANTOM Consortium has provided the most comprehensive mouse cDNA collection for functional studies and extensive maps of the human and mouse transcriptome comprising promoters, enhancers, as well as the network of their regulatory interactions. More importantly, serendipitous observations of the FANTOM dataset led us to realize that the mammalian genome is pervasively transcribed, even from retrotransposon elements, which were previously considered junk DNA. The majority of products from the mammalian genome are long non-coding RNAs (lncRNAs), including sense-antisense, intergenic, and enhancer RNAs. While the biological function has been elucidated for some lncRNAs, more than 98 % of them remain without a known function. We argue that large-scale studies are urgently needed to address the functional role of lncRNAs.

Promoter-level expression clustering identifies time development of transcriptional regulatory cascades initiated by ErbB receptors in breast cancer cells.

The analysis of CAGE (Cap Analysis of Gene Expression) time-course has been proposed by the FANTOM5 Consortium to extend the understanding of the sequence of events facilitating cell state transition at the level of promoter regulation. To identify the most prominent transcriptional regulations induced by growth factors in human breast cancer, we apply here the Complexity Invariant Dynamic Time Warping motif EnRichment (CIDER) analysis approach to the CAGE time-course datasets of MCF-7 cells stimulated by epidermal growth factor (EGF) or heregulin (HRG). We identify a multi-level cascade of regulations rooted by the Serum Response Factor (SRF) transcription factor, connecting the MAPK-mediated transduction of the HRG stimulus to the negative regulation of the MAPK pathway by the members of the DUSP family phosphatases. The finding confirms the known primary role of FOS and FOSL1, members of AP-1 family, in shaping gene expression in response to HRG induction. Moreover, we identify a new potential regulation of DUSP5 and RARA (known to antagonize the transcriptional regulation induced by the estrogen receptors) by the activity of the AP-1 complex, specific to HRG response. The results indicate that a divergence in AP-1 regulation determines cellular changes of breast cancer cells stimulated by ErbB receptors.

Position-specific binding of FUS to nascent RNA regulates mRNA length.

More than half of all human genes produce prematurely terminated polyadenylated short mRNAs. However, the underlying mechanisms remain largely elusive. CLIP-seq (cross-linking immunoprecipitation [CLIP] combined with deep sequencing) of FUS (fused in sarcoma) in neuronal cells showed that FUS is frequently clustered around an alternative polyadenylation (APA) site of nascent RNA. ChIP-seq (chromatin immunoprecipitation [ChIP] combined with deep sequencing) of RNA polymerase II (RNAP II) demonstrated that FUS stalls RNAP II and prematurely terminates transcription. When an APA site is located upstream of an FUS cluster, FUS enhances polyadenylation by recruiting CPSF160 and up-regulates the alternative short transcript. In contrast, when an APA site is located downstream from an FUS cluster, polyadenylation is not activated, and the RNAP II-suppressing effect of FUS leads to down-regulation of the alternative short transcript. CAGE-seq (cap analysis of gene expression [CAGE] combined with deep sequencing) and PolyA-seq (a strand-specific and quantitative method for high-throughput sequencing of 3' ends of polyadenylated transcripts) revealed that position-specific regulation of mRNA lengths by FUS is operational in two-thirds of transcripts in neuronal cells, with enrichment in genes involved in synaptic activities.

Novel splice variants in the 5'UTR of Gtf2i expressed in the rat brain: alternative 5'UTRs and differential expression in the neuronal dendrites.

General transcription factor II-I (Gtf2i) is a transcription factor and one of the genes implicated in Willams-Beuren syndrome, an autism spectrum disorder. In this study, we investigated splice variants of the Gtf2i gene in both the 5'untranslated region (5'UTR) and the coding region. To search for novel 5'UTRs of Gtf2i, we utilized the cap analysis gene expression database of the mouse. We identified seven novel Gtf2i transcripts with alternatively spliced 5'UTRs in the rat brain. We also identified four novel splice variants in the coding sequence of Gtf2i. Furthermore, we identified a selective usage of certain types of 5'UTR by coding variants. In situ hybridization demonstrated a differential pattern of expression of Gtf2i mRNAs with alternatively spliced 5'UTRs among neuronal cells, and the localization of one of the variants in neuronal dendrites in the rat brain. Immunohistochemistry also demonstrated a distribution of Gtf2i-immunoreactivity in the dendrites. These results suggest multiple pathways of expression of Gtf2i gene in the brain. The expression patterns may be under the control of alternative promoters coupled to the alternative splicing in the coding region. Differential localization of mRNA to neuronal dendrites suggests spatiotemporal-specific translation at the post-synaptic sites that is involved in transfer of synaptic activity to expression of specific sets of genes in the nucleus. Gtf2i is a transcription factor and implicated in Willams-Beuren syndrome. We identified seven novel Gtf2i transcripts with alternatively spliced 5'UTRs in the rat brain. In situ hybridization demonstrated a differential expression of Gtf2i mRNAs with different 5'UTRs in somas and dendrites of neuronal cells. Differential localization of mRNA to neuronal dendrites suggests spatiotemporal-specific translation at the postsynaptic sites. The scheme shows genomic structure showing the positions of the potential transcription start tags (rDEC695, rDEC3D7, rDEC1D3, rDEC104, rDEC072 and rDEBE25). Newly identified exons (1.1-1.6) are shown with the white boxes. The distances from rDEC695-5'end are indicated in bp.

Expression of mouse Dab2ip transcript variants and gene methylation during brain development

Dab2ip (DOC-2/DAB2 interacting protein) is a RasGAP protein which shows a growth-inhibitory effect in human prostate cancer cell lines. Recent studies have shown that Dab2ip also plays an important role in regulating dendrite development and neuronal migration during brain development. In this study, we provide a more complete description of the mouse Dab2ip (mDab2ip) gene locus and examined DNA methylation and expression of Dab2ip during cerebellar development. Analysis of cDNA sequences in public databases revealed a total of 20 possible exons for mDab2ip gene, spanning over 172kb. Using Cap Analysis of Gene Expression (CAGE) data available through FANTOM5 project, we deduced five different transcription start sites for mDab2ip. Here, we characterized three different mDab2ip transcript variants beginning with exon 1. These transcripts varied by the presence or absence of exons 3 and 5, which encode a putative nuclear localization signal and the N-terminal region of a PH-domain, respectively. The 5′ region of the mDab2ip gene contains three putative CpG islands (CpG131, CpG54, and CpG85). Interestingly, CpG54 and CpG85 are localized on exons 3 and 5. Bisulfate DNA sequencing showed that methylation level of CpG54 remained constant whereas methylation of CpG85 increased during cerebellar development. Real-time PCR analysis showed that the proportion of PH-domain containing mDab2ip transcripts increased during cerebellar development, in correlation with the increase in CpG85 methylation. These data suggest that site-specific methylation of mDab2ip gene during cerebellar development may play a role in inclusion of exon 5, resulting in a Dab2ip transcript variant that encodes a full pleckstrin homology (PH) domain.

Occupancy by key transcription factors is a more accurate predictor of enhancer activity than histone modifications or chromatin accessibility.

BackgroundRegulated gene expression controls organismal development, and variation in regulatory patterns has been implicated in complex traits. Thus accurate prediction of enhancers is important for further understanding of these processes. Genome-wide measurement of epigenetic features, such as histone modifications and occupancy by transcription factors, is improving enhancer predictions, but the contribution of these features to prediction accuracy is not known. Given the importance of the hematopoietic transcription factor TAL1 for erythroid gene activation, we predicted candidate enhancers based on genomic occupancy by TAL1 and measured their activity. Contributions of multiple features to enhancer prediction were evaluated based on the results of these and other studies.ResultsTAL1-bound DNA segments were active enhancers at a high rate both in transient transfections of cultured cells (39 of 79, or 56%) and transgenic mice (43 of 66, or 65%). The level of binding signal for TAL1 or GATA1 did not help distinguish TAL1-bound DNA segments as active versus inactive enhancers, nor did the density of regulation-related histone modifications. A meta-analysis of results from this and other studies (273 tested predicted enhancers) showed that the presence of TAL1, GATA1, EP300, SMAD1, H3K4 methylation, H3K27ac, and CAGE tags at DNase hypersensitive sites gave the most accurate predictors of enhancer activity, with a success rate over 80% and a median threefold increase in activity. Chromatin accessibility assays and the histone modifications H3K4me1 and H3K27ac were sensitive for finding enhancers, but they have high false positive rates unless transcription factor occupancy is also included.ConclusionsOccupancy by key transcription factors such as TAL1, GATA1, SMAD1, and EP300, along with evidence of transcription, improves the accuracy of enhancer predictions based on epigenetic features.Electronic supplementary materialThe online version of this article (doi:10.1186/s13072-015-0009-5) contains supplementary material, which is available to authorized users.

Transcriptional dynamics reveal critical roles for non-coding RNAs in the immediate-early response.

The immediate-early response mediates cell fate in response to a variety of extracellular stimuli and is dysregulated in many cancers. However, the specificity of the response across stimuli and cell types, and the roles of non-coding RNAs are not well understood. Using a large collection of densely-sampled time series expression data we have examined the induction of the immediate-early response in unparalleled detail, across cell types and stimuli. We exploit cap analysis of gene expression (CAGE) time series datasets to directly measure promoter activities over time. Using a novel analysis method for time series data we identify transcripts with expression patterns that closely resemble the dynamics of known immediate-early genes (IEGs) and this enables a comprehensive comparative study of these genes and their chromatin state. Surprisingly, these data suggest that the earliest transcriptional responses often involve promoters generating non-coding RNAs, many of which are produced in advance of canonical protein-coding IEGs. IEGs are known to be capable of induction without de novo protein synthesis. Consistent with this, we find that the response of both protein-coding and non-coding RNA IEGs can be explained by their transcriptionally poised, permissive chromatin state prior to stimulation. We also explore the function of non-coding RNAs in the attenuation of the immediate early response in a small RNA sequencing dataset matched to the CAGE data: We identify a novel set of microRNAs responsible for the attenuation of the IEG response in an estrogen receptor positive cancer cell line. Our computational statistical method is well suited to meta-analyses as there is no requirement for transcripts to pass thresholds for significant differential expression between time points, and it is agnostic to the number of time points per dataset.

Human endogenous retroviruses sustain complex and cooperative regulation of gene-containing loci and unannotated megabase-sized regions.

BackgroundEvidence suggests that some human endogenous retroviruses and endogenous retrovirus-like repeats (here collectively ERVs) regulate the expression of neighboring genes in normal and disease states; e.g. the human globin locus is regulated by an ERV9 that coordinates long-range gene switching during hematopoiesis and activates also intergenic transcripts. While complex transcription regulation is associated with integration of certain exogenous retroviruses, comparable regulation sustained by ERVs is less understood.FindingsWe analyzed ERV transcription using ERV9 consensus sequences and publically available RNA-sequencing, chromatin immunoprecipitation with sequencing (ChIP-seq) and cap analysis gene expression (CAGE) data from ENCODE. We discovered previously undescribed and advanced transcription regulation mechanisms in several human reference cell lines. We show that regulation by ERVs involves long-ranging activations including complex RNA splicing patterns, and transcription of large unannotated regions ranging in size from several hundred kb to around 1 Mb. Moreover, regulation was found to be cooperatively sustained in some loci by multiple ERVs and also non-LTR repeats.ConclusionOur analyses show that endogenous retroviruses sustain advanced transcription regulation in human cell lines, which shows similarities to complex insertional mutagenesis effects exerted by exogenous retroviruses. By exposing previously undescribed regulation effects, this study should prove useful for understanding fundamental transcription mechanisms resulting from evolutionary acquisition of retroviral sequence in the human genome.Electronic supplementary materialThe online version of this article (doi:10.1186/s12977-015-0161-9) contains supplementary material, which is available to authorized users.

Transcriptome analysis of controlled and therapy-resistant childhood asthma reveals distinct gene expression profiles

Children with problematic severe asthma have poor disease control despite high doses of inhaled corticosteroids and additional therapy, leading to personal suffering, early deterioration of lung function, and significant consumption of health care resources. If no exacerbating factors, such as smoking or allergies, are found after extensive investigation, these children are given a diagnosis of therapy-resistant (or therapy-refractory) asthma (SA). We sought to deepen our understanding of childhood SA by analyzing gene expression and modeling the underlying regulatory transcription factor networks in peripheral blood leukocytes. Gene expression was analyzed by using Cap Analysis of Gene Expression in children with SA (n= 13), children with controlled persistent asthma (n= 15), and age-matched healthy control subjects (n= 9). Cap Analysis of Gene Expression sequencing detects the transcription start sites of known and novel mRNAs and noncoding RNAs. Sample groups could be separated by hierarchical clustering on 1305 differentially expressed transcription start sites, including 816 known genes and several novel transcripts. Ten of 13 tested novel transcripts were validated by means of RT-PCR and Sanger sequencing. Expression of RAR-related orphan receptor A(RORA), which has been linked to asthma in genome-wide association studies, was significantly upregulated in patients with SA. Gene network modeling revealed decreased glucocorticoid receptor signaling and increased activity of the mitogen-activated protein kinase and Jun kinase cascades in patients with SA. Circulating leukocytes from children with controlled asthma and those with SA have distinct gene expression profiles, demonstrating the possible development of specific molecular biomarkers and supporting the need for novel therapeutic approaches.