Search

To establish a quantitative method to evaluate the DNA methylation level of an immediate upstream region of major BRCA1 transcriptional start sites (TSSs), and to investigate whether methylation of the region is a prognostic factor in high-grade serous ovarian cancer patients after neoadjuvant chemotherapy.Ninety-two FFPE samples of advanced high-grade serous ovarian cancers after neoadjuvant chemotherapy between 2011 and 2018 were used for mutation and methylation analysis. DNA methylation levels were assessed by pyrosequencing and DNA methylation microarray. An association between methylation level (or a mutation) and progression-free survival was assessed by Kaplan-Meier analysis.Major BRCA1 transcripts and CpG sites immediately upstream of their TSSs were identified, and a pyrosequencing method was developed. BRCA1 methylation, BRCA1/2 mutations, and a RAD51C mutation were detected in 17/79 (21.5%), 17/92 (18.5%), and 1/92 (1.1%) high-grade serious ovarian cancer samples. In univariate analysis, BRCA1 methylation and no residual tumor were associated with progression-free survival (BRCA1 methylation: P = 0.025, no residual tumor: P = 0.0026). Multivariate analysis showed that both BRCA1 methylation (P = 0.038, HR = 0.47, 95% CI: 0.21-0.96) and no residual tumor (P = 0.012, HR = 0.49, 95% CI: 0.28-0.85) were significant favorable prognostic factors.A quantitative method to estimate the methylation level of the immediate upstream region of major BRCA1 TSSs was established. Methylation of the region of was an independent favorable prognostic factor in high-grade serous ovarian cancer patients.

RNA polymerase II–mediated eukaryotic transcription starts with the assembly of the preinitiation complex (PIC) on core promoters. The +1 nucleosome is well positioned about 40 base pairs downstream of the transcription start site (TSS) and is commonly known as a barrier of transcription. The +1 nucleosome–bound PIC-Mediator structures show that PIC-Mediator prefers binding to T40N nucleosome located 40 base pairs downstream of TSS and contacts T50N but not the T70N nucleosome. The nucleosome facilitates the organization of PIC-Mediator on the promoter by binding TFIIH subunit p52 and Mediator subunits MED19 and MED26 and may contribute to transcription initiation. PIC-Mediator exhibits multiple nucleosome-binding patterns, supporting a structural role of the +1 nucleosome in the coordination of PIC-Mediator assembly. Our study reveals the molecular mechanism of PIC-Mediator organization on chromatin and underscores the significance of the +1 nucleosome in regulating transcription initiation.

Eukaryotes have evolved multiple ATP-dependent chromatin remodelers to shape the nucleosome landscape. We recently uncovered an evolutionarily conserved SWItch/Sucrose Non-Fermentable (SWI/SNF) chromatin remodeler complex in plants reminiscent of the mammalian BAF subclass, which specifically incorporates the MINUSCULE (MINU) catalytic subunits and the TRIPLE PHD FINGERS (TPF) signature subunits. Here we report experimental evidence that establishes the functional relevance of TPF proteins for the complex activity. Our results show that depletion of TPF triggers similar pleiotropic phenotypes and molecular defects to those found in minu mutants. Moreover, we report the genomic location of MINU2 and TPF proteins as representative members of this SWI/SNF complex and their impact on nucleosome positioning and transcription. These analyses unravel the binding of the complex to thousands of genes where it modulates the position of the +1 nucleosome. These targets tend to produce 5′-shifted transcripts in the tpf and minu mutants pointing to the participation of the complex in alternative transcription start site usage. Interestingly, there is a remarkable correlation between +1 nucleosome shift and 5′ transcript length change suggesting their functional connection. In summary, this study unravels the function of a plant SWI/SNF complex involved in +1 nucleosome positioning and transcription start site determination., MCIN/AEI/10.13039/501100011033 [RYC2018-024108-I, PID2019-108577GA-I00 to J.G.B.]. Funding for open access charge: CSIC Open Access Publication Support Initiative through its Unit of Information Resources for Research (URICI).

The precise and accurate identification and quantification of transcriptional start sites (TSSs) is key to understanding the control of transcription. The core promoter consists of the TSS and proximal non-coding sequences, which are critical in transcriptional regulation. Therefore, the accurate identification of TSSs is important for understanding the molecular regulation of transcription. Existing protocols for TSS identification are challenging and expensive, leaving high-quality data available for a small subset of organisms. This sparsity of data impairs study of TSS usage across tissues or in an evolutionary context. To address these shortcomings, we developed Smart-Seq2 Rolling Circle to Concatemeric Consensus (Smar2C2), which identifies and quantifies TSSs and transcription termination sites. Smar2C2 incorporates unique molecular identifiers that allowed for the identification of as many as 70 million sites, with no known upper limit. We have also generated TSS data sets from as little as 40 pg of total RNA, which was the smallest input tested. In this study, we used Smar2C2 to identify TSSs in Glycine max (soybean), Oryza sativa (rice), Sorghum bicolor (sorghum), Triticum aestivum (wheat) and Zea mays (maize) across multiple tissues. This wide panel of plant TSSs facilitated the identification of evolutionarily conserved features, such as novel patterns in the dinucleotides that compose the initiator element (Inr), that correlated with promoter expression levels across all species examined. We also discovered sequence variations in known promoter motifs that are positioned reliably close to the TSS, such as differences in the TATA box and in the Inr that may prove significant to our understanding and control of transcription initiation. Smar2C2 allows for the easy study of these critical sequences, providing a tool to facilitate discovery.

The prevalent one-dimensional alignment of genomic signals to a reference landmark is a cornerstone of current methods to study transcription and its DNA-dependent processes but it is prone to mask potential relations among multiple DNA elements. We developed a systematic approach to align genomic signals to multiple locations simultaneously by expanding the dimensionality of the genomic-coordinate space. We analyzed transcription in human and uncovered a complex dependence on the relative position of neighboring transcription start sites (TSSs) that is consistently conserved among cell types. The dependence ranges from enhancement to suppression of transcription depending on the relative distances to the TSSs, their intragenic position, and the transcriptional activity of the gene. Our results reveal a conserved hierarchy of alternative TSS usage within a previously unrecognized level of genomic organization and provide a general methodology to analyze complex functional relationships among multiple types of DNA elements.

