Non-B DNA structure-induced genetic instability and evolution

Intragenomic mutational heterogeneity: structural and functional insights from gene evolution

Variation of mutation rates between species has been documented over decades, but the variation between different regions of a genome has been less often discussed. Recent studies using high-quality sequence data have revealed previously unknown levels of intragenomic heterogeneity of mutation rates and their association with other structural and functional features of DNA sequences. This article reviews accumulating evidence of this intragenomic heterogeneity and speculates its cause and influence on organismal phenotypes.

Non-canonical DNA in human and other ape telomere-to-telomere genomes

Non-canonical (non-B) DNA structures-e.g. bent DNA, hairpins, G-quadruplexes (G4s), Z-DNA, etc.-which form at certain sequence motifs (e.g. A-phased repeats, inverted repeats, etc.), have emerged as important regulators of cellular processes and drivers of genome evolution. Yet, they have been understudied due to their repetitive nature and potentially inaccurate sequences generated with short-read technologies. Here we comprehensively characterize such motifs in the long-read telomere-to-telomere (T2T) genomes of human, bonobo, chimpanzee, gorilla, Bornean orangutan, Sumatran orangutan, and siamang. Non-B DNA motifs are enriched at the genomic regions added to T2T assemblies and occupy 9%-15%, 9%-11%, and 12%-38% of autosomes and chromosomes X and Y, respectively. G4s and Z-DNA are enriched at promoters and enhancers, as well as at origins of replication. Repetitive sequences harbor more non-B DNA motifs than non-repetitive sequences, especially in the short arms of acrocentric chromosomes. Most centromeres and/or their flanking regions are enriched in at least one non-B DNA motif type, consistent with a potential role of non-B structures in determining centromeres. Our results highlight the uneven distribution of predicted non-B DNA structures across ape genomes and suggest their novel functions in previously inaccessible genomic regions.

DNA supercoiling-mediated G4/R-loop formation tunes transcription by controlling the access of RNA polymerase

RNA polymerase (RNAP) is a processive motor that modulates DNA supercoiling and reshapes DNA structures. The feedback loop between the DNA topology and transcription remains elusive. Here, we investigate the impact of potential G-quadruplex forming sequences (PQS) on transcription in response to DNA supercoiling. We find that supercoiled DNA increases transcription frequency 10-fold higher than relaxed DNA, which lead to an abrupt formation of G-quadruplex (G4) and R-loop structures. Moreover, the stable R-loop relieves topological strain, facilitated by G4 formation. The cooperative formation of G4/R-loop effectively alters the DNA topology around the promoter and suppresses transcriptional activity by impeding RNAP loading. These findings highlight negative supercoiling as a built-in spring that triggers a transcriptional burst followed by a rapid suppression upon G4/R-loop formation. This study sheds light on the intricate interplay between DNA topology and structural change in transcriptional regulation, with implications for understanding gene expression dynamics.

DNA2 and MSH2 activity collectively mediate chemically stabilized G4 for efficient telomere replication

G-quadruplexes (G4s) are widely existing stable DNA secondary structures in mammalian cells. A long-standing hypothesis is that timely resolution of G4s is needed for efficient and faithful DNA replication. In vitro , G4s may be unwound by helicases or alternatively resolved via DNA2 nuclease mediated G4 cleavage. However, little is known about the biological significance and regulatory mechanism of the DNA2-mediated G4 removal pathway. Here, we report that DNA2 deficiency or its chemical inhibition leads to a significant accumulation of G4s and stalled replication forks at telomeres, which is demonstrated by a high-resolution technology: Single molecular analysis of replicating DNA (SMARD). We further identify that the DNA repair complex MutSα (MSH2-MSH6) binds G4s and stimulates G4 resolution via DNA2-mediated G4 excision. MSH2 deficiency, like DNA2 deficiency or inhibition, causes G4 accumulation and defective telomere replication. Meanwhile, G4-stabilizing environmental compounds block G4 unwinding by helicases but not G4 cleavage by DNA2. Consequently, G4 stabilizers impair telomere replication and cause telomere instabilities, especially in cells deficient in DNA2 or MSH2.

Non-canonical DNA in human and other ape telomere-to-telomere genomes

Non-canonical (non-B) DNA structures-e.g., bent DNA, hairpins, G-quadruplexes (G4s), Z-DNA, etc.-which form at certain sequence motifs (e.g., A-phased repeats, inverted repeats, etc.), have emerged as important regulators of cellular processes and drivers of genome evolution. Yet, they have been understudied due to their repetitive nature and potentially inaccurate sequences generated with short-read technologies. Here we comprehensively characterize such motifs in the long-read telomere-to-telomere (T2T) genomes of human, bonobo, chimpanzee, gorilla, Bornean orangutan, Sumatran orangutan, and siamang. Non-B DNA motifs are enriched at the genomic regions added to T2T assemblies, and occupy 9-15%, 9-11%, and 12-38% of autosomes, and chromosomes X and Y, respectively. G4s and Z-DNA are enriched at promoters and enhancers, as well as at origins of replication. Repetitive sequences harbor more non-B DNA motifs than non-repetitive sequences, especially in the short arms of acrocentric chromosomes. Most centromeres and/or their flanking regions are enriched in at least one non-B DNA motif type, consistent with a potential role of non-B structures in determining centromeres. Our results highlight the uneven distribution of predicted non-B DNA structures across ape genomes and suggest their novel functions in previously inaccessible genomic regions.

Z-DNA at the crossroads: untangling its role in genome dynamics

DNA can fold into noncanonical left-handed Z-DNA conformation beyond the right-handed B-DNA. While its crystal structure was discovered nearly four decades ago, it was predominantly considered a structural curiosity. Recent evidence suggests that Z-DNA formation occurs in nuclear and mitochondrial DNA (mtDNA), with significant biological implications. However, our understanding of its roles remains in its infancy, primarily due to a lack of study tools. In this review we summarize the structure and function of Z-DNA within the genome while addressing the difficulties associated with identifying and investigating its role(s). We then critically evaluate several intracellular factors that can modulate and regulate Z-DNA. Additionally, we discuss the recent technological and methodological advances that may overcome the challenges and enhance our understanding of Z-DNA.

Ribosomal DNA arrays are the most H-DNA rich element in the human genome

Repetitive DNA sequences can form noncanonical structures such as H-DNA. The new telomere-to-telomere genome assembly for the human genome has eliminated gaps, enabling examination of highly repetitive regions including centromeric and pericentromeric repeats and ribosomal DNA arrays. We find that H-DNA appears once every 25 000 base pairs in the human genome. Its distribution is highly inhomogeneous with H-DNA motif hotspots being detectable in acrocentric chromosomes. Ribosomal DNA arrays are the genomic element with a 40.94-fold H-DNA enrichment. Across acrocentric chromosomes, we report that 54.82% of H-DNA motifs found in these chromosomes are in rDNA array loci. We discover that binding sites for the PRDM9-B allele, a variant of the PRDM9 protein, are enriched for H-DNA motifs. We further investigate these findings through an analysis of PRDM-9 ChIP-seq data across various PRDM-9 alleles, observing an enrichment of H-DNA motifs in the binding sites of A-like alleles (including A, B, and N alleles), but not C-like alleles (including C and L4 alleles). The enrichment of H-DNA motifs at ribosomal DNA arrays is consistent in nonhuman great ape genomes. We conclude that ribosomal DNA arrays are the most enriched genomic loci for H-DNA sequences in human and other great ape genomes.

