Mutagenesis at methylated CpG sequences

The evolutionary consequences of interactions between the epigenome, the genome and the environment

The epigenome is the suite of interacting chemical marks and molecules that helps to shape patterns of development, phenotypic plasticity and gene regulation, in part due to its responsiveness to environmental stimuli. There is increasing interest in understanding the functional and evolutionary importance of this sensitivity under ecologically realistic conditions. Observations that epigenetic variation abounds in natural populations have prompted speculation that it may facilitate evolutionary responses to rapid environmental perturbations, such as those occurring under climate change. A frequent point of contention is whether epigenetic variants reflect genetic variation or are independent of it. The genome and epigenome often appear tightly linked and interdependent. While many epigenetic changes are genetically determined, the converse is also true, with DNA sequence changes influenced by the presence of epigenetic marks. Understanding how the epigenome, genome and environment interact with one another is therefore an essential step in explaining the broader evolutionary consequences of epigenomic variation. Drawing on results from experimental and comparative studies carried out in diverse plant and animal species, we synthesize our current understanding of how these factors interact to shape phenotypic variation in natural populations, with a focus on identifying similarities and differences between taxonomic groups. We describe the main components of the epigenome and how they vary within and between taxa. We review how variation in the epigenome interacts with genetic features and environmental determinants, with a focus on the role of transposable elements (TEs) in integrating the epigenome, genome and environment. And we look at recent studies investigating the functional and evolutionary consequences of these interactions. Although epigenetic differentiation in nature is likely often a result of drift or selection on stochastic epimutations, there is growing evidence that a significant fraction of it can be stably inherited and could therefore contribute to evolution independently of genetic change.

Methylation and hydroxymethylation of cytosine alter activity and fidelity of translesion DNA polymerases

Epigenetic cytosine methylation covers most of genomic CpG dinucleotides in human cells. In addition to common deamination-mediated mutagenesis at CpG sites, an alternative deamination-independent pathway associated with DNA polymerase activity was previously described. This mutagenesis is characterized by the TCG→TTG mutational signature and is believed to arise from dAMP misincorporation opposite 5-methylcytosine (mC) or its oxidized derivative 5-hydroxymethylcytosine (hmC) by B-family replicative DNA polymerases with disrupted proofreading 3→5'-exonuclease activity. In addition to being less stable and pro-mutagenic themselves, cytosine modifications also increase the risk of adjacent nucleotides damage, including the formation of 8-oxo-2'-deoxyguanosine (8-oxoG), a well-known mutagenic lesion. The effect of cytosine methylation on error-prone DNA polymerases lacking proofreading activity and involved in repair and DNA translesion synthesis remains unexplored. Here we analyze the efficiency and fidelity of translesion Y-family polymerases (Pol κ, Pol η, Pol ι and REV1) and primase-polymerase PrimPol opposite mC and hmC as well as opposite 8-oxoG adjacent to mC in the TCG context. We demonstrate that epigenetic cytosine modifications suppress Pol ι and REV1 activities and lead to increasing dAMP misincorporation by PrimPol, Pol κ and Pol ι in vitro. Cytosine methylation also increases misincorporation of dAMP opposite the adjacent 8-oxoG by PrimPol, decreases the TLS activity of Pol η opposite the lesion but increases dCMP incorporation opposite 8-oxoG by REV1. Altogether, these data suggest that methylation and hydroxymethylation of cytosine alter activity and fidelity of translesion DNA polymerases.

Biochemical reconstitution of a major age-related cancer mutational signature by heat-induced spontaneous deamination of 5-methylcytosine residues, repair of uracil residues, and DNA replication

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 ε. These results provide biochemical evidence to support that the majority of mutations in cancers are produced by ordinary DNA replication on spontaneously damaged DNA.

Analysis of gliomas DNA methylation: Assessment of pre-analytical variables

Precision oncology is driven by molecular biomarkers. For glioblastoma multiforme (GBM), the most common malignant adult primary brain tumor, O6-methylguanine-DNA methyltransferase ( MGMT ) gene DNA promoter methylation is an important prognostic and treatment clinical biomarker. Time consuming pre-analytical steps such as biospecimen storage before fixing, sampling, and processing are major sources of errors and batch effects, that are further confounded by intra-tumor heterogeneity of MGMT promoter methylation. To assess the effect of pre-analytical variables on GBM DNA methylation, tissue storage/sampling (CryoGrid), sample preparation multi-sonicator (PIXUL) and 5-methylcytosine (5mC) DNA immunoprecipitation (Matrix MeDIP-qPCR/seq) platforms were used. MGMT promoter CpG methylation was examined in 173 surgical samples from 90 individuals, 50 of these were used for intra-tumor heterogeneity studies. MGMT promoter methylation levels in paired frozen and formalin fixed paraffin embedded (FFPE) samples were very close, confirming suitability of FFPE for MGMT promoter methylation analysis in clinical settings. Matrix MeDIP-qPCR yielded similar results to methylation specific PCR (MS-PCR). Warm ex-vivo ischemia (37°C up to 4hrs) and 3 cycles of repeated sample thawing and freezing did not alter 5mC levels at MGMT promoter, exon and upstream enhancer regions, demonstrating the resistance of DNA methylation to the most common variations in sample processing conditions that might be encountered in research and clinical settings. 20-30% of specimens exhibited intratumor heterogeneity in the MGMT DNA promoter methylation. Collectively these data demonstrate that variations in sample fixation, ischemia duration and temperature, and DNA methylation assay technique do not have significant impact on assessment of MGMT promoter methylation status. However, intratumor methylation heterogeneity underscores the need for histologic verification and value of multiple biopsies at different GBM geographic tumor sites in assessment of MGMT promoter methylation. Matrix-MeDIP-seq analysis revealed that MGMT promoter methylation status clustered with other differentially methylated genomic loci (e.g. HOXA and lncRNAs), that are likewise resilient to variation in above post-resection pre-analytical conditions. These MGMT -associated global DNA methylation patterns offer new opportunities to validate more granular data-based epigenetic GBM clinical biomarkers where the CryoGrid-PIXUL-Matrix toolbox could prove to be useful.

BRCA1, BRCA2, and TP53 germline and somatic variants and clinicopathological characteristics of Brazilian patients with epithelial ovarian cancer

BACKGROUND

Approximately 3/4 of ovarian cancers are diagnosed in advanced stages, with the high-grade epithelial ovarian carcinoma (EOC) accounting for 90% of the cases. EOC present high genomic instability and somatic loss-of-function variants in genes associated with homologous recombination mutational repair pathway (HR), such as BRCA1 and BRCA2, and in TP53. The identification of germline variants in HR genes in EOC is relevant for treatment of platinum resistant tumors and relapsed tumors with therapies based in synthetic lethality such as PARP inhibitors. Patients with somatic variants in HR genes may also benefit from these therapies. In this work was analyzed the frequency of somatic variants in BRCA1, BRCA2, and TP53 in an EOC cohort of Brazilian patients, estimating the proportion of variants in tumoral tissue and their association with progression-free survival and overall survival.

METHODS

The study was conducted with paired blood/tumor samples from 56 patients. Germline and tumoral sequences of BRCA1, BRCA2, and TP53 were obtained by massive parallel sequencing. The HaplotypeCaller method was used for calling germline variants, and somatic variants were called with Mutect2.

RESULTS

A total of 26 germline variants were found, and seven patients presented germline pathogenic or likely pathogenic variants in BRCA1 or BRCA2. The analysis of tumoral tissue identified 52 somatic variants in 41 patients, being 43 somatic variants affecting or likely affecting protein functionality. Survival analyses showed that tumor staging was associated with overall survival (OS), while the presence of somatic mutation in TP53 was not associated with OS or progression-free survival.

CONCLUSION

Frequency of pathogenic or likely pathogenic germline variants in BRCA1 and BRCA2 (12.5%) was lower in comparison with other studies. TP53 was the most altered gene in tumors, with 62.5% presenting likely non-functional or non-functional somatic variants, while eight 14.2% presented likely non-functional or non-functional somatic variants in BRCA1 or BRCA2.

Different gene alterations in patients with non-small-cell lung cancer between the eastern and southern China

