Search

Microdeletions and gross deletions are important causes (~20%) of human inherited disease and their genomic locations are strongly influenced by the local DNA sequence environment. This notwithstanding, no study has systematically examined their underlying generative mechanisms. Here, we obtained 42,098 pathogenic microdeletions and gross deletions from the Human Gene Mutation Database (HGMD) that together form a continuum of germline deletions ranging in size from 1bp to 28,394,429bp. We analyzed the DNA sequence within 1-kb of the breakpoint junctions and found that the frequencies of non-B DNA-forming repeats, GC-content, and the presence of seven of 78 specific sequence motifs in the vicinity of pathogenic deletions correlated with deletion length for deletions of length ≤30 bp. Further, we found that the presence of DR, GQ and STR repeats is important for the formation of longer deletions (>30 bp) but not for the formation of shorter deletions (≤30 bp) whilst significantly (Chi-square test P-value30 bp). We provide evidence to support a functional distinction between microdeletions and gross deletions. Finally, we propose that a deletion length cut-off of 25-30bp may serve as an objective means to functionally distinguish microdeletions from gross deletions.

SummaryThe pre-mRNA life cycle requires intron processing; yet, how intron processing defects influence splicing and gene expression is unclear. Here, we find TTDN1, which is frequently mutated in non-photosensitive trichothiodystrophy (NP-TTD), functionally links intron lariat processing to the spliceosome. The conserved TTDN1 C-terminal region directly binds lariat debranching enzyme DBR1, while its N-terminal intrinsically disordered region (IDR) binds the intron binding complex (IBC). The IDR forms condensatesin vitroand is needed for IBC interaction. TTDN1 loss causes significant intron lariat accumulation, as well as splicing and gene expression defects, mirroring phenotypes observed in NP-TTD patient cells.Ttdn1Δ/Δmice recapitulate intron processing defects and neurodevelopmental phenotypes seen in NP-TTD. A DBR1-IDR fusion recruits DBR1 to the IBC and circumvents the requirement for TTDN1, indicating this tethering role as its major molecular function. Collectively, our findings unveil key functional connections between lariat processing, splicing outcomes, and NP-TTD molecular pathology.

