Genome-wide analysis of alternative promoters of human genes using a custom promoter tiling array

Evaluating the regulatory function of non-coding autism-associated single nucleotide polymorphisms on gene expression in human brain tissue

Common variants account for most of the estimated heritability associated with autism spectrum disorder (autism). Although several replicable single nucleotide polymorphisms (SNPs) for the condition have been detected using genome-wide association study (GWAS) methodologies, their pathophysiological relevance remains elusive. Examining this is complicated, however, as all detected loci are situated within non-coding regions of the genome. It is therefore likely that they possess roles of regulatory function as opposed to directly affecting gene coding sequences. To bridge the gap between SNP discovery and mechanistic insight, we applied a comprehensive bioinformatic pipeline to functionally annotate autism-associated polymorphisms and their non-coding linkage disequilibrium (i.e., non-randomly associated) partners. We identified 82 DNA variants of probable regulatory function that may contribute to autism pathogenesis. To validate these predictions, we measured the impact of 11 high-confidence candidates and their GWAS linkage disequilibrium partners on gene expression in human brain tissue from Autistic and non-Autistic donors. Although a small number of the surveyed variants exhibited measurable influence on gene expression as determined via quantitative polymerase chain reaction, these did not survive correction for multiple comparisons. Additionally, no significant genotype-by-diagnosis effects were observed for any of the SNP-gene associations. We contend that this may reflect an inability to effectively capture the modest, neurodevelopmental-specific impact of individual variants on biological dysregulation in available post-mortem tissue samples, as well as limitations in the existing autism GWAS data.

Alternative promoters in CpG depleted regions are prevalently associated with epigenetic misregulation of liver cancer transcriptomes

Transcriptional regulation is commonly governed by alternative promoters. However, the regulatory architecture in alternative and reference promoters, and how they differ, remains elusive. In 100 CAGE-seq libraries from hepatocellular carcinoma patients, here we annotate 4083 alternative promoters in 2926 multi-promoter genes, which are largely undetected in normal livers. These genes are enriched in oncogenic processes and predominantly show association with overall survival. Alternative promoters are narrow nucleosome depleted regions, CpG island depleted, and enriched for tissue-specific transcription factors. Globally tumors lose DNA methylation. We show hierarchical retention of intragenic DNA methylation with CG-poor regions rapidly losing methylation, while CG-rich regions retain it, a process mediated by differential SETD2, H3K36me3, DNMT3B, and TET1 binding. This mechanism is validated in SETD2 knockdown cells and SETD2-mutated patients. Selective DNA methylation loss in CG-poor regions makes the chromatin accessible for alternative transcription. We show alternative promoters can control tumor transcriptomes and their regulatory architecture.

Integrative analysis of Iso-Seq and RNA-seq reveals dynamic changes of alternative promoter, alternative splicing and alternative polyadenylation during Angiotensin II-induced senescence in rat primary aortic endothelial cells

In eukaryotes, alternative promoter (AP), alternative splicing (AS), and alternative polyadenylation (APA) are three crucial regulatory mechanisms that modulate message RNA (mRNA) diversity. Although AP, AS and APA are involved in diverse biological processess, whether they have dynamic changes in Angiotensin II (Ang II) induced senescence in rat primary aortic endothelial cells (RAECs), an important cellular model for studying cardiovascular disease, remains unclear. Here we integrated both PacBio single-molecule long-read isoform sequencing (Iso-Seq) and Illumina short-read RNA sequencing (RNA-seq) to analyze the changes of AP, AS and APA in Ang II-induced senescent RAECs. Iso-Seq generated 36,278 isoforms from 10,145 gene loci and 65.81% of these isoforms are novel, which were further cross-validated by public data obtained by other techonologies such as CAGE, PolyA-Seq and 3'READS. APA contributed most to novel isoforms, followed by AS and AP. Further investigation showed that AP, AS and APA could all contribute to the regulation of isoform, but AS has more dynamic changes compared to AP and APA upon Ang II stimulation. Genes undergoing AP, AS and APA in Ang II-treated cells are enriched in various pathways related to aging or senescence, suggesting that these molecular changes are involved in functional alterations during Ang II-induced senescence. Together, the present study largely improved the annotation of rat genome and revealed gene expression changes at isoform level, extending the understanding of the complexity of gene regulation in Ang II-treated RAECs, and also provided novel clues for discovering the regulatory mechanism undelying Ang II caused vascular senescence and diseases.

Hippocampal ceRNA networks from chronic intermittent ethanol vapor-exposed male mice and functional analysis of top-ranked lncRNA genes for ethanol drinking phenotypes

The molecular mechanisms regulating the development and progression of alcohol use disorder (AUD) are largely unknown. While noncoding RNAs have previously been implicated as playing key roles in AUD, long-noncoding RNA (lncRNA) remains understudied in relation to AUD. In this study, we first identified ethanol-responsive lncRNAs in the mouse hippocampus that are transcriptional network hub genes. Microarray analysis of lncRNA, miRNA, circular RNA, and protein coding gene expression in the hippocampus from chronic intermittent ethanol vapor- or air- (control) exposed mice was used to identify ethanol-responsive competing endogenous RNA (ceRNA) networks. Highly interconnected lncRNAs (genes that had the strongest overall correlation to all other dysregulated genes identified) were ranked. The top four lncRNAs were novel, previously uncharacterized genes named Gm42575, 4930413E15Rik, Gm15767, and Gm33447, hereafter referred to as Pitt1, Pitt2, Pitt3, and Pitt4, respectively. We subsequently tested the hypothesis that CRISPR/Cas9 mutagenesis of the putative promoter and first exon of these lncRNAs in C57BL/6J mice would alter ethanol drinking behavior. The Drinking in the Dark (DID) assay was used to examine binge-like drinking behavior, and the Every-Other-Day Two-Bottle Choice (EOD-2BC) assay was used to examine intermittent ethanol consumption and preference. No significant differences between control and mutant mice were observed in the DID assay. Female-specific reductions in ethanol consumption were observed in the EOD-2BC assay for Pitt1, Pitt3, and Pitt4 mutant mice compared to controls. Male-specific alterations in ethanol preference were observed for Pitt1 and Pitt2. Female-specific increases in ethanol preference were observed for Pitt3 and Pitt4. Total fluid consumption was reduced in Pitt1 and Pitt2 mutants at 15% v/v ethanol and in Pitt3 and Pitt4 at 20% v/v ethanol in females only. We conclude that all lncRNAs targeted altered ethanol drinking behavior, and that lncRNAs Pitt1, Pitt3, and Pitt4 influenced ethanol consumption in a sex-specific manner. Further research is necessary to elucidate the biological mechanisms for these effects. These findings add to the literature implicating noncoding RNAs in AUD and suggest lncRNAs also play an important regulatory role in the disease.