Promoters are consensus DNA sequences located near the transcription start sites and they play an important role in transcription initiation. Due to their importance in biological processes, the identification of promoters is significantly important for characterizing the expression of the genes. Numerous computational methods have been proposed to predict promoters. However, it is difficult for these methods to achieve satisfactory performance in multiple species. In this study, we propose a novel weighted average ensemble learning model, termed iPro-WAEL, for identifying promoters in multiple species, including Human, Mouse, E.coli, Arabidopsis, B.amyloliquefaciens, B.subtilis and R.capsulatus. Extensive benchmarking experiments illustrate that iPro-WAEL has optimal performance and is superior to the current methods in promoter prediction. The experimental results also demonstrate a satisfactory prediction ability of iPro-WAEL on cross-cell lines, promoters annotated by other methods and distinguishing between promoters and enhancers. Moreover, we identify the most important transcription factor binding site (TFBS) motif in promoter regions to facilitate the study of identifying important motifs in the promoter regions. The source code of iPro-WAEL is freely available at https://github.com/HaoWuLab-Bioinformatics/iPro-WAEL.

CMTR1 (cap methyltransferase 1) catalyses methylation of the first transcribed nucleotide of RNAPII transcripts (N1 2′-O-Me), creating part of the mammalian RNA cap structure. In addition to marking RNA as self, N1 2′-O-Me has ill-defined roles in RNA expression and translation. Here, we investigated the gene specificity of CMTR1 and its impact on RNA expression in embryonic stem cells. Using chromatin immunoprecipitation, CMTR1 was found to bind to transcription start sites (TSS) correlating with RNAPII levels, predominantly binding at histone genes and ribosomal protein (RP) genes. Repression of CMTR1 expression resulted in repression of RNAPII binding at the TSS and repression of RNA expression, particularly of histone and RP genes. In correlation with regulation of histones and RP genes, CMTR1 repression resulted in repression of translation and induction of DNA replication stress and damage. Indicating a direct role for CMTR1 in transcription, addition of recombinant CMTR1 to purified nuclei increased transcription of the histone and RP genes. CMTR1 was found to be upregulated during neural differentiation and there was an enhanced requirement for CMTR1 for gene expression and proliferation during this process. We highlight the distinct roles of the cap methyltransferases RNMT and CMTR1 in target gene expression and differentiation.

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.

The genomic sequence of the horse has been available since 2009, providing critical resources for discovering important genomic variants regarding both animal health and population structures. However, to fully understand the functional implications of these variants, detailed annotation of the horse genome is required. Currently, the horse genome is annotated using limited available RNA-seq data, as well as through comparative genomics by translating human and mouse genome annotation. While this approach has served the equine researchers well and led to various discoveries that would improve care and management of horses, many important questions remain unanswered. The limitation of the current annotation is two pronged. First, the heavy reliance on short-read sequencing-based RNA-seq data meant that alternate isoforms could not be accurately resolved. Second, epigenomic regulatory elements are crucial to detailed understanding of gene expression network but are yet to be systemically identified in the horse. Many of these regulatory elements, including enhancers, promoters, and insulators, are either not transcribed or transcribed at a very low level, necessitating alternate approaches to identify them. To solve above problems, the Functional Annotation of the Animal Genomes (FAANG) project proposed a systemic approach to tissue collection, phenotyping, and data generation, adopting the blueprint laid out by the Encyclopedia of DNA Elements (ENCODE) project. Here we detail the first comprehensive overview of gene expression and regulation in the horse, providing the equine research community an expanded set of genomics resources for studies of complex traits in the horse.Author SummaryFunctional annotation of a reference genome provides critical information that pertains the tissue-specific gene expression and regulation. Non-model organisms often rely on existing annotations of human and mouse genomes and the conservation between species for their genome annotation. This approach has limited power in annotating transcripts and regulatory elements that are less evolutionarily conserved. Such are the cases of alternatively spliced isoforms and enhancer elements. In a large-scale collaborated effort, Functional Annotation of Animal Genome (FAANG) aims to generate species-specific and tissue-aware functional annotation for farm animals. In this study, we present the overall annotation efforts and findings from the equine FAANG group. This integrated annotation for the horse genome provides, for the first time, a comprehensive overview of gene expression and regulation landscape in nine major equine tissues, as well as an analytical framework for further inclusion of other important tissues.

The conserved transcription elongation factor Spt6 makes several contacts with the RNA Polymerase II (RNAPII) complex, including a high-affinity interaction between the Spt6 tandem SH2 domain (Spt6-tSH2) and phosphorylated residues of the Rpb1 subunit in the linker between the catalytic core and the C-terminal domain (CTD) heptad repeats. This interaction contributes to generic localization of Spt6, but we show here that it also has gene-specific roles. Disrupting the interface affected transcription start site selection at a subset of genes whose expression is regulated by this choice, and this was accompanied by changes in a distinct pattern of Spt6 accumulation at these sites. Splicing efficiency was also diminished, as was apparent progression through introns that encode snoRNAs. Chromatin-mediated repression was impaired, and a distinct role in maintaining +1 nucleosomes was identified, especially at ribosomal protein genes. The Spt6-tSH2:Rpb1 interface therefore has both genome-wide functions and local roles at subsets of genes where dynamic decisions regarding initiation, transcript processing, or termination are made. We propose that the interaction modulates the availability or activity of the core elongation and histone chaperone functions of Spt6, contributing to coordination between RNAPII and its accessory factors as varying local conditions call for dynamic responses., Graphical Abstract Graphical AbstractRNAPII is recruited to a promoter and forms a pre-initiation complex. Conversion of the PIC to an elongation complex involves loading multiple factors that promote elongation, modify chromatin, and process the emerging transcript. Many of these factors are associated with the heptad repeats of the CTD of Rpb1 (gray shapes bound to green line). The Spt6 core domain has conserved functions in RNAPII elongation, and we propose that the Spt6-tSH2 domain coordinates these elongation functions with the activities of the factors associated with the Rpb1 CTD by orienting the CTD and mediating interactions among the factors, the other domains of Spt6, and the catalytic engine of RNAPII.

