Search

Patterns of transcriptional activity are encoded in our genome through regulatory elements such as promoters or enhancers that, paradoxically, contain similar assortments of sequence-specific transcription factor (TF) binding sites1-3. Knowledge of how these sequence motifs encode multiple, often overlapping, gene expression programs is central to understanding gene regulation and how mutations in non-coding DNA manifest in disease4,5. Here, by studying gene regulation from the perspective of individual transcription start sites (TSSs), using natural genetic variation, perturbation of endogenous TF protein levels and massively parallel analysis of natural and synthetic regulatory elements, we show that the effect of TF binding on transcription initiation is position dependent. Analysing TF-binding-site occurrences relative to the TSS, we identified several motifs with highly preferential positioning. We show that these patterns are a combination of a TF's distinct functional profiles-many TFs, including canonical activators such as NRF1, NFY and Sp1, activate or repress transcription initiation depending on their precise position relative to the TSS. As such, TFs and their spacing collectively guide the site and frequency of transcription initiation. More broadly, these findings reveal how similar assortments of TF binding sites can generate distinct gene regulatory outcomes depending on their spatial configuration and how DNA sequence polymorphisms may contribute to transcription variation and disease and underscore a critical role for TSS data in decoding the regulatory information of our genome.

Cis-regulatory elements (CREs) can be classified by the shapes of their transcription start site (TSS) profiles, which are indicative of distinct regulatory mechanisms. Massively parallel reporter assays (MPRAs) are increasingly being used to study CRE regulatory mechanisms, yet the degree to which MPRAs replicate individual endogenous TSS profiles has not been determined. Here, we present a new low-input MPRA protocol (TSS-MPRA) that enables measuring TSS profiles of episomal reporters as well as after lentiviral reporter chromatinization. To sensitively compare MPRA and endogenous TSS profiles, we developed a novel dissimilarity scoring algorithm (WIP score) that outperforms the frequently used earth mover's distance on experimental data. Using TSS-MPRA and WIP scoring on 500 unique reporter inserts, we found that short (153 bp) MPRA promoter inserts replicate the endogenous TSS patterns of ∼60% of promoters. Lentiviral reporter chromatinization did not improve fidelity of TSS-MPRA initiation patterns, and increasing insert size frequently led to activation of extraneous TSS in the MPRA that are not active in vivo. We discuss the implications of our findings, which highlight important caveats when using MPRAs to study transcription mechanisms. Finally, we illustrate how TSS-MPRA and WIP scoring can provide novel insights into the impact of transcription factor motif mutations and genetic variants on TSS patterns and transcription levels.

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.

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. Due to the limited availability of functional data for the equine genome, as well as the technical limitations of short-read RNA-seq, existing annotation of the equine genome contains limited information about important aspects of gene regulation, such as alternate isoforms and regulatory elements, which are either not transcribed or transcribed at a very low level. 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, presenting 39,625 novel transcripts, 84,613 candidate cis-regulatory elements (CRE) and their target genes, 332,115 open chromatin regions genome wide across a diverse set of tissues. We showed substantial concordance between chromatin accessibility, chromatin states in different genic features and gene expression. This comprehensive and expanded set of genomics resources will provide the equine research community ample opportunities for studies of complex traits in the horse.

Accurate transcription start site (TSS) annotations are essential for understanding transcriptional regulation and its role in human disease. Gene collections such as GENCODE contain annotations for tens of thousands of TSSs, but not all of these annotations are experimentally validated nor do they contain information on cell type-specific usage. Therefore, we sought to generate a collection of experimentally validated TSSs by integrating RNA Annotation and Mapping of Promoters for the Analysis of Gene Expression (RAMPAGE) data from 115 cell and tissue types, which resulted in a collection of approximately 50 thousand representative RAMPAGE peaks. These peaks are primarily proximal to GENCODE-annotated TSSs and are concordant with other transcription assays. Because RAMPAGE uses paired-end reads, we were then able to connect peaks to transcripts by analyzing the genomic positions of the 3' ends of read mates. Using this paired-end information, we classified the vast majority (37 thousand) of our RAMPAGE peaks as verified TSSs, updating TSS annotations for 20% of GENCODE genes. We also found that these updated TSS annotations are supported by epigenomic and other transcriptomic data sets. To show the utility of this RAMPAGE rPeak collection, we intersected it with the NHGRI/EBI genome-wide association study (GWAS) catalog and identified new candidate GWAS genes. Overall, our work shows the importance of integrating experimental data to further refine TSS annotations and provides a valuable resource for the biological community.

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 rise in throughput and quality of long-read sequencing should allow unambiguous identification of full-length transcript isoforms. However, its application to single-cell RNA-seq has been limited by throughput and expense. Here we develop and characterize long-read Split-seq (LR-Split-seq), which uses combinatorial barcoding to sequence single cells with long reads. Applied to the C2C12 myogenic system, LR-split-seq associates isoforms to cell types with relative economy and design flexibility. We find widespread evidence of changing isoform expression during differentiation including alternative transcription start sites (TSS) and/or alternative internal exon usage. LR-Split-seq provides an affordable method for identifying cluster-specific isoforms in single cells.

The general transcription factor TFIIH contains three ATP-dependent catalytic activities. TFIIH functions in nucleotide excision repair primarily as a DNA helicase and in Pol II transcription initiation as a dsDNA translocase and protein kinase. During initiation, the XPB/Ssl2 subunit of TFIIH couples ATP hydrolysis to dsDNA translocation facilitating promoter opening and the kinase module phosphorylates Pol II to facilitate the transition to elongation. These functions are conserved between metazoans and yeast; however, yeast TFIIH also drives transcription start-site scanning in which Pol II scans downstream DNA to locate productive start-sites. The ten-subunit holo-TFIIH from S. cerevisiae has a processive dsDNA translocase activity required for scanning and a structural role in scanning has been ascribed to the three-subunit TFIIH kinase module. Here, we assess the dsDNA translocase activity of ten-subunit holo- and core-TFIIH complexes (i.e. seven subunits, lacking the kinase module) from both S. cerevisiae and H. sapiens. We find that neither holo nor core human TFIIH exhibit processive translocation, consistent with the lack of start-site scanning in humans. Furthermore, in contrast to holo-TFIIH, the S. cerevisiae core-TFIIH also lacks processive translocation and its dsDNA-stimulated ATPase activity was reduced ~5-fold to a level comparable to the human complexes, potentially explaining the reported upstream shift in start-site observed in vitro in the absence of the S. cerevisiae kinase module. These results suggest that neither human nor S. cerevisiae core-TFIIH can translocate efficiently, and that the S. cerevisiae kinase module functions as a processivity factor to allow for robust transcription start-site scanning.