Quest for space: Tenacity of DNA, Protein, and Lipid macromolecules in intracellular crowded environment

The biochemical processes in the cellular milieu involving biomacromolecular interaction usually occur in crowded and heterogeneous environments, impacting their structure, stability, and reactivity. The crowded environment in vivo is typically ignored for experimental investigations since the studies get complex due to intracellular biophysical interactions between nucleic acids, proteins, cellular membranes, and various cations/anions present in the cell. Thus, being a ubiquitous property of all cells, studying those biophysical aspects affecting biochemical processes under realistically crowded conditions is of prime importance. Crowders or crowding agents are usually exploited to mimic the in vivo conditions on interacting with such genomic species, revealing structural and functional changes resulting from excluded volume and soft interactions. In the last few years, studies including crowders of varied sizes have gained attention concerning the consequences of crowding agents on biomolecular structural transitions and stability. This review comprehensively summarizes macromolecular crowding, emphasizing the biophysical effects and contribution of soft interactions in the heterogeneous cellular environment.

Dissecting non-B DNA structural motifs in untranslated regions of eukaryotic genomes

The untranslated regions (UTRs) of genes significantly impact various biological processes, including transcription, posttranscriptional control, mRNA stability, localization, and translation efficiency. In functional areas of genomes, non-B DNA structures such as cruciform, curved, triplex, G-quadruplex, and Z-DNA structures are common and have an impact on cellular physiology. Although the role of these structures in cis-regulatory regions such as promoters is well established in eukaryotic genomes, their prevalence within UTRs across different eukaryotic classes has not been extensively documented. Our study investigated the prevalence of various non-B DNA motifs within the 5' and 3' UTRs across diverse eukaryotic species. Our comparative analysis encompassed the 5'-UTRs and 3'UTRs of 360 species representing diverse eukaryotic domains of life, including Arthropoda (Diptera, Hemiptera, and Hymenoptera), Chordata (Artiodactyla, Carnivora, Galliformes, Passeriformes, Primates, Rodentia, Squamata, Testudines), Magnoliophyta (Brassicales), Fabales (Poales), and Nematoda (Rhabditida), on the basis of datasets derived from the UTRdb. We observed that species belonging to taxonomic orders such as Rhabditida, Diptera, Brassicales, and Hemiptera present a prevalence of curved DNA motifs in their UTRs, whereas orders such as Testudines, Galliformes, and Rodentia present a preponderance of G-quadruplexes in both UTRs. The distribution of motifs is conserved across different taxonomic classes, although species-specific variations in motif preferences were also observed. Our research unequivocally illuminates the prevalence and potential functional implications of non-B DNA motifs, offering invaluable insights into the evolutionary and biological significance of these structures.

Polymorphic potential of SRF binding site of c-Fos gene promoter: in vitro study

Recently published in vivo observations have highlighted the presence of cruciform structures within the genome, suggesting their potential significance in the rapid recognition of the target sequence for transcription factor binding. In this in vitro study, we investigate the organization and stability of the sense (coding) strand within the Serum Response Element of the c-Fos gene promoter (c-Fos SRE), specifically focusing on segments spanning 12 to 36 nucleotides, centered around the CArG-box. Through a thorough examination of UV absorption patterns with varying temperatures, we identified the emergence of a remarkably stable structure, which we conclusively characterized as a hairpin using complementary 1H NMR experiments. Our research decisively ruled out the formation of homoduplexes, as confirmed by supplementary fluorescence experiments. Utilizing molecular dynamics simulations with atomic distance constraints derived from NMR data, we explored the structural intricacies of the compact hairpin. Notably, the loop consisting of the six-membered A/T sequence demonstrated substantial stabilization through extensive stacking, non-canonical inter-base hydrogen bonding, and hydrophobic clustering of thymine methyl groups. These findings suggest the potential of the c-Fos SRE to adopt a cruciform structure (consisting of two opposing hairpins), potentially providing a topological recognition site for the SRF transcription factor under cellular conditions. Our results should inspire further biochemical and in vivo studies to explore the functional implications of these non-canonical DNA structures.

A comprehensive study of Z-DNA density and its evolutionary implications in birds

BACKGROUND

Z-DNA, a left-handed helical form of DNA, plays a significant role in genomic stability and gene regulation. Its formation, associated with high GC content and repetitive sequences, is linked to genomic instability, potentially leading to large-scale deletions and contributing to phenotypic diversity and evolutionary adaptation.

RESULTS

In this study, we analyzed the density of Z-DNA-prone motifs of 154 avian genomes using the non-B DNA Motif Search Tool (nBMST). Our findings indicate a higher prevalence of Z-DNA motifs in promoter regions across all avian species compared to other genomic regions. A negative correlation was observed between Z-DNA density and developmental time in birds, suggesting that species with shorter developmental periods tend to have higher Z-DNA densities. This relationship implies that Z-DNA may influence the timing and regulation of development in avian species. Furthermore, Z-DNA density showed associations with traits such as body mass, egg mass, and genome size, highlighting the complex interactions between genome architecture and phenotypic characteristics. Gene Ontology (GO) analysis revealed that Z-DNA motifs are enriched in genes involved in nucleic acid binding, kinase activity, and translation regulation, suggesting a role in fine-tuning gene expression essential for cellular functions and responses to environmental changes. Additionally, the potential of Z-DNA to drive genomic instability and facilitate adaptive evolution underscores its importance in shaping phenotypic diversity.

CONCLUSIONS

This study emphasizes the role of Z-DNA as a dynamic genomic element contributing to gene regulation, genomic stability, and phenotypic diversity in avian species. Future research should experimentally validate these associations and explore the molecular mechanisms by which Z-DNA influences avian biology.

Mutations at BCL11B Exon 4 Associated with T Cell Acute Lymphoblastic Leukemia Are Facilitated by AID and Formation of Non-B DNA Conformations

