Distinct DNA repair pathways cause genomic instability at alternative DNA structures

Formation of left-handed helices by C2'-fluorinated nucleic acids under physiological salt conditions

Recent findings in cell biology have rekindled interest in Z-DNA, the left-handed helical form of DNA. We report here that two minimally modified nucleosides, 2'F-araC and 2'F-riboG, induce the formation of the Z-form under low ionic strength. We show that oligomers entirely made of these two nucleosides exclusively produce left-handed duplexes that bind to the Zα domain of ADAR1. The effect of the two nucleotides is so dramatic that Z-form duplexes are the only species observed in 10 mM sodium phosphate buffer and neutral pH, and no B-form is observed at any temperature. Hence, in contrast to other studies reporting formation of Z/B-form equilibria by a preference for purine glycosidic angles in syn, our NMR and computational work revealed that sequential 2'F…H2N and intramolecular 3'H…N3' interactions stabilize the left-handed helix. The equilibrium between B- and Z- forms is slow in the 19F NMR time scale (≥ms), and each conformation exhibited unprecedented chemical shift differences in the 19F signals. This observation led to a reliable estimation of the relative population of B and Z species and enabled us to monitor B-Z transitions under different conditions. The unique features of 2'F-modified DNA should thus be a valuable addition to existing techniques for specific detection of new Z-binding proteins and ligands.

Elucidation of how the Mir-23-27-24 cluster regulates development and aging

MicroRNAs (miRNAs) are pivotal regulators of gene expression and are involved in biological processes spanning from early developmental stages to the intricate process of aging. Extensive research has underscored the fundamental role of miRNAs in orchestrating eukaryotic development, with disruptions in miRNA biogenesis resulting in early lethality. Moreover, perturbations in miRNA function have been implicated in the aging process, particularly in model organisms such as nematodes and flies. miRNAs tend to be clustered in vertebrate genomes, finely modulating an array of biological pathways through clustering within a single transcript. Although extensive research of their developmental roles has been conducted, the potential implications of miRNA clusters in regulating aging remain largely unclear. In this review, we use the Mir-23-27-24 cluster as a paradigm, shedding light on the nuanced physiological functions of miRNA clusters during embryonic development and exploring their potential involvement in the aging process. Moreover, we advocate further research into the intricate interplay among miRNA clusters, particularly the Mir-23-27-24 cluster, in shaping the regulatory landscape of aging.

TMPRSS2 is a tumor suppressor and its downregulation promotes antitumor immunity and immunotherapy response in lung adenocarcinoma

BACKGROUND

TMPRSS2, a key molecule for SARS-CoV-2 invading human host cells, has an association with cancer. However, its association with lung cancer remains insufficiently unexplored.

METHODS

In five bulk transcriptomics datasets, one single-cell RNA sequencing (scRNA-seq) dataset and one proteomics dataset for lung adenocarcinoma (LUAD), we explored associations between TMPRSS2 expression and immune signatures, tumor progression phenotypes, genomic features, and clinical prognosis in LUAD by the bioinformatics approach. Furthermore, we performed experimental validation of the bioinformatics findings.

RESULTS

TMPRSS2 expression levels correlated negatively with the enrichment levels of both immune-stimulatory and immune-inhibitory signatures, while they correlated positively with the ratios of immune-stimulatory/immune-inhibitory signatures. It indicated that TMPRSS2 levels had a stronger negative correlation with immune-inhibitory than with immune-stimulatory signatures. TMPRSS2 downregulation correlated with increased proliferation, stemness, genomic instability, tumor progression, and worse survival in LUAD. We further validated that TMPRSS2 was downregulated with tumor progression in the LUAD cohort we collected from Jiangsu Cancer Hospital, China. In vitro and in vivo experiments verified the association of TMPRSS2 deficiency with increased tumor cell proliferation and invasion and antitumor immunity in LUAD. Moreover, in vivo experiments demonstrated that TMPRSS2-knockdown tumors were more sensitive to BMS-1, an inhibitor of PD-1/PD-L1.

CONCLUSIONS

TMPRSS2 is a tumor suppressor, while its downregulation is a positive biomarker of immunotherapy in LUAD. Our data provide a potential link between lung cancer and pneumonia caused by SARS-CoV-2 infection.

Microsatellite break-induced replication generates highly mutagenized extrachromosomal circular DNAs