Female mosquitoes transmit numerous devastating human diseases because they require vertebrate blood meal for egg development. MicroRNAs (miRNAs) play critical roles across multiple reproductive processes in female Aedes aegypti mosquitoes. However, how miRNAs are controlled to coordinate their activity with the demands of mosquito reproduction remains largely unknown. We report that the ecdysone receptor (EcR)-mediated 20-hydroxyecdysone (20E) signaling regulates miRNA expression in female mosquitoes. EcR RNA-interference silencing linked to small RNA-sequencing analysis reveals that EcR not only activates but also represses miRNA expression in the female mosquito fat body, a functional analog of the vertebrate liver. EcR directly represses the expression of clustered miR-275 and miR-305 before blood feeding when the 20E titer is low, whereas it activates their expression in response to the increased 20E titer after a blood meal. Furthermore, we find that SMRTER, an insect analog of the vertebrate nuclear receptor corepressors SMRT and N-CoR, interacts with EcR in a 20E-sensitive manner and is required for EcR-mediated repression of miRNA expression in Ae. aegypti mosquitoes. In addition, we demonstrate that miR-275 and miR-305 directly target glutamate semialdehyde dehydrogenase and AAEL009899, respectively, to facilitate egg development. This study reveals a mechanism for how miRNAs are controlled by the 20E signaling pathway to coordinate their activity with the demands of mosquito reproduction.

Using the Cap Analysis of Gene Expression (CAGE) technology, the FANTOM5 consortium provided one of the most comprehensive maps of transcription start sites (TSSs) in several species. Strikingly, ~72% of them could not be assigned to a specific gene and initiate at unconventional regions, outside promoters or enhancers. Here, we probe these unassigned TSSs and show that, in all species studied, a significant fraction of CAGE peaks initiate at microsatellites, also called short tandem repeats (STRs). To confirm this transcription, we develop Cap Trap RNA-seq, a technology which combines cap trapping and long read MinION sequencing. We train sequence-based deep learning models able to predict CAGE signal at STRs with high accuracy. These models unveil the importance of STR surrounding sequences not only to distinguish STR classes, but also to predict the level of transcription initiation. Importantly, genetic variants linked to human diseases are preferentially found at STRs with high transcription initiation level, supporting the biological and clinical relevance of transcription initiation at STRs. Together, our results extend the repertoire of non-coding transcription associated with DNA tandem repeats and complexify STR polymorphism.

Characterizing transcription start sites is essential for understanding the regulatory mechanisms that control gene expression. Recently, a new bovine genome assembly (ARS-UCD1.2) with high continuity, accuracy, and completeness was released; however, the functional annotation of the bovine genome lacks precise transcription start sites and contains a low number of transcripts in comparison to human and mouse. By using the RAMPAGE approach, this study identified transcription start sites at high resolution in a large collection of bovine tissues. We found several known and novel transcription start sites attributed to promoters of protein-coding and lncRNA genes that were validated through experimental and in silico evidence. With these findings, the annotation of transcription start sites in cattle reached a level comparable to the mouse and human genome annotations. In addition, we identified and characterized transcription start sites for antisense transcripts derived from bidirectional promoters, potential lncRNAs, mRNAs, and pre-miRNAs. We also analyzed the quantitative aspects of RAMPAGE to produce a promoter activity atlas, reaching highly reproducible results comparable to traditional RNA-seq. Coexpression networks revealed considerable use of tissue-specific promoters, especially between brain and testicle, which expressed several genes in common from alternate loci. Furthermore, regions surrounding coexpressed modules were enriched in binding factor motifs representative of each tissue. The comprehensive annotation of promoters in such a large collection of tissues will substantially contribute to our understanding of gene expression in cattle and other mammalian species, shortening the gap between genotypes and phenotypes.

Human genome stability requires efficient repair of oxidized bases, which is initiated via damage recognition and excision by NEIL1 and other base excision repair (BER) pathway DNA glycosylases (DGs). However, the biological mechanisms underlying detection of damaged bases among the million-fold excess of undamaged bases remain enigmatic. Indeed, mutation rates vary greatly within individual genomes, and lesion recognition by purified DGs in the chromatin context is inefficient. Employing super-resolution microscopy and co-immunoprecipitation assays, we find that acetylated NEIL1 (AcNEIL1), but not its non-acetylated form, is predominantly localized in the nucleus in association with epigenetic marks of uncondensed chromatin. Furthermore, chromatin immunoprecipitation followed by high-throughput sequencing (ChIP-seq) revealed non-random AcNEIL1 binding near transcription start sites of weakly transcribed genes and along highly transcribed chromatin domains. Bioinformatic analyses revealed a striking correspondence between AcNEIL1 occupancy along the genome and mutation rates, with AcNEIL1-occupied sites exhibiting fewer mutations compared to AcNEIL1-free domains, both in cancer genomes and in population variation. Intriguingly, from the evolutionarily conserved unstructured domain that targets NEIL1 to open chromatin, its damage surveillance of highly oxidation-susceptible sites to preserve essential gene function and to limit instability and cancer likely originated ∼500 million years ago during the buildup of free atmospheric oxygen.

Sustained, drug-free control of HIV-1 replication is naturally achieved in less than 0.5% of infected individuals (here termed 'elite controllers'), despite the presence of a replication-competent viral reservoir1. Inducing such an ability to spontaneously maintain undetectable plasma viraemia is a major objective of HIV-1 cure research, but the characteristics of proviral reservoirs in elite controllers remain to be determined. Here, using next-generation sequencing of near-full-length single HIV-1 genomes and corresponding chromosomal integration sites, we show that the proviral reservoirs of elite controllers frequently consist of oligoclonal to near-monoclonal clusters of intact proviral sequences. In contrast to individuals treated with long-term antiretroviral therapy, intact proviral sequences from elite controllers were integrated at highly distinct sites in the human genome and were preferentially located in centromeric satellite DNA or in Krüppel-associated box domain-containing zinc finger genes on chromosome 19, both of which are associated with heterochromatin features. Moreover, the integration sites of intact proviral sequences from elite controllers showed an increased distance to transcriptional start sites and accessible chromatin of the host genome and were enriched in repressive chromatin marks. These data suggest that a distinct configuration of the proviral reservoir represents a structural correlate of natural viral control, and that the quality, rather than the quantity, of viral reservoirs can be an important distinguishing feature for a functional cure of HIV-1 infection. Moreover, in one elite controller, we were unable to detect intact proviral sequences despite analysing more than 1.5 billion peripheral blood mononuclear cells, which raises the possibility that a sterilizing cure of HIV-1 infection, which has previously been observed only following allogeneic haematopoietic stem cell transplantation2,3, may be feasible in rare instances.

