Search

Mounting evidence supports the idea that transcriptional patterns serve as more specific identifiers of active enhancers than histone marks; however, the optimal strategy to identify active enhancers both experimentally and computationally has not been determined. Here, we compared 13 genome-wide RNA sequencing (RNA-seq) assays in K562 cells and show that nuclear run-on followed by cap-selection assay (GRO/PRO-cap) has advantages in enhancer RNA detection and active enhancer identification. We also introduce a tool, peak identifier for nascent transcript starts (PINTS), to identify active promoters and enhancers genome wide and pinpoint the precise location of 5' transcription start sites. Finally, we compiled a comprehensive enhancer candidate compendium based on the detected enhancer RNA (eRNA) transcription start sites (TSSs) available in 120 cell and tissue types, which can be accessed at https://pints.yulab.org . With knowledge of the best available assays and pipelines, this large-scale annotation of candidate enhancers will pave the way for selection and characterization of their functions in a time- and labor-efficient manner.

Background Promoters, non-coding DNA sequences located at upstream regions of the transcription start site of genes/gene clusters, are essential regulatory elements for the initiation and regulation of transcriptional processes. Furthermore, identifying promoters in DNA sequences and genomes significantly contributes to discovering entire structures of genes of interest. Therefore, exploration of promoter regions is one of the most imperative topics in molecular genetics and biology. Besides experimental techniques, computational methods have been developed to predict promoters. In this study, we propose iPromoter-Seqvec – an efficient computational model to predict TATA and non-TATA promoters in human and mouse genomes using bidirectional long short-term memory neural networks in combination with sequence-embedded features extracted from input sequences. The promoter and non-promoter sequences were retrieved from the Eukaryotic Promoter database and then were refined to create four benchmark datasets. Results The area under the receiver operating characteristic curve (AUCROC) and the area under the precision-recall curve (AUCPR) were used as two key metrics to evaluate model performance. Results on independent test sets showed that iPromoter-Seqvec outperformed other state-of-the-art methods with AUCROC values ranging from 0.85 to 0.99 and AUCPR values ranging from 0.86 to 0.99. Models predicting TATA promoters in both species had slightly higher predictive power compared to those predicting non-TATA promoters. With a novel idea of constructing artificial non-promoter sequences based on promoter sequences, our models were able to learn highly specific characteristics discriminating promoters from non-promoters to improve predictive efficiency. Conclusions iPromoter-Seqvec is a stable and robust model for predicting both TATA and non-TATA promoters in human and mouse genomes. Our proposed method was also deployed as an online web server with a user-friendly interface to support research communities. Links to our source codes and web server are available at https://github.com/mldlproject/2022-iPromoter-Seqvec.

Background More than half of human protein-coding genes have an alternative transcription start site (TSS). We aimed to investigate the contribution of alternative TSSs to the acute-stress-induced transcriptome response in human tissue (skeletal muscle) using the cap analysis of gene expression approach. TSSs were examined at baseline and during recovery after acute stress (a cycling exercise). Results We identified 44,680 CAGE TSS clusters (including 3764 first defined) belonging to 12,268 genes and annotated for the first time 290 TSSs belonging to 163 genes. The transcriptome dynamically changes during the first hours after acute stress; the change in the expression of 10% of genes was associated with the activation of alternative TSSs, indicating differential TSSs usage. The majority of the alternative TSSs do not increase proteome complexity suggesting that the function of thousands of alternative TSSs is associated with the fine regulation of mRNA isoform expression from a gene due to the transcription factor-specific activation of various alternative TSSs. We identified individual muscle promoter regions for each TSS using muscle open chromatin data (ATAC-seq and DNase-seq). Then, using the positional weight matrix approach we predicted time course activation of “classic” transcription factors involved in response of skeletal muscle to contractile activity, as well as diversity of less/un-investigated factors. Conclusions Transcriptome response induced by acute stress related to activation of the alternative TSSs indicates that differential TSSs usage is an essential mechanism of fine regulation of gene response to stress stimulus. A comprehensive resource of accurate TSSs and individual promoter regions for each TSS in muscle was created. This resource together with the positional weight matrix approach can be used to accurate prediction of TFs in any gene(s) of interest involved in the response to various stimuli, interventions or pathological conditions in human skeletal muscle.

Background New genes continuously emerge from non-coding DNA or by diverging from existing genes, but most of them are rapidly lost and only a few become fixed within the population. We hypothesized that young genes are subject to transcriptional and post-transcriptional regulation to limit their expression and minimize their exposure to purifying selection. Results We performed a protein-based homology search across the tree of life to determine the evolutionary age of protein-coding genes present in the rice genome. We found that young genes in rice have relatively low expression levels, which can be attributed to distal enhancers, and closed chromatin conformation at their transcription start sites (TSS). The chromatin in TSS regions can be re-modeled in response to abiotic stress, indicating conditional expression of young genes. Furthermore, transcripts of young genes in Arabidopsis tend to be targeted by nonsense-mediated RNA decay, presenting another layer of regulation limiting their expression. Conclusions These data suggest that transcriptional and post-transcriptional mechanisms contribute to the conditional expression of young genes, which may alleviate purging selection while providing an opportunity for phenotypic exposure and functionalization.

Transcription initiates at promoters, DNA regions recognized by a DNA-dependent RNA polymerase. We previously identified horizontally acquired Escherichia coli promoters from which the direction of transcription was unclear. In the present study, we show that more than half of these promoters are bidirectional and drive divergent transcription. Using genome-scale approaches, we demonstrate that 19% of all transcription start sites detected in E. coli are associated with a bidirectional promoter. Bidirectional promoters are similarly common in diverse bacteria and archaea, and have inherent symmetry: specific bases required for transcription initiation are reciprocally co-located on opposite DNA strands. Bidirectional promoters enable co-regulation of divergent genes and are enriched in both intergenic and horizontally acquired regions. Divergent transcription is conserved among bacteria, archaea and eukaryotes, but the underlying mechanisms for bidirectionality are different.

Background The zone adjacent to a transcription start site (TSS), namely, the promoter, is primarily involved in the process of DNA transcription initiation and regulation. As a result, proper promoter identification is critical for further understanding the mechanism of the networks controlling genomic regulation. A number of methodologies for the identification of promoters have been proposed. Nonetheless, due to the great heterogeneity existing in promoters, the results of these procedures are still unsatisfactory. In order to establish additional discriminative characteristics and properly recognize promoters, we developed the hybrid model for promoter identification (HMPI), a hybrid deep learning model that can characterize both the native sequences of promoters and the morphological outline of promoters at the same time. We developed the HMPI to combine a method called the PSFN (promoter sequence features network), which characterizes native promoter sequences and deduces sequence features, with a technique referred to as the DSPN (deep structural profiles network), which is specially structured to model the promoters in terms of their structural profile and to deduce their structural attributes. Results The HMPI was applied to human, plant and Escherichia coli K-12 strain datasets, and the findings showed that the HMPI was successful at extracting the features of the promoter while greatly enhancing the promoter identification performance. In addition, after the improvements of synthetic sampling, transfer learning and label smoothing regularization, the improved HMPI models achieved good results in identifying subtypes of promoters on prokaryotic promoter datasets. Conclusions The results showed that the HMPI was successful at extracting the features of promoters while greatly enhancing the performance of identifying promoters on both eukaryotic and prokaryotic datasets, and the improved HMPI models are good at identifying subtypes of promoters on prokaryotic promoter datasets. The HMPI is additionally adaptable to different biological functional sequences, allowing for the addition of new features or models.

