High-resolution, genome-wide mapping of positive supercoiling in chromosomes

Depletion of CpG dinucleotides in bacterial genomes may represent an adaptation to high temperatures

Dinucleotide biases have been widely investigated in the genomes of eukaryotes and viruses, but not in bacteria. We assembled a dataset of bacterial genomes (>15 000), which are representative of the genetic diversity in the kingdom Eubacteria, and we analyzed dinucleotide biases in relation to different traits. We found that TpA dinucleotides are the most depleted and that CpG dinucleotides show the widest dispersion. The abundances of both dinucleotides vary with genomic G + C content and show a very strong phylogenetic signal. After accounting for G + C content and phylogenetic inertia, we analyzed different bacterial lifestyle traits. We found that temperature preferences associate with the abundance of CpG dinucleotides, with thermophiles/hyperthemophiles being particularly depleted. Conversely, the TpA dinucleotide displays a bias that only depends on genomic G + C composition. Using predictions of intrinsic cyclizability we also show that CpG depletion may associate with higher DNA bendability in both thermophiles/hyperthermophiles and mesophiles, and that the former are predicted to have significantly more flexible genomes than the latter. We suggest that higher bendability is advantageous at high temperatures because it facilitates DNA positive supercoiling and that, through modulation of DNA mechanical properties, local or global CpG depletion controls genome organization, most likely not only in bacteria.

Nucleosomes play a dual role in regulating transcription dynamics

Transcription has a mechanical component, as the translocation of the transcription machinery or RNA polymerase (RNAP) on DNA or chromatin is dynamically coupled to the chromatin torsion. This posits chromatin mechanics as a possible regulator of eukaryotic transcription, however, the modes and mechanisms of this regulation are elusive. Here, we first take a statistical mechanics approach to model the torsional response of topology-constrained chromatin. Our model recapitulates the experimentally observed weaker torsional stiffness of chromatin compared to bare DNA and proposes structural transitions of nucleosomes into chirally distinct states as the driver of the contrasting torsional mechanics. Coupling chromatin mechanics with RNAP translocation in stochastic simulations, we reveal a complex interplay of DNA supercoiling and nucleosome dynamics in governing RNAP velocity. Nucleosomes play a dual role in controlling the transcription dynamics. The steric barrier aspect of nucleosomes in the gene body counteracts transcription via hindering RNAP motion, whereas the chiral transitions facilitate RNAP motion via driving a low restoring torque upon twisting the DNA. While nucleosomes with low dissociation rates are typically transcriptionally repressive, highly dynamic nucleosomes offer less of a steric barrier and enhance the transcription elongation dynamics of weakly transcribed genes via buffering DNA twist. We use the model to predict transcription-dependent levels of DNA supercoiling in segments of the budding yeast genome that are in accord with available experimental data. The model unveils a paradigm of DNA supercoiling-mediated interaction between genes and makes testable predictions that will guide experimental design.

Structure of the E. coli nucleoid-associated protein YejK reveals a novel DNA binding clamp

Nucleoid-associated proteins (NAPs) play central roles in bacterial chromosome organization and DNA processes. The Escherichia coli YejK protein is a highly abundant, yet poorly understood NAP. YejK proteins are conserved among Gram-negative bacteria but show no homology to any previously characterized DNA-binding protein. Hence, how YejK binds DNA is unknown. To gain insight into YejK structure and its DNA binding mechanism we performed biochemical and structural analyses on the E. coli YejK protein. Biochemical assays demonstrate that, unlike many NAPs, YejK does not show a preference for AT-rich DNA and binds non-sequence specifically. A crystal structure revealed YejK adopts a novel fold comprised of two domains. Strikingly, each of the domains harbors an extended arm that mediates dimerization, creating an asymmetric clamp with a 30 Å diameter pore. The lining of the pore is electropositive and mutagenesis combined with fluorescence polarization assays support DNA binding within the pore. Finally, our biochemical analyses on truncated YejK proteins suggest a mechanism for YejK clamp loading. Thus, these data reveal YejK contains a newly described DNA-binding motif that functions as a novel clamp.

Bacterial chromatin proteins, transcription, and DNA topology: Inseparable partners in the control of gene expression

DNA in bacterial chromosomes is organized into higher-order structures by DNA-binding proteins called nucleoid-associated proteins (NAPs) or bacterial chromatin proteins (BCPs). BCPs often bind to or near DNA loci transcribed by RNA polymerase (RNAP) and can either increase or decrease gene expression. To understand the mechanisms by which BCPs alter transcription, one must consider both steric effects and the topological forces that arise when DNA deviates from its fully relaxed double-helical structure. Transcribing RNAP creates DNA negative (-) supercoils upstream and positive (+) supercoils downstream whenever RNAP and DNA are unable to rotate freely. This (-) and (+) supercoiling generates topological forces that resist forward translocation of DNA through RNAP unless the supercoiling is constrained by BCPs or relieved by topoisomerases. BCPs also may enhance topological stress and overall can either inhibit or aid transcription. Here, we review current understanding of how RNAP, BCPs, and DNA topology interplay to control gene expression.