G-quadruplexes are nucleic acid secondary structures formed by a stack of square planar G-quartets. G-quadruplexes were implicated in many biological functions including telomere maintenance, replication, transcription, and translation, in many species including humans and plants. For wheat, however, though it is one of the world's most important staple food, no G-quadruplex studies have been reported to date. Here, we computationally identify putative G4 structures (G4s) in wheat genome for the first time and compare its distribution across the genome against five other genomes (human, maize, Arabidopsis, rice, and sorghum). We identified close to 1 million G4 motifs with a density of 76 G4s/Mb across the whole genome and 93 G4s/Mb over genic regions. Remarkably, G4s were enriched around three regions, two located on the antisense and one on the sense strand at the following positions: 1) the transcription start site (TSS) (antisense), 2) the first coding domain sequence (CDS) (antisense), and 3) the start codon (sense). Functional enrichment analysis revealed that the gene models containing G4 motifs within these peaks were associated with specific gene ontology (GO) terms, such as developmental process, localization, and cellular component organization or biogenesis. We investigated genes encoding MADS-box transcription factors and showed examples of G4 motifs within critical regulatory regions in the VRN-1 genes in wheat. Furthermore, comparison with other plants showed that monocots share a similar distribution of G4s, but Arabidopsis shows a unique G4 distribution. Our study shows for the first time the prevalence and possible functional roles of G4s in wheat.

Metazoan development requires the robust proliferation of progenitor cells, the identities of which are established by tightly controlled transcriptional networks1. As gene expression is globally inhibited during mitosis, the transcriptional programs that define cell identity must be restarted in each cell cycle2-5 but how this is accomplished is poorly understood. Here we identify a ubiquitin-dependent mechanism that integrates gene expression with cell division to preserve cell identity. We found that WDR5 and TBP, which bind active interphase promoters6,7, recruit the anaphase-promoting complex (APC/C) to specific transcription start sites during mitosis. This allows APC/C to decorate histones with ubiquitin chains branched at Lys11 and Lys48 (K11/K48-branched ubiquitin chains) that recruit p97 (also known as VCP) and the proteasome, which ensures the rapid expression of pluripotency genes in the next cell cycle. Mitotic exit and the re-initiation of transcription are thus controlled by a single regulator (APC/C), which provides a robust mechanism for maintaining cell identity throughout cell division.

MeCP2 plays a multifaceted role in gene expression regulation and chromatin organization. Interaction between MeCP2 and methylated DNA in the regulation of gene expression is well established. However, the widespread distribution of MeCP2 suggests it has additional interactions with chromatin. Here we demonstrate, by both biochemical and genomic analyses, that MeCP2 directly interacts with nucleosomes and its genomic distribution correlates with that of H3K27me3. In particular, the methyl-CpG-binding domain of MeCP2 shows preferential interactions with H3K27me3. We further observe that the impact of MeCP2 on transcriptional changes correlates with histone post-translational modification patterns. Our findings indicate that MeCP2 interacts with genomic loci via binding to DNA as well as histones, and that interaction between MeCP2 and histone proteins plays a key role in gene expression regulation.

Retroelement integration into host genomes affects chromosome structure and function. A goal of a considerable number of investigations is to elucidate features influencing insertion site selection. The Saccharomyces cerevisiae Ty3 retrotransposon inserts proximal to the transcription start sites (TSS) of genes transcribed by RNA polymerase III (RNAP3). In this study, differential patterns of insertion were profiled genome-wide using a random barcode-tagged Ty3. Saturation transposition showed that tRNA genes (tDNAs) are targeted at widely different frequencies even within isoacceptor families. Ectopic expression of Ty3 integrase (IN) showed that it localized to targets independent of other Ty3 proteins and cDNA. IN, RNAP3, and transcription factor Brf1 were enriched at tDNA targets with high frequencies of transposition. To examine potential effects of cis-acting DNA features on transposition, targeting was tested on high-copy plasmids with restricted amounts of 5' flanking sequence plus tDNA. Relative activity of targets was reconstituted in these constructions. Weighting of genomic insertions according to frequency identified an A/T-rich sequence followed by C as the dominant site of strand transfer. This site lies immediately adjacent to the adenines previously implicated in the RNAP3 TSS motif (CAA). In silico DNA structural analysis upstream of this motif showed that targets with elevated DNA curvature coincide with reduced integration. We propose that integration mediated by the Ty3 intasome complex (IN and cDNA) is subject to inputs from a combination of host factor occupancy and insertion site architecture, and that this results in the wide range of Ty3 targeting frequencies.

Determining transcriptional and translational regulatory elements in GC-rich Streptomyces genomes is essential to elucidating the complex regulatory networks that govern secondary metabolite biosynthetic gene cluster (BGC) expression. However, information about such regulatory elements has been limited for Streptomyces genomes. To address this limitation, a high-quality genome sequence of β-lactam antibiotic-producing Streptomyces clavuligerus ATCC 27 064 is completed, which contains 7163 newly annotated genes. This provides a fundamental reference genome sequence to integrate multiple genome-scale data types, including dRNA-Seq, RNA-Seq and ribosome profiling. Data integration results in the precise determination of 2659 transcription start sites which reveal transcriptional and translational regulatory elements, including -10 and -35 promoter components specific to sigma (σ) factors, and 5'-untranslated region as a determinant for translation efficiency regulation. Particularly, sequence analysis of a wide diversity of the -35 components enables us to predict potential σ-factor regulons, along with various spacer lengths between the -10 and -35 elements. At last, the primary transcriptome landscape of the β-lactam biosynthetic pathway is analyzed, suggesting temporal changes in metabolism for the synthesis of secondary metabolites driven by transcriptional regulation. This comprehensive genetic information provides a versatile genetic resource for rational engineering of secondary metabolite BGCs in Streptomyces.

Although promoter-associated CpG islands have been established as targets of DNA methylation changes in cancer, previous studies suggest that epigenetic dysregulation outside the promoter region may be more closely associated with transcriptional changes. Here we examine DNA methylation, chromatin marks, and transcriptional alterations to define the relationship between transcriptional modulation and spatial changes in chromatin structure. Using human papillomavirus-related oropharyngeal carcinoma as a model, we show aberrant enrichment of repressive H3K9me3 at the transcriptional start site (TSS) with methylation-associated, tumor-specific gene silencing. Further analysis identifies a hypermethylated subtype which shows a functional convergence on MYC targets and association with CREBBP/EP300 mutation. The tumor-specific shift to transcriptional repression associated with DNA methylation at TSSs was confirmed in multiple tumor types. Our data may show a common underlying epigenetic dysregulation in cancer associated with broad enrichment of repressive chromatin marks and aberrant DNA hypermethylation at TSSs in combination with MYC network activation.

The potential impact of structural variants includes not only the duplication or deletion of coding sequences, but also the perturbation of noncoding DNA regulatory elements and structural chromatin features, including topological domains (TADs). Structural variants disrupting TAD boundaries have been implicated both in cancer and developmental disease; this likely occurs via "enhancer hijacking," whereby removal of the TAD boundary exposes enhancers to new target transcription start sites (TSSs). With this functional role, we hypothesized that boundaries would display evidence for negative selection. Here we demonstrate that the chromatin landscape constrains structural variation both within healthy humans and across primate evolution. In contrast, in patients with developmental delay, variants occur remarkably uniformly across genomic features, suggesting a potentially broad role for enhancer hijacking in human disease.