Genes are transcribed from transcription start sites (TSSs), and their position in a genome is strictly controlled to avoid mis-expression of undesired regions. In this study, we designed and developed a methodology for the evaluation of promoter context, which detects proximal promoter regions from − 200 to − 60 bp relative to a TSS, in Arabidopsis and rice genomes. The method positively evaluates spacer sequences and Regulatory Element Groups, but not core promoter elements like TATA boxes, and is able to predict the position of a TSS within a width of 200 bp. An important feature of the evaluation/prediction method is its independence of the core promoter elements, which was demonstrated by successful prediction of all the TATA, GA, and coreless types of promoters without notable differences in the accuracy of prediction. The positive relationship identified between the evaluation scores and gene expression levels suggests that this method is useful for the evaluation of promoter maturity.

miRNAs form a class of noncoding RNAs, involved in post-transcriptional regulation of gene expression, broadly studied for their involvement in physiological and pathological context. Inhibition of mature miRNA transcripts, commonly used in miRNA loss-of-function experiments, may not be specific in case of miRNAs with high sequence homology, e.g. miRNAs from the same seed family. Phenotypic effects of miRNA repression might be biased by the repression of highly similar miRNAs. Another challenge is simultaneous inhibition of multiple miRNAs encoded within policistronic clusters, potentially co-regulating common biological processes. To elucidate roles of miRNA clusters and miRNAs with high sequence homology, it is of key importance to selectively repress only the miRNAs of interest. Targeting miRNAs on genomic level with CRISPR/dCas9-based methods is an attractive alternative to blocking mature miRNAs. Yet, so far no clear guidelines on the design of CRISPR inhibition (CRISPRi) experiments, specifically for miRNA repression, have been proposed. To address this need, here we propose a strategy for effective inhibition of miRNAs and miRNA clusters using CRISPRi. We provide clues on how to approach the challenges in using CRISPR/dCas in miRNA studies, which include prediction of miRNA transcription start sites (TSSs) and the design of single guide RNAs (sgRNAs). The strategy implements three TSS prediction online tools, dedicated specifically for miRNAs: miRStart, FANTOM 5 miRNA atlas, DIANA-miRGen, and CRISPOR tool for sgRNAs design; it includes testing and selection of optimal sgRNAs. We demonstrate that compared to siRNA/shRNA-based miRNA silencing, CRISPRi improves the repression specificity for miRNAs with highly similar sequence and contribute to higher uniformity of the effects of silencing the whole miRNA clusters. This strategy may be adapted for CRISPR-mediated activation (CRISPRa) of miRNA expression.

Distal enhancers play pivotal roles in development and disease yet remain one of the least understood regulatory elements. We used massively parallel reporter assays to perform functional comparisons of two leading enhancer models and find that gene-distal transcription start sites are robust predictors of active enhancers with higher resolution than histone modifications. We show that active enhancer units are precisely delineated by active transcription start sites, validate that these boundaries are sufficient for capturing enhancer function, and confirm that core promoter sequences are necessary for this activity. We assay adjacent enhancers and find that their joint activity is often driven by the stronger unit within the cluster. Finally, we validate these results through functional dissection of a distal enhancer cluster using CRISPR–Cas9 deletions. In summary, definition of high-resolution enhancer boundaries enables deconvolution of complex regulatory loci into modular units. A massively parallel reporter assay (eSTARR-seq) shows that gene-distal transcription start sites can delineate active enhancers with higher resolution than histone modifications.

Animal chromosomes are partitioned into contact domains. Pathogenic domain disruptions can result from chromosomal rearrangements or perturbation of architectural factors. However, such broad-scale alterations are insufficient to define the minimal requirements for domain formation. Moreover, to what extent domains can be engineered is just beginning to be explored. In an attempt to create contact domains, we inserted a 2-kb DNA sequence underlying a tissue-invariant domain boundary-containing a CTCF-binding site (CBS) and a transcription start site (TSS)-into 16 ectopic loci across 11 chromosomes, and characterized its architectural impact. Depending on local constraints, this fragment variably formed new domains, partitioned existing ones, altered compartmentalization and initiated contacts reflecting chromatin loop extrusion. Deletions of the CBS or the TSS individually or in combination within inserts revealed its distinct contributions to genome folding. Altogether, short DNA insertions can suffice to shape the spatial genome in a manner influenced by chromatin context.

While CpG dinucleotides are significantly reduced compared to other dinucleotides in mammalian genomes, they can congregate and form CpG islands, which localize around the 5ʹ regions of genes, where they function as promoters. CpG-island promoters are generally unmethylated and are often found in housekeeping genes. However, their nucleotide sequences and existence per se are not conserved between humans and mice, which may be due to evolutionary gain and loss of the regulatory regions. In this study, human and rhesus monkey genomes, with moderately conserved sequences, were compared at base resolution. Using transcription start site data, we first validated our methods’ ability to identify orthologous promoters and indicated a limitation using the 5ʹ end of curated gene models, such as NCBI RefSeq, as their transcription start sites. We found that, in addition to deamination mutations, insertions and deletions of bases, repeats, and long fragments contributed to the mutations of CpG dinucleotides. We also observed that the G + C contents tended to change in CpG-poor environments, while CpG content was altered in G + C-rich environments. While loss of CpG islands can be caused by gradual decreases in CpG sites, gain of these islands appear to require two distinct nucleotide altering steps. Taken together, our findings provide novel insights into the process of acquisition and diversification of CpG-island promoters in vertebrates. Electronic supplementary material The online version of this article (10.1007/s00335-020-09844-2) contains supplementary material, which is available to authorized users.

While horizontal gene transfer is prevalent across the biosphere, the regulatory features that enable expression and functionalization of foreign DNA remain poorly understood. Here, we combine high-throughput promoter activity measurements and large-scale genomic analysis of regulatory regions to investigate the cross-compatibility of regulatory elements (REs) in bacteria. Functional characterization of thousands of natural REs in three distinct bacterial species revealed distinct expression patterns according to RE and recipient phylogeny. Host capacity to activate foreign promoters was proportional to their genomic GC content, while many low GC regulatory elements were both broadly active and had more transcription start sites across hosts. The difference in expression capabilities could be explained by the influence of the host GC content on the stringency of the AT-rich canonical σ70 motif necessary for transcription initiation. We further confirm the generalizability of this model and find widespread GC content adaptation of the σ70 motif in a set of 1,545 genomes from all major bacterial phyla. Our analysis identifies a key mechanism by which the strength of the AT-rich σ70 motif relative to a host’s genomic GC content governs the capacity for expression of acquired DNA. These findings shed light on regulatory adaptation in the context of evolving genomic composition.

Gene expression using alternative transcription start sites (TSSs) is an important transcriptional regulatory mechanism for environmental responses in eukaryotes. Here, we identify two alternative TSSs in the enolase-encoding gene (enoA) in Aspergillus oryzae, an industrially important filamentous fungus. TSS use in enoA is strictly dependent on the difference in glycolytic and gluconeogenic carbon sources. Transcription from the upstream TSS (uTSS) or downstream TSS (dTSS) predominantly occurs under gluconeogenic or glycolytic conditions, respectively. In addition to enoA, most glycolytic genes involved in reversible reactions possess alternative TSSs. The fbaA gene, which encodes fructose-bisphosphate aldolase, also shows stringent alternative TSS selection, similar to enoA. Alignment of promoter sequences of enolase-encoding genes in Aspergillus predicted two conserved regions that contain a putative cis-element required for enoA transcription from each TSS. However, uTSS-mediated transcription of the acuN gene, an enoA ortholog in Aspergillus nidulans, is not strictly dependent on carbon source, unlike enoA. Furthermore, enoA transcript levels in glycolytic conditions are higher than in gluconeogenic conditions. Conversely, acuN is more highly transcribed in gluconeogenic conditions. This suggests that the stringent usage of alternative TSSs and higher transcription in glycolytic conditions in enoA may reflect that the A. oryzae evolutionary genetic background was domesticated by exclusive growth in starch-rich environments. These findings provide novel insights into the complexity and diversity of transcriptional regulation of glycolytic/gluconeogenic genes among Aspergilli.

