Deamination Of 5-methylcytosine Research Articles

Base-excision repair pathway shapes 5-methylcytosine deamination signatures in pan-cancer genomes.

Transition of cytosine to thymine in CpG dinucleotides is the most frequent type of mutation in cancer. This increased mutability is commonly attributed to the spontaneous deamination of 5-methylcytosine (5mC), which is normally repaired by the base-excision repair (BER) pathway. However, the contribution of 5mC deamination in the increasing diversity of cancer mutational signatures remains poorly explored. We integrate mutational signatures analysis in a large series of tumor whole genomes with lineage-specific epigenomic data to draw a detailed view of 5mC deamination in cancer. We uncover tumor type-specific patterns of 5mC deamination signatures in CpG and non-CpG contexts. We demonstrate that the BER glycosylase MBD4 preferentially binds to active chromatin and early replicating DNA, which correlates with lower mutational burden in these domains. We validate our findings by modeling BER deficiencies in isogenic cell models. Here, we establish MBD4 as the main actor responsible for 5mC deamination repair in humans.

EPID-22. MUTATIONAL SIGNATURE ASSOCIATED WITH HALOALKANE EXPOSURE ENRICHED IN FIREFIGHTERS IN THE ADULT GLIOMA STUDY

Abstract Using data from The Cancer Genome Atlas (TCGA) and Glioma Longitudinal Analysis (GLASS), we identified specific environmentally related mutational signatures in glioma including a signature associated with haloalkane exposure (SBS42) from the Catalogue of Somatic Mutational Signatures (COSMIC). Haloalkanes are widely used commercially including in flame retardants and fire extinguishants and are an intriguing finding given an observed increased glioma risk in firefighters. To validate our findings, we use new study participants with data on occupation (not available in either TCGA/GLASS) from the University of California, San Francisco Adult Glioma Study (AGS) to compare mutational signatures in tumors from 17 persons with glioma and documented occupational exposure to firefighting with those of 18 matched glioma patients without history of firefighting. Using whole exome sequencing, we identified environmental exposure signatures in tumor tissue from firefighters as well as in non-firefighters who had other relevant occupational exposure (e.g. auto/truck mechanic work). In these data, glioma is largely associated with aging and mutational signatures relating to endogenous mutational processes that correlate with age, such as spontaneous or enzymatic deamination of 5-methylcytosine. However, some gliomas have detectable signatures associated with exogenous mutational processes, such as SBS42 haloalkanes. Thus, we confirm detection of haloalkane signatures in persons likely highly exposed, i.e. long-term firefighters. Identifying exogenous mutational processes in cancers are important as they may inform public health intervention strategy to reduce mutagenesis and prevent cancer inception. Identifying occupational correlates with SBS42, associated with occupational exposure to haloalkanes, will pinpoint occupational hazards that may be avoidable. This is especially important for glioma and other cancers where exogenous mutagenesis is not well established.

Effect of C-to-T transition at CpG sites on tumor suppressor genes in tumor development in cattle evaluated by somatic mutation analysis in enzootic bovine leukosis.

Oncogenic transformation of normal cells is caused by mutations and chromosomal abnormalities in cancer-related genes. Enzootic bovine leukosis (EBL) is a malignant B-cell lymphoma caused by bovine leukemia virus (BLV) infection in cattle. Although a small fraction of BLV-infected cattle develops EBL after a long latent period, the mechanisms for oncogenesis in EBL cattle remain largely unknown. In this study, we analyzed the types and patterns of somatic mutations in cancer cells from 36 EBL cases, targeting 21 cancer-related genes. Various somatic mutations were identified in eight genes, TP53, KMT2D, CREBBP, KRAS, PTEN, NOTCH1, MYD88, and CARD11. In addition, TP53 gene was found to be mutated in 69.4% of EBL cases, with most being biallelic mutations. In some cases, associations were observed between the ages at which cattle had developed EBL and somatic mutation patterns; young onset of EBL possibly occurs due to high impact mutations affecting protein translation and biallelic mutations. Furthermore, nucleotide substitution patterns indicated that cytosine at CpG sites tended to be converted to thymine in many EBL cases, which was considered to be the result of spontaneous deamination of 5-methylcytosine. These results demonstrate how somatic mutations have occurred in cancer cells leading to EBL development, thereby explaining its pathogenic mechanism. These findings will contribute to a better understanding and future elucidation of disease progression in BLV infection.IMPORTANCEEnzootic bovine leukosis (EBL) is a malignant and lethal disease in cattle. Currently, there are no effective vaccines or therapeutic methods against bovine leukemia virus (BLV) infection, resulting in severe economic losses in livestock industry. This study provides a renewed hypothesis to explain the general mechanisms of EBL onset by combining the previous finding that several integration sites of BLV provirus can affect the increase in survival and proliferation of infected cells. We demonstrate that two additional random events are necessary for oncogenic transformation in infected cell clones, elucidating the reason why only few infected cattle develop EBL. Further exploration of somatic mutation and BLV integration sites could support this hypothesis more firmly, potentially contributing to the development of novel control methods for EBL onset.