CCCTC-binding factor (CTCF) plays a key role in the formation of topologically associating domains (TADs) and loops in interphase. During mitosis TADs are absent, but how TAD formation is dynamically controlled during the cell cycle is not known. Several contradicting observations have been made regarding CTCF binding to mitotic chromatin using both genomics- and microscopy-based techniques. Here, we have used four different assays to address this debate. First, using 5C, we confirmed that TADs and CTCF loops are readily detected in interphase, but absent during prometaphase. Second, ATAC-seq analysis showed that CTCF sites display greatly reduced accessibility and lose the CTCF footprint in prometaphase, suggesting loss of CTCF binding and rearrangement of the nucleosomal array around the binding motif. In contrast, transcription start sites remain accessible in prometaphase, although adjacent nucleosomes can also become repositioned and occupy at least a subset of start sites during mitosis. Third, loss of site-specific CTCF binding was directly demonstrated using CUT&RUN. Histone modifications and histone variants are maintained in mitosis, suggesting a role in bookmarking of active CTCF sites. Finally, live-cell imaging, fluorescence recovery after photobleaching, and single molecule tracking showed that almost all CTCF chromatin binding is lost in prometaphase. Combined, our results demonstrate loss of CTCF binding to CTCF sites during prometaphase and rearrangement of the chromatin landscape around CTCF motifs. This, combined with loss of cohesin, would contribute to the observed loss of TADs and CTCF loops during mitosis and reveals that CTCF sites, key architectural cis-elements, display cell cycle stage-dependent dynamics in factor binding and nucleosome positioning.

The regulation of transcription initiation is critical for developmental and cellular processes. RNA polymerase II (Pol II) is recruited by the basal transcription machinery to the core promoter where Pol II initiates transcription. The core promoter encompasses the region from -40 to +40 bp relative to the +1 transcription start site (TSS). Core promoters may contain one or more core promoter motifs that confer specific properties to the core promoter, such as the TATA box, initiator (Inr) and motifs that are located downstream of the TSS, namely, motif 10 element (MTE), the downstream core promoter element (DPE) and the Bridge, a bipartite core promoter element. We had previously shown that Caudal, an enhancer-binding homeodomain transcription factor and a key regulator of the Hox gene network, is a DPE-specific activator. Interestingly, pair-rule proteins have been implicated in enhancer-promoter communication at the engrailed locus. Fushi tarazu (Ftz) is an enhancer-binding homeodomain transcription factor encoded by the ftz pair-rule gene. Ftz works in concert with its co-factor, Ftz-F1, to activate transcription. Here, we examined whether Ftz and Ftz-F1 activate transcription with a preference for a specific core promoter motif. Our analysis revealed that similarly to Caudal, Ftz and Ftz-F1 activate the promoter containing a TATA box mutation to significantly higher levels than the promoter containing a DPE mutation, thus demonstrating a preference for the DPE motif. We further discovered that Ftz target genes are enriched for a combination of functional downstream core promoter elements that are conserved among Drosophila species. Thus, the unique combination (Inr, Bridge and DPE) of functional downstream core promoter elements within Ftz target genes highlights the complexity of transcriptional regulation via the core promoter in the transcription of different developmental gene regulatory networks.

Although recent studies provide evidence for a common genetic basis between complex traits and Mendelian disorders, a thorough quantification of their overlap in a phenotype-specific manner remains elusive. Here, we have quantified the overlap of genes identified through large-scale genome-wide association studies (GWASs) for 62 complex traits and diseases with genes containing mutations known to cause 20 broad categories of Mendelian disorders. We identified a significant enrichment of genes linked to phenotypically matched Mendelian disorders in GWAS gene sets; of the total 1,240 comparisons, a higher proportion of phenotypically matched or related pairs (n = 50 of 92 [54%]) than phenotypically unmatched pairs (n = 27 of 1,148 [2%]) demonstrated significant overlap, confirming a phenotype-specific enrichment pattern. Further, we observed elevated GWAS effect sizes near genes linked to phenotypically matched Mendelian disorders. Finally, we report examples of GWAS variants localized at the transcription start site or physically interacting with the promoters of genes linked to phenotypically matched Mendelian disorders. Our results are consistent with the hypothesis that genes that are disrupted in Mendelian disorders are dysregulated by non-coding variants in complex traits and demonstrate how leveraging findings from related Mendelian disorders and functional genomic datasets can prioritize genes that are putatively dysregulated by local and distal non-coding GWAS variants.

RNA-sequencing (RNA-seq) is a powerful technique to investigate and quantify entire transcriptomes. Recent advances in the field have made it possible to explore the transcriptomes of single cells. However, most widely used RNA-seq protocols fail to provide crucial information regarding transcription start sites. Here we present a protocol, Tn5Prime, that takes advantage of the Tn5 transposase-based Smart-seq2 protocol to create RNA-seq libraries that capture the 5' end of transcripts. The Tn5Prime method dramatically streamlines the 5' capture process and is both cost effective and reliable. By applying Tn5Prime to bulk RNA and single cell samples, we were able to define transcription start sites as well as quantify transcriptomes at high accuracy and reproducibility. Additionally, similar to 3' end-based high-throughput methods like Drop-seq and 10× Genomics Chromium, the 5' capture Tn5Prime method allows the introduction of cellular identifiers during reverse transcription, simplifying the analysis of large numbers of single cells. In contrast to 3' end-based methods, Tn5Prime also enables the assembly of the variable 5' ends of the antibody sequences present in single B-cell data. Therefore, Tn5Prime presents a robust tool for both basic and applied research into the adaptive immune system and beyond.

Staphylococcus aureus infection is a rising public health care threat. S. aureus is believed to have elaborate regulatory networks that orchestrate its virulence. Despite its importance, the systematic understanding of the transcriptional landscape of S. aureus is limited. Here, we describe the primary transcriptome landscape of an epidemic USA300 isolate of community-acquired methicillin-resistant S. aureus. We experimentally determined 1,861 transcription start sites with their principal promoter elements, including well-conserved -35 and -10 elements and weakly conserved -16 element and 5' untranslated regions containing AG-rich Shine-Dalgarno sequence. In addition, we identified 225 genes whose transcription was initiated from multiple transcription start sites, suggesting potential regulatory functions at transcription level. Along with the transcription unit architecture derived by integrating the primary transcriptome analysis with operon prediction, the measurement of differential gene expression revealed the regulatory framework of the virulence regulator Agr, the SarA-family transcriptional regulators, and β-lactam resistance regulators. Interestingly, we observed a complex interplay between virulence regulation, β-lactam resistance, and metabolism, suggesting a possible tradeoff between pathogenesis and drug resistance in the USA300 strain. Our results provide platform resource for the location of transcription initiation and an in-depth understanding of transcriptional regulation of pathogenesis, virulence, and antibiotic resistance in S. aureus.