Regulators of rDNA array morphology in fission yeast

Nucleolar morphology is a well-established indicator of ribosome biogenesis activity that has served as the foundation of many screens investigating ribosome production. Missing from this field of study is a broad-scale investigation of the regulation of ribosomal DNA morphology, despite the essential role of rRNA gene transcription in modulating ribosome output. We hypothesized that the morphology of rDNA arrays reflects ribosome biogenesis activity. We established GapR-GFP, a prokaryotic DNA-binding protein that recognizes transcriptionally-induced overtwisted DNA, as a live visual fluorescent marker for quantitative analysis of rDNA organization in Schizosaccharomyces pombe. We found that the morphology-which we refer to as spatial organization-of the rDNA arrays is dynamic throughout the cell cycle, under glucose starvation, RNA pol I inhibition, and TOR activation. Screening the haploid S. pombe Bioneer deletion collection for spatial organization phenotypes revealed large ribosomal protein (RPL) gene deletions that alter rDNA organization. Further work revealed RPL gene deletion mutants with altered rDNA organization also demonstrate resistance to the TOR inhibitor Torin1. A genetic analysis of signaling pathways essential for this resistance phenotype implicated many factors including a conserved MAPK, Pmk1, previously linked to extracellular stress responses. We propose RPL gene deletion triggers altered rDNA morphology due to compensatory changes in ribosome biogenesis via multiple signaling pathways, and we further suggest compensatory responses may contribute to human diseases such as ribosomopathies. Altogether, GapR-GFP is a powerful tool for live visual reporting on rDNA morphology under myriad conditions.

Spatio-temporal organization of the E. coli chromosome from base to cellular length scales

Escherichia coli has been a vital model organism for studying chromosomal structure, thanks, in part, to its small and circular genome (4.6 million base pairs) and well-characterized biochemical pathways. Over the last several decades, we have made considerable progress in understanding the intricacies of the structure and subsequent function of the E. coli nucleoid. At the smallest scale, DNA, with no physical constraints, takes on a shape reminiscent of a randomly twisted cable, forming mostly random coils but partly affected by its stiffness. This ball-of-spaghetti-like shape forms a structure several times too large to fit into the cell. Once the physiological constraints of the cell are added, the DNA takes on overtwisted (negatively supercoiled) structures, which are shaped by an intricate interplay of many proteins carrying out essential biological processes. At shorter length scales (up to about 1 kb), nucleoid-associated proteins organize and condense the chromosome by inducing loops, bends, and forming bridges. Zooming out further and including cellular processes, topological domains are formed, which are flanked by supercoiling barriers. At the megabase-scale both large, highly self-interacting regions (macrodomains) and strong contacts between distant but co-regulated genes have been observed. At the largest scale, the nucleoid forms a helical ellipsoid. In this review, we will explore the history and recent advances that pave the way for a better understanding of E. coli chromosome organization and structure, discussing the cellular processes that drive changes in DNA shape, and what contributes to compaction and formation of dynamic structures, and in turn how bacterial chromatin affects key processes such as transcription and replication.

Physical models of bacterial chromosomes

The interplay between bacterial chromosome organization and functions such as transcription and replication can be studied in increasing detail using novel experimental techniques. Interpreting the resulting quantitative data, however, can be theoretically challenging. In this minireview, we discuss how connecting experimental observations to biophysical theory and modeling can give rise to new insights on bacterial chromosome organization. We consider three flavors of models of increasing complexity: simple polymer models that explore how physical constraints, such as confinement or plectoneme branching, can affect bacterial chromosome organization; bottom-up mechanistic models that connect these constraints to their underlying causes, for instance, chromosome compaction to macromolecular crowding, or supercoiling to transcription; and finally, data-driven methods for inferring interpretable and quantitative models directly from complex experimental data. Using recent examples, we discuss how biophysical models can both deepen our understanding of how bacterial chromosomes are structured and give rise to novel predictions about bacterial chromosome organization.

Enricherator: A Bayesian Method for Inferring Regularized Genome-wide Enrichments from Sequencing Count Data