Geographical differences are conspicuous in lung cancer, and the distinct molecular features of lung tumor between Western patients and Asian patients have been demonstrated. However, the etiology of non-small-cell lung cancer (NSCLC) and the distribution of associated molecular aberrations in China have not been fully elucidated. The mutational profiles of 12 lung cancer-related genes were investigated in 85 patients from eastern China and 88 patients from southern China who had been histologically confirmed NSCLC. Overall, 93.6% (162/173) of tumor samples harbored at least one somatic alteration. The most frequently mutated genes were TP53 (56.1%), EGFR (50.3%), and KRAS (14.5%). We found that EGFR mutated much more frequently (60.0% vs 40.9%, P = 0.012) and TP53 mutations had significantly lower incidence (47.1% vs 64.8%, P = 0.019) in eastern cohort than that in southern cohort. Mutational signature analysis revealed a region-related mutagenesis mechanism characterized by a high prevalence of C to T transitions mainly occurring at CpG dinucleotides in southern patients. This study reveals the difference in the mutational features between NSCLC patients in eastern and southern China. The distinct patterns of gene mutation could provide clues for the mechanism of carcinogenesis of lung cancer, offering opportunities to stratify patients into optimal treatment plans based on genomic profiles.

Environmental carcinogens disproportionally mutate genes implicated in neurodevelopmental disorders

Introduction

De novo mutations contribute to a large proportion of sporadic psychiatric and developmental disorders, yet the potential role of environmental carcinogens as drivers of causal de novo mutations in neurodevelopmental disorders is poorly studied.

Methods

To explore environmental mutation vulnerability of disease-associated gene sets, we analyzed publicly available whole genome sequencing datasets of mutations in human induced pluripotent stem cell clonal lines exposed to 12 classes of environmental carcinogens, and human lung cancers from individuals living in highly polluted regions. We compared observed rates of exposure-induced mutations in disease-related gene sets with the expected rates of mutations based on control genes randomly sampled from the genome using exact binomial tests. To explore the role of sequence characteristics in mutation vulnerability, we modeled the effects of sequence length, gene expression, and percent GC content on mutation rates of entire genes and gene coding sequences using multivariate Quasi-Poisson regressions.

Results

We demonstrate that several mutagens, including radiation and polycyclic aromatic hydrocarbons, disproportionately mutate genes related to neurodevelopmental disorders including autism spectrum disorders, schizophrenia, and attention deficit hyperactivity disorder. Other disease genes including amyotrophic lateral sclerosis, Alzheimer's disease, congenital heart disease, orofacial clefts, and coronary artery disease were generally not mutated more than expected. Longer sequence length was more strongly associated with elevated mutations in entire genes compared with mutations in coding sequences. Increased expression was associated with decreased coding sequence mutation rate, but not with the mutability of entire genes. Increased GC content was associated with increased coding sequence mutation rates but decreased mutation rates in entire genes.

Discussion

Our findings support the possibility that neurodevelopmental disorder genetic etiology is partially driven by a contribution of environment-induced germ line and somatic mutations.

The Mutagenic Consequences of DNA Methylation within and across Generations

DNA methylation is an epigenetic modification with wide-ranging consequences across the life of an organism. This modification can be stable, persisting through development despite changing environmental conditions. However, in other contexts, DNA methylation can also be flexible, underlying organismal phenotypic plasticity. One underappreciated aspect of DNA methylation is that it is a potent mutagen; methylated cytosines mutate at a much faster rate than other genetic motifs. This mutagenic property of DNA methylation has been largely ignored in eco-evolutionary literature, despite its prevalence. Here, we explore how DNA methylation induced by environmental and other factors could promote mutation and lead to evolutionary change at a more rapid rate and in a more directed manner than through stochastic genetic mutations alone. We argue for future research on the evolutionary implications of DNA methylation driven mutations both within the lifetime of organisms, as well as across timescales.

Single Amino Acid Substitution the DNA Repairing Gene Radiation-Sensitive 4 Contributes to Ultraviolet Tolerance of a Plant Pathogen

To successfully survive and reproduce, all species constantly modify the structure and expression of their genomes to cope with changing environmental conditions including ultraviolet (UV) radiation. Thus, knowledge of species adaptation to environmental changes is a central theme of evolutionary studies which could have important implication for disease management and social-ecological sustainability in the future but is generally insufficient. Here, we investigated the evolution of UV adaptation in organisms by population genetic analysis of sequence structure, physiochemistry, transcription, and fitness variation in the radiation-sensitive 4 (RAD4) gene of the Irish potato famine pathogen Phytophthora infestans sampled from various altitudes. We found that RAD4 is a key gene determining the resistance of the pathogen to UV stress as indicated by strong phenotype-genotype-geography associations and upregulated transcription after UV exposure. We also found conserved evolution in the RAD4 gene. Only five nucleotide haplotypes corresponding to three protein isoforms generated by point mutations were detected in the 140 sequences analyzed and the mutations were constrained to the N-terminal domain of the protein. Physiochemical changes associated with non-synonymous mutations generate severe fitness penalty to mutants, which are purged out by natural selection, leading to the conserved evolution observed in the gene.

Genomic alterations of dermatofibrosarcoma protuberans revealed by whole-genome sequencing

Background

Dermatofibrosarcoma protuberans (DFSP) is a rare and marginal cutaneous sarcoma of intermediate‐grade malignancy, for which the genomic landscape remains unclear. Understanding the landscape of DFSP will help to further classify the genomic pathway of malignant development in soft tissue.

Objectives

To identify the comprehensive molecular pathogenesis of DFSP.

Methods

In this study, the comprehensive genomic features, with 53 tumour‐normal pairs of DFSP, were revealed by whole‐genome sequencing.

Results

The mutational signature 1 (C > T mutation at CpG dinucleotides) is featured in DFSP, resulting in higher mutations in DNA replication. Interestingly, the recurrence of DFSP is correlated with low tumour mutation burden. Novel mutation genes in DFSP were identified, including MUC4/6, KMT2C and BRCA1, and subsequently, three molecular subtypes of DFSP were classified on the basis of MUC4 and MUC6 mutations. Various structural aberrations including genomic rearrangements were identified in DSFPs, particularly in 17q and 22q, which cause oncogene amplification (AKT1, SPHK1, COL1A1, PDGFβ) or tumour suppressor deletion (CDKN2A/B). In addition to gene fusion of COL1A1‐PDGFβ [t(17;22)], we identified gene fusion of SLC2A5‐BTBD7 [t(1;14)] in DFSP through whole‐genome sequencing, and verified it experimentally. Enrichment analysis of altered molecules revealed that DNA repair, cell cycle, phosphoinositide 3‐kinase and Janus kinase pathways were primarily involved in DFSP.

Conclusions

This is the first large‐scale whole‐genome sequencing for DFSP, and our findings describe the comprehensive genomic landscape, highlighting the molecular complexity and genomic aberrations of DFSP. Our findings also provide novel potential diagnostic and therapeutic targets for this disease.

What is already known about this topic?

Chromosomal translocation between chromosome 17 and chromosome 22 is the main feature in the pathogenesis of dermatofibrosarcoma protuberans (DFSP).

What does this study add?

We describe the comprehensive genomic landscape of DFSP, highlighting the molecular complexity and genomic aberrations.

Our findings provide novel potential diagnostic and therapeutic targets for this disease.

What is the translational message?

Our study revealed novel molecular subtypes of DFSP based on genetic mutations, which benefits precision diagnosis.

We also found oncogene amplification, including AKT1 and SPHK1, which provides novel potential target molecules for further DFSP treatment.

In addition to gene fusion of COL1A1‐PDGFβ, we identified a novel gene fusion of SLC2A5‐BTBD7 in DFSP, which is a novel potential diagnostic and therapeutic target for this disease.

Species-Specific Paternal Age Effects and Sperm Methylation Levels of Developmentally Important Genes