Determination of eukaryotic transcription start sites (TSSs) has been based on methods that require the cap structure at the 5′ end of transcripts derived from Pol II RNA polymerase. Consequently, these methods do not reveal TSSs derived from the other RNA polymerases that also play critical roles in various cell functions. To address this limitation, we developed ReCappable-seq, which comprehensively identifies TSS for both Pol II and non–Pol II transcripts at single-nucleotide resolution. The method relies on specific enzymatic exchange of 5′ m7G caps and 5′ triphosphates with a selectable tag. When applied to human transcriptomes, ReCappable-seq identifies Pol II TSSs that are in agreement with orthogonal methods such as CAGE. Additionally, ReCappable-seq reveals a rich landscape of TSSs associated with Pol III transcripts that have not previously been amenable to study at genome-wide scale. Novel TSS from non–Pol II transcription can be located in the nuclear and mitochondrial genomes. ReCappable-seq interrogates the regulatory landscape of coding and noncoding RNA concurrently and enables the classification of epigenetic profiles associated with Pol II and non–Pol II TSS.

The Saccharomyces cerevisiae gene deletion collection is widely used for functional gene annotation and genetic interaction analyses. However, the standard G418-resistance cassette used to produce knockout mutants delivers strong regulatory elements into the target genetic loci. To date, its side effects on the expression of neighboring genes have never been systematically assessed. Here, using ribosome profiling data, RT-qPCR, and reporter expression, we investigated perturbations induced by the KanMX module. Our analysis revealed significant alterations in the transcription efficiency of neighboring genes and, more importantly, severe impairment of their mRNA translation, leading to changes in protein abundance. In the ‘head-to-head’ orientation of the deleted and neighboring genes, knockout often led to a shift of the transcription start site of the latter, introducing new uAUG codon(s) into the expanded 5′ untranslated region (5′ UTR). In the ‘tail-to-tail’ arrangement, knockout led to activation of alternative polyadenylation signals in the neighboring gene, thus altering its 3′ UTR. These events may explain the so-called neighboring gene effect (NGE), i.e. false genetic interactions of the deleted genes. We estimate that in as much as ∼1/5 of knockout strains the expression of neighboring genes may be substantially (>2-fold) deregulated at the level of translation.

Cutaneous squamous cell carcinomas (SCCs) are frequent heterogeneous tumors arising from sun-exposed regions of the skin and characterized by complex pathogenesis. HOPX is a member of the homeodomain-containing superfamily of proteins holding an atypical homeodomain unable to bind to DNA. First discovered in the heart as a regulator of cardiac development, in the skin, HOPX modulates the terminal differentiation of keratinocytes. There is a particular interest in studying HOPX in squamous skin carcinogenesis because it has the atypical structure and the functional duality as an oncogene and a tumor suppressor gene, reported in different malignancies. In this study, we analyzed the effects of HOPX knockdown and overexpression on SCC tumorigenicity in vitro and in vivo. Our data show that HOPX knockdown in SCC cells inhibits their proliferative and invasive activity through the acceleration of apoptosis. We established that methylation of two alternative HOPX promoters leads to differential expression of HOPX transcripts in normal keratinocytes and SCC cells. Importantly, we report that HOPX acts as an oncogene in the pathogenesis of SCC probably through the activation of the second alternative promoter and the modulation of apoptosis.