One of the primary reasons behind the pathogenesis of T cell acute lymphoblastic leukemia (T-ALL) is the deregulation of the transcription factor BCL11B. The exon 4 of BCL11B harbors several driver mutations, which abolishes its DNA-binding ability. The high frequency of C > T or G > A conversion in close vicinity of AID (Activation-induced cytidine deaminase)-hotspot motifs in the deregulated gene prompted us to investigate the role of AID in BCL11B mutagenesis. Our results reveal that AID is expressed in T-ALL patient-derived cells, binds to BCL11B fragile region (FR) in exon 4 of T cells in vivo, and generates a signature mutation pattern in this region. The mutation frequency in BCL11B FR could be modulated upon overexpression of the AID gene in the knockout background, further suggesting the involvement of AID in BCL11B mutagenesis. Importantly, various lines of experimentation reveal that BCL11B FR could fold into parallel G-quadruplex, triplex, and hairpin structures, which could act as a replication/transcription block, causing mutagenesis. Thus, our results suggest that AID binds to BCL11B exon 4 due to non-B DNA formation, causing U:G mismatches or replication blocks, which, when repaired erroneously, generates deleterious mutations, resulting in loss of functionality of BCL11B, and thus becomes the cause of T-ALL.

RAG1/2 induces double-stranded DNA breaks at non-Ig loci in the proximity of single sequence repeats in developing B cells

In developing B cells, V(D)J gene recombination is initiated by the RAG1/2 endonuclease complex, introducing double-stranded DNA breaks (DSBs) in V, D, and J genes and resulting in the formation of the hypervariable parts of immunoglobulins (Ig). Persistent or aberrant RAG1/2 targeting is a potential threat to genome integrity. While RAG1 and RAG2 have been shown to bind various regions genome-wide, the in vivo off-target DNA damage instigated by RAG1/2 endonuclease remains less well understood. In the current study, we identified regions containing RAG1/2-induced DNA breaks in mouse pre-B cells on a genome-wide scale using a global DNA DSB detection strategy. We detected 1489 putative RAG1/2-dependent DSBs, most of which were located outside the Ig loci. DNA sequence motif analysis showed a specific enrichment of RAG1/2-induced DNA DSBs at GA- and CA-repeats and GC-rich motifs. These findings provide further insights into RAG1/2 off-target activity. The ability of RAG1/2 to introduce DSBs on the non-Ig loci during the endogenous V(D)J recombination emphasizes its genotoxic potential in developing lymphocytes.

Effects of interrupting residues on DNA dumbbell structures formed by CCTG tetranucleotide repeats associated with myotonic dystrophy type 2

Myotonic dystrophy type 2 (DM2) is a neurogenerative disease caused by caprylic/capric triglyceride (CCTG) tetranucleotide repeat expansions in intron 1 of the cellular nucleic acid-binding protein (CNBP) gene. Non-B DNA structures formed by CCTG repeats can promote genetic instability, whereas interrupting motifs of NCTG (N = A/T/G) within CCTG repeats help to maintain genomic stability. However, whether the interrupting motifs can affect DNA structures of CCTG repeats remains unclear. Here, we report that four CCTG repeats with an interrupting 3'-A/T/G residue formed dumbbell structures, whereas a non-interrupting 3'-C residue resulted in a multi-loop structure exhibiting conformational dynamics that may contribute to a higher tendency of escaping from DNA mismatch repair and causing repeat expansions. The results provide new structural insights into the genetic instability of CCTG repeats in DM2.

Enriched G4 forming repeats in the human genome are associated with robust well-coordinated transcription and reduced cancer transcriptome variation

Non-B DNA G-quadruplex (G4) structures with guanine (G) runs of 2 to 4 repeats can trigger opposing experimental transcriptional impacts. Here, we used bioinformatic algorithms to comprehensively assess correlations of steady-state RNA transcript levels with all putative G4 sequence (pG4) locations genome-wide in three mammalian genomes and in normal and tumor human tissues. The human pG4-containing gene set displays higher expression levels than the set without pG4, supporting and extending some prior observations. pG4 enrichment at transcription start sites (TSSs) in human, but not chimpanzee and mouse genomes, suggests possible positive selection pressure for pG4 at human TSS, potentially driving genome rewiring and gene expression divergence between human and chimpanzee. Comprehensive bioinformatic analyses revealed lower pG4-containing gene set variability in humans and among different pG4 genes in tumors. As G4 stabilizers are under therapeutic consideration for cancer and pathogens, such distinctions between human normal and tumor G4s along with other species merit attention. Furthermore, in germline and cancer sequences, the most mutagenic pG4 mapped to regions promoting alternative DNA structures. Overall findings establish high pG4 at TSS as a human genome attribute statistically associated with robust well-coordinated transcription and reduced cancer transcriptome variation with implications for biology, model organisms, and medicine.

R-loop functions in Brca1-associated mammary tumorigenesis

Deleterious accumulation of R-loops, a DNA-RNA hybrid structure, contributes to genome instability. They are associated with BRCA1 mutation-related breast cancer, an estrogen receptor α negative (ERα-) tumor type originating from luminal progenitor cells. However, a presumed causality of R-loops in tumorigenesis has not been established in vivo. Here, we overexpress mouse Rnaseh1 (Rh1-OE) in vivo to remove accumulated R-loops in Brca1-deficient mouse mammary epithelium (BKO). R-loop removal exacerbates DNA replication stress in proliferating BKO mammary epithelial cells, with little effect on homology-directed repair of double-strand breaks following ionizing radiation. Compared to their BKO counterparts, BKO-Rh1-OE mammary glands contain fewer luminal progenitor cells but more mature luminal cells. Despite a similar incidence of spontaneous mammary tumors in BKO and BKO-Rh1-OE mice, a significant percentage of BKO-Rh1-OE tumors express ERα and progesterone receptor. Our results suggest that rather than directly elevating the overall tumor incidence, R-loops influence the mammary tumor subtype by shaping the cell of origin for Brca1 tumors.

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.

Novel B-DNA dermatophyte assay for demonstration of canonical DNA in dermatophytes: Histopathologic characterization by artificial intelligence

We describe a novel assay and artificial intelligence-driven histopathologic approach identifying dermatophytes in human skin tissue sections (ie, B-DNA dermatophyte assay) and demonstrate, for the first time, the presence of dermatophytes in tissue using immunohistochemistry to detect canonical right-handed double-stranded (ds) B-DNA. Immunohistochemistry was performed using anti-ds-B-DNA monoclonal antibodies with formalin-fixed paraffin-embedded tissues to determine the presence of dermatophytes. The B-DNA assay resulted in a more accurate identification of dermatophytes, nuclear morphology, dimensions, and gene expression of dermatophytes (ie, optical density values) than periodic acid-Schiff (PAS), Grocott methenamine silver (GMS), or hematoxylin and eosin (H&E) stains. The novel assay guided by artificial intelligence allowed for efficient identification of different types of dermatophytes (eg, hyphae, microconidia, macroconidia, and arthroconidia). Using the B-DNA dermatophyte assay as a clinical tool for diagnosing dermatophytes is an alternative to PAS, GMS, and H&E as a fast and inexpensive way to accurately detect dermatophytosis and reduce the number of false negatives. Our assay resulted in superior identification, sensitivity, life cycle stages, and morphology compared to H&E, PAS, and GMS stains. This method detects a specific structural marker (ie, ds-B-DNA), which can assist with diagnosis of dermatophytes. It represents a significant advantage over methods currently in use.