Human DNA polymerase ε is a source of C>T mutations at CpG dinucleotides

C-to-T transitions in CpG dinucleotides are the most prevalent mutations in human cancers and genetic diseases. These mutations have been attributed to deamination of 5-methylcytosine (5mC), an epigenetic modification found on CpGs. We recently linked CpG>TpG mutations to replication and hypothesized that errors introduced by polymerase ε (Pol ε) may represent an alternative source of mutations. Here we present a new method called polymerase error rate sequencing (PER-seq) to measure the error spectrum of DNA polymerases in isolation. We find that the most common human cancer-associated Pol ε mutant (P286R) produces an excess of CpG>TpG errors, phenocopying the mutation spectrum of tumors carrying this mutation and deficiencies in mismatch repair. Notably, we also discover that wild-type Pol ε has a sevenfold higher error rate when replicating 5mCpG compared to C in other contexts. Together, our results from PER-seq and human cancers demonstrate that replication errors are a major contributor to CpG>TpG mutagenesis in replicating cells, fundamentally changing our understanding of this important disease-causing mutational mechanism.

Biochemical reconstitution of heat-induced mutational processes.

Non-enzymatic spontaneous deamination of 5-methylcytosine, producing thymine, is the proposed etiology of cancer mutational signature 1, which is the most predominant signature in all cancers. Here, the proposed mutational process was reconstituted using synthetic DNA and purified proteins. First, single-stranded DNA containing 5-methylcytosine at CpG context was incubated at an elevated temperature to accelerate spontaneous DNA damage. Then, the DNA was treated with uracil DNA glycosylase to remove uracil residues that were formed by deamination of cytosine. The resulting DNA was then used as a template for DNA synthesis by yeast DNA polymerase δ. The DNA products were analyzed by next-generation DNA sequencing, and mutation frequencies were quantified. The observed mutations after this process were exclusively C>T mutations at CpG context, which was very similar to signature 1. When 5-methylcytosine modification and uracil DNA glycosylase were both omitted, C>T mutations were produced on C residues in all sequence contexts, but these mutations were diminished by uracil DNA glycosylase-treatment. These results indicate that the CpG>TpG mutations were produced by the deamination of 5-methylcytosine. Additional mutations, mainly C>G, were introduced by yeast DNA polymerase ζ on the heat-damaged DNA, indicating that G residues of the templates were also damaged. However, the damage on G residues was not converted to mutations with DNA polymerase δ or ε.

Somatic CpG hypermutation is associated with mismatch repair deficiency in cancer

Somatic hypermutation in cancer has gained momentum with the increased use of tumour mutation burden as a biomarker for immune checkpoint inhibitors. Spontaneous deamination of 5-methylcytosine to thymine at CpG dinucleotides is one of the most ubiquitous endogenous mutational processes in normal and cancer cells. Here, we performed a systematic investigation of somatic CpG hypermutation at a pan-cancer level. We studied 30,191 cancer patients and 103 cancer types and developed an algorithm to identify somatic CpG hypermutation. Across cancer types, we observed the highest prevalence in paediatric leukaemia (3.5%), paediatric high-grade glioma (1.7%), and colorectal cancer (1%). We discovered germline variants and somatic mutations in the mismatch repair complex MutSα (MSH2-MSH6) as genetic drivers of somatic CpG hypermutation in cancer, which frequently converged on CpG sites and TP53 driver mutations. We further observe an association between somatic CpG hypermutation and response to immune checkpoint inhibitors. Overall, our study identified novel cancer types that display somatic CpG hypermutation, strong association with MutSα-deficiency, and potential utility in cancer immunotherapy.