Extrachromosomal circular DNAs (eccDNAs) are produced from all regions of the eucaryotic genome. We used inverse PCR of non-B microsatellites capable of forming hairpin, triplex, quadruplex and AT-rich structures integrated at a common ectopic chromosomal site to show that these non-B DNAs generate highly mutagenized eccDNAs by replication-dependent mechanisms. Mutagenesis occurs within the non-B DNAs and extends several kilobases bidirectionally into flanking and nonallelic DNA. Each non-B DNA exhibits a different pattern of mutagenesis, while sister clones containing the same non-B DNA also display distinct patterns of recombination, microhomology-mediated template switching and base substitutions. Mutations include mismatches, short duplications, long nontemplated insertions, large deletions and template switches to sister chromatids and nonallelic chromosomes. Drug-induced replication stress or the depletion of DNA repair factors Rad51, the COPS2 signalosome subunit or POLη change the pattern of template switching and alter the eccDNA mutagenic profiles. We propose an asynchronous capture model based on break-induced replication from microsatellite-induced DNA double strand breaks to account for the generation and circularization of mutagenized eccDNAs and the appearance of genomic homologous recombination deficiency (HRD) scars. These results may help to explain the appearance of tumor eccDNAS and their roles in neoantigen production, oncogenesis and resistance to chemotherapy.

Microsatellite break-induced replication generates highly mutagenized extrachromosomal circular DNAs

Extrachromosomal circular DNAs (eccDNAs) are produced from all regions of the eucaryotic genome. In tumors, highly transcribed eccDNAs have been implicated in oncogenesis, neoantigen production and resistance to chemotherapy. Here we show that unstable microsatellites capable of forming hairpin, triplex, quadruplex and AT-rich structures generate eccDNAs when integrated at a common ectopic site in human cells. These non-B DNA prone microsatellites form eccDNAs by replication-dependent mechanisms. The microsatellite-based eccDNAs are highly mutagenized and display template switches to sister chromatids and to nonallelic chromosomal sites. High frequency mutagenesis occurs within the eccDNA microsatellites and extends bidirectionally for several kilobases into flanking DNA and nonallelic DNA. Mutations include mismatches, short duplications, longer nontemplated insertions and large deletions. Template switching leads to recurrent deletions and recombination domains within the eccDNAs. Template switching events are microhomology-mediated, but do not occur at all potential sites of complementarity. Each microsatellite exhibits a distinct pattern of recombination, microhomology choice and base substitution signature. Depletion of Rad51, the COPS2 signalosome subunit or POLη alter the eccDNA mutagenic profiles. We propose an asynchronous capture model based on break-induced replication from microsatellite-induced DNA breaks for the generation and circularization of mutagenized eccDNAs and genomic homologous recombination deficiency (HRD) scars.

The evolutionary loss of the Eh1 motif in FoxE1 in the lineage of placental mammals

Forkhead box E1 (FoxE1) protein is a transcriptional regulator known to play a major role in the development of the thyroid gland. By performing sequence alignments, we detected a deletion in FoxE1, which occurred in the evolution of mammals, near the point of divergence of placental mammals. This deletion led to the loss of the majority of the Eh1 motif, which was important for interactions with transcriptional corepressors. To investigate a potential mechanism for this deletion, we analyzed replication through the deletion area in mammalian cells with two-dimensional gel electrophoresis, and in vitro, using a primer extension reaction. We demonstrated that the area of the deletion presented an obstacle for replication in both assays. The exact position of polymerization arrest in primer extension indicated that it was most likely caused by a quadruplex DNA structure. The quadruplex structure hypothesis is also consistent with the exact borders of the deletion. The exact roles of these evolutionary changes in FoxE1 family proteins are still to be determined.

Elevated MSH2 MSH3 expression interferes with DNA metabolism in vivo

The Msh2-Msh3 mismatch repair (MMR) complex in Saccharomyces cerevisiae recognizes and directs repair of insertion/deletion loops (IDLs) up to ∼17 nucleotides. Msh2-Msh3 also recognizes and binds distinct looped and branched DNA structures with varying affinities, thereby contributing to genome stability outside post-replicative MMR through homologous recombination, double-strand break repair (DSBR) and the DNA damage response. In contrast, Msh2-Msh3 promotes genome instability through trinucleotide repeat (TNR) expansions, presumably by binding structures that form from single-stranded (ss) TNR sequences. We previously demonstrated that Msh2-Msh3 binding to 5' ssDNA flap structures interfered with Rad27 (Fen1 in humans)-mediated Okazaki fragment maturation (OFM) in vitro. Here we demonstrate that elevated Msh2-Msh3 levels interfere with DNA replication and base excision repair in vivo. Elevated Msh2-Msh3 also induced a cell cycle arrest that was dependent on RAD9 and ELG1 and led to PCNA modification. These phenotypes also required Msh2-Msh3 ATPase activity and downstream MMR proteins, indicating an active mechanism that is not simply a result of Msh2-Msh3 DNA-binding activity. This study provides new mechanistic details regarding how excess Msh2-Msh3 can disrupt DNA replication and repair and highlights the role of Msh2-Msh3 protein abundance in Msh2-Msh3-mediated genomic instability.

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.

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.

Dynamic alternative DNA structures in biology and disease

