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.

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.

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 p53 family of proteins evolved from a common ancestor into three separate genes encoding proteins that act as transcription factors with distinct cellular roles. Isoforms of each member that lack specific regions or domains are suggested to result from alternative transcription start sites, alternative splicing or alternative translation initiation, and have the potential to exponentially increase the functional repertoire of each gene. However, evidence supporting the presence of individual protein variants at functional levels is often limited and is inferred by mRNA detection using highly sensitive amplification techniques. We provide a critical appraisal of the current evidence for the origins, expression, functions and regulation of p53-family isoforms. We conclude that despite the wealth of publications, several putative isoforms remain poorly established. Future research with improved technical approaches and the generation of isoform-specific protein detection reagents is required to establish the physiological relevance of p53-family isoforms in health and disease. In addition, our analyses suggest that p53-family variants evolved partly through convergent rather than divergent evolution from the ancestral gene., (© The Author(s) 2024. Published by Oxford University Press on behalf of Nucleic Acids Research.)

RNA polymerase II (pol II) initiates transcription from transcription start sites (TSSs) located ∼30-35 bp downstream of the TATA box in metazoans, whereas in the yeast Saccharomyces cerevisiae, pol II scans further downstream TSSs located ∼40-120 bp downstream of the TATA box. Previously, we found that removal of the kinase module TFIIK (Kin28-Ccl1-Tfb3) from TFIIH shifts the TSS in a yeast in vitro system upstream to the location observed in metazoans and that addition of recombinant Tfb3 back to TFIIH-ΔTFIIK restores the downstream TSS usage. Here, we report that this biochemical activity of yeast TFIIK in TSS scanning is attributable to the Tfb3 RING domain at the interface with pol II in the pre-initiation complex (PIC): especially, swapping Tfb3 Pro51-a residue conserved among all fungi-with Ala or Ser as in MAT1, the metazoan homolog of Tfb3, confers an upstream TSS shift in vitro in a similar manner to the removal of TFIIK. Yeast genetic analysis suggests that both Pro51 and Arg64 of Tfb3 are required to maintain the stability of the Tfb3-pol II interface in the PIC. Cryo-electron microscopy analysis of a yeast PIC lacking TFIIK reveals considerable variability in the orientation of TFIIH, which impairs TSS scanning after promoter opening., (© The Author(s) 2024. Published by Oxford University Press on behalf of Nucleic Acids Research.)

Fanconi anemia (FA) is a rare genetic disease characterized by congenital abnormalities and increased risk for bone marrow failure and cancer. Central nervous system defects, including acute and irreversible loss of neurological function and white matter lesions with calcifications, have become increasingly recognized among FA patients, and are collectively referred to as Fanconi Anemia Neurological Syndrome or FANS. The molecular etiology of FANS is poorly understood. In this study, we have used a functional integrative genomics approach to further define the function of the FANCD2 protein and FA pathway. Combined analysis of new and existing FANCD2 ChIP-seq datasets demonstrates that FANCD2 binds nonrandomly throughout the genome with binding enriched at transcription start sites and in broad regions spanning protein-coding gene bodies. FANCD2 demonstrates a strong preference for large neural genes involved in neuronal differentiation, synapse function, and cell adhesion, with many of these genes implicated in neurodevelopmental and neuropsychiatric disorders. Furthermore, FANCD2 binds to regions of the genome that replicate late, undergo mitotic DNA synthesis (MiDAS) under conditions of replication stress, and are hotspots for copy number variation. Our analysis describes an important targeted role for FANCD2 and the FA pathway in the maintenance of large neural gene stability., (© 2024. The Author(s).)

Much of what we know about eukaryotic transcription stems from animals and yeast; however, plants evolved separately for over a billion years, leaving ample time for divergence in transcriptional regulation. Here we set out to elucidate fundamental properties of cis-regulatory sequences in plants. Using massively parallel reporter assays across four plant species, we demonstrate the central role of sequences downstream of the transcription start site (TSS) in transcriptional regulation. Unlike animal enhancers that are position independent, plant regulatory elements depend on their position, as altering their location relative to the TSS significantly affects transcription. We highlight the importance of the region downstream of the TSS in regulating transcription by identifying a DNA motif that is conserved across vascular plants and is sufficient to enhance gene expression in a dose-dependent manner. The identification of a large number of position-dependent enhancers points to fundamental differences in gene regulation between plants and animals., (© 2024. The Author(s).)

HIV-1, the causative agent of AIDS, is a retrovirus that packages two copies of unspliced viral RNA as a dimer into newly budding virions. The unspliced viral RNA also serves as an mRNA template for translation of two polyproteins. Recent studies suggest that the fate of the viral RNA (genome or mRNA) is determined at the level of transcription. RNA polymerase II uses heterogeneous transcription start sites to generate major transcripts that differ in only two guanosines at the 5' end. Remarkably, this two-nucleotide difference is sufficient to alter the structure of the 5'-untranslated region and generate two pools of RNA with distinct functions. The presence of both RNA species is needed for optimal viral replication and fitness., Competing Interests: Declaration of competing interest None., (Published by Elsevier Ltd.)

Alternative transcription start sites (TSS) are widespread in eukaryotes and can alter the 5' UTR length and coding potential of transcripts. Here we show that inorganic phosphate (Pi) availability regulates the usage of several alternative TSS in Arabidopsis (Arabidopsis thaliana). In comparison to phytohormone treatment, Pi had a pronounced and specific effect on the usage of many alternative TSS. By combining short-read RNA sequencing with long-read sequencing of full-length mRNAs, we identified a set of 45 genes showing alternative TSS under Pi deficiency. Alternative TSS affected several processes, such as translation via the exclusion of upstream open reading frames present in the 5' UTR of RETICULAN LIKE PROTEIN B1 mRNA, and subcellular localization via removal of the plastid transit peptide coding region from the mRNAs of HEME OXYGENASE 1 and SULFOQUINOVOSYLDIACYLGLYCEROL 2. Several alternative TSS also generated shorter transcripts lacking the coding potential for important domains. For example, the EVOLUTIONARILY CONSERVED C-TERMINAL REGION 4 (ECT4) locus, which encodes an N 6 -methyladenosine (m 6 A) reader, strongly expressed under Pi deficiency a short noncoding transcript (named ALT ECT4 ) ~550 nt long with a TSS in the penultimate intron. The specific and robust induction of ALT ECT4 production by Pi deficiency led to the identification of a role for m 6 A readers in primary root growth in response to low phosphate that is dependent on iron and is involved in modulating cell division in the root meristem. Our results identify alternative TSS usage as an important process in the plant response to Pi deficiency., (© 2024 The Author(s). The Plant Journal published by Society for Experimental Biology and John Wiley & Sons Ltd.)

In this study, we employed short- and long-read sequencing technologies to delineate the transcriptional architecture of the human monkeypox virus and to identify key regulatory elements that govern its gene expression. Specifically, we conducted a transcriptomic analysis to annotate the transcription start sites (TSSs) and transcription end sites (TESs) of the virus by utilizing Cap Analysis of gene expression sequencing on the Illumina platform and direct RNA sequencing on the Oxford Nanopore technology device. Our investigations uncovered significant complexity in the use of alternative TSSs and TESs in viral genes. In this research, we also detected the promoter elements and poly(A) signals associated with the viral genes. Additionally, we identified novel genes in both the left and right variable regions of the viral genome.IMPORTANCEGenerally, gaining insight into how the transcription of a virus is regulated offers insights into the key mechanisms that control its life cycle. The recent outbreak of the human monkeypox virus has underscored the necessity of understanding the basic biology of its causative agent. Our results are pivotal for constructing a comprehensive transcriptomic atlas of the human monkeypox virus, providing valuable resources for future studies., Competing Interests: The authors declare no conflict of interest.

Differentiation of female germline stem cells into a mature oocyte includes the expression of RNAs and proteins that drive early embryonic development in Drosophila . We have little insight into what activates the expression of these maternal factors. One candidate is the zinc-finger protein OVO. OVO is required for female germline viability and has been shown to positively regulate its own expression, as well as a downstream target, ovarian tumor , by binding to the transcriptional start site (TSS). To find additional OVO targets in the female germline and further elucidate OVO's role in oocyte development, we performed ChIP-seq to determine genome-wide OVO occupancy, as well as RNA-seq comparing hypomorphic and wild type rescue ovo alleles. OVO preferentially binds in close proximity to target TSSs genome-wide, is associated with open chromatin, transcriptionally active histone marks, and OVO-dependent expression. Motif enrichment analysis on OVO ChIP peaks identified a 5'-TAACNGT-3' OVO DNA binding motif spatially enriched near TSSs. However, the OVO DNA binding motif does not exhibit precise motif spacing relative to the TSS characteristic of RNA polymerase II complex binding core promoter elements. Integrated genomics analysis showed that 525 genes that are bound and increase in expression downstream of OVO are known to be essential maternally expressed genes. These include genes involved in anterior/posterior/germ plasm specification ( bcd, exu, swa, osk, nos, aub, pgc, gcl ), egg activation ( png, plu, gnu, wisp, C(3)g, mtrm ), translational regulation ( cup , orb , bru1, me31B ), and vitelline membrane formation ( fs(1)N , fs(1)M3 , clos ). This suggests that OVO is a master transcriptional regulator of oocyte development and is responsible for the expression of structural components of the egg as well as maternally provided RNAs that are required for early embryonic development., Competing Interests: LB, SM, JG, BO No competing interests declared

Mycobacterium tuberculosis (MTB) is a pathogen that is known for its ability to persist in harsh environments and cause chronic infections. Understanding the regulatory networks of MTB is crucial for developing effective treatments. Small regulatory RNAs (sRNAs) play important roles in gene expression regulation in all kingdoms of life, and their classification based solely on genomic location can be imprecise due to the computational-based prediction of protein-coding genes in bacteria, which often neglects segments of mRNA such as 5'UTRs, 3'UTRs, and intercistronic regions of operons. To address this issue, our study simultaneously discovered genomic features such as TSSs, UTRs, and operons together with sRNAs in the M. tuberculosis H37Rv strain (ATCC 27294) across multiple stress conditions. Our analysis identified 1,376 sRNA candidates and 8,173 TSSs in MTB, providing valuable insights into its complex regulatory landscape. TSS mapping enabled us to classify these sRNAs into more specific categories, including promoter-associated sRNAs, 5'UTR-derived sRNAs, 3'UTR-derived sRNAs, true intergenic sRNAs, and antisense sRNAs. Three of these sRNA candidates were experimentally validated using 3'-RACE-PCR: predictedRNA_0240, predictedRNA_0325, and predictedRNA_0578. Future characterization and validation are necessary to fully elucidate the functions and roles of these sRNAs in MTB. Our study is the first to simultaneously unravel TSSs and sRNAs in MTB and demonstrate that the identification of other genomic features, such as TSSs, UTRs, and operons, allows for more accurate and specific classification of sRNAs., (© 2024. The Author(s), under exclusive licence to Springer-Verlag GmbH Germany, part of Springer Nature.)

The H3K4 methyltransferase SETD1A plays an essential role in both development and cancer. However, essential components involved in SETD1A chromatin binding remain unclear. Here, we discovered that BOD1L exhibits the highest correlated SETD1A co-dependency in human cancer cell lines. BOD1L knockout reduces leukemia cells in vitro and in vivo, and mimics the transcriptional profiles observed in SETD1A knockout cells. The loss of BOD1L immediately reduced SETD1A distribution at transcriptional start sites (TSS), induced transcriptional elongation defect, and increased the RNA polymerase II content at TSS; however, it did not reduce H3K4me3. The Shg1 domain of BOD1L has a DNA binding ability, and a tryptophan residue (W104) in the domain recruits SETD1A to chromatin through the association with SETD1A FLOS domain. In addition, the BOD1L-SETD1A complex associates with transcriptional regulators, including E2Fs. These results reveal that BOD1L mediates chromatin and SETD1A, and regulates the non-canonical function of SETD1A in transcription., (© The Author(s) 2024. Published by Oxford University Press on behalf of Nucleic Acids Research.)

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.

Syntenic long non-coding RNAs (lncRNAs) often show limited sequence conservation across species, prompting concern in the field. This study delves into functional signatures of syntenic lncRNAs between humans and zebrafish. Syntenic lncRNAs are highly expressed in zebrafish, with ∼90 % located near protein-coding genes, either in sense or antisense orientation. During early zebrafish development and in human embryonic stem cells (H1-hESC), syntenic lncRNA loci are enriched with cis-regulatory repressor signatures, influencing the expression of development-associated genes. In later zebrafish developmental stages and specific human cell lines, these syntenic lncRNA loci function as enhancers or transcription start sites (TSS) for protein-coding genes. Analysis of transposable elements (TEs) in syntenic lncRNA sequences revealed intriguing patterns: human lncRNAs are enriched in simple repeat elements, while their zebrafish counterparts show enrichment in LTR elements. This sequence evolution likely arises from post-rearrangement mutations that enhance DNA elements or cis-regulatory functions. It may also contribute to vertebrate innovation by creating novel transcription factor binding sites within the locus. This study highlights the conserved functionality of syntenic lncRNA loci through DNA elements, emphasizing their conserved roles across species despite sequence divergence., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2024 Elsevier B.V. All rights reserved.)

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.

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.

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.

We report that the mammalian 5-methylcytosine (5mC) oxidase Tet3 exists as three major isoforms and characterized the full-length isoform containing an N-terminal CXXC domain (Tet3FL). This CXXC domain binds to unmethylated CpGs, but, unexpectedly, its highest affinity is toward 5-carboxylcytosine (5caC). We determined the crystal structure of the CXXC domain-5caC-DNA complex, revealing the structural basis of the binding specificity of this domain as a reader of CcaCG sequences. Mapping of Tet3FL in neuronal cells shows that Tet3FL is localized precisely at the transcription start sites (TSSs) of genes involved in lysosome function, mRNA processing, and key genes of the base excision repair pathway. Therefore, Tet3FL may function as a regulator of 5caC removal by base excision repair. Active removal of accumulating 5mC from the TSSs of genes coding for lysosomal proteins by Tet3FL in postmitotic neurons of the brain may be important for preventing neurodegenerative diseases.

BackgroundTranscription initiation regulation is mediated by sequence-specific interactions between DNA-binding proteins (transcription factors) and cis-elements, where BRE, TATA, INR, DPE and MTE motifs constitute canonical core motifs for basal transcription initiation of genes. Accurate identification of transcription start site (TSS) and their corresponding promoter regions is critical for delineation of these motifs. To this end, the genome scale analysis of core promoter architecture in insects has been confined to Drosophila. The recently sequenced Tsetse fly genome provides a unique opportunity to analyze transcription initiation regulation machinery in blood-feeding insects.ResultsA computational method for identification of TSS in newly sequenced Tsetse fly genome was evaluated, using TSS seq tags sampled from two developmental stages namely; larvae and pupae. There were 3134 tag clusters among which 45.4% (1424) of the tag clusters mapped to first coding exons or their proximal predicted 5'UTR regions and 1.0% (31) tag clusters mapping to transposons, within a threshold of 100 tags per cluster. These 1393 non transposon-derived core promoters had propensity for AT nucleotides. The -1/+1 and 1/+1 positions in D. melanogaster, and G. m. morsitans had propensity for CA and AA dinucleotides respectively. The 1393 tag clusters comprised narrow promoters (5%), broad with peak promoters (23%) and broad without peak promoters (72%). Two-way motif co-occurrence analysis showed that the MTE-DPE pair is over-represented in broad core promoters. The frequently occurring triplet motifs in all promoter classes are the INR-MTE-DPE, TATA-MTE-DPE and TATA-INR-DPE. Promoters without the TATA motif had higher frequency of the MTE and INR motifs than those observed in Drosophila, where the DPE motif occur more frequently in promoters without TATA motif. Gene ontology terms associated with developmental processes were overrepresented in the narrow and broad with peak promoters.ConclusionsThe study has identified different motif combinations associated with broad promoters in a blood-feeding insect. In the case of TATA-less core promoters, G.m. morsitans uses the MTE to compensate for the lack of a TATA motif. The increasing availability of TSS seq data allows for revision of existing gene annotation datasets with the potential of identifying new transcriptional units.

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.)