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.

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.

O-GlcNAcylation of High Mobility Group Box 1 (HMGB1) Alters Its DNA Binding and DNA Damage Processing Activities.

Protein O-GlcNAcylation is an essential and dynamic regulator of myriad cellular processes, including DNA replication and repair. Proteomic studies have identified the multifunctional nuclear protein HMGB1 as O-GlcNAcylated, providing a potential link between this modification and DNA damage responses. Here, we verify the protein's endogenous modification at S100 and S107 and found that the major modification site is S100, a residue that can potentially influence HMGB1-DNA interactions. Using synthetic protein chemistry, we generated site-specifically O-GlcNAc-modified HMGB1 at S100 and characterized biochemically the effect of the sugar modification on its DNA binding activity. We found that O-GlcNAc alters HMGB1 binding to linear, nucleosomal, supercoiled, cruciform, and interstrand cross-linked damaged DNA, generally resulting in enhanced oligomerization on these DNA structures. Using cell-free extracts, we also found that O-GlcNAc reduces the ability of HMGB1 to facilitate DNA repair, resulting in error-prone processing of damaged DNA. Our results expand our understanding of the molecular consequences of O-GlcNAc and how it affects protein-DNA interfaces. Importantly, our work may also support a link between upregulated O-GlcNAc levels and increased rates of mutations in certain cancer states.

FANCM and FAAP24 maintain genome stability via cooperative as well as unique functions.

The DNA remodeling enzyme FANCM and its DNA-binding partner, FAAP24, constitute a complex involved in the activation of Fanconi anemia (FA) DNA damage response mechanism, but neither gene has distinct patient mutants. In this study, we created isogenic models for both FANCM and FAAP24 and investigated their integrated functions in DNA damage response. We found that FANCM and FAAP24 coordinately facilitate FA pathway activation and suppress sister chromatid exchange. Importantly, we show that FANCM and FAAP24 possess nonoverlapping functions such that FAAP24 promotes ATR-mediated checkpoint activation particularly in response to DNA crosslinking agents, whereas FANCM participates in recombination-independent interstrand crosslink repair by facilitating recruitment of lesion incision activities, which requires its translocase activity. Our data suggest that FANCM and FAAP24 play multiple, while not fully epistatic, roles in maintaining genomic integrity.

A tunable assay for modulators of genome-destabilizing DNA structures.

Regions of genomic instability are not random and often co-localize with DNA sequences that can adopt alternative DNA structures (i.e. non-B DNA, such as H-DNA). Non-B DNA-forming sequences are highly enriched at translocation breakpoints in human cancer genomes, representing an endogenous source of genetic instability. However, a further understanding of the mechanisms involved in non-B DNA-induced genetic instability is needed. Small molecules that can modulate the formation/stability of non-B DNA structures, and therefore the subsequent mutagenic outcome, represent valuable tools to study DNA structure-induced genetic instability. To this end, we have developed a tunable Förster resonance energy transfer (FRET)-based assay to detect triplex/H-DNA-destabilizing and -stabilizing ligands. The assay was designed by incorporating a fluorophore-quencher pair in a naturally-occurring H-DNA-forming sequence from a chromosomal breakpoint hotspot in the human c-MYC oncogene. By tuning triplex stability via buffer composition, the assay functions as a dual-reporter that can identify stabilizers and destabilizers, simultaneously. The assay principle was demonstrated using known triplex stabilizers, BePI and coralyne, and a complementary oligonucleotide to mimic a destabilizer, MCRa2. The potential of the assay was validated in a 384-well plate with 320 custom-assembled compounds. The discovery of novel triplex stabilizers/destabilizers may allow the regulation of genetic instability in human genomes.

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.

Circulating microRNA-590-5p functions as a liquid biopsy marker in non-small cell lung cancer.