Inflammation drives alternative first exon usage to regulate immune genes including a novel iron-regulated isoform of Aim2

Determining the layers of gene regulation within the innate immune response is critical to our understanding of the cellular responses to infection and dysregulation in disease. We identified a conserved mechanism of gene regulation in human and mouse via changes in alternative first exon (AFE) usage following inflammation, resulting in changes to the isoforms produced. Of these AFE events, we identified 95 unannotated transcription start sites in mice using a de novo transcriptome generated by long-read native RNA-sequencing, one of which is in the cytosolic receptor for dsDNA and known inflammatory inducible gene, Aim2. We show that this unannotated AFE isoform of Aim2 is the predominant isoform expressed during inflammation and contains an iron-responsive element in its 5′UTR enabling mRNA translation to be regulated by iron levels. This work highlights the importance of examining alternative isoform changes and translational regulation in the innate immune response and uncovers novel regulatory mechanisms of Aim2.

Human Endogenous Retroviruses as Biomedicine Markers

Human endogenous retroviruses (HERVs) were formed via ancient integration of exogenous retroviruses into the human genome and are considered to be viral "fossils". The human genome is embedded with a considerable amount of HERVs, witnessing the long-term evolutionary history of the viruses and the host. Most HERVs have lost coding capability during selection but still function in terms of HERV-mediated regulation of host gene expression. In this review, we summarize the roles of HERV activation in response to viral infections and diseases, and emphasize the potential use of HERVs as biomedicine markers in the early diagnosis of diseases such as cancer, which provides a new perspective for the clinical application of HERVs.

Promoter switching in response to changing environment and elevated expression of protein-coding genes overlapping at their 5' ends

Despite the number of studies focused on sense-antisense transcription, the key question of whether such organization evolved as a regulator of gene expression or if this is only a byproduct of other regulatory processes has not been elucidated to date. In this study, protein-coding sense-antisense gene pairs were analyzed with a particular focus on pairs overlapping at their 5' ends. Analyses were performed in 73 human transcription start site libraries. The results of our studies showed that the overlap between genes is not a stable feature and depends on which TSSs are utilized in a given cell type. An analysis of gene expression did not confirm that overlap between genes causes downregulation of their expression. This observation contradicts earlier findings. In addition, we showed that the switch from one promoter to another, leading to genes overlap, may occur in response to changing environment of a cell or tissue. We also demonstrated that in transfected and cancerous cells genes overlap is observed more often in comparison with normal tissues. Moreover, utilization of overlapping promoters depends on particular state of a cell and, at least in some groups of genes, is not merely coincidental.

Cancer Epigenomics and Beyond: Advancing the Precision Oncology Paradigm

How cancers are characterized and treated has evolved over the past few decades. Major advances in genomics tools and techniques have revealed interlinked regulatory pathways of cancers with unprecedented detail. Early discoveries led to success with rationally targeted small molecules and more recently with immunomodulatory agents, setting the stage for precision oncology. However, drug resistance to every agent has thus far proven intractable, sending us back to fill the gaps in our rudimentary knowledge of tumor biology. Epigenetics is emerging as a fundamental process in every hallmark of cancer. Large-scale interrogation of the cancer epigenome continues to reveal new mechanisms of astounding complexity. In this review, I present selected experimental and clinical examples that have shaped our understanding of cancer at the molecular level. Translation of our collective erudition into revolutionary diagnostic and treatment strategies will advance the precision oncology paradigm.

Tumor-Based Genetic Testing and Familial Cancer Risk

As genetic testing on somatic tumor tissue becomes a more routine part of personalized cancer treatment, a growing opportunity arises to identify hereditary germline variants within those results. These germline results can affect future cancer screening for both patients and their family members. Finding this germline information can be complicated as a result of differences between somatic and germline testing processes, nomenclature, and outcome goals (e.g., treatment impact). The goal of this review is to highlight differences between somatic and germline testing and outline a potential guide to allow for appropriate clinical interpretation of somatic testing results in order to better facilitate genetic counseling referrals and confirmatory germline testing.

Exon-Mediated Activation of Transcription Starts

The processing of RNA transcripts from mammalian genes occurs in proximity to their transcription. Here, we describe a phenomenon affecting thousands of genes that we call exon-mediated activation of transcription starts (EMATS), in which the splicing of internal exons impacts promoter choice and the expression level of the gene. We observed that evolutionary gain of internal exons is associated with gain of new transcription start sites (TSSs) nearby and increased gene expression. Inhibiting exon splicing reduced transcription from nearby promoters, and creation of new spliced exons activated transcription from cryptic promoters. The strongest effects occurred for weak promoters located proximal and upstream of efficiently spliced exons. Together, our findings support a model in which splicing recruits transcription machinery locally to influence TSS choice and identify exon gain, loss, and regulatory change as major contributors to the evolution of alternative promoters and gene expression in mammals.

Expressional activation and functional roles of human endogenous retroviruses in cancers

Human endogenous retroviruses (HERVs) are widely believed to be remnants of ancestral germ line infections by exogenous retroviruses. Although HERVs are deemed as “nonfunctional DNAs” due to loss of most of their viral protein coding capacity during evolution as part of the human genome, cumulative evidences are showing the expressional activation and potential roles of HERVs in diseases especially cancers. Work by other researchers and us has observed the dysregulation of HERVs in cancers, identified new HERV‐related genes, and revealed their potential importance in cancer development. Here, we summarized the current knowledge on the mechanisms of the expressional activation and functional roles of HERVs, with a focus on the H family HERV (HERV‐H), in carcinogenesis. HERV expression is regulated by external chemical or physical substances and exogenous virus infection, as well as host factors such as epigenetic DNA methylation, transcription factors, cytokines, and small RNAs. Diverse roles of HERVs have been proposed by acting in the forms of noncoding RNAs, proteins, and transcriptional regulators during carcinogenesis. However, much remains to be learnt about the contributions of HERVs to human cancers. More investigation is warranted to elucidate the functions of these “fossil remnants” yet important viral DNAs in the human genome.