Replication of [AT/TA]25 Microsatellite Sequences by Human DNA Polymerase δ Holoenzymes Is Dependent on dNTP and RPA Levels

Fragile sites are unstable genomic regions that are prone to breakage during stressed DNA replication. Several common fragile sites (CFS) contain A+T-rich regions including perfect [AT/TA] microsatellite repeats that may collapse into hairpins when in single-stranded DNA (ssDNA) form and coincide with chromosomal hotspots for breakage and rearrangements. While many factors contribute to CFS instability, evidence exists for replication stalling within [AT/TA] microsatellite repeats. Currently, it is unknown how stress causes replication stalling within [AT/TA] microsatellite repeats. To investigate this, we utilized FRET to characterize the structures of [AT/TA]25 sequences and also reconstituted lagging strand replication to characterize the progression of pol δ holoenzymes through A+T-rich sequences. The results indicate that [AT/TA]25 sequences adopt hairpins that are unwound by the major ssDNA-binding complex, RPA, and the progression of pol δ holoenzymes through A+T-rich sequences saturated with RPA is dependent on the template sequence and dNTP concentration. Importantly, the effects of RPA on the replication of [AT/TA]25 sequences are dependent on dNTP concentration, whereas the effects of RPA on the replication of A+T-rich, nonstructure-forming sequences are independent of dNTP concentration. Collectively, these results reveal complexities in lagging strand replication and provide novel insights into how [AT/TA] microsatellite repeats contribute to genome instability.

Genomic insights into the evolution and adaptation of secondary metabolite gene clusters in fungicolous species Cladobotryum mycophilum ATHUM6906

Mycophilic or fungicolous fungi can be found wherever fungi exist since they are able to colonize other fungi, which occupy a diverse range of habitats. Some fungicolous species cause important diseases on Basidiomycetes, and thus, they are the main reason for the destruction of mushroom cultivations. Nonetheless, despite their ecological significance, their genomic data remain limited. Cladobotryum mycophilum is one of the most aggressive species of the genus, destroying the economically important Agaricus bisporus cultivations. The 40.7 Mb whole genome of the Greek isolate ATHUM6906 is assembled in 16 fragments, including the mitochondrial genome and 2 small circular mitochondrial plasmids, in this study. This genome includes a comprehensive set of 12,282 protein coding, 56 rRNA, and 273 tRNA genes. Transposable elements, CAZymes, and pathogenicity related genes were also examined. The genome of C. mycophilum contained a diverse arsenal of genes involved in secondary metabolism, forming 106 biosynthetic gene clusters, which renders this genome as one of the most BGC abundant among fungicolous species. Comparative analyses were performed for genomes of species of the family Hypocreaceae. Some BGCs identified in C. mycophilum genome exhibited similarities to clusters found in the family Hypocreaceae, suggesting vertical heritage. In contrast, certain BGCs showed a scattered distribution among Hypocreaceae species or were solely found in Cladobotryum genomes. This work provides evidence of extensive BGC losses, horizontal gene transfer events, and formation of novel BGCs during evolution, potentially driven by neutral or even positive selection pressures. These events may increase Cladobotryum fitness under various environmental conditions and potentially during host-fungus interaction.

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.

DNA strand breaks at centromeres: Friend or foe?

Centromeres are large structural regions in the genomic DNA, which are essential for accurately transmitting a complete set of chromosomes to daughter cells during cell division. In humans, centromeres consist of highly repetitive α-satellite DNA sequences and unique epigenetic components, forming large proteinaceous structures required for chromosome segregation. Despite their biological importance, there is a growing body of evidence for centromere breakage across the cell cycle, including periods of quiescence. In this review, we provide an up-to-date examination of the distinct centromere environments at different stages of the cell cycle, highlighting their plausible contribution to centromere breakage. Additionally, we explore the implications of these breaks on centromere function, both in terms of negative consequences and potential positive effects.

R-Loop Functions in Brca1 -Associated Mammary Tumorigenesis

Excessive R-loops, a DNA-RNA hybrid structure, are associated with genome instability and BRCA1 mutation-related breast cancer. Yet the causality of R-loops in tumorigenesis remains unclear. Here we show that R-loop removal by Rnaseh1 overexpression (Rh1-OE) in Brca1 -knockout (BKO) mouse mammary epithelium exacerbates DNA replication stress without affecting homology-directed DNA repair. R-loop removal also diminishes luminal progenitors, the cell of origin for estrogen receptor α (ERα)-negative BKO tumors. However, R-loop reduction does not dampen spontaneous BKO tumor incidence. Rather, it gives rise to a significant percentage of ERα-expressing BKO tumors. Thus, R-loops reshape mammary tumor subtype rather than promoting tumorigenesis.

ZSCAN4-binding motif-TGCACAC is conserved and enriched in CA/TG microsatellites in both mouse and human genomes

The Zinc finger and SCAN domain containing 4 (ZSCAN4) protein, expressed transiently in pluripotent stem cells, gametes, and early embryos, extends telomeres, enhances genome stability, and improves karyotypes in mouse embryonic stem (mES) cells. To gain insights into the mechanism of ZSCAN4 function, we identified genome-wide binding sites of endogenous ZSCAN4 protein using ChIP-seq technology in mouse and human ES cells, where the expression of endogenous ZSCAN4 was induced by treating cells with retinoic acids or by overexpressing DUX4. We revealed that both mouse and human ZSCAN4 bind to the TGCACAC motif located in CA/TG microsatellite repeats, which are known to form unstable left-handed duplexes called Z-DNA that can induce double-strand DNA breaks and mutations. These ZSCAN4 binding sites are mostly located in intergenic and intronic regions of the genomes. By generating ZSCAN4 knockout in human ES cells, we showed that ZSCAN4 does not seem to be involved in transcriptional regulation. We also found that ectopic expression of mouse ZSCAN4 enhances the suppression of chromatin at ZSCAN4-binding sites. These results together suggest that some of the ZSCAN4 functions are mediated by binding to the error-prone regions in mouse and human genomes.

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.

The impact of chromatin on double-strand break repair: Imaging tools and discoveries

Eukaryotic nuclei are constantly being exposed to factors that break or chemically modify the DNA. Accurate repair of this DNA damage is crucial to prevent DNA mutations and maintain optimal cell function. To overcome the detrimental effects of DNA damage, a multitude of repair pathways has evolved. These pathways need to function properly within the different chromatin domains present in the nucleus. Each of these domains exhibit distinct molecular- and bio-physical characteristics that can influence the response to DNA damage. In particular, chromatin domains highly enriched for repetitive DNA sequences, such as nucleoli, centromeres and pericentromeric heterochromatin require tailored repair mechanisms to safeguard genome stability. Work from the past decades has led to the development of innovative imaging tools as well as inducible DNA damage techniques to gain new insights into the impact of these repetitive chromatin domains on the DNA repair process. Here we summarize these tools with a particular focus on Double-Strand Break (DSB) repair, and discuss the insights gained into our understanding of the influence of chromatin domains on DSB -dynamics and -repair pathway choice.