Simultaneous assessment of human genome and methylome data in a single experiment using limited deamination of methylated cytosine.

Multiomics require concerted recording of independent information, ideally from a single experiment. In this study, we introduce RIMS-seq2, a high-throughput technique to simultaneously sequence genomes and overlay methylation information while requiring only a small modification of the experimental protocol for high-throughput DNA sequencing to include a controlled deamination step. Importantly, the rate of deamination of 5-methylcytosine is negligible and thus does not interfere with standard DNA sequencing and data processing. Thus, RIMS-seq2 libraries from whole- or targeted-genome sequencing show the same germline variation calling accuracy and sensitivity compared with standard DNA-seq. Additionally, regional methylation levels provide an accurate map of the human methylome.

Whole‑exome sequencing reveals Lewis lung carcinoma is a hypermutated Kras/Nras–mutant cancer with extensive regional mutation clusters in its genome

Lewis lung carcinoma (LLC), as a widely used preclinical cancer model, has still not been genetically and genomically characterized. Here, we performed a whole–exome sequencing analysis on the LLC cell line to elucidate its molecular characteristics and etiologies. Our data showed that LLC originated from a male mouse belonging to C57BL/6L (a transitional strain between C57BL/6J and C57BL/6N) and contains substantial somatic SNV and InDel mutations (> 20,000). Extensive regional mutation clusters are present in its genome, which were caused mainly by the mutational processes underlying the SBS1, SBS5, SBS15, SBS17a, and SBS21 signatures during frequent structural rearrangements. Thirty three deleterious mutations are present in 30 cancer genes including Kras, Nras, Trp53, Dcc, and Cacna1d. Cdkn2a and Cdkn2b are biallelically deleted from the genome. Five pathways (RTK/RAS, p53, cell cycle, TGFB, and Hippo) are oncogenically deregulated or affected. The major mutational processes in LLC include chromosomal instability, exposure to metabolic mutagens, spontaneous 5–methylcytosine deamination, defective DNA mismatch repair, and reactive oxygen species. Our data also suggest that LLC is a lung cancer similar to human lung adenocarcinoma. This study lays a molecular basis for the more targeted application of LLC in preclinical research.

Comprehensive Genomic and Transcriptomic Analysis of Sclerosing Pneumocytoma

Sclerosing pneumocytoma is a rare and distinct lung neoplasm whose histogenesis and molecular alterations are the subject of ongoing research. Our recent study revealed that AKT1 internal tandem duplications (ITD), point mutations, and short indels were present in almost all tested sclerosing pneumocytomas, suggesting that AKT1 mutations are a major driving oncogenic event in this tumor. Although the pathogenic role of AKT1 point mutations is well established, the significance of AKT1 ITD in oncogenesis remains largely unexplored. We conducted comprehensive genomic and transcriptomic analyses of sclerosing pneumocytoma to address this knowledge gap. RNA-sequencing data from 23 tumors and whole-exome sequencing data from 44 tumors were used to obtain insights into their genetic and transcriptomic profiles. Our analysis revealed a high degree of genetic and transcriptomic similarity between tumors carrying AKT1 ITD and those with AKT1 point mutations. Mutational signature analysis revealed COSMIC signatures 1 and 5 as the prevailing signatures of sclerosing pneumocytoma, associated with the spontaneous deamination of 5-methylcytosine and an unknown etiology, respectively. RNA-sequencing data analysis revealed that the sclerosing pneumocytoma gene expression profile is characterized by activation of the PI3K/AKT/mTOR pathway, which exhibits significant similarity between tumors harboring AKT1 ITD and those with AKT1 point mutations. Notably, an upregulation of SOX9, a transcription factor known for its involvement in fetal lung development, was observed in sclerosing pneumocytoma. Specifically, SOX9 expression was prominent in the round cell component, whereas it was relatively lower in the surface cell component of the tumor. To the best of our knowledge, this is the first comprehensive investigation of the genomic and transcriptomic characteristics of sclerosing pneumocytoma. Results of the present study provide insights into the molecular attributes of sclerosing pneumocytoma and a basis for future studies of this enigmatic tumor.

Whole-genome high-fidelity sequencing: A novel approach to detecting and characterization of mutagenicity in vivo.