Selective Activation of Alternative MYC Core Promoters by Wnt-Responsive Enhancers

In Metazoans, transcription of most genes is driven by the use of multiple alternative promoters. Although the precise regulation of alternative promoters is important for proper gene expression, the mechanisms that mediates their differential utilization remains unclear. Here, we investigate how the two alternative promoters (P1, P2) that drive MYC expression are regulated. We find that P1 and P2 can be differentially regulated across cell-types and that their selective usage is largely mediated by distal regulatory sequences. Moreover, we show that in colon carcinoma cells, Wnt-responsive enhancers preferentially upregulate transcription from the P1 promoter using reporter assays and in the context of the endogenous Wnt induction. In addition, multiple enhancer deletions using CRISPR/Cas9 corroborate the regulatory specificity of P1. Finally, we show that preferential activation between Wnt-responsive enhancers and the P1 promoter is influenced by the distinct core promoter elements that are present in the MYC promoters. Taken together, our results provide new insight into how enhancers can specifically target alternative promoters and suggest that formation of these selective interactions could allow more precise combinatorial regulation of transcription initiation.

Cancer transcriptome profiling at the juncture of clinical translation

Methodological breakthroughs over the past four decades have repeatedly revolutionized transcriptome profiling. Using RNA sequencing (RNA-seq), it has now become possible to sequence and quantify the transcriptional outputs of individual cells or thousands of samples. These transcriptomes provide a link between cellular phenotypes and their molecular underpinnings, such as mutations. In the context of cancer, this link represents an opportunity to dissect the complexity and heterogeneity of tumours and to discover new biomarkers or therapeutic strategies. Here, we review the rationale, methodology and translational impact of transcriptome profiling in cancer.

Robo3A and Robo3B expression is regulated via alternative promoters and mRNA stability

Background

The transmembrane receptor family Roundabout (Robo) was described to have an essential role in the developing nervous system. Recent studies demonstrated that Robo3 shows an altered expression in rheumatoid arthritis as well as in melanoma.

Context and purpose of the study

Until today no detailed studies of the two Robo3 isoforms (Robo3A and Robo3B) and their roles in rheumatoid arthritis synovial fibroblasts, respectively malignant melanoma are available. To get a better understanding regarding the role of Robo3A and Robo3B in the molecular process of rheumatoid arthritis and melanoma the exact characterization of expression and regulation is object of this study.

Results

mRNA and protein expression of the transcriptional variants were analyzed by quantitative RT-PCR respectively western blotting and revealed particularly enhanced expression of Robo3B in rheumatoid arthritis and melanoma. Promoter assays and inhibitor studies also disclosed that there is apparently a cell- and isoform-specific regulation of the Robo3 expression. Finally, dissimilar mRNA stabilities of Robo3A and Robo3B are identified as decisive posttranscriptional gene expression control.

Conclusion

In summary, this study supported an isotype specific role of Robo3B in disease hinting to different functional roles of each isoform.

What Is the Impact of mRNA 5' TL Heterogeneity on Translational Start Site Selection and the Mammalian Cellular Phenotype?

A major determinant in the efficiency of ribosome loading onto mRNAs is the 5′ TL (transcript leader or 5′ UTR). In addition, elements within this region also impact on start site selection demonstrating that it can modulate the protein readout at both quantitative and qualitative levels. With the increasing wealth of data generated by the mining of the mammalian transcriptome, it has become evident that a genes 5′ TL is not homogeneous but actually exhibits significant heterogeneity. This arises due to the utilization of alternative promoters, and is further compounded by significant variability with regards to the precise transcriptional start sites of each (not to mention alternative splicing). Consequently, the transcript for a protein coding gene is not a unique mRNA, but in-fact a complexed quasi-species of variants whose composition may respond to the changing physiological environment of the cell. Here we examine the potential impact of these events with regards to the protein readout.

Transcriptome-Wide Analysis Reveals Modulation of Human Macrophage Inflammatory Phenotype Through Alternative Splicing

OBJECTIVE

Human macrophages can shift phenotype across the inflammatory M1 and reparative M2 spectrum in response to environmental challenges, but the mechanisms promoting inflammatory and cardiometabolic disease-associated M1 phenotypes remain incompletely understood. Alternative splicing (AS) is emerging as an important regulator of cellular function, yet its role in macrophage activation is largely unknown. We investigated the extent to which AS occurs in M1 activation within the cardiometabolic disease context and validated a functional genomic cell model for studying human macrophage-related AS events.

APPROACH AND RESULTS

From deep RNA-sequencing of resting, M1, and M2 primary human monocyte-derived macrophages, we found 3860 differentially expressed genes in M1 activation and detected 233 M1-induced AS events; the majority of AS events were cell- and M1-specific with enrichment for pathways relevant to macrophage inflammation. Using genetic variant data for 10 cardiometabolic traits, we identified 28 trait-associated variants within the genomic loci of 21 alternatively spliced genes and 15 variants within 7 differentially expressed regulatory splicing factors in M1 activation. Knockdown of 1 such splicing factor, CELF1, in primary human macrophages led to increased inflammatory response to M1 stimulation, demonstrating CELF1's potential modulation of the M1 phenotype. Finally, we demonstrated that an induced pluripotent stem cell-derived macrophage system recapitulates M1-associated AS events and provides a high-fidelity macrophage AS model.

CONCLUSIONS

AS plays a role in defining macrophage phenotype in a cell- and stimulus-specific fashion. Alternatively spliced genes and splicing factors with trait-associated variants may reveal novel pathways and targets in cardiometabolic diseases.

Unravelling the RNA-Binding Properties of SAFB Proteins in Breast Cancer Cells

