Dynamic alternative DNA structures in biology and disease

kmerDB: A database encompassing the set of genomic and proteomic sequence information for each species

The decrease in sequencing expenses has facilitated the creation of reference genomes and proteomes for an expanding array of organisms. Nevertheless, no established repository that details organism-specific genomic and proteomic sequences of specific lengths, referred to as kmers, exists to our knowledge. In this article, we present kmerDB, a database accessible through an interactive web interface that provides kmer-based information from genomic and proteomic sequences in a systematic way. kmerDB currently contains 202,340,859,107 base pairs and 19,304,903,356 amino acids, spanning 54,039 and 21,865 reference genomes and proteomes, respectively, as well as 6,905,362 and 149,305,183 genomic and proteomic species-specific sequences, termed quasi-primes. Additionally, we provide access to 5,186,757 nucleic and 214,904,089 peptide sequences absent from every genome and proteome, termed primes. kmerDB features a user-friendly interface offering various search options and filters for easy parsing and searching. The service is available at: www.kmerdb.com.

Methods and applications of genome-wide profiling of DNA damage and rare mutations

DNA damage is a threat to genome integrity and can be a cause of many human diseases, owing to either changes in the chemical structure of DNA or conversion of the damage into a mutation, that is, a permanent change in DNA sequence. Determining the exact positions of DNA damage and ensuing mutations in the genome are important for identifying mechanisms of disease aetiology when characteristic mutations are prevalent and probably causative in a particular disease. However, this approach is challenging particularly when levels of DNA damage are low, for example, as a result of chronic exposure to environmental agents or certain endogenous processes, such as the generation of reactive oxygen species. Over the past few years, a comprehensive toolbox of genome-wide methods has been developed for the detection of DNA damage and rare mutations at single-nucleotide resolution in mammalian cells. Here, we review and compare these methods, describe their current applications and discuss future research questions that can now be addressed.

Controlling genome topology with sequences that trigger post-replication gap formation during replisome passage: the E. coli RRS elements

We report that the Escherichia coli chromosome includes novel GC-rich genomic structural elements that trigger formation of post-replication gaps upon replisome passage. The two nearly perfect 222 bp repeats, designated Replication Risk Sequences or RRS, are each 650 kb from the terminus sequence dif and flank the Ter macrodomain. RRS sequence and positioning is highly conserved in enterobacteria. At least one RRS appears to be essential unless a 200 kb region encompassing one of them is amplified. The RRS contain a G-quadruplex on the lagging strand which impedes DNA polymerase extension producing lagging strand ssDNA gaps, $ \le$2000 bp long, upon replisome passage. Deletion of both RRS elements has substantial effects on global genome structure and topology. We hypothesize that RRS elements serve as topological relief valves during chromosome replication and segregation. There have been no screens for genomic sequences that trigger transient gap formation. Functional analogs of RRS could be widespread, possibly including some enigmatic G-quadruplexes in eukaryotes.

i-Motif DNA: identification, formation, and cellular functions

An i-motif (iM) is a four-stranded (quadruplex) DNA structure that folds from cytosine (C)-rich sequences. iMs can fold under many different conditions in vitro, which paves the way for their formation in living cells. iMs are thought to play key roles in various DNA transactions, notably in the regulation of genome stability, gene transcription, mRNA translation, DNA replication, telomere and centromere functions, and human diseases. We summarize the different techniques used to assess the folding of iMs in vitro and provide an overview of the internal and external factors that affect their formation and stability in vivo. We describe the possible biological relevance of iMs and propose directions towards their use as target in biology.

YY1 is involved in homologous recombination inhibition at guanine quadruplex sites in human cells

Homologous recombination (HR) is a key process for repairing DNA double strand breaks and for promoting genetic diversity. However, HR occurs unevenly across the genome, and certain genomic features can influence its activity. One such feature is the presence of guanine quadruplexes (G4s), stable secondary structures widely distributed throughout the genome. These G4s play essential roles in gene transcription and genome stability regulation. Especially, elevated G4 levels in cells deficient in the Bloom syndrome helicase (BLM) significantly enhance HR at G4 sites, potentially threatening genome stability. Here, we investigated the role of G4-binding protein Yin Yang-1 (YY1) in modulating HR at G4 sites in human cells. Our results show that YY1's binding to G4 structures suppresses sister chromatid exchange after BLM knockdown, and YY1's chromatin occupancy negatively correlates with the overall HR rate observed across the genome. By limiting RAD51 homolog 1 (RAD51) access, YY1 preferentially binds to essential genomic regions, shielding them from excessive HR. Our findings unveil a novel role of YY1-G4 interaction, revealing novel insights into cellular mechanisms involved in HR regulation.

The determinants of the rarity of nucleic and peptide short sequences in nature

The prevalence of nucleic and peptide short sequences across organismal genomes and proteomes has not been thoroughly investigated. We examined 45 785 reference genomes and 21 871 reference proteomes, spanning archaea, bacteria, eukaryotes and viruses to calculate the rarity of short sequences in them. To capture this, we developed a metric of the rarity of each sequence in nature, the rarity index. We find that the frequency of certain dipeptides in rare oligopeptide sequences is hundreds of times lower than expected, which is not the case for any dinucleotides. We also generate predictive regression models that infer the rarity of nucleic and proteomic sequences across nature or within each domain of life and viruses separately. When examining each of the three domains of life and viruses separately, the R² performance of the model predicting rarity for 5-mer peptides from mono- and dipeptides ranged between 0.814 and 0.932. A separate model predicting rarity for 10-mer oligonucleotides from mono- and dinucleotides achieved R² performance between 0.408 and 0.606. Our results indicate that the mono- and dinucleotide composition of nucleic sequences and the mono- and dipeptide composition of peptide sequences can explain a significant proportion of the variance in their frequencies in nature.

The complete sequence and comparative analysis of ape sex chromosomes