Repetitive elements in the human genome, once considered 'junk DNA', are now known to adopt more than a dozen alternative (that is, non-B) DNA structures, such as self-annealed hairpins, left-handed Z-DNA, three-stranded triplexes (H-DNA) or four-stranded guanine quadruplex structures (G4 DNA). These dynamic conformations can act as functional genomic elements involved in DNA replication and transcription, chromatin organization and genome stability. In addition, recent studies have revealed a role for these alternative structures in triggering error-generating DNA repair processes, thereby actively enabling genome plasticity. As a driving force for genetic variation, non-B DNA structures thus contribute to both disease aetiology and evolution.

Revisiting mutagenesis at non-B DNA motifs in the human genome

Non-B DNA structures formed by repetitive sequence motifs are known instigators of mutagenesis in experimental systems. Analyzing this phenomenon computationally in the human genome requires careful disentangling of intrinsic confounding factors, including overlapping and interrupted motifs and recurrent sequencing errors. Here, we show that accounting for these factors eliminates all signals of repeat-induced mutagenesis that extend beyond the motif boundary, and eliminates or dramatically shrinks the magnitude of mutagenesis within some motifs, contradicting previous reports. Mutagenesis not attributable to artifacts revealed several biological mechanisms. Polymerase slippage generates frequent indels within every variety of short tandem repeat motif, implicating slipped-strand structures. Interruption-correcting single nucleotide variants within short tandem repeats may originate from error-prone polymerases. Secondary-structure formation promotes single nucleotide variants within palindromic repeats and duplications within direct repeats. G-quadruplex motifs cause recurrent sequencing errors, whereas mutagenesis at Z-DNAs is conspicuously absent.

Z-DNA and Z-RNA: Methods-Past and Future

A quote attributed to Yogi Berra makes the observation that "It's tough to make predictions, especially about the future," highlighting the difficulties posed to an author writing a manuscript like the present. The history of Z-DNA shows that earlier postulates about its biology have failed the test of time, both those from proponents who were wildly enthusiastic in enunciating roles that till this day still remain elusive to experimental validation and those from skeptics within the larger community who considered the field a folly, presumably because of the limitations in the methods available at that time. If anything, the biological roles we now know for Z-DNA and Z-RNA were not anticipated by anyone, even when those early predictions are interpreted in the most favorable way possible. The breakthroughs in the field were made using a combination of methods, especially those based on human and mouse genetic approaches informed by the biochemical and biophysical characterization of the Zα family of proteins. The first success was with the p150 Zα isoform of ADAR1 (adenosine deaminase RNA specific), with insights into the functions of ZBP1 (Z-DNA-binding protein 1) following soon after from the cell death community. Just as the replacement of mechanical clocks by more accurate designs changed expectations about navigation, the discovery of the roles assigned by nature to alternative conformations like Z-DNA has forever altered our view of how the genome operates. These recent advances have been driven by better methodology and by better analytical approaches. This article will briefly describe the methods that were key to these discoveries and highlight areas where new method development is likely to further advance our knowledge.

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.

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.

The involvement of homeobox-C 4 in predicting prognosis and unraveling immune landscape across multiple cancers via integrated analysis

Background: There has been growing evidence that the aberrantly expressed Homeobox-C 4 (HOXC4) plays crucial roles in the development of some cancer types. However, it remains unclear as far as its expression patterns and prognostic significance are concerned, as is tumor immunity.

Methods: To investigate the expression levels and prognostic implications of HOXC4, multiple data sources were used in conjunction with quantitative real-time polymerase chain reaction (qRT-PCR) verification. Afterward, diverse immunological-related analyses, along with anti-cancer drug sensitivity, were performed in a number of cancer types. A further exploration of the underlying mechanisms of HOXC4 in tumorigenesis and immunity was carried out using the Gene Set Enrichment Analysis (GSEA) and the Gene Set Variation Analysis (GSVA).

Results: Based on extensive database mining, HOXC4 was ubiquitously expressed across 21 tumor cell lines and significantly higher than that of normal tissues in 21 tumor types. The outcome of survival analysis including overall survival (OS), disease-free interval (DFI), disease-specific survival (DSS) and progression-free interval (PFI) revealed that upregulation of HOXC4 expression in several cancers was associated with worse prognosis. Additionally, HOXC4 was observed to correlate closely with colon adenocarcinoma (COAD), head and neck squamous cell carcinoma (HNSC), lower grade glioma (LGG), liver hepatocellular carcinoma (LIHC), rectum adenocarcinoma (READ), and thyroid carcinoma (THCA) in terms of tumor immune cells infiltration. As a result of our comprehensive pan-cancer study, we have identified a significant link between the expression of HOXC4 and the efficacy of immunotherapy-related treatments, together with anti-cancer drug sensitivity. As a final note, HOXC4 was found to modulate multiple signaling pathways involved in tumorigenesis and immunity.

Conclusion: HOXC4 has been implicated in our study for the first time as an oncogene in cancers with a poor prognosis, potentially laying the groundwork for promising clinical biomarkers and immunotherapy approaches.

Creation and resolution of non-B-DNA structural impediments during replication