Scaffold attachment factor B1 (SAFB1) and SAFB2 proteins are oestrogen (ER) corepressors that bind to and modulate ER activity through chromatin remodelling or interaction with the basal transcription machinery. SAFB proteins also have an internal RNA-recognition motif but little is known about the RNA-binding properties of SAFB1 or SAFB2. We utilised crosslinking and immunoprecipitation (iCLIP) coupled with high-throughput sequencing to enable a transcriptome-wide mapping of SAFB1 protein-RNA interactions in breast cancer MCF-7 cells. Analysis of crosslinking frequency mapped to transcript regions revealed that SAFB1 binds to coding and noncoding RNAs (ncRNAs). The highest proportion of SAFB1 crosslink sites mapped to ncRNAs, followed by intergenic regions, open reading frames (ORFs), introns, and 3′ or 5′ untranslated regions (UTR). Furthermore, we reveal that SAFB1 binds directly to RNA and its binding is particularly enriched at purine-rich sequences not dissimilar to the RNA-binding motifs for SR proteins. Using RNAi, we also show, for the first time, that single depletion of either SAFB1 or SAFB2 leads to an increase in expression of the other SAFB protein in both MCF-7 and MDA-MD231 breast cancer cells.

An expanding universe of the non-coding genome in cancer biology

Neoplastic transformation is caused by accumulation of genetic and epigenetic alterations that ultimately convert normal cells into tumor cells with uncontrolled proliferation and survival, unlimited replicative potential and invasive growth [Hanahan,D. et al. (2011) Hallmarks of cancer: the next generation. Cell, 144, 646-674]. Although the majority of the cancer studies have focused on the functions of protein-coding genes, emerging evidence has started to reveal the importance of the vast non-coding genome, which constitutes more than 98% of the human genome. A number of non-coding RNAs (ncRNAs) derived from the 'dark matter' of the human genome exhibit cancer-specific differential expression and/or genomic alterations, and it is increasingly clear that ncRNAs, including small ncRNAs and long ncRNAs (lncRNAs), play an important role in cancer development by regulating protein-coding gene expression through diverse mechanisms. In addition to ncRNAs, nearly half of the mammalian genomes consist of transposable elements, particularly retrotransposons. Once depicted as selfish genomic parasites that propagate at the expense of host fitness, retrotransposon elements could also confer regulatory complexity to the host genomes during development and disease. Reactivation of retrotransposons in cancer, while capable of causing insertional mutagenesis and genome rearrangements to promote oncogenesis, could also alter host gene expression networks to favor tumor development. Taken together, the functional significance of non-coding genome in tumorigenesis has been previously underestimated, and diverse transcripts derived from the non-coding genome could act as integral functional components of the oncogene and tumor suppressor network.

Epidermal growth-factor-induced transcript isoform variation drives mammary cell migration

Signal-induced transcript isoform variation (TIV) includes alternative promoter usage as well as alternative splicing and alternative polyadenylation of mRNA. To assess the phenotypic relevance of signal-induced TIV, we employed exon arrays and breast epithelial cells, which migrate in response to the epidermal growth factor (EGF). We show that EGF rapidly – within one hour – induces widespread TIV in a significant fraction of the transcriptome. Importantly, TIV characterizes many genes that display no differential expression upon stimulus. In addition, similar EGF-dependent changes are shared by a panel of mammary cell lines. A functional screen, which utilized isoform-specific siRNA oligonucleotides, indicated that several isoforms play essential, non-redundant roles in EGF-induced mammary cell migration. Taken together, our findings highlight the importance of TIV in the rapid evolvement of a phenotypic response to extracellular signals.

Isoform level expression profiles provide better cancer signatures than gene level expression profiles

BACKGROUND

The majority of mammalian genes generate multiple transcript variants and protein isoforms through alternative transcription and/or alternative splicing, and the dynamic changes at the transcript/isoform level between non-oncogenic and cancer cells remain largely unexplored. We hypothesized that isoform level expression profiles would be better than gene level expression profiles at discriminating between non-oncogenic and cancer cellsgene level.

METHODS

We analyzed 160 Affymetrix exon-array datasets, comprising cell lines of non-oncogenic or oncogenic tissue origins. We obtained the transcript-level and gene level expression estimates, and used unsupervised and supervised clustering algorithms to study the profile similarity between the samples at both gene and isoform levels.

RESULTS

Hierarchical clustering, based on isoform level expressions, effectively grouped the non-oncogenic and oncogenic cell lines with a virtually perfect homogeneity-grouping rate (97.5%), regardless of the tissue origin of the cell lines. However, gene levelthis rate was much lower, being 75% at best based on the gene level expressions. Statistical analyses of the difference between cancer and non-oncogenic samples identified the existence of numerous genes with differentially expressed isoforms, which otherwise were not significant at the gene level. We also found that canonical pathways of protein ubiquitination, purine metabolism, and breast-cancer regulation by stathmin1 were significantly enriched among genes thatshow differential expression at isoform level but not at gene level.

CONCLUSIONS

In summary, cancer cell lines, regardless of their tissue of origin, can be effectively discriminated from non-cancer cell lines at isoform level, but not at gene level. This study suggests the existence of an isoform signature, rather than a gene signature, which could be used to distinguish cancer cells from normal cells.

Cloning full-length transcripts and transcript variants using 5' and 3' RACE

