TATA and paused promoters active in differentiated tissues have distinct expression characteristics

Single-cell nascent RNA sequencing unveils coordinated global transcription

Transcription is the primary regulatory step in gene expression. Divergent transcription initiation from promoters and enhancers produces stable RNAs from genes and unstable RNAs from enhancers1,2. Nascent RNA capture and sequencing assays simultaneously measure gene and enhancer activity in cell populations3. However, fundamental questions about the temporal regulation of transcription and enhancer-gene coordination remain unanswered, primarily because of the absence of a single-cell perspective on active transcription. In this study, we present scGRO-seq-a new single-cell nascent RNA sequencing assay that uses click chemistry-and unveil coordinated transcription throughout the genome. We demonstrate the episodic nature of transcription and the co-transcription of functionally related genes. scGRO-seq can estimate burst size and frequency by directly quantifying transcribing RNA polymerases in individual cells and can leverage replication-dependent non-polyadenylated histone gene transcription to elucidate cell cycle dynamics. The single-nucleotide spatial and temporal resolution of scGRO-seq enables the identification of networks of enhancers and genes. Our results suggest that the bursting of transcription at super-enhancers precedes bursting from associated genes. By imparting insights into the dynamic nature of global transcription and the origin and propagation of transcription signals, we demonstrate the ability of scGRO-seq to investigate the mechanisms of transcription regulation and the role of enhancers in gene expression.

Tissue-specific RNA Polymerase II promoter-proximal pause release and burst kinetics in a Drosophila embryonic patterning network

BACKGROUND

Formation of tissue-specific transcriptional programs underlies multicellular development, including dorsoventral (DV) patterning of the Drosophila embryo. This involves interactions between transcriptional enhancers and promoters in a chromatin context, but how the chromatin landscape influences transcription is not fully understood.

RESULTS

Here we comprehensively resolve differential transcriptional and chromatin states during Drosophila DV patterning. We find that RNA Polymerase II pausing is established at DV promoters prior to zygotic genome activation (ZGA), that pausing persists irrespective of cell fate, but that release into productive elongation is tightly regulated and accompanied by tissue-specific P-TEFb recruitment. DV enhancers acquire distinct tissue-specific chromatin states through CBP-mediated histone acetylation that predict the transcriptional output of target genes, whereas promoter states are more tissue-invariant. Transcriptome-wide inference of burst kinetics in different cell types revealed that while DV genes are generally characterized by a high burst size, either burst size or frequency can differ between tissues.

CONCLUSIONS

The data suggest that pausing is established by pioneer transcription factors prior to ZGA and that release from pausing is imparted by enhancer chromatin state to regulate bursting in a tissue-specific manner in the early embryo. Our results uncover how developmental patterning is orchestrated by tissue-specific bursts of transcription from Pol II primed promoters in response to enhancer regulatory cues.

Lola-I is a promoter pioneer factor that establishes de novo Pol II pausing during development

While the accessibility of enhancers is dynamically regulated during development, promoters tend to be constitutively accessible and poised for activation by paused Pol II. By studying Lola-I, a Drosophila zinc finger transcription factor, we show here that the promoter state can also be subject to developmental regulation independently of gene activation. Lola-I is ubiquitously expressed at the end of embryogenesis and causes its target promoters to become accessible and acquire paused Pol II throughout the embryo. This promoter transition is required but not sufficient for tissue-specific target gene activation. Lola-I mediates this function by depleting promoter nucleosomes, similar to the action of pioneer factors at enhancers. These results uncover a level of regulation for promoters that is normally found at enhancers and reveal a mechanism for the de novo establishment of paused Pol II at promoters.

A bioinformatics screen reveals hox and chromatin remodeling factors at the Drosophila histone locus

BACKGROUND

Cells orchestrate histone biogenesis with strict temporal and quantitative control. To efficiently regulate histone biogenesis, the repetitive Drosophila melanogaster replication-dependent histone genes are arrayed and clustered at a single locus. Regulatory factors concentrate in a nuclear body known as the histone locus body (HLB), which forms around the locus. Historically, HLB factors are largely discovered by chance, and few are known to interact directly with DNA. It is therefore unclear how the histone genes are specifically targeted for unique and coordinated regulation.

RESULTS

To expand the list of known HLB factors, we performed a candidate-based screen by mapping 30 publicly available ChIP datasets of 27 unique factors to the Drosophila histone gene array. We identified novel transcription factor candidates, including the Drosophila Hox proteins Ultrabithorax (Ubx), Abdominal-A (Abd-A), and Abdominal-B (Abd-B), suggesting a new pathway for these factors in influencing body plan morphogenesis. Additionally, we identified six other factors that target the histone gene array: JIL-1, hormone-like receptor 78 (Hr78), the long isoform of female sterile homeotic (1) (fs(1)h) as well as the general transcription factors TBP associated factor 1 (TAF-1), Transcription Factor IIB (TFIIB), and Transcription Factor IIF (TFIIF).