Changes in the 5' leader of an mRNA can have profound effects on its translational efficiency with little effect on abundance. Sequencing-based methods to accurately map the 5' leader by identifying the first transcribed nucleotide rely on enzymatic removal of the 5' eukaryotic cap structure by tobacco acid pyrophosphatase (TAP). However, commercial TAP production has been problematic and has now been discontinued. RppH, a bacterial enzyme that can also cleave the 5' cap, and Cap-Clip, a plant-derived enzyme, have been marketed as TAP replacements. We have engineered a Schizosaccharomyces pombe Edc1-fused Dcp1-Dcp2 decapping enzyme that functions as a superior TAP replacement. It can be purified from E. coli overexpression in high yields using standard biochemical methods. This constitutively active enzyme is four orders of magnitude more catalytically efficient than RppH at 5' cap removal, compares favorably to Cap-Clip, and the 5' monophosphorylated RNA product is suitable for standard RNA cloning methods. This engineered enzyme is a better replacement for TAP treatment than the current marketed use of RppH and can be produced cost-effectively in a general laboratory setting, unlike Cap-Clip.

The growth factor Neuregulin-1 (NRG1) has pleiotropic roles in proliferation and differentiation of the stem cell niche in different tissues. It has been implicated in gut, brain and muscle development and repair. Six isoform classes of NRG1 and over 28 protein isoforms have been previously described. Here we report a new class of NRG1, designated NRG1-VII to denote that these NRG1 isoforms arise from a myeloid-specific transcriptional start site (TSS) previously uncharacterized. Long-read sequencing was used to identify eight high-confidence NRG1-VII transcripts. These transcripts presented major structural differences from one another, through the use of cassette exons and alternative stop codons. Expression of NRG1-VII was confirmed in primary human monocytes and tissue resident macrophages and induced pluripotent stem cell-derived macrophages (iPSC-derived macrophages). Isoform switching via cassette exon usage and alternate polyadenylation was apparent during monocyte maturation and macrophage differentiation. NRG1-VII is the major class expressed by the myeloid lineage, including tissue-resident macrophages. Analysis of public gene expression data indicates that monocytes and macrophages are a primary source of NRG1. The size and structure of class VII isoforms suggests that they may be more diffusible through tissues than other NRG1 classes. However, the specific roles of class VII variants in tissue homeostasis and repair have not yet been determined., (© 2024. The Author(s).)

Transcription regulation in cestodes has been little studied. Here, we characterize the Taenia solium TATA-binding protein (TBP) gene. We found binding sites for transcription factors such as NF1, YY1, and AP-1 in the proximal promoter. We also identified two TATA-like elements in the promoter; however, neither could bind TBP. Additionally, we mapped the transcription start site (A+1) within an initiator and identified a putative downstream promoter element (DPE) located at +27 bp relative to the transcription start site. These two elements are important and functional for gene expression. Moreover, we identified the genes encoding T. solium TBP-Associated Factor 6 (TsTAF6) and 9 (TsTAF9). A Western blot assay revealed that both factors are expressed in the parasite; electrophoretic mobility shift assays and super-shift assays revealed interactions between the DPE probe and TsTAF6-TsTAF9. Finally, we used molecular dynamics simulations to formulate an interaction model among TsTAF6, TsTAF9, and the DPE probe; we stabilized the model with interactions between the histone fold domain pair in TAFs and several pairs of nucleotides in the DPE probe. We discuss novel and interesting features of the TsTAF6-TsTAF9 complex for interaction with DPE on T. solium promoters., Competing Interests: The authors declare that they have no conflict of interest., (Copyright: © 2024 Rodríguez-Lima et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.)

Alternative transcription start sites can affect transcript isoform diversity and translation levels. In a recently described form of gene regulation, coordinated transcriptional and translational interference results in transcript isoform-dependent changes in protein expression. Specifically, a long undecoded transcript isoform (LUTI) is transcribed from a gene-distal promoter, interfering with expression of the gene-proximal promoter. Although transcriptional and chromatin features associated with LUTI expression have been described, the mechanism underlying LUTI-based transcriptional interference is not well understood. Using an unbiased genetic approach followed by functional genomics, we uncovered that the Swi/Snf chromatin remodeling complex is required for co-transcriptional nucleosome remodeling that leads to LUTI-based repression. We identified genes with tandem promoters that rely on Swi/Snf function for transcriptional interference during protein folding stress, including LUTI-regulated genes. This study provides clear evidence for Swi/Snf playing a direct role in gene repression via a cis transcriptional interference mechanism., Competing Interests: Declaration of interests The authors declare no competing interests., (Copyright © 2024 The Authors. Published by Elsevier Inc. All rights reserved.)

Background: In vertebrates, most protein-coding genes have a peak of GC-content near their 5' transcriptional start site (TSS). This feature promotes both the efficient nuclear export and translation of mRNAs. Despite the importance of GC-content for RNA metabolism, its general features, origin, and maintenance remain mysterious. We investigate the evolutionary forces shaping GC-content at the transcriptional start site (TSS) of genes through both comparative genomic analysis of nucleotide substitution rates between different species and by examining human de novo mutations., Results: Our data suggests that GC-peaks at TSSs were present in the last common ancestor of amniotes, and likely that of vertebrates. We observe that in apes and rodents, where recombination is directed away from TSSs by PRDM9, GC-content at the 5' end of protein-coding gene is currently undergoing mutational decay. In canids, which lack PRDM9 and perform recombination at TSSs, GC-content at the 5' end of protein-coding is increasing. We show that these patterns extend into the 5' end of the open reading frame, thus impacting synonymous codon position choices., Conclusions: Our results indicate that the dynamics of this GC-peak in amniotes is largely shaped by historic patterns of recombination. Since decay of GC-content towards the mutation rate equilibrium is the default state for non-functional DNA, the observed decrease in GC-content at TSSs in apes and rodents indicates that the GC-peak is not being maintained by selection on most protein-coding genes in those species., (© 2024. The Author(s).)

Mammalian RNA polymerase II preinitiation complexes assemble adjacent to a nucleosome whose proximal edge (NPE) is typically 40 to 50 bp downstream of the transcription start site. At active promoters, that +1 nucleosome is universally modified by trimethylation on lysine 4 of histone H3 (H3K4me3). The Pol II preinitiation complex only extends 35 bp beyond the transcription start site, but nucleosomal templates with an NPE at +51 are nearly inactive in vitro with promoters that lack a TATA element and thus depend on TFIID for promoter recognition. Significantly, this inhibition is relieved when the +1 nucleosome contains H3K4me3, which can interact with TFIID subunits. Here, we show that H3K4me3 templates with both TATA and TATA-less promoters are active with +35 NPEs when transcription is driven by TFIID. Templates with +20 NPE are also active but at reduced levels compared to +35 and +51 NPEs, consistent with a general inhibition of promoter function when the proximal nucleosome encroaches on the preinitiation complex. Remarkably, dinucleosome templates support transcription when H3K4me3 is only present in the distal nucleosome, suggesting that TFIID-H3K4me3 interaction does not require modification of the +1 nucleosome. Transcription reactions performed with an alternative protocol retaining most nuclear factors results primarily in early termination, with a minority of complexes successfully traversing the first nucleosome. In such reactions, the +1 nucleosome does not substantially affect the level of termination even with an NPE of +20, indicating that a nucleosome barrier is not a major driver of early termination by Pol II., Competing Interests: Conflict of interest The authors declare that they have no conflicts of interest with the contents of this article., (Copyright © 2024 The Authors. Published by Elsevier Inc. All rights reserved.)