Gene transcripts and transcript variants must be cloned to characterize gene function and regulation. However, obtaining full-length cDNAs with accurate sequences from the 5' end through to the 3' end can be challenging. Here we describe a reverse-transcriptase-based method for obtaining full-length cDNAs using the SMARTer ("Switching Mechanism At RNA Termini") RACE technology developed by Clontech. RNA is isolated from the tissue of interest and annealed to a primer (a modified oligo(dT) primer for polyA+ transcripts; random hexamers or a gene-specific primer for polyA- transcripts). A modified MMLV-reverse transcriptase uses the primer to initiate cDNA synthesis from RNA transcript(s) annealed to the primer and continues cDNA synthesis (reverse transcription) towards the 5' end of the transcript(s). Importantly, this reverse transcriptase possesses terminal transferase activity, so when it reaches the 5' end of a transcript it adds a 3-5 residue "tail" to the newly synthesized cDNA strand. Included in the reverse transcriptase reaction mix is an oligonucleotide containing a sequence tag as well as a terminal series of modified bases that anneal to the 3-5 residue tail on the newly synthesized cDNA. The reverse transcriptase proceeds from the end of the transcript onwards into the modified bases and the rest of the sequence-tagged oligo. The newly synthesized cDNA now has a sequence tag attached to it and can be used as a template for PCR, with one primer complementary to the sequence tag and the second primer specific to the gene of interest. The fragment can be cloned and sequenced or just sequenced directly. If high-quality, undegraded RNA is used, obtaining the true 5' end of a transcript is greatly enhanced. In combination with 3' RACE, full-length transcripts are easily cloned. This method provides sequence information on important regulatory regions, such as 5' and 3' UTRs and flanking regions, and is ideal for detecting transcript variants, including those with alternative transcriptional start sites, alternative splicing, and/or alternative polyadenylation.

Integrative genome-wide chromatin signature analysis using finite mixture models

Regulation of gene expression has been shown to involve not only the binding of transcription factor at target gene promoters but also the characterization of histone around which DNA is wrapped around. Some histone modification, for example di-methylated histone H3 at lysine 4 (H3K4me2), has been shown to bind to promoters and activate target genes. However, no clear pattern has been shown to predict human promoters. This paper proposed a novel quantitative approach to characterize patterns of promoter regions and predict novel and alternative promoters. We utilized high-throughput data generated using chromatin immunoprecipitation methods followed by massively parallel sequencing (ChIP-seq) technology on RNA Polymerase II (Pol-II) and H3K4me2. Common patterns of promoter regions are modeled using a mixture model involving double-exponential and uniform distributions. The fitted model obtained were then used to search for regions displaying similar patterns over the entire genome to find novel and alternative promoters. Regions with high correlations with the common patterns are identified as putative novel promoters. We used this proposed algorithm, RNA-seq data and several transcripts databases to find alternative promoters in MCF7 (normal breast cancer) cell line. We found 7,235 high-confidence regions that display the identified promoter patterns. Of these, 4,167 regions (58%) can be mapped to RefSeq regions. 2,444 regions are in a gene body or overlap with transcripts (non-coding RNAs, ESTs, and transcripts that are predicted by RNA-seq data). Some of these maybe potential alternative promoters. We also found 193 regions that map to enhancer regions (represented by androgen and estrogen receptor binding sites) and other regulatory regions such as CTCF (CCCTC binding factor) and CpG island. Around 5% (431 regions) of these correlated regions do not overlap with any transcripts or regulatory regions suggesting that these might be potential new promoters or markers for other annotation which are currently undiscovered.

Use of alternative promoters may hide genetic effects on phenotypic traits

Genome-wide association studies have identified a multitude of single-nucleotide polymorphisms (SNPs) associated with a wide spectrum of human phenotypic traits. However, the SNPs identified so far do not explain much of the expected genetic variation and they are poor predictors of the occurrence of disease. I recently advanced the hypothesis that there is person-to-person variation in the use of alternative regulatory elements (for example, gene promoters) and this new source of variation may explain in part the low genetic variation accounted for known genetic variants. In the present report a simple mathematical model is developed to explore the biological consequences of the proposed hypothesis. The model predicts that in presence of person-to-person variation in the use of alternative promoters the observable effects of genetic variants located inside promoters will be smaller than their actual effects. As a consequence, genetic variation because of those observed polymorphisms will be reduced. The present report suggests new paths of research to elucidate the genetic basis of human complex traits.

Cooperativity of stress-responsive transcription factors in core hypoxia-inducible factor binding regions

The transcriptional response driven by Hypoxia-inducible factor (HIF) is central to the adaptation to oxygen restriction. Despite recent characterization of genome-wide HIF DNA binding locations and hypoxia-regulated transcripts in different cell types, the molecular bases of HIF target selection remain unresolved. Herein, we combined multi-level experimental data and computational predictions to identify sequence motifs that may contribute to HIF target selectivity. We obtained a core set of bona fide HIF binding regions by integrating multiple HIF1 DNA binding and hypoxia expression profiling datasets. This core set exhibits evolutionarily conserved binding regions and is enriched in functional responses to hypoxia. Computational prediction of enriched transcription factor binding sites identified sequence motifs corresponding to several stress-responsive transcription factors, such as activator protein 1 (AP1), cAMP response element-binding (CREB), or CCAAT-enhancer binding protein (CEBP). Experimental validations on HIF-regulated promoters suggest a functional role of the identified motifs in modulating HIF-mediated transcription. Accordingly, transcriptional targets of these factors are over-represented in a sorted list of hypoxia-regulated genes. Altogether, our results implicate cooperativity among stress-responsive transcription factors in fine-tuning the HIF transcriptional response.

Identification of a novel, tissue-specific ABCG2 promoter expressed in pediatric acute megakaryoblastic leukemia

ABCG2 encodes a transporter protein that is associated with multidrug-resistant phenotypes in many cancers, including acute myeloid leukemia (AML); high levels of expression are generally associated with a poor prognosis. To better understand how expression of ABCG2 is controlled in pediatric AML, we performed a detailed analysis of the ABCG2 transcript isoforms from a variety of tissue sources, including 85 pediatric AML samples. These studies revealed a complex 5' untranslated region (UTR) with 6 novel exons and multiple splice variants. Samples from children with acute megakaryoblastic leukemia (AML FAB-M7) not associated with Down syndrome showed uniformly higher levels of ABCG2 transcripts than samples from children with other AML subtypes. A novel 5' UTR identified 90kb upstream of the exon 2 translation initiation site was expressed only in M7 AML subtypes. An associated upstream promoter fragment was shown to be selectively expressed in megakaryoblastic leukemia cells but not in human epithelial cell lines. These findings identify a new tissue-specific ABCG2 promoter that is selectively expressed in pediatric M7 AML. We also show a relatively high incidence of ABCG2 mRNA expression in non-Down associated M7 AML, which may contribute to the relatively poor prognosis of the M7 AML subtype.

Epigenetic Control of MicroRNA Expression and Aging