Despite the availability of various diagnostic procedures, a tissue biopsy is still indispensable for the routine diagnosis of lung cancer. However, inaccurate diagnoses can occur, leading to inefficient cancer management. In this context, use of circulating microRNAs (miRNAs) may serve as diagnostic tools as liquid biopsies, and as biomarkers to better understand the molecular mechanisms involved in the progression of cancer. We identified miR-590-5p as a potential prognostic marker in the progression of non-small cell lung cancer (NSCLC). We were able to detect this miRNA in blood plasma samples of NSCLC patients through quantitative real-time PCR. Our data showed an ~7.5-fold downregulation of miR-590-5p in NSCLC patients compared to healthy controls, which correlated with several clinicopathological features. Further, overexpression of miR-590-5p led to decreased cell viability, proliferation, colony formation, migration, and invasion potential of lung cancer cells, whereas its knockdown showed the opposite effect. In addition, the levels of several proteins involved in the epithelial-to-mesenchymal transition negatively correlated with miR-590-5p levels in lung adenocarcinoma cells and tumors of NSCLC patients. Further, dual-luciferase reporter assays identified STAT3 as a direct target of miR-590-5p, which negatively regulated STAT3 activation and its downstream signaling molecules (eg, Cyclin D1, c-Myc, Vimentin, and ß-catenin) involved in tumorigenesis. Taken together, our study suggests that miR-590-5p functions as a tumor suppressor in NSCLC through regulating the STAT3 pathway, and may serve as a useful biomarker for the diagnosis/prognosis of NSCLC, and as a potential therapeutic target for the treatment of NSCLC.

Obesity and cancer: A mechanistic overview of metabolic changes in obesity that impact genetic instability.

Obesity, defined as a state of positive energy balance with a body mass index exceeding 30 kg/m2 in adults and 95th percentile in children, is an increasing global concern. Approximately one-third of the world's population is overweight or obese, and in the United States alone, obesity affects one in six children. Meta-analysis studies suggest that obesity increases the likelihood of developing several types of cancer, and with poorer outcomes, especially in children. The contribution of obesity to cancer risk requires a better understanding of the association between obesity-induced metabolic changes and its impact on genomic instability, which is a major driving force of tumorigenesis. In this review, we discuss how molecular changes during adipose tissue dysregulation can result in oxidative stress and subsequent DNA damage. This represents one of the many critical steps connecting obesity and cancer since oxidative DNA lesions can result in cancer-associated genetic instability. In addition, the by-products of the oxidative degradation of lipids (e.g., malondialdehyde, 4-hydroxynonenal, and acrolein), and gut microbiota-mediated secondary bile acid metabolites (e.g., deoxycholic acid and lithocholic acid), can function as genotoxic agents and tumor promoters. We also discuss how obesity can impact DNA repair efficiency, potentially contributing to cancer initiation and progression. Finally, we outline obesity-related epigenetic changes and identify the gaps in knowledge to be addressed for the development of better therapeutic strategies for the prevention and treatment of obesity-related cancers.

The homologous recombination protein RAD51D protects the genome from large deletions.

Homologous recombination (HR) is a DNA double-strand break (DSB) repair pathway that protects the genome from chromosomal instability. RAD51 mediator proteins (i.e. paralogs) are critical for efficient HR in mammalian cells. However, how HR-deficient cells process DSBs is not clear. Here, we utilized a loss-of-function HR-reporter substrate to simultaneously monitor HR-mediated gene conversion and non-conservative mutation events. The assay is designed around a heteroallelic duplication of the Aprt gene at its endogenous locus in isogenic Chinese hamster ovary cell lines. We found that RAD51D-deficient cells had a reduced capacity for HR-mediated gene conversion both spontaneously and in response to I-SceI-induced DSBs. Further, RAD51D-deficiency shifted DSB repair toward highly deleterious single-strand annealing (SSA) and end-joining processes that led to the loss of large chromosomal segments surrounding site-specific DSBs at an exceptionally high frequency. Deletions in the proximity of the break were due to a non-homologous end-joining pathway, while larger deletions were processed via a SSA pathway. Overall, our data revealed that, in addition to leading to chromosomal abnormalities, RAD51D-deficiency resulted in a high frequency of deletions advancing our understanding of how a RAD51 paralog is involved in maintaining genomic stability and how its deficiency may predispose cells to tumorigenesis.

Alternative DNA structure formation in the mutagenic human c-MYC promoter.