A growing number of sperm methylome analyses have identified genomic loci that are susceptible to paternal age effects in a variety of mammalian species, including human, bovine, and mouse. However, there is little overlap between different data sets. Here, we studied whether or not paternal age effects on the sperm epigenome have been conserved in mammalian evolution and compared methylation patterns of orthologous regulatory regions (mainly gene promoters) containing both conserved and non-conserved CpG sites in 94 human, 36 bovine, and 94 mouse sperm samples, using bisulfite pyrosequencing. We discovered three (NFKB2, RASGEF1C, and RPL6) age-related differentially methylated regions (ageDMRs) in humans, four (CHD7, HDAC11, PAK1, and PTK2B) in bovines, and three (Def6, Nrxn2, and Tbx19) in mice. Remarkably, the identified sperm ageDMRs were all species-specific. Most ageDMRs were in genomic regions with medium methylation levels and large methylation variation. Orthologous regions in species not showing this age effect were either hypermethylated (>80%) or hypomethylated (<20%). In humans and mice, ageDMRs lost methylation, whereas bovine ageDMRs gained methylation with age. Our results are in line with the hypothesis that sperm ageDMRs are in regions under epigenomic evolution and may be part of an epigenetic mechanism(s) for lineage-specific environmental adaptations and provide a solid basis for studies on downstream effects in the genes analyzed here.

Somatic mutational profiles and germline polygenic risk scores in human cancer

BACKGROUND

The mutational profile of cancer reflects the activity of the mutagenic processes which have been operative throughout the lineage of the cancer cell. These processes leave characteristic profiles of somatic mutations called mutational signatures. Mutational signatures, including single-base substitution (SBS) signatures, may reflect the effects of exogenous or endogenous exposures.

METHODS

We used polygenic risk scores (PRS) to summarize common germline variation associated with cancer risk and other cancer-related traits and examined the association between somatic mutational profiles and germline PRS in 12 cancer types from The Cancer Genome Atlas. Somatic mutational profiles were constructed from whole-exome sequencing data of primary tumors. PRS were calculated for the 12 selected cancer types and 9 non-cancer traits, including cancer risk determinants, hormonal factors, and immune-mediated inflammatory diseases, using germline genetic data and published summary statistics from genome-wide association studies.

RESULTS

We found 17 statistically significant associations between somatic mutational profiles and germline PRS after Bonferroni correction (p < 3.15 × 10^-5), including positive associations between germline inflammatory bowel disease PRS and number of somatic mutations attributed to signature SBS1 in prostate cancer and APOBEC-related signatures in breast cancer. Positive associations were also found between age at menarche PRS and mutation counts of SBS1 in overall and estrogen receptor-positive breast cancer. Consistent with prior studies that found an inverse association between the pubertal development PRS and risk of prostate cancer, likely reflecting hormone-related mechanisms, we found an inverse association between age at menarche PRS and mutation counts of SBS1 in prostate cancer. Inverse associations were also found between several cancer PRS and tumor mutation counts.

CONCLUSIONS

Our analysis suggests that there are robust associations between tumor somatic mutational profiles and germline PRS. These may reflect the mechanisms through hormone regulation and immune responses that contribute to cancer etiology and drive cancer progression.

LONG/TERM GENETIC AND EPIGENETIC DISORDERS IN PERSONS EXPOSED TO IONIZING RADIATION AND THEIR DESCENDANTS (review)

The review is devoted to long-term genetic and epigenetic disorders in exposed individuals and their descendants,namely to cytogenetic effects in the Chornobyl NPP accident clean-up workers and their children, DNA methylation as an epigenetic modification of human genome. Data presented in review expand the understanding of risk of the prolonged exposure for the present and future generations, which is one of key problems posed by fundamental radiation genetics and human radiobiology.

Lack of genotype-by-environment interaction suggests limited potential for evolutionary changes in plasticity in the eastern oyster, Crassostrea virginica

Eastern oysters in the northern Gulf of Mexico are facing rapid environmental changes and can respond to this change via plasticity or evolution. Plasticity can act as an immediate buffer against environmental change, but this buffering could impact the organism's ability to evolve in subsequent generations. While plasticity and evolution are not mutually exclusive, the relative contribution and interaction between them remains unclear. In this study, we investigate the roles of plastic and evolved responses to environmental variation and Perkinsus marinus infection in Crassostrea virginica by using a common garden experiment with 80 oysters from six families outplanted at two field sites naturally differing in salinity. We use growth data, P. marinus infection intensities, 3' RNA sequencing (TagSeq) and low-coverage whole-genome sequencing to identify the effect of genotype, environment and genotype-by-environment interaction on the oyster's response to site. As one of first studies to characterize the joint effects of genotype and environment on transcriptomic and morphological profiles in a natural setting, we demonstrate that C. virginica has a highly plastic response to environment and that this response is parallel among genotypes. We also find that genes responding to genotype have distinct and opposing profiles compared to genes responding to environment with regard to expression levels, Ka/Ks ratios and nucleotide diversity. Our findings suggest that C. virginica may be able to buffer the immediate impacts of future environmental changes by altering gene expression and physiology, but the lack of genetic variation in plasticity suggests limited capacity for evolved responses.

Characteristic mutations induced in the small intestine of Msh2-knockout gpt delta mice

Background

Base pair mismatches in genomic DNA can result in mutagenesis, and consequently in tumorigenesis. To investigate how mismatch repair deficiency increases mutagenicity under oxidative stress, we examined the type and frequency of mutations arising in the mucosa of the small intestine of mice carrying a reporter gene encoding guanine phosphoribosyltransferase (gpt) and in which the Msh2 gene, which encodes a component of the mismatch repair system, was either intact (Msh2+/+::gpt/0; Msh2-bearing) or homozygously knockout (KO) (Msh2−/−::gpt/0; Msh2-KO).

Results

Gpt mutant frequency in the small intestine of Msh2-KO mice was about 10 times that in Msh2-bearing mice. Mutant frequency in the Msh2-KO mice was not further enhanced by administration of potassium bromate, an oxidative stress inducer, in the drinking water at a dose of 1.5 g/L for 28 days. Mutation analysis showed that the characteristic mutation in the small intestine of the Msh2-KO mice was G-to-A transition, irrespective of whether potassium bromate was administered. Furthermore, administration of potassium bromate induced mutations at specific sites in gpt in the Msh2-KO mice: G-to-A transition was frequently induced at two known sites of spontaneous mutation (nucleotides 110 and 115, CpG sites) and at nucleotides 92 and 113 (3′-side of 5′-GpG-3′), and these sites were confirmed to be mutation hotspots in potassium bromate-administered Msh2-KO mice. Administration of potassium bromate also induced characteristic mutations, mainly single-base deletion and insertion of an adenine residue, in sequences of three to five adenine nucleotides (A-runs) in Msh2-KO mice, and elevated the overall proportion of single-base deletions plus insertions in Msh2-KO mice.

Conclusions

Our previous study revealed that administration of potassium bromate enhanced tumorigenesis in the small intestine of Msh2-KO mice and induced G-to-A transition in the Ctnnb1 gene. Based on our present and previous observations, we propose that oxidative stress under conditions of mismatch repair deficiency accelerates the induction of single-adenine deletions at specific sites in oncogenes, which enhances tumorigenesis in a synergistic manner with G-to-A transition in other oncogenes (e.g., Ctnnb1).

Supplementary Information

The online version contains supplementary material available at 10.1186/s41021-021-00196-0.

Whole Exome Sequencing reveals NOTCH1 mutations in anaplastic large cell lymphoma and points to Notch both as a key pathway and a potential therapeutic target

Patients diagnosed with Anaplastic Large Cell Lymphoma (ALCL) are still treated with toxic multi-agent chemotherapy and as many as 25-50% of patients relapse. To understand disease pathology and to uncover novel targets for therapy, Whole-Exome Sequencing (WES) of Anaplastic Lymphoma Kinase (ALK)+ ALCL was performed as well as Gene-Set Enrichment Analysis. This revealed that the T-cell receptor (TCR) and Notch pathways were the most enriched in mutations. In particular, variant T349P of NOTCH1, which confers a growth advantage to cells in which it is expressed, was detected in 12% of ALK+ and ALK- ALCL patient samples. Furthermore, we demonstrate that NPM-ALK promotes NOTCH1 expression through binding of STAT3 upstream of NOTCH1. Moreover, inhibition of NOTCH1 with γ-secretase inhibitors (GSIs) or silencing by shRNA leads to apoptosis; co-treatment in vitro with the ALK inhibitor Crizotinib led to additive/synergistic anti-tumour activity suggesting this may be an appropriate combination therapy for future use in the circumvention of ALK inhibitor resistance. Indeed, Crizotinib-resistant and sensitive ALCL were equally sensitive to GSIs. In conclusion, we show a variant in the extracellular domain of NOTCH1 that provides a growth advantage to cells and confirm the suitability of the Notch pathway as a second-line druggable target in ALK+ ALCL.