MicroRNAs are a major category among the noncoding RNA fraction that negatively regulate gene expression at the post-transcriptional level, by either degrading the target messages or inhibiting their translation. MicroRNAs may be referred to as ‘dimmer switches’ of gene expression, because of their ability to repress gene expression without completely silencing it. Whether through up-regulating specific groups of microRNAs to suppress unwanted gene expressions, or by down-regulating other microRNAs whose target genes’ expression is necessary for cellular function, such as cell proliferation, apoptosis, or differentiation, these regulatory RNAs play pivotal roles in a wide variety of cellular processes. The equilibrium between these two groups of microRNA expressions largely determines the function of particular cell types. Our recent results with several model systems show that upon aging, there is a trend of up-regulation of microRNA expression, with concomitant inverse down-regulation of target genes. This review addresses molecular mechanisms that may provide the underlying control for this up-regulating trend, focusing on activation by various microRNAs’ own promoters, through binding with pivotal transcription factors, stress response, methylation of clustered DNA domains, etc. Thus, epigenomic control of aging may be due in part to heightened promoter activation of unwanted microRNA expressions, which in turn down-regulate their target gene products. Overriding and dampening the activation of these noncoding RNAs may prove to be a new frontier for future research, to delay aging and extend healthy life-span.

Genome-wide approaches for cancer gene discovery

One of the central aims of cancer research is to identify and characterize cancer-causing alterations in cancer genomes. In recent years, unprecedented advances in genome-wide sequencing, functional genomics technologies for RNA interference screens and methods for evaluating three-dimensional chromatin organization in vivo have resulted in important discoveries regarding human cancer. The cancer-causing genes identified from these new genome-wide technologies have also provided opportunities for effective and personalized cancer therapy. In this review, we describe some of the most recent technologies for cancer gene discovery. We also provide specific examples in which these technologies have proven remarkably successful in uncovering important cancer-causing alterations.

Computational identification of transcriptional regulators in human endotoxemia

One of the great challenges in the post-genomic era is to decipher the underlying principles governing the dynamics of biological responses. As modulating gene expression levels is among the key regulatory responses of an organism to changes in its environment, identifying biologically relevant transcriptional regulators and their putative regulatory interactions with target genes is an essential step towards studying the complex dynamics of transcriptional regulation. We present an analysis that integrates various computational and biological aspects to explore the transcriptional regulation of systemic inflammatory responses through a human endotoxemia model. Given a high-dimensional transcriptional profiling dataset from human blood leukocytes, an elementary set of temporal dynamic responses which capture the essence of a pro-inflammatory phase, a counter-regulatory response and a dysregulation in leukocyte bioenergetics has been extracted. Upon identification of these expression patterns, fourteen inflammation-specific gene batteries that represent groups of hypothetically ‘coregulated’ genes are proposed. Subsequently, statistically significant cis-regulatory modules (CRMs) are identified and decomposed into a list of critical transcription factors (34) that are validated largely on primary literature. Finally, our analysis further allows for the construction of a dynamic representation of the temporal transcriptional regulatory program across the host, deciphering possible combinatorial interactions among factors under which they might be active. Although much remains to be explored, this study has computationally identified key transcription factors and proposed a putative time-dependent transcriptional regulatory program associated with critical transcriptional inflammatory responses. These results provide a solid foundation for future investigations to elucidate the underlying transcriptional regulatory mechanisms under the host inflammatory response. Also, the assumption that coexpressed genes that are functionally relevant are more likely to share some common transcriptional regulatory mechanism seems to be promising, making the proposed framework become essential in unravelling context-specific transcriptional regulatory interactions underlying diverse mammalian biological processes.

Comparative analysis of acute and chronic corticosteroid pharmacogenomic effects in rat liver: transcriptional dynamics and regulatory structures

Background

Comprehensively understanding corticosteroid pharmacogenomic effects is an essential step towards an insight into the underlying molecular mechanisms for both beneficial and detrimental clinical effects. Nevertheless, even in a single tissue different methods of corticosteroid administration can induce different patterns of expression and regulatory control structures. Therefore, rich in vivo datasets of pharmacological time-series with two dosing regimens sampled from rat liver are examined for temporal patterns of changes in gene expression and their regulatory commonalities.

Results

The study addresses two issues, including (1) identifying significant transcriptional modules coupled with dynamic expression patterns and (2) predicting relevant common transcriptional controls to better understand the underlying mechanisms of corticosteroid adverse effects. Following the orientation of meta-analysis, an extended computational approach that explores the concept of agreement matrix from consensus clustering has been proposed with the aims of identifying gene clusters that share common expression patterns across multiple dosing regimens as well as handling challenges in the analysis of microarray data from heterogeneous sources, e.g. different platforms and time-grids in this study. Six significant transcriptional modules coupled with typical patterns of expression have been identified. Functional analysis reveals that virtually all enriched functions (gene ontologies, pathways) in these modules are shown to be related to metabolic processes, implying the importance of these modules in adverse effects under the administration of corticosteroids. Relevant putative transcriptional regulators (e.g. RXRF, FKHD, SP1F) are also predicted to provide another source of information towards better understanding the complexities of expression patterns and the underlying regulatory mechanisms of those modules.

Conclusions

We have proposed a framework to identify significant coexpressed clusters of genes across multiple conditions experimented from different microarray platforms, time-grids, and also tissues if applicable. Analysis on rich in vivo datasets of corticosteroid time-series yielded significant insights into the pharmacogenomic effects of corticosteroids, especially the relevance to metabolic side-effects. This has been illustrated through enriched metabolic functions in those transcriptional modules and the presence of GRE binding motifs in those enriched pathways, providing significant modules for further analysis on pharmacogenomic corticosteroid effects.

Genome-wide mapping of RNA Pol-II promoter usage in mouse tissues by ChIP-seq