Direct DNA sequencing can be used for characterizing mutagenicity in simple and complex biological models. Recently we described a method of whole-genome sequencing for detecting mutations in simple models of cultured bacteria, mammalian cells, and nematode. In the current proof-of-concept study, we expand and improve our method for evaluating a more complex mammalian biological model in outbred mice. We detail the method by applying it to a small set of animals treated with a mutagen with known mutagenicity profiles, N-ethyl-N-nitrosourea (ENU), for consistency with the known data. Whole-genome high-fidelity sequencing (HiFi Sequencing) showed frequencies and spectra of background mutations in tissues of untreated mice that were consistent with normal ageing and characterized by spontaneous or enzymatic deamination of 5-methylcytosine. In mice treated with a single 40mg/kg dose of ENU, the frequency of mutations in the genomic DNA of solid tissues increased up to 7-fold, with the greatest increase observed in the spleen and the smallest increase in the liver. The most common mutations detected in ENU-treated mice were T>A transitions and T>C transversions, consistent with the types of mutations caused by alkylating agents. The data suggest that HiFi Sequencing may be useful for characterizing mutagenicity of novel compounds in various biological models.

Deficiency of the Arabidopsis mismatch repair MSH6 attenuates Pseudomonas syringae invasion

The MutS homolog 6 (MSH6) is a nuclear DNA mismatch repair (MMR) gene that encodes the MSH6 protein. MSH6 interacts with MSH2 to form the MutSα heterodimer. MutSα corrects DNA mismatches and unpaired nucleotides arising during DNA replication, deamination of 5-methylcytosine, and recombination between non-identical DNA sequences. In addition to correcting DNA biosynthetic errors, MutSα also recognizes chemically damaged DNA bases. Here, we show that inactivation of MSH6 affects the basal susceptibility of Arabidopsis thaliana to Pseudomonas syringae pv tomato DC3000. The msh6 T-DNA insertional mutant exhibited a reduced susceptibility to the bacterial invasion. This heightened basal resistance of msh6 mutants appears to be dependent on an increased stomatal closure, an accumulation of H2O2 and double-strand breaks (DSBs) and a constitutive expression of pathogenesis-related (NPR1 and PR1) and DNA damage response (RAD51D and SOG1) genes. Complementation of this mutant with the MSH6 wild type allele under the control of its own promoter resulted in reversal of the basal bacterial resistance phenotype and the stomatal closure back to wild type levels. Taken together, these results demonstrate that inactivation of MSH6 increases Arabidopsis basal susceptibility to the bacterial pathogen and suggests a link between DNA repair and stress signaling in plants.

EBS-seq: enrichment-based method for accurate analysis of 5-hydroxymethylcytosine at single-base resolution

BackgroundA growing body of research has emphasized 5-hydroxymethylcytosine (5hmC) as an important epigenetic mark. High-resolution methods to detect 5hmC require high sequencing depth and are therefore expensive. Many studies have used enrichment-based methods to detect 5hmC; however, conventional enrichment-based methods have limited resolution. To overcome these limitations, we developed EBS-seq, a cost-efficient method for 5hmC detection with single-base resolution that combines the advantages of high-resolution methods and enrichment-based methods.ResultsEBS-seq uses selective labeling of 5hmC, deamination of cytosine and 5-methylcytosine, pull-down of labeled 5hmC, and C-to-T conversion during DNA amplification. Using this method, we profiled 5hmC in HEK293T cells and two colorectal cancer samples. Compared with conventional enrichment-based 5hmC detection, EBS-seq improved 5hmC signals by localizing them at single-base resolution. Furthermore, EBS-seq was able to determine 5hmC levels in CpG-dense regions where distortion of signals can occur, such as CpG islands and CpG shores. Comparing EBS-seq and conventional high-resolution 5hmC detection by ACE-seq, we showed that EBS-seq is more effective at finding 5hmC sites. Using EBS-seq, we found strong associations between gene expression and gene-body 5hmC content in both HEK293T cells and colorectal cancer samples.ConclusionsEBS-seq is a reliable and cost-efficient method for 5hmC detection because it simultaneously enriches 5hmC-containing DNA fragments and localizes 5hmC signals at single-base resolution. This method is a promising choice for 5hmC detection in challenging clinical samples with low 5hmC levels, such as cancer tissues.

Cytidine analogs are synthetic lethal with base excision repair default due to MBD4 deficiency