A pervasive question in biological research studying gene regulation, chromatin structure, or genomics is where, and to what extent, does a signal of interest arise genome-wide? This question is addressed using a variety of methods relying on high-throughput sequencing data as their final output, including ChIP-seq for protein-DNA interactions,1 GapR-seq for measuring supercoiling,2 and HBD-seq or DRIP-seq for R-loop positioning.3,4 Current computational methods to calculate genome-wide enrichment of the signal of interest usually do not properly handle the count-based nature of sequencing data, they often do not make use of the local correlation structure of sequencing data, and they do not apply any regularization of enrichment estimates. This can result in unrealistic estimates of the true underlying biological enrichment of interest, unrealistically low estimates of confidence in point estimates of enrichment (or no estimates of confidence at all), unrealistic gyrations in enrichment estimates at very close (<10 bp) genomic loci due to noise inherent in sequencing data, and in a multiple-hypothesis testing problem during interpretation of genome-wide enrichment estimates. We developed a tool called Enricherator to infer genome-wide enrichments from sequencing count data. Enricherator uses the variational Bayes algorithm to fit a generalized linear model to sequencing count data and to sample from the approximate posterior distribution of enrichment estimates (https://github.com/jwschroeder3/enricherator). Enrichments inferred by Enricherator more precisely identify known binding sites in cases where low coverage between binding sites leads to false-positive peak calls in these noisy regions of the genome; these benefits extend to published datasets.

Loop-extruding Smc5/6 organizes transcription-induced positive DNA supercoils

The structural maintenance of chromosomes (SMC) protein complexes-cohesin, condensin, and the Smc5/6 complex (Smc5/6)-are essential for chromosome function. At the molecular level, these complexes fold DNA by loop extrusion. Accordingly, cohesin creates chromosome loops in interphase, and condensin compacts mitotic chromosomes. However, the role of Smc5/6's recently discovered DNA loop extrusion activity is unknown. Here, we uncover that Smc5/6 associates with transcription-induced positively supercoiled DNA at cohesin-dependent loop boundaries on budding yeast (Saccharomyces cerevisiae) chromosomes. Mechanistically, single-molecule imaging reveals that dimers of Smc5/6 specifically recognize the tip of positively supercoiled DNA plectonemes and efficiently initiate loop extrusion to gather the supercoiled DNA into a large plectonemic loop. Finally, Hi-C analysis shows that Smc5/6 links chromosomal regions containing transcription-induced positive supercoiling in cis. Altogether, our findings indicate that Smc5/6 controls the three-dimensional organization of chromosomes by recognizing and initiating loop extrusion on positively supercoiled DNA.

Transcription-induced domains form the elementary constraining building blocks of bacterial chromosomes

Transcription generates local topological and mechanical constraints on the DNA fiber, leading to the generation of supercoiled chromosome domains in bacteria. However, the global impact of transcription on chromosome organization remains elusive, as the scale of genes and operons in bacteria remains well below the resolution of chromosomal contact maps generated using Hi-C (~5-10 kb). Here we combined sub-kb Hi-C contact maps and chromosome engineering to visualize individual transcriptional units. We show that transcriptional units form discrete three-dimensional transcription-induced domains that impose mechanical and topological constraints on their neighboring sequences at larger scales, modifying their localization and dynamics. These results show that transcriptional domains constitute primary building blocks of bacterial chromosome folding and locally impose structural and dynamic constraints.

Differential roles of positive and negative supercoiling in organizing the E. coli genome

This study aims to explore whether and how positive and negative supercoiling contribute to the three-dimensional (3D) organization of the bacterial genome. We used recently published Escherichia coli GapR ChIP-seq and TopoI ChIP-seq (also called EcTopoI-seq) data, which marks positive and negative supercoiling sites, respectively, to study how supercoiling correlates with the spatial contact maps obtained from chromosome conformation capture sequencing (Hi-C and 5C). We find that supercoiled chromosomal loci have overall higher Hi-C contact frequencies than sites that are not supercoiled. Surprisingly, positive supercoiling corresponds to higher spatial contact than negative supercoiling. Additionally, positive, but not negative, supercoiling could be identified from Hi-C data with high accuracy. We further find that the majority of positive and negative supercoils coincide with highly active transcription units, with a minor group likely associated with replication and other genomic processes. Our results show that both positive and negative supercoiling enhance spatial contact, with positive supercoiling playing a larger role in bringing genomic loci closer in space. Based on our results, we propose new physical models of how the E. coli chromosome is organized by positive and negative supercoils.

The nucleolar shell provides anchoring sites for DNA untwisting

DNA underwinding (untwisting) is a crucial step in transcriptional activation. DNA underwinding occurs between the site where torque is generated by RNA polymerase (RNAP) and the site where the axial rotation of DNA is constrained. However, what constrains DNA axial rotation in the nucleus is yet unknown. Here, we show that the anchorage to the nuclear protein condensates constrains DNA axial rotation for DNA underwinding in the nucleolus. In situ super-resolution imaging of underwound DNA reveal that underwound DNA accumulates in the nucleolus, a nuclear condensate with a core-shell structure. Specifically, underwound DNA is distributed in the nucleolar core owing to RNA polymerase I (RNAPI) activities. Furthermore, underwound DNA in the core decreases when nucleolar shell components are prevented from binding to their recognition structure, G-quadruplex (G4). Taken together, these results suggest that the nucleolar shell provides anchoring sites that constrain DNA axial rotation for RNAPI-driven DNA underwinding in the core. Our findings will contribute to understanding how nuclear protein condensates make up constraints for the site-specific regulation of DNA underwinding and transcription.

How Do Thermophiles Organize Their Genomes?

All cells must maintain the structural and functional integrity of the genome under a wide range of environments. High temperatures pose a formidable challenge to cells by denaturing the DNA double helix, causing chemical damage to DNA, and increasing the random thermal motion of chromosomes. Thermophiles, predominantly classified as bacteria or archaea, exhibit an exceptional capacity to mitigate these detrimental effects and prosper under extreme thermal conditions, with some species tolerating temperatures higher than 100°C. Their genomes are mainly characterized by the presence of reverse gyrase, a unique topoisomerase that introduces positive supercoils into DNA. This enzyme has been suggested to maintain the genome integrity of thermophiles by limiting DNA melting and mediating DNA repair. Previous studies provided significant insights into the mechanisms by which NAPs, histones, SMC superfamily proteins, and polyamines affect the 3D genomes of thermophiles across different scales. Here, I discuss current knowledge of the genome organization in thermophiles and pertinent research questions for future investigations.

ERK2-topoisomerase II regulatory axis is important for gene activation in immediate early genes

The function of the mitogen-activated protein kinase signaling pathway is required for the activation of immediate early genes (IEGs), including EGR1 and FOS, for cell growth and proliferation. Recent studies have identified topoisomerase II (TOP2) as one of the important regulators of the transcriptional activation of IEGs. However, the mechanism underlying transcriptional regulation involving TOP2 in IEG activation has remained unknown. Here, we demonstrate that ERK2, but not ERK1, is important for IEG transcriptional activation and report a critical ELK1 binding sequence for ERK2 function at the EGR1 gene. Our data indicate that both ERK1 and ERK2 extensively phosphorylate the C-terminal domain of TOP2B at mutual and distinctive residues. Although both ERK1 and ERK2 enhance the catalytic rate of TOP2B required to relax positive DNA supercoiling, ERK2 delays TOP2B catalysis of negative DNA supercoiling. In addition, ERK1 may relax DNA supercoiling by itself. ERK2 catalytic inhibition or knock-down interferes with transcription and deregulates TOP2B in IEGs. Furthermore, we present the first cryo-EM structure of the human cell-purified TOP2B and etoposide together with the EGR1 transcriptional start site (-30 to +20) that has the strongest affinity to TOP2B within -423 to +332. The structure shows TOP2B-mediated breakage and dramatic bending of the DNA. Transcription is activated by etoposide, while it is inhibited by ICRF193 at EGR1 and FOS, suggesting that TOP2B-mediated DNA break to favor transcriptional activation. Taken together, this study suggests that activated ERK2 phosphorylates TOP2B to regulate TOP2-DNA interactions and favor transcriptional activation in IEGs. We propose that TOP2B association, catalysis, and dissociation on its substrate DNA are important processes for regulating transcription and that ERK2-mediated TOP2B phosphorylation may be key for the catalysis and dissociation steps.

The abscission checkpoint senses chromatin bridges through Top2α recruitment to DNA knots

In response to chromatin bridges, the abscission checkpoint delays completion of cytokinesis to prevent chromosome breakage or tetraploidization. Here, we show that spontaneous or replication stress-induced chromatin bridges exhibit "knots" of catenated and overtwisted DNA next to the midbody. Topoisomerase IIα (Top2α) forms abortive Top2-DNA cleavage complexes (Top2ccs) on DNA knots; furthermore, impaired Top2α-DNA cleavage activity correlates with chromatin bridge breakage in cytokinesis. Proteasomal degradation of Top2ccs is required for Rad17 localization to Top2-generated double-strand DNA ends on DNA knots; in turn, Rad17 promotes local recruitment of the MRN complex and downstream ATM-Chk2-INCENP signaling to delay abscission and prevent chromatin breakage. In contrast, dicentric chromosomes that do not exhibit knotted DNA fail to activate the abscission checkpoint in human cells. These findings are the first to describe a mechanism by which the abscission checkpoint detects chromatin bridges, through generation of abortive Top2ccs on DNA knots, to preserve genome integrity.

Insights in bacterial genome folding

Chromosomes in all domains of life are well-defined structural entities with complex hierarchical organization. The regulation of this hierarchical organization and its functional interplay with gene expression or other chromosome metabolic processes such as repair, replication, or segregation is actively investigated in a variety of species, including prokaryotes. Bacterial chromosomes are typically gene-dense with few non-coding sequences and are organized into the nucleoid, a membrane-less compartment composed of DNA, RNA, and proteins (nucleoid-associated proteins or NAPs). The continuous improvement of imaging and genomic methods has put the organization of these Mb-long molecules at reach, allowing to disambiguate some of their highly dynamic properties and intertwined structural features. Here we review and discuss some of the recent advances in the field of bacterial chromosome organization.

Direct observation of coordinated assembly of individual native centromeric nucleosomes

Eukaryotic chromosome segregation requires the kinetochore, a megadalton-sized machine that forms on specialized centromeric chromatin containing CENP-A, a histone H3 variant. CENP-A deposition requires a chaperone protein HJURP that targets it to the centromere, but it has remained unclear whether HJURP has additional functions beyond CENP-A targeting and why high AT DNA content, which disfavors nucleosome assembly, is widely conserved at centromeres. To overcome the difficulties of studying nucleosome formation in vivo, we developed a microscopy assay that enables direct observation of de novo centromeric nucleosome recruitment and maintenance with single molecule resolution. Using this assay, we discover that CENP-A can arrive at centromeres without its dedicated centromere-specific chaperone HJURP, but stable incorporation depends on HJURP and additional DNA-binding proteins of the inner kinetochore. We also show that homopolymer AT runs in the yeast centromeres are essential for efficient CENP-A deposition. Together, our findings reveal requirements for stable nucleosome formation and provide a foundation for further studies of the assembly and dynamics of native kinetochore complexes.

Dynamics of chromosome organization in a minimal bacterial cell

Computational models of cells cannot be considered complete unless they include the most fundamental process of life, the replication and inheritance of genetic material. By creating a computational framework to model systems of replicating bacterial chromosomes as polymers at 10 bp resolution with Brownian dynamics, we investigate changes in chromosome organization during replication and extend the applicability of an existing whole-cell model (WCM) for a genetically minimal bacterium, JCVI-syn3A, to the entire cell-cycle. To achieve cell-scale chromosome structures that are realistic, we model the chromosome as a self-avoiding homopolymer with bending and torsional stiffnesses that capture the essential mechanical properties of dsDNA in Syn3A. In addition, the conformations of the circular DNA must avoid overlapping with ribosomes identitied in cryo-electron tomograms. While Syn3A lacks the complex regulatory systems known to orchestrate chromosome segregation in other bacteria, its minimized genome retains essential loop-extruding structural maintenance of chromosomes (SMC) protein complexes (SMC-scpAB) and topoisomerases. Through implementing the effects of these proteins in our simulations of replicating chromosomes, we find that they alone are sufficient for simultaneous chromosome segregation across all generations within nested theta structures. This supports previous studies suggesting loop-extrusion serves as a near-universal mechanism for chromosome organization within bacterial and eukaryotic cells. Furthermore, we analyze ribosome diffusion under the influence of the chromosome and calculate in silico chromosome contact maps that capture inter-daughter interactions. Finally, we present a methodology to map the polymer model of the chromosome to a Martini coarse-grained representation to prepare molecular dynamics models of entire Syn3A cells, which serves as an ultimate means of validation for cell states predicted by the WCM.

In Vivo Sampling of Intracellular Heterogeneity of Pseudomonas putida Enables Multiobjective Optimization of Genetic Devices

The inner physicochemical heterogeneity of bacterial cells generates three-dimensional (3D)-dependent variations of resources for effective expression of given chromosomally located genes. This fact has been exploited for adjusting the most favorable parameters for implanting a complex device for optogenetic control of biofilm formation in the soil bacterium Pseudomonas putida. To this end, a DNA segment encoding a superactive variant of the Caulobacter crescendus diguanylate cyclase PleD expressed under the control of the cyanobacterial light-responsive CcaSR system was placed in a mini-Tn5 transposon vector and randomly inserted through the chromosome of wild-type and biofilm-deficient variants of P. putida lacking the wsp gene cluster. This operation delivered a collection of clones covering a whole range of biofilm-building capacities and dynamic ranges in response to green light. Since the phenotypic output of the device depends on a large number of parameters (multiple promoters, RNA stability, translational efficacy, metabolic precursors, protein folding, etc.), we argue that random chromosomal insertions enable sampling the intracellular milieu for an optimal set of resources that deliver a preset phenotypic specification. Results thus support the notion that the context dependency can be exploited as a tool for multiobjective optimization, rather than a foe to be suppressed in Synthetic Biology constructs.

DNA supercoiling restricts the transcriptional bursting of neighboring eukaryotic genes

DNA supercoiling has emerged as a major contributor to gene regulation in bacteria, but how DNA supercoiling impacts transcription dynamics in eukaryotes is unclear. Here, using single-molecule dual-color nascent transcription imaging in budding yeast, we show that transcriptional bursting of divergent and tandem GAL genes is coupled. Temporal coupling of neighboring genes requires rapid release of DNA supercoils by topoisomerases. When DNA supercoils accumulate, transcription of one gene inhibits transcription at its adjacent genes. Transcription inhibition of the GAL genes results from destabilized binding of the transcription factor Gal4. Moreover, wild-type yeast minimizes supercoiling-mediated inhibition by maintaining sufficient levels of topoisomerases. Overall, we discover fundamental differences in transcriptional control by DNA supercoiling between bacteria and yeast and show that rapid supercoiling release in eukaryotes ensures proper gene expression of neighboring genes.

Direct observation of coordinated assembly of individual native centromeric nucleosomes

Eukaryotic chromosome segregation requires the kinetochore, a megadalton-sized machine that forms on specialized centromeric chromatin containing CENP-A, a histone H3 variant. CENP-A deposition requires a chaperone protein HJURP that targets it to the centromere, but it has remained unclear whether HJURP has additional functions beyond CENP-A targeting and why high AT DNA content, which disfavors nucleosome assembly, is widely conserved at centromeres. To overcome the difficulties of studying nucleosome formation in vivo, we developed a microscopy assay that enables direct observation of de novo centromeric nucleosome recruitment and maintenance with single molecule resolution. Using this assay, we discover that CENP-A can arrive at centromeres without its dedicated centromere-specific chaperone HJURP, but stable incorporation depends on HJURP and additional DNA-binding proteins of the inner kinetochore. We also show that homopolymer AT runs in the yeast centromeres are essential for efficient CENP-A deposition. Together, our findings reveal requirements for stable nucleosome formation and provide a foundation for further studies of the assembly and dynamics of native kinetochore complexes.

Single-nucleotide resolution detection of Topo IV cleavage activity in the Escherichia coli genome with Topo-Seq

Topoisomerase IV (Topo IV) is the main decatenation enzyme in Escherichia coli; it removes catenation links that are formed during DNA replication. Topo IV binding and cleavage sites were previously identified in the E. coli genome with ChIP-Seq and NorfIP. Here, we used a more sensitive, single-nucleotide resolution Topo-Seq procedure to identify Topo IV cleavage sites (TCSs) genome-wide. We detected thousands of TCSs scattered in the bacterial genome. The determined cleavage motif of Topo IV contained previously known cleavage determinants (−4G/+8C, −2A/+6 T, −1 T/+5A) and additional, not observed previously, positions −7C/+11G and −6C/+10G. TCSs were depleted in the Ter macrodomain except for two exceptionally strong non-canonical cleavage sites located in 33 and 38 bp from the XerC-box of the dif-site. Topo IV cleavage activity was increased in Left and Right macrodomains flanking the Ter macrodomain and was especially high in the 50–60 kb region containing the oriC origin of replication. Topo IV enrichment was also increased downstream of highly active transcription units, indicating that the enzyme is involved in relaxation of transcription-induced positive supercoiling.

Women’s contribution in understanding how topoisomerases, supercoiling, and transcription control genome organization

One of the biggest paradoxes in biology is that human genome is roughly 2 m long, while the nucleus containing it is almost one million times smaller. To fit into the nucleus, DNA twists, bends and folds into several hierarchical levels of compaction. Still, DNA has to maintain a high degree of accessibility to be readily replicated and transcribed by proteins. How compaction and accessibility co-exist functionally in human cells is still a matter of debate. Here, we discuss how the torsional stress of the DNA helix acts as a buffer, regulating both chromatin compaction and accessibility. We will focus on chromatin supercoiling and on the emerging role of topoisomerases as pivotal regulators of genome organization. We will mainly highlight the major breakthrough studies led by women, with the intention of celebrating the work of this group that remains a minority within the scientific community.

Variable DNA topology is an epigenetic generator of physiological heterogeneity in bacterial populations

Transcription is a noisy and stochastic process that produces sibling-to-sibling variations in physiology across a population of genetically identical cells. This pattern of diversity reflects, in part, the burst-like nature of transcription. Transcription bursting has many causes and a failure to remove the supercoils that accumulate in DNA during transcription elongation is an important contributor. Positive supercoiling of the DNA ahead of the transcription elongation complex can result in RNA polymerase stalling if this DNA topological roadblock is not removed. The relaxation of these positive supercoils is performed by the ATP-dependent type II topoisomerases DNA gyrase and topoisomerase IV. Interference with the action of these topoisomerases involving, inter alia, topoisomerase poisons, fluctuations in the [ATP]/[ADP] ratio, and/or the intervention of nucleoid-associated proteins with GapR-like or YejK-like activities, may have consequences for the smooth operation of the transcriptional machinery. Antibiotic-tolerant (but not resistant) persister cells are among the phenotypic outliers that may emerge. However, interference with type II topoisomerase activity can have much broader consequences, making it an important epigenetic driver of physiological diversity in the bacterial population.

Supercoiling-mediated feedback rapidly couples and tunes transcription

Transcription induces a wave of DNA supercoiling, altering the binding affinity of RNA polymerases and reshaping the biochemical landscape of gene regulation. As supercoiling rapidly diffuses, transcription dynamically reshapes the regulation of proximal genes, forming a complex feedback loop. However, a theoretical framework is needed to integrate biophysical regulation with biochemical transcriptional regulation. To investigate the role of supercoiling-mediated feedback within multi-gene systems, we model transcriptional regulation under the influence of supercoiling-mediated polymerase dynamics, allowing us to identify patterns of expression that result from physical inter-gene coupling. We find that gene syntax—the relative ordering and orientation of genes—defines the expression profiles, variance, burst dynamics, and inter-gene correlation of two-gene systems. Furthermore, supercoiling can enhance or weaken biochemical regulation. Our results suggest that supercoiling couples behavior between neighboring genes, providing a regulatory mechanism that tunes transcriptional variance in engineered gene networks and explains the behavior of co-localized native circuits.

Supercoiling-mediated feedback couples the transcription of proximal genes. Here, Johnstone and Galloway provide a framework for integrating biochemical gene regulation with the biophysical effects of DNA supercoiling. This unified model provides design principles for improving the performance of gene networks, developing novel regulatory functions, and accessing previously inaccessible regulatory dynamics.

TOP1 CAD-seq: A protocol to map catalytically engaged topoisomerase 1 in human cells

TOP1 CAD-seq enables mapping of TOP1 sites of covalent engagement with DNA. The procedure depends upon enrichment of DNA-covalent adducts using chaotropic salts and immunoprecipitation with an antibody specific for TOP1. Here, we describe a step-by-step protocol compatible with Illumina sequencing and bioinformatic pipeline for preliminary data analysis. Compared to other approaches for the genomic study of topoisomerases, TOP1 CAD-seq provides information about active TOP1 engaged on the DNA, taking advantage of low background due to absence of crosslinking.

For complete details on the use and execution of this protocol, please refer to Das et al. (2022).

•

TOP1 CAD-seq maps sites of topoisomerase 1 catalytic engagement with the DNA

•

The absence of a cross-linking step reduces the background signal

•

The method is high throughput and takes no longer than a regular ChIP-seq protocol

•

The method can be used as a proxy to assess the level of supercoiling in the genome

Publisher’s note: Undertaking any experimental protocol requires adherence to local institutional guidelines for laboratory safety and ethics.

Interaction between transcribing RNA polymerase and topoisomerase I prevents R-loop formation in E. coli

Bacterial topoisomerase I (TopoI) removes excessive negative supercoiling and is thought to relax DNA molecules during transcription, replication and other processes. Using ChIP-Seq, we show that TopoI of Escherichia coli (EcTopoI) is colocalized, genome-wide, with transcribing RNA polymerase (RNAP). Treatment with transcription elongation inhibitor rifampicin leads to EcTopoI relocation to promoter regions, where RNAP also accumulates. When a 14 kDa RNAP-binding EcTopoI C-terminal domain (CTD) is overexpressed, colocalization of EcTopoI and RNAP along the transcription units is reduced. Pull-down experiments directly show that the two enzymes interact in vivo. Using ChIP-Seq and Topo-Seq, we demonstrate that EcTopoI is enriched upstream (within up to 12-15 kb) of highly-active transcription units, indicating that EcTopoI relaxes negative supercoiling generated by transcription. Uncoupling of the RNAP:EcTopoI interaction by either overexpression of EcTopoI competitor (CTD or inactive EcTopoI Y319F mutant) or deletion of EcTopoI domains involved in the interaction is toxic for cells and leads to excessive negative plasmid supercoiling. Moreover, uncoupling of the RNAP:EcTopoI interaction leads to R-loops accumulation genome-wide, indicating that this interaction is required for prevention of R-loops formation.

In E. coli, disruption of TopoI and RNAP interaction decreases cells viability and leads to hypernegative DNA supercoiling and R loops accumulation. TopoI and DNA gyrase bind around transcription units and TopoI recognizes cleavage sites by a specific motif and negative supercoiling.

Condensin-driven loop extrusion on supercoiled DNA

Condensin, a structural maintenance of chromosomes (SMC) complex, has been shown to be a molecular motor protein that organizes chromosomes by extruding loops of DNA. In cells, such loop extrusion is challenged by many potential conflicts, for example, the torsional stresses that are generated by other DNA-processing enzymes. It has so far remained unclear how DNA supercoiling affects loop extrusion. Here, we use time-lapse single-molecule imaging to study condensin-driven DNA loop extrusion on supercoiled DNA. We find that condensin binding and DNA looping are stimulated by positively supercoiled DNA, and condensin preferentially binds near the tips of supercoiled plectonemes. Upon loop extrusion, condensin collects nearby plectonemes into a single supercoiled loop that is highly stable. Atomic force microscopy imaging shows that condensin generates supercoils in the presence of ATP. Our findings provide insight into the topology-regulated loading and formation of supercoiled loops by SMC complexes and clarify the interplay of loop extrusion and supercoiling.

Integrating cellular and molecular structures and dynamics into whole-cell models

A complete description of the state of the cell requires knowledge of its size, shape, components, intracellular reactions, and interactions with its environment-all of these as a function of time and cell growth. Adding to this list is the need for theoretical models and simulations that integrate and help to interpret this daunting amount of experimental data. It seems like an overwhelming list of requirements, but progress is being made on many fronts. In this review, we discuss the current challenges and problems in obtaining sufficient information about each aspect of a dynamical whole-cell model (DWCM) for simple and well-studied bacterial systems.

Psoralen mapping reveals a bacterial genome supercoiling landscape dominated by transcription

DNA supercoiling is a key regulator of all DNA metabolic processes including replication, transcription, and recombination, yet a reliable genomic assay for supercoiling is lacking. Here, we present a robust and flexible method (Psora-seq) to measure whole-genome supercoiling at high resolution. Using this tool in Escherichia coli, we observe a supercoiling landscape that is well correlated to transcription. Supercoiling twin-domains generated by RNA polymerase complexes span 25 kb in each direction - an order of magnitude farther than previous measurements in any organism. Thus, ribosomal and many other highly expressed genes strongly affect the topology of about 40 neighboring genes each, creating highly integrated gene circuits. Genomic patterns of supercoiling revealed by Psora-seq could be aptly predicted from modeling based on gene expression levels alone, indicating that transcription is the major determinant of chromosome supercoiling. Large-scale supercoiling patterns were highly symmetrical between left and right chromosome arms (replichores), indicating that DNA replication also strongly influences supercoiling. Skew in the axis of symmetry from the natural ori-ter axis supports previous indications that the rightward replication fork is delayed several minutes after initiation. Implications of supercoiling on DNA replication and chromosome domain structure are discussed.

DNA supercoiling-mediated collective behavior of co-transcribing RNA polymerases

Multiple RNA polymerases (RNAPs) transcribing a gene have been known to exhibit collective group behavior, causing the transcription elongation rate to increase with the rate of transcription initiation. Such behavior has long been believed to be driven by a physical interaction or ‘push’ between closely spaced RNAPs. However, recent studies have posited that RNAPs separated by longer distances may cooperate by modifying the DNA segment under transcription. Here, we present a theoretical model incorporating the mechanical coupling between RNAP translocation and the DNA torsional response. Using stochastic simulations, we demonstrate DNA supercoiling-mediated long-range cooperation between co-transcribing RNAPs. We find that inhibiting transcription initiation can slow down the already recruited RNAPs, in agreement with recent experimental observations, and predict that the average transcription elongation rate varies non-monotonically with the rate of transcription initiation. We further show that while RNAPs transcribing neighboring genes oriented in tandem can cooperate, those transcribing genes in divergent or convergent orientations can act antagonistically, and that such behavior holds over a large range of intergenic separations. Our model makes testable predictions, revealing how the mechanical interplay between RNAPs and the DNA they transcribe can govern transcriptional dynamics.

R-Loop Tracker: Web Access-Based Tool for R-Loop Detection and Analysis in Genomic DNA Sequences

R-loops are common non-B nucleic acid structures formed by a three-stranded nucleic acid composed of an RNA–DNA hybrid and a displaced single-stranded DNA (ssDNA) loop. Because the aberrant R-loop formation leads to increased mutagenesis, hyper-recombination, rearrangements, and transcription-replication collisions, it is regarded as important in human diseases. Therefore, its prevalence and distribution in genomes are studied intensively. However, in silico tools for R-loop prediction are limited, and therefore, we have developed the R-loop tracker tool, which was implemented as a part of the DNA Analyser web server. This new tool is focused upon (1) prediction of R-loops in genomic DNA without length and sequence limitations; (2) integration of R-loop tracker results with other tools for nucleic acids analyses, including Genome Browser; (3) internal cross-evaluation of in silico results with experimental data, where available; (4) easy export and correlation analyses with other genome features and markers; and (5) enhanced visualization outputs. Our new R-loop tracker tool is freely accessible on the web pages of DNA Analyser tools, and its implementation on the web-based server allows effective analyses not only for DNA segments but also for full chromosomes and genomes.

Bacterial transcription during growth arrest

Bacteria in most natural environments spend substantial periods of time limited for essential nutrients and not actively dividing. While transcriptional activity under these conditions is substantially reduced compared to that occurring during active growth, observations from diverse organisms and experimental approaches have shown that new transcription still occurs and is important for survival. Much of our understanding of transcription regulation has come from measuring transcripts in exponentially growing cells, or from in vitro experiments focused on transcription from highly active promoters by the housekeeping RNA polymerase holoenzyme. The fact that transcription during growth arrest occurs at low levels and is highly heterogeneous has posed challenges for its study. However, new methods of measuring low levels of gene expression activity, even in single cells, offer exciting opportunities for directly investigating transcriptional activity and its regulation during growth arrest. Furthermore, much of the rich structural and biochemical data from decades of work on the bacterial transcriptional machinery is also relevant to growth arrest. In this review, the physiological changes likely affecting transcription during growth arrest are first considered. Next, possible adaptations to help facilitate ongoing transcription during growth arrest are discussed. Finally, new insights from several recently published datasets investigating mRNA transcripts in single bacterial cells at various growth phases will be explored. Keywords: Growth arrest, stationary phase, RNA polymerase, nucleoid condensation, population heterogeneity.

High-resolution, genome-wide mapping of positive supercoiling in chromosomes

Supercoiling impacts DNA replication, transcription, protein binding to DNA, and the three-dimensional organization of chromosomes. However, there are currently no methods to directly interrogate or map positive supercoils, so their distribution in genomes remains unknown. Here, we describe a method, GapR-seq, based on the chromatin immunoprecipitation of GapR, a bacterial protein that preferentially recognizes overtwisted DNA, for generating high-resolution maps of positive supercoiling. Applying this method to Escherichia coli and Saccharomyces cerevisiae, we find that positive supercoiling is widespread, associated with transcription, and particularly enriched between convergently oriented genes, consistent with the 'twin-domain' model of supercoiling. In yeast, we also find positive supercoils associated with centromeres, cohesin-binding sites, autonomously replicating sites, and the borders of R-loops (DNA-RNA hybrids). Our results suggest that GapR-seq is a powerful approach, likely applicable in any organism, to investigate aspects of chromosome structure and organization not accessible by Hi-C or other existing methods.