Apes possess two sex chromosomes-the male-specific Y chromosome and the X chromosome, which is present in both males and females. The Y chromosome is crucial for male reproduction, with deletions being linked to infertility1. The X chromosome is vital for reproduction and cognition2. Variation in mating patterns and brain function among apes suggests corresponding differences in their sex chromosomes. However, owing to their repetitive nature and incomplete reference assemblies, ape sex chromosomes have been challenging to study. Here, using the methodology developed for the telomere-to-telomere (T2T) human genome, we produced gapless assemblies of the X and Y chromosomes for five great apes (bonobo (Pan paniscus), chimpanzee (Pan troglodytes), western lowland gorilla (Gorilla gorilla gorilla), Bornean orangutan (Pongo pygmaeus) and Sumatran orangutan (Pongo abelii)) and a lesser ape (the siamang gibbon (Symphalangus syndactylus)), and untangled the intricacies of their evolution. Compared with the X chromosomes, the ape Y chromosomes vary greatly in size and have low alignability and high levels of structural rearrangements-owing to the accumulation of lineage-specific ampliconic regions, palindromes, transposable elements and satellites. Many Y chromosome genes expand in multi-copy families and some evolve under purifying selection. Thus, the Y chromosome exhibits dynamic evolution, whereas the X chromosome is more stable. Mapping short-read sequencing data to these assemblies revealed diversity and selection patterns on sex chromosomes of more than 100 individual great apes. These reference assemblies are expected to inform human evolution and conservation genetics of non-human apes, all of which are endangered species.

Chromatin damage generated by DNA intercalators leads to degradation of RNA Polymerase II

In cancer therapy, DNA intercalators are mainly known for their capacity to kill cells by inducing DNA damage. Recently, several DNA intercalators have attracted much interest given their ability to inhibit RNA Polymerase I transcription (BMH-21), evict histones (Aclarubicin) or induce chromatin trapping of FACT (Curaxin CBL0137). Interestingly, these DNA intercalators lack the capacity to induce DNA damage while still retaining cytotoxic effects and stabilize p53. Herein, we report that these DNA intercalators impact chromatin biology by interfering with the chromatin stability of RNA polymerases I, II and III. These three compounds have the capacity to induce degradation of RNA polymerase II and they simultaneously enable the trapping of Topoisomerases TOP2A and TOP2B on the chromatin. In addition, BMH-21 also acts as a catalytic inhibitor of Topoisomerase II, resembling Aclarubicin. Moreover, BMH-21 induces chromatin trapping of the histone chaperone FACT and propels accumulation of Z-DNA and histone eviction, similarly to Aclarubicin and CBL0137. These DNA intercalators have a cumulative impact on general transcription machinery by inducing accumulation of topological defects and impacting nuclear chromatin. Therefore, their cytotoxic capabilities may be the result of compounding deleterious effects on chromatin homeostasis.

Pathogenic CANVAS (AAGGG)n repeats stall DNA replication due to the formation of alternative DNA structures

CANVAS is a recently characterized repeat expansion disease, most commonly caused by homozygous expansions of an intronic (A2G3)n repeat in the RFC1 gene. There are a multitude of repeat motifs found in the human population at this locus, some of which are pathogenic and others benign. In this study, we conducted structure-functional analyses of the pathogenic (A2G3)n and nonpathogenic (A4G)n repeats. We found that the pathogenic, but not the nonpathogenic, repeat presents a potent, orientation-dependent impediment to DNA polymerization in vitro. The pattern of the polymerization blockage is consistent with triplex or quadruplex formation in the presence of magnesium or potassium ions, respectively. Chemical probing of both repeats in vitro reveals triplex H-DNA formation by only the pathogenic repeat. Consistently, bioinformatic analysis of S1-END-seq data from human cell lines shows preferential H-DNA formation genome-wide by (A2G3)n motifs over (A4G)n motifs. Finally, the pathogenic, but not the nonpathogenic, repeat stalls replication fork progression in yeast and human cells. We hypothesize that the CANVAS-causing (A2G3)n repeat represents a challenge to genome stability by folding into alternative DNA structures that stall DNA replication.

Synthesis and antiproliferative evaluation of novel 3,5,8-trisubstituted coumarins against breast cancer

Aim: This study focused on designing and synthesizing novel derivatives of 3,5,8-trisubstituted coumarin. Results: The synthesized compounds, particularly compound 5, exhibited significant cytotoxic effects on MCF-7 cells, surpassing staurosporine, and reduced toxicity toward MCF-10A cells, highlighting potential pharmacological advantages. Further, compound 5 altered the cell cycle and significantly increased apoptosis in MCF-7 cells, involving both early (41.7-fold) and late stages (33-fold), while moderately affecting necrotic signaling. The antitumor activity was linked to a notable reduction (4.78-fold) in topoisomerase IIβ expression. Molecular modeling indicated compound 5's strong affinity for EGFR, human EGF2 and topoisomerase II proteins. Conclusion: These findings highlight compound 5 as a multifaceted antitumor agent for breast cancer.

Controlling Genome Topology with Sequences that Trigger Post-replication Gap Formation During Replisome Passage: The E. coli RRS Elements

We report that the Escherichia coli chromosome includes novel GC-rich genomic structural elements that trigger formation of post-replication gaps upon replisome passage. The two nearly perfect 222 bp repeats, designated Replication Risk Sequences or RRS, are each 650 kb from the terminus sequence dif and flank the Ter macrodomain. RRS sequence and positioning is highly conserved in enterobacteria. At least one RRS appears to be essential unless a 200 kb region encompassing one of them is amplified. The RRS contain a G-quadruplex on the lagging strand which impedes DNA polymerase extension producing lagging strand ssDNA gaps, ≤2000 bp long, upon replisome passage. Deletion of both RRS elements has substantial effects on global genome structure and topology. We hypothesize that RRS elements serve as topological relief valves during chromosome replication and segregation. There have been no screens for genomic sequences that trigger transient gap formation. Functional analogs of RRS could be widespread, possibly including some enigmatic G-quadruplexes in eukaryotes.

