Alternative DNA secondary structure formation affects RNA polymerase II promoter-proximal pausing in human

DNA fragility at topologically associated domain boundaries is promoted by alternative DNA secondary structure and topoisomerase II activity

CCCTC-binding factor (CTCF) binding sites are hotspots of genome instability. Although many factors have been associated with CTCF binding site fragility, no study has integrated all fragility-related factors to understand the mechanism(s) of how they work together. Using an unbiased, genome-wide approach, we found that DNA double-strand breaks (DSBs) are enriched at strong, but not weak, CTCF binding sites in five human cell types. Energetically favorable alternative DNA secondary structures underlie strong CTCF binding sites. These structures coincided with the location of topoisomerase II (TOP2) cleavage complex, suggesting that DNA secondary structure acts as a recognition sequence for TOP2 binding and cleavage at CTCF binding sites. Furthermore, CTCF knockdown significantly increased DSBs at strong CTCF binding sites and at CTCF sites that are located at topologically associated domain (TAD) boundaries. TAD boundary-associated CTCF sites that lost CTCF upon knockdown displayed increased DSBs when compared to the gained sites, and those lost sites are overrepresented with G-quadruplexes, suggesting that the structures act as boundary insulators in the absence of CTCF, and contribute to increased DSBs. These results model how alternative DNA secondary structures facilitate recruitment of TOP2 to CTCF binding sites, providing mechanistic insight into DNA fragility at CTCF binding sites.

Non-B-form DNA is associated with centromere stability in newly-formed polyploid wheat

Non-B-form DNA differs from the classic B-DNA double helix structure and plays a crucial regulatory role in replication and transcription. However, the role of non-B-form DNA in centromeres, especially in polyploid wheat, remains elusive. Here, we systematically analyzed seven non-B-form DNA motif profiles (A-phased DNA repeat, direct repeat, G-quadruplex, inverted repeat, mirror repeat, short tandem repeat, and Z-DNA) in hexaploid wheat. We found that three of these non-B-form DNA motifs were enriched at centromeric regions, especially at the CENH3-binding sites, suggesting that non-B-form DNA may create a favorable loading environment for the CENH3 nucleosome. To investigate the dynamics of centromeric non-B form DNA during the alloploidization process, we analyzed DNA secondary structure using CENH3 ChIP-seq data from newly formed allotetraploid wheat and its two diploid ancestors. We found that newly formed allotetraploid wheat formed more non-B-form DNA in centromeric regions compared with their parents, suggesting that non-B-form DNA is related to the localization of the centromeric regions in newly formed wheat. Furthermore, non-B-form DNA enriched in the centromeric regions was found to preferentially form on young LTR retrotransposons, explaining CENH3's tendency to bind to younger LTR. Collectively, our study describes the landscape of non-B-form DNA in the wheat genome, and sheds light on its potential role in the evolution of polyploid centromeres.

Bacillus subtilis stress-associated mutagenesis and developmental DNA repair

SUMMARYThe metabolic conditions that prevail during bacterial growth have evolved with the faithful operation of repair systems that recognize and eliminate DNA lesions caused by intracellular and exogenous agents. This idea is supported by the low rate of spontaneous mutations (10-9) that occur in replicating cells, maintaining genome integrity. In contrast, when growth and/or replication cease, bacteria frequently process DNA lesions in an error-prone manner. DNA repairs provide cells with the tools needed for maintaining homeostasis during stressful conditions and depend on the developmental context in which repair events occur. Thus, different physiological scenarios can be anticipated. In nutritionally stressed bacteria, different components of the base excision repair pathway may process damaged DNA in an error-prone approach, promoting genetic variability. Interestingly, suppressing the mismatch repair machinery and activating specific DNA glycosylases promote stationary-phase mutations. Current evidence also suggests that in resting cells, coupling repair processes to actively transcribed genes may promote multiple genetic transactions that are advantageous for stressed cells. DNA repair during sporulation is of interest as a model to understand how transcriptional processes influence the formation of mutations in conditions where replication is halted. Current reports indicate that transcriptional coupling repair-dependent and -independent processes operate in differentiating cells to process spontaneous and induced DNA damage and that error-prone synthesis of DNA is involved in these events. These and other noncanonical ways of DNA repair that contribute to mutagenesis, survival, and evolution are reviewed in this manuscript.

Tyrosine 1-phosphorylated RNA polymerase II transcribes PROMPTs to facilitate proximal promoter pausing and induce global transcriptional repression in response to DNA damage

DNA damage triggers a complex transcriptional response that involves both activation and repression of gene expression. In this study, we investigated global changes in transcription in response to ionizing irradiation (IR), which induces double-strand breaks in DNA. We used mNET-seq to profile nascent transcripts bound to different phosphorylated forms of the RNA polymerase II (RNA Pol II) C-terminal domain (CTD). We found that IR leads to global transcriptional repression of protein-coding genes, accompanied by an increase in antisense transcripts near promoters, called PROMPTs, transcribed by RNA Pol II phosphorylated on tyrosine 1 (Y1P) residue of the CTD. These Y1P-transcribed PROMPTs are enriched for PRC2 binding sites and associated with RNA Pol II proximal promoter pausing. We show the interaction between Y1P RNA Pol II and PRC2, as well as PRC2 binding to PROMPTs. Inhibition of PROMPTs or depletion of PRC2 leads to loss of transcriptional repression. Our results reveal a novel function of Y1P-dependent PROMPTs in mediating PRC2 recruitment to chromatin and RNA Pol II promoter pausing in response to DNA damage.

Genetic Diversity in Bronchial Asthma Susceptibility: Exploring the Role of Vitamin D Receptor Gene Polymorphisms in Varied Geographic Contexts