During replication, folding of the DNA template into non-B-form secondary structures provides one of the most abundant impediments to the smooth progression of the replisome. The core replisome collaborates with multiple accessory factors to ensure timely and accurate duplication of the genome and epigenome. Here, we discuss the forces that drive non-B structure formation and the evidence that secondary structures are a significant and frequent source of replication stress that must be actively countered. Taking advantage of recent advances in the molecular and structural biology of the yeast and human replisomes, we examine how structures form and how they may be sensed and resolved during replication.

DNA Bending Force Facilitates Z-DNA Formation under Physiological Salt Conditions

Z-DNA, a noncanonical helical structure of double-stranded DNA (dsDNA), plays pivotal roles in various biological processes, including transcription regulation. Mechanical stresses on dsDNA, such as twisting and stretching, help to form Z-DNA. However, the effect of DNA bending, one of the most common dsDNA deformations, on Z-DNA formation is utterly unknown. Here, we show that DNA bending induces the formation of Z-DNA, that is, more Z-DNA is formed as the bending force becomes stronger. We regulated the bending force on dsDNA by using D-shaped DNA nanostructures. The B-Z transition was observed by single-molecule fluorescence resonance energy transfer. We found that as the bending force became stronger, Z-DNA was formed at lower Mg²⁺ concentrations. When dsDNA contained cytosine methylations, the B-Z transition occurred at 78 mM Mg²⁺ (midpoint) in the absence of the bending force. However, the B-Z transition occurred at a 28-fold lower Mg²⁺ concentration (2.8 mM) in the presence of the bending force. Monte Carlo simulation suggested that the B-Z transition stabilizes the bent form via the formation of the B-Z junction with base extrusion, which effectively releases the bending stress on DNA. Our results clearly show that the bending force facilitates the B-Z transition under physiological salt conditions.

Z-DNA is remodelled by ZBTB43 in prospermatogonia to safeguard the germline genome and epigenome

Mutagenic purine–pyrimidine repeats can adopt the left-handed Z-DNA conformation. DNA breaks at potential Z-DNA sites can lead to somatic mutations in cancer or to germline mutations that are transmitted to the next generation. It is not known whether any mechanism exists in the germ line to control Z-DNA structure and DNA breaks at purine–pyrimidine repeats. Here we provide genetic, epigenomic and biochemical evidence for the existence of a biological process that erases Z-DNA specifically in germ cells of the mouse male foetus. We show that a previously uncharacterized zinc finger protein, ZBTB43, binds to and removes Z-DNA, preventing the formation of DNA double-strand breaks. By removing Z-DNA, ZBTB43 also promotes de novo DNA methylation at CG-containing purine–pyrimidine repeats in prospermatogonia. Therefore, the genomic and epigenomic integrity of the species is safeguarded by remodelling DNA structure in the mammalian germ line during a critical window of germline epigenome reprogramming.

Meng et al. show that ZBTB43 alters Z-DNA structures to prevent deleterious double-strand breaks and promote DNA methylation at purine–pyrimidine repeats in the mouse germ line.

APOBEC mutagenesis is low in most types of non-B DNA structures

While somatic mutations are known to be enriched in genome regions with non-canonical DNA secondary structure, the impact of particular mutagens still needs to be elucidated. Here, we demonstrate that in human cancers, the APOBEC mutagenesis is not enriched in direct repeats, mirror repeats, short tandem repeats, and G-quadruplexes, and even decreased below its level in B-DNA for cancer samples with very high APOBEC activity. In contrast, we observe that the APOBEC-induced mutational density is positively associated with APOBEC activity in inverted repeats (cruciform structures), where the impact of cytosine at the 3’-end of the hairpin loop is substantial. Surprisingly, the APOBEC-signature mutation density per TC motif in the single-stranded DNA of a G-quadruplex (G4) is lower than in the four-stranded part of G4 and in B-DNA. The APOBEC mutagenesis, as well as the UV-mutagenesis in melanoma samples, are absent in Z-DNA regions, owing to the depletion of their mutational signature motifs.

•

APOBEC mutagenesis is not enriched in most non-canonical DNA structures

•

Inverted repeats (cruciform structures) show increased APOBEC mutagenesis

•

G-quadruplex’s unstructured strand has low APOBEC-induced mutation density

•

Decrease of APOBEC mutagenesis in non-B DNA possibly associated with PrimPol

The nuclease activity of DNA2 promotes exonuclease 1-independent mismatch repair

The DNA mismatch repair (MMR) system is a major DNA repair system that corrects DNA replication errors. In eukaryotes, the MMR system functions via mechanisms both dependent on and independent of exonuclease 1 (EXO1), an enzyme that has multiple roles in DNA metabolism. Although the mechanism of EXO1-dependent MMR is well understood, less is known about EXO1-independent MMR. Here, we provide genetic and biochemical evidence that the DNA2 nuclease/helicase has a role in EXO1-independent MMR. Biochemical reactions reconstituted with purified human proteins demonstrated that the nuclease activity of DNA2 promotes an EXO1-independent MMR reaction via a mismatch excision-independent mechanism that involves DNA polymerase δ. We show that DNA polymerase ε is not able to replace DNA polymerase δ in the DNA2-promoted MMR reaction. Unlike its nuclease activity, the helicase activity of DNA2 is dispensable for the ability of the protein to enhance the MMR reaction. Further examination established that DNA2 acts in the EXO1-independent MMR reaction by increasing the strand-displacement activity of DNA polymerase δ. These data reveal a mechanism for EXO1-independent mismatch repair.