CONCLUSIONS

Our foundational screen provides several candidates for future studies into factors that may influence histone biogenesis. Further, our study emphasizes the powerful reservoir of publicly available datasets, which can be mined as a primary screening technique.

Single-cell nascent RNA sequencing using click-chemistry unveils coordinated transcription

Transcription is the primary regulatory step in gene expression. Divergent transcription initiation from promoters and enhancers produces stable RNAs from genes and unstable RNAs from enhancers1-5. Nascent RNA capture and sequencing assays simultaneously measure gene and enhancer activity in cell populations6-9. However, fundamental questions in the temporal regulation of transcription and enhancer-gene synchrony remain unanswered primarily due to the absence of a single-cell perspective on active transcription. In this study, we present scGRO-seq - a novel single-cell nascent RNA sequencing assay using click-chemistry - and unveil the coordinated transcription throughout the genome. scGRO-seq demonstrates the episodic nature of transcription, and estimates burst size and frequency by directly quantifying transcribing RNA polymerases in individual cells. It reveals the co-transcription of functionally related genes and leverages the replication-dependent non-polyadenylated histone genes transcription to elucidate cell-cycle dynamics. The single-nucleotide spatial and temporal resolution of scGRO-seq identifies networks of enhancers and genes and indicates that the bursting of transcription at super-enhancers precedes the burst from associated genes. By imparting insights into the dynamic nature of transcription and the origin and propagation of transcription signals, scGRO-seq demonstrates its unique ability to investigate the mechanisms of transcription regulation and the role of enhancers in gene expression.

Genetic variation in chromatin state across multiple tissues in Drosophila melanogaster

We use ATAC-seq to examine chromatin accessibility for four different tissues in Drosophila melanogaster: adult female brain, ovaries, and both wing and eye-antennal imaginal discs from males. Each tissue is assayed in eight different inbred strain genetic backgrounds, seven associated with a reference quality genome assembly. We develop a method for the quantile normalization of ATAC-seq fragments and test for differences in coverage among genotypes, tissues, and their interaction at 44099 peaks throughout the euchromatic genome. For the strains with reference quality genome assemblies, we correct ATAC-seq profiles for read mis-mapping due to nearby polymorphic structural variants (SVs). Comparing coverage among genotypes without accounting for SVs results in a highly elevated rate (55%) of identifying false positive differences in chromatin state between genotypes. After SV correction, we identify 1050, 30383, and 4508 regions whose peak heights are polymorphic among genotypes, among tissues, or exhibit genotype-by-tissue interactions, respectively. Finally, we identify 3988 candidate causative variants that explain at least 80% of the variance in chromatin state at nearby ATAC-seq peaks.

CBP and Gcn5 drive zygotic genome activation independently of their catalytic activity

Zygotic genome activation (ZGA) is a crucial step of embryonic development. So far, little is known about the role of chromatin factors during this process. Here, we used an in vivo RNA interference reverse genetic screen to identify chromatin factors necessary for embryonic development in Drosophila melanogaster. Our screen reveals that histone acetyltransferases (HATs) and histone deacetylases are crucial ZGA regulators. We demonstrate that Nejire (CBP/EP300 ortholog) is essential for the acetylation of histone H3 lysine-18 and lysine-27, whereas Gcn5 (GCN5/PCAF ortholog) for lysine-9 of H3 at ZGA, with these marks being enriched at all actively transcribed genes. Nonetheless, these HATs activate distinct sets of genes. Unexpectedly, individual catalytic dead mutants of either Nejire or Gcn5 can activate zygotic transcription (ZGA) and transactivate a reporter gene in vitro. Together, our data identify Nejire and Gcn5 as key regulators of ZGA.

Genome-wide inference reveals that feedback regulations constrain promoter-dependent transcriptional burst kinetics