The sperm epigenome is thought to affect the developmental programming of the resulting embryo, influencing health and disease in later life. Age-related methylation changes in the sperm of old fathers may mediate the increased risks for reproductive and offspring medical problems. The impact of paternal age on sperm methylation has been extensively studied in humans and, to a lesser extent, in rodents and cattle. Here, we performed a comparative analysis of paternal age effects on protein-coding genes in the human and marmoset sperm methylomes. The marmoset has gained growing importance as a non-human primate model of aging and age-related diseases. Using reduced representation bisulfite sequencing, we identified age-related differentially methylated transcription start site (ageTSS) regions in 204 marmoset and 27 human genes. The direction of methylation changes was the opposite, increasing with age in marmosets and decreasing in humans. None of the identified ageTSS was differentially methylated in both species. Although the average methylation levels of all TSS regions were highly correlated between marmosets and humans, with the majority of TSS being hypomethylated in sperm, more than 300 protein-coding genes were endowed with species-specifically (hypo)methylated TSS. Several genes of the glycosphingolipid (GSL) biosynthesis pathway, which plays a role in embryonic stem cell differentiation and regulation of development, were hypomethylated (<5%) in human and fully methylated (>95%) in marmoset sperm. The expression levels and patterns of defined sets of GSL genes differed considerably between human and marmoset pre-implantation embryo stages and blastocyst tissues, respectively., (© 2024 The Authors. Aging Cell published by Anatomical Society and John Wiley & Sons Ltd.)

Drought stress (DS) is one of the major constraints limiting yield in crop plants including rice. Gene regulation under DS is largely governed by accessibility of the transcription factors (TFs) to their cognate cis-regulatory elements (CREs). In this study, we used DNase I hypersensitive assays followed by sequencing to identify the accessible chromatin regions under DS in a drought-sensitive (IR64) and a drought-tolerant (N22) rice cultivar. Our results indicated that DNase I hypersensitive sites (DHSs) were highly enriched at transcription start sites (TSSs) and numerous DHSs were detected in the promoter regions. DHSs were concurrent with epigenetic marks and the genes harboring DHSs in their TSS and promoter regions were highly expressed. In addition, DS induced changes in DHSs (∆DHSs) in TSS and promoter regions were positively correlated with upregulation of several genes involved in drought/abiotic stress response, those encoding TFs and located within drought-associated quantitative trait loci, much preferentially in the drought-tolerant cultivar. The CREs representing the binding sites of TFs involved in DS response were detected within the ∆DHSs, suggesting differential accessibility of TFs to their cognate sites under DS in different rice cultivars, which may be further deployed for enhancing drought tolerance in rice., (© 2024 Society for Experimental Biology and John Wiley & Sons Ltd.)

A universal characteristic of eukaryotic transcription is that the promoter recruits RNA polymerase II (RNAPII) to produce both precursor mRNAs (pre-mRNAs) and short unstable promoter upstream transcripts (PROMPTs) toward the opposite direction. However, how the transcription machinery selects the correct direction to produce pre-mRNAs is largely unknown. Here, through multiple acute auxin-inducible degradation systems, we show that rapid depletion of an RNAPII-binding protein complex, Integrator, results in robust PROMPT accumulation throughout the genome. Interestingly, the accumulation of PROMPTs is compensated by the reduction of pre-mRNA transcripts in actively transcribed genes. Consistently, Integrator depletion alters the distribution of polymerase between the sense and antisense directions, which is marked by increased RNAPII-carboxy-terminal domain Tyr1 phosphorylation at PROMPT regions and a reduced Ser2 phosphorylation level at transcription start sites. Mechanistically, the endonuclease activity of Integrator is critical to suppress PROMPT production. Furthermore, our data indicate that the presence of U1 binding sites on nascent transcripts could counteract the cleavage activity of Integrator. In this process, the absence of robust U1 signal at most PROMPTs allows Integrator to suppress the antisense transcription and shift the transcriptional balance in favor of the sense direction., (© 2024. The Author(s), under exclusive licence to Springer Nature America, Inc.)

The three-dimensional (3D) organization of chromatin within the nucleus is crucial for gene regulation. However, the 3D architectural features that coordinate the activation of an entire chromosome remain largely unknown. We introduce an omics method, RNA-associated chromatin DNA-DNA interactions, that integrates RNA polymerase II (RNAPII)-mediated regulome with stochastic optical reconstruction microscopy to investigate the landscape of noncoding RNA roX2 -associated chromatin topology for gene equalization to achieve dosage compensation. Our findings reveal that roX2 anchors to the target gene transcription end sites (TESs) and spreads in a distinctive boot-shaped configuration, promoting a more open chromatin state for hyperactivation. Furthermore, roX2 arches TES to transcription start sites to enhance transcriptional loops, potentially facilitating RNAPII convoying and connecting proximal promoter-promoter transcriptional hubs for synergistic gene regulation. These TESs cluster as roX2 compartments, surrounded by inactive domains for coactivation of multiple genes within the roX2 territory. In addition, roX2 structures gradually form and scaffold for stepwise coactivation in dosage compensation.

Alternative transcription start site (TSS) usage regulation has been identified as a major means of gene expression regulation in metazoans. However, in fungi, its impact remains elusive as its study has thus far been restricted to model yeasts. Here, we first re-analyzed TSS-seq data to define genuine TSS clusters in 2 species of pathogenic Cryptococcus. We identified 2 types of TSS clusters associated with specific DNA sequence motifs. Our analysis also revealed that alternative TSS usage regulation in response to environmental cues is widespread in Cryptococcus, altering gene expression and protein targeting. Importantly, we performed a forward genetic screen to identify a unique transcription factor (TF) named Tur1, which regulates alternative TSS (altTSS) usage genome-wide when cells switch from exponential phase to stationary phase. ChiP-Seq and DamID-Seq analyses suggest that at some loci, the role of Tur1 might be direct. Tur1 has been previously shown to be essential for virulence in C. neoformans. We demonstrated here that a tur1Δ mutant strain is more sensitive to superoxide stress and phagocytosed more efficiently by macrophages than the wild-type (WT) strain., Competing Interests: The authors have declared that no competing interests exist., (Copyright: © 2024 Dang et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.)

Pausing of RNA polymerase II (Pol II) at transcription start sites (TSSs) primes target genes for productive elongation. Coincidentally, DNA double-strand breaks (DSBs) enrich at highly transcribed and Pol II-paused genes, although their interplay remains undefined. Using androgen receptor (AR) signaling as a model, we have uncovered AR-interacting protein 4 (ARIP4) helicase as a driver of androgen-dependent transcription induction. Chromatin immunoprecipitation sequencing analysis revealed that ARIP4 preferentially co-occupies TSSs with paused Pol II. Moreover, we found that ARIP4 complexes with topoisomerase II beta and mediates transient DSB formation upon hormone stimulation. Accordingly, ARIP4 deficiency compromised release of paused Pol II and resulted in R-loop accumulation at a panel of highly transcribed AR target genes. Last, we showed that ARIP4 binds and unwinds R-loops in vitro and that its expression positively correlates with prostate cancer progression. We propose that androgen stimulation triggers ARIP4-mediated unwinding of R-loops at TSSs, enforcing Pol II pause release to effectively drive an androgen-dependent expression program.