Gerontology through the Eyes of 21st Century Toxicology

Toxicants can cause cells to experience DNA damage, leading them to cellular senescence. Discovering mechanisms of cellular aging from birth to death will ease the process of understanding aging.

G-quadruplexes in helminth parasites

Parasitic helminths infecting humans are highly prevalent infecting ∼2 billion people worldwide, causing inflammatory responses, malnutrition and anemia that are the primary cause of morbidity. In addition, helminth infections of cattle have a significant economic impact on livestock production, milk yield and fertility. The etiological agents of helminth infections are mainly Nematodes (roundworms) and Platyhelminths (flatworms). G-quadruplexes (G4) are unusual nucleic acid structures formed by G-rich sequences that can be recognized by specific G4 ligands. Here we used the G4Hunter Web Tool to identify and compare potential G4 sequences (PQS) in the nuclear and mitochondrial genomes of various helminths to identify G4 ligand targets. PQS are nonrandomly distributed in these genomes and often located in the proximity of genes. Unexpectedly, a Nematode, Ascaris lumbricoides, was found to be highly enriched in stable PQS. This species can tolerate high-stability G4 structures, which are not counter selected at all, in stark contrast to most other species. We experimentally confirmed G4 formation for sequences found in four different parasitic helminths. Small molecules able to selectively recognize G4 were found to bind to Schistosoma mansoni G4 motifs. Two of these ligands demonstrated potent activity both against larval and adult stages of this parasite.

LncRNA MT1DP promotes cadmium-induced DNA replication stress by inhibiting chromatin recruitment of SMARCAL1

Cadmium (Cd) is a well-known carcinogenic metal and widespread environmental pollutant. The effect of Cd-induced carcinogenesis is partly due to accumulated DNA damage and chromosomal aberrations, but the exact mechanisms underlying the genotoxicity of Cd have not been clearly understood. Here, we found that one long non-coding RNA MT1DP is participated in Cd-induced DNA damage and replication stress. Through analyzing the residents from Cd-contaminated area in Southern China, we found that blood DNA repair genes are down-regulated in individuals with high urine Cd values compared to those with low urine Cd values, which contrast to the blood MT1DP levels. Through in vitro experiments, we found that MT1DP promotes Cd-induced DNA damage response, genome instability and replication fork stalling. Mechanically, upon Cd treatment, ATR is activated to enhance HIF-1α expression, which in turn promotes the transcription level of MT1DP. Subsequently MT1DP is recruited on the chromatin and binds to SMARCAL1 to competitive inhibit latter's interaction with RPA complexes, finally leading to increased replication stress and DNA damage. In summary, this study provides clear evidence for the role of epigenetic regulation on the genotoxic effect of Cd, and MT1DP-mediated replication stress may represent a novel mechanism for Cd-induced carcinogenesis.

Non-Canonical Helical Structure of Nucleic Acids Containing Base-Modified Nucleotides

Chemically modified nucleobases are thought to be important for therapeutic purposes as well as diagnosing genetic diseases and have been widely involved in research fields such as molecular biology and biochemical studies. Many artificially modified nucleobases, such as methyl, halogen, and aryl modifications of purines at the C8 position and pyrimidines at the C5 position, are widely studied for their biological functions. DNA containing these modified nucleobases can form non-canonical helical structures such as Z-DNA, G-quadruplex, i-motif, and triplex. This review summarizes the synthesis of chemically modified nucleotides: (i) methylation, bromination, and arylation of purine at the C8 position and (ii) methylation, bromination, and arylation of pyrimidine at the C5 position. Additionally, we introduce the non-canonical structures of nucleic acids containing these modifications.

Homologous recombination within repetitive DNA

Many microsatellite DNA sequences are able to form non-B form DNA secondary structures, such as hairpin loops, cruciforms, triplex DNA or G-quadruplexes. These DNA structures can form a significant impediment to DNA replication and repair, leading to DNA nicks, gaps, and breaks, which can be repaired by homologous recombination (HR). Recent work understanding HR at structure-forming repeats has focused on genetic requirements for replication fork restart, break induced replication (BIR) at broken forks, recombination during and after relocalization of breaks or stalled forks to the nuclear periphery, and how repair pathway choice and kinetics are navigated in the presence of a repeat tract. In this review, we summarize recent developments that illuminate the role of recombination in repairing DNA damage or causing tract length changes within repetitive DNA and its role in maintaining genome stability.