Bronchial asthma (BA) exhibits varying prevalence across global populations, prompting a comprehensive investigation into genetic and environmental determinants. Vitamin D is a potent immunomodulator capable of suppressing inflammatory signals in several cell types involved in the asthmatic response; it exerts effects on the immune system by binding to the nuclear vitamin D receptor (VDR). VDR gene genetic variations are affecting serum vitamin D levels with a possible role in the BA risk. The current study aimed to examine the complex interaction of various factors (genetic background, serum vitamin D levels, and geographic location) to identify differences in the influence of these factors on the susceptibility to asthma between populations at different latitudes. Focusing on Eastern European cohorts from Latvia and Lithuania and comparing them with published data on East Asian populations, we explore the impact of VDR gene polymorphisms on BA susceptibility. Genotyping four key VDR SNPs and assessing their association with 25-hydroxyvitamin D levels, our study unveils significant associations of the studied loci with the risk of asthma-both risk-reducing and increasing effects, differently distributed between Baltic and East Asian populations. The functional effects of in silico VDR gene genetic variations are also identified and discussed.

Combining a prioritization strategy and functional studies nominates 5'UTR variants underlying inherited retinal disease

BACKGROUND

5' untranslated regions (5'UTRs) are essential modulators of protein translation. Predicting the impact of 5'UTR variants is challenging and rarely performed in routine diagnostics. Here, we present a combined approach of a comprehensive prioritization strategy and functional assays to evaluate 5'UTR variation in two large cohorts of patients with inherited retinal diseases (IRDs).

METHODS

We performed an isoform-level re-analysis of retinal RNA-seq data to identify the protein-coding transcripts of 378 IRD genes with highest expression in retina. We evaluated the coverage of their 5'UTRs by different whole exome sequencing (WES) kits. The selected 5'UTRs were analyzed in whole genome sequencing (WGS) and WES data from IRD sub-cohorts from the 100,000 Genomes Project (n = 2397 WGS) and an in-house database (n = 1682 WES), respectively. Identified variants were annotated for 5'UTR-relevant features and classified into seven categories based on their predicted functional consequence. We developed a variant prioritization strategy by integrating population frequency, specific criteria for each category, and family and phenotypic data. A selection of candidate variants underwent functional validation using diverse approaches.

RESULTS

Isoform-level re-quantification of retinal gene expression revealed 76 IRD genes with a non-canonical retina-enriched isoform, of which 20 display a fully distinct 5'UTR compared to that of their canonical isoform. Depending on the probe design, 3-20% of IRD genes have 5'UTRs fully captured by WES. After analyzing these regions in both cohorts, we prioritized 11 (likely) pathogenic variants in 10 genes (ARL3, MERTK, NDP, NMNAT1, NPHP4, PAX6, PRPF31, PRPF4, RDH12, RD3), of which 7 were novel. Functional analyses further supported the pathogenicity of three variants. Mis-splicing was demonstrated for the PRPF31:c.-9+1G>T variant. The MERTK:c.-125G>A variant, overlapping a transcriptional start site, was shown to significantly reduce both luciferase mRNA levels and activity. The RDH12:c.-123C>T variant was found in cis with the hypomorphic RDH12:c.701G>A (p.Arg234His) variant in 11 patients. This 5'UTR variant, predicted to introduce an upstream open reading frame, was shown to result in reduced RDH12 protein but unaltered mRNA levels.

CONCLUSIONS

This study demonstrates the importance of 5'UTR variants implicated in IRDs and provides a systematic approach for 5'UTR annotation and validation that is applicable to other inherited diseases.

Assessing the Quality of Cotranscriptional Folding Simulations

Structural changes in RNAs are an important contributor to controlling gene expression not only at the posttranscriptional stage but also during transcription. A subclass of riboswitches and RNA thermometers located in the 5' region of the primary transcript regulates the downstream functional unit - usually an ORF - through premature termination of transcription. Not only such elements occur naturally, but they are also attractive devices in synthetic biology. The possibility to design such riboswitches or RNA thermometers is thus of considerable practical interest. Since these functional RNA elements act already during transcription, it is important to model and understand the dynamics of folding and, in particular, the formation of intermediate structures concurrently with transcription. Cotranscriptional folding simulations are therefore an important step to verify the functionality of design constructs before conducting expensive and labor-intensive wet lab experiments. For RNAs, full-fledged molecular dynamics simulations are far beyond practical reach because of both the size of the molecules and the timescales of interest. Even at the simplified level of secondary structures, further approximations are necessary. The BarMap approach is based on representing the secondary structure landscape for each individual transcription step by a coarse-grained representation that only retains a small set of low-energy local minima and the energy barriers between them. The folding dynamics between two transcriptional elongation steps is modeled as a Markov process on this representation. Maps between pairs of consecutive coarse-grained landscapes make it possible to follow the folding process as it changes in response to transcription elongation. In its original implementation, the BarMap software provides a general framework to investigate RNA folding dynamics on temporally changing landscapes. It is, however, difficult to use in particular for specific scenarios such as cotranscriptional folding. To overcome this limitation, we developed the user-friendly BarMap-QA pipeline described in detail in this contribution. It is illustrated here by an elaborate example that emphasizes the careful monitoring of several quality measures. Using an iterative workflow, a reliable and complete kinetics simulation of a synthetic, transcription-regulating riboswitch is obtained using minimal computational resources. All programs and scripts used in this contribution are free software and available for download as a source distribution for Linux® or as a platform-independent Docker® image including support for Apple macOS® and Microsoft Windows®.

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.

CHEX-seq detects single-cell genomic single-stranded DNA with catalytical potential

Genomic DNA (gDNA) undergoes structural interconversion between single- and double-stranded states during transcription, DNA repair and replication, which is critical for cellular homeostasis. We describe "CHEX-seq" which identifies the single-stranded DNA (ssDNA) in situ in individual cells. CHEX-seq uses 3'-terminal blocked, light-activatable probes to prime the copying of ssDNA into complementary DNA that is sequenced, thereby reporting the genome-wide single-stranded chromatin landscape. CHEX-seq is benchmarked in human K562 cells, and its utilities are demonstrated in cultures of mouse and human brain cells as well as immunostained spatially localized neurons in brain sections. The amount of ssDNA is dynamically regulated in response to perturbation. CHEX-seq also identifies single-stranded regions of mitochondrial DNA in single cells. Surprisingly, CHEX-seq identifies single-stranded loci in mouse and human gDNA that catalyze porphyrin metalation in vitro, suggesting a catalytic activity for genomic ssDNA. We posit that endogenous DNA enzymatic activity is a function of genomic ssDNA.