Understanding the regulatory mechanisms of gene expression is a crucial objective in genomics. Although the DNA sequence near the transcription start site (TSS) offers valuable insights, recent methods suggest that analyzing only the surrounding DNA may not suffice to accurately predict gene expression levels. We developed GENet (Gene Expression Network from Histone and Transcription Factor Integration), a novel approach that integrates essential regulatory signals from transcription factors and histone modifications into a graph-based model. GENet extends beyond simple DNA sequence analysis by incorporating additional layers of genetic control, which are vital for determining gene expression. Our method markedly enhances the prediction of mRNA levels compared to previous models that depend solely on DNA sequence data. The results underscore the significance of including comprehensive regulatory information in gene expression studies. GENet emerges as a promising tool for researchers, with potential applications extending from fundamental biological research to the development of medical therapies.

The long noncoding RNA TERRA is transcribed from telomeres in virtually all eukaryotes with linear chromosomes. In humans, TERRA transcription is driven in part by promoters comprising CpG dinucleotide-rich repeats of 29 bp repeats, believed to be present in half of the subtelomeres. Thus far, TERRA expression has been analyzed mainly using molecular biology-based approaches that only generate partial and somehow biased results. Here, we present a novel experimental pipeline to study human TERRA based on long-read sequencing (TERRA ONTseq). By applying TERRA ONTseq to different cell lines, we show that the vast majority of human telomeres produce TERRA and that the cellular levels of TERRA transcripts vary according to their chromosomes of origin. Using TERRA ONTseq, we also identified regions containing TERRA transcription start sites (TSSs) in more than half of human subtelomeres. TERRA TSS regions are generally found immediately downstream from 29 bp repeat-related sequences, which appear to be more widespread than previously estimated. Finally, we isolated a novel TERRA promoter from the highly expressed subtelomere of the long arm of Chromosome 7. With the development of TERRA ONTseq, we provide a refined picture of human TERRA biogenesis and expression and we equip the scientific community with an invaluable tool for future studies., (© 2024 Rodrigues et al.; Published by Cold Spring Harbor Laboratory Press for the RNA Society.)

Our study employs pooled CRISPR screens, integrating 2D and 3D culture models, to identify miRNAs critical in Breast Cancer (BC) tumoursphere formation. These screens combine with RNA-seq experiments allowing identification of miRNA signatures and targets essential for tumoursphere growth. miR-4787-3p exhibits significant up-regulation in BC, particularly in basal-like BCs, suggesting its association with aggressive disease. Surprisingly, despite its location within the 5'UTR of a protein coding gene, which defines DROSHA-independent transcription start site (TSS)-miRNAs, we find it dependant on both DROSHA and DICER1 for maturation. Inhibition of miR-4787-3p hinders tumoursphere formation, highlighting its potential as a therapeutic target in BC. Our study proposes elevated miR-4787-3p expression as a potential prognostic biomarker for adverse outcomes in BC. We find that protein-coding genes positively selected in the CRISPR screens are enriched of miR-4787-3p targets. Of these targets, we select ARHGAP17, FOXO3A, and PDCD4 as known tumour suppressors in cancer and experimentally validate the interaction of miR-4787-3p with their 3'UTRs. Our work illuminates the molecular mechanisms underpinning miR-4787-3p's oncogenic role in BC. These findings advocate for clinical investigations targeting miR-4787-3p and underscore its prognostic significance, offering promising avenues for tailored therapeutic interventions and prognostic assessments in BC., (© 2024. Crown.)

Abnormal transcription initiation from alternative first exon has been reported to promote tumorigenesis. However, the prevalence and impact of gene expression regulation mediated by alternative tandem transcription initiation were mostly unknown in cancer. Here, we developed a robust computational method to analyze alternative tandem transcription start site (TSS) usage from standard RNA sequencing data. Applying this method to pan-cancer RNA sequencing datasets, we observed widespread dysregulation of tandem TSS usage in tumors, many of which were independent of changes in overall expression level or alternative first exon usage. We showed that the dynamics of tandem TSS usage was associated with epigenomic modulation. We found that significant 5' untranslated region shortening of gene TIMM13 contributed to increased protein production, and up-regulation of TIMM13 by CRISPR-mediated transcriptional activation promoted proliferation and migration of lung cancer cells. Our findings suggest that dysregulated tandem TSS usage represents an addtional layer of cancer-associated transcriptome alterations.

Transcribed enhancer maps can reveal nuclear interactions underpinning each cell type and connect specific cell types to diseases. Using a 5' single-cell RNA sequencing approach, we defined transcription start sites of enhancer RNAs and other classes of coding and noncoding RNAs in human CD4 + T cells, revealing cellular heterogeneity and differentiation trajectories. Integration of these datasets with single-cell chromatin profiles showed that active enhancers with bidirectional RNA transcription are highly cell type-specific and that disease heritability is strongly enriched in these enhancers. The resulting cell type-resolved multimodal atlas of bidirectionally transcribed enhancers, which we linked with promoters using fine-scale chromatin contact maps, enabled us to systematically interpret genetic variants associated with a range of immune-mediated diseases.

During spermatogenesis, a large number of germline genes essential for male fertility are coordinately activated. However, it remains unknown how timely activation of this group of germline genes is accomplished. Here we show that Polycomb-repressive complex 1 (PRC1) directs timely activation of germline genes during spermatogenesis. Inactivation of PRC1 in male germ cells results in the gradual loss of a stem cell population and severe differentiation defects, leading to male infertility. In the stem cell population, RNF2, the dominant catalytic subunit of PRC1, activates transcription of Sall4, which codes for a transcription factor essential for subsequent spermatogenic differentiation. Furthermore, RNF2 and SALL4 together occupy transcription start sites of germline genes in the stem cell population. Once differentiation commences, these germline genes are activated to enable the progression of spermatogenesis. Our study identifies a novel mechanism by which Polycomb directs the developmental process by activating a group of lineage-specific genes.

N6-methyldeoxyadenine (6mA) is a noncanonical DNA base modification present at low levels in plant and animal genomes, but its prevalence and association with genome function in other eukaryotic lineages remains poorly understood. Here we report that abundant 6mA is associated with transcriptionally active genes in early-diverging fungal lineages. Using single-molecule long-read sequencing of 16 diverse fungal genomes, we observed that up to 2.8% of all adenines were methylated in early-diverging fungi, far exceeding levels observed in other eukaryotes and more derived fungi. 6mA occurred symmetrically at ApT dinucleotides and was concentrated in dense methylated adenine clusters surrounding the transcriptional start sites of expressed genes; its distribution was inversely correlated with that of 5-methylcytosine. Our results show a striking contrast in the genomic distributions of 6mA and 5-methylcytosine and reinforce a distinct role for 6mA as a gene-expression-associated epigenomic mark in eukaryotes.