Alternative promoters that are differentially used in various cellular contexts and tissue types add to the transcriptional complexity in mammalian genome. Identification of alternative promoters and the annotation of their activity in different tissues is one of the major challenges in understanding the transcriptional regulation of the mammalian genes and their isoforms. To determine the use of alternative promoters in different tissues, we performed ChIP-seq experiments using antibody against RNA Pol-II, in five adult mouse tissues (brain, liver, lung, spleen and kidney). Our analysis identified 38 639 Pol-II promoters, including 12 270 novel promoters, for both protein coding and non-coding mouse genes. Of these, 6384 promoters are tissue specific which are CpG poor and we find that only 34% of the novel promoters are located in CpG-rich regions, suggesting that novel promoters are mostly tissue specific. By identifying the Pol-II bound promoter(s) of each annotated gene in a given tissue, we found that 37% of the protein coding genes use alternative promoters in the five mouse tissues. The promoter annotations and ChIP-seq data presented here will aid ongoing efforts of characterizing gene regulatory regions in mammalian genomes.

Bioinformatic analysis of the human DHRS4 gene cluster and a proposed mechanism for its transcriptional regulation

BACKGROUND

The human DHRS4 gene cluster consists of three genes, DHRS4, DHRS4L2 and DHRS4L1. Among them, DHRS4 encodes NADP(H)-dependent retinol dehydrogenase/reductase. In a previous study, we investigated the alternative splicing of DHRS4 and DHRS4L2. DHRS4L1 was added to the gene cluster recently, but little is known about its structure and expression. To reveal the regulatory mechanism of the DHRS4 gene cluster expression, we studied the structure and transcription of DHRS4L1 in the context of the transcriptional behaviors of the human DHRS4 gene cluster. Based on the results of bioinformatics analysis, we propose a possible mechanism for the transcriptional regulation of the human DHRS4 gene cluster.

RESULTS

The homologous comparison analysis suggests that DHRS4, DHRS4L2 and DHRS4L1 are three homologous genes in human. DHRS4L1 and DHRS4L2 are paralogues of DHRS4, and DHRS4L2 is the most recent member of the DHRS4 gene cluster. In the minus strand of the human DHRS4 gene cluster, a gene transcribed in an antisense direction was found containing a 5' sequence overlapping the region of exon 1 and promoter of DHRS4. By cloning the full length of RNA variants through 5'RACE and 3'RACE, we identified two transcription start sites, within exon a2 and exon 1, of this newly named gene DHRS4L1 using neuroblastoma cell line BE(2)-M17. Analysis of exon composition in the transcripts of DHRS4 gene cluster revealed that exon 1 was absent in all the transcripts initiated from exon a1 of DHRS4L2 and exon a2 of DHRS4L1.

CONCLUSIONS

Alternatively spliced RNA variants are prevalent in the human DHRS4 gene cluster. Based on the analysis of gene transcripts and bioinformatic prediction, we propose here that antisense transcription may be involved in the transcriptional initiation regulation of DHRS4 and in the posttranscriptional splicing of DHRS4L2 and DRHS4L1 for the homologous identity of DHRS4 gene cluster. Beside the alternative transcriptional start sites, the antisense RNA is novel possible factor serving to remove exon 1 from the transcripts initiated from exon a1 and exon a2.

Estrogen regulation and physiopathologic significance of alternative promoters in breast cancer

Alternative promoters (AP) occur in >30% protein-coding genes and contribute to proteome diversity. However, large-scale analyses of AP regulation are lacking, and little is known about their potential physiopathologic significance. To better understand the transcriptomic effect of estrogens, which play a major role in breast cancer, we analyzed gene and AP regulation by estradiol in MCF7 cells using pan-genomic exon arrays. We thereby identified novel estrogen-regulated genes (ERG) and determined the regulation of AP-encoded transcripts in 150 regulated genes. In <30% cases, APs were regulated in a similar manner by estradiol, whereas in >70% cases, they were regulated differentially. The patterns of AP regulation correlated with the patterns of estrogen receptor alpha (ERalpha) and CCCTC-binding factor (CTCF) binding sites at regulated gene loci. Interestingly, among genes with differentially regulated (DR) APs, we identified cases where estradiol regulated APs in an opposite manner, sometimes without affecting global gene expression levels. This promoter switch was mediated by the DDX5/DDX17 family of ERalpha coregulators. Finally, genes with DR promoters were preferentially involved in specific processes (e.g., cell structure and motility, and cell cycle). We show, in particular, that isoforms encoded by the NET1 gene APs, which are inversely regulated by estradiol, play distinct roles in cell adhesion and cell cycle regulation and that their expression is differentially associated with prognosis in ER(+) breast cancer. Altogether, this study identifies the patterns of AP regulation in ERGs and shows the contribution of AP-encoded isoforms to the estradiol-regulated transcriptome as well as their physiopathologic significance in breast cancer.

Increasing alternative promoter repertories is positively associated with differential expression and disease susceptibility

Background

Alternative Promoter (AP) usages have been shown to enable diversified transcriptional regulation of individual gene in a context-specific (e.g., pathway, cell lineage, tissue type, and development stage et. ac.) way. Aberrant uses of APs have been directly linked to mechanism of certain human diseases. However, whether or not there exists a general link between a gene's AP repertoire and its expression diversity is currently unknown. The general relation between a gene's AP repertoire and its disease susceptibility also remains largely unexplored.

Methodology/Principal Findings

Based on the differential expression ratio inferred from all human microarray data in NCBI GEO and the list of disease genes curated in public repositories, we systemically analyzed the general relation of AP repertoire with expression diversity and disease susceptibility. We found that genes with APs are more likely to be differentially expressed and/or disease associated than those with Single Promoter (SP), and genes with more APs are more likely differentially expressed and disease susceptible than those with less APs. Further analysis showed that genes with increased number of APs tend to have increased length in all aspects of gene structure including 3′ UTR, be associated with increased duplicability, and have increased connectivity in protein-protein interaction network.

Conclusions

Our genome-wide analysis provided evidences that increasing alternative promoter repertories is positively associated with differential expression and disease susceptibility.

Annotation of gene promoters by integrative data-mining of ChIP-seq Pol-II enrichment data

Background

Use of alternative gene promoters that drive widespread cell-type, tissue-type or developmental gene regulation in mammalian genomes is a common phenomenon. Chromatin immunoprecipitation methods coupled with DNA microarray (ChIP-chip) or massive parallel sequencing (ChIP-seq) are enabling genome-wide identification of active promoters in different cellular conditions using antibodies against Pol-II. However, these methods produce enrichment not only near the gene promoters but also inside the genes and other genomic regions due to the non-specificity of the antibodies used in ChIP. Further, the use of these methods is limited by their high cost and strong dependence on cellular type and context.