The chromatin architecture in promoters is thought to regulate gene expression, but it remains uncertain how most transcription factors (TFs) impact nucleosome position. The MuvB TF complex regulates cell-cycle dependent gene-expression and is critical for differentiation and proliferation during development and cancer. MuvB can both positively and negatively regulate expression, but the structure of MuvB and its biochemical function are poorly understood. Here we determine the overall architecture of MuvB assembly and the crystal structure of a subcomplex critical for MuvB function in gene repression. We find that the MuvB subunits LIN9 and LIN37 function as scaffolding proteins that arrange the other subunits LIN52, LIN54 and RBAP48 for TF, DNA, and histone binding, respectively. Biochemical and structural data demonstrate that MuvB binds nucleosomes through an interface that is distinct from LIN54-DNA consensus site recognition and that MuvB increases nucleosome occupancy in a reconstituted promoter. We find in arrested cells that MuvB primarily associates with a tightly positioned +1 nucleosome near the transcription start site (TSS) of MuvB-regulated genes. These results support a model that MuvB binds and stabilizes nucleosomes just downstream of the TSS on its target promoters to repress gene expression., The MuvB protein complex regulates genes that are differentially expressed through the cell cycle, yet its precise molecular function has remained unclear. Here the authors reveal MuvB associates with the nucleosome adjacent to the transcription start site of cell-cycle genes and that the tight positioning of this nucleosome correlates with MuvB-dependent gene repression.

CUX1, encoding a homeodomain-containing transcription factor, is recurrently deleted or mutated in multiple tumor types. In myeloid neoplasms, CUX1 deletion or mutation carries a poor prognosis. We have previously established that CUX1 functions as a tumor suppressor in hematopoietic cells across multiple organisms. Others, however, have described oncogenic functions of CUX1 in solid tumors, often attributed to truncated CUX1 isoforms, p75 and p110, generated by an alternative transcriptional start site or post-translational cleavage, respectively. Given the clinical relevance, it is imperative to clarify these discrepant activities. Herein, we sought to determine the CUX1 isoforms expressed in hematopoietic cells and find that they express the full-length p200 isoform. Through the course of this analysis, we found no evidence of the p75 alternative transcript in any cell type examined. Using an array of orthogonal approaches, including biochemistry, proteomics, CRISPR/Cas9 genomic editing, and analysis of functional genomics datasets across a spectrum of normal and malignant tissue types, we found no data to support the existence of the CUX1 p75 isoform as previously described. Based on these results, prior studies of p75 require reevaluation, including the interpretation of oncogenic roles attributed to CUX1.

[Background] Variation in gene expression underlies interindividual variability in relevant traits including immune response. However, the genetic variation responsible for these gene expression changes remains largely unknown. Among the non-coding variants that could be relevant, transposable element insertions are promising candidates as they have been shown to be a rich and diverse source of cis-regulatory elements., [Results] In this work, we use a population genetics approach to identify transposable element insertions likely to increase the tolerance of Drosophila melanogaster to bacterial infection by affecting the expression of immune-related genes. We identify 12 insertions associated with allele-specific expression changes in immune-related genes. We experimentally validate three of these insertions including one likely to be acting as a silencer, one as an enhancer, and one with a dual role as enhancer and promoter. The direction in the change of gene expression associated with the presence of several of these insertions is consistent with an increased survival to infection. Indeed, for one of the insertions, we show that this is the case by analyzing both natural populations and CRISPR/Cas9 mutants in which the insertion is deleted from its native genomic context., [Conclusions] We show that transposable elements contribute to gene expression variation in response to infection in D. melanogaster and that this variation is likely to affect their survival capacity. Because the role of transposable elements as regulatory elements is not restricted to Drosophila, transposable elements are likely to play a role in immune response in other organisms as well., This project has received funding from the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation program (H2020-ERC-2014-CoG-647900) and under the FP7 program (FP7-PEOPLE-2011-CIG-293860) and by the MEC/FEDER (BFU2014-57779-P). A.U. was a FPI fellow (BES-2012-052999) and JG was a Ramon y Cajal fellow (RYC-2010-07306).

To further the understanding of the evolution of transcriptional regulation, we profiled genome-wide transcriptional start sites (TSSs) in two sub-species, Bos taurus taurus and Bos taurus indicus, that diverged approximately 500,000 years ago. Evolutionary and developmental-stage differences in TSSs were detected across the sub-species, including translocation of dominant TSS and changes in TSS distribution. The 16% of all SNPs located in significant differentially used TSS clusters across sub-species had significant shifts in allele frequency (472 SNPs), indicating they may have been subject to selection. In spleen and muscle, a higher relative TSS expression was observed in Bos indicus than Bos taurus for all heat shock protein genes, which may be responsible for the tropical adaptation of Bos indicus., Forutan et al. measure RNA expression and map transcription start sites (TSSs) between two sub-species of cattle in adult and fetal tissues and demonstrate translocation and changes in TSS distribution between the sub-species. This study provides insight into cattle adaptation and provides support of rapidly evolving TSS between two recently diverged subspecies.

Transcriptional cofactors (COFs) communicate regulatory cues from enhancers to promoters and are central effectors of transcription activation and gene expression1. Although some COFs have been shown to prefer certain promoter types2–5 over others (for example, see refs 6,7), the extent to which different COFs display intrinsic specificities for distinct promoters is unclear. Here we use a high-throughput promoter-activity assay in Drosophila melanogaster S2 cells to screen 23 COFs for their ability to activate 72,000 candidate core promoters (CPs). We observe differential activation of CPs, indicating distinct regulatory preferences or ‘compatibilities’8,9 between COFs and specific types of CPs. These functionally distinct CP types are differentially enriched for known sequence elements2,4, such as the TATA box, downstream promoter element (DPE) or TCT motif, and display distinct chromatin properties at endogenous loci. Notably, the CP types differ in their relative abundance of H3K4me3 and H3K4me1 marks (see also refs 10–12), suggesting that these histone modifications might distinguish trans-regulatory factors rather than promoter- versus enhancer-type cis-regulatory elements. We confirm the existence of distinct COF–CP compatibilities in two additional Drosophila cell lines and in human cells, for which we find COFs that prefer TATA-box or CpG-island promoters, respectively. Distinct compatibilities between COFs and promoters can explain how different enhancers specifically activate distinct sets of genes9, alternative promoters within the same genes, and distinct transcription start sites within the same promoter13. Thus, COF–promoter compatibilities may underlie distinct transcriptional programs in species as divergent as flies and humans. A screen of 23 transcriptional cofactors for their ability to activate 72,000 candidate core promoters in Drosophila melanogaster identified distinct compatibility groups, providing insight into mechanisms that underlie the selective activation of transcriptional programs.

Although DNA replication is a fundamental aspect of biology, it is not known what determines where DNA replication starts and stops in the human genome. Here we directly identify and quantitatively compare sites of replication initiation and termination in untransformed human cells. We report that replication preferentially initiates at the transcription start site of genes occupied by high levels of RNA polymerase II, and terminates at their polyadenylation sites, thus ensuring global co-directionality of transcription and replication, particularly at gene 5’ ends. During replication stress, replication initiation is stimulated downstream of genes and termination is redistributed to gene bodies; this globally re-orients replication relative to transcription around gene 3’ ends. These data suggest that replication initiation and termination are coupled to transcription in human cells, and propose a model for the impact of replication stress on genome integrity.