Background and Motivation: Whilst DNA repeat expansions cause numerous heritable human disorders, their origins and underlying pathological mechanisms are often unclear.Method: We collated a dataset comprising 224 human repeat expansions encompassing 203 different genes, and performed a systematic analysis with respect to key features at the DNA-, RNA- and protein-levels. Comparison with controls without known pathogenicity and genomic regions lacking repeats, allowed the construction of the first model to discriminate repeat regions harboring pathogenic repeat expansions (DPREx).Results: At the DNA level, pathogenic repeat expansions exhibited stronger signals for DNA regulatory factors (e.g. H3K4me3, transcription factor-binding sites) in exons, promoters, 5’UTRs, and 5’genes but not significantly different in introns, 3’UTRs and 3’genes than controls. At the RNA-level, pathogenic repeat expansions showed lower free energy for forming RNA secondary structure and were closer to splice sites in introns, exons, promoters and 5’genes than controls. At the protein level, pathogenic repeat expansions preferred to form coil than other types of secondary structures, and tended to encode surface-located protein domains. Additionally, pathogenic repeat expansions were also enriched in non-B DNA structures. Guided by these features, DPREx (http://biomed.nscc-gz.cn/zhaolab/geneprediction/#/) achieved an Area Under the Curve (AUC) value of 0.88 in an independent dataset test. Conclusion: Pathogenic repeat expansions are located so as to exert a synergistic, multi-level influence on stress responses and inter-molecular connections involving DNA, RNA and proteins, thereby impacting the relationship between genotype and clinical phenotype.

Significance DNA repair is essential to life and to avoidance of genome instability and cancer. Xeroderma pigmentosum group G (XPG) protein acts in multiple DNA repair pathways, both as an active enzyme and as a scaffold for coordinating with other repair proteins. We present here the structure of the catalytic domain responsible for its DNA binding and nuclease activity. Our analysis provides structure-based hypotheses for how XPG recognizes its bubble DNA substrate and predictions of the structural impacts of XPG disease mutations associated with two phenotypically distinct diseases: xeroderma pigmentosum (XP, skin cancer prone) or Cockayne syndrome (XP/CS, severe progressive developmental defects)., Xeroderma pigmentosum group G (XPG) protein is both a functional partner in multiple DNA damage responses (DDR) and a pathway coordinator and structure-specific endonuclease in nucleotide excision repair (NER). Different mutations in the XPG gene ERCC5 lead to either of two distinct human diseases: Cancer-prone xeroderma pigmentosum (XP-G) or the fatal neurodevelopmental disorder Cockayne syndrome (XP-G/CS). To address the enigmatic structural mechanism for these differing disease phenotypes and for XPG’s role in multiple DDRs, here we determined the crystal structure of human XPG catalytic domain (XPGcat), revealing XPG-specific features for its activities and regulation. Furthermore, XPG DNA binding elements conserved with FEN1 superfamily members enable insights on DNA interactions. Notably, all but one of the known pathogenic point mutations map to XPGcat, and both XP-G and XP-G/CS mutations destabilize XPG and reduce its cellular protein levels. Mapping the distinct mutation classes provides structure-based predictions for disease phenotypes: Residues mutated in XP-G are positioned to reduce local stability and NER activity, whereas residues mutated in XP-G/CS have implied long-range structural defects that would likely disrupt stability of the whole protein, and thus interfere with its functional interactions. Combined data from crystallography, biochemistry, small angle X-ray scattering, and electron microscopy unveil an XPG homodimer that binds, unstacks, and sculpts duplex DNA at internal unpaired regions (bubbles) into strongly bent structures, and suggest how XPG complexes may bind both NER bubble junctions and replication forks. Collective results support XPG scaffolding and DNA sculpting functions in multiple DDR processes to maintain genome stability.

The massive amount of experimental DNA and RNA sequence information provides an encyclopedia for cell biology that requires computational tools for efficient interpretation. The ability to write and apply simple computing scripts propels the investigator beyond the boundaries of online analysis tools to more broadly interrogate laboratory experimental data and to integrate them with all available datasets to test and challenge hypotheses. Here we describe robust prototypic bash and C++ scripts with metrics and methods for validation that we have made publicly available to address the roles of non-B DNA-forming motifs in eliciting genetic instability and to query The Cancer Genome Atlas. Importantly, the methods presented provide practical data interpretation tools to examine fundamental relationships and to enable insights and correlations between alterations in gene expression patterns and patient outcome. The exemplary source codes described are simple and can be efficiently modified, elaborated, and applied to other relationships and areas of investigation.

Significance Common fragile sites (CFSs) are normal loci that are genetically unstable under normal and oncogenic replication stress. Pol eta has been proposed to play a key role in CFS replication. Here, we show that in the absence of Pol eta, replication at five specific CFS loci is perturbed, with fork pausing observed at several sites. Sequence analysis showed that certain pause sites are associated with the presence of non-B DNA motifs, while others are not. Importantly, pause sites are located within regions of increased genetic variation in healthy human populations that could be attributed to Pol eta activity. Our data unveil a role for Pol eta in overcoming replication stress, reducing DNA breakage, and promoting genetic variation at CFSs., Common fragile sites (CFSs) are difficult-to-replicate genomic regions that form gaps and breaks on metaphase chromosomes under replication stress. They are hotspots for chromosomal instability in cancer. Repetitive sequences located at CFS loci are inefficiently copied by replicative DNA polymerase (Pol) delta. However, translesion synthesis Pol eta has been shown to efficiently polymerize CFS-associated repetitive sequences in vitro and facilitate CFS stability by a mechanism that is not fully understood. Here, by locus-specific, single-molecule replication analysis, we identified a crucial role for Pol eta (encoded by the gene POLH) in the in vivo replication of CFSs, even without exogenous stress. We find that Pol eta deficiency induces replication pausing, increases initiation events, and alters the direction of replication-fork progression at CFS-FRA16D in both lymphoblasts and fibroblasts. Furthermore, certain replication pause sites at CFS-FRA16D were associated with the presence of non-B DNA-forming motifs, implying that non-B DNA structures could increase replication hindrance in the absence of Pol eta. Further, in Pol eta-deficient fibroblasts, there was an increase in fork pausing at fibroblast-specific CFSs. Importantly, while not all pause sites were associated with non-B DNA structures, they were embedded within regions of increased genetic variation in the healthy human population, with mutational spectra consistent with Pol eta activity. From these findings, we propose that Pol eta replicating through CFSs may result in genetic variations found in the human population at these sites.

International audience; Summary Central to genotoxic responses is their ability to sense highly specific signals to activate the appropriate repair response. We previously reported that the activation of the ASCC-ALKBH3 repair pathway is exquisitely specific to alkylation damage in human cells. Yet the mechanistic basis for the selectivity of this pathway was not immediately obvious. Here, we demonstrate that RNA but not DNA alkylation is the initiating signal for this process. Aberrantly methylated RNA is sufficient to recruit ASCC, while an RNA dealkylase suppresses ASCC recruitment during chemical alkylation. In turn, recruitment of ASCC during alkylation damage, which is mediated by the E3 ubiquitin ligase RNF113A, suppresses transcription and R-loop formation. We further show that alkylated pre-mRNA is sufficient to activate RNF113A E3 ligase in vitro in a manner dependent on its RNA binding Zn-finger domain. Together, our work identifies an unexpected role for RNA damage in eliciting a specific response to genotoxins.

Alveolar capillary dysplasia with misalignment of pulmonary veins (ACDMPV) is a rare lethal lung developmental disease. Affected infants manifest with severe respiratory distress and refractory pulmonary hypertension and uniformly die in the first month of life. Heterozygous point mutations or copy-number variant deletions involving FOXF1 and/or its upstream lung-specific enhancer on 16q24.1 have been identified in the vast majority of ACDMPV patients. We have previously described two unrelated families with a de novo pathogenic frameshift variant c.691_698del (p.Ala231Argfs*61) in the exon 1 of FOXF1. Here, we present a third unrelated ACDMPV family with the same de novo variant and propose that a direct tandem repeat of eight consecutive nucleotides GCGGCGGC within the ~ 4 kb CpG island in FOXF1 exon 1 is a novel mutation hotspot causative for ACDMPV.

Vigilin (Vgl1) is essential for heterochromatin formation, chromosome segregation, and mRNA stability and is associated with autism spectrum disorders and cancer: vigilin, for example, can suppress proto-oncogene c-fms expression in breast cancer. Conserved from yeast to humans, vigilin is an RNA-binding protein with 14 tandemly arranged nonidentical hnRNP K-type homology (KH) domains. Here, we report that vigilin depletion increased cell sensitivity to cisplatin- or ionizing radiation (IR)-induced cell death and genomic instability due to defective DNA repair. Vigilin depletion delayed dephosphorylation of IR-induced γ-H2AX and elevated levels of residual 53BP1 and RIF1 foci, while reducing Rad51 and BRCA1 focus formation, DNA end resection, and double-strand break (DSB) repair. We show that vigilin interacts with the DNA damage response (DDR) proteins RAD51 and BRCA1, and vigilin depletion impairs their recruitment to DSB sites. Transient hydroxyurea (HU)-induced replicative stress in vigilin-depleted cells increased replication fork stalling and blocked restart of DNA synthesis. Furthermore, histone acetylation promoted vigilin recruitment to DSBs preferentially in the transcriptionally active genome. These findings uncover a novel vigilin role in DNA damage repair with implications for autism and cancer-related disorders.

Arrival of the novel SARS-CoV-2 has launched a worldwide effort to identify both pre-approved and novel therapeutics targeting the viral proteome, highlighting the urgent need for efficient drug discovery strategies. Even with effective vaccines, infection is possible, and at-risk populations would benefit from effective drug compounds that reduce the lethality and lasting damage of COVID-19 infection. The CoV-2 MacroD-like macrodomain (Mac1) is implicated in viral pathogenicity by disrupting host innate immunity through its mono (ADP-ribosyl) hydrolase activity, making it a prime target for antiviral therapy. We therefore solved the structure of CoV-2 Mac1 from non-structural protein 3 (Nsp3) and applied structural and sequence-based genetic tracing, including newly determined A. pompejana MacroD2 and GDAP2 amino acid sequences, to compare and contrast CoV-2 Mac1 with the functionally related human DNA-damage signaling factor poly (ADP-ribose) glycohydrolase (PARG). Previously, identified targetable features of the PARG active site allowed us to develop a pharmacologically useful PARG inhibitor (PARGi). Here, we developed a focused chemical library and determined 6 novel PARGi X-ray crystal structures for comparative analysis. We applied this knowledge to discovery of CoV-2 Mac1 inhibitors by combining computation and structural analysis to identify PARGi fragments with potential to bind the distal-ribose and adenosyl pockets of the CoV-2 Mac1 active site. Scaffold development of these PARGi fragments has yielded two novel compounds, PARG-345 and PARG-329, that crystallize within the Mac1 active site, providing critical structure-activity data and a pathway for inhibitor optimization. The reported structural findings demonstrate ways to harness our PARGi synthesis and characterization pipeline to develop CoV-2 Mac1 inhibitors targeting the ADP-ribose active site. Together, these structural and computational analyses reveal a path for accelerating development of antiviral therapeutics from pre-existing drug optimization pipelines.

GRB2 controls radiation-induced DNA damage repair—a predictive biomarker for PARPi and radiation therapy outcome., DNA double-strand break (DSB) repair is initiated by MRE11 nuclease for both homology-directed repair (HDR) and alternative end joining (Alt-EJ). Here, we found that GRB2, crucial to timely proliferative RAS/MAPK pathway activation, unexpectedly forms a biophysically validated GRB2-MRE11 (GM) complex for efficient HDR initiation. GRB2-SH2 domain targets the GM complex to phosphorylated H2AX at DSBs. GRB2 K109 ubiquitination by E3 ubiquitin ligase RBBP6 releases MRE11 promoting HDR. RBBP6 depletion results in prolonged GM complex and HDR defects. GRB2 knockout increased MRE11-XRCC1 complex and Alt-EJ. Reconstitution with separation-of-function GRB2 mutant caused HDR deficiency and synthetic lethality with PARP inhibitor. Cell and cancer genome analyses suggest biomarkers of low GRB2 for noncanonical HDR deficiency and high MRE11 and GRB2 expression for worse survival in HDR-proficient patients. These findings establish GRB2’s role in binding, targeting, and releasing MRE11 to promote efficient HDR over Alt-EJ DSB repair, with implications for genome stability and cancer biology.

SummaryA critical question in genome stability is the nature of the chemical damage responsible for repair activation. We previously reported a novel pathway specifically activated during alkylation damage in human cells, where the E3 ubiquitin ligase RNF113A mediates the recruitment of the ASCC repair complex. Yet the mechanistic basis for the alkylation damage selectivity of this pathway remains unclear. Here, we demonstrate that RNA but not DNA alkylation is the initiating signal for this process. Aberrantly methylated RNA is sufficient to recruit ASCC, while an RNA dealkylase suppresses ASCC recruitment during chemical alkylation. This aberrant RNA methylation causes transcriptional repression in a manner dependent on the ASCC complex. We show that an alkylated pre-mRNA, or an RNA containing a single damaged base, is sufficient to activate RNF113A E3 activity in a phosphorylation-dependent manner. Together, our work identifies an unexpected role for RNA damage in eliciting a DNA repair response, and suggests that RNA may serve as the “canary in the coal mine” for sensing alkylation damage.

The DNA and protein complex known as chromatin is subject to posttranslational modifications (PTMs) that regulate cellular functions such that PTM dysregulation can lead to disease, including cancer. One critical PTM is acetylation/deacetylation, which is being investigated as a means to develop targeted cancer therapies. The histone acetyltransferase (HAT) family of proteins performs histone acetylation. In humans, MOF (hMOF), a member of the MYST family of HATs, acetylates histone H4 at lysine 16 (H4K16ac). MOF-mediated acetylation plays a critical role in the DNA damage response (DDR) and embryonic stem cell development. Functionally, MOF is found in two distinct complexes: NSL (nonspecific lethal) in humans and MSL (male-specific lethal) in flies. The NSL complex is also able to acetylate additional histone H4 sites. Dysregulation of MOF activity occurs in multiple cancers, including ovarian cancer, medulloblastoma, breast cancer, colorectal cancer, and lung cancer. Bioinformatics analysis of KAT8, the gene encoding hMOF, indicated that it is highly overexpressed in kidney tumors as part of a concerted gene coexpression program that can support high levels of chromosome segregation and cell proliferation. The linkage between MOF and tumor proliferation suggests that there are additional functions of MOF that remain to be discovered.

Homologous recombination/end joining (HR/HEJ)-deficient cancers with BRCA mutations utilize alternative DNA double-strand break repair pathways, particularly alternative non-homologous end joining or microhomology-mediated end joining (alt-EJ/MMEJ) during S and G2 cell cycle phases. Depletion of alt-EJ factors, including XRCC1, PARP1 and POLQ, is synthetically lethal with BRCA2 deficiency; yet, XRCC1 roles in HR-deficient cancers and replication stress are enigmatic. Here, we show that after replication stress, XRCC1 forms an active repair complex with POLQ and MRE11 that supports alt-EJ activity in vitro. BRCA2 limits XRCC1 recruitment and repair complex formation to suppress alt-EJ at stalled forks. Without BRCA2 fork protection, XRCC1 enables cells to complete DNA replication at the expense of increased genome instability by promoting MRE11-dependent fork resection and restart. High XRCC1 and MRE11 gene expression negatively impacts Kaplan–Meier survival curves and hazard ratios for HR-deficient breast cancer patients in The Cancer Genome Atlas. The additive effects of depleting both BRCA2 and XRCC1 indicate distinct pathways for replication restart. Our collective data show that XRCC1-mediated processing contributes to replication fork degradation, replication restart and chromosome aberrations in BRCA2-deficient cells, uncovering new roles of XRCC1 and microhomology-mediated repair mechanisms in HR-deficient cancers, with implications for chemotherapeutic strategies targeting POLQ and PARP activities.

DNA-dependent protein kinase catalytic subunit (DNA-PKcs) is a key member of the phosphatidylinositol-3 kinase-like (PIKK) family of protein kinases with critical roles in DNA-double strand break repair, transcription, metastasis, mitosis, RNA processing, and innate and adaptive immunity. The absence of DNA-PKcs from many model organisms has led to the assumption that DNA-PKcs is a vertebrate-specific PIKK. Here, we find that DNA-PKcs is widely distributed in invertebrates, fungi, plants, and protists, and that threonines 2609, 2638, and 2647 of the ABCDE cluster of phosphorylation sites are highly conserved amongst most Eukaryotes. Furthermore, we identify highly conserved amino acid sequence motifs and domains that are characteristic of DNA-PKcs relative to other PIKKs. These include residues in the Forehead domain and a novel motif we have termed YRPD, located in an α helix C-terminal to the ABCDE phosphorylation site loop. Combining sequence with biochemistry plus structural data on human DNA-PKcs unveils conserved sequence and conformational features with functional insights and implications. The defined generally progressive DNA-PKcs sequence diversification uncovers conserved functionality supported by Evolutionary Trace analysis, suggesting that for many organisms both functional sites and evolutionary pressures remain identical due to fundamental cell biology. The mining of cancer genomic data and germline mutations causing human inherited disease reveal that robust DNA-PKcs activity in tumors is detrimental to patient survival, whereas germline mutations compromising function are linked to severe immunodeficiency and neuronal degeneration. We anticipate that these collective results will enable ongoing DNA-PKcs functional analyses with biological and medical implications.

Summary: Sequences with the capacity to adopt alternative DNA structures have been implicated in cancer etiology; however, the mechanisms are unclear. For example, H-DNA-forming sequences within oncogenes have been shown to stimulate genetic instability in mammals. Here, we report that H-DNA-forming sequences are enriched at translocation breakpoints in human cancer genomes, further implicating them in cancer etiology. H-DNA-induced mutations were suppressed in human cells deficient in the nucleotide excision repair nucleases, ERCC1-XPF and XPG, but were stimulated in cells deficient in FEN1, a replication-related endonuclease. Further, we found that these nucleases cleaved H-DNA conformations, and the interactions of modeled H-DNA with ERCC1-XPF, XPG, and FEN1 proteins were explored at the sub-molecular level. The results suggest mechanisms of genetic instability triggered by H-DNA through distinct structure-specific, cleavage-based replication-independent and replication-dependent pathways, providing critical evidence for a role of the DNA structure itself in the etiology of cancer and other human diseases. : DNA sequences that can adopt alternative structures, such as H-DNA, have been implicated in cancer etiology. Zhao et al. found that such sequences are enriched at translocation breakpoints in human cancer genomes and that repair- and replication-related nucleases cleave H-DNA in both error-free “replication-related” and mutagenic “replication-independent” mechanisms. Keywords: DNA structure, H-DNA, chromosome translocation, mutation hotspot, DNA repair, nuclease, structural simulation

Background SNCA multiplication is a genomic cause of familial PD, showing dosage-dependent toxicity. Until now, nonallelic homologous recombination was suggested as the mechanism of SNCA duplication, based on various types of repetitive elements found in the spanning region of the breakpoints. However, the sequence at the breakpoint was analyzed only for 1 case. Objectives We have analyzed the breakpoint sequences of 6 patients with PD who had duplicated SNCA using whole-genome sequencing data to elucidate the mechanism of SNCA duplication. Methods Six patient samples with SNCA duplication underwent whole-genome sequencing. The duplicated regions were defined with nucleotide-resolution breakpoints, which were confirmed by junction polymerase chain reaction and Sanger sequencing. The search for potential non-B DNA-forming sequences and stem-loop structure predictions was conducted. Results Duplicated regions ranged from the smallest region of 718.3 kb to the largest one of 4,162 kb. Repetitive elements were found at 8 of the 12 breakpoint sequences on each side of the junction, but none of the pairs shared overt homologies. Five of these six junctions had microhomologies (2-4 bp) at the breakpoint, and a short stretch of sequences was inserted in 3 cases. All except one junction were located within or next to stem-loop structures. Conclusion Our study has determined that homologous recombination mechanisms involving repetitive elements are not the main cause of the duplication of SNCA. The presence of microhomology at the junctions and their position within stem-loop structures suggest that replication-based rearrangements may be a common mechanism for SNCA amplification. © 2020 International Parkinson and Movement Disorder Society.

Stalled DNA replication fork restart after stress as orchestrated by ATR kinase, BLM helicase, and structure-specific nucleases enables replication, cell survival, and genome stability. Here we unveil human exonuclease V (EXO5) as an ATR-regulated DNA structure-specific nuclease and BLM partner for replication fork restart. We find that elevated EXO5 in tumors correlates with increased mutation loads and poor patient survival, suggesting that EXO5 upregulation has oncogenic potential. Structural, mechanistic, and mutational analyses of EXO5 and EXO5-DNA complexes reveal a single-stranded DNA binding channel with an adjacent ATR phosphorylation motif (T88Q89) that regulates EXO5 nuclease activity and BLM binding identified by mass spectrometric analysis. EXO5 phospho-mimetic mutant rescues the restart defect from EXO5 depletion that decreases fork progression, DNA damage repair, and cell survival. EXO5 depletion furthermore rescues survival of FANCA-deficient cells and indicates EXO5 functions epistatically with SMARCAL1 and BLM. Thus, an EXO5 axis connects ATR and BLM in directing replication fork restart.

A hallmark of cancer is genomic instability, which can enable cancer cells to evade therapeutic strategies. Here we employed a computational approach to uncover mechanisms underlying cancer mutational burden by focusing upon relationships between 1) translocation breakpoints and the thousands of G4 DNA-forming sequences within retrotransposons impacting transcription and exemplifying probable non-B DNA structures and 2) transcriptome profiling and cancer mutations. We determined the location and number of G4 DNA-forming sequences in the Genome Reference Consortium Human Build 38 and found a total of 358,605 covering ∼13.4 million bases. By analyzing >97,000 unique translocation breakpoints from the Catalogue Of Somatic Mutations In Cancer (COSMIC), we found that breakpoints are overrepresented at G4 DNA-forming sequences within hominid-specific SVA retrotransposons, and generally occur in tumors with mutations in tumor suppressor genes, such as TP53. Furthermore, correlation analyses between mRNA levels and exome mutational loads from The Cancer Genome Atlas (TCGA) encompassing >450,000 gene-mutation regressions revealed strong positive and negative associations, which depended upon tissue of origin. The strongest positive correlations originated from genes not listed as cancer genes in COSMIC; yet, these show strong predictive power for survival in most tumor types by Kaplan-Meier estimation. Thus, correlation analyses of DNA structure and gene expression with mutation loads complement and extend more traditional approaches to elucidate processes shaping genomic instability in cancer. The combined results point to G4 DNA, activation of cell cycle/DNA repair pathways, and mitochondrial dysfunction as three major factors driving the accumulation of somatic mutations in cancer cells.