Mutation 'hotspot' regions in the genome are susceptible to genetic instability, implicating them in diseases. These hotspots are not random and often co-localize with DNA sequences potentially capable of adopting alternative DNA structures (non-B DNA, e.g. H-DNA and G4-DNA), which have been identified as endogenous sources of genomic instability. There are regions that contain overlapping sequences that may form more than one non-B DNA structure. The extent to which one structure impacts the formation/stability of another, within the sequence, is not fully understood. To address this issue, we investigated the folding preferences of oligonucleotides from a chromosomal breakpoint hotspot in the human c-MYC oncogene containing both potential G4-forming and H-DNA-forming elements. We characterized the structures formed in the presence of G4-DNA-stabilizing K+ ions or H-DNA-stabilizing Mg2+ ions using multiple techniques. We found that under conditions favorable for H-DNA formation, a stable intramolecular triplex DNA structure predominated; whereas, under K+-rich, G4-DNA-forming conditions, a plurality of unfolded and folded species were present. Thus, within a limited region containing sequences with the potential to adopt multiple structures, only one structure predominates under a given condition. The predominance of H-DNA implicates this structure in the instability associated with the human c-MYC oncogene.

Low dose radiation primed iNOS + M1macrophages modulate angiogenic programming of tumor derived endothelium.

Solid tumors are covered by stroma, which is hypoxic in nature and composed of various non-malignant components such as endothelial cells, fibroblasts, and pericytes that support tumor growth. Tumor stroma represents a mechanical barrier for tumor infiltration of CD8+ effector T cells in particular. In this context, our previous studies have demonstrated the therapeutic impact of Low-Dose Radiation (LDR)-primed and M1-retuned (iNOS+) peritumoral macrophages that produce inducible nitric oxide, have immunological roles on tumor infiltration of effector T cells, cancer-related inflammation, and subsequent tumor immune rejection in a mouse model of pancreatic cancer. These findings suggested a possible modification of tumor endothelium by LDR-primed macrophages. In line with these observations, here we demonstrate the influence of LDR in down-modulating HIF-1 in irradiated tumors in the course of polarization of irradiated tumor-associated macrophages toward an M1 phenotype. Furthermore, we demonstrate that M1 macrophages which are primed by LDR can directly influence angiogenic responses in eNOS+ endothelial cells which produce nitric oxide having both vascular and physiological roles. Furthermore, we demonstrate that naïve macrophages, upon differentiating to an M1 phenotype either by Th1 stimuli or LDR, potentially modify sphingosine-1-phosphate/VEGF-induced angiogenic signaling in tumor-derived endothelial cells with tumorigenic potential, thus indicating the significance of iNOS+ macrophages in modulating signaling in eNOS+ tumor-derived endothelium. Our study suggests that iNOS+ macrophages can activate tumor endothelium which may contribute to cancer-directed immunotherapy in particular.

HMGB1 interacts with XPA to facilitate the processing of DNA interstrand crosslinks in human cells.

Many effective agents used in cancer chemotherapy cause DNA interstrand crosslinks (ICLs), which covalently link both strands of the double helix together resulting in cytotoxicity. ICLs are thought to be processed by proteins from a variety of DNA repair pathways; however, a clear understanding of ICL recognition and repair processing in human cells is lacking. Previously, we found that the high mobility group box 1 (HMGB1) protein bound to triplex-directed psoralen ICLs (TFO-ICLs) in vitro, cooperatively with NER damage recognition proteins, promoted removal of UVC-induced lesions and facilitated error-free repair of TFO-ICLs in mouse fibroblasts. Here, we demonstrate that HMGB1 recognizes TFO-ICLs in human cells, and its depletion increases ICL-induced mutagenesis in human cells without altering the mutation spectra. In contrast, HMGB1 depletion in XPA-deficient human cells significantly altered the ICL-induced mutation spectrum from predominantly T→A to T→G transversions. Moreover, the recruitment of XPA and HMGB1 to the ICLs is co-dependent. Finally, we show that HMGB1 specifically introduces negative supercoils in ICL-containing plasmids in HeLa cell extracts. Taken together, our data suggest that in human cells, HMGB1 functions in association with XPA on ICLs and facilitates the formation of a favorable architectural environment for ICL repair processing.

Translocation and deletion breakpoints in cancer genomes are associated with potential non-B DNA-forming sequences.