Altitudinal Heterogeneity of UV Adaptation in Phytophthorainfestans Is Associated with the Spatial Distribution of a DNA Repair Gene

Climate change is considered a major threat to society and nature. UV irradiation is the most important environmental genotoxic agent. Thus, how elevated UV irradiation may influence human health and ecosystems has generated wide concern in the scientific community, as well as with policy makers and the public in general. In this study, we investigated patterns and mechanisms of UV adaptation in natural ecosystems by studying a gene-specific variation in the potato late blight pathogen, Phytophthora infestans. We compared the sequence characteristics of radiation sensitive 23 (RAD23), a gene involved in the nucleotide excision repair (NER) pathway and UV tolerance, in P. infestans isolates sampled from various altitudes. We found that lower genetic variation in the RAD23 gene was caused by natural selection. The hypothesis that UV irradiation drives this selection was supported by strong correlations between the genomic characteristics and altitudinal origin (historic UV irradiation) of the RAD23 sequences with UV tolerance of the P. infestans isolates. These results indicate that the RAD23 gene plays an important role in the adaptation of P. infestans to UV stress. We also found that different climatic factors could work synergistically to determine the evolutionary adaptation of species, making the influence of climate change on ecological functions and resilience more difficult to predict. Future attention should aim at understanding the collective impact generated by simultaneous change in several climate factors on species adaptation and ecological sustainability, using state of the art technologies such as experimental evolution, genome-wide scanning, and proteomics.

The emergence of the brain non-CpG methylation system in vertebrates

Mammalian brains feature exceptionally high levels of non-CpG DNA methylation alongside the canonical form of CpG methylation. Non-CpG methylation plays a critical regulatory role in cognitive function, which is mediated by the binding of MeCP2, the transcriptional regulator that when mutated causes Rett syndrome. However, it is unclear whether the non-CpG neural methylation system is restricted to mammalian species with complex cognitive abilities or has deeper evolutionary origins. To test this, we investigated brain DNA methylation across 12 distantly related animal lineages, revealing that non-CpG methylation is restricted to vertebrates. We discovered that in vertebrates, non-CpG methylation is enriched within a highly conserved set of developmental genes transcriptionally repressed in adult brains, indicating that it demarcates a deeply conserved regulatory program. We also found that the writer of non-CpG methylation, DNMT3A, and the reader, MeCP2, originated at the onset of vertebrates as a result of the ancestral vertebrate whole-genome duplication. Together, we demonstrate how this novel layer of epigenetic information assembled at the root of vertebrates and gained new regulatory roles independent of the ancestral form of the canonical CpG methylation. This suggests that the emergence of non-CpG methylation may have fostered the evolution of sophisticated cognitive abilities found in the vertebrate lineage.

Supervised mutational signatures for obesity and other tissue-specific etiological factors in cancer

Determining the etiologic basis of the mutations that are responsible for cancer is one of the fundamental challenges in modern cancer research. Different mutational processes induce different types of DNA mutations, providing 'mutational signatures' that have led to key insights into cancer etiology. The most widely used signatures for assessing genomic data are based on unsupervised patterns that are then retrospectively correlated with certain features of cancer. We show here that supervised machine-learning techniques can identify signatures, called SuperSigs, that are more predictive than those currently available. Surprisingly, we found that aging yields different SuperSigs in different tissues, and the same is true for environmental exposures. We were able to discover SuperSigs associated with obesity, the most important lifestyle factor contributing to cancer in Western populations.

Prenatal and postnatal diagnosis of Schuurs-Hoeijmakers syndrome: Case series and review of the literature

Schuurs-Hoeijmakers syndrome (SHS) is a rare syndrome involving a de novo variant in the PACS1 gene on chromosome 11q13. There are 36 individuals published in the literature so far, mostly diagnosed postnatally (34/36) after recognizing the typical facial features co-occurring with developmental delay, intellectual disability, and multiple malformations. Herein, we present one prenatal and 15 postnatal cases with the recurrent heterozygous pathogenic variant NM_018026.3:c.607C>T p.(Arg203Trp) in the PACS1 gene detected by exome sequencing. These 16 cases were identified by mining Centogene and the Hong Kong clinical genetic service databases. Collectively, the 49 postnatally diagnosed individuals present with typical facial features and developmental delay, while the three prenatally diagnosed individuals present with multiple congenital anomalies. In the current study, the use of exome sequencing as an unbiased diagnostic tool aided the diagnosis of SHS (pre- and postnatally). The identification of additional cases with SHS add to the current understanding of the clinical phenotype associated with pathogenic PACS1 variants. Databases combining clinical and genetic information are helpful for the study of rare diseases.

Expanding the genetic spectrum of choroideremia in an Australian cohort: report of five novel CHM variants

Choroideremia is an X-linked chorioretinal dystrophy caused by mutations in the CHM gene. Several CHM gene replacement clinical trials are in advanced stages. In this study, we report the molecular confirmation of choroideremia in 14 Australian families sourced from the Australian Inherited Retinal Disease Registry and DNA Bank. Sixteen males (14 symptomatic) and 18 females (4 symptomatic; 14 obligate carriers) were identified for analysis. Participants' DNA was analyzed for disease-causing CHM variants by Sanger sequencing, TaqMan qPCR and targeted NGS. We report phenotypic and genotypic data for the 14 symptomatic males and four females manifesting disease symptoms. A pathogenic or likely pathogenic CHM variant was detected in all families. Eight variants were previously reported, and five were novel. Two de novo variants were identified. We previously reported the molecular confirmation of choroideremia in 11 Australian families. This study expands the CHM genetically confirmed Australian cohort to 32 males and four affected carrier females.

[Mutation signatures in head and neck squamous cell carcinoma : Pathogenesis and therapeutic potential]

BACKGROUND

The pathogenesis of head and neck squamous cell carcinoma (HNSCC) is a complex and multistage process which results from the interaction of exogenous and endogenous cellular processes. Each of these processes leaves a characteristic pattern of mutations on the tumor genome, a so-called mutational signature.

STATE OF THE ART

The subject of current studies is to decipher specific signatures of mutational processes operating during HNSCC pathogenesis and to address their prognostic value. Computational analysis of genomic sequencing data by The Cancer Genome Atlas (TCGA) revealed mutational signatures 1, 2, 4, 5, 7, and 13 as the main players in HNSCC pathogenesis. Signature 16 was first discovered in human papillomavirus (HPV)-negative oral and oropharyngeal tumors. In many studies, an association of signature 16 with alcohol and tobacco consumption as well as with an unfavorable prognosis was described.

Origins and characterization of variants shared between databases of somatic and germline human mutations

Background

Mutations arise in the human genome in two major settings: the germline and the soma. These settings involve different inheritance patterns, time scales, chromatin structures, and environmental exposures, all of which impact the resulting distribution of substitutions. Nonetheless, many of the same single nucleotide variants (SNVs) are shared between germline and somatic mutation databases, such as between the gnomAD database of 120,000 germline exomes and the TCGA database of 10,000 somatic exomes. Here, we sought to explain this overlap.

Results

After strict filtering to exclude common germline polymorphisms and sites with poor coverage or mappability, we found 336,987 variants shared between the somatic and germline databases. A uniform statistical model explains 34% of these shared variants; a model that incorporates the varying mutation rates of the basic mutation types explains another 50% of shared variants; and a model that includes extended nucleotide contexts (e.g. surrounding 3 bases on either side) explains an additional 4% of shared variants. Analysis of read depth finds mixed evidence that up to 4% of the shared variants may represent germline variants leaked into somatic call sets. 9% of the shared variants are not explained by any model. Sequencing errors and convergent evolution did not account for these. We surveyed other factors as well: Cancers driven by endogenous mutational processes share a greater fraction of variants with the germline, and recently derived germline variants were more likely to be somatically shared than were ancient germline ones.

Conclusions

Overall, we find that shared variants largely represent bona fide biological occurrences of the same variant in the germline and somatic setting and arise primarily because DNA has some of the same basic chemical vulnerabilities in either setting. Moreover, we find mixed evidence that somatic call-sets leak appreciable numbers of germline variants, which is relevant to genomic privacy regulations. In future studies, the similar chemical vulnerability of DNA between the somatic and germline settings might be used to help identify disease-related genes by guiding the development of background-mutation models that are informed by both somatic and germline patterns of variation.