Replication of [AT/TA]25 microsatellite sequences by human DNA polymerase δ holoenzymes is dependent on dNTP and RPA levels

Difficult-to-Replicate Sequences (DiToRS) are natural impediments in the human genome that inhibit DNA replication under endogenous replication. Some of the most widely-studied DiToRS are A+T-rich, high "flexibility regions," including long stretches of perfect [AT/TA] microsatellite repeats that have the potential to collapse into hairpin structures when in single-stranded DNA (ssDNA) form and are sites of recurrent structural variation and double-stranded DNA (dsDNA) breaks. Currently, it is unclear how these flexibility regions impact DNA replication, greatly limiting our fundamental understanding of human genome stability. To investigate replication through flexibility regions, we utilized FRET to characterize the effects of the major ssDNA-binding complex, RPA, on the structure of perfect [AT/TA]25 microsatellite repeats and also re-constituted human lagging strand replication to quantitatively characterize initial encounters of pol δ holoenzymes with A+T-rich DNA template sequences. The results indicate that [AT/TA]25 sequences adopt hairpin structures that are unwound by RPA and pol δ holoenzymes support dNTP incorporation through the [AT/TA]25 sequences as well as an A+T-rich, non-structure forming sequence. Furthermore, the extent of dNTP incorporation is dependent on the sequence of the DNA template and the concentration of dNTPs. Importantly, the effects of RPA on the replication of [AT/TA]25 sequences are dependent on the concentration of dNTPs, whereas the effects of RPA on the replication of an A+T-rich, non-structure forming sequence are independent of dNTP concentration. Collectively, these results reveal complexities in lagging strand replication and provide novel insights into how flexibility regions contribute to genome instability.

Computational analysis on the dissemination of non-B DNA structural motifs in promoter regions of 1180 cellular genomes

The promoter regions of gene regulation are under evolutionary constraints and earlier studies uncovered that they are characterized by enrichment of functional non-B DNA structural signatures like curved DNA, cruciform DNA, G-quadruplex, triple-helical DNA, slipped DNA structures, and Z-DNA. However, these studies are restricted to a few model organisms, single non-B DNA motif types, or whole genomic sequences, and their comparative accumulation in promoter regions of different domains of life has not been reported comprehensively. In this study, for the first time, we investigated the preponderance of non-B DNA-prone motifs in promoter regions in 1180 genomes belonging to 28 taxonomic groups using the non-B DNA Motif Search Tool (nBMST). The trends suggest that they are predominant in promoters compared to the upstream and downstream regions of all three domains of life and variably linked to taxonomic groups. Cruciform DNA motif is the most abundant form of non-B DNA, spanning from archaea to lower eukaryotes. Curved DNA motifs are prominent in host-associated bacteria, and suppressed in mammals. Triplex-DNA and slipped DNA structure repeats are discretely dispersed in all lineages. G-quadruplex motifs are significantly enriched in mammals. We also observed that the unique enrichment of non-B DNA in promoters is strongly linked to genome GC, size, evolutionary time divergence, and ecological adaptations. Overall, our work systematically reports the unique non-B DNA structural landscape of cellular organisms from the perspective of the cis-regulatory code of genomes.

Non-B DNA structures as a booster of genome instability

Non-canonical secondary structures (NCSs) are alternative nucleic acid structures that differ from the canonical B-DNA conformation. NCSs often occur in repetitive DNA sequences and can adopt different conformations depending on the sequence. The majority of these structures form in the context of physiological processes, such as transcription-associated R-loops, G4s, as well as hairpins and slipped-strand DNA, whose formation can be dependent on DNA replication. It is therefore not surprising that NCSs play important roles in the regulation of key biological processes. In the last years, increasing published data have supported their biological role thanks to genome-wide studies and the development of bioinformatic prediction tools. Data have also highlighted the pathological role of these secondary structures. Indeed, the alteration or stabilization of NCSs can cause the impairment of transcription and DNA replication, modification in chromatin structure and DNA damage. These events lead to a wide range of recombination events, deletions, mutations and chromosomal aberrations, well-known hallmarks of genome instability which are strongly associated with human diseases. In this review, we summarize molecular processes through which NCSs trigger genome instability, with a focus on G-quadruplex, i-motif, R-loop, Z-DNA, hairpin, cruciform and multi-stranded structures known as triplexes.

Large-scale expansions of Friedreich's ataxia GAA•TTC repeats in an experimental human system: role of DNA replication and prevention by LNA-DNA oligonucleotides and PNA oligomers

Friedreich's ataxia (FRDA) is caused by expansions of GAA•TTC repeats in the first intron of the human FXN gene that occur during both intergenerational transmissions and in somatic cells. Here we describe an experimental system to analyze large-scale repeat expansions in cultured human cells. It employs a shuttle plasmid that can replicate from the SV40 origin in human cells or be stably maintained in S. cerevisiae utilizing ARS4-CEN6. It also contains a selectable cassette allowing us to detect repeat expansions that accumulated in human cells upon plasmid transformation into yeast. We indeed observed massive expansions of GAA•TTC repeats, making it the first genetically tractable experimental system to study large-scale repeat expansions in human cells. Further, GAA•TTC repeats stall replication fork progression, while the frequency of repeat expansions appears to depend on proteins implicated in replication fork stalling, reversal, and restart. Locked nucleic acid (LNA)-DNA mixmer oligonucleotides and peptide nucleic acid (PNA) oligomers, which interfere with triplex formation at GAA•TTC repeats in vitro, prevented the expansion of these repeats in human cells. We hypothesize, therefore, that triplex formation by GAA•TTC repeats stall replication fork progression, ultimately leading to repeat expansions during replication fork restart.

Deep statistical modelling of nanopore sequencing translocation times reveals latent non-B DNA structures

MOTIVATION

Non-canonical (or non-B) DNA are genomic regions whose three-dimensional conformation deviates from the canonical double helix. Non-B DNA play an important role in basic cellular processes and are associated with genomic instability, gene regulation, and oncogenesis. Experimental methods are low-throughput and can detect only a limited set of non-B DNA structures, while computational methods rely on non-B DNA base motifs, which are necessary but not sufficient indicators of non-B structures. Oxford Nanopore sequencing is an efficient and low-cost platform, but it is currently unknown whether nanopore reads can be used for identifying non-B structures.

RESULTS

We build the first computational pipeline to predict non-B DNA structures from nanopore sequencing. We formalize non-B detection as a novelty detection problem and develop the GoFAE-DND, an autoencoder that uses goodness-of-fit (GoF) tests as a regularizer. A discriminative loss encourages non-B DNA to be poorly reconstructed and optimizing Gaussian GoF tests allows for the computation of P-values that indicate non-B structures. Based on whole genome nanopore sequencing of NA12878, we show that there exist significant differences between the timing of DNA translocation for non-B DNA bases compared with B-DNA. We demonstrate the efficacy of our approach through comparisons with novelty detection methods using experimental data and data synthesized from a new translocation time simulator. Experimental validations suggest that reliable detection of non-B DNA from nanopore sequencing is achievable.