Methods

We trained and tested different state-of-art ensemble and meta classification methods for identification of Pol-II enriched promoter and Pol-II enriched non-promoter sequences, each of length 500 bp. The classification models were trained and tested on a bench-mark dataset, using a set of 39 different feature variables that are based on chromatin modification signatures and various DNA sequence features. The best performing model was applied on seven published ChIP-seq Pol-II datasets to provide genome wide annotation of mouse gene promoters.

Results

We present a novel algorithm based on supervised learning methods to discriminate promoter associated Pol-II enrichment from enrichment elsewhere in the genome in ChIP-chip/seq profiles. We accumulated a dataset of 11,773 promoter and 46,167 non-promoter sequences, each of length 500 bp, generated from RNA Pol-II ChIP-seq data of five tissues (Brain, Kidney, Liver, Lung and Spleen). We evaluated the classification models in building the best predictor and found that Bagging and Random Forest based approaches give the best accuracy. We implemented the algorithm on seven different published ChIP-seq datasets to provide a comprehensive set of promoter annotations for both protein-coding and non-coding genes in the mouse genome. The resulting annotations contain 13,413 (4,747) protein-coding (non-coding) genes with single promoters and 9,929 (1,858) protein-coding (non-coding) genes with two or more alternative promoters, and a significant number of unassigned novel promoters.

Conclusion

Our new algorithm can successfully predict the promoters from the genome wide profile of Pol-II bound regions. In addition, our algorithm performs significantly better than existing promoter prediction methods and can be applied for genome-wide predictions of Pol-II promoters.

NF-YC complexity is generated by dual promoters and alternative splicing

The CCAAT box is a DNA element present in the majority of human promoters, bound by the trimeric NF-Y, composed of NF-YA, NF-YB, and NF-YC subunits. We describe and characterize novel isoforms of one of the two histone-like subunits, NF-YC. The locus generates a minimum of four splicing products, mainly located within the Q-rich activation domain. The abundance of each isoform is cell-dependent; only one major NF-YC isoform is present in a given cell type. The 37- and 50-kDa isoforms are mutually exclusive, and preferential pairings with NF-YA isoforms possess different transcriptional activities, with specific combinations being more active on selected promoters. The transcriptional regulation of the NF-YC locus is also complex, and mRNAs arise from the two promoters P1 and P2. Transient transfections, chromatin immunoprecipitations, and reverse transcription-PCRs indicate that P1 has a robust housekeeping activity; P2 possesses a lower basal activity, but it is induced in response to DNA damage in a p53-dependent way. Alternative promoter usage directly affects NF-YC splicing, with the 50-kDa transcript being excluded from P2. Specific functional inactivation of the 37-kDa isoform affects the basal levels of G(1)/S blocking and pro-apoptotic genes but not G(2)/M promoters. In summary, our data highlight an unexpected degree of complexity and regulation of the NF-YC gene, demonstrating the existence of a discrete cohort of NF-Y trimer subtypes resulting from the functional diversification of Q-rich transactivating subunits and a specific role of the 37-kDa isoform in suppression of the DNA damage-response under growing conditions.

Comparative study on ChIP-seq data: normalization and binding pattern characterization

MOTIVATION

Antibody-based Chromatin Immunoprecipitation assay followed by high-throughput sequencing technology (ChIP-seq) is a relatively new method to study the binding patterns of specific protein molecules over the entire genome. ChIP-seq technology allows scientist to get more comprehensive results in shorter time. Here, we present a non-linear normalization algorithm and a mixture modeling method for comparing ChIP-seq data from multiple samples and characterizing genes based on their RNA polymerase II (Pol II) binding patterns.

RESULTS

We apply a two-step non-linear normalization method based on locally weighted regression (LOESS) approach to compare ChIP-seq data across multiple samples and model the difference using an Exponential-Normal(K) mixture model. Fitted model is used to identify genes associated with differential binding sites based on local false discovery rate (fdr). These genes are then standardized and hierarchically clustered to characterize their Pol II binding patterns. As a case study, we apply the analysis procedure comparing normal breast cancer (MCF7) to tamoxifen-resistant (OHT) cell line. We find enriched regions that are associated with cancer (P < 0.0001). Our findings also imply that there may be a dysregulation of cell cycle and gene expression control pathways in the tamoxifen-resistant cells. These results show that the non-linear normalization method can be used to analyze ChIP-seq data across multiple samples.

AVAILABILITY

Data are available at http://www.bmi.osu.edu/~khuang/Data/ChIP/RNAPII/.

Novel genes in cell cycle control and lipid metabolism with dynamically regulated binding sites for sterol regulatory element-binding protein 1 and RNA polymerase II in HepG2 cells detected by chromatin immunoprecipitation with microarray detection

Sterol regulatory element-binding proteins 1 and 2 (SREBP-1 and SREBP-2) are important regulators of genes involved in cholesterol and fatty acid metabolism, but have also been implicated in the regulation of the cell cycle and have been associated with the pathogenesis of type 2 diabetes, atherosclerosis and obesity, among others. In this study, we aimed to characterize the binding sites of SREBP-1 and RNA polymerase II through chromatin immunoprecipitation and microarray analysis in 1% of the human genome, as defined by the Encyclopaedia of DNA Elements consortium, in a hepatocellular carcinoma cell line (HepG2). Our data identified novel binding sites for SREBP-1 in genes directly or indirectly involved in cholesterol metabolism, e.g. apolipoprotein C-III (APOC3). The most interesting biological findings were the binding sites for SREBP-1 in genes for host cell factor C1 (HCFC1), involved in cell cycle regulation, and for filamin A (FLNA). For RNA polymerase II, we found binding sites at classical promoters, but also in intergenic and intragenic regions. Furthermore, we found evidence of sterol-regulated binding of SREBP-1 and RNA polymerase II to HCFC1 and FLNA. From the results of this work, we infer that SREBP-1 may be involved in processes other than lipid metabolism.