Progesterone and its receptor, PR, are essential for uterine leiomyoma (LM, a.k.a., fibroid) tumorigenesis, but the underlying cellular and molecular mechanisms remain unclear. The receptor activator of NF-κB (RANKL) was recently identified as a novel progesterone/PR-responsive gene that plays an important role in promoting LM growth. Here, we used RANKL as a representative gene to investigate how steroid hormone, genetic, and epigenetic signals are integrated to regulate LM stem cell (LSC) function. We demonstrated that RANKL specifically upregulates LSC proliferation through activation of Cyclin D1. RANKL gene transcription was robustly induced by the progesterone agonist R5020, leading to a dramatically higher RANKL expression in LM compared to adjacent myometrial (MM) tissue. MethylCap-Seq revealed a differentially methylated region (DMR) adjacent to the distal PR-binding site (PRBS) 87 kb upstream of the RANKL transcription start site. Hypermethylation of the DMR inhibited recruitment of PR to the adjacent PRBS. Luciferase assays indicated that the DMR and distal PRBS constitute a novel RANKL distal regulatory element that actively regulates RANKL expression. Furthermore, MED12 physically interacts with PR in LM tissue. The interaction between MED12 and PR, binding of PR and MED12 to PRBS, and RANKL gene expression are significantly higher in LM containing a distinct MED12 mutation (G44D) than in LM with wild-type MED12. In summary, our findings suggest that DNA methylation and MED12 mutation together constitute a complex regulatory network that affects progesterone/PR-mediated RANKL gene expression, with an important role in activating stem cell proliferation and fibroid tumor development.

Background Cancer genomes are shaped by mutational processes with complex spatial variation at multiple scales. Entire classes of regulatory elements are affected by local variations in mutation frequency. However, the underlying mechanisms with functional and genetic determinants remain poorly understood. Results We characterise the mutational landscape of 1.3 million gene-regulatory and chromatin architectural elements in 2419 whole cancer genomes with transcriptional and pathway activity, functional conservation and recurrent driver events. We develop RM2, a statistical model that quantifies mutational enrichment or depletion in classes of genomic elements through genetic, trinucleotide and megabase-scale effects. We report a map of localised mutational processes affecting CTCF binding sites, transcription start sites (TSS) and tissue-specific open-chromatin regions. Increased mutation frequency in TSSs associates with mRNA abundance in most cancer types, while open-chromatin regions are generally enriched in mutations. We identify ~ 10,000 CTCF binding sites with core DNA motifs and constitutive binding in 66 cell types that represent focal points of mutagenesis. We detect site-specific mutational signature enrichments, such as SBS40 in open-chromatin regions in prostate cancer and SBS17b in CTCF binding sites in gastrointestinal cancers. Candidate drivers of localised mutagenesis are also apparent: BRAF mutations associate with mutational enrichments at CTCF binding sites in melanoma, and ARID1A mutations with TSS-specific mutagenesis in pancreatic cancer. Conclusions Our method and catalogue of localised mutational processes provide novel perspectives to cancer genome evolution, mutagenesis, DNA repair and driver gene discovery. The functional and genetic correlates of mutational processes suggest mechanistic hypotheses for future studies.