The resolution and formation of facultative heterochromatin are essential for development, reprogramming, and oncogenesis. The mechanisms underlying these changes are poorly understood owing to the difficulty of studying heterochromatin dynamics and structure in vivo. We devised an in vivo approach to investigate these mechanisms and found that topoisomerase II (TOP2), but not TOP1, synergizes with BAF (mSWI/SNF) ATP-dependent chromatin remodeling complexes genome-wide to resolve facultative heterochromatin to accessible chromatin independent of transcription. This indicates that changes in DNA topology that take place through (de-)catenation rather than the release of torsional stress through swiveling are necessary for heterochromatin resolution. TOP2 and BAF cooperate to recruit pluripotency factors, which explains some of the instructive roles of BAF complexes. Unexpectedly, we found that TOP2 also plays a role in the re-formation of facultative heterochromatin; this finding suggests that facultative heterochromatin and accessible chromatin exist at different states of catenation or other topologies, which might be critical to their structures.

DNA sequence signals in the core promoter, such as the initiator (Inr), direct transcription initiation by RNA polymerase II. Here we show that the human Inr has the consensus of BBCA+1BW at focused promoters in which transcription initiates at a single site or a narrow cluster of sites. The analysis of 7678 focused transcription start sites revealed 40% with a perfect match to the Inr and 16% with a single mismatch outside of the CA+1 core. TATA-like sequences are underrepresented in Inr promoters. This consensus is a key component of the DNA sequence rules that specify transcription initiation in humans.

DNA sequence signals in the core promoter, such as the initiator (Inr), direct transcription initiation by RNA polymerase II. Here we show that the human Inr has the consensus of BBCA+1BW at focused promoters in which transcription initiates at a single site or a narrow cluster of sites. The analysis of 7678 focused transcription start sites revealed 40% with a perfect match to the Inr and 16% with a single mismatch outside of the CA+1 core. TATA-like sequences are underrepresented in Inr promoters. This consensus is a key component of the DNA sequence rules that specify transcription initiation in humans.

The Cancer Genome Atlas (TCGA) network study of 12 cancer types (PanCancer 12) revealed frequent mutation of TP53, and amplification and expression of related TP63 isoform ΔNp63 in squamous cancers. Further, aberrant expression of inflammatory genes and TP53/p63/p73 targets were detected in the PanCancer 12 project, reminiscent of gene programs comodulated by cREL/ΔNp63/TAp73 transcription factors we uncovered in head and neck squamous cell carcinomas (HNSCCs). However, how inflammatory gene signatures and cREL/p63/p73 targets are comodulated genome wide is unclear. Here, we examined how the inflammatory factor tumor necrosis factor-α (TNF-α) broadly modulates redistribution of cREL with ΔNp63α/TAp73 complexes and signatures genome wide in the HNSCC model UM-SCC46 using chromatin immunoprecipitation sequencing (ChIP-seq). TNF-α enhanced genome-wide co-occupancy of cREL with ΔNp63α on TP53/p63 sites, while unexpectedly promoting redistribution of TAp73 from TP53 to activator protein-1 (AP-1) sites. cREL, ΔNp63α and TAp73 binding and oligomerization on NF-κB-, TP53- or AP-1-specific sequences were independently validated by ChIP-qPCR (quantitative PCR), oligonucleotide-binding assays and analytical ultracentrifugation. Function of the binding activity was confirmed using TP53-, AP-1- and NF-κB-specific REs or p21, SERPINE1 and IL-6 promoter luciferase reporter activities. Concurrently, TNF-α regulated a broad gene network with cobinding activities for cREL, ΔNp63α and TAp73 observed upon array profiling and reverse transcription-PCR. Overlapping target gene signatures were observed in squamous cancer subsets and in inflamed skin of transgenic mice overexpressing ΔNp63α. Furthermore, multiple target genes identified in this study were linked to TP63 and TP73 activity and increased gene expression in large squamous cancer samples from PanCancer 12 TCGA by CircleMap. PARADIGM inferred pathway analysis revealed the network connection of TP63 and NF-κB complexes through an AP-1 hub, further supporting our findings. Thus, inflammatory cytokine TNF-α mediates genome-wide redistribution of the cREL/p63/p73, and AP-1 interactome, to diminish TAp73 tumor suppressor function and reciprocally activate NF-κB and AP-1 gene programs implicated in malignancy.

The canonical view of DNA methylation, a pivotal epigenetic regulation mechanism in eukaryotes, dictates its role as a suppressor of gene activity, particularly within promoter regions. However, this view is being challenged as it is becoming increasingly evident that the connection between DNA methylation and gene expression varies depending on the genomic location and is therefore more complex than initially thought. We examined DNA methylation levels in the gut epithelium of Atlantic salmon ( Salmo salar ) using whole-genome bisulfite sequencing, which we correlated with gene expression data from RNA sequencing of the same gut tissue sample (RNA-seq). Assuming epigenetic signals might be pronounced between distinctive phenotypes, we compared large and small fish, finding 22 significant associations between 22 differentially methylated regions and 21 genes. We did not detect significant methylation differences between large and small fish. However, we observed a consistent signal of methylation levels around the transcription start sites (TSS), being negatively correlated with the expression levels of those genes. We found both negative and positive associations of methylation levels with gene expression further upstream or downstream of the TSS, revealing a more unpredictable pattern. The 21 genes showing significant methylation-expression correlations were involved in biological processes related to salmon health, such as growth and immune responses. Deciphering how DNA methylation affects the expression of such genes holds great potential for future applications. For instance, our results suggest the importance of genomic context in targeting epigenetic modifications to improve the welfare of aquaculture species like Atlantic salmon.

Nuclear pore complexes (NPCs) emerged as nuclear transport channels in eukaryotic cells ∼1.5 billion years ago. While the primary role of NPCs is to regulate nucleo-cytoplasmic transport, recent research suggests that certain NPC proteins have additionally acquired the role of affecting gene expression at the nuclear periphery and in the nucleoplasm in metazoans. Here we identify a widely expressed variant of the transmembrane nucleoporin (Nup) Pom121 (named sPom121, for "soluble Pom121") that arose by genomic rearrangement before the divergence of hominoids. sPom121 lacks the nuclear membrane-anchoring domain and thus does not localize to the NPC. Instead, sPom121 colocalizes and interacts with nucleoplasmic Nup98, a previously identified transcriptional regulator, at gene promoters to control transcription of its target genes in human cells. Interestingly, sPom121 transcripts appear independently in several mammalian species, suggesting convergent innovation of Nup-mediated transcription regulation during mammalian evolution. Our findings implicate alternate transcription initiation as a mechanism to increase the functional diversity of NPC components.