DNA language models are powerful predictors of genome-wide variant effects

The expanding catalog of genome-wide association studies (GWAS) provides biological insights across a variety of species, but identifying the causal variants behind these associations remains a significant challenge. Experimental validation is both labor-intensive and costly, highlighting the need for accurate, scalable computational methods to predict the effects of genetic variants across the entire genome. Inspired by recent progress in natural language processing, unsupervised pretraining on large protein sequence databases has proven successful in extracting complex information related to proteins. These models showcase their ability to learn variant effects in coding regions using an unsupervised approach. Expanding on this idea, we here introduce the Genomic Pre-trained Network (GPN), a model designed to learn genome-wide variant effects through unsupervised pretraining on genomic DNA sequences. Our model also successfully learns gene structure and DNA motifs without any supervision. To demonstrate its utility, we train GPN on unaligned reference genomes of Arabidopsis thaliana and seven related species within the Brassicales order and evaluate its ability to predict the functional impact of genetic variants in A. thaliana by utilizing allele frequencies from the 1001 Genomes Project and a comprehensive database of GWAS. Notably, GPN outperforms predictors based on popular conservation scores such as phyloP and phastCons. Our predictions for A. thaliana can be visualized as sequence logos in the UCSC Genome Browser (https://genome.ucsc.edu/s/gbenegas/gpn-arabidopsis). We provide code (https://github.com/songlab-cal/gpn) to train GPN for any given species using its DNA sequence alone, enabling unsupervised prediction of variant effects across the entire genome.

How does precursor RNA structure influence RNA processing and gene expression?

RNA is a fundamental biomolecule that has many purposes within cells. Due to its single-stranded and flexible nature, RNA naturally folds into complex and dynamic structures. Recent technological and computational advances have produced an explosion of RNA structural data. Many RNA structures have regulatory and functional properties. Studying the structure of nascent RNAs is particularly challenging due to their low abundance and long length, but their structures are important because they can influence RNA processing. Precursor RNA processing is a nexus of pathways that determines mature isoform composition and that controls gene expression. In this review, we examine what is known about human nascent RNA structure and the influence of RNA structure on processing of precursor RNAs. These known structures provide examples of how other nascent RNAs may be structured and show how novel RNA structures may influence RNA processing including splicing and polyadenylation. RNA structures can be targeted therapeutically to treat disease.

Roles of G4-DNA and G4-RNA in Class Switch Recombination and Additional Regulations in B-Lymphocytes

Mature B cells notably diversify immunoglobulin (Ig) production through class switch recombination (CSR), allowing the junction of distant "switch" (S) regions. CSR is initiated by activation-induced deaminase (AID), which targets cytosines adequately exposed within single-stranded DNA of transcribed targeted S regions, with a specific affinity for WRCY motifs. In mammals, G-rich sequences are additionally present in S regions, forming canonical G-quadruplexes (G4s) DNA structures, which favor CSR. Small molecules interacting with G4-DNA (G4 ligands), proved able to regulate CSR in B lymphocytes, either positively (such as for nucleoside diphosphate kinase isoforms) or negatively (such as for RHPS4). G4-DNA is also implicated in the control of transcription, and due to their impact on both CSR and transcriptional regulation, G4-rich sequences likely play a role in the natural history of B cell malignancies. Since G4-DNA stands at multiple locations in the genome, notably within oncogene promoters, it remains to be clarified how it can more specifically promote legitimate CSR in physiology, rather than pathogenic translocation. The specific regulatory role of G4 structures in transcribed DNA and/or in corresponding transcripts and recombination hereby appears as a major issue for understanding immune responses and lymphomagenesis.

Studying RNA Polymerase II Promoter-Proximal Pausing by In Vitro Immobilized Template and Transcription Assays

The immobilized template assay is a versatile biochemical method for studying protein-nucleic acid interactions. Using this method, immobilized nucleic acid-associated or specific proteins can be identified and quantified by techniques such as mass spectrometry and immunoblotting. Here, a modified immobilized template assay combined with in vitro transcription assay to study the function of transcription factors and transcriptional activities at the human heat shock protein 70 (HSP70) gene is described. Notably, this method can be applied to study other important genes and transcription factors in vitro.

DeePromClass: Delineator for Eukaryotic Core Promoters Employing Deep Neural Networks

Computational promoter identification in eukaryotes is a classical biological problem that should be refurbished with the availability of an avalanche of experimental data and emerging deep learning technologies. The current knowledge indicates that eukaryotic core promoters display multifarious signals such as TATA-Box, Inr element, TCT, and Pause-button, etc., and structural motifs such as G-quadruplexes. In the present study, we combined the power of deep learning with a plethora of promoter motifs to delineate promoter and non-promoters gleaned from the statistical properties of DNA sequence arrangement. To this end, we implemented convolutional neural network (CNN) and long short-term memory (LSTM) recurrent neural network architecture for five model systems with [-100 to +50] segments relative to the transcription start site being the core promoter. Unlike previous state-of-the-art tools, which furnish a binary decision of promoter or non-promoter, we classify a chunk of 151mer sequence into a promoter along with the consensus signal type or a non-promoter. The combined CNN-LSTM model; we call "DeePromClass", achieved testing accuracy of 90.6%, 93.6%, 91.8%, 86.5%, and 84.0% for S. cerevisiae, C. elegans, D. melanogaster, Mus musculus, and Homo sapiens respectively. In total, our tool provides an insightful update on next-generation promoter prediction tools for promoter biologists.

Dysregulation of Plasma miR-146a and miR-155 Expression Profile in Mycosis Fungoides Is Associated with rs2910164 and rs767649 Polymorphisms

Diagnosis of Mycosis Fungoides (MF) may be challenging, due to its polymorphic nature. The use of miRNAs as biomarkers to assist in diagnosis has been investigated, mainly in skin lesion biopsies. The purpose of this study is to evaluate the plasma levels of miR-146a and miR-155 in MF patients and to investigate their association with SNPs of their genes. Plasma miRNAs were quantified by RT-qPCR. Genomic DNA was used for SNPs’ genotyping by Sanger sequencing. Plasma levels of miR-146a and miR-155 were significantly higher in patients vs. controls, in early MF patients vs. controls, and in advanced vs. early MF patients. Both miRNAs’ levels were significantly higher in stage IIB vs. early-stage patients. miR-155 plasma levels were significantly higher in patients with skin tumors or erythroderma. CC genotype (rs2910164 C>G) was significantly more frequent in healthy controls and associated with lower MF risk and lower miR-146a levels. The AA genotype (rs767649 T>A) was significantly more frequent in patients and correlated with increased MF risk and increased miR-155 levels. The combination of GG+AA was only detected in patients and was correlated with higher MF susceptibility. Increased mir-146a and mir-155 plasma levels in MF is an important finding to establish putative noninvasive biomarkers. The presence of SNPs is closely associated with miRs’ expression, and possibly with disease susceptibility.

High-throughput techniques enable advances in the roles of DNA and RNA secondary structures in transcriptional and post-transcriptional gene regulation

The most stable structure of DNA is the canonical right-handed double helix termed B DNA. However, certain environments and sequence motifs favor alternative conformations, termed non-canonical secondary structures. The roles of DNA and RNA secondary structures in transcriptional regulation remain incompletely understood. However, advances in high-throughput assays have enabled genome wide characterization of some secondary structures. Here, we describe their regulatory functions in promoters and 3’UTRs, providing insights into key mechanisms through which they regulate gene expression. We discuss their implication in human disease, and how advances in molecular technologies and emerging high-throughput experimental methods could provide additional insights.

DSIF modulates RNA polymerase II occupancy according to template G + C content

The DSIF complex comprising the Supt4h and Supt5h transcription elongation proteins clamps RNA polymerase II (RNAPII) onto DNA templates, facilitating polymerase processivity. Lowering DSIF components can differentially decrease expression of alleles containing nucleotide repeat expansions, suggesting that RNAPII transit through repeat expansions is dependent on DSIF functions. To globally identify sequence features that affect dependence of the polymerase on DSIF in human cells, we used ultra-deep ChIP-seq analysis and RNA-seq to investigate and quantify the genome-wide effects of Supt4h loss on template occupancy and transcript production. Our results indicate that RNAPII dependence on Supt4h varies according to G + C content. Effects of DSIF knockdown were prominent during transcription of sequences high in G + C but minimal for sequences low in G + C and were particularly evident for G + C-rich segments of long genes. Reanalysis of previously published ChIP-seq data obtained from mouse cells showed similar effects of template G + C composition on Supt5h actions. Our evidence that DSIF dependency varies globally in different template regions according to template sequence composition suggests that G + C content may have a role in the selectivity of Supt4h knockdown and Supt5h knockdown during transcription of gene alleles containing expansions of G + C-rich repeats.

Delineation of the DNA Structural Features of Eukaryotic Core Promoter Classes

The eukaryotic transcription is orchestrated from a chunk of the DNA region stated as the core promoter. Multifarious and punctilious core promoter signals, viz., TATA-box, Inr, BREs, and Pause Button, are associated with a subset of genes and regulate their spatiotemporal expression. However, the core promoter architecture linked with these signals has not been investigated exhaustively for several species. In this study, we attempted to envisage the adaptive binding landscape of the transcription initiation machinery as a function of DNA structure. To this end, we deployed a set of k-mer based DNA structural estimates and regular expression models derived from experiments, molecular dynamic simulations, and theoretical frameworks, and high-throughout promoter data sets retrieved from the eukaryotic promoter database. We categorized protein-coding gene core promoters based on characteristic motifs at precise locations and analyzed the B-DNA structural properties and non-B-DNA structural motifs for 15 different eukaryotic genomes. We observed that Inr, BREd, and no-motif classes display common patterns of DNA sequence and structural environment. TATA-containing, BREu, and Pause Button classes show a deviant behavior with the TATA class displaying varied axial and twisting flexibility while BREu and Pause Button leaned toward G-quadruplex motif enrichment. Intriguingly, DNA meltability and shape signals are conserved irrespective of the presence or absence of distinct core promoter motifs in the majority of species. Altogether, here we delineated the conserved DNA structural signals associated with several promoter classes that may contribute to the chromatin configuration, orchestration of transcription machinery, and DNA duplex melting during the transcription process.

Unique properties of spacer acquisition by the type III-A CRISPR-Cas system

Type III CRISPR-Cas systems have a unique mode of interference, involving crRNA-guided recognition of nascent RNA and leading to DNA and RNA degradation. How type III systems acquire new CRISPR spacers is currently not well understood. Here, we characterize CRISPR spacer uptake by a type III-A system within its native host, Streptococcus thermophilus. Adaptation by the type II-A system in the same host provided a basis for comparison. Cas1 and Cas2 proteins were critical for type III adaptation but deletion of genes responsible for crRNA biogenesis or interference did not detectably change spacer uptake patterns, except those related to host counter-selection. Unlike the type II-A system, type III spacers are acquired in a PAM- and orientation-independent manner. Interestingly, certain regions of plasmids and the host genome were particularly well-sampled during type III-A, but not type II-A, spacer uptake. These regions included the single-stranded origins of rolling-circle replicating plasmids, rRNA and tRNA encoding gene clusters, promoter regions of expressed genes and 5′ UTR regions involved in transcription attenuation. These features share the potential to form DNA secondary structures, suggesting a preferred substrate for type III adaptation. Lastly, the type III-A system adapted to and protected host cells from lytic phage infection.

Timing RNA polymerase pausing with TV-PRO-seq

Transcription of many genes in metazoans is subject to polymerase pausing, which is the transient stop of transcriptionally engaged polymerases. This is known to mainly occur in promoter-proximal regions but it is not well understood. In particular, a genome-wide measurement of pausing times at high resolution has been lacking. We present here the time-variant precision nuclear run-on and sequencing (TV-PRO-seq) assay, an extension of the standard PRO-seq that allows us to estimate genome-wide pausing times at single-base resolution. Its application to human cells demonstrates that, proximal to promoters, polymerases pause more frequently but for shorter times than in other genomic regions. Comparison with single-cell gene expression data reveals that the polymerase pausing times are longer in highly expressed genes, while transcriptionally noisier genes have higher pausing frequencies and slightly longer pausing times. Analyses of histone modifications suggest that the marker H3K36me3 is related to the polymerase pausing.

A Bayesian model based computational analysis of the relationship between bisulfite accessible single-stranded DNA in chromatin and somatic hypermutation of immunoglobulin genes

The B cells in our body generate protective antibodies by introducing somatic hypermutations (SHM) into the variable region of immunoglobulin genes (IgVs). The mutations are generated by activation induced deaminase (AID) that converts cytosine to uracil in single stranded DNA (ssDNA) generated during transcription. Attempts have been made to correlate SHM with ssDNA using bisulfite to chemically convert cytosines that are accessible in the intact chromatin of mutating B cells. These studies have been complicated by using different definitions of "bisulfite accessible regions" (BARs). Recently, deep-sequencing has provided much larger datasets of such regions but computational methods are needed to enable this analysis. Here we leveraged the deep-sequencing approach with unique molecular identifiers and developed a novel Hidden Markov Model based Bayesian Segmentation algorithm to characterize the ssDNA regions in the IGHV4-34 gene of the human Ramos B cell line. Combining hierarchical clustering and our new Bayesian model, we identified recurrent BARs in certain subregions of both top and bottom strands of this gene. Using this new system, the average size of BARs is about 15 bp. We also identified potential G-quadruplex DNA structures in this gene and found that the BARs co-locate with G-quadruplex structures in the opposite strand. Using various correlation analyses, there is not a direct site-to-site relationship between the bisulfite accessible ssDNA and all sites of SHM but most of the highly AID mutated sites are within 15 bp of a BAR. In summary, we developed a novel platform to study single stranded DNA in chromatin at a base pair resolution that reveals potential relationships among BARs, SHM and G-quadruplexes. This platform could be applied to genome wide studies in the future.

G-Quadruplex in Gene Encoding Large Subunit of Plant RNA Polymerase II: A Billion-Year-Old Story

G-quadruplexes have long been perceived as rare and physiologically unimportant nucleic acid structures. However, several studies have revealed their importance in molecular processes, suggesting their possible role in replication and gene expression regulation. Pathways involving G-quadruplexes are intensively studied, especially in the context of human diseases, while their involvement in gene expression regulation in plants remains largely unexplored. Here, we conducted a bioinformatic study and performed a complex circular dichroism measurement to identify a stable G-quadruplex in the gene RPB1, coding for the RNA polymerase II large subunit. We found that this G-quadruplex-forming locus is highly evolutionarily conserved amongst plants sensu lato (Archaeplastida) that share a common ancestor more than one billion years old. Finally, we discussed a new hypothesis regarding G-quadruplexes interacting with UV light in plants to potentially form an additional layer of the regulatory network.

Genome-wide profiling reveals functional interplay of DNA sequence composition, transcriptional activity, and nucleosome positioning in driving DNA supercoiling and helix destabilization in C. elegans

DNA topology and alternative DNA structures are implicated in regulating diverse biological processes. Although biomechanical properties of these structures have been studied extensively in vitro, characterization in vivo, particularly in multicellular organisms, is limited. We devised new methods to map DNA supercoiling and single-stranded DNA in Caenorhabditis elegans embryos and diapause larvae. To map supercoiling, we quantified the incorporation of biotinylated psoralen into DNA using high-throughput sequencing. To map single-stranded DNA, we combined permanganate treatment with genome-wide sequencing of induced double-stranded breaks. We found high levels of negative supercoiling at transcription start sites (TSSs) in embryos. GC-rich regions flanked by a sharp GC-to-AT transition delineate boundaries of supercoil propagation. In contrast to TSSs in embryos, TSSs in diapause larvae showed dramatic reductions in negative supercoiling without concomitant attenuation of transcription, suggesting developmental-stage-specific regulation. To assess whether alternative DNA structures control chromosome architecture and gene expression, we examined DNA supercoiling in the context of X-Chromosome dosage compensation. We showed that the condensin dosage compensation complex creates negative supercoils locally at its highest-occupancy binding sites but found no evidence for large-scale supercoiling domains along X Chromosomes. In contrast to transcription-coupled negative supercoiling, single-strandedness, which is most pronounced at transcript end sites, is dependent on high AT content and symmetrically positioned nucleosomes. We propose that sharp transitions in sequence composition at functional genomic elements constitute a common regulatory code and that DNA structure and propagation of torsional stress at regulatory elements are critical parameters in shaping important developmental events.

Constrained G4 structures unveil topology specificity of known and new G4 binding proteins

G-quadruplexes (G4) are non-canonical secondary structures consisting in stacked tetrads of hydrogen-bonded guanines bases. An essential feature of G4 is their intrinsic polymorphic nature, which is characterized by the equilibrium between several conformations (also called topologies) and the presence of different types of loops with variable lengths. In cells, G4 functions rely on protein or enzymatic factors that recognize and promote or resolve these structures. In order to characterize new G4-dependent mechanisms, extensive researches aimed at identifying new G4 binding proteins. Using G-rich single-stranded oligonucleotides that adopt non-controlled G4 conformations, a large number of G4-binding proteins have been identified in vitro, but their specificity towards G4 topology remained unknown. Constrained G4 structures are biomolecular objects based on the use of a rigid cyclic peptide scaffold as a template for directing the intramolecular assembly of the anchored oligonucleotides into a single and stabilized G4 topology. Here, using various constrained RNA or DNA G4 as baits in human cell extracts, we establish the topology preference of several well-known G4-interacting factors. Moreover, we identify new G4-interacting proteins such as the NELF complex involved in the RNA-Pol II pausing mechanism, and we show that it impacts the clastogenic effect of the G4-ligand pyridostatin.

Non-B DNA-Forming Motifs Promote Mfd-Dependent Stationary-Phase Mutagenesis in Bacillus subtilis

Transcription-induced mutagenic mechanisms limit genetic changes to times when expression happens and to coding DNA. It has been hypothesized that intrinsic sequences that have the potential to form alternate DNA structures, such as non-B DNA structures, influence these mechanisms. Non-B DNA structures are promoted by transcription and induce genome instability in eukaryotic cells, but their impact in bacterial genomes is less known. Here, we investigated if G4 DNA- and hairpin-forming motifs influence stationary-phase mutagenesis in Bacillus subtilis. We developed a system to measure the influence of non-B DNA on B. subtilis stationary-phase mutagenesis by deleting the wild-type argF at its chromosomal position and introducing IPTG-inducible argF alleles differing in their ability to form hairpin and G4 DNA structures into an ectopic locus. Using this system, we found that sequences predicted to form non-B DNA structures promoted mutagenesis in B. subtilis stationary-phase cells; such a response did not occur in growing conditions. We also found that the transcription-coupled repair factor Mfd promoted mutagenesis at these predicted structures. In summary, we showed that non-B DNA-forming motifs promote genetic instability, particularly in coding regions in stressed cells; therefore, non-B DNA structures may have a spatial and temporal mutagenic effect in bacteria. This study provides insights into mechanisms that prevent or promote mutagenesis and advances our understanding of processes underlying bacterial evolution.

Conserved DNA sequence features underlie pervasive RNA polymerase pausing

Pausing of transcribing RNA polymerase is regulated and creates opportunities to control gene expression. Research in metazoans has so far mainly focused on RNA polymerase II (Pol II) promoter-proximal pausing leaving the pervasive nature of pausing and its regulatory potential in mammalian cells unclear. Here, we developed a pause detecting algorithm (PDA) for nucleotide-resolution occupancy data and a new native elongating transcript sequencing approach, termed nested NET-seq, that strongly reduces artifactual peaks commonly misinterpreted as pausing sites. Leveraging PDA and nested NET-seq reveal widespread genome-wide Pol II pausing at single-nucleotide resolution in human cells. Notably, the majority of Pol II pauses occur outside of promoter-proximal gene regions primarily along the gene-body of transcribed genes. Sequence analysis combined with machine learning modeling reveals DNA sequence properties underlying widespread transcriptional pausing including a new pause motif. Interestingly, key sequence determinants of RNA polymerase pausing are conserved between human cells and bacteria. These studies indicate pervasive sequence-induced transcriptional pausing in human cells and the knowledge of exact pause locations implies potential functional roles in gene expression.

A dual role for H2A.Z.1 in modulating the dynamics of RNA polymerase II initiation and elongation

RNA polymerase II (RNAPII) pausing immediately downstream of the transcription start site is a critical rate-limiting step for the expression of most metazoan genes. During pause release, RNAPII encounters a highly conserved +1 H2A.Z nucleosome, yet how this histone variant contributes to transcription is poorly understood. Here, using an inducible protein degron system combined with genomic approaches and live cell super-resolution microscopy, we show that H2A.Z.1 modulates RNAPII dynamics across most genes in murine embryonic stem cells. Our quantitative analysis shows that H2A.Z.1 slows the rate of RNAPII pause release and consequently impacts negative elongation factor dynamics as well as nascent transcription. Consequently, H2A.Z.1 also impacts re-loading of the pre-initiation complex components TFIIB and TBP. Altogether, this work provides a critical mechanistic link between H2A.Z.1 and the proper induction of mammalian gene expression programs through the regulation of RNAPII dynamics and pause release.

Evolution of Diverse Strategies for Promoter Regulation

DNA is fundamentally important for all cellular organisms due to its role as a store of hereditary genetic information. The precise and accurate regulation of gene transcription depends primarily on promoters, which vary significantly within and between genomes. Some promoters are rich in specific types of bases, while others have more varied, complex sequence characteristics. However, it is not only base sequence but also epigenetic modifications and altered DNA structure that regulate promoter activity. Significantly, many promoters across all organisms contain sequences that can form intrastrand hairpins (cruciforms) or four-stranded structures (G-quadruplex or i-motif). In this review we integrate recent studies on promoter regulation that highlight the importance of DNA structure in the evolutionary adaptation of promoter sequences.

Short-range template switching in great ape genomes explored using pair hidden Markov models

Many complex genomic rearrangements arise through template switch errors, which occur in DNA replication when there is a transient polymerase switch to an alternate template nearby in three-dimensional space. While typically investigated at kilobase-to-megabase scales, the genomic and evolutionary consequences of this mutational process are not well characterised at smaller scales, where they are often interpreted as clusters of independent substitutions, insertions and deletions. Here we present an improved statistical approach using pair hidden Markov models, and use it to detect and describe short-range template switches underlying clusters of mutations in the multi-way alignment of hominid genomes. Using robust statistics derived from evolutionary genomic simulations, we show that template switch events have been widespread in the evolution of the great apes’ genomes and provide a parsimonious explanation for the presence of many complex mutation clusters in their phylogenetic context. Larger-scale mechanisms of genome rearrangement are typically associated with structural features around breakpoints, and accordingly we show that atypical patterns of secondary structure formation and DNA bending are present at the initial template switch loci. Our methods improve on previous non-probabilistic approaches for computational detection of template switch mutations, allowing the statistical significance of events to be assessed. By specifying realistic evolutionary parameters based on the genomes and taxa involved, our methods can be readily adapted to other intra- or inter-species comparisons.

DNA replication is an imperfect process which causes the mutations that give rise to genetic diversity during the evolution of genomes. While many mutations are independent, single-nucleotide substitutions or small insertions and deletions, some mutations arise as nonindependent clusters of substitutions and larger scale chromosomal rearrangements. Large-scale rearrangements (also called structural variants) in particular can have a profound impact on genome evolution and contribute to both germline and somatic disease in humans. The replication-based mechanisms underlying structural variation typically involve a polymerase switch event in which a large segment of DNA is copied using a template from an alternate location in the genome. Methods for identifying these template switch mutations lack the power to detect smaller scale rearrangements which can arise through the same replication-based pathways. Here we outline a model which can detect and assess the statistical significance of such small-scale template switches within their evolutionary context. We show that these events are widespread in the evolution of great apes and that the genomic features associated with these small-scale rearrangements are similar to those of large-scale structural variants.

Distinct properties and functions of CTCF revealed by a rapidly inducible degron system

CCCTC-binding factor (CTCF) is a conserved zinc finger transcription factor implicated in a wide range of functions, including genome organization, transcription activation, and elongation. To explore the basis for CTCF functional diversity, we coupled an auxin-induced degron system with precision nuclear run-on. Unexpectedly, oriented CTCF motifs in gene bodies are associated with transcriptional stalling in a manner independent of bound CTCF. Moreover, CTCF at different binding sites (CBSs) displays highly variable resistance to degradation. Motif sequence does not significantly predict degradation behavior, but location at chromatin boundaries and chromatin loop anchors, as well as co-occupancy with cohesin, are associated with delayed degradation. Single-molecule tracking experiments link chromatin residence time to CTCF degradation kinetics, which has ramifications regarding architectural CTCF functions. Our study highlights the heterogeneity of CBSs, uncovers properties specific to architecturally important CBSs, and provides insights into the basic processes of genome organization and transcription regulation.

Causes and consequences of RNA polymerase II stalling during transcript elongation

The journey of RNA polymerase II (Pol II) as it transcribes a gene is anything but a smooth ride. Transcript elongation is discontinuous and can be perturbed by intrinsic regulatory barriers, such as promoter-proximal pausing, nucleosomes, RNA secondary structures and the underlying DNA sequence. More substantial blocking of Pol II translocation can be caused by other physiological circumstances and extrinsic obstacles, including other transcribing polymerases, the replication machinery and several types of DNA damage, such as bulky lesions and DNA double-strand breaks. Although numerous different obstacles cause Pol II stalling or arrest, the cell somehow distinguishes between them and invokes different mechanisms to resolve each roadblock. Resolution of Pol II blocking can be as straightforward as temporary backtracking and transcription elongation factor S-II (TFIIS)-dependent RNA cleavage, or as drastic as premature transcription termination or degradation of polyubiquitylated Pol II and its associated nascent RNA. In this Review, we discuss the current knowledge of how these different Pol II stalling contexts are distinguished by the cell, how they overlap with each other, how they are resolved and how, when unresolved, they can cause genome instability.

Topoisomerase II contributes to DNA secondary structure-mediated double-stranded breaks

DNA double-stranded breaks (DSBs) trigger human genome instability, therefore identifying what factors contribute to DSB induction is critical for our understanding of human disease etiology. Using an unbiased, genome-wide approach, we found that genomic regions with the ability to form highly stable DNA secondary structures are enriched for endogenous DSBs in human cells. Human genomic regions predicted to form non-B-form DNA induced gross chromosomal rearrangements in yeast and displayed high indel frequency in human genomes. The extent of instability in both analyses is in concordance with the structure forming ability of these regions. We also observed an enrichment of DNA secondary structure-prone sites overlapping transcription start sites (TSSs) and CCCTC-binding factor (CTCF) binding sites, and uncovered an increase in DSBs at highly stable DNA secondary structure regions, in response to etoposide, an inhibitor of topoisomerase II (TOP2) re-ligation activity. Importantly, we found that TOP2 deficiency in both yeast and human leads to a significant reduction in DSBs at structure-prone loci, and that sites of TOP2 cleavage have a greater ability to form highly stable DNA secondary structures. This study reveals a direct role for TOP2 in generating secondary structure-mediated DNA fragility, advancing our understanding of mechanisms underlying human genome instability.

Oncogenic seRNA functional activation: a novel mechanism of tumorigenesis

seRNA is a noncoding RNA (ncRNA) transcribed from active super-enhancer (SE), through which SE exerts biological functions and participates in various physiological and pathological processes. seRNA recruits cofactor, RNA polymerase II and mediator to constitute and stabilize chromatin loop SE and promoter region, which regulates target genes transcription. In tumorigenesis, DNA insertion, deletion, translocation, focal amplification and carcinogen factor mediate oncogenic SE generation, meanwhile, oncogenic SE transcribes into tumor-related seRNA, termed as oncogenic seRNA. Oncogenic seRNA participates in tumorigenesis through activating various signal-pathways. The recent reports showed that oncogenic seRNA implicates in a widespread range of cytopathological processes in cancer progression including cell proliferation, apoptosis, autophagy, epithelial-mesenchymal transition, extracellular matrix stiffness and angiogenesis. In this article, we comprehensively summarized seRNA’s characteristics and functions, and emphatically introduced inducible formation of oncogenic seRNA and its functional mechanisms. Lastly, some research strategies on oncogenic seRNA were introduced, and the perspectives on cancer therapy that targets oncogenic seRNA were also discussed.

Pausing sites of RNA polymerase II on actively transcribed genes are enriched in DNA double-stranded breaks

DNA double-stranded breaks (DSBs) are strongly associated with active transcription, and promoter-proximal pausing of RNA polymerase II (Pol II) is a critical step in transcriptional regulation. Mapping the distribution of DSBs along actively expressed genes and identifying the location of DSBs relative to pausing sites can provide mechanistic insights into transcriptional regulation. Using genome-wide DNA break mapping/sequencing techniques at single-nucleotide resolution in human cells, we found that DSBs are preferentially located around transcription start sites of highly transcribed and paused genes and that Pol II promoter-proximal pausing sites are enriched in DSBs. We observed that DSB frequency at pausing sites increases as the strength of pausing increases, regardless of whether the pausing sites are near or far from annotated transcription start sites. Inhibition of topoisomerase I and II by camptothecin and etoposide treatment, respectively, increased DSBs at the pausing sites as the concentrations of drugs increased, demonstrating the involvement of topoisomerases in DSB generation at the pausing sites. DNA breaks generated by topoisomerases are short-lived because of the religation activity of these enzymes, which these drugs inhibit; therefore, the observation of increased DSBs with increasing drug doses at pausing sites indicated active recruitment of topoisomerases to these sites. Furthermore, the enrichment and locations of DSBs at pausing sites were shared among different cell types, suggesting that Pol II promoter-proximal pausing is a common regulatory mechanism. Our findings support a model in which topoisomerases participate in Pol II promoter-proximal pausing and indicated that DSBs at pausing sites contribute to transcriptional activation.

Broken by the Cut: A Journey into the Role of Topoisomerase II in DNA Fragility

DNA topoisomerase II (TOP2) plays a critical role in many processes such as replication and transcription, where it resolves DNA structures and relieves torsional stress. Recent evidence demonstrated the association of TOP2 with topologically associated domains (TAD) boundaries and CCCTC-binding factor (CTCF) binding sites. At these sites, TOP2 promotes interactions between enhancers and gene promoters, and relieves torsional stress that accumulates at these physical barriers. Interestingly, in executing its enzymatic function, TOP2 contributes to DNA fragility through re-ligation failure, which results in persistent DNA breaks when unrepaired or illegitimately repaired. Here, we discuss the biological processes for which TOP2 is required and the steps at which it can introduce DNA breaks. We describe the repair processes that follow removal of TOP2 adducts and the resultant broken DNA ends, and present how these processes can contribute to disease-associated mutations. Furthermore, we examine the involvement of TOP2-induced breaks in the formation of oncogenic translocations of leukemia and papillary thyroid cancer, as well as the role of TOP2 and proteins which repair TOP2 adducts in other diseases. The participation of TOP2 in generating persistent DNA breaks and leading to diseases such as cancer, could have an impact on disease treatment and prevention.

Promoter-proximal pausing of RNA polymerase II: a nexus of gene regulation

In this review, Core et al. discuss the recent advances in our understanding of the early steps in Pol II transcription, highlighting the events and factors involved in the establishment and release of paused Pol II. They also discuss a number of unanswered questions about the regulation and function of Pol II pausing.

Precise spatio–temporal control of gene activity is essential for organismal development, growth, and survival in a changing environment. Decisive steps in gene regulation involve the pausing of RNA polymerase II (Pol II) in early elongation, and the controlled release of paused polymerase into productive RNA synthesis. Here we describe the factors that enable pausing and the events that trigger Pol II release into the gene. We also discuss open questions in the field concerning the stability of paused Pol II, nucleosomes as obstacles to elongation, and potential roles of pausing in defining the precision and dynamics of gene expression.

The Presence and Localization of G-Quadruplex Forming Sequences in the Domain of Bacteria

The role of local DNA structures in the regulation of basic cellular processes is an emerging field of research. Amongst local non-B DNA structures, the significance of G-quadruplexes was demonstrated in the last decade, and their presence and functional relevance has been demonstrated in many genomes, including humans. In this study, we analyzed the presence and locations of G-quadruplex-forming sequences by G4Hunter in all complete bacterial genomes available in the NCBI database. G-quadruplex-forming sequences were identified in all species, however the frequency differed significantly across evolutionary groups. The highest frequency of G-quadruplex forming sequences was detected in the subgroup Deinococcus-Thermus, and the lowest frequency in Thermotogae. G-quadruplex forming sequences are non-randomly distributed and are favored in various evolutionary groups. G-quadruplex-forming sequences are enriched in ncRNA segments followed by mRNAs. Analyses of surrounding sequences showed G-quadruplex-forming sequences around tRNA and regulatory sequences. These data point to the unique and non-random localization of G-quadruplex-forming sequences in bacterial genomes.

The role of fork stalling and DNA structures in causing chromosome fragility

Alternative non-B form DNA structures, also called secondary structures, can form in certain DNA sequences under conditions that produce single-stranded DNA, such as during replication, transcription, and repair. Direct links between secondary structure formation, replication fork stalling, and genomic instability have been found for many repeated DNA sequences that cause disease when they expand. Common fragile sites (CFSs) are known to be AT-rich and break under replication stress, yet the molecular basis for their fragility is still being investigated. Over the past several years, new evidence has linked both the formation of secondary structures and transcription to fork stalling and fragility of CFSs. How these two events may synergize to cause fragility and the role of nuclease cleavage at secondary structures in rare and CFSs are discussed here. We also highlight evidence for a new hypothesis that secondary structures at CFSs not only initiate fragility but also inhibit healing, resulting in their characteristic appearance.

Developing Novel G-Quadruplex Ligands: from Interaction with Nucleic Acids to Interfering with Nucleic Acid⁻Protein Interaction

G-quadruplex is a special secondary structure of nucleic acids in guanine-rich sequences of genome. G-quadruplexes have been proved to be involved in the regulation of replication, DNA damage repair, and transcription and translation of oncogenes or other cancer-related genes. Therefore, targeting G-quadruplexes has become a novel promising anti-tumor strategy. Different kinds of small molecules targeting the G-quadruplexes have been designed, synthesized, and identified as potential anti-tumor agents, including molecules directly bind to the G-quadruplex and molecules interfering with the binding between the G-quadruplex structures and related binding proteins. This review will explore the feasibility of G-quadruplex ligands acting as anti-tumor drugs, from basis to application. Meanwhile, since helicase is the most well-defined G-quadruplex-related protein, the most extensive research on the relationship between helicase and G-quadruplexes, and its meaning in drug design, is emphasized.