Gross chromosomal rearrangements (including translocations, deletions, insertions and duplications) are a hallmark of cancer genomes and often create oncogenic fusion genes. An obligate step in the generation of such gross rearrangements is the formation of DNA double-strand breaks (DSBs). Since the genomic distribution of rearrangement breakpoints is non-random, intrinsic cellular factors may predispose certain genomic regions to breakage. Notably, certain DNA sequences with the potential to fold into secondary structures [potential non-B DNA structures (PONDS); e.g. triplexes, quadruplexes, hairpin/cruciforms, Z-DNA and single-stranded looped-out structures with implications in DNA replication and transcription] can stimulate the formation of DNA DSBs. Here, we tested the postulate that these DNA sequences might be found at, or in close proximity to, rearrangement breakpoints. By analyzing the distribution of PONDS-forming sequences within ±500 bases of 19 947 translocation and 46 365 sequence-characterized deletion breakpoints in cancer genomes, we find significant association between PONDS-forming repeats and cancer breakpoints. Specifically, (AT)n, (GAA)n and (GAAA)n constitute the most frequent repeats at translocation breakpoints, whereas A-tracts occur preferentially at deletion breakpoints. Translocation breakpoints near PONDS-forming repeats also recur in different individuals and patient tumor samples. Hence, PONDS-forming sequences represent an intrinsic risk factor for genomic rearrangements in cancer genomes.

New Perspectives on DNA and RNA Triplexes As Effectors of Biological Activity.

Since the first description of the canonical B-form DNA double helix, it has been suggested that alternative DNA, DNA-RNA, and RNA structures exist and act as functional genomic elements. Indeed, over the past few years it has become clear that, in addition to serving as a repository for genetic information, genomic DNA elicits biological responses by adopting conformations that differ from the canonical right-handed double helix, and by interacting with RNA molecules to form complex secondary structures. This review focuses on recent advances on three-stranded (triplex) nucleic acids, with an emphasis on DNA-RNA and RNA-RNA interactions. Emerging work reveals that triplex interactions between noncoding RNAs and duplex DNA serve as platforms for delivering site-specific epigenetic marks critical for the regulation of gene expression. Additionally, an increasing body of genetic and structural studies demonstrates that triplex RNA-RNA interactions are essential for performing catalytic and regulatory functions in cellular nucleoprotein complexes, including spliceosomes and telomerases, and for enabling protein recoding during programmed ribosomal frameshifting. Thus, evidence is mounting that DNA and RNA triplex interactions are implemented to perform a range of diverse biological activities in the cell, some of which will be discussed in this review.

Local DNA dynamics shape mutational patterns of mononucleotide repeats in human genomes.

Single base substitutions (SBSs) and insertions/deletions are critical for generating population diversity and can lead both to inherited disease and cancer. Whereas on a genome-wide scale SBSs are influenced by cellular factors, on a fine scale SBSs are influenced by the local DNA sequence-context, although the role of flanking sequence is often unclear. Herein, we used bioinformatics, molecular dynamics and hybrid quantum mechanics/molecular mechanics to analyze sequence context-dependent mutagenesis at mononucleotide repeats (A-tracts and G-tracts) in human population variation and in cancer genomes. SBSs and insertions/deletions occur predominantly at the first and last base-pairs of A-tracts, whereas they are concentrated at the second and third base-pairs in G-tracts. These positions correspond to the most flexible sites along A-tracts, and to sites where a 'hole', generated by the loss of an electron through oxidation, is most likely to be localized in G-tracts. For A-tracts, most SBSs occur in the direction of the base-pair flanking the tracts. We conclude that intrinsic features of local DNA structure, i.e. base-pair flexibility and charge transfer, render specific nucleotides along mononucleotide runs susceptible to base modification, which then yields mutations. Thus, local DNA dynamics contributes to phenotypic variation and disease in the human population.

Guanine holes are prominent targets for mutation in cancer and inherited disease.

Single base substitutions constitute the most frequent type of human gene mutation and are a leading cause of cancer and inherited disease. These alterations occur non-randomly in DNA, being strongly influenced by the local nucleotide sequence context. However, the molecular mechanisms underlying such sequence context-dependent mutagenesis are not fully understood. Using bioinformatics, computational and molecular modeling analyses, we have determined the frequencies of mutation at G • C bp in the context of all 64 5'-NGNN-3' motifs that contain the mutation at the second position. Twenty-four datasets were employed, comprising >530,000 somatic single base substitutions from 21 cancer genomes, >77,000 germline single-base substitutions causing or associated with human inherited disease and 16.7 million benign germline single-nucleotide variants. In several cancer types, the number of mutated motifs correlated both with the free energies of base stacking and the energies required for abstracting an electron from the target guanines (ionization potentials). Similar correlations were also evident for the pathological missense and nonsense germline mutations, but only when the target guanines were located on the non-transcribed DNA strand. Likewise, pathogenic splicing mutations predominantly affected positions in which a purine was located on the non-transcribed DNA strand. Novel candidate driver mutations and tissue-specific mutational patterns were also identified in the cancer datasets. We conclude that electron transfer reactions within the DNA molecule contribute to sequence context-dependent mutagenesis, involving both somatic driver and passenger mutations in cancer, as well as germline alterations causing or associated with inherited disease.