Inactivating mutations of MBD4 have been reported in subsets of various tumors. A deficiency of this DNA glycosylase, recognizing specifically T:G mismatch resulting from the deamination of methyl-cytosine, results in a hypermutated phenotype due to the accumulation of CpG>TpG transitions. Here, we hypothesize that the difference in DNA metabolism consecutive to MBD4 deficiency may result in specific cytotoxicities in MBD4-deficient tumor cells in a synthetic lethality fashion. After a large-scale drug repurposing screen, we show in two isogenic MBD4 knock-out cell models that the inactivation of MBD4 sensitizes cancer cells to cytidine analogs. We further confirm the exquisite activity of gemcitabine in an MBD4-deficient co-clinical model as (i) it completely prevented the development of an MBD4-deficient uveal melanoma patient-derived xenograft and (ii) treatment in the corresponding patient resulted in an exceptional tumor response. These data suggest that patients harboring MBD4-deficient tumors may be treated efficiently by cytidine analogs.

Modulation of DNA Methylation/Demethylation Reactions Induced by Nutraceuticals and Pollutants of Exposome Can Promote a C > T Mutation in the Breast Cancer Predisposing Gene PALB2

Background: Deregulation of DNA methylation/demethylation reactions may be the source of C > T mutation via active deamination of 5-methylcytosine to thymine. Exposome, that is to say, the totality of exposures to which an individual is subjected during their life, can deregulate these reactions. Thus, one may wonder whether the exposome can induce C > T mutations in the breast cancer-predisposing gene PALB2. Methods: Our work is based on the exposure of MCF10A mammary epithelial cells to seven compounds of our exposome (folate, Diuron, glyphosate, PFOA, iron, zinc, and ascorbic acid) alone or in cocktail. The qMSRE and RMS techniques were used to study the impact of these exposures on the level of methylation and mutation of the PALB2 gene. Results: Here, we have found that exposome compounds (nutriments, ions, pollutants) promoting the cytosine methylation and the 5-methylcytosine deamination have the ability to promote a specific C > T mutation in the PALB2 gene. Interestingly, we also noted that the addition of exposome compounds promoting the TET-mediated conversion of 5-methylcytosine (Ascorbic acid and iron) abrogates the presence of C > T mutation in the PALB2 gene. Conclusions: Our study provides a proof of concept supporting the idea that exposomes can generate genetic mutation by affecting DNA methylation/demethylation.

The mutational signatures of formalin fixation on the human genome

Clinical archives of patient material near-exclusively consist of formalin-fixed and paraffin-embedded (FFPE) blocks. The ability to precisely characterise mutational signatures from FFPE-derived DNA has tremendous translational potential. However, sequencing of DNA derived from FFPE material is known to be riddled with artefacts. Here we derive genome-wide mutational signatures caused by formalin fixation. We show that the FFPE-signature is highly similar to signature 30 (the signature of Base Excision Repair deficiency due to NTHL1 mutations), and chemical repair of DNA lesions leads to a signature highly similar to signature 1 (clock-like signature due to spontaneous deamination of methylcytosine). We demonstrate that using uncorrected mutational catalogues of FFPE samples leads to major mis-assignment of signature activities. To correct for this, we introduce FFPEsig, a computational algorithm to rectify the formalin-induced artefacts in the mutational catalogue. We demonstrate that FFPEsig enables accurate mutational signature analysis both in simulated and whole-genome sequenced FFPE cancer samples. FFPEsig thus provides an opportunity to unlock additional clinical potential of archival patient tissues.

Comprehensive Analysis of Genomic Alterations in Hepatoid Adenocarcinoma of the Stomach and Identification of Clinically Actionable Alterations.