Characterizing Mutational Signatures in Human Cancer Cell Lines Reveals Episodic APOBEC Mutagenesis

Multiple signatures of somatic mutations have been identified in cancer genomes. Exome sequences of 1,001 human cancer cell lines and 577 xenografts revealed most common mutational signatures, indicating past activity of the underlying processes, usually in appropriate cancer types. To investigate ongoing patterns of mutational-signature generation, cell lines were cultured for extended periods and subsequently DNA sequenced. Signatures of discontinued exposures, including tobacco smoke and ultraviolet light, were not generated in vitro. Signatures of normal and defective DNA repair and replication continued to be generated at roughly stable mutation rates. Signatures of APOBEC cytidine deaminase DNA-editing exhibited substantial fluctuations in mutation rate over time with episodic bursts of mutations. The initiating factors for the bursts are unclear, although retrotransposon mobilization may contribute. The examined cell lines constitute a resource of live experimental models of mutational processes, which potentially retain patterns of activity and regulation operative in primary human cancers.

•

Annotation of mutational signatures across 1,001 cancer cell lines and 577 PDXs

•

Activities of mutational processes determined over time in cancer cell lines

•

APOBEC-associated mutagenesis is often ongoing and can be episodic

•

Detection of mutational signatures by single-cell sequencing

Molecular Profiles and Metastasis Markers in Chinese Patients with Gastric Carcinoma

The goal of this work was to investigate the molecular profiles and metastasis markers in Chinese patients with gastric carcinoma (GC). In total, we performed whole exome sequencing (WES) on 74 GC patients with tumor and adjacent normal formalin-fixed, paraffin-embedded (FFPE) tissue samples. The mutation spectrum of these samples showed a high concordance with TCGA and other studies on GC. PTPRT is significantly associated with metastasis of GC, suggesting its predictive role in metastasis of GC. Patients carrying BRCA2 mutations tend not to metastasize, which may be related to their sensitivity to chemotherapy. Mutations in MACF1, CDC27, HMCN1, CDH1 and PDZD2 were moderately enriched in peritoneal metastasis (PM) samples. Furthermore, we found two genomic regions (1p36.21 and Xq26.3) were associated with PM of GC, and patients with amplification of 1p36.21 and Xq26.3 have a worse prognosis (P = 0.002, 0.01, respectively). Our analysis provides GC patients with potential markers for single and combination therapies.

Phylogenetic Shifts in Gene Body Methylation Correlate with Gene Expression and Reflect Trait Conservation

A subset of genes in plant genomes are labeled with DNA methylation specifically at CG residues. These genes, known as gene-body methylated (gbM), have a number of associated characteristics. They tend to have longer sequences, to be enriched for intermediate expression levels, and to be associated with slower rates of molecular evolution. Most importantly, gbM genes tend to maintain their level of DNA methylation between species, suggesting that this trait is under evolutionary constraint. Given the degree of conservation in gbM, we still know surprisingly little about its function in plant genomes or whether gbM is itself a target of selection. To address these questions, we surveyed DNA methylation across eight grass (Poaceae) species that span a gradient of genome sizes. We first established that genome size correlates with genome-wide DNA methylation levels, but less so for genic levels. We then leveraged genomic data to identify a set of 2,982 putative orthologs among the eight species and examined shifts of methylation status for each ortholog in a phylogenetic context. A total of 55% of orthologs exhibited a shift in gbM, but these shifts occurred predominantly on terminal branches, indicating that shifts in gbM are rarely conveyed over time. Finally, we found that the degree of conservation of gbM across species is associated with increased gene length, reduced rates of molecular evolution, and increased gene expression level, but reduced gene expression variation across species. Overall, these observations suggest a basis for evolutionary pressure to maintain gbM status over evolutionary time.

Protein Interactions at Oxidized 5-Methylcytosine Bases

5-Methylcytosine (5mC), the major modified DNA base in mammalian cells, can be oxidized enzymatically to 5-hydroxymethylcytosine (5hmC), 5-formylcytosine (5fC), and 5-carboxylcytosine (5caC) by the Ten-Eleven-Translocation (TET) family of proteins. Whereas 5fC and 5caC are recognized and removed by base excision repair proteins, the 5hmC base accumulates to substantial levels in certain cell types such as brain-derived neurons and is viewed as a relatively stable DNA base. As such, the existence of "reader" proteins that recognize 5hmC would be a logical assumption, and various searches have been undertaken to identify proteins that specifically bind to 5hmC and the other oxidized 5mC bases. However, the existence of definitive 5hmC "readers" has remained unclear and proteins interacting specifically with 5fC or 5caC are also very few. On the other hand, 5hmC is incapable of interacting with a number of proteins that recognize 5mC at CpG sequences, suggesting that 5hmC is an anti-reader modification that may serve to displace 5mC readers from DNA. In this review article, we discuss candidate proteins that may interact with oxidized 5mC bases.

Hidden Markov models lead to higher resolution maps of mutation signature activity in cancer

Knowing the activity of the mutational processes shaping a cancer genome may provide insight into tumorigenesis and personalized therapy. It is thus important to characterize the signatures of active mutational processes in patients from their patterns of single base substitutions. However, mutational processes do not act uniformly on the genome, leading to statistical dependencies among neighboring mutations. To account for such dependencies, we develop the first sequence-dependent model, SigMa, for mutation signatures. We apply SigMa to characterize genomic and other factors that influence the activity of mutation signatures in breast cancer. We show that SigMa outperforms previous approaches, revealing novel insights on signature etiology. The source code for SigMa is publicly available at https://github.com/lrgr/sigma.

Canine osteosarcoma genome sequencing identifies recurrent mutations in DMD and the histone methyltransferase gene SETD2

Osteosarcoma (OS) is a rare, metastatic, human adolescent cancer that also occurs in pet dogs. To define the genomic underpinnings of canine OS, we performed multi-platform analysis of OS tumors from 59 dogs, including whole genome sequencing (n = 24) and whole exome sequencing (WES; n = 13) of primary tumors and matched normal tissue, WES (n = 10) of matched primary/metastatic/normal samples and RNA sequencing (n = 54) of primary tumors. We found that canine OS recapitulates features of human OS including low point mutation burden (median 1.98 per Mb) with a trend towards higher burden in metastases, high structural complexity, frequent TP53 (71%), PI3K pathway (37%), and MAPK pathway mutations (17%), and low expression of immune-associated genes. We also identified novel features of canine OS including putatively inactivating somatic SETD2 (42%) and DMD (50%) aberrations. These findings set the stage for understanding OS development in dogs and humans, and establish genomic contexts for future comparative analyses.

Germline Predisposition and Copy Number Alteration in Pre-stage Lung Adenocarcinomas Presenting as Ground-Glass Nodules

Objective: Synchronous multiple ground-glass nodules (SM-GGNs) are a distinct entity of lung cancer which has been emerging increasingly in recent years in China. The oncogenesis molecular mechanisms of SM-GGNs remain elusive. Methods: We investigated single nucleotide variations (SNV), insertions and deletions (INDEL), somatic copy number variations (CNV), and germline mutations of 69 SM-GGN samples collected from 31 patients, using target sequencing (TRS) and whole exome sequencing (WES). Results: In the entire cohort, many known driver mutations were found, including EGFR (21.7%), BRAF (14.5%), and KRAS (6%). However, only one out of the 31 patients had the same somatic missense or truncated events within SM-GGNs, indicating the independent origins for almost all of these SM-GGNs. Many germline mutations with a low frequency in the Chinese population, and genes harboring both germline and somatic variations, were discovered in these pre-stage GGNs. These GGNs also bore large segments of copy number gains and/or losses. The CNV segment number tended to be positively correlated with the germline mutations (r = 0.57). The CNV sizes were correlated with the somatic mutations (r = 0.55). A moderate correlation (r = 0.54) was also shown between the somatic and germline mutations. Conclusion: Our data suggests that the precancerous unstable CNVs with potentially predisposing genetic backgrounds may foster the onset of driver mutations and the development of independent SM-GGNs during the local stimulation of mutagens.

DNA Methylation Patterns in the Social Spider, Stegodyphus dumicola