The somatic autosomal mutation matrix in cancer genomes.

DNA damage in somatic cells originates from both environmental and endogenous sources, giving rise to mutations through multiple mechanisms. When these mutations affect the function of critical genes, cancer may ensue. Although identifying genomic subsets of mutated genes may inform therapeutic options, a systematic survey of tumor mutational spectra is required to improve our understanding of the underlying mechanisms of mutagenesis involved in cancer etiology. Recent studies have presented genome-wide sets of somatic mutations as a 96-element vector, a procedure that only captures the immediate neighbors of the mutated nucleotide. Herein, we present a 32 × 12 mutation matrix that captures the nucleotide pattern two nucleotides upstream and downstream of the mutation. A somatic autosomal mutation matrix (SAMM) was constructed from tumor-specific mutations derived from each of 909 individual cancer genomes harboring a total of 10,681,843 single-base substitutions. In addition, mechanistic template mutation matrices (MTMMs) representing oxidative DNA damage, ultraviolet-induced DNA damage, (5m)CpG deamination, and APOBEC-mediated cytosine mutation, are presented. MTMMs were mapped to the individual tumor SAMMs to determine the maximum contribution of each mutational mechanism to the overall mutation pattern. A Manhattan distance across all SAMM elements between any two tumor genomes was used to determine their relative distance. Employing this metric, 89.5% of all tumor genomes were found to have a nearest neighbor from the same tissue of origin. When a distance-dependent 6-nearest neighbor classifier was used, 10.4% of the SAMMs had an Undetermined tissue of origin, and 92.2% of the remaining SAMMs were assigned to the correct tissue of origin. [corrected]. Thus, although tumors from different tissues may have similar mutation patterns, their SAMMs often display signatures that are characteristic of specific tissues.

Long non-coding RNA: A new paradigm for lung cancer.

Lung cancer is the leading cause of cancer-related deaths worldwide. Recent advances in whole genome transcriptome analysis have enabled the identification of numerous members of a novel class of non-coding RNAs, i.e., long non-coding RNAs (lncRNAs), which play important roles in a wide range of biological processes and whose deregulation causes human disease, including cancer. Herein we provide a comprehensive survey of lncRNAs associated with lung cancer, with particular focus on the functions that either facilitate or inhibit the progression of lung cancer and the pathways involved. Emerging data on the use of lncRNAs as biomarkers for the diagnosis and prognosis of cancer are also discussed. We cast this information within the wider perspective of lncRNA biogenesis and molecular functions in the cell. Relationships that exist between lncRNAs, genome-wide transcription, and lung cancer are discussed. Deepening our understanding on these processes is critical not only from a mechanistic standpoint, but also for the development of novel biomarkers and effective therapeutic targets for cancer patients.

DHX9 helicase is involved in preventing genomic instability induced by alternatively structured DNA in human cells.

Sequences that have the capacity to adopt alternative (i.e. non-B) DNA structures in the human genome have been implicated in stimulating genomic instability. Previously, we found that a naturally occurring intra-molecular triplex (H-DNA) caused genetic instability in mammals largely in the form of DNA double-strand breaks. Thus, it is of interest to determine the mechanism(s) involved in processing H-DNA. Recently, we demonstrated that human DHX9 helicase preferentially unwinds inter-molecular triplex DNA in vitro. Herein, we used a mutation-reporter system containing H-DNA to examine the relevance of DHX9 activity on naturally occurring H-DNA structures in human cells. We found that H-DNA significantly increased mutagenesis in small-interfering siRNA-treated, DHX9-depleted cells, affecting mostly deletions. Moreover, DHX9 associated with H-DNA in the context of supercoiled plasmids. To further investigate the role of DHX9 in the recognition/processing of H-DNA, we performed binding assays in vitro and chromatin immunoprecipitation assays in U2OS cells. DHX9 recognized H-DNA, as evidenced by its binding to the H-DNA structure and enrichment at the H-DNA region compared with a control region in human cells. These composite data implicate DHX9 in processing H-DNA structures in vivo and support its role in the overall maintenance of genomic stability at sites of alternatively structured DNA.