AVAILABILITY AND IMPLEMENTATION

Source code is available at https://github.com/bayesomicslab/ONT-nonb-GoFAE-DND.

Accurate sequencing of DNA motifs able to form alternative (non-B) structures

Approximately 13% of the human genome at certain motifs have the potential to form noncanonical (non-B) DNA structures (e.g., G-quadruplexes, cruciforms, and Z-DNA), which regulate many cellular processes but also affect the activity of polymerases and helicases. Because sequencing technologies use these enzymes, they might possess increased errors at non-B structures. To evaluate this, we analyzed error rates, read depth, and base quality of Illumina, Pacific Biosciences (PacBio) HiFi, and Oxford Nanopore Technologies (ONT) sequencing at non-B motifs. All technologies showed altered sequencing success for most non-B motif types, although this could be owing to several factors, including structure formation, biased GC content, and the presence of homopolymers. Single-nucleotide mismatch errors had low biases in HiFi and ONT for all non-B motif types but were increased for G-quadruplexes and Z-DNA in all three technologies. Deletion errors were increased for all non-B types but Z-DNA in Illumina and HiFi, as well as only for G-quadruplexes in ONT. Insertion errors for non-B motifs were highly, moderately, and slightly elevated in Illumina, HiFi, and ONT, respectively. Additionally, we developed a probabilistic approach to determine the number of false positives at non-B motifs depending on sample size and variant frequency, and applied it to publicly available data sets (1000 Genomes, Simons Genome Diversity Project, and gnomAD). We conclude that elevated sequencing errors at non-B DNA motifs should be considered in low-read-depth studies (single-cell, ancient DNA, and pooled-sample population sequencing) and in scoring rare variants. Combining technologies should maximize sequencing accuracy in future studies of non-B DNA.

Evolutionary honing in and mutational replacement: how long-term directed mutational responses to specific environmental pressures are possible

Recent results have shown that the human malaria-resistant hemoglobin S mutation originates de novo more frequently in the gene and in the population where it is of adaptive significance, namely, in the hemoglobin subunit beta gene compared to the nonresistant but otherwise identical 20A[Formula: see text]T mutation in the hemoglobin subunit delta gene, and in sub-Saharan Africans, who have been subject to intense malarial pressure for many generations, compared to northern Europeans, who have not. This finding raises a fundamental challenge to the traditional notion of accidental mutation. Here, we address this finding with the replacement hypothesis, according to which preexisting genetic interactions can lead directly and mechanistically to mutations that simplify and replace them. Thus, an evolutionary process under selection can gradually hone in on interactions of importance for the currently evolving adaptations, from which large-effect mutations follow that are relevant to these adaptations. We exemplify this hypothesis using multiple types of mutation, including gene fusion mutations, gene duplication mutations, A[Formula: see text]G mutations in RNA-edited sites and transcription-associated mutations, and place it in the broader context of a system-level view of mutation origination called interaction-based evolution. Potential consequences include that similarity of mutation pressures may contribute to parallel evolution in genetically related species, that the evolution of genome organization may be driven by mutational mechanisms, that transposable element movements may also be explained by replacement, and that long-term directed mutational responses to specific environmental pressures are possible. Such mutational phenomena need to be further tested by future studies in natural and artificial settings.

NBBC: a non-B DNA burden explorer in cancer

Alternate (non-B) DNA-forming structures, such as Z-DNA, G-quadruplex, triplex have demonstrated a potential role in cancer etiology. It has been found that non-B DNA-forming sequences can stimulate genetic instability in human cancer genomes, implicating them in the development of cancer and other genetic diseases. While there exist several non-B prediction tools and databases, they lack the ability to both analyze and visualize non-B data within a cancer context. Herein, we introduce NBBC, a non-B DNA burden explorer in cancer, that offers analyses and visualizations for non-B DNA forming motifs. To do so, we introduce 'non-B burden' as a metric to summarize the prevalence of non-B DNA motifs at the gene-, signature- and genomic site-levels. Using our non-B burden metric, we developed two analyses modules within a cancer context to assist in exploring both gene- and motif-level non-B type heterogeneity among gene signatures. NBBC is designed to serve as a new analysis and visualization platform for the exploration of non-B DNA, guided by non-B burden as a novel marker.

Stable bulged G-quadruplexes in the human genome: identification, experimental validation and functionalization

DNA sequence composition determines the topology and stability of G-quadruplexes (G4s). Bulged G-quadruplex structures (G4-Bs) are a subset of G4s characterized by 3D conformations with bulges. Current search algorithms fail to capture stable G4-B, making their genome-wide study infeasible. Here, we introduced a large family of computationally defined and experimentally verified potential G4-B forming sequences (pG4-BS). We found 478 263 pG4-BS regions that do not overlap 'canonical' G4-forming sequences in the human genome and are preferentially localized in transcription regulatory regions including R-loops and open chromatin. Over 90% of protein-coding genes contain pG4-BS in their promoter or gene body. We observed generally higher pG4-BS content in R-loops and their flanks, longer genes that are associated with brain tissue, immune and developmental processes. Also, the presence of pG4-BS on both template and non-template strands in promoters is associated with oncogenesis, cardiovascular disease and stemness. Our G4-BS models predicted G4-forming ability in vitro with 91.5% accuracy. Analysis of G4-seq and CUT&Tag data strongly supports the existence of G4-BS conformations genome-wide. We reconstructed a novel G4-B 3D structure located in the E2F8 promoter. This study defines a large family of G4-like sequences, offering new insights into the essential biological functions and potential future therapeutic uses of G4-B.

G-QINDER Tool: Bioinformatically Predicted Formation of Different Four-Stranded DNA Motifs from (GT)n and (GA)n Repeats

The recently introduced semi-orthogonal system of nucleic acid imaging offers a greatly improved method of identifying DNA sequences that are capable of adopting noncanonical structures. This paper uses our newly developed G-QINDER tool to identify specific repeat sequences that adopt unique structural motifs in DNA: TG and AG repeats. The structures were found to adopt a left-handed G-quadruplex form under extreme crowding conditions and a unique tetrahelical motif under certain other conditions. The tetrahelical structure likely consists of stacked AGAG-tetrads but, unlike G-quadruplexes, their stability does not appear to be dependent on the type of monovalent cation present. The occurrence of TG and AG repeats in genomes is not rare, and they are also found frequently in the regulatory regions of nucleic acids, so it is reasonable to assume that putative structural motifs, like other noncanonical forms, could play an important regulatory role in cells. This hypothesis is supported by the structural stability of the AGAG motif; its unfolding can occur even at physiological temperatures since the melting temperature is primarily dependent on the number of AG repeats in the sequence.

Breakpoints characterisation of the genomic deletions identified by MLPA in alkaptonuria patients