Variation in DNA methylation patterns among genes, individuals, and populations appears to be highly variable among taxa, but our understanding of the functional significance of this variation is still incomplete. We here present the first whole genome bisulfite sequencing of a chelicerate species, the social spider Stegodyphus dumicola. We show that DNA methylation occurs mainly in CpG context and is concentrated in genes. This is a pattern also documented in other invertebrates. We present RNA sequence data to investigate the role of DNA methylation in gene regulation and show that, within individuals, methylated genes are more expressed than genes that are not methylated and that methylated genes are more stably expressed across individuals than unmethylated genes. Although no causal association is shown, this lends support for the implication of DNA CpG methylation in regulating gene expression in invertebrates. Differential DNA methylation between populations showed a small but significant correlation with differential gene expression. This is consistent with a possible role of DNA methylation in local adaptation. Based on indirect inference of the presence and pattern of DNA methylation in chelicerate species whose genomes have been sequenced, we performed a comparative phylogenetic analysis. We found strong evidence for exon DNA methylation in the horseshoe crab Limulus polyphemus and in all spider and scorpion species, while most Parasitiformes and Acariformes species seem to have lost DNA methylation.

Hypermutagenesis in untreated adult gliomas due to inherited mismatch mutations

Hypermutagenesis refers to marked increase in the number of mutations due to continuous mutagenic process. Hypermutated tumors, have being found in several tumor types, are associated with inherited or acquired alterations in the DNA repair pathways. Hypermutation has been observed in a subset of adult glioma patients as a direct result of temozolomide(TMZ)-induced mutagenesis. In our study, we have identified a rare subset of treatment-naïve adult gliomas with de novo hypermutator phenotype and explored the evolution of spontaneous and treatment-induced hypermutagenesis. We conducted Whole-Exome Sequencing (WES), Whole-Transcriptome Sequencing (WTS), and Single-Cell Sequencing (SCS) of TMZ-naïve and post-TMZ-treated hypermutated tumors to identify distinct clinical or genomic manifestations that contribute to the development of hypermutation in untreated adult gliomas. TMZ-naïve hypermutated tumors were marked by absence of IDH1 somatic mutation and MGMT promoter (pMGMT) methylation, two genomic traits that were significantly associated with the TMZ-induced hypermutagenic event in glioblastoma, and harbored inherited alterations in the mismatch repair (MMR) machinery. The immediate family members of the TMZ-naive hypermutated glioma patients were also previous diagnosed with cancer development history, suggesting that germline dysfunction of the MMR pathway could potentially pose hereditary risk to genetic predisposition of carcinogenesis in gliomas. Lastly, both TMZ-naïve and post-TMZ-treated hypermutated tumors exhibited a significant accumulation of neoantigen loads, suggesting immunotherapeutic alternatives. Our results present new and unique understanding of hypermutagenic process in adult gliomas and an important step towards clinical implication of immunotherapy in glioma treatment.

High prevalence of focal and multi-focal somatic genetic variants in the human brain

Somatic mutations during stem cell division are responsible for several cancers. In principle, a similar process could occur during the intense cell proliferation accompanying human brain development, leading to the accumulation of regionally distributed foci of mutations. Using dual platform >5000-fold depth sequencing of 102 genes in 173 adult human brain samples, we detect and validate somatic mutations in 27 of 54 brains. Using a mathematical model of neurodevelopment and approximate Bayesian inference, we predict that macroscopic islands of pathologically mutated neurons are likely to be common in the general population. The detected mutation spectrum also includes DNMT3A and TET2 which are likely to have originated from blood cell lineages. Together, these findings establish developmental mutagenesis as a potential mechanism for neurodegenerative disorders, and provide a novel mechanism for the regional onset and focal pathology in sporadic cases.

Spectrum of MECP2 mutations in Vietnamese patients with RETT syndrome

Background

Rett syndrome (RTT) is a severe neurodevelopmental disorder in children characterized by a normal neurodevelopmental process in the first 6–18 months followed by a period of motor and vocal deterioration with stereotypic hand movements. Incidence of RTT is mostly due to de novo mutation in the MECP2 gene (methyl-CpG-binding protein 2).

Methods

The study assessed 27 female patients presented with classical RTT phenotype age range from 18 months to 48 months. Specialist carried out the clinical evaluation and diagnosis according to RTT diagnosis criteria. Blood samples from patients were then collected for genomic DNA extraction. We next performed MECP2 gene amplification and sequencing of the whole coding region to screen for mutations.

Result

MECP2 mutation was found in 20 patients (74%) including: 2 missense, 4 nonsense, 6 frameshift and 2 deletion mutation. The study identified 14 pathogenic mutations which we found 4 mutation, to our knowledge and extensive search, not priory reported in any mutation database or publication: c.1384-1385DelGT, c.1205insT, c.717delC and c.1132_1207del77. High percentage of C > T (70%) in CpG sites mutation was found.

Conclusion

Our result reveals a high percentage of C > T mutation in CpG hot spot, which is more prone to modification and more likely to be detected in RTT as a disorder is strictly due to de novo mutations. The study is the first to identify the mutation spectrum of MECP2 gene in Vietnamese patients and also an important step toward better diagnosis and care for RTT patients in Vietnam.

Excision of the doubly methylated base N4,5-dimethylcytosine from DNA by Escherichia coli Nei and Fpg proteins

Cytosine (C) in DNA is often modified to 5-methylcytosine (m5C) to execute important cellular functions. Despite the significance of m5C for epigenetic regulation in mammals, damage to m5C has received little attention. For instance, almost no studies exist on erroneous methylation of m5C by alkylating agents to doubly or triply methylated bases. Owing to chemical evidence, and because many prokaryotes express methyltransferases able to convert m5C into N4,5-dimethylcytosine (m N4,5C) in DNA, m N4,5C is probably present in vivo We screened a series of glycosylases from prokaryotic to human and found significant DNA incision activity of the Escherichia coli Nei and Fpg proteins at m N4,5C residues in vitro The activity of Nei was highest opposite cognate guanine followed by adenine, thymine (T) and C. Fpg-complemented Nei by exhibiting the highest activity opposite C followed by lower activity opposite T. To our knowledge, this is the first description of a repair enzyme activity at a further methylated m5C in DNA, as well as the first alkylated base allocated as a Nei or Fpg substrate. Based on our observed high sensitivity to nuclease S1 digestion, we suggest that m N4,5C occurs as a disturbing lesion in DNA and that Nei may serve as a major DNA glycosylase in E. coli to initiate its repair.This article is part of a discussion meeting issue 'Frontiers in epigenetic chemical biology'.

Defining Driver DNA Methylation Changes in Human Cancer

Human malignant tumors are characterized by pervasive changes in the patterns of DNA methylation. These changes include a globally hypomethylated tumor cell genome and the focal hypermethylation of numerous 5′-cytosine-phosphate-guanine-3′ (CpG) islands, many of them associated with gene promoters. It has been challenging to link specific DNA methylation changes with tumorigenesis in a cause-and-effect relationship. Some evidence suggests that cancer-associated DNA hypomethylation may increase genomic instability. Promoter hypermethylation events can lead to silencing of genes functioning in pathways reflecting hallmarks of cancer, including DNA repair, cell cycle regulation, promotion of apoptosis or control of key tumor-relevant signaling networks. A convincing argument for a tumor-driving role of DNA methylation can be made when the same genes are also frequently mutated in cancer. Many of the most commonly hypermethylated genes encode developmental transcription factors, the methylation of which may lead to permanent gene silencing. Inactivation of such genes will deprive the cells in which the tumor may initiate from the option of undergoing or maintaining lineage differentiation and will lock them into a perpetuated stem cell-like state thus providing an additional window for cell transformation.

H3K36me3-mediated mismatch repair preferentially protects actively transcribed genes from mutation