GATA3 is a transcription factor that plays an important role in T cell lineage differentiation and T-helper 2 (Th2) type immune responses. In this study, we developed two rat antibodies against Atlantic salmon GATA-3 (anti-r

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.

Germline mutation of PTEN is causally observed in Cowden syndrome (CS) and is one of the most common, penetrant risk genes for autism spectrum disorder (ASD). However, the majority of individuals who present with CS-like clinical features are PTEN-mutation negative. Reassessment of PTEN promoter regulation may help explain abnormal PTEN dosage, as only the minimal promoter and coding regions are currently included in diagnostic PTEN mutation analysis. Therefore, we reanalyzed the architecture of the PTEN promoter using next-generation sequencing datasets. Specifically, run-on sequencing assays identified two additional transcription start regions (TSRs) at −2053 and −1906 basepairs from the canonical start of PTEN, thus extending the PTEN 5'UTR and redefining the PTEN promoter. We show that these novel upstream TSRs are active in cancer cell lines, human cancer and normal tissue. Furthermore, these TSRs can produce novel PTEN transcripts due to the introduction of new splice donors at −2041, −1826 and −1355, which may allow for splicing out of the PTEN 5'UTR or the first and second exon in upstream-initiated transcripts. Combining ENCODE ChIP-seq and pertinent literature, we also compile and analyze all transcription factors (TFs) binding at the redefined PTEN locus. Enrichment analyses suggest that TFs bind specifically to the upstream TSRs may be implicated in inflammatory processes. Altogether, these data redefine the architecture of the PTEN promoter, an important step toward a comprehensive model of PTEN transcription regulation, a basis for future investigations into the new promoters’ role in disease pathogenesis.

Background The widespread usage of Cap Analysis of Gene Expression (CAGE) has led to numerous breakthroughs in understanding the transcription mechanisms. Recent evidence in the literature, however, suggests that CAGE suffers from transcriptional and technical noise. Regardless of the sample quality, there is a significant number of CAGE peaks that are not associated with transcription initiation events. This type of signal is typically attributed to technical noise and more frequently to random five-prime capping or transcription bioproducts. Thus, the need for computational methods emerges, that can accurately increase the signal-to-noise ratio in CAGE data, resulting in error-free transcription start site (TSS) annotation and quantification of regulatory region usage. In this study, we present DeepTSS, a novel computational method for processing CAGE samples, that combines genomic signal processing (GSP), structural DNA features, evolutionary conservation evidence and raw DNA sequence with Deep Learning (DL) to provide single-nucleotide TSS predictions with unprecedented levels of performance. Results To evaluate DeepTSS, we utilized experimental data, protein-coding gene annotations and computationally-derived genome segmentations by chromatin states. DeepTSS was found to outperform existing algorithms on all benchmarks, achieving 98% precision and 96% sensitivity (accuracy 95.4%) on the protein-coding gene strategy, with 96.66% of its positive predictions overlapping active chromatin, 98.27% and 92.04% co-localized with at least one transcription factor and H3K4me3 peak. Conclusions CAGE is a key protocol in deciphering the language of transcription, however, as every experimental protocol, it suffers from biological and technical noise that can severely affect downstream analyses. DeepTSS is a novel DL-based method for effectively removing noisy CAGE signal. In contrast to existing software, DeepTSS does not require feature selection since the embedded convolutional layers can readily identify patterns and only utilize the important ones for the classification task. This study highlights the key role that DL can play in Molecular Biology, by removing the inherent flaws of experimental protocols, that form the backbone of contemporary research. Here, we show how DeepTSS can unleash the full potential of an already popular and mature method such as CAGE, and push the boundaries of coding and non-coding gene expression regulator research even further.

HIV-1 must package its RNA genome to generate infectious viruses. Recent studies have revealed that during genome packaging, HIV-1 not only excludes cellular mRNAs, but also distinguishes among full-length viral RNAs. Using NL4-3 and MAL molecular clones, multiple transcription start sites (TSS) were identified, which generate full-length RNAs that differ by only a few nucleotides at the 5' end. However, HIV-1 selectively packages RNAs containing one guanosine (1G RNA) over RNAs with three guanosines (3G RNA) at the 5' end. Thus, the 5' context of HIV-1 full-length RNA can affect its function. To determine whether the regulation of genome packaging by TSS usage is unique to NL4-3 and MAL, we examined 15 primate lentiviruses including transmitted founder viruses of HIV-1, HIV-2, and several simian immunodeficiency viruses (SIVs). We found that all 15 viruses used multiple TSS to some extent. However, the level of TSS heterogeneity in infected cells varied greatly, even among closely related viruses belonging to the same subtype. Most viruses also exhibited selective packaging of specific full-length viral RNA species into particles. These findings demonstrate that TSS heterogeneity and selective packaging of certain full-length viral RNA species are conserved features of primate lentiviruses. In addition, an SIV strain closely related to the progenitor virus that gave rise to HIV-1 group M, the pandemic pathogen, exhibited TSS usage similar to some HIV-1 strains and preferentially packaged 1G RNA. These findings indicate that multiple TSS usage and selective packaging of a particular unspliced RNA species predate the emergence of HIV-1.

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.

How DNA sequence affects the dynamics and position of RNA Polymerase II (Pol II) during transcription remains poorly understood. Here, we used naturally occurring genetic variation in F1 hybrid mice to explore how DNA sequence differences affect the genome-wide distribution of Pol II. We measured the position and orientation of Pol II in eight organs collected from heterozygous F1 hybrid mice using ChRO-seq. Our data revealed a strong genetic basis for the precise coordinates of transcription initiation and promoter proximal pause, allowing us to redefine molecular models of core transcriptional processes. Our results implicate DNA sequence, including both known and novel DNA sequence motifs, as key determinants of the position of Pol II initiation and pause. We report evidence that initiation site selection follows a stochastic process similar to Brownian motion along the DNA template. We found widespread differences in the position of transcription termination, which impact the primary structure and stability of mature mRNA. Finally, we report evidence that allelic changes in transcription often affect mRNA and ncRNA expression across broad genomic domains. Collectively, we reveal how DNA sequences shape core transcriptional processes at single nucleotide resolution in mammals.

Serum/glucocorticoid-regulated kinase 1 (SGK1) is predominantly expressed in endothelial cells of mouse embryos, and Sgk1 null mice show embryonic lethality due to impaired vascular formation. However, how the SGK1 expression is controlled in developing vasculature remains unknown. In this study, we first identified a proximal endothelial enhancer through lacZ reporter mouse analyses. The mouse Sgk1 proximal enhancer was narrowed down to the 5' region of the major transcription initiation site, while a human corresponding region possessed relatively weak activity. We then searched for distal enhancer candidates using in silico analyses of publicly available databases for DNase accessibility, RNA polymerase association and chromatin modification. A region approximately 500 kb distant from the human SGK1 gene was conserved in the mouse, and the mouse and human genomic fragments drove transcription restricted to embryonic endothelial cells. Minimal fragments of both proximal and distal enhancers had consensus binding elements for the ETS transcription factors, which were essential for the responsiveness to ERG, FLI1 and ETS1 proteins in luciferase assays and the endothelial lacZ reporter expression in mouse embryos. These results suggest that endothelial SGK1 expression in embryonic vasculature is maintained through at least two ETS-regulated enhancers located in the proximal and distal regions.

Engineering microorganisms into biological factories that convert renewable feedstocks into valuable materials is a major goal of synthetic biology; however, for many nonmodel organisms, we do not yet have the genetic tools, such as suites of strong promoters, necessary to effectively engineer them. In this work, we developed a computational framework that can leverage standard RNA-seq data sets to identify sets of constitutive, strongly expressed genes and predict strong promoter signals within their upstream regions. The framework was applied to a diverse collection of RNA-seq data measured for the methanotroph Methylotuvimicrobium buryatense 5GB1 and identified 25 genes that were constitutively, strongly expressed across 12 experimental conditions. For each gene, the framework predicted short (27-30 nucleotide) sequences as candidate promoters and derived -35 and -10 consensus promoter motifs (TTGACA and TATAAT, respectively) for strong expression in M. buryatense. This consensus closely matches the canonical E. coli sigma-70 motif and was found to be enriched in promoter regions of the genome. A subset of promoter predictions was experimentally validated in a XylE reporter assay, including the consensus promoter, which showed high expression. The pmoC, pqqA, and ssrA promoter predictions were additionally screened in an experiment that scrambled the -35 and -10 signal sequences, confirming that transcription initiation was disrupted when these specific regions of the predicted sequence were altered. These results indicate that the computational framework can make biologically meaningful promoter predictions and identify key pieces of regulatory systems that can serve as foundational tools for engineering diverse microorganisms for biomolecule production.

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.

Identifying the tissue-specific molecular signatures of active regulatory elements is critical to understand gene regulatory mechanisms. Here, we identify transcription start sites (TSS) using cap analysis of gene expression (CAGE) across 57 human pancreatic islet samples. We identify 9,954 reproducible CAGE tag clusters (TCs), ~20% of which are islet-specific and occur mostly distal to known gene TSSs. We integrated islet CAGE data with histone modification and chromatin accessibility profiles to identify epigenomic signatures of transcription initiation. Using a massively parallel reporter assay, we validate transcriptional enhancer activity (5% FDR) for 2,279 of 3,378 (~68%) tested islet CAGE elements. TCs within accessible enhancers show higher enrichment to overlap type 2 diabetes genome-wide association study (GWAS) signals than existing islet annotations, which emphasizes the utility of mapping CAGE profiles in disease-relevant tissue. This work provides a high-resolution map of transcriptional initiation in human pancreatic islets with utility for dissecting functional enhancers at GWAS loci.

Enhancers are well established as critical regulators of gene expression, but the mechanisms underlying the molecular basis of their specificity and activity are only partly understood. One of the most exciting recent observations is the discovery of enhancer RNA (eRNA), a class of noncoding RNAs derived from enhancer regions. Transcription of developmentally regulated enhancers has been observed to be associated with their active state. The nature of transcripts (eRNA) and their functional attributes are diverse and context dependent. The majority of eRNA are nonpolyadenylated and present in low abundance owing to their low stability, and may represent transcriptional noise. However, some eRNAs have been reported to be reasonably long and stable, are enriched in nuclei, exhibit tissue‐specific expression and may contribute to enhancer function. Transcription of enhancers has been postulated to mediate enhancer function through either the act of transcription or via the transcribed RNA per se and is a useful feature to be analysed to understand mechanisms underlying enhancer activity. Enhancer Eβ at the murine TCRβ locus has been reported to exhibit enhanced occupancy of RNA polymerase II in developing thymocytes. Here, we investigated the transcriptional potential of Eβ in developing thymocytes and detected overlapping bidirectional transcripts at Eβ ranging between 0.7 and 1.7 kb. These noncoding transcripts are capped, polyadenylated, nuclear and expressed specifically in thymocytes. Delineation of these characteristics is important to further investigate their functional roles in mediating enhancer activity., Several enhancers in their active state exhibit enrichment of H3K4‐me3 and H3K27‐ac, accumulate RNA polymerase II and undergo unidirectional or bidirectional transcription. In this study, enhancer Eβ of the TCRβ locus was found to produce noncoding, convergent and bidirectional transcripts, consistent with its critical regulatory role, in developing thymocytes. It will be interesting to evaluate the functional and mechanistic significance of transcription at enhancer Eβ.

The manner in which newborn coding sequences and their transcriptional competency emerge during the process of gene evolution remains unclear. Here, we experimentally simulated eukaryotic gene origination processes by mimicking horizontal gene transfer events in the plant genome. We mapped the precise position of the transcription start sites (TSSs) of hundreds of newly introduced promoterless firefly luciferase (LUC) coding sequences in the genome of Arabidopsis thaliana cultured cells. The systematic characterization of the LUC-TSSs revealed that 80% of them occurred under the influence of endogenous promoters, while the remainder underwent de novo activation in the intergenic regions, starting from pyrimidine-purine dinucleotides. These de novo TSSs obeyed unexpected rules; they predominantly occurred ∼100 bp upstream of the LUC inserts and did not overlap with Kozak-containing putative open reading frames (ORFs). These features were the output of the immediate responses to the sequence insertions, rather than a bias in the screening of the LUC gene function. Regarding the wild-type genic TSSs, they appeared to have evolved to lack any ORFs in their vicinities. Therefore, the repulsion by the de novo TSSs of Kozak-containing ORFs described above might be the first selection gate for the occurrence and evolution of TSSs in the plant genome. Based on these results, we characterized the de novo type of TSS identified in the plant genome and discuss its significance in genome evolution.

The bacterial σ factor plays the central role in promoter recognition by RNA polymerase (RNAP). The primary σ factor, involved in transcription of housekeeping genes, was also shown to participate in the initiation of RNA synthesis and promoter escape by RNAP. In the open promoter complex, the σ finger formed by σ region 3.2 directly interacts with the template DNA strand upstream of the transcription start site. Here, we analysed the role of the σ finger in transcription initiation by four alternative σ factors in Escherichia coli, σ38, σ32, σ28 and σ24. We found that deletions of the σ finger to various extent compromise the activity of RNAP holoenzymes containing alternative σ factors, especially at low NTP concentrations. All four σs are able to utilize NADH as a noncanonical priming substrate but it has only mild effects on the efficiency of transcription initiation. The mediators of the stringent response, transcription factor DksA and the alarmone ppGpp decrease RNAP activity and promoter complex stability for all four σ factors on tested promoters. For all σs except σ38, deletions of the σ finger conversely increase the stability of promoter complexes and decrease their sensitivity to DksA and ppGpp. The result suggests that the σ finger plays a universal role in transcription initiation by alternative σ factors and sensitizes promoter complexes to the action of global transcription regulators DksA and ppGpp by modulating promoter complex stability.

Paenibacillus durus strain ATCC 35681T is a Gram-positive diazotroph that displayed capability of fixing nitrogen even in the presence of nitrate or ammonium. However, the nitrogen fixation activity was detected only at day 1 of growth when cultured in liquid nitrogen-enriched medium. The transcripts of all the nifH homologues were present throughout the 9-day study. When grown in nitrogen-depleted medium, nitrogenase activities occurred from day 1 until day 6 and the nifH transcripts were also present during the course of the study albeit at different levels. In both studies, the absence of nitrogen fixation activity regardless of the presence of the nifH transcripts raised the possibility of a post-transcriptional or post-translational regulation of the system. A putative SigA box sequence was found upstream of the transcription start site of nifB1, the first gene in the major nitrogen fixation cluster. The upstream region of nifB2 showed a promoter recognisable by SigE, a sigma factor normally involved in sporulation.

Motivation Cell type-specific activities of cis-regulatory elements (CRE) are central to understanding gene regulation and disease predisposition. Single-cell RNA 5′end sequencing (sc-end5-seq) captures the transcription start sites (TSS) which can be used as a proxy to measure the activity of transcribed CREs (tCREs). However, a substantial fraction of TSS identified from sc-end5-seq data may not be genuine due to various artifacts, hindering the use of sc-end5-seq for de novo discovery of tCREs. Results We developed SCAFE—Single-Cell Analysis of Five-prime Ends—a software suite that processes sc-end5-seq data to de novo identify TSS clusters based on multiple logistic regression. It annotates tCREs based on the identified TSS clusters and generates a tCRE-by-cell count matrix for downstream analyses. The software suite consists of a set of flexible tools that could either be run independently or as pre-configured workflows. Availability and implementation SCAFE is implemented in Perl and R. The source code and documentation are freely available for download under the MIT License from https://github.com/chung-lab/SCAFE. Docker images are available from https://hub.docker.com/r/cchon/scafe. The submitted software version and test data are archived at https://doi.org/10.5281/zenodo.7023163 and https://doi.org/10.5281/zenodo.7024060, respectively. Supplementary information Supplementary data are available at Bioinformatics online.

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.

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.

Transcription start site (TSS) usage is a critical factor in the regulation of gene expression. A number of methods for global TSS mapping have been developed, but barriers of expense, technical difficulty, time, and/or cost have limited their broader adoption. To address these issues, we developed Survey of TRanscription Initiation at Promoter Elements with high-throughput sequencing (STRIPE-seq). Requiring only three enzymatic steps with intervening bead cleanups, a STRIPE-seq library can be prepared from as little as 50 ng total RNA in ~5 h at a cost of ~$12 (US). In addition to profiling TSS usage, STRIPE-seq provides information on transcript levels that can be used for differential expression analysis. Thanks to its simplicity and low cost, we envision that STRIPE-seq could be employed by any molecular biology laboratory interested in profiling transcription initiation.

Recognition of mammalian mitochondrial promoters requires the concerted action of mitochondrial RNA polymerase (mtRNAP) and transcription initiation factors TFAM and TFB2M. In this work, we found that transcript slippage results in heterogeneity of the human mitochondrial transcripts in vivo and in vitro. This allowed us to correctly interpret the RNAseq data, identify the bona fide transcription start sites (TSS), and assign mitochondrial promoters for > 50% of mammalian species and some other vertebrates. The divergent structure of the mammalian promoters reveals previously unappreciated aspects of mtDNA evolution. The correct assignment of TSS also enabled us to establish the precise register of the DNA in the initiation complex and permitted investigation of the sequence-specific protein-DNA interactions. We determined the molecular basis of promoter recognition by mtRNAP and TFB2M, which cooperatively recognize bases near TSS in a species-specific manner. Our findings reveal a role of mitochondrial transcription machinery in mitonuclear coevolution and speciation.

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.

Brucella melitensis (B. melitensis) is an important facultative intracellular bacterium that causes global zoonotic diseases. Continuous intracellular survival and replication are the main obstruction responsible for the accessibility of prevention and treatment of brucellosis. Bacteria respond to complex environment by regulating gene expression. Many regulatory factors function at loci where RNA polymerase initiates messenger RNA synthesis. However, limited gene annotation is a current obstacle for the research on expression regulation in bacteria. To improve annotation and explore potential functional sites, we proposed a novel genome-wide method called Capping-seq for transcription start site (TSS) mapping in B. melitensis. This technique combines capture of capped primary transcripts with Single Molecule Real-Time (SMRT) sequencing technology. We identified 2,369 TSSs at single nucleotide resolution by Capping-seq. TSSs analysis of Brucella transcripts showed a preference of purine on the TSS positions. Our results revealed that -35 and -10 elements of promoter contained consensus sequences of TTGNNN and TATNNN, respectively. The 5' ends analysis showed that 57% genes are associated with more than one TSS and 47% genes contain long leader regions, suggested potential complex regulation at the 5' ends of genes in B. melitensis. Moreover, we identified 52 leaderless genes that are mainly involved in the metabolic processes. Overall, Capping-seq technology provides a unique solution for TSS determination in prokaryotes. Our findings develop a systematic insight into the primary transcriptome characterization of B. melitensis. This study represents a critical basis for investigating gene regulation and pathogenesis of Brucella.

Mitochondria are specialized compartments that produce requisite ATP to fuel cellular functions and serve as centers of metabolite processing, cellular signaling, and apoptosis. To accomplish these roles, mitochondria rely on the genetic information in their small genome (mitochondrial DNA) and the nucleus. A growing appreciation for mitochondria's role in a myriad of human diseases, including inherited genetic disorders, degenerative diseases, inflammation, and cancer, has fueled the study of biochemical mechanisms that control mitochondrial function. The mitochondrial transcriptional machinery is different from nuclear machinery. The in vitro re-constituted transcriptional complexes of Saccharomyces cerevisiae (yeast) and humans, aided with high-resolution structures and biochemical characterizations, have provided a deeper understanding of the mechanism and regulation of mitochondrial DNA transcription. In this review, we will discuss recent advances in the structure and mechanism of mitochondrial transcription initiation. We will follow up with recent discoveries and formative findings regarding the regulatory events that control mitochondrial DNA transcription, focusing on those involved in cross-talk between the mitochondria and nucleus. ispartof: JOURNAL OF BIOLOGICAL CHEMISTRY vol:295 issue:52 pages:18406-18425 ispartof: location:United States status: published

Deregulation of microRNA (miRNA) expression plays a critical role in the transition from a physiological to a pathological state. The accurate miRNA promoter identification in multiple cell types is a fundamental endeavor towards understanding and characterizing the underlying mechanisms of both physiological as well as pathological conditions. DIANA-miRGen v4 (www.microrna.gr/mirgenv4) provides cell type specific miRNA transcription start sites (TSSs) for over 1500 miRNAs retrieved from the analysis of >1000 cap analysis of gene expression (CAGE) samples corresponding to 133 tissues, cell lines and primary cells available in FANTOM repository. MiRNA TSS locations were associated with transcription factor binding site (TFBSs) annotation, for >280 TFs, derived from analyzing the majority of ENCODE ChIP-Seq datasets. For the first time, clusters of cell types having common miRNA TSSs are characterized and provided through a user friendly interface with multiple layers of customization. DIANA-miRGen v4 significantly improves our understanding of miRNA biogenesis regulation at the transcriptional level by providing a unique integration of high-quality annotations for hundreds of cell specific miRNA promoters with experimentally derived TFBSs.

The molecular process of transcription by RNA Polymerase II is highly conserved among eukaryotes (“classic model”). Intriguingly, a distinct way of locating transcription start sites (TSSs) was found in a budding yeast Saccharomyces cerevisiae (“scanning model”). The origin of the “scanning model” and its underlying genetic mechanisms remain unsolved. Herein, we applied genomic approaches to address these questions. We first identified TSSs at a single-nucleotide resolution for 12 yeast species using the nAnT-iCAGE technique, which significantly improved the annotations of these genomes by providing accurate 5’boundaries of protein-coding genes. We then infer the initiation mechanism of a species based on its TSS maps and genome sequences. We found that the “scanning model” had originated after the split of Yarrowia lipolytica and the rest of budding yeasts. An adenine-rich region immediately upstream of TSS had appeared during the evolution of the “scanning model” species, which might facilitate TSS selection in these species. Both initiation mechanisms share a strong preference for pyrimidine-purine dinucleotides surrounding the TSS. Our results suggested that the purine is required for accurately recruiting the first nucleotide, increasing the chance of being capped during mRNA maturation, which is critical for efficient translation initiation. Based on our findings, we proposed a model of TSS selection for the “scanning model” species. Besides, our study also demonstrated that the intrinsic sequence feature primarily determines the distribution of initiation activities within a core promoter (core promoter shape).

Bcr‐abl1 oncogene causes a shift in the transcription start site of the SMS1 gene (SGMS1) encoding the sphingomyelin (SM) synthesizing enzyme, sphingomyelin synthase 1 (SMS1). This results in an mRNA with a significantly shorter 5′‐UTR, called 7‐SGMS1, which is translated more efficiently than another transcript (IIb‐SGMS1) with a longer 5′UTR in Bcr‐abl1‐positive cells. Here, we determine the effects of these alternative 5′UTRs on SMS1 translation and investigate the key features underlying such regulation. First, the presence of the longer IIb 5′UTR is sufficient to greatly impair translation of a reporter gene. Deletion of the upstream open reading frame (−164 nt) or of the predicted stem‐loops in the 5′UTR of IIb‐SGMS1 has minimal effects on SGMS1 translation. Conversely, deletion of nucleotides −310 to −132 enhanced transcription of IIb‐SGMS1 to reach that of 7‐SGMS1. We thus suggest that regulatory features within nucleotides −310 and −132 modulate IIb‐SGMS1 translation efficiency.

Ebola virus (EBOV) RNA has the potential to form hairpin structures at the transcription start sequence (TSS) and reinitiation sites of internal genes, both on the genomic and antigenomic/mRNA level. Hairpin formation involving the TSS and the spacer sequence between promotor elements (PE) 1 and 2 was suggested to regulate viral transcription. Here, we provide evidence that such RNA structures form during RNA synthesis by the viral polymerase and affect its activity. This was analysed using monocistronic minigenomes carrying hairpin structure variants in the TSS-spacer region that differ in length and stability. Transcription and replication were measured via reporter activity and by qRT-PCR quantification of the distinct viral RNA species. We demonstrate that viral RNA synthesis is remarkably tolerant to spacer extensions of up to ~54 nt, but declines beyond this length limit (~25% residual activity for a 66-nt extension). Minor incremental stabilizations of hairpin structures in the TSS-spacer region and on the mRNA/antigenomic level were found to rapidly abolish viral polymerase activity, which may be exploited for antisense strategies to inhibit viral RNA synthesis. Finally, balanced viral transcription and replication can still occur when any RNA structure formation potential at the TSS is eliminated, provided that hexamer phasing in the promoter region is maintained. Altogether, the findings deepen and refine our insight into structure and length constraints within the EBOV transcription and replication promoter and suggest a remarkable flexibility of the viral polymerase in recognition of PE1 and PE2.

Estrogen receptor alpha (ERα) is a ligand-dependent transcription regulator, containing two transactivation functional domains, AF-1 and AF-2. The selective estrogen receptor modulators (SERMs), including 4-hydroxytamoxifen (4OHT), activate AF-1 preferentially rather than AF-2. However, it is unclear whether this specific function is related to the tissue-selective functionality of SERMs. Moreover, there is no information determining AF-1 dependent estrogenic-genes existing in tissues. We sought to identify AF-1 dependent estrogenic-genes using the AF-2 mutated knock-in (KI) mouse model, AF2ERKI. AF2ER is an AF-2 disrupted estradiol (E2)-insensitive mutant ERα, but AF-1 dependent transcription can be activated by the estrogen-antagonists, fulvestrant (ICI) and 4OHT. Gene profiling and ChIP-Seq analysis identified Klk1b21 as an ICI-inducible gene in AF2ERKI uterus. The regulatory activity was analyzed further using a cell-based reporter assay. The 5’-flanking 0.4k bp region of Klk1b21 gene responded as an ERα AF-1 dependent estrogen-responsive promoter. The 150 bp minimum ERα binding element (EBE) consists of three direct repeats. These three half-site sequences were essential for the ERα dependent transactivation and were differentially recognized by E2 and 4OHT for the gene activation. This response was impaired when the minimum EBE was fused with a thymidine-kinase promoter but could be restored by fusion with the 100 bp minimum transcription initiation element (TIE) of Klk1b21, suggesting that the cooperative function of EBE and TIE is essential for mediating AF-1 dependent transactivation. These findings provide the first in vivo evidence that endogenous ERα AF-1 dominant estrogenic-genes exist in estrogen-responsive organs. Such findings will aid in understanding the mechanism of ERα dependent tissue-selective activity of SERMs.

The triple-negative breast cancer (TNBC), a subtype of breast cancer which lacks of targeted therapies, exhibits a poor prognosis. It was shown recently that the PIM1 oncogene is highly related to the proliferation of TNBC cells. A quadruplex–duplex hybrid (QDH) forming sequence was recently found to exist near the transcription start site of PIM1. This structure could be an attractive target for regulation of the PIM1 gene expression and thus the treatment of TNBC. Here, we present the solution structures of two QDHs that could coexist in the human PIM1 gene. Form 1 is a three-G-tetrad-layered (3+1) G-quadruplex containing a propeller loop, a lateral loop and a stem-loop made up of three G•C Watson–Crick base pairs. On the other hand, Form 2 is an anti-parallel G-quadruplex comprising two G-tetrads and a G•C•G•C tetrad; the structure has three lateral loops with the middle stem-loop made up of two Watson-Crick G•C base pairs. These structures provide valuable information for the design of G-quadruplex-specific ligands for PIM1 transcription regulation.

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.

The aggregation and joint analysis of large numbers of exome sequences has recently made it possible to de-rive estimates of intolerance to loss-of-function (LoF) variation for human genes. Here, we demonstrate strong and widespread coupling between genic LoF-intolerance and promoter CpG density across the human genome. Genes downstream of the most CpG-rich pro-moters (top 10% CpG density) have a 67.2% probability of being highly LoF-intolerant, using the LOEUF metric from gnomAD. This is in contrast to 7.4% of genes downstream of the most CpG-poor (bottom 10% CpG density) promoters. Combining promoter CpG density with exonic and promoter conservation explains 33.4% of the variation in LOEUF, and the contribution of CpG density exceeds the individual contributions of exonic and promoter conservation. We leverage this to train a simple and easily interpretable predictive model that out-performs other existing predictors and allows us to classify 1,760 genes – which currently lack reliable LOEUF estimates – as highly LoF-intolerant or not. These predictions have the potential to aid in the interpretation of novel patient variants. Moreover, our results reveal that high CpG density is not merely a generic feature of human promoters, but is preferentially encountered at the promoters of the most selectively constrained genes, calling into question the prevailing view that CpG islands are not subject to selection.

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.

mRNA levels are determined by the balance between mRNA synthesis and decay. Protein factors that mediate both processes, including the 5'-3' exonuclease Xrn1, are responsible for a cross-talk between the two processes that buffers steady-state mRNA levels. However, the roles of these proteins in transcription remain elusive and controversial. Applying native elongating transcript sequencing (NET-seq) to yeast cells, we show that Xrn1 functions mainly as a transcriptional activator and that its disruption manifests as a reduction of RNA polymerase II (Pol II) occupancy downstream of transcription start sites. By combining our sequencing data and mathematical modeling of transcription, we found that Xrn1 modulates transcription initiation and elongation of its target genes. Furthermore, Pol II occupancy markedly increased near cleavage and polyadenylation sites in xrn1Δ cells, whereas its activity decreased, a characteristic feature of backtracked Pol II. We also provide indirect evidence that Xrn1 is involved in transcription termination downstream of polyadenylation sites. We noted that two additional decay factors, Dhh1 and Lsm1, seem to function similarly to Xrn1 in transcription, perhaps as a complex, and that the decay factors Ccr4 and Rpb4 also perturb transcription in other ways. Interestingly, the decay factors could differentiate between SAGA- and TFIID-dominated promoters. These two classes of genes responded differently to XRN1 deletion in mRNA synthesis and were differentially regulated by mRNA decay pathways, raising the possibility that one distinction between these two gene classes lies in the mechanisms that balance mRNA synthesis with mRNA decay.

Plants have evolved efficient mechanisms for adapting to temperature fluctuations, known as heat stress response and heat stress memory. Although the transcriptional regulatory network of plant heat stress response has been established, little is known about the genome-wide transcriptional changes occurring within the first several minutes after heat shock. Here, we investigated the nascent RNA and mature messenger RNA (mRNA) from plant leaf tissues exposed to 5 min of heat shock treatment using global run-on sequencing and RNA sequencing methods. Only a small group of genes were up- or downregulated at both the nascent RNA and mRNA levels. Primed plants that were already exposed to mild heat stress exhibited a more drastic alteration at multiple transcriptional steps than naïve plants that had not experienced heat stress. Upon heat shock, we also observed the following: (i) engaged RNA polymerase II accumulated downstream of transcription start sites; (ii) 5' pausing release was a rate-limiting step for the induction of some heat shock protein genes; (iii) numerous genes switched transcription modes; (iv) pervasive read-through was induced at terminators; and (v) heat stress memory occurs at multiple steps of the transcription cycle, such as at Pol II recruitment, 5' pausing, elongation, and termination.

Porcine reproductive and respiratory syndrome (PRRS) is the most economically important infectious disease of pigs worldwide. Our previous study revealed that Tongcheng (TC) pigs display higher resistance to PRRS than Largewhite (LW) pigs, but the genetic mechanism remains unknown. Here, we first confirmed that CXCL14 was downregulated in lungs and porcine alveolar macrophages (PAMs) responding to PRRS virus (PRRSV) infection, but the decline in LW pigs was more obvious than that in TC pigs. Then, we found that the overexpression of CXCL14 activated type-I interferon (IFN-I) signaling by upregulating interferon beta (IFNB), which plays a major role in the antiviral effect. To further decipher the mechanism underlying its differential expression, we characterized the core promoter of CXCL14 as being located from -145 to 276 bp of the transcription start site (TSS) and identified two main haplotypes that displayed significant differential transcriptional activities. We further identified two coupled point mutations that altered the binding status of CEBPB and were responsible for the differential expression in TC and LW pigs. The regulatory effect of CEBPB on CXCL14 was further confirmed by RNA interference (RNAi) and chromatin immunoprecipitation (ChIP), providing crucial clues for deciphering the mechanism of CXCL14 downregulation in unusual conditions. The present study revealed the potential antiviral effect of CXCL14, occurring via activation of interferon signaling, and suggested that CXCL14 contributes to the PRRS resistance of TC pigs.

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.