An emerging molecular understanding of RNA alkylation and its removal is transforming our knowledge of RNA biology and its interplay with cancer chemotherapy responses. DNA modifications are known to perform critical functions depending on the genome template, including gene expression, DNA replication timing, and DNA damage protection, yet current results suggest that the chemical diversity of DNA modifications pales in comparison to those on RNA. More than 150 RNA modifications have been identified to date, and their complete functional implications are still being unveiled. These include intrinsic roles such as proper processing and RNA maturation; emerging evidence has furthermore uncovered RNA modification "readers", seemingly analogous to those identified for histone modifications. These modification recognition factors may regulate mRNA stability, localization, and interaction with translation machinery, affecting gene expression. Not surprisingly, tumors differentially modulate factors involved in expressing these marks, contributing to both tumorigenesis and responses to alkylating chemotherapy. Here we describe the current understanding of RNA modifications and their removal, with a focus primarily on methylation and alkylation as functionally relevant changes to the transcriptome. Intriguingly, some of the same RNA modifications elicited by physiological processes are also produced by alkylating agents, thus blurring the lines between what is a physiological mark and a damage-induced modification. Furthermore, we find that a high level of gene expression of enzymes with RNA dealkylation activity is a sensitive readout for poor survival in four different cancer types, underscoring the likely importance of examining RNA dealkylation mechanisms to cancer biology and for cancer treatment and prognosis.