Histone H3 trimethylation at lysine 36 (H3K36me3) is an important histone mark involved in both transcription elongation and DNA mismatch repair (MMR). It is known that H3K36me3 recruits the mismatch-recognition protein MutSα to replicating chromatin via its physical interaction with MutSα's PWWP domain, but the exact role of H3K36me3 in transcription is undefined. Using ChIP combined with whole-genome DNA sequencing analysis, we demonstrate here that H3K36me3, together with MutSα, is involved in protecting against mutation, preferentially in actively transcribed genomic regions. We found that H3K36me3 and MutSα are much more co-enriched in exons and actively transcribed regions than in introns and nontranscribed regions. The H3K36me3-MutSα co-enrichment correlated with a much lower mutation frequency in exons and actively transcribed regions than in introns and nontranscribed regions. Correspondingly, depleting H3K36me3 or disrupting the H3K36me3-MutSα interaction elevated the spontaneous mutation frequency in actively transcribed genes, but it had little influence on the mutation frequency in nontranscribed or transcriptionally inactive regions. Similarly, H₂O₂-induced mutations, which mainly cause base oxidations, preferentially occurred in actively transcribed genes in MMR-deficient cells. The data presented here suggest that H3K36me3-mediated MMR preferentially safeguards actively transcribed genes not only during replication by efficiently correcting mispairs in early replicating chromatin but also during transcription by directly or indirectly removing DNA lesions associated with a persistently open chromatin structure.

Alu exaptation enriches the human transcriptome by introducing new gene ends

In mammals, transposable elements are largely silenced, but under fortuitous circumstances may be co-opted to play a functional role. Here, we show that when Alu elements are inserted within or nearby genes in sense orientation, they may contribute to the transcriptome diversity by forming new cleavage and polyadenylation sites. We mapped these new gene ends in human onto the Alu sequence and identified three hotspots of cleavage and polyadenylation site formation. Interestingly, the native Alu sequence does not contain any canonical polyadenylation signal. We therefore studied what evolutionary processes might explain the formation of these specific hotspots of novel gene ends. We show that two of the three hotspots might have emerged from mutational processes that turned sequences that resemble polyadenylation signals into full-blown canonical signals, whereas one hotspot is tightly linked to the process of Alu insertion into the genome. Overall, Alu elements may lie behind the formation of 302 new gene end variants, affecting a total of 243 genes. Intergenic Alu elements may elongate genes by creating a downstream cleavage site, intronic Alu elements may lead to gene variants which code for truncated proteins, and 3'UTR Alu elements may result in gene variants with alternative 3'UTR.

Genomic and transcriptomic characterization of skull base chordoma

Skull base chordoma is a primary rare malignant bone-origin tumor showing relatively slow growth pattern and locally destructive lesions, which can only be characterized by histologic components. There is no available prognostic or therapeutic biomarker to predict clinical outcome or treatment response and the molecular mechanisms underlying chordoma development still remain unexplored. Therefore, we sought out to identify novel somatic variations that are associated with chordoma progression and potentially employed as therapeutic targets. Thirteen skull base chordomas were subjected for whole-exome and/or whole-transcriptome sequencing. In process, we have identified chromosomal aberration in 1p, 7, 10, 13 and 17q, high frequency of functional germline SNP of the T gene, rs2305089 (P = 0.0038) and several recurrent alterations including MUC4, NBPF1, NPIPB15 mutations and novel gene fusion of SAMD5-SASH1 for the first time in skull base chordoma.

Mutational signatures and chromosome alteration profiles of squamous cell carcinomas of the vulva

Vulvar squamous cell carcinoma (SCC) consists of two different etiologic categories: human papilloma virus (HPV)-associated (HPV (+)) and HPV-non-associated (HPV (−)). There have been no genome-wide studies on the genetic alterations of vulvar SCCs or on the differences between HPV (+) and HPV (−) vulvar SCCs. In this study, we performed whole-exome sequencing and copy number profiling of 6 HPV (+) and 9 HPV (−) vulvar SCCs and found known mutations (TP53, CDKN2A and HRAS) and copy number alterations (CNAs) (7p and 8q gains and 2q loss) in HPV (−) SCCs. In HPV (+), we found novel mutations in PIK3CA, BRCA2 and FBXW7 that had not been reported in vulvar SCCs. HPV (−) SCCs exhibited more mutational loads (numbers of nonsilent mutations and driver mutations) than HPV (+) SCCs, but the CNA loads and mutation signatures between HPV (+) and HPV (−) SCCs did not differ. Of note, 40% and 40% of the 15 vulvar SCCs harbored PIK3CA and FAT1 alterations, respectively. In addition, we found that the SCCs harbored kataegis (a localized hypermutation) in 2 HPV (+) SCCs and copy-neutral losses of heterozygosity in 4 (one HPV (+) and 3 HPV (−)) SCCs. Our data indicate that HPV (+) and HPV (−) vulvar SCCs may have different mutation and CNA profiles but that there are genomic features common to SCCs. Our data provide useful information for both HPV (+) and HPV (−) vulvar SCCs and may aid in the development of clinical treatment strategies.

High glucose induces podocyte epithelial‑to‑mesenchymal transition by demethylation‑mediated enhancement of MMP9 expression

Abnormal expression of matrix metalloproteinase 9 (MMP9) is correlated with podocyte epithelial-to-mesenchymal transition (EMT) in diabetic nephropathy (DN). However, the mechanisms underlying this process are not well defined. Site-specific demethylation may sustain high expression levels of target genes. In the present study, in order to investigate the association between DNA demethylation of MMP9 promoter and podocyte EMT in DN, human podocytes were cultured in high-glucose (HG) medium and a rat model of DN was established by intraperitoneal injection of streptozotocin (STZ) to determine whether site-specific demethylation of the MMP9 promoter was involved in regulating podocyte EMT in DN. The MTT assay was used to assess the effects of HG culture on the growth of podocytes, and the demethylation status of the MMP9 promoter was assessed by bisulfite sequencing polymerase chain reaction. mRNA and protein expression levels of MMP9, α-smooth muscle actin (α-SMA), podocalyxin and fibronectin-1 in podocytes were assessed by reverse transcription-quantitative PCR (RT-qPCR) and western blot analyses. The results demonstrated that HG treatment up regulated the expression of MMP9, α-SMA and fibronectin-1, but down regulated the expression of podocalyxin in podocytes. The MMP9 promoter region was revealed to contain a variety of demethylated CpG sites, and HG treatment reduced the rate of MMP9 promotermethylation, which, in turn, enhanced its promoter activity. In summary, these data suggested that demethylation of the MMP9 promoter may serve an important role in podocyte EMT in DN. The demethylation status of the MMP9 promoter maybe used as an important prognostic marker of DN in clinic.

BeWith: A Between-Within method to discover relationships between cancer modules via integrated analysis of mutual exclusivity, co-occurrence and functional interactions

The analysis of the mutational landscape of cancer, including mutual exclusivity and co-occurrence of mutations, has been instrumental in studying the disease. We hypothesized that exploring the interplay between co-occurrence, mutual exclusivity, and functional interactions between genes will further improve our understanding of the disease and help to uncover new relations between cancer driving genes and pathways. To this end, we designed a general framework, BeWith, for identifying modules with different combinations of mutation and interaction patterns. We focused on three different settings of the BeWith schema: (i) BeME-WithFun, in which the relations between modules are enriched with mutual exclusivity, while genes within each module are functionally related; (ii) BeME-WithCo, which combines mutual exclusivity between modules with co-occurrence within modules; and (iii) BeCo-WithMEFun, which ensures co-occurrence between modules, while the within module relations combine mutual exclusivity and functional interactions. We formulated the BeWith framework using Integer Linear Programming (ILP), enabling us to find optimally scoring sets of modules. Our results demonstrate the utility of BeWith in providing novel information about mutational patterns, driver genes, and pathways. In particular, BeME-WithFun helped identify functionally coherent modules that might be relevant for cancer progression. In addition to finding previously well-known drivers, the identified modules pointed to other novel findings such as the interaction between NCOR2 and NCOA3 in breast cancer. Additionally, an application of the BeME-WithCo setting revealed that gene groups differ with respect to their vulnerability to different mutagenic processes, and helped us to uncover pairs of genes with potentially synergistic effects, including a potential synergy between mutations in TP53 and the metastasis related DCC gene. Overall, BeWith not only helped us uncover relations between potential driver genes and pathways, but also provided additional insights on patterns of the mutational landscape, going beyond cancer driving mutations. Implementation is available at https://www.ncbi.nlm.nih.gov/CBBresearch/Przytycka/software/bewith.html