Simple SummaryHepatoid adenocarcinoma of the stomach (HAS) is a subset of gastric cancer (GC) histologically characterized by hepatocellular carcinoma-like foci with or without alpha-fetoprotein secretion, which is easily misdiagnosed. Genomic alterations and potential targets for this population are still largely unknown. Additionally, treatment regimens of HAS are mainly based on GC guidelines, which is not reasonable for diseases with great heterogeneity. The present study comprehensively depicts the genomic features of HAS, and they are significantly different from GC, AFP-producing GC (AFPGC), and liver hepatocellular carcinoma (LIHC). Multiple aggressive behavior-related amplificated or deleted regions in HAS are firstly reported. Moreover, reliable and practicable clinically actionable alterations for HAS are identified, providing evidence for making personalized therapy based on the genomic characteristics of HAS instead of GC.Hepatoid adenocarcinoma of the stomach (HAS) is a rare malignancy with aggressive biological behavior. This study aimed to compare the genetic landscape of HAS with liver hepatocellular carcinoma (LIHC), gastric cancer (GC), and AFP-producing GC (AFPGC) and identify clinically actionable alterations. Thirty-eight cases of HAS were collected for whole-exome sequencing. Significantly mutated genes were identified. TP53 was the most frequently mutated gene (66%). Hypoxia, TNF-α/NFκB, mitotic spindle assembly, DNA repair, and p53 signaling pathways mutated frequently. Mutagenesis mechanisms in HAS were associated with spontaneous or enzymatic deamination of 5-methylcytosine to thymine and defective homologous recombination-related DNA damage repair. However, LIHC was characteristic of exposure to aflatoxin and aristolochic acid. The copy number variants (CNVs) in HAS was significantly different compared to LIHC, GC, and AFPGC. Aggressive behavior-related CNVs were identified, including local vascular invasion, advanced stages, and adverse prognosis. In 55.26% of HAS patients there existed at least one clinically actionable alteration, including ERBB2, FGFR1, CDK4, EGFR, MET, and MDM2 amplifications and BRCA1/2 mutations. MDM2 amplification with functional TP53 was detected in 5% of HAS patients, which was proved sensitive to MDM2 inhibitors. A total of 10.53% of HAS patients harbored TMB > 10 muts/Mb. These findings improve our understanding of the genomic features of HAS and provide potential therapeutic targets.

The mismatch-repair proteins MSH2 and MSH6 interact with the imprinting control regions through the ZFP57-KAP1 complex

BackgroundImprinting Control Regions (ICRs) are CpG-rich sequences acquiring differential methylation in the female and male germline and maintaining it in a parental origin-specific manner in somatic cells. Despite their expected high mutation rate due to spontaneous deamination of methylated cytosines, ICRs show conservation of CpG-richness and CpG-containing transcription factor binding sites in mammalian species. However, little is known about the mechanisms contributing to the maintenance of a high density of methyl CpGs at these loci.ResultsTo gain functional insights into the mechanisms for maintaining CpG methylation, we sought to identify the proteins binding the methylated allele of the ICRs by determining the interactors of ZFP57 that recognizes a methylated hexanucleotide motif of these DNA regions in mouse ESCs. By using a tagged approach coupled to LC–MS/MS analysis, we identified several proteins, including factors involved in mRNA processing/splicing, chromosome organization, transcription and DNA repair processes. The presence of the post-replicative mismatch-repair (MMR) complex components MSH2 and MSH6 among the identified ZFP57 interactors prompted us to investigate their DNA binding profile by chromatin immunoprecipitation and sequencing. We demonstrated that MSH2 was enriched at gene promoters overlapping unmethylated CpG islands and at repeats. We also found that both MSH2 and MSH6 interacted with the methylated allele of the ICRs, where their binding to DNA was mediated by the ZFP57/KAP1 complex.ConclusionsOur findings show that the MMR complex is concentrated on gene promoters and repeats in mouse ESCs, suggesting that maintaining the integrity of these regions is a primary function of highly proliferating cells. Furthermore, the demonstration that MSH2/MSH6 are recruited to the methylated allele of the ICRs through interaction with ZFP57/KAP1 suggests a role of the MMR complex in the maintenance of the integrity of these regulatory regions and evolution of genomic imprinting in mammalian species.

Characterization of mutational signatures of homologous recombination deficiency (HRD) in Chinese breast cancer.