MicroRNAs (miRNAs) play important roles in post-transcriptional gene regulation and phenotype development. Understanding the regulation of miRNA genes is critical to understand gene regulation. One of the challenges to study miRNA gene regulation is the lack of condition-specific annotation of miRNA transcription start sites (TSSs). Unlike protein-coding genes, miRNA TSSs can be tens of thousands of nucleotides away from the precursor miRNAs and they are hard to be detected by conventional RNA-Seq experiments. A number of studies have been attempted to computationally predict miRNA TSSs. However, high-resolution condition-specific miRNA TSS prediction remains a challenging problem. Recently, deep learning models have been successfully applied to various bioinformatics problems but have not been effectively created for condition-specific miRNA TSS prediction. Here we created a two-stream deep learning model called D-miRT for computational prediction of condition-specific miRNA TSSs (http://hulab.ucf.edu/research/projects/DmiRT/). D-miRT is a natural fit for the integration of low-resolution epigenetic features (DNase-Seq and histone modification data) and high-resolution sequence features. Compared with alternative computational models on different sets of training data, D-miRT outperformed all baseline models and demonstrated high accuracy for condition-specific miRNA TSS prediction tasks. Comparing with the most recent approaches on cell-specific miRNA TSS identification using cell lines that were unseen to the model training processes, D-miRT also showed superior performance.

Ribosomal RNA (rRNA) is most highly expressed in rapidly growing bacteria and is drastically downregulated under stress conditions by the global transcriptional regulator DksA and the alarmone ppGpp. Here, we determined cryo-electron microscopy structures of the Escherichia coli RNA polymerase (RNAP) σ70 holoenzyme during rRNA promoter recognition with and without DksA/ppGpp. RNAP contacts the UP element using dimerized α subunit carboxyl-terminal domains and scrunches the template DNA with the σ finger and β’ lid to select the transcription start site favorable for rapid promoter escape. Promoter binding induces conformational change of σ domain 2 that opens a gate for DNA loading and ejects σ1.1 from the RNAP cleft to facilitate open complex formation. DksA/ppGpp binding also opens the DNA loading gate, which is not coupled to σ1.1 ejection and impedes open complex formation. These results provide a molecular basis for the exceptionally active rRNA transcription and its vulnerability to DksA/ppGpp., Ribosomal RNA (rRNA) expression is regulated at the initiation stage of RNA synthesis. Here, the authors report cryo-EM structures of E. coli RNA polymerase and rRNA promoter complex with DksA/ppGpp on the way to open complex formation, identifying key steps in promoter recognition and opening.

The Origin Recognition Complex (ORC) is an evolutionarily conserved six-subunit protein complex that binds specific sites at many locations to coordinately replicate the entire eukaryote genome. Though highly conserved in structure, ORC’s selectivity for replication origins has diverged tremendously between yeasts and humans to adapt to vastly different life cycles. In this work, we demonstrate that the selectivity determinant of ORC for DNA binding lies in a 19-amino acid insertion helix in the Orc4 subunit, which is present in yeast but absent in human. Removal of this motif from Orc4 transforms the yeast ORC, which selects origins based on base-specific binding at defined locations, into one whose selectivity is dictated by chromatin landscape and afforded with plasticity, as reported for human. Notably, the altered yeast ORC has acquired an affinity for regions near transcriptional start sites (TSSs), which the human ORC also favors., In most model yeast species the Origin Recognition Complex (ORC) binds defined and species-specific base sequences while in humans what determines the binding appears to be more complex. Here the authors reveal that the yeast’s ORC complex binding specificity is dependent on a 19-amino acid insertion helix in the Orc4 subunit which is lost in human.

Long-read sequencing (LRS) has become a standard approach for transcriptome analysis in recent years. Bovine alphaherpesvirus 1 (BoHV-1) is an important pathogen of cattle worldwide. This study reports the profiling of the dynamic lytic transcriptome of BoHV-1 using two long-read sequencing (LRS) techniques, the Oxford Nanopore Technologies MinION, and the LoopSeq synthetic LRS methods, using multiple library preparation protocols. In this work, we annotated viral mRNAs and non-coding transcripts, and a large number of transcript isoforms, including transcription start and end sites, as well as splice variants of BoHV-1. Our analysis demonstrated an extremely complex pattern of transcriptional overlaps.

Background The majority of eukaryotic promoters utilize multiple transcription start sites (TSSs). How multiple TSSs are specified at individual promoters across eukaryotes is not understood for most species. In Saccharomyces cerevisiae, a pre-initiation complex (PIC) comprised of Pol II and conserved general transcription factors (GTFs) assembles and opens DNA upstream of TSSs. Evidence from model promoters indicates that the PIC scans from upstream to downstream to identify TSSs. Prior results suggest that TSS distributions at promoters where scanning occurs shift in a polar fashion upon alteration in Pol II catalytic activity or GTF function. Results To determine the extent of promoter scanning across promoter classes in S. cerevisiae, we perturb Pol II catalytic activity and GTF function and analyze their effects on TSS usage genome-wide. We find that alterations to Pol II, TFIIB, or TFIIF function widely alter the initiation landscape consistent with promoter scanning operating at all yeast promoters, regardless of promoter class. Promoter architecture, however, can determine the extent of promoter sensitivity to altered Pol II activity in ways that are predicted by a scanning model. Conclusions Our observations coupled with previous data validate key predictions of the scanning model for Pol II initiation in yeast, which we term the shooting gallery. In this model, Pol II catalytic activity and the rate and processivity of Pol II scanning together with promoter sequence determine the distribution of TSSs and their usage.

Characterization of the genomic distances over which transcription factor (TF) binding influences gene expression is important for inferring target genes from TF chromatin immunoprecipitation followed by sequencing (ChIP-seq) data. Here we systematically examine the relationship between thousands of TF and histone modification ChIP-seq data sets with thousands of gene expression profiles. We develop a model for integrating these data, which reveals two classes of TFs with distinct ranges of regulatory influence, chromatin-binding preferences, and auto-regulatory properties. We find that the regulatory range of the same TF bound within different topologically associating domains (TADs) depend on intrinsic TAD properties such as local gene density and G/C content, but also on the TAD chromatin states. Our results suggest that considering TF type, binding distance to gene locus, as well as chromatin context is important in identifying implicated TFs from GWAS SNPs., Characterization of the distance over which TF binding influences gene expression is important for inferring target genes. Here the authors study this relationship using thousands of genomic data sets, finding two classes of TFs with distinct ranges of regulatory influence modulated by chromatin states of topologically associated domains.

Snorc (Small NOvel Rich in Cartilage) has been identified as a chondrocyte-specific gene in the mouse. Yet little is known about the SNORC protein biochemical properties, and mechanistically how the gene is regulated transcriptionally in a tissue-specific manner. The goals of the present study were to shed light on those important aspects. The chondrocyte nature of Snorc expression was confirmed in mouse and rat tissues, in differentiated (day 7) ATDC5, and in RCS cells where it was constitutive. Topological mapping and biochemical analysis brought experimental evidences that SNORC is a type I protein carrying a chondroitin sulfate (CS) attached to serine 44. The anomalous migration of SNORC on SDS-PAGE was due to its primary polypeptide features, suggesting no additional post-translational modifications apart from the CS glycosaminoglycan. A highly conserved SOX9-binding enhancer located in intron 1 was necessary to drive transcription of Snorc in the mouse, rat, and human. The enhancer was active independently of orientation and whether located in a heterologous promoter or intron. Crispr-mediated inactivation of the enhancer in RCS cells caused reduction of Snorc. Transgenic mice carrying the intronic multimerized enhancer drove high expression of a βGeo reporter in chondrocytes, but not in the hypertrophic zone. Altogether these data confirmed the chondrocyte-specific nature of Snorc and revealed dependency on the intronic enhancer binding of SOX9 for transcription.

Background DNA methylation is an epigenetic modification that mainly repress expression of genes essential during embryogenesis and development. There are key ATPase-dependent enzymes that read or write DNA methylation to remodel chromatin and regulate gene expression. Structural maintenance of chromosome hinge domain containing 1 (SMCHD1) is an architectural protein that regulates expression of numerous genes, some of which are imprinted, that are sensitive to DNA methylation. In addition, SMCHD1 germline mutations lead to developmental diseases; facioscapulohumoral muscular dystrophy (FSHD), bosma arhinia and micropthalmia (BAMS). Current evidence suggests that SMCHD1 functions through maintenance or de novo DNA methylation required for chromatin compaction. However, it is unclear if DNA methylation is also essential for genomic recruitment of SMCHD1 and its role as an architectural protein. We previously isolated SMCHD1 using a methylated DNA region from mouse pituitary growth hormone (Gh1) promoter, suggesting that methylation is required for SMCHD1 DNA binding. The goal of this study was to further understand DNA methylation directed role of SMCHD1 in regulating gene expression. Therefore, we profiled SMCHD1 genome wide occupancy in human neuroblastoma SH-SY5Y cells and evaluated if DNA methylation is required for SMCHD1 genomic binding by treating cells with the DNA demethylating reagent, 5-azacytidine (5-azaC). Results Our data suggest that the majority of SMCHD1 binding occurs in intron and intergenic regions. Gene ontology analysis of genes associated with SMCHD1 genomic occupancy that is sensitive to 5-azaC treatment suggests SMCHD1 involvement in central nervous system development. The potassium voltage-gated channel subfamily Q member1 (KCNQ1) gene that associates with central nervous system is a known SMCHD1 target. We showed SMCHD1 binding to an intronic region of KCNQ1 that is lost following 5-azaC treatment suggesting DNA methylation facilitated binding of SMCHD1. Indeed, deletion of SMCHD1 by CRISPR- Cas9 increases KCNQ1 gene expression confirming its role in regulating KCNQ1 gene expression. Conclusion These findings provide novel insights on DNA methylation directed function of SMCHD1 in regulating expression of genes associated with central nervous system development that impact future drug development strategies.

Nucleosome positioning is critical to chromatin accessibility and is associated with gene expression programs in cells1–3. Previous nucleosome mapping methods assemble profiles from cell populations and reveal a cell-averaged pattern: nucleosomes are positioned and form a phased array that surrounds the transcription start sites of active genes3–6 and DNase I hypersensitive sites7. However, even in a homogenous population of cells, cells exhibit heterogeneity in expression in response to active signalling8,9 that may be related to heterogeneity in chromatin accessibility10–12. Here we report a technique, termed single-cell micrococcal nuclease sequencing (scMNase-seq), that can be used to simultaneously measure genome-wide nucleosome positioning and chromatin accessibility in single cells. Application of scMNase-seq to NIH3T3 cells, mouse primary naive CD4 T cells and mouse embryonic stem cells reveals two principles of nucleosome organization: first, nucleosomes in heterochromatin regions, or that surround the transcription start sites of silent genes, show large variation in positioning across different cells but are highly uniformly spaced along the nucleosome array; and second, nucleosomes that surround the transcription start sites of active genes and DNase I hypersensitive sites show little variation in positioning across different cells but are relatively heterogeneously spaced along the nucleosome array. We found a bimodal distribution of nucleosome spacing at DNase I hypersensitive sites, which corresponds to inaccessible and accessible states and is associated with nucleosome variation and variation in accessibility across cells. Nucleosome variation is smaller within single cells than across cells, and smaller within the same cell type than across cell types. A large fraction of naive CD4 T cells and mouse embryonic stem cells shows depleted nucleosome occupancy at the de novo enhancers detected in their respective differentiated lineages, revealing the existence of cells primed for differentiation to specific lineages in undifferentiated cell populations.

Background Most statistical methods used to identify cancer driver genes are either biased due to choice of assumed parametric models or insensitive to directional relationships important for causal inference. To overcome modeling biases and directional insensitivity, a recent statistical functional chi-squared test (FunChisq) detects directional association via model-free functional dependency. FunChisq examines patterns pointing from independent to dependent variables arising from linear, non-linear, or many-to-one functional relationships. Meanwhile, the Functional Annotation of Mammalian Genome 5 (FANTOM5) project surveyed gene expression at over 200,000 transcription start sites (TSSs) in nearly all human tissue types, primary cell types, and cancer cell lines. The data cover TSSs originated from both coding and noncoding genes. For the vast uncharacterized human TSSs that may exhibit complex patterns in cancer versus normal tissues, the model-free property of FunChisq provides us an unprecedented opportunity to assess the evidence for a gene’s directional effect on human cancer. Results We first evaluated FunChisq and six other methods using 719 curated cancer genes on the FANTOM5 data. FunChisq performed best in detecting known cancer driver genes from non-cancer genes. We also show the capacity of FunChisq to reveal non-monotonic patterns of functional association, to which typical differential analysis methods such as t-test are insensitive. Further applying FunChisq to screen unannotated TSSs in FANTOM5, we predicted 1108 putative cancer driver noncoding RNAs, stronger than 90% of curated cancer driver genes. Next, we compared leukemia samples against other samples in FANTOM5 and FunChisq predicted 332/79 potential biomarkers for lymphoid/myeloid leukemia, stronger than the TSSs of all 87/100 known driver genes in lymphoid/myeloid leukemia. Conclusions This study demonstrated the advantage of FunChisq in revealing directional association, especially in detecting non-monotonic patterns. Here, we also provide the most comprehensive catalog of high-quality biomarkers that may play a causative role in human cancers, including putative cancer driver noncoding RNAs and lymphoid/myeloid leukemia specific biomarkers.

The Drosophila nuclear receptor-binding SET domain protein (NSD) gene encodes the Drosophila ortholog of mammalian NSD family members that are important in many aspects of development and disease in humans. In this study, we observed that overexpression of Drosophila NSD in imaginal discs induces organ atrophy. Thus, to gain an understanding of the transcriptional regulation of the gene, we analyzed the NSD promoter region. First, we identified the presence of three putative DNA replication-related element (DRE) sequences in its promoter region, where DRE-binding factor (DREF) could bind for transcriptional activation. In the experiments with the fly GAL4-UAS system, we demonstrated that overexpressed DREF increased the endogenous NSD transcription. To confirm the role of DREF as a transcriptional activator on the NSD expression, we generated a series of luciferase reporter gene constructs containing deleted portions of the 5'-flanking regions as well as point mutations in the putative DRE sites. When transiently transfected into S2 cells, the deletion construct containing no DRE sites showed dramatic decrease in the NSD promoter activity, but only two sites near the transcriptional start site were important. Furthermore, we verified the direct interaction of DREF with the two positively cis-acting sequences on the NSD promoter by chromatin immunoprecipitation assay. Taken together, these results demonstrated that NSD is one of the downstream targets of the DRE/DREF pathway that is associated with various cellular processes in Drosophila, indicating that our findings may contribute to the understanding of molecular mechanisms in complex disorders associated with NSD family members in humans.

Histone modifications, such as the frequently occurring lysine succinylation, are central to the regulation of chromatin-based processes. However, the mechanism and functional consequences of histone succinylation are unknown. Here we show that the α-ketoglutarate dehydrogenase (α-KGDH) complex is localized in the nucleus in human cell lines and binds to lysine acetyltransferase 2A (KAT2A, also known as GCN5) in the promoter regions of genes. We show that succinyl-coenzyme A (succinyl-CoA) binds to KAT2A. The crystal structure of the catalytic domain of KAT2A in complex with succinyl-CoA at 2.3 Å resolution shows that succinyl-CoA binds to a deep cleft of KAT2A with the succinyl moiety pointing towards the end of a flexible loop 3, which adopts different structural conformations in succinyl-CoA-bound and acetyl-CoA-bound forms. Site-directed mutagenesis indicates that tyrosine 645 in this loop has an important role in the selective binding of succinyl-CoA over acetyl-CoA. KAT2A acts as a succinyltransferase and succinylates histone H3 on lysine 79, with a maximum frequency around the transcription start sites of genes. Preventing the α-KGDH complex from entering the nucleus, or expression of KAT2A(Tyr645Ala), reduces gene expression and inhibits tumour cell proliferation and tumour growth. These findings reveal an important mechanism of histone modification and demonstrate that local generation of succinyl-CoA by the nuclear α-KGDH complex coupled with the succinyltransferase activity of KAT2A is instrumental in histone succinylation, tumour cell proliferation, and tumour development.

This month’s Genome Watch highlights approaches that enable high-resolution mapping of transcription start sites in apicomplexan parasites.

Increased chromatin accessibility is a feature of cell-type-specific cis-regulatory elements; therefore, mapping of DNase I hypersensitive sites (DHSs) enables the detection of active regulatory elements of transcription, including promoters, enhancers, insulators and locus-control regions. Single-cell DNase sequencing (scDNase-seq) is a method of detecting genome-wide DHSs when starting with either single cells or

Xanthomonas campestris pv. campestris (Xcc) is the phytopathogen that causes black rot disease in cruciferous plants. The XCC2366 gene product is annotated as a protein belonging to the TetR family of transcriptional regulators. In this study, we evaluated the function and expression of the XCC2366 gene. Mutational analysis demonstrated that XCC2366 is involved in the resistance to acriflavin and is necessary for virulence in Xcc. In addition, the XCC2366 transcription initiation site was mapped at nucleotide A, 63 nucleotide upstream of the XCC2366 translation start codon. Furthermore, transcriptional analysis revealed that the expression of XCC2366 is induced in the presence of acriflavin. Reporter assay also showed that XCC2366 regulates its own expression under acriflavin-supplemented condition. To the best of our knowledge, acriflavin resistance-related gene in the crucifer pathogen Xcc was characterized for the first time.

Prochlorococcus is a globally abundant marine cyanobacterium with many adaptations that reduce cellular nutrient requirements, facilitating growth in its nutrient-poor environment. One such genomic adaptation is the preferential utilization of amino acids containing fewer N-atoms, which minimizes cellular nitrogen requirements. We predicted that transcriptional regulation might further reduce cellular N budgets during transient N limitation. To explore this, we compared transcription start sites (TSSs) in Prochlorococcus MED4 under N-deprived and N-replete conditions. Of 64 genes with primary and internal TSSs in both conditions, N-deprived cells initiated transcription downstream of primary TSSs more frequently than N-replete cells. Additionally, 117 genes with only an internal TSS demonstrated increased internal transcription under N-deprivation. These shortened transcripts encode predicted proteins with an average of 21% less N content compared to full-length transcripts. We hypothesized that low translation rates, which afford greater control over protein abundances, would be beneficial to relatively slow-growing organisms like Prochlorococcus. Consistent with this idea, we found that Prochlorococcus exhibits greater usage of glycine-glycine motifs, which causes translational pausing, when compared to faster growing microbes. Our findings indicate that structural changes occur within the Prochlorococcus MED4 transcriptome during N-deprivation, potentially altering the size and structure of proteins expressed under nutrient limitation., Gordon and Betty Moore Foundation (Grant GBMF495), Simons Foundation (Award 329108), National Science Foundation (U.S.) (Grant DBI-0424599)

Intergenic regions of divergent gene pairs show bidirectional promoter activity but whether regulatory sequences for gene expression in opposite directions are shared is not established. In this study, promoters of divergently arranged gene pair At4g35640-At4g35650 (SERAT3;2-IDH-III) of Arabidopsis thaliana were analyzed to identify overlapping regulatory regions. Both genes showed the highest expression in flower buds and flowers. 5' RACE experiments extended the intergenic region from 161 bp shown in TAIR annotation to 512 bp. GUS analysis of transgenic A. thaliana plants carrying the 691 bp fragment (512 bp intergenic region plus 5' UTR of both the genes) linked to uidA gene revealed that SERAT3;2 promoter drives gene expression in the tapetum, whereas IDH-III promoter functions specifically in microspores/pollen. Serial 5' deletion of the 691 bp fragment showed SERAT3;2 promoter extends up to -355 position, whereas IDH-III promoter encompasses the 512 bp intergenic region. In transgenics, uidA transcript levels were lower than native SERAT3;2 and IDH-III transcripts indicating presence of additional cis regulatory elements beyond the 691 bp fragment. The present study demonstrated for the first time occurrence of a nested promoter in plants and identified a novel bidirectional promoter capable of driving gene expression in tapetum and microspores/pollen.