The genomic landscape of cancer obtained from large scale studies revealed mutational patterns such as mutual exclusivity and co-occurrences. Despite significant efforts, the understanding of the patterns harbored by specific genes and their interplay with functional interactions are still limited. Both mutual exclusivity and co-occurrence can arise due to several reasons. For example, two functionally interacting genes may dysregulate the same cancer related pathway when either of them is mutated, leading to mutually exclusive mutations. Alternatively, mutual exclusivity might reflect mutations specific to two different cancer types. Methods for joint analysis of co-occurrence, mutual exclusivity, and functional interaction relationships can lead to a better understanding of the causes and impacts of mutations in cancer. We report a new computational approach, BeWith, which identifies groups of genes (or gene modules) with coherent patterns within modules, but distinct properties among genes in different modules. The general formulation of our method allows us to investigate various aspects of the cancer mutational landscape, leading to uncovering relationships between mutated gene modules, cancer subtypes, and mutational signatures.

The interaction between cytosine methylation and processes of DNA replication and repair shape the mutational landscape of cancer genomes

Methylated cytosines (5mCs) are frequently mutated in the genome. However, no studies have yet comprehensively analysed mutation–methylation associations across cancer types. Here we analyse 916 cancer genomes, together with tissue type-specific methylation and replication timing data. We describe a strong mutation–methylation association across colorectal cancer subtypes, most interestingly in samples with microsatellite instability (MSI) or Polymerase epsilon (POLE) exonuclease domain mutations. By analysing genomic regions with differential mismatch repair (MMR) efficiency, we suggest a possible role for MMR in the correction of 5mC deamination events, potentially accounting for the high rate of 5mC mutation accumulation in MSI tumours. Additionally, we propose that mutant POLE asserts a mutator phenotype specifically at 5mCs, and we find coding mutation hotspots in POLE-mutant cancers at highly-methylated CpGs in the tumour-suppressor genes APC and TP53. Finally, using multivariable regression models, we demonstrate that different cancers exhibit distinct mutation–methylation associations, with DNA repair influencing such associations in certain cancer genomes. Taken together, we find differential associations with methylation that are vital for accurately predicting expected mutation loads across cancer types. Our findings reveal links between methylation and common mutation and repair processes, with these mechanisms defining a key part of the mutational landscape of cancer genomes.

Construction of a combinatorial pipeline using two somatic variant calling methods for whole exome sequence data of gastric cancer

High-throughput next-generation sequencing is a powerful tool to identify the genotypic landscapes of somatic variants and therapeutic targets in various cancers including gastric cancer, forming the basis for personalized medicine in the clinical setting. Although the advent of many computational algorithms leads to higher accuracy in somatic variant calling, no standard method exists due to the limitations of each method. Here, we constructed a new pipeline. We combined two different somatic variant callers with different algorithms, Strelka and VarScan 2, and evaluated performance using whole exome sequencing data obtained from 19 Japanese cases with gastric cancer (GC); then, we characterized these tumors based on identified driver molecular alterations. More single nucleotide variants (SNVs) and small insertions/deletions were detected by Strelka and VarScan 2, respectively. SNVs detected by both tools showed higher accuracy for estimating somatic variants compared with those detected by only one of the two tools and accurately showed the mutation signature and mutations of driver genes reported for GC. Our combinatorial pipeline may have an advantage in detection of somatic mutations in GC and may be useful for further genomic characterization of Japanese patients with GC to improve the efficacy of GC treatments. J. Med. Invest. 64: 233-240, August, 2017.

Equivalent missense variant in the FOXP2 and FOXP1 transcription factors causes distinct neurodevelopmental disorders

The closely related paralogues FOXP2 and FOXP1 encode transcription factors with shared functions in the development of many tissues, including the brain. However, while mutations in FOXP2 lead to a speech/language disorder characterized by childhood apraxia of speech (CAS), the clinical profile of FOXP1 variants includes a broader neurodevelopmental phenotype with global developmental delay, intellectual disability, and speech/language impairment. Using clinical whole-exome sequencing, we report an identical de novo missense FOXP1 variant identified in three unrelated patients. The variant, p.R514H, is located in the forkhead-box DNA-binding domain and is equivalent to the well-studied p.R553H FOXP2 variant that cosegregates with CAS in a large UK family. We present here for the first time a direct comparison of the molecular and clinical consequences of the same mutation affecting the equivalent residue in FOXP1 and FOXP2. Detailed functional characterization of the two variants in cell model systems revealed very similar molecular consequences, including aberrant subcellular localization, disruption of transcription factor activity, and deleterious effects on protein interactions. Nonetheless, clinical manifestations were broader and more severe in the three cases carrying the p.R514H FOXP1 variant than in individuals with the p.R553H variant related to CAS, highlighting divergent roles of FOXP2 and FOXP1 in neurodevelopment.

Postzygotic single-nucleotide mosaicisms contribute to the etiology of autism spectrum disorder and autistic traits and the origin of mutations

The roles and characteristics of postzygotic single‐nucleotide mosaicisms (pSNMs) in autism spectrum disorders (ASDs) remain unclear. In this study of the whole exomes of 2,361 families in the Simons Simplex Collection, we identified 1,248 putative pSNMs in children and 285 de novo SNPs in children with detectable parental mosaicism. Ultra‐deep amplicon resequencing suggested a validation rate of 51%. Analyses of validated pSNMs revealed that missense/loss‐of‐function (LoF) pSNMs with a high mutant allele fraction (MAF≥ 0.2) contributed to ASD diagnoses (P = 0.022, odds ratio [OR] = 5.25), whereas missense/LoF pSNMs with a low MAF (MAF<0.2) contributed to autistic traits in male non‐ASD siblings (P = 0.033). LoF pSNMs in parents were less likely to be transmitted to offspring than neutral pSNMs (P = 0.037), and missense/LoF pSNMs in parents with a low MAF were transmitted more to probands than to siblings (P = 0.016, OR = 1.45). We estimated that pSNMs in probands or de novo mutations inherited from parental pSNMs increased the risk of ASD by approximately 6%. Adding pSNMs into the transmission and de novo association test model revealed 13 new ASD risk genes. These results expand the existing repertoire of genes involved in ASD and shed new light on the contribution of genomic mosaicisms to ASD diagnoses and autistic traits.

An evolutionary case for functional gene body methylation in plants and animals

Methylation in the bodies of active genes is common in animals and vascular plants. Evolutionary patterns indicate homeostatic functions for this type of methylation.

Stress and the Emerging Roles of Chromatin Remodeling in Signal Integration and Stable Transmission of Reversible Phenotypes

The influence of early life experience and degree of parental-infant attachment on emotional development in children and adolescents has been comprehensively studied. Structural and mechanistic insight into the biological foundation and maintenance of mammalian defensive systems (metabolic, immune, nervous and behavioral) is slowly advancing through the emerging field of developmental molecular (epi)genetics. Initial evidence revealed that differential nurture early in life generates stable differences in offspring hypothalamic-pituitary-adrenal (HPA) regulation, in part, through chromatin remodeling and changes in DNA methylation of specific genes expressed in the brain, revealing physical, biochemical and molecular paths for the epidemiological concept of gene-environment interactions. Herein, a primary molecular mechanism underpinning the early developmental programming and lifelong maintenance of defensive (emotional) responses in the offspring is the alteration of chromatin domains of specific genomic regions from a condensed state (heterochromatin) to a transcriptionally accessible state (euchromatin). Conversely, DNA methylation promotes the formation of heterochromatin, which is essential for gene silencing, genomic integrity and chromosome segregation. Therefore, inter-individual differences in chromatin modifications and DNA methylation marks hold great potential for assessing the impact of both early life experience and effectiveness of intervention programs—from guided psychosocial strategies focused on changing behavior to pharmacological treatments that target chromatin remodeling and DNA methylation enzymes to dietary approaches that alter cellular pools of metabolic intermediates and methyl donors to affect nutrient bioavailability and metabolism. In this review article, we discuss the potential molecular mechanism(s) of gene regulation associated with chromatin modeling and programming of endocrine (e.g., HPA and metabolic or cardiovascular) and behavioral (e.g., fearfulness, vigilance) responses to stress, including alterations in DNA methylation and the role of DNA repair machinery. From parental history (e.g., drugs, housing, illness, nutrition, socialization) to maternal-offspring exchanges of nutrition, microbiota, antibodies and stimulation, the nature of nurture provides not only mechanistic insight into how experiences propagate from external to internal variables, but also identifies a composite therapeutic target, chromatin modeling, for gestational/prenatal stress, adolescent anxiety/depression and adult-onset neuropsychiatric disease.