De Novo Mutation Rate Variation and Its Determinants in Chlamydomonas

De novo mutations are central for evolution, since they provide the raw material for natural selection by regenerating genetic variation. However, studying de novo mutations is challenging and is generally restricted to model species, so we have a limited understanding of the evolution of the mutation rate and spectrum between closely related species. Here, we present a mutation accumulation (MA) experiment to study de novo mutation in the unicellular green alga Chlamydomonas incerta and perform comparative analyses with its closest known relative, Chlamydomonas reinhardtii. Using whole-genome sequencing data, we estimate that the median single nucleotide mutation (SNM) rate in C. incerta is μ = 7.6 × 10-10, and is highly variable between MA lines, ranging from μ = 0.35 × 10-10 to μ = 131.7 × 10-10. The SNM rate is strongly positively correlated with the mutation rate for insertions and deletions between lines (r > 0.97). We infer that the genomic factors associated with variation in the mutation rate are similar to those in C. reinhardtii, allowing for cross-prediction between species. Among these genomic factors, sequence context and complexity are more important than GC content. With the exception of a remarkably high C→T bias, the SNM spectrum differs markedly between the two Chlamydomonas species. Our results suggest that similar genomic and biological characteristics may result in a similar mutation rate in the two species, whereas the SNM spectrum has more freedom to diverge.

The Simple Biology of Flipons and Condensates Enhances the Evolution of Complexity

The classical genetic code maps nucleotide triplets to amino acids. The associated sequence composition is complex, representing many elaborations during evolution of form and function. Other genomic elements code for the expression and processing of RNA transcripts. However, over 50% of the human genome consists of widely dispersed repetitive sequences. Among these are simple sequence repeats (SSRs), representing a class of flipons, that under physiological conditions, form alternative nucleic acid conformations such as Z-DNA, G4 quartets, I-motifs, and triplexes. Proteins that bind in a structure-specific manner enable the seeding of condensates with the potential to regulate a wide range of biological processes. SSRs also encode the low complexity peptide repeats to patch condensates together, increasing the number of combinations possible. In situations where SSRs are transcribed, SSR-specific, single-stranded binding proteins may further impact condensate formation. Jointly, flipons and patches speed evolution by enhancing the functionality of condensates. Here, the focus is on the selection of SSR flipons and peptide patches that solve for survival under a wide range of environmental contexts, generating complexity with simple parts.

Special Issue: A, B and Z: The Structure, Function and Genetics of Z-DNA and Z-RNA

It is now difficult to believe that a biological function for the left-handed Z-DNA and Z-RNA conformations was once controversial. The papers in this Special Issue, "Z-DNA and Z-RNA: from Physical Structure to Biological Function", are based on presentations at the ABZ2021 meeting that was held virtually on 19 May 2021 and provide evidence for several biological functions of these structures. The first of its kind, this international conference gathered over 200 scientists from many disciplines to specifically address progress in research involving Z-DNA and Z-RNA. These high-energy left-handed conformers of B-DNA and A-RNA are associated with biological functions and disease outcomes, as evidenced from both mouse and human genetic studies. These alternative structures, referred to as "flipons", form under physiological conditions, regulate type I interferon responses and induce necroptosis during viral infection. They can also stimulate genetic instability, resulting in adaptive evolution and diseases such as cancer. The meeting featured cutting-edge science that was, for the most part, unpublished. We plan for the ABZ meeting to reconvene in 2022.

Beyond the double helix: DNA structural diversity and the PDB

The determination of the double helical structure of DNA in 1953 remains the landmark event in the development of modern biological and biomedical science. This structure has also been the starting point for the determination of some 2000 DNA crystal structures in the subsequent 68 years. Their structural diversity has extended to the demonstration of sequence-dependent local structure in duplex DNA, to DNA bending in short and long sequences and in the DNA wound round the nucleosome, and to left-handed duplex DNAs. Beyond the double helix itself, in circumstances where DNA sequences are or can be induced to unwind from being duplex, a wide variety of topologies and forms can exist. Quadruplex structures, based on four-stranded cores of stacked G-quartets, are prevalent though not randomly distributed in the human and other genomes and can play roles in transcription, translation, and replication. Yet more complex folds can result in DNAs with extended tertiary structures and enzymatic/catalytic activity. The Protein Data Bank is the depository of all these structures, and the resource where structures can be critically examined and validated, as well as compared one with another to facilitate analysis of conformational and base morphology features. This review will briefly survey the major structural classes of DNAs and illustrate their significance, together with some examples of how the use of the Protein Data Bank by for example, data mining, has illuminated DNA structural concepts.

ERCC1 mutations impede DNA damage repair and cause liver and kidney dysfunction in patients