e15129 Background: Homologous recombination deficiency (HRD) is a common molecular characteristic of genomic instability and has been proved to be an important biomarker for personalized therapy decisions and stratified management in breast cancer. In this study, we explored the mutational features of HRD tumors identified by HRD score in Chinese breast cancer. Methods: Tumor tissue DNA from 68 breast cancer patients was isolated and sequenced by next-generation sequencing (NGS) technologies with AcornMed 808-genes panel. HRD score was calculated as an unweighted sum of loss of heterozygosity (LOH) score, telomeric allelic imbalance (TAI) score and large-scale state transition (LST) score by AcornHRD algorithm, and mutations in 39 homologous recombination (HR) related genes including BRCA1/2 were identified for each case. Results: Based on the mutation alternation profiles from our cohort, somatic BRCA1/2 mutations were present in 15 patients (22.1%) while 29 patients (42.6%) harbored one or more mutations in HR-related genes excluding BRCA1/2. Germline mutations of BRCA1/2 and other HR-related genes occurred in 11.8% and 4.4% of cases detected, respectively.In our cohort, the median HRD score was 40 (range from 3 to 92). Sixty-eight cases were divided into HRD-high or HRD-low group based on a cut-off value of ≥ 42 or < 42, with HRD-high group accounting for approximately 47% of our cohort, indicating the combination with HRD score yielded an increment of 26.5% in coverage of HRD within BRCA1/2 wildtype cases. The HRD-high group covered 82.4% of the BRCA-mutated cases and also showed higher mutation level of BRCA1/2-excluded HR genes than HRD-low group (p = 0.001). Prevalence of mutations in a total of 12 genes, including BRCA1(31.3% vs 5.6%, p = 0.009) and NBN (25% vs 5.6%, p = 0.038), were significant difference between HRD-high and HRD-low group. Besides, Cosmic Signature 1 (spontaneous deamination of 5-methylcytosine) and 5 (unknown) were existed in both groups, whereas Signature 6 (defective DNA mismatch repair) was only discovered in HRD-high group. Compared with HRD-low ones, copy number variations ratio was higher in HRD-high group (p = 0.031). Meanwhile, the patients with ERBB2 amplification (HER2+) showed a lower level of HRD score than HER2- ones. Conclusions: This study comprehensively delineated the genetic characteristics of HRD tumors, which suggests that HRD status marked by HRD score is strongly associated with BRCA1/2 and HR-related genes mutation as well as other genomic instability signatures, thus indicating its predictive value as a promising biomarker for DNA repair deficiency related drugs or therapies in Chinese breast cancer.

Somatic mutation rates scale with lifespan across mammals

The rates and patterns of somatic mutation in normal tissues are largely unknown outside of humans1–7. Comparative analyses can shed light on the diversity of mutagenesis across species, and on long-standing hypotheses about the evolution of somatic mutation rates and their role in cancer and ageing. Here we performed whole-genome sequencing of 208 intestinal crypts from 56 individuals to study the landscape of somatic mutation across 16 mammalian species. We found that somatic mutagenesis was dominated by seemingly endogenous mutational processes in all species, including 5-methylcytosine deamination and oxidative damage. With some differences, mutational signatures in other species resembled those described in humans8, although the relative contribution of each signature varied across species. Notably, the somatic mutation rate per year varied greatly across species and exhibited a strong inverse relationship with species lifespan, with no other life-history trait studied showing a comparable association. Despite widely different life histories among the species we examined—including variation of around 30-fold in lifespan and around 40,000-fold in body mass—the somatic mutation burden at the end of lifespan varied only by a factor of around 3. These data unveil common mutational processes across mammals, and suggest that somatic mutation rates are evolutionarily constrained and may be a contributing factor in ageing.

Gigantic genomes of salamanders indicate that body temperature, not genome size, is the driver of global methylation and 5-methylcytosine deamination in vertebrates.

Transposable elements (TEs) are sequences that replicate and move throughout genomes, and they can be silenced through methylation of cytosines at CpG dinucleotides. TE abundance contributes to genome size, but TE silencing variation across genomes of different sizes remains underexplored. Salamanders include most of the largest C-values - 9 to 120 Gb. We measured CpG methylation levels in salamanders with genomes ranging from 2N = ∼58 Gb to 4N = ∼116 Gb. We compared these levels to results from endo- and ectothermic vertebrates with more typical genomes. Salamander methylation levels are approximately 90%, higher than all endotherms. However, salamander methylation does not differ from other ectotherms, despite an approximately 100-fold difference in nuclear DNA content. Because methylation affects the nucleotide compositional landscape through 5-methylcytosine deamination to thymine, we quantified salamander CpG dinucleotide levels and compared them to other vertebrates. Salamanders and other ectotherms have comparable CpG levels, and ectotherm levels are higher than endotherms. These data show no shift in global methylation at the base of salamanders, despite a dramatic increase in TE load and genome size. This result is reconcilable with previous studies that considered endothermy and ectothermy, which may be more important drivers of methylation in vertebrates than genome size.