Until recently, mainly DNA sequencing has been used to identify variants within the gene coding for homogentisate dioxygenase (HGD, 3q13.33) that cause alkaptonuria (AKU), an autosomal recessive inborn error of metabolism of tyrosine. In order to identify possible larger genomic deletions we have developed a novel Multiplex Ligation-dependent Probe Amplification (MLPA) assay specific for this gene (HGD-MLPA) and tested it successfully in healthy controls and in patients carrying two known previously identified HGD deletions. Subsequently, we analysed 22 AKU patients in whom only one or none classical HGD variant was found by sequencing. Using HGD-MLPA and sequencing, we identified four larger deletions encompassing from 1 to 4 exons of this gene and we defined their exact breakpoints: deletion of exons 1-4 (c.1-8460_282 + 6727del), deletion of exons 5 and 6 (c.283-9199_434 + 1688del), deletion of exon 11 (c.775-1915_879 + 1293del), and deletion of exon 13 (c.1007-1709_1188 + 1121del). We suggest including MLPA in the DNA diagnostic protocols for AKU in cases where DNA sequencing does not lead to identification of both HGD variants.

DNA Polymerase Delta Exhibits Altered Catalytic Properties on Lysine Acetylation

DNA polymerase delta is the primary polymerase that is involved in undamaged nuclear lagging strand DNA replication. Our mass-spectroscopic analysis has revealed that the human DNA polymerase δ is acetylated on subunits p125, p68, and p12. Using substrates that simulate Okazaki fragment intermediates, we studied alterations in the catalytic properties of acetylated polymerase and compared it to the unmodified form. The current data show that the acetylated form of human pol δ displays a higher polymerization activity compared to the unmodified form of the enzyme. Additionally, acetylation enhances the ability of the polymerase to resolve complex structures such as G-quadruplexes and other secondary structures that might be present on the template strand. More importantly, the ability of pol δ to displace a downstream DNA fragment is enhanced upon acetylation. Our current results suggest that acetylation has a profound effect on the activity of pol δ and supports the hypothesis that acetylation may promote higher-fidelity DNA replication.

Aberrant HMGA2 Expression Sustains Genome Instability That Promotes Metastasis and Therapeutic Resistance in Colorectal Cancer

Colorectal cancer (CRC) is one of the most lethal cancers worldwide, accounting for nearly ~10% of all cancer diagnoses and deaths. Current therapeutic approaches have considerably increased survival for patients diagnosed at early stages; however, ~20% of CRC patients are diagnosed with late-stage, metastatic CRC, where 5-year survival rates drop to 6–13% and treatment options are limited. Genome instability is an enabling hallmark of cancer that confers increased acquisition of genetic alterations, mutations, copy number variations and chromosomal rearrangements. In that regard, research has shown a clear association between genome instability and CRC, as the accumulation of aberrations in cancer-related genes provides subpopulations of cells with several advantages, such as increased proliferation rates, metastatic potential and therapeutic resistance. Although numerous genes have been associated with CRC, few have been validated as predictive biomarkers of metastasis or therapeutic resistance. A growing body of evidence suggests a member of the High-Mobility Group A (HMGA) gene family, HMGA2, is a potential biomarker of metastatic spread and therapeutic resistance. HMGA2 is expressed in embryonic tissues and is frequently upregulated in aggressively growing cancers, including CRC. As an architectural, non-histone chromatin binding factor, it initiates chromatin decompaction to facilitate transcriptional regulation. HMGA2 maintains the capacity for stem cell renewal in embryonic and cancer tissues and is a known promoter of epithelial-to-mesenchymal transition in tumor cells. This review will focus on the known molecular mechanisms by which HMGA2 exerts genome protective functions that contribute to cancer cell survival and chemoresistance in CRC.

Suppressors of Break-Induced Replication in Human Cells

Short tandem DNA repeats are drivers of genome instability. To identify suppressors of break-induced mutagenesis human cells, unbiased genetic screens were conducted using a lentiviral shRNA library. The recipient cells possessed fragile non-B DNA that could induce DNA double-strand breaks (DSBs), integrated at an ectopic chromosomal site adjacent to a thymidine kinase marker gene. Mutagenesis of the thymidine kinase gene rendered cells resistant to the nucleoside analog ganciclovir (GCV). The screen identified genes that have established roles in DNA replication and repair, chromatin modification, responses to ionizing radiation, and genes encoding proteins enriched at replication forks. Novel loci implicated in BIR included olfactory receptors, the G0S2 oncogene/tumor suppressor axis, the EIF3H-METTL3 translational regulator, and the SUDS3 subunit of the Sin3A corepressor. Consistent with a role in suppressing BIR, siRNA knockdown of selected candidates increased the frequency of the GCV^r phenotype and increased DNA rearrangements near the ectopic non-B DNA. Inverse PCR and DNA sequence analyses showed that hits identified in the screen increased genome instability. Further analysis quantitated repeat-induced hypermutagenesis at the ectopic site and showed that knockdown of a primary hit, COPS2, induced mutagenic hotspots, remodeled the replication fork, and increased nonallelic chromosome template switches.

Sequence-specific binding behavior of coralyne toward triplex DNA: An ultrafast time-resolved fluorescence spectroscopy study

Triplex DNA structure has potential therapeutic application in inhibiting the expression of genes involved in cancer and other diseases. As a DNA-targeting antitumor and antibiotic drug, coralyne shows a remarkable binding propensity to triplex over canonical duplex and thus can modulate the stability of triplex structure, providing a prospective gene targeting strategy. Much less is known, however, about coralyne-binding interactions with triplex. By combining multiple steady-state spectroscopy with ultrafast fluorescence spectroscopy, we have investigated the binding behaviors of coralyne with typical triplexes. Upon binding with a G-containing triplex, the fluorescence of coralyne is markedly quenched owing to the photoinduced electron transfer (PET) of coralyne with the G base. Systematic studies show that the PET rates are sensitive to the binding configuration and local microenvironment, from which the coexisting binding modes of monomeric (full and partial) intercalation and aggregate stacking along the sugar-phosphate backbone are distinguished and their respective contributions are determined. It shows that coralyne has preferences for monomeric intercalation within CGG triplex and pure TAT triplex, whereas CGC⁺ triplex adopts mainly backbone binding of coralyne aggregates due to charge repulsion, revealing the sequence-specific binding selectivity. The triplex-DNA-induced aggregation of coralyne could be used as a probe for recognizing the water content in local DNA structures. The strong π-π stacking of intercalated coralyne monomer with base-triplets plays an important role in stabilizing the triplex structure. These results provide mechanistic insights for understanding the remarkable propensity of coralyne in selective binding to triplex DNA and shed light on the prospective applications of coralyne-triplex targeted anti-gene therapeutics.

Beyond the Primary Structure of Nucleic Acids: Potential Roles of Epigenetics and Noncanonical Structures in the Regulations of Plant Growth and Stress Responses