Missense/nonsense mutations and microdeletions/microinsertions (

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.

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, 5mCpG 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, 86.9 % of all SAMMs were assigned to the correct tissue of origin. Thus, although tumors from different tissues may have similar mutation patterns, their SAMMs often display signatures that are characteristic of specific tissues. Electronic supplementary material The online version of this article (doi:10.1007/s00439-015-1566-1) contains supplementary material, which is available to authorized users.

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.

SummaryAnalyses of chromosomal aberrations in human genetic disorders have revealed that inverted repeat sequences (IRs) often co-localize with endogenous chromosomal instability and breakage hotspots. Approximately 80% of all IRs in the human genome are short (

Cancer genome sequence data provide an invaluable resource for inferring the key mechanisms by which mutations arise in cancer cells, favoring their survival, proliferation and invasiveness. Here we examine recent advances in understanding the molecular mechanisms responsible for the predominant type of genetic alteration found in cancer cells, somatic single base substitutions (SBSs). Cytosine methylation, demethylation and deamination, charge transfer reactions in DNA, DNA replication timing, chromatin status and altered DNA proofreading activities are all now known to contribute to the mechanisms leading to base substitution mutagenesis. We review current hypotheses as to the major processes that give rise to SBSs and evaluate their relative relevance in the light of knowledge acquired from cancer genome sequencing projects and the study of base modifications, DNA repair and lesion bypass. Although gene expression data on APOBEC3B enzymes provide support for a role in cancer mutagenesis through U:G mismatch intermediates, the enzyme preference for single-stranded DNA may limit its activity genome-wide. For SBSs at both CG:CG and YC:GR sites, we outline evidence for a prominent role of damage by charge transfer reactions that follow interactions of the DNA with reactive oxygen species (ROS) and other endogenous or exogenous electron-abstracting molecules.

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.

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.

Although alternative DNA secondary structures (non-B DNA) can induce genomic rearrangements, their associated mutational spectra remain largely unknown. The helicase activity of WRN, which is absent in the human progeroid Werner syndrome, is thought to counteract this genomic instability. We determined non-B DNA-induced mutation frequencies and spectra in human U2OS osteosarcoma cells and assessed the role of WRN in isogenic knockdown (WRN-KD) cells using a supF gene mutation reporter system flanked by triplex- or Z-DNA-forming sequences. Although both non-B DNA and WRN-KD served to increase the mutation frequency, the increase afforded by WRN-KD was independent of DNA structure despite the fact that purified WRN helicase was found to resolve these structures in vitro. In U2OS cells, ∼70% of mutations comprised single-base substitutions, mostly at G·C base-pairs, with the remaining ∼30% being microdeletions. The number of mutations at G·C base-pairs in the context of NGNN/NNCN sequences correlated well with predicted free energies of base stacking and ionization potentials, suggesting a possible origin via oxidation reactions involving electron loss and subsequent electron transfer (hole migration) between neighboring bases. A set of ∼40,000 somatic mutations at G·C base pairs identified in a lung cancer genome exhibited similar correlations, implying that hole migration may also be involved. We conclude that alternative DNA conformations, WRN deficiency and lung tumorigenesis may all serve to increase the mutation rate by promoting, through diverse pathways, oxidation reactions that perturb the electron orbitals of neighboring bases. It follows that such “hole migration” is likely to play a much more widespread role in mutagenesis than previously anticipated.

Repetitive DNA motifs may fold into non-B DNA structures, including cruciforms/hairpins, triplexes, slipped conformations, quadruplexes, and left-handed Z-DNA, thereby representing chromosomal targets for DNA repair, recombination, and aberrant DNA synthesis leading to repeat expansion or genomic rearrangements associated with neurodegenerative and genomic disorders. Hairpins and quadruplexes also determined the relative abundances of simple sequence repeats (SSR) in vertebrate genomes, whereas strong base stacking has permitted the expansion of purine·pyrimidine-rich SSR during evolutionary time. SSR are enriched in regulatory and cancer-related gene classes, where they have been actively recruited to participate in both gene and protein functions. SSR polymorphic alleles in the population are associated with cancer susceptibility, including within genes that appear to share regulatory circuits involving reactive oxygen species. © 2009 Wiley-Liss, Inc.

Microsatellites are abundant in vertebrate genomes, but their sequence representation and length distributions vary greatly within each family of repeats (e.g., tetranucleotides). Biophysical studies of 82 synthetic single-stranded oligonucleotides comprising all tetra- and trinucleotide repeats revealed an inverse correlation between the stability of folded-back hairpin and quadruplex structures and the sequence representation for repeats ≥30 bp in length in nine vertebrate genomes. Alternatively, the predicted energies of base-stacking interactions correlated directly with the longest length distributions in vertebrate genomes. Genome-wide analyses indicated that unstable sequences, such as CAG:CTG and CCG:CGG, were over-represented in coding regions and that micro/minisatellites were recruited in genes involved in transcription and signaling pathways, particularly in the nervous system. Microsatellite instability (MSI) is a hallmark of cancer, and length polymorphism within genes can confer susceptibility to inherited disease. Sequences that manifest the highest MSI values also displayed the strongest base-stacking interactions; analyses of 62 tri- and tetranucleotide repeat-containing genes associated with human genetic disease revealed enrichments similar to those noted for micro/minisatellite-containing genes. We conclude that DNA structure and base-stacking determined the number and length distributions of microsatellite repeats in vertebrate genomes over evolutionary time and that micro/minisatellites have been recruited to participate in both gene and protein function.

Missense/nonsense mutations and microdeletions/microinsertions (21 bp) represent ∼ 76% of all mutations causing human inherited disease, and their occurrence has been associated with sequence motifs (direct, inverted, and mirror repeats; G-quartets) capable of adopting non-B DNA structures. We found that a significant proportion (∼ 21%) of both microdeletions and microinsertions occur within direct repeats, and are explicable by slipped misalignment. A novel mutational mechanism, DNA triplex formation followed by DNA repair, may explain ∼ 5% of microdeletions and microinsertions at mirror repeats. Further, G-quartets, direct, and inverted repeats also appear to play a prominent role in mediating missense mutations, whereas only direct and inverted repeats mediate nonsense mutations. We suggest a mutational mechanism involving slipped strand mispairing, slipped structure formation, and DNA repair, to explain ∼ 15% of missense and ∼ 12% of nonsense mutations yielding perfect direct repeats from imperfect repeats, or the extension of existing direct repeats. Similar proportions of missense and nonsense mutations were explicable by hairpin/loop formation and DNA repair, yielding perfect inverted repeats from imperfect repeats. We also propose a model for single base-pair substitution based on one-electron oxidation reactions at G-quadruplex DNA. Overall, the proposed mechanisms provide support for a role for non-B DNA structures in human gene mutagenesis.

Non-B DNA conformations adopted by certain types of DNA sequences promote genetic instabilities, especially gross rearrangements including translocations. We conclude the following: (a) slipped (hairpin) structures, cruciforms, triplexes, tetraplexes and i-motifs, and left-handed Z-DNA are formed in chromosomes and elicit profound genetic consequences via recombination-repair, (b) repeating sequences, probably in their non-B conformations, cause gross genomic rearrangements (translocations, deletions, insertions, inversions, and duplications), and (c) these rearrangements are the genetic basis for numerous human diseases including polycystic kidney disease, adrenoleukodystrophy, follicular lymphomas, and spermatogenic failure.

Homo(purine*pyrimidine) sequences (R*Y tracts) with mirror repeat symmetries form stable triplexes that block replication and transcription and promote genetic rearrangements. A systematic search was conducted to map the location of the longest R*Y tracts in the human genome in order to assess their potential function(s). The 814 R*Y tracts withor =250 uninterrupted base pairs were preferentially clustered in the pseudoautosomal region of the sex chromosomes and located in the introns of 228 annotated genes whose protein products were associated with functions at the cell membrane. These genes were highly expressed in the brain and particularly in genes associated with susceptibility to mental disorders, such as schizophrenia. The set of 1957 genes harboring the 2886 R*Y tracts withor =100 uninterrupted base pairs was additionally enriched in proteins associated with phosphorylation, signal transduction, development and morphogenesis. Comparisons of theor =250 bp R*Y tracts in the mouse and chimpanzee genomes indicated that these sequences have mutated faster than the surrounding regions and are longer in humans than in chimpanzees. These results support a role for long R*Y tracts in promoting recombination and genome diversity during evolution through destabilization of chromosomal DNA, thereby inducing repair and mutation.

Genomic rearrangements are a frequent source of instability, but the mechanisms involved are poorly understood. A 2.5-kbp poly(purine·pyrimidine) sequence from the human PKD1 gene, known to form non-B DNA structures, induced long deletions and other instabilities in plasmids that were mediated by mismatch repair and, in some cases, transcription. The breakpoints occurred at predicted non-B DNA structures. Distance measurements also indicated a significant proximity of alternating purine-pyrimidine and oligo(purine·pyrimidine) tracts to breakpoint junctions in 222 gross deletions and translocations, respectively, involved in human diseases. In 11 deletions analyzed, breakpoints were explicable by non-B DNA structure formation. We conclude that alternative DNA conformations trigger genomic rearrangements through recombination-repair activities.