ERCC1-XPF is a multifunctional endonuclease involved in nucleotide excision repair (NER), interstrand cross-link (ICL) repair, and DNA double-strand break (DSB) repair. Only two patients with bi-allelic ERCC1 mutations have been reported, both of whom had features of Cockayne syndrome and died in infancy. Here, we describe two siblings with bi-allelic ERCC1 mutations in their teenage years. Genomic sequencing identified a deletion and a missense variant (R156W) within ERCC1 that disrupts a salt bridge below the XPA-binding pocket. Patient-derived fibroblasts and knock-in epithelial cells carrying the R156W substitution show dramatically reduced protein levels of ERCC1 and XPF. Moreover, mutant ERCC1 weakly interacts with NER and ICL repair proteins, resulting in diminished recruitment to DNA damage. Consequently, patient cells show strongly reduced NER activity and increased chromosome breakage induced by DNA cross-linkers, while DSB repair was relatively normal. We report a new case of ERCC1 deficiency that severely affects NER and considerably impacts ICL repair, which together result in a unique phenotype combining short stature, photosensitivity, and progressive liver and kidney dysfunction.

Non-B DNA: a major contributor to small- and large-scale variation in nucleotide substitution frequencies across the genome

Approximately 13% of the human genome can fold into non-canonical (non-B) DNA structures (e.g. G-quadruplexes, Z-DNA, etc.), which have been implicated in vital cellular processes. Non-B DNA also hinders replication, increasing errors and facilitating mutagenesis, yet its contribution to genome-wide variation in mutation rates remains unexplored. Here, we conducted a comprehensive analysis of nucleotide substitution frequencies at non-B DNA loci within noncoding, non-repetitive genome regions, their ±2 kb flanking regions, and 1-Megabase windows, using human-orangutan divergence and human single-nucleotide polymorphisms. Functional data analysis at single-base resolution demonstrated that substitution frequencies are usually elevated at non-B DNA, with patterns specific to each non-B DNA type. Mirror, direct and inverted repeats have higher substitution frequencies in spacers than in repeat arms, whereas G-quadruplexes, particularly stable ones, have higher substitution frequencies in loops than in stems. Several non-B DNA types also affect substitution frequencies in their flanking regions. Finally, non-B DNA explains more variation than any other predictor in multiple regression models for diversity or divergence at 1-Megabase scale. Thus, non-B DNA substantially contributes to variation in substitution frequencies at small and large scales. Our results highlight the role of non-B DNA in germline mutagenesis with implications to evolution and genetic diseases.

Structure-forming repeats and their impact on genome stability

Repetitive sequences throughout the genome are a major source of endogenous DNA damage, due to the propensity of many of them to form alternative non-B DNA structures that can interfere with replication, transcription, and DNA repair. These repetitive sequences are prone to breakage (fragility) and instability (changes in repeat number). Repeat fragility and expansions are linked to several diseases, including many cancers and neurodegenerative diseases, hence the importance of understanding the mechanisms that cause genome instability and contribute to these diseases. This review focuses on recent findings of mechanisms causing repeat fragility and instability, new associations between repeat expansions and genetic diseases, and potential therapeutic options to target repeat expansions.

BRCA1-Associated Protein Is a Potential Prognostic Biomarker and Is Correlated With Immune Infiltration in Liver Hepatocellular Carcinoma: A Pan-Cancer Analysis

Background

BRCA1-associated protein (BRAP) is a critical gene that regulates inflammation-related signaling pathway and affects patients’ prognosis in esophageal squamous cell carcinoma (ESCC). However, its roles in different cancers remain largely unknown.

Methods

BRAP expression in human pan-cancer was analyzed via the Genotype-Tissue Expression (GTEx) and The Cancer Genome Atlas (TCGA) database. Pearson correlation analysis was used to evaluate the association between BRAP expression with mismatch repair (MMR) gene mutation and DNA methyltransferase. We evaluated the influence of BRAP on clinical prognosis by univariate survival analysis. Moreover, the correlation between BRAP and tumor immune infiltration was analyzed via the Tumor Immune Evaluation Resource (TIMER) database. Pearson correlation analysis was used to investigate the correlation between BRAP expression and immune checkpoint genes expression.

Results

BRAP is abnormally overexpressed and significantly correlated with MMR gene mutation level and DNA methyltransferase expression in human pan-cancer. Univariate survival analysis showed that BRAP was significant with patients’ overall survival (OS) in six cancer types, disease-free interval (DFI) in three cancer types, and progression-free interval (PFI) in two cancer types. Remarkably, increased BRAP expression was strongly correlated with patients’ poor prognosis in liver hepatocellular carcinoma (LIHC), whether OS (P < 0.0001, hazard ratio (HR) = 1.1), DFI (P = 0.00099, HR = 1.06), or PFI (P = 0.00025, HR = 1.07). Moreover, a positive relationship was found between BRAP expression and immune infiltrating cells including B cell, CD4 + T cell, CD8 + T cell, dendritic cell, macrophage cell, and neutrophil cell in colon adenocarcinoma (COAD), kidney renal clear cell carcinoma (KIRC), and LIHC. Additionally, BRAP expression showed strong correlations with immune checkpoint genes in LIHC.