RNA polymerase II (Pol II) pauses downstream of the transcription initiation site before beginning productive elongation. This pause is a key component of metazoan gene expression regulation. Some promoters have a strong disposition for Pol II pausing and often mediate faster, more synchronous changes in expression. This requires multiple rounds of transcription and thus cannot rely solely on pause release. However, it is unclear how pausing affects the initiation of new transcripts during consecutive rounds of transcription. Using our recently developed ChIP-nexus method, we find that Pol II pausing inhibits new initiation. We propose that paused Pol II helps prevent new initiation between transcription bursts, which may reduce noise.

Facioscapulohumeral dystrophy (FSHD; OMIM #158900, #158901) is caused by mis-expression of the DUX4 transcription factor in skeletal muscle1. Animal models of FSHD are hampered by incomplete knowledge of the conservation of the DUX4 transcriptional program in other species2–5. Despite divergence of their binding motifs, both mouse Dux and human DUX4 activate genes associated with cleavage-stage embryos, including MERV-L and ERVL-MaLR retrotransposons, in mouse and human muscle cells respectively. When expressed in mouse cells, human DUX4 maintained modest activation of cleavage-stage genes driven by conventional promoters, but did not activate MERV-L-promoted genes. These findings indicate that the ancestral DUX4-factor regulated genes characteristic of cleavage-stage embryos driven by conventional promoters, whereas divergence of the DUX4/Dux homeodomains correlates with retrotransposon specificity. These results provide insight into how species balance conservation of a core transcriptional program with innovation at retrotransposon promoters and provide a basis for animal models that recreate the FSHD transcriptome.

Animal promoters initiate transcription either at precise positions (narrow promoters) or dispersed regions (broad promoters), a distinction referred to as promoter shape. Although highly conserved, the functional properties of promoters with different shapes and the genetic basis of their evolution remain unclear. Here we used natural genetic variation across a panel of 81 Drosophila lines to measure changes in transcriptional start site (TSS) usage, identifying thousands of genetic variants affecting transcript levels (strength) or the distribution of TSSs within a promoter (shape). Our results identify promoter shape as a molecular trait that can evolve independently of promoter strength. Broad promoters typically harbor shape-associated variants, with signatures of adaptive selection. Single-cell measurements demonstrate that variants modulating promoter shape often increase expression noise, whereas heteroallelic interactions with other promoter variants alleviate these effects. These results uncover new functional properties of natural promoters and suggest the minimization of expression noise as an important factor in promoter evolution. Fil: Schor, Ignacio Esteban. European Molecular Biology Laboratory Heidelberg; Alemania. Consejo Nacional de Investigaciones Científicas y Técnicas. Oficina de Coordinación Administrativa Ciudad Universitaria. Instituto de Fisiología, Biología Molecular y Neurociencias. Universidad de Buenos Aires. Facultad de Ciencias Exactas y Naturales. Instituto de Fisiología, Biología Molecular y Neurociencias; Argentina Fil: Degner, Jacob F.. European Molecular Biology Laboratory; Alemania Fil: Harnett, Dermot. European Molecular Biology Laboratory; Alemania Fil: Cannavò, Enrico. European Molecular Biology Laboratory; Alemania Fil: Casale, Francesco P.. European Bioinformatics Institute; Reino Unido Fil: Shim, Heejung. Purdue University; Estados Unidos Fil: Garfield, David A.. European Molecular Biology Laboratory; Alemania Fil: Birney, Ewan. European Bioinformatics Institute; Reino Unido Fil: Stephens, Matthew. University of Chicago; Estados Unidos Fil: Stegle, Oliver. European Bioinformatics Institute; Reino Unido Fil: Furlong, Eileen E. M.. European Molecular Biology Laboratory; Alemania

The Arabidopsis thaliana promoter trap mutant Bitrap-112 expressing green fluorescent protein (GFP) gene in the ovules was found to carry transferred DNA (T-DNA) insertion at -309 position of the APETALA2 (AP2) gene. Bitrap-112 line did not show phenotype associated with the AP2 mutation, suggesting that T-DNA insertion did not interrupt the AP2 promoter. Further, head-to-head orientation of GFP and AP2 genes indicated that the AP2 promoter could be bidirectional. A detailed deletion analysis of the upstream sequences of the AP2 gene was done to identify the promoter. GUS assay of transgenic A. thaliana plants carrying various AP2 upstream fragments fused to the uidA gene showed that ~200-bp 5' UTR sequences are capable of driving gene expression at low levels in vegetative tissues whereas inclusion of further upstream sequences (~300 bp) enhanced uidA expression comparable to native AP2 expression levels in various tissues including ovules. In the reverse orientation, the 519-bp AP2 upstream fragment was found to drive gene expression in immature ovules and pollen. Absence of antisense transcripts corresponding to the sequences upstream of AP2 gene in wild-type A. thaliana plants suggests that promoter trapping has uncovered a cryptic promoter, which in reverse orientation is capable of driving gene expression in ovules and anthers.

The 22-nucleotide non-coding microRNAs (miRNAs) are mostly transcribed by RNA polymerase II and are similar to protein-coding genes. Unlike the clear process from stem-loop precursors to mature miRNAs, the primary transcriptional regulation of miRNA, especially in plants, still needs to be further clarified, including the original transcription start site, functional cis-elements and primary transcript structures. Due to several well-characterized transcription signals in the promoter region, we proposed a systemic approach integrating multidimensional "omics" (including genomics, transcriptomics, and epigenomics) data to improve the genome-wide identification of primary miRNA transcripts. Here, we used the model plant Arabidopsis thaliana to improve the ability to identify candidate promoter locations in intergenic miRNAs and to determine rules for identifying primary transcription start sites of miRNAs by integrating high-throughput omics data, such as the DNase I hypersensitive sites, chromatin immunoprecipitation-sequencing of polymerase II and H3K4me3, as well as high throughput transcriptomic data. As a result, 93% of refined primary transcripts could be confirmed by the primer pairs from a previous study. Cis-element and secondary structure analyses also supported the feasibility of our results. This work will contribute to the primary transcriptional regulatory analysis of miRNAs, and the conserved regulatory pattern may be a suitable miRNA characteristic in other plant species.

The switch from fetal to adult hemoglobin relies on repression or silencing of the upstream γ-globin gene, but identification of the transcriptional repressors that bind to the sites at which a cluster of naturally occurring variants associated with HPFH (hereditary persistence of fetal hemoglobin) are found has been elusive. A new study provides mechanistic evidence for the direct binding of BCL11A and ZBTB7A, two previously identified γ-globin gene repressors.

Background Global RNA sequencing technologies have revealed widespread RNA polymerase II (Pol II) transcription outside of gene promoters. Small 5′-capped RNA sequencing (Start-seq) originally developed for the detection of promoter-proximal Pol II pausing has helped improve annotation of Transcription Start Sites (TSSs) of genes as well as identification of non-genic regulatory elements. However, apart from the most well studied genomes of human and mouse, mammalian transcription has not been profiled with sufficiently high precision. Results We prepared and sequenced Start-seq libraries from rat (Rattus norgevicus) primary neural progenitor cells. Over 48 million uniquely mappable reads from two independent biological replicates allowed us to define the TSSs of 7365 known genes in the rn6 genome, reannotating 2503 TSSs by more than 5 base pairs, characterize promoter-associated antisense transcription, and profile Pol II pausing. By combining TSS data with polyA-selected RNA sequencing, we also identified thousands of potential new genes producing stable RNA as well as non-genic transcripts representing possible regulatory elements. Conclusions Our study has produced the first Start-seq dataset for the rat. Apart from profiling transcription initiation, our data reaffirm the prevalence of Pol II pausing across the rat genome and indicate conservation of pausing mechanisms across metazoan genomes. We suggest that pausing location, at least in mammals, is constrained by a distance from initiation of transcription, whether it occurs at or outside of a gene promoter. Abundant antisense transcription initiation around protein coding genes indicates that Pol II recruited to the vicinity of a promoter is distributed to available start sites of transcription at either DNA strand. Transcriptome profiling of neural progenitors presented here will facilitate further studies of other rat cell types as well as other organisms. Electronic supplementary material The online version of this article (10.1186/s12864-019-5829-4) contains supplementary material, which is available to authorized users.

Background In mammals, fine-tuned regulation of gene expression leads to transcription initiation from diverse transcription start sites (TSSs) and multiple core promoters. Although polysome association is a critical step in translation, whether polysome selectively uses TSSs and core promoters and how this could impact translation remains elusive. Results In this study, we used CAGE followed by deep sequencing to globally profile the transcript 5′ isoforms in the translatome and transcriptome of human HEK293 cells at single-nucleotide resolution. By comparing the two profiles, we identified the 5′ isoforms preferentially used in translatome and revealed a widespread selective usage of TSSs (32.0%) and core promoters (48.7%) by polysome. We discovered the transcription initiation patterns and the sequence characteristics that were highly correlated with polysome selection. We further identified 5804 genes significantly enriched or depleted in translatome and showed that polysome selection was an important contributing factor to the abundance of related gene products. Moreover, after comparison with public transcriptome CAGE data from 180 human tissues and primary cells, we raised a question on whether it is a widely adopted mechanism to regulate translation efficiency by changing the transcription initiation sites on the transcription level in cells of different conditions. Conclusions Using HEK293 cells as a model, we delineated an indirect selection toward TSSs and core promoters by the translation machinery. Our findings lend additional evidence for a much closer coordination between transcription and translation, warranting future translatome studies in more cell types and conditions to develop a more intricate regulatory model for gene expression. Electronic supplementary material The online version of this article (10.1186/s12864-019-5650-0) contains supplementary material, which is available to authorized users.

Hypoxia-inducible factor 2α (HIF2α) plays critical roles in cancer progression. Although the mechanisms of HIF2α translation and degradation have been well studied, the mechanism for HIF2α regulation at transcriptional level is still not fully understood. Here, we present evidence that DNA methylation in promoter contributes to transcription of EPAS1 coding HIF2α. Methylated CpG binding protein 3 (MBD3) contributes to the intricate regulatory mechanism. We showed that MBD3 bound to the EPAS1 promoter in breast cancer cells and amplified EPAS1 transcription through demethylating CpG located around transcriptional start site in MDA-MB-468 cells. This enabled MDA-MB-468 cells to activate HIF2α-mediated angiogenesis. However, in 7860 cells, the demethylation function of MBD3 on EPAS1 was not observed because of the poor methylated-CpG promoter. Nevertheless, depletion of MBD3 induced by shRNA decreased EPAS1 transcription and therefore decreased HIF2α-mediated cellular response in both MDA-MB-468 and 7860 cancer cells. These results indicated that the endogenous MBD3 was involved in regulating the transcription and therefore the transcriptional activities of HIF2α, suggesting that MBD3 may be a potential therapeutic target of tumor.

Galanin receptor 2 (GAL2R) is a G protein-coupled receptor for the neuropeptide galanin that regulates many important physiological functions and pathological processes. To investigate the molecular mechanism governing GAL2R gene transcription, the rat GAL2R promoter was isolated and analyzed. We found that the region from -320 to -300 of the GAL2R promoter contains two putative ETS-1 elements and plays an important role in regulating GAL2R promoter activity. We also showed that transcription factor ETS-1 bound to this region in vitro and in vivo. Overexpression of ETS-1 significantly increased GAL2R promoter activity and transcription of the GAL2R gene, whereas knockdown of ETS-1 produced the opposite effects. In addition, we showed that ETS-1 recruited co-activator p300 to the GAL2R promoter. These data indicate a role for ETS-1 in the control of the GAL2R gene expression and provide a basis for understanding the transcriptional regulation of the GAL2R gene.

Fatty acid binding protein 3 (FABP3) is a member of the FABP family which bind fatty acids and have an important role in fatty acid metabolism. A large number of studies have shown that the genetic polymorphisms of FABP3 are positively correlated with intramuscular fat (IMF) content in domestic animals, however, the function and transcriptional characteristics of FABP3 in cattle remain unclear. Real-time PCR analysis revealed that bovine FABP3 was highly expressed in cardiac tissue. The 5'-regulatory region of bovine FABP3 was cloned and its transcription initiation sites were identified. Sequence analysis showed that many transcriptional factor binding sites including TATA-box and CCAAT-box were present on the 5'-flanking region of bovine FABP3, and four CpG islands were found on nucleotides from -891 to +118. Seven serial deletion constructs of the 5'-regulatory region evaluated in dual-luciferase reporter assay indicated that its core promoter was 384 base pairs upstream from the transcription initiation site. The transcriptional factor binding sites RXRα, KLF15, CREB and Sp1 were conserved in the core promoter of cattle, sheep, pigs and dogs. These results provide further understanding of the function and regulation mechanism of bovine FABP3.

Low phosphorus (P) tolerance in rice is a biologically and agronomically important character. Low P tolerant Indica-type rice genotypes, Sahbhagi Dhan (SD) and Chakhao Poreiton (CP), are adapted to acidic soils and show variable response to low P levels. Using RNAseq approach, transcriptome data was generated from roots of SD and CP after 15 days of low P treatment to understand differences and similarities at molecular level. In response to low P, number of genes up-regulated (1318) was more when compared with down-regulated genes (761). Eight hundred twenty-one genes found to be significantly regulated between SD and CP in response to low P. De novo assembly using plant database led to further identification of 1535 novel transcripts. Functional annotation of significantly expressed genes suggests two distinct methods of low P tolerance. While root system architecture in SD works through serine-threonine kinase PSTOL1, suberin-mediated cell wall modification seems to be key in CP. The transcription data indicated that CP relies more on releasing its internally bound Pi and coping with low P levels by transcriptional and translational modifications and using dehydration response-based signals. Role of P transporters seems to be vital in response to low P in CP while sugar- and auxin-mediated pathway seems to be preferred in SD. At least six small RNA clusters overlap with transcripts highly expressed under low P, suggesting role of RNA super clusters in nutrient response in plants. These results help us to understand and thereby devise better strategy to enhance low P tolerance in Indica-type rice.