Oxidative DNA damage on the VEGF G-quadruplex forming promoter is repaired via long-patch BER

In response to oxidative damage, base excision repair (BER) enzymes perturb the structural equilibrium of the VEGF promoter between B-form and G4 DNA conformations, resulting in epigenetic-like modifications of gene expression. However, the mechanistic details remain enigmatic, including the activity and coordination of BER enzymes on the damaged G4 promoter. To address this, we investigated the ability of each BER factor to conduct its repair activity on VEGF promoter G4 DNA substrates by employing pre-steady-state kinetics assays and in vitro coupled BER assays. OGG1 was able to initiate BER on double-stranded VEGF promoter G4 DNA substrates. Moreover, pre-steady-state kinetics revealed that compared to B-form DNA, APE1 repair activity on the G4 was decreased ~two-fold and is the result of slower product release as opposed to inefficient strand cleavage. Interestingly, Pol β performs multiple insertions on G4 substrates via strand displacement DNA synthesis in contrast to a single insertion on B-form DNA. The multiple insertions inhibit ligation of the Pol β products, and hence BER is not completed on the VEGF G4 promoter substrates through canonical short-patch BER. Instead, repair requires the long-patch BER flap-endonuclease activity of FEN1 in response to the multiple insertions by Pol β prior to ligation. Because the BER proteins and their repair activities are a key part of the VEGF transcriptional enhancement in response to oxidative DNA damage of the G4 VEGF promoter, the new insights reported here on BER activity in the context of this promoter are relevant toward understanding the mechanism of transcriptional regulation.

Structural investigation of pathogenic RFC1 AAGGG pentanucleotide repeats reveals a role of G-quadruplex in dysregulated gene expression in CANVAS

An expansion of AAGGG pentanucleotide repeats in the replication factor C subunit 1 (RFC1) gene is the genetic cause of cerebellar ataxia, neuropathy, and vestibular areflexia syndrome (CANVAS), and it also links to several other neurodegenerative diseases including the Parkinson's disease. However, the pathogenic mechanism of RFC1 AAGGG repeat expansion remains enigmatic. Here, we report that the pathogenic RFC1 AAGGG repeats form DNA and RNA parallel G-quadruplex (G4) structures that play a role in impairing biological processes. We determine the first high-resolution nuclear magnetic resonance (NMR) structure of a bimolecular parallel G4 formed by d(AAGGG)2AA and reveal how AAGGG repeats fold into a higher-order structure composed of three G-tetrad layers, and further demonstrate the formation of intramolecular G4s in longer DNA and RNA repeats. The pathogenic AAGGG repeats, but not the nonpathogenic AAAAG repeats, form G4 structures to stall DNA replication and reduce gene expression via impairing the translation process in a repeat-length-dependent manner. Our results provide an unprecedented structural basis for understanding the pathogenic mechanism of AAGGG repeat expansion associated with CANVAS. In addition, the high-resolution structures resolved in this study will facilitate rational design of small-molecule ligands and helicases targeting G4s formed by AAGGG repeats for therapeutic interventions.

G-quadruplex-mediated genomic instability drives SNVs in cancer

G-quadruplex (G4s) DNA structures have been implicated in inducing genomic instability and contributing to cancer development. However, the relationship between G4s and cancer-related single nucleotide variants (cSNVs) in clinical settings remains unclear. In this large-scale study, we integrated experimentally validated G4s with genomic cSNVs from 13480 cancer patients to investigate the spatial association of G4s with the cellular cSNV landscape. Our findings demonstrate an increase in local genomic instability with increasing local G4 content in cancer patients, suggesting a potential role for G4s in driving cSNVs. Notably, we observed distinct spatial patterns of cSNVs and common single nucleotide variants (dbSNVs) in relation to G4s, implying different mechanisms for their generation and accumulation. We further demonstrate large, cancer-specific differences in the relationship of G4s and cSNVs, which could have important implications for a new class of G4-stabilizing cancer therapeutics. Moreover, we show that high G4-content can serve as a prognostic marker for local cSNV density and patient survival rates. Our findings underscore the importance of considering G4s in cancer research and highlight the need for further investigation into the underlying molecular mechanisms of G4-mediated genomic instability, especially in the context of cancer.

Dimeric structures of DNA ATTTC repeats promoted by divalent cations

Structural studies of repetitive DNA sequences may provide insights why and how certain repeat instabilities in their number and nucleotide sequence are managed or even required for normal cell physiology, while genomic variability associated with repeat expansions may also be disease-causing. The pentanucleotide ATTTC repeats occur in hundreds of genes important for various cellular processes, while their insertion and expansion in noncoding regions are associated with neurodegeneration, particularly with subtypes of spinocerebellar ataxia and familial adult myoclonic epilepsy. We describe a new striking domain-swapped DNA-DNA interaction triggered by the addition of divalent cations, including Mg2+ and Ca2+. The results of NMR characterization of d(ATTTC)3 in solution show that the oligonucleotide folds into a novel 3D architecture with two central C:C+ base pairs sandwiched between a couple of T:T base pairs. This structural element, referred to here as the TCCTzip, is characterized by intercalative hydrogen-bonding, while the nucleobase moieties are poorly stacked. The 5'- and 3'-ends of TCCTzip motif are connected by stem-loop segments characterized by A:T base pairs and stacking interactions. Insights embodied in the non-canonical DNA structure are expected to advance our understanding of why only certain pyrimidine-rich DNA repeats appear to be pathogenic, while others can occur in the human genome without any harmful consequences.

Distinguishing G-Quadruplexes Stabilizer and Chaperone for c-MYC Promoter G-Quadruplexes through Single-Molecule Manipulation

G-quadruplex (G4) selective stabilizing ligands can regulate c-MYC gene expression, but the kinetic basis remains unclear. Determining the effects of ligands on c-MYC promoter G4s' folding/unfolding kinetics is challenging due to the polymorphic nature of G4s and the high energy barrier to unfold c-MYC promoter G4s. Here, we used single-molecule magnetic tweezers to manipulate a duplex hairpin containing a c-MYC promoter sequence to mimic the transiently denatured duplex during transcription. We measured the effects of six commonly used G4s binding ligands on the competition between quadruplex and duplex structures, as well as the folding/unfolding kinetics of G4s. Our results revealed two distinct roles for G4s selective stabilization: CX-5461 is mainly acting as c-MYC G4s stabilizer, reducing the unfolding rate (ku) of c-MYC G4s, whereas PDS and 360A also act as G4s chaperone, accelerating the folding rates (kf) of c-MYC G4s. qRT-PCR results obtained from CA46 and Raji cell lines demonstrated that G4s stabilizing ligands can downregulate c-MYC expression, while G4s stabilizer CX-5461 exhibited the strongest c-MYC gene suppression. These results shed light on the potential of manipulating G4s' folding/unfolding kinetics by ligands for precise regulation of promoter G4-associated biological activities.

mEnrich-seq: methylation-guided enrichment sequencing of bacterial taxa of interest from microbiome

Metagenomics has enabled the comprehensive study of microbiomes. However, many applications would benefit from a method that sequences specific bacterial taxa of interest, but not most background taxa. We developed mEnrich-seq (in which 'm' stands for methylation and seq for sequencing) for enriching taxa of interest from metagenomic DNA before sequencing. The core idea is to exploit the self versus nonself differentiation by natural bacterial DNA methylation and rationally choose methylation-sensitive restriction enzymes, individually or in combination, to deplete host and background taxa while enriching targeted taxa. This idea is integrated with library preparation procedures and applied in several applications to enrich (up to 117-fold) pathogenic or beneficial bacteria from human urine and fecal samples, including species that are hard to culture or of low abundance. We assessed 4,601 bacterial strains with mapped methylomes so far and showed broad applicability of mEnrich-seq. mEnrich-seq provides microbiome researchers with a versatile and cost-effective approach for selective sequencing of diverse taxa of interest.

DNA at the center of mammalian innate immune recognition of bacterial biofilms

Historically, the study of innate immune detection of bacterial infections has focused on the recognition of pathogen-associated molecular patterns (PAMPs) from bacteria growing as single cells in planktonic phase. However, over the past two decades, studies have highlighted an adaptive advantage of bacteria: the formation of biofilms. These structures are complex fortresses that stand against a hostile environment, including antibiotics and immune responses. Extracellular DNA (eDNA) is a crucial component of the matrix of most known biofilms. In this opinion article, I propose that eDNA is a universal PAMP that the immune system uses to recognize biofilms. Outstanding questions concern the discrimination between biofilm-associated eDNA and DNA from planktonic bacteria, the innate receptors involved, and the immune response to biofilms.

In vivo detection of DNA secondary structures using permanganate/S1 footprinting with direct adapter ligation and sequencing (PDAL-Seq)

DNA secondary structures are essential elements of the genomic landscape, playing a critical role in regulating various cellular processes. These structures refer to G-quadruplexes, cruciforms, Z-DNA or H-DNA structures, amongst others (collectively called 'non-B DNA'), which DNA molecules can adopt beyond the B conformation. DNA secondary structures have significant biological roles, and their landscape is dynamic and can rearrange due to various factors, including changes in cellular conditions, temperature, and DNA-binding proteins. Understanding this dynamic nature is crucial for unraveling their functions in cellular processes. Detecting DNA secondary structures remains a challenge. Conventional methods, such as gel electrophoresis and chemical probing, have limitations in terms of sensitivity and specificity. Emerging techniques, including next-generation sequencing and single-molecule approaches, offer promise but face challenges since these techniques are mostly limited to only one type of secondary structure. Here we describe an updated version of a technique permanganate/S1 nuclease footprinting, which uses potassium permanganate to trap single-stranded DNA regions as found in many non-B structures, in combination with S1 nuclease digest and adapter ligation to detect genome-wide non-B formation. To overcome technical hurdles, we combined this method with direct adapter ligation and sequencing (PDAL-Seq). Furthermore, we established a user-friendly pipeline available on Galaxy to standardize PDAL-Seq data analysis. This optimized method allows the analysis of many types of DNA secondary structures that form in a living cell and will advance our knowledge of their roles in health and disease.

MoCoLo: a testing framework for motif co-localization

Sequence-level data offers insights into biological processes through the interaction of two or more genomic features from the same or different molecular data types. Within motifs, this interaction is often explored via the co-occurrence of feature genomic tracks using fixed-segments or analytical tests that respectively require window size determination and risk of false positives from over-simplified models. Moreover, methods for robustly examining the co-localization of genomic features, and thereby understanding their spatial interaction, have been elusive. We present a new analytical method for examining feature interaction by introducing the notion of reciprocal co-occurrence, define statistics to estimate it and hypotheses to test for it. Our approach leverages conditional motif co-occurrence events between features to infer their co-localization. Using reverse conditional probabilities and introducing a novel simulation approach that retains motif properties (e.g. length, guanine-content), our method further accounts for potential confounders in testing. As a proof-of-concept, motif co-localization (MoCoLo) confirmed the co-occurrence of histone markers in a breast cancer cell line. As a novel analysis, MoCoLo identified significant co-localization of oxidative DNA damage within non-B DNA-forming regions that significantly differed between non-B DNA structures. Altogether, these findings demonstrate the potential utility of MoCoLo for testing spatial interactions between genomic features via their co-localization.

DNA superhelicity

Closing each strand of a DNA duplex upon itself fixes its linking number L. This topological condition couples together the secondary and tertiary structures of the resulting ccDNA topoisomer, a constraint that is not present in otherwise identical nicked or linear DNAs. Fixing L has a range of structural, energetic and functional consequences. Here we consider how L having different integer values (that is, different superhelicities) affects ccDNA molecules. The approaches used are primarily theoretical, and are developed from a historical perspective. In brief, processes that either relax or increase superhelicity, or repartition what is there, may either release or require free energy. The energies involved can be substantial, sufficient to influence many events, directly or indirectly. Here two examples are developed. The changes of unconstrained superhelicity that occur during nucleosome attachment and release are examined. And a simple theoretical model of superhelically driven DNA structural transitions is described that calculates equilibrium distributions for populations of identical topoisomers. This model is used to examine how these distributions change with superhelicity and other factors, and applied to analyze several situations of biological interest.

Purification and biochemical characterization of the G4 resolvase and DNA helicase FANCJ

G-quadruplex (G4) DNA or RNA poses a unique nucleic acid structure in genomic transactions. Because of the unique topology presented by G4, cells have exquisite mechanisms and pathways to metabolize G4 that arise in guanine-rich regions of the genome such as telomeres, promoter regions, ribosomal DNA, and other chromosomal elements. G4 resolvases are often represented by a class of molecular motors known as helicases that disrupt the Hoogsteen hydrogen bonds in G4 by harnessing the chemical energy of nucleoside triphosphate hydrolysis. Of special interest to researchers in the field, including us, is the human FANCJ DNA helicase that efficiently resolves G4 DNA structures. Notably, FANCJ mutations are linked to Fanconi Anemia and are prominent in breast and ovarian cancer. Since our discovery that FANCJ efficiently resolves G4 DNA structures 15 years ago, we and other labs have characterized mechanistic aspects of FANCJ-catalyzed G4 resolution and its biological importance in genomic integrity and cellular DNA replication. In addition to its G4 resolvase function, FANCJ is also a classic DNA helicase that acts on conventional duplex DNA structures, which are relevant to the enzyme's role in interstrand cross link repair, double-strand break repair via homologous recombination, and response to replication stress. Here, we describe detailed procedures for the purification of recombinant FANCJ protein and characterization of its G4 resolvase and duplex DNA helicase activity.

EndoQuad: a comprehensive genome-wide experimentally validated endogenous G-quadruplex database

G-quadruplexes (G4s) are non-canonical four-stranded structures and are emerging as novel genetic regulatory elements. However, a comprehensive genomic annotation of endogenous G4s (eG4s) and systematic characterization of their regulatory network are still lacking, posing major challenges for eG4 research. Here, we present EndoQuad (https://EndoQuad.chenzxlab.cn/) to address these pressing issues by integrating high-throughput experimental data. First, based on high-quality genome-wide eG4s mapping datasets (human: 1181; mouse: 24; chicken: 2) generated by G4 ChIP-seq/CUT&Tag, we generate a reference set of genome-wide eG4s. Our multi-omics analyses show that most eG4s are identified in one or a few cell types. The eG4s with higher occurrences across samples are more structurally stable, evolutionarily conserved, enriched in promoter regions, mark highly expressed genes and associate with complex regulatory programs, demonstrating higher confidence level for further experiments. Finally, we integrate millions of functional genomic variants and prioritize eG4s with regulatory functions in disease and cancer contexts. These efforts have culminated in the comprehensive and interactive database of experimentally validated DNA eG4s. As such, EndoQuad enables users to easily access, download and repurpose these data for their own research. EndoQuad will become a one-stop resource for eG4 research and lay the foundation for future functional studies.

Synchrotron-based FTIR microspectroscopy reveals DNA methylation profile in DNA-HALO structure

Fourier transform infrared (FTIR) spectroscopy is a rapid, non-destructive and label-free technique for identifying subtle changes in all bio-macromolecules, and has been used as a method of choice for studying DNA conformation, secondary DNA structure transition and DNA damage. In addition, the specific level of chromatin complexity is introduced via epigenetic modifications forcing the technological upgrade in the analysis of such an intricacy. As the most studied epigenetic mechanism, DNA methylation is a major regulator of transcriptional activity, involved in the suppression of a broad spectrum of genes and its deregulation is involved in all non-communicable diseases. The present study was designed to explore the use of synchrotron-based FTIR analysis to monitor the subtle changes in molecule bases regarding the DNA methylation status of cytosine in the whole genome. In order to reveal the conformation-related best sample for FTIR-based DNA methylation analysis in situ, we used methodology for nuclear HALO preparations and slightly modified it to isolated DNA in HALO formations. Nuclear DNA-HALOs represent samples with preserved higher-order chromatin structure liberated of any protein residues that are closer to native DNA conformation than genomic DNA (gDNA) isolated by the standard batch procedure. Using FTIR spectroscopy we analyzed the DNA methylation profile of isolated gDNA and compared it with the DNA-HALOs. This study demonstrated the potential of FTIR microspectroscopy to detect DNA methylation marks in analyzed DNA-HALO specimens more precisely in comparison with classical DNA extraction procedures that yield unstructured whole genomic DNA. In addition, we used different cell types to assess their global DNA methylation profile, as well as defined specific infrared peaks that can be used for screening DNA methylation.

A beginner's handbook to identify and characterize i-motif DNA

Genomic DNA exhibits an innate ability to manifest diverse sequence-dependent secondary structures, serving crucial functions in gene regulation and cellular equilibrium. While extensive research has confirmed the formation of G-quadruplex structures by guanine-rich sequences in vitro and in cells, recent investigations have turned the quadruplex community's attention to the cytosine (C)-rich complementary strands that can adopt unique tetra-stranded conformation, termed as intercalated motif or i-motif. I-motifs are stabilized by hemi-protonated C:CH+ base pairs under acidic conditions. Initially, the in vivo occurrence of i-motifs was underestimated because their formation is favored at non-physiological pH. However, groundbreaking research utilizing the structure-specific iMab antibody and high-throughput sequencing have recently detected their conserved dispersion throughout the genome, challenging previous assumptions. Given the evolving nature of this research field, it becomes imperative to conduct independent in vitro experiments aimed at identifying potential i-motif formation in C-rich sequences and consolidating the findings to address the properties of i-motifs. This chapter serves as an introductory guide for the swift identification of novel i-motifs, where we present an experimental framework for investigating and characterizing i-motif sequences in vitro. In this chapter, we selected a synthetic oligonucleotide (C7T3) sequence and outlined appropriate methodologies for annealing the i-motif structure into suitable buffers. Then, we validated its formation by CD (Circular Dichroism) and NMR (Nuclear Magnetic Resonance) spectroscopy. Finally, we provided a thorough account of the step-by-step procedures to investigate the effect of i-motif formation on the stalling or retardation of DNA replication using high resolution primer extension assays.

The Complete Sequence and Comparative Analysis of Ape Sex Chromosomes

Apes possess two sex chromosomes-the male-specific Y and the X shared by males and females. The Y chromosome is crucial for male reproduction, with deletions linked to infertility. The X chromosome carries genes vital for reproduction and cognition. Variation in mating patterns and brain function among great apes suggests corresponding differences in their sex chromosome structure and evolution. However, due to their highly repetitive nature and incomplete reference assemblies, ape sex chromosomes have been challenging to study. Here, using the state-of-the-art experimental and computational methods developed for the telomere-to-telomere (T2T) human genome, we produced gapless, complete assemblies of the X and Y chromosomes for five great apes (chimpanzee, bonobo, gorilla, Bornean and Sumatran orangutans) and a lesser ape, the siamang gibbon. These assemblies completely resolved ampliconic, palindromic, and satellite sequences, including the entire centromeres, allowing us to untangle the intricacies of ape sex chromosome evolution. We found that, compared to the X, ape Y chromosomes vary greatly in size and have low alignability and high levels of structural rearrangements. This divergence on the Y arises from the accumulation of lineage-specific ampliconic regions and palindromes (which are shared more broadly among species on the X) and from the abundance of transposable elements and satellites (which have a lower representation on the X). Our analysis of Y chromosome genes revealed lineage-specific expansions of multi-copy gene families and signatures of purifying selection. In summary, the Y exhibits dynamic evolution, while the X is more stable. Finally, mapping short-read sequencing data from >100 great ape individuals revealed the patterns of diversity and selection on their sex chromosomes, demonstrating the utility of these reference assemblies for studies of great ape evolution. These complete sex chromosome assemblies are expected to further inform conservation genetics of nonhuman apes, all of which are endangered species.

Non-canonical DNA structures in the human ribosomal DNA

Non-canonical structures (NCS) refer to the various forms of DNA that differ from the B-conformation described by Watson and Crick. It has been found that these structures are usual components of the genome, actively participating in its essential functions. The present review is focused on the nine kinds of NCS appearing or likely to appear in human ribosomal DNA (rDNA): supercoiling structures, R-loops, G-quadruplexes, i-motifs, DNA triplexes, cruciform structures, DNA bubbles, and A and Z DNA conformations. We discuss the conditions of their generation, including their sequence specificity, distribution within the locus, dynamics, and beneficial and detrimental role in the cell.

Sensing mitochondrial DNA stress in cardiotoxicity

Cytoplasmic mitochondrial DNA (mtDNA) can trigger the interferon response to promote disease progression, but mtDNA sensing mechanisms remain elusive. Lei et al. have shown that Z-DNA binding protein1 (ZBP1) cooperates with cyclic GMP-AMP synthase (cGAS) to sense Z-form mtDNA and transmit mtDNA stress signals to promote diseases such as cardiotoxicity, providing an important piece of the mtDNA stress landscape.

Epstein-Barr virus and host cell 3D genome organization

The human genome is organized in an extremely complexed yet ordered way within the nucleus. Genome organization plays a critical role in the regulation of gene expression. Viruses manipulate the host machinery to influence host genome organization to favor their survival and promote disease development. Epstein-Barr virus (EBV) is a common human virus, whose infection is associated with various diseases, including infectious mononucleosis, cancer, and autoimmune disorders. This review summarizes our current knowledge of how EBV uses different strategies to control the cellular 3D genome organization to affect cell gene expression to transform normal cells into lymphoblasts.

HLTF Prevents G4 Accumulation and Promotes G4-induced Fork Slowing to Maintain Genome Stability

G-quadruplexes (G4s) form throughout the genome and influence important cellular processes, but their deregulation can challenge DNA replication fork progression and threaten genome stability. Here, we demonstrate an unexpected, dual role for the dsDNA translocase HLTF in G4 metabolism. First, we find that HLTF is enriched at G4s in the human genome and suppresses G4 accumulation throughout the cell cycle using its ATPase activity. This function of HLTF affects telomere maintenance by restricting alternative lengthening of telomeres, a process stimulated by G4s. We also show that HLTF and MSH2, a mismatch repair factor that binds G4s, act in independent pathways to suppress G4s and to promote resistance to G4 stabilization. In a second, distinct role, HLTF restrains DNA synthesis upon G4 stabilization by suppressing PrimPol-dependent repriming. Together, the dual functions of HLTF in the G4 response prevent DNA damage and potentially mutagenic replication to safeguard genome stability.

Detection of alternative DNA structures and its implications for human disease

Around 3% of the genome consists of simple DNA repeats that are prone to forming alternative (non-B) DNA structures, such as hairpins, cruciforms, triplexes (H-DNA), four-stranded guanine quadruplexes (G4-DNA), and others, as well as composite RNA:DNA structures (e.g., R-loops, G-loops, and H-loops). These DNA structures are dynamic and favored by the unwinding of duplex DNA. For many years, the association of alternative DNA structures with genome function was limited by the lack of methods to detect them in vivo. Here, we review the recent advancements in the field and present state-of-the-art technologies and methods to study alternative DNA structures. We discuss the limitations of these methods as well as how they are beginning to provide insights into causal relationships between alternative DNA structures, genome function and stability, and human disease.

Yeast transcription factor Msn2 binds to G4 DNA

Sequences capable of forming quadruplex or G4 DNA are prevalent in the promoter regions. The transformation from canonical to non-canonical secondary structure apparently regulates transcription of a number of human genes. In the budding yeast Saccharomyces cerevisiae, we identified 37 genes with a G4 motif in the promoters including 20 genes that contain stress response element (STRE) overlapping a G4 motif. STRE is the binding site of stress response regulators Msn2 and Msn4, transcription factors belonging to the C2H2 zinc-finger protein family. We show here that Msn2 binds directly to the G4 DNA structure through its zinc-finger domain with a dissociation constant similar to that of STRE-binding and that, in a stress condition, Msn2 is enriched at G4 DNA-forming loci in the yeast genome. For a large fraction of genes with G4/STRE-containing promoters, treating with G4-ligands led to significant elevations in transcription levels. Such transcriptional elevation was greatly diminished in a msn2Δ msn4Δ background and was partly muted when the G4 motif was disrupted. Taken together, our data suggest that G4 DNA could be an alternative binding site of Msn2 in addition to STRE, and that G4 DNA formation could be an important element of transcriptional regulation in yeast.

Folding Double-Stranded DNA into Designed Shapes with Triplex-Forming Oligonucleotides

The compaction and organization of genomic DNA is a central mechanism in eukaryotic cells, but engineered architectural control over double-stranded DNA (dsDNA) is notably challenging. Here, long dsDNA templates are folded into designed shapes via triplex-mediated self-assembly. Triplex-forming oligonucleotides (TFOs) bind purines in dsDNA via normal or reverse Hoogsteen interactions. In the triplex origami methodology, these non-canonical interactions are programmed to compact dsDNA (linear or plasmid) into well-defined objects, which demonstrate a variety of structural features: hollow and raster-filled, single- and multi-layered, with custom curvatures and geometries, and featuring lattice-free, square-, or honeycomb-pleated internal arrangements. Surprisingly, the length of integrated and free-standing dsDNA loops can be modulated with near-perfect efficiency; from hundreds down to only six bp (2 nm). The inherent rigidity of dsDNA promotes structural robustness and non-periodic structures of almost 25.000 nt are therefore formed with fewer unique starting materials, compared to other DNA-based self-assembly methods. Densely triplexed structures also resist degradation by DNase I. Triplex-mediated dsDNA folding is methodologically straightforward and orthogonal to Watson-Crick-based methods. Moreover, it enables unprecedented spatial control over dsDNA templates.

Genomic Instability of G-Quadruplex Sequences in Escherichia coli: Roles of DinG, RecG, and RecQ Helicases

Guanine-rich DNA can fold into highly stable four-stranded DNA structures called G-quadruplexes (G4). Originally identified in sequences from telomeres and oncogene promoters, they can alter DNA metabolism. Indeed, G4-forming sequences represent obstacles for the DNA polymerase, with important consequences for cell life as they may lead to genomic instability. To understand their role in bacterial genomic instability, different G-quadruplex-forming repeats were cloned into an Escherichia coli genetic system that reports frameshifts and complete or partial deletions of the repeat when the G-tract comprises either the leading or lagging template strand during replication. These repeats formed stable G-quadruplexes in single-stranded DNA but not naturally supercoiled double-stranded DNA. Nevertheless, transcription promoted G-quadruplex formation in the resulting R-loop for (G3T)4 and (G3T)8 repeats. Depending on genetic background and sequence propensity for structure formation, mutation rates varied by five orders of magnitude. Furthermore, while in vitro approaches have shown that bacterial helicases can resolve G4, it is still unclear whether G4 unwinding is important in vivo. Here, we show that a mutation in recG decreased mutation rates, while deficiencies in the structure-specific helicases DinG and RecQ increased mutation rates. These results suggest that G-quadruplex formation promotes genetic instability in bacteria and that helicases play an important role in controlling this process in vivo.

Pathogenic CANVAS (AAGGG)n repeats stall DNA replication due to the formation of alternative DNA structures

CANVAS is a recently characterized repeat expansion disease, most commonly caused by homozygous expansions of an intronic (A2G3)n repeat in the RFC1 gene. There are a multitude of repeat motifs found in the human population at this locus, some of which are pathogenic and others benign. In this study, we conducted structure-functional analyses of the main pathogenic (A2G3)n and the main nonpathogenic (A4G)n repeats. We found that the pathogenic, but not the nonpathogenic, repeat presents a potent, orientation-dependent impediment to DNA polymerization in vitro. The pattern of the polymerization blockage is consistent with triplex or quadruplex formation in the presence of magnesium or potassium ions, respectively. Chemical probing of both repeats in supercoiled DNA reveals triplex H-DNA formation by the pathogenic repeat. Consistently, bioinformatic analysis of the S1-END-seq data from human cell lines shows preferential H-DNA formation genome-wide by (A2G3)n motifs over (A4G)n motifs in vivo. Finally, the pathogenic, but not the non-pathogenic, repeat stalls replication fork progression in yeast and human cells. We hypothesize that CANVAS-causing (A2G3)n repeat represents a challenge to genome stability by folding into alternative DNA structures that stall DNA replication.

G-quadruplexes from non-coding RNAs

Non-coding RNAs (ncRNAs) are significant regulators of gene expression in a wide range of biological processes, such as transcription, RNA maturation, or translation. ncRNAs interplay with proteins or other RNAs through not only classical sequence-based mechanisms but also unique higher-order structures such as RNA G-quadruplexes (rG4s). rG4s are predictably formed in guanine-rich sequences and are closely related to various human diseases, such as tumors, neurodegenerative diseases, and infections. This review focuses on the vital role of rG4s in ncRNAs, particularly lncRNAs and miRNAs. We outline the dynamic balance between rG4s and RNA stem-loop/hairpin structures and the interplay between ncRNAs and interactors, thereby modulating gene expression and disease progression. A complete understanding of the biological regulatory role and mechanism of rG4s in ncRNAs affirms the critical importance of folding into the appropriate three-dimensional structure in maintaining or modulating the functions of ncRNAs. It makes them novel therapeutic targets for adjusting potential-G4-containing-ncRNAs-associated diseases.

Chromosomal fragile site breakage by EBV-encoded EBNA1 at clustered repeats

Epstein-Barr virus (EBV) is an oncogenic herpesvirus associated with several cancers of lymphocytic and epithelial origin1-3. EBV encodes EBNA1, which binds to a cluster of 20 copies of an 18-base-pair palindromic sequence in the EBV genome4-6. EBNA1 also associates with host chromosomes at non-sequence-specific sites7, thereby enabling viral persistence. Here we show that the sequence-specific DNA-binding domain of EBNA1 binds to a cluster of tandemly repeated copies of an EBV-like, 18-base-pair imperfect palindromic sequence encompassing a region of about 21 kilobases at human chromosome 11q23. In situ visualization of the repetitive EBNA1-binding site reveals aberrant structures on mitotic chromosomes characteristic of inherently fragile DNA. We demonstrate that increasing levels of EBNA1 binding trigger dose-dependent breakage at 11q23, producing a fusogenic centromere-containing fragment and an acentric distal fragment, with both mis-segregated into micronuclei in the next cell cycles. In cells latently infected with EBV, elevating EBNA1 abundance by as little as twofold was sufficient to trigger breakage at 11q23. Examination of whole-genome sequencing of EBV-associated nasopharyngeal carcinomas revealed that structural variants are highly enriched on chromosome 11. Presence of EBV is also shown to be associated with an enrichment of chromosome 11 rearrangements across 2,439 tumours from 38 cancer types. Our results identify a previously unappreciated link between EBV and genomic instability, wherein EBNA1-induced breakage at 11q23 triggers acquisition of structural variations in chromosome 11.

Studies of the Mechanism of Nucleosome Dynamics: A Review on Multifactorial Regulation from Computational and Experimental Cases

The nucleosome, which organizes the long coil of genomic DNA in a highly condensed, polymeric way, is thought to be the basic unit of chromosomal structure. As the most important protein–DNA complex, its structural and dynamic features have been successively revealed in recent years. However, its regulatory mechanism, which is modulated by multiple factors, still requires systemic discussion. This study summarizes the regulatory factors of the nucleosome’s dynamic features from the perspective of histone modification, DNA methylation, and the nucleosome-interacting factors (transcription factors and nucleosome-remodeling proteins and cations) and focuses on the research exploring the molecular mechanism through both computational and experimental approaches. The regulatory factors that affect the dynamic features of nucleosomes are also discussed in detail, such as unwrapping, wrapping, sliding, and stacking. Due to the complexity of the high-order topological structures of nucleosomes and the comprehensive effects of regulatory factors, the research on the functional modulation mechanism of nucleosomes has encountered great challenges. The integration of computational and experimental approaches, the construction of physical modes for nucleosomes, and the application of deep learning techniques will provide promising opportunities for further exploration.

Template-Assembled Synthetic G-Quartets (TASQs): multiTASQing Molecular Tools for Investigating DNA and RNA G-Quadruplex Biology

Biomimetics is defined as a "practice of making technological design that copies natural processes", with the idea that "nature has already solved the challenges we are trying to solve" (Cambridge Dictionary). The challenge we decided to address several years ago was the selective targeting of G quadruplexes (G4s) by small molecules (G4 ligands). Why? Because G4s, which are four-stranded DNA and RNA structures that fold from guanine (G)-rich sequences, are suspected to play key biological roles in human cells and diseases. Selective G4 ligands can thus be used as small-molecule modulators to gain a deep understanding of cell circuitry where G4s are involved, thus complying with the very definition of chemical biology (Stuart Schreiber) applied here to G4 biology. How? Following a biomimetic approach that hinges on the observation that G4s are stable secondary structures owing to the ability of Gs to self-associate to form G quartets, and then of G quartets to self-stack to form the columnar core of G4s. Therefore, using a synthetic G quartet as a G4 ligand represents a unique example of biomimetic recognition of G4s.We formulated this hypothesis more than a decade ago, stepping on years of research on Gs, G4s, and G4 ligands. Our approach led to the design, synthesis, and use of a broad family of synthetic G quartets, also referred to as TASQs for template-assembled synthetic G quartets (John Sherman). This quest led us across various chemical lands (organic and supramolecular chemistry, chemical biology, and genetics), along a route on which every new generation of TASQ was a milestone in the growing portfolio of ever smarter molecular tools to decipher G4 biology. As discussed in this Account, we detail how and why we successively develop the very first prototypes of (i) biomimetic ligands, which interact with G4s according to a bioinspired, like-likes-like interaction between two G quartets, one from the ligand, the other from the G4; (ii) smart ligands, which adopt their active conformation only in the presence of their G4 targets; (iii) twice-as-smart ligands, which act as both smart ligands and smart fluorescent probes, whose fluorescence is triggered (turned on) upon interaction with their G4 targets; and (iv) multivalent ligands, which display additional functionalities enabling the detection, isolation, and identification of G4s both in vitro and in vivo. This quest led us to gather a panel of 14 molecular tools which were used to investigate the biology of G4s at a cellular level, from basic optical imaging to multiomics studies.

Impacts of Molecular Structure on Nucleic Acid-Protein Interactions

Interactions between nucleic acids and proteins are some of the most important interactions in biology because they are the cornerstones for fundamental biological processes, such as replication, transcription, and recombination [...].