Conclusion

BRAP expression is increased in human pan-cancer samples compared with normal tissues. Overexpression of BRAP is correlated with poor prognosis and immune infiltration in multiple cancers, especially in LIHC. These findings suggest that BRAP may be used as a potential molecular biomarker for determining prognosis and immune infiltration in LIHC.

The Cyclically Seasonal Drosophila subobscura Inversion O7 Originated From Fragile Genomic Sites and Relocated Immunity and Metabolic Genes

Chromosome inversions are important contributors to standing genetic variation in Drosophila subobscura. Presently, the species is experiencing a rapid replacement of high-latitude by low-latitude inversions associated with global warming. Yet not all low-latitude inversions are correlated with the ongoing warming trend. This is particularly unexpected in the case of O7 because it shows a regular seasonal cycle that peaks in summer and rose with a heatwave. The inconsistent behavior of O7 across components of the ambient temperature suggests that is causally more complex than simply due to temperature alone. In order to understand the dynamics of O7, high-quality genomic data are needed to determine both the breakpoints and the genetic content. To fill this gap, here we generated a PacBio long read-based chromosome-scale genome assembly, from a highly homozygous line made isogenic for an O3 + 4 + 7 chromosome. Then we isolated the complete continuous sequence of O7 by conserved synteny analysis with the available reference genome. Main findings include the following: (i) the assembled O7 inversion stretches 9.936 Mb, containing > 1,000 annotated genes; (ii) O7 had a complex origin, involving multiple breaks associated with non-B DNA-forming motifs, formation of a microinversion, and ectopic repair in trans with the two homologous chromosomes; (iii) the O7 breakpoints carry a pre-inversion record of fragility, including a sequence insertion, and transposition with later inverted duplication of an Attacin immunity gene; and (iv) the O7 inversion relocated the major insulin signaling forkhead box subgroup O (foxo) gene in tight linkage with its antagonistic regulatory partner serine/threonine-protein kinase B (Akt1) and disrupted concerted evolution of the two inverted Attacin duplicates, reattaching them to dFOXO metabolic enhancers. Our findings suggest that O7 exerts antagonistic pleiotropic effects on reproduction and immunity, setting a framework to understand its relationship with climate change. Furthermore, they are relevant for fragility in genome rearrangement evolution and for current views on the contribution of breakage versus repair in shaping inversion-breakpoint junctions.

Somatic mutations in neurodegeneration: An update

Mosaicism, the presence of genomic differences between cells due to post-zygotic somatic mutations, is widespread in the human body, including within the brain. A role for this in neurodegenerative diseases has long been hypothesised, and technical developments are now allowing the question to be addressed in detail. The rapidly accumulating evidence is discussed in this review, with a focus on recent developments. Somatic mutations of numerous types may occur, including single nucleotide variants (SNVs), copy number variants (CNVs), and retrotransposon insertions. They could act as initiators or risk factors, especially if they arise in development, although they could also result from the disease process, potentially contributing to progression. In common sporadic neurodegenerative disorders, relevant mutations have been reported in synucleinopathies, comprising somatic gains of SNCA in Parkinson's disease and multiple system atrophy, and in Alzheimer's disease, where a novel recombination mechanism leading to somatic variants of APP, as well as an excess of somatic SNVs affecting tau phosphorylation, have been reported. In Mendelian repeat expansion disorders, mosaicism due to somatic instability, first detected 25 years ago, has come to the forefront. Brain somatic SNVs occur in DNA repair disorders, and there is evidence for a role of several ALS genes in DNA repair. While numerous challenges, and need for further validation, remain, this new, or perhaps rediscovered, area of research has the potential to transform our understanding of neurodegeneration.

Simple Repeats as Building Blocks for Genetic Computers

Processing of RNA involves heterogeneous nuclear ribonucleoproteins. The simple sequence repeats (SSRs) they bind can also adopt alternative DNA structures, like Z DNA, triplexes, G quadruplexes, and I motifs. Those SSRs capable of switching conformation under physiological conditions (called flipons) are genetic elements that can encode alternative RNA processing by their effects on RNA processivity, most likely as DNA:RNA hybrids. Flipons are elements of a binary, instructive genetic code directing how genomic sequences are compiled into transcripts. The combinatorial nature of this code provides a rich set of options for creating genetic computers able to reproduce themselves and use a heritable and evolvable code to optimize survival. The underlying computational logic potentiates a diverse set of genetic programs that modify cis-mediated heritability and disease risk.

Distinct mechanisms of mutagenic processing of alternative DNA structures by repair proteins

Repetitive sequences can form a variety of alternative DNA structures (non-B DNA) that can modulate transcription, replication, and repair. However, non-B DNA-forming sequences can also stimulate mutagenesis, and are enriched at mutation hotspots in human cancer genomes. Interestingly, different types of non-B DNA stimulate mutagenesis via distinct repair processing mechanisms.