Epigenetics deals with changes in gene expression that are not caused by modifications in the primary sequence of nucleic acids. These changes beyond primary structures of nucleic acids not only include DNA/RNA methylation, but also other reversible conversions, together with histone modifications or RNA interference. In addition, under particular conditions (such as specific ion concentrations or protein-induced stabilization), the right-handed double-stranded DNA helix (B-DNA) can form noncanonical structures commonly described as "non-B DNA" structures. These structures comprise, for example, cruciforms, i-motifs, triplexes, and G-quadruplexes. Their formation often leads to significant differences in replication and transcription rates. Noncanonical RNA structures have also been documented to play important roles in translation regulation and the biology of noncoding RNAs. In human and animal studies, the frequency and dynamics of noncanonical DNA and RNA structures are intensively investigated, especially in the field of cancer research and neurodegenerative diseases. In contrast, noncanonical DNA and RNA structures in plants have been on the fringes of interest for a long time and only a few studies deal with their formation, regulation, and physiological importance for plant stress responses. Herein, we present a review focused on the main fields of epigenetics in plants and their possible roles in stress responses and signaling, with special attention dedicated to noncanonical DNA and RNA structures.

Detection of Z-DNA Structures in Supercoiled Genome

Z-DNAs are nucleic acid secondary structures that form within a special pattern of nucleotides and are promoted by DNA supercoiling. Through Z-DNA formation, DNA encodes information by dynamic changes in its secondary structure. A growing body of evidence indicates that Z-DNA formation can play a role in gene regulation; it can affect chromatin architecture and demonstrates its association with genomic instability, genetic diseases, and genome evolution. Many functional roles of Z-DNA are yet to be discovered highlighting the need for techniques to detect genome-wide folding of DNA into this structure. Here, we describe an approach to convert linear genome into supercoiled genome sponsoring Z-DNA formation. Applying permanganate-based methodology and high-throughput sequencing to supercoiled genome allows genome-wide detection of single-stranded DNA. Single-stranded DNA is characteristic of the junctions between the classical B-form of DNA and Z-DNA. Consequently, analysis of single-stranded DNA map provides snapshots of the Z-DNA conformation over the whole genome.

Methods to Study Z-DNA-Induced Genetic Instability

Alternative DNA structures that differ from the canonical B-DNA double helix, including Z-DNA, have received much attention recently due to their impact on DNA metabolic processes, including replication, transcription, and genome maintenance. Non-B-DNA-forming sequences can also stimulate genetic instability associated with disease development and evolution. Z-DNA can stimulate different types of genetic instability events in different species, and several different assays have been established to detect Z-DNA-induced DNA strand breaks and mutagenesis in prokaryotic and eukaryotic systems. In this chapter, we will introduce some of these methods including Z-DNA-induced mutation screening and detection of Z-DNA-induced strand breaks in mammalian cells, yeast, and mammalian cell extracts. Results from these assays should provide better insight into the mechanisms of Z-DNA-related genetic instability in different eukaryotic model systems.

Mitochondria and G-quadruplex evolution: an intertwined relationship

G-quadruplexes (G4s) are non-canonical structures formed in guanine (G)-rich sequences through stacked G tetrads by Hoogsteen hydrogen bonding. Several studies have demonstrated the existence of G4s in the genome of various organisms, including humans, and have proposed that G4s have a regulatory role in various cellular functions. However, little is known regarding the dissemination of G4s in mitochondria. In this review, we report the observation that the number of potential G4-forming sequences in the mitochondrial genome increases with the evolutionary complexity of different species, suggesting that G4s have a beneficial role in higher-order organisms. We also discuss the possible function of G4s in mitochondrial (mt)DNA and long noncoding (lnc)RNA and their role in various biological processes.

The remodeling of Z-DNA in the mammalian germ line

We recently discovered a novel biological process, the scheduled remodeling of Z-DNA structures in the developing fetal mouse male germ cells [Nat. Cell Biol. 24, 1141-1153]. This process affects purine/pyrimidine dinucleotide repeat (PPR) rich sequences, which can form stable left-handed Z-DNA structures. The protein that carries out this function is identified as ZBTB43, member of a large family of ZBTB proteins. Z-DNA remodeling by ZBTB43 not only coincides with global remodeling of DNA methylation and chromatin events in the male germ line, but it also is a prerequisite for de novo DNA methylation. When ZBTB43 changes DNA structure from the left-handed zigzag shaped Z-DNA to the regular smooth right-handed B-DNA, it also generates a suitable substrate for the de novo DNA methyltransferase, DNMT3A. By instructing de novo DNA methylation at PPRs in prospermatogonia, ZBTB43 safeguards epigenomic integrity of the male gamete. PPRs are fragile sequences, sites of large deletions and rearrangements in mammalian cells, and this fragility is thought to be due to Z-DNA structure formation rather than the sequence itself. This idea is now supported by the in vivo finding that DNA double strand breaks accumulate in mutant prospermatogonia which lack ZBTB43-dependent Z-DNA remodeling. If unrepaired, double stranded DNA breaks can lead to germ line mutations. Therefore, by preventing such breaks ZBTB43 is critical for guarding genome stability between generations. Here, we discuss the significance and implications of these findings in more detail.

Revisiting recent unusual drug-DNA complex structures: Implications for cancer and neurological disease diagnostics and therapeutics

DNA plays a crucial role in various biological processes such as protein production, replication, recombination etc. by adopting different conformations. Targeting these conformations by small molecules is not only important for disease therapy, but also improves our understanding of the mechanisms of disease development. In this review, we provide an overview of some of the most recent ligand-DNA complexes that have diagnostic and therapeutic applications in neurological diseases caused by abnormal repeat expansions and in cancer associated with mismatches. In addition, we have discussed important implications of ligands targeting higher-order structures, such as four-way junctions, G-quadruplexes and triplexes for drug discovery and DNA nanotechnology. We provide an overview of the results and perspectives of such structural studies on ligand-DNA interactions.

An updated overview of experimental and computational approaches to identify non-canonical DNA/RNA structures with emphasis on G-quadruplexes and R-loops

Multiple types of non-canonical nucleic acid structures play essential roles in DNA recombination and replication, transcription, and genomic instability and have been associated with several human diseases. Thus, an increasing number of experimental and bioinformatics methods have been developed to identify these structures. To date, most reviews have focused on the features of non-canonical DNA/RNA structure formation, experimental approaches to mapping these structures, and the association of these structures with diseases. In addition, two reviews of computational algorithms for the prediction of non-canonical nucleic acid structures have been published. One of these reviews focused only on computational approaches for G4 detection until 2020. The other mainly summarized the computational tools for predicting cruciform, H-DNA and Z-DNA, in which the algorithms discussed were published before 2012. Since then, several experimental and computational methods have been developed. However, a systematic review including the conformation, sequencing mapping methods and computational prediction strategies for these structures has not yet been published. The purpose of this review is to provide an updated overview of conformation, current sequencing technologies and computational identification methods for non-canonical nucleic acid structures, as well as their strengths and weaknesses. We expect that this review will aid in understanding how these structures are characterised and how they contribute to related biological processes and diseases.