Gene expression in mammalian cells is highly variable and episodic, resulting in a series of discontinuous bursts of mRNAs. A challenge is to understand how static promoter architecture and dynamic feedback regulations dictate bursting on a genome-wide scale. Although single-cell RNA sequencing (scRNA-seq) provides an opportunity to address this challenge, effective analytical methods are scarce. We developed an interpretable and scalable inference framework, which combined experimental data with a mechanistic model to infer transcriptional burst kinetics (sizes and frequencies) and feedback regulations. Applying this framework to scRNA-seq data generated from embryonic mouse fibroblast cells, we found Simpson's paradoxes, i.e. genome-wide burst kinetics exhibit different characteristics in two cases without and with distinguishing feedback regulations. We also showed that feedbacks differently modulate burst frequencies and sizes and conceal the effects of transcription start site distributions on burst kinetics. Notably, only in the presence of positive feedback, TATA genes are expressed with high burst frequencies and enhancer-promoter interactions mainly modulate burst frequencies. The developed inference method provided a flexible and efficient way to investigate transcriptional burst kinetics and the obtained results would be helpful for understanding cell development and fate decision.

Genome-Wide Identification and Characterization of CPR5 Genes in Gossypium Reveals Their Potential Role in Trichome Development

Trichomes protect plants against insects, microbes, herbivores, and abiotic damages and assist seed dispersal. The function of CPR5 genes have been found to be involved in the trichome development but the research on the underlying genetic and molecular mechanisms are extremely limited. Herein, genome wide identification and characterization of CPR5 genes was performed. In total, 26 CPR5 family members were identified in Gossypium species. Phylogenetic analysis, structural characteristics, and synteny analysis of CPR5s showed the conserved evolution relationships of CPR5. The promoter analysis of CPR5 genes revealed hormone, stress, and development-related cis-elements. Gene ontology (GO) enrichment analysis showed that the CPR5 genes were largely related to biological regulation, developmental process, multicellular organismal process. Protein-protein interaction analysis predicted several trichome development related proteins (SIM, LGO, and GRL) directly interacting with CPR5 genes. Further, nine putative Gossypium-miRNAs were also identified, targeting Gossypium CPR5 genes. RNA-Seq data of G. arboreum (with trichomes) and G. herbaceum (with no trichomes) was used to perform the co-expression network analysis. GheCPR5.1 was identified as a hub gene in a co-expression network analysis. RT-qPCR of GheCPR5.1 gene in different tissues suggests that this gene has higher expressions in the petiole and might be a key candidate involved in the trichome development. Virus induced gene silencing of GheCPR5.1 (Ghe02G17590) confirms its role in trichome development and elongation. Current results provide proofs of the possible role of CPR5 genes and provide preliminary information for further studies of GheCPR5.1 functions in trichome development.

Identification and Somatic Characterization of the Germline PTEN Promoter Variant rs34149102 in a Family with Gastrointestinal and Breast Tumors

Genetic variants located in non-coding regions can affect processes that regulate protein expression, functionally contributing to human disease. Germline heterozygous mutations in the non-coding region of the PTEN gene have been previously identified in patients with PTEN hamartoma tumor syndrome (PHTS) diagnosed with breast, thyroid, and/or endometrial cancer. In this study, we report a PTEN promoter variant (rs34149102 A allele) that was identified by direct sequencing in an Italian family with a history of gastroesophageal junction (GEJ) adenocarcinoma and breast cancer. In order to investigate the putative functional role of the rs34149102 A allele variant, we evaluated the status of PTEN alterations at the somatic level. We found that PTEN protein expression was absent in the GEJ adenocarcinoma tissue of the index case. Moreover, we detected the occurrence of copy number loss involving the PTEN rs34149102 major C allele in tumor tissue, revealing that the second allele was somatically inactivated. This variant is located within an active regulatory region of the PTEN core promoter, and in silico analysis suggests that it may affect the binding of the nuclear transcription factor MAZ and hence PTEN expression. Overall, these results reveal the functional role of the PTEN promoter rs34149102 A allele variant in the modulation of PTEN protein expression and highlight its contribution to hereditary cancer risk.

Two promoters integrate multiple enhancer inputs to drive wild-type knirps expression in the Drosophila melanogaster embryo

Proper development depends on precise spatiotemporal gene expression patterns. Most developmental genes are regulated by multiple enhancers and often by multiple core promoters that generate similar transcripts. We hypothesize that multiple promoters may be required either because enhancers prefer a specific promoter or because multiple promoters serve as a redundancy mechanism. To test these hypotheses, we studied the expression of the knirps locus in the early Drosophila melanogaster embryo, which is mediated by multiple enhancers and core promoters. We found that one of these promoters resembles a typical "sharp" developmental promoter, while the other resembles a "broad" promoter usually associated with housekeeping genes. Using synthetic reporter constructs, we found that some, but not all, enhancers in the locus show a preference for one promoter, indicating that promoters provide both redundancy and specificity. By analyzing the reporter dynamics, we identified specific burst properties during the transcription process, namely burst size and frequency, that are most strongly tuned by the combination of promoter and enhancer. Using locus-sized reporters, we discovered that enhancers with no promoter preference in a synthetic setting have a preference in the locus context. Our results suggest that the presence of multiple promoters in a locus is due both to enhancer preference and a need for redundancy and that "broad" promoters with dispersed transcription start sites are common among developmental genes. They also imply that it can be difficult to extrapolate expression measurements from synthetic reporters to the locus context, where other variables shape a gene's overall expression pattern.

Histone variant H2A.Z regulates zygotic genome activation

During embryogenesis, the genome shifts from transcriptionally quiescent to extensively active in a process known as Zygotic Genome Activation (ZGA). In Drosophila, the pioneer factor Zelda is known to be essential for the progression of development; still, it regulates the activation of only a small subset of genes at ZGA. However, thousands of genes do not require Zelda, suggesting that other mechanisms exist. By conducting GRO-seq, HiC and ChIP-seq in Drosophila embryos, we demonstrate that up to 65% of zygotically activated genes are enriched for the histone variant H2A.Z. H2A.Z enrichment precedes ZGA and RNA Polymerase II loading onto chromatin. In vivo knockdown of maternally contributed Domino, a histone chaperone and ATPase, reduces H2A.Z deposition at transcription start sites, causes global downregulation of housekeeping genes at ZGA, and compromises the establishment of the 3D chromatin structure. We infer that H2A.Z is essential for the de novo establishment of transcriptional programs during ZGA via chromatin reorganization.

During embryogenesis, the genome becomes transcriptionally active in a process known as zygotic genome activation (ZGA); how ZGA is initiated is still an open question. Here the authors show histone variant H2A.Z deposition precedes RNA polymerase II binding on chromatin, before ZGA. H2A.Z loss causes transcriptional downregulation of ZGA genes and leads to changes in the 3D genome organization.

Engineering a genome-reduced bacterium to eliminate Staphylococcus aureus biofilms in vivo

Bacteria present a promising delivery system for treating human diseases. Here, we engineered the genome-reduced human lung pathogen Mycoplasma pneumoniae as a live biotherapeutic to treat biofilm-associated bacterial infections. This strain has a unique genetic code, which hinders gene transfer to most other bacterial genera, and it lacks a cell wall, which allows it to express proteins that target peptidoglycans of pathogenic bacteria. We first determined that removal of the pathogenic factors fully attenuated the chassis strain in vivo. We then designed synthetic promoters and identified an endogenous peptide signal sequence that, when fused to heterologous proteins, promotes efficient secretion. Based on this, we equipped the chassis strain with a genetic platform designed to secrete antibiofilm and bactericidal enzymes, resulting in a strain capable of dissolving Staphylococcus aureus biofilms preformed on catheters in vitro, ex vivo, and in vivo. To our knowledge, this is the first engineered genome-reduced bacterium that can fight against clinically relevant biofilm-associated bacterial infections.

Single-cell imaging of genome organization and dynamics

Probing the architecture, mechanism, and dynamics of genome folding is fundamental to our understanding of genome function in homeostasis and disease. Most chromosome conformation capture studies dissect the genome architecture with population- and time-averaged snapshots and thus have limited capabilities to reveal 3D nuclear organization and dynamics at the single-cell level. Here, we discuss emerging imaging techniques ranging from light microscopy to electron microscopy that enable investigation of genome folding and dynamics at high spatial and temporal resolution. Results from these studies complement genomic data, unveiling principles underlying the spatial arrangement of the genome and its potential functional links to diverse biological activities in the nucleus.

Sequence analysis of tyrosine recombinases allows annotation of mobile genetic elements in prokaryotic genomes

Mobile genetic elements (MGEs) sequester and mobilize antibiotic resistance genes across bacterial genomes. Efficient and reliable identification of such elements is necessary to follow resistance spreading. However, automated tools for MGE identification are missing. Tyrosine recombinase (YR) proteins drive MGE mobilization and could provide markers for MGE detection, but they constitute a diverse family also involved in housekeeping functions. Here, we conducted a comprehensive survey of YRs from bacterial, archaeal, and phage genomes and developed a sequence-based classification system that dissects the characteristics of MGE-borne YRs. We revealed that MGE-related YRs evolved from non-mobile YRs by acquisition of a regulatory arm-binding domain that is essential for their mobility function. Based on these results, we further identified numerous unknown MGEs. This work provides a resource for comparative analysis and functional annotation of YRs and aids the development of computational tools for MGE annotation. Additionally, we reveal how YRs adapted to drive gene transfer across species and provide a tool to better characterize antibiotic resistance dissemination.