Pancreatic cancer acquires resistance to MAPK pathway inhibition by clonal expansion and adaptive DNA hypermethylation.

BACKGROUND: Pancreatic ductal adenocarcinoma (PDAC) is an aggressive cancer with poor prognosis. It is marked by extraordinary resistance to conventional therapies including chemotherapy and radiation, as well as to essentially all targeted therapies evaluated so far. More than 90% of PDAC cases harbor an activating KRAS mutation. As the most common KRAS variants in PDAC remain undruggable so far, it seemed promising to inhibit a downstream target in the MAPK pathway such as MEK1/2, but up to now preclinical and clinical evaluation of MEK inhibitors (MEKi) failed due to inherent and acquired resistance mechanisms. To gain insights into molecular changes during the formation of resistance to oncogenic MAPK pathway inhibition, we utilized short-term passaged primary tumor cells from ten PDACs of genetically engineered mice. We followed gain and loss of resistance upon MEKi exposure and withdrawal by longitudinal integrative analysis of whole genome sequencing, whole genome bisulfite sequencing, RNA-sequencing and mass spectrometry data. RESULTS: We found that resistant cell populations under increasing MEKi treatment evolved by the expansion of a single clone but were not a direct consequence of known resistance-conferring mutations. Rather, resistant cells showed adaptive DNA hypermethylation of 209 and hypomethylation of 8 genomic sites, most of which overlap with regulatory elements known to be active in murine PDAC cells. Both DNA methylation changes and MEKi resistance were transient and reversible upon drug withdrawal. Furthermore, MEKi resistance could be reversed by DNA methyltransferase inhibition with remarkable sensitivity exclusively in the resistant cells. CONCLUSION: Overall, the concept of acquired therapy resistance as a result of the expansion of a single cell clone with epigenetic plasticity sheds light on genetic, epigenetic and phenotypic patterns during evolvement of treatment resistance in a tumor with high adaptive capabilities and provides potential for reversion through epigenetic targeting.

Early-set POMC methylation variability is accompanied by increased risk for obesity and is addressable by MC4R agonist treatment.

Increasing evidence points toward epigenetic variants as a risk factor for developing obesity. We analyzed DNA methylation of the POMC (pro-opiomelanocortin) gene, which is pivotal for satiety regulation. We identified sex-specific and nongenetically determined POMC hypermethylation associated with a 1.4-fold (confidence interval, 1.03 to 2.04) increased individual risk of developing obesity. To investigate the early embryonic establishment of POMC methylation states, we established a human embryonic stem cell (hESC) model. Here, hESCs (WA01) were transferred into a naïve state, which was associated with a reduction of DNA methylation. Naïve hESCs were differentiated via a formative state into POMC-expressing hypothalamic neurons, which was accompanied by re-establishment of DNA methylation patterning. We observed that reduced POMC gene expression was associated with increased POMC methylation in POMC-expressing neurons. On the basis of these findings, we treated POMC-hypermethylated obese individuals (n = 5) with an MC4R agonist and observed a body weight reduction of 4.66 ± 2.16% (means ± SD) over a mean treatment duration of 38.4 ± 26.0 weeks. In summary, we identified an epigenetic obesity risk variant at the POMC gene fulfilling the criteria for a metastable epiallele established in early embryonic development that may be addressable by MC4R agonist treatment to reduce body weight.

Prevalence of Pathogenic Germline Variants in Women with Non-Familial Unilateral Triple-Negative Breast Cancer.

Introduction: International guidelines recommend genetic testing for women with familial breast cancer at an expected prevalence of pathogenic germline variants (PVs) of at least 10%. In a study sample of the German Consortium for Hereditary Breast and Ovarian Cancer (GC-HBOC), we have previously shown that women with TNBC diagnosed before the age of 50 years but without a family history of breast or ovarian cancer (sTNBC) meet this criterion. The present study investigates the PV prevalence in BRCA1, BRCA2, and nine additional cancer predisposition genes in an extended sTNBC study sample including a cohort of women with a later age at sTNBC diagnosis. Patients and Methods: In 1,600 women with sTNBC (median age at diagnosis: 41 years, range 19-78 years), we investigated the association between age at diagnosis and PV occurrence in cancer predisposition genes using logistic regression. Results: 260 sTNBC patients (16.2%) were found to have a PV in cancer predisposition genes (BRCA1: n = 170 [10.6%]; BRCA2: n = 46 [2.9%], other: n = 44 [2.8%]). The PV prevalence in women diagnosed between 50 and 59 years (n = 194) was 11.3% (22/194). Logistic regression showed a significant increase in PV prevalence with decreasing age at diagnosis (OR 1.41 per 10 years younger age at diagnosis; 95% confidence interval: 1.21-1.65; p < 0.001). The PV prevalence predicted by the model was above 10% for diagnoses before the age of 56.8 years. Conclusion: Based on the data presented, we recommend genetic testing by gene panel analysis for sTNBC patients diagnosed before the age of 60 years. Due to the still uncertain estimate for women with sTNBC diagnosed above the age of 60 years, further studies are needed.

Bi-allelic variants of FILIP1 cause congenital myopathy, dysmorphism and neurological defects.

Filamin-A-interacting protein 1 (FILIP1) is a structural protein that is involved in neuronal and muscle function and integrity and interacts with FLNa and FLNc. Pathogenic variants in filamin-encoding genes have been linked to neurological disorders (FLNA) and muscle diseases characterized by myofibrillar perturbations (FLNC), but human diseases associated with FILIP1 variants have not yet been described. Here, we report on five patients from four unrelated consanguineous families with homozygous FILIP1 variants (two nonsense and two missense). Functional studies indicated altered stability of the FILIP1 protein carrying the p.[Pro1133Leu] variant. Patients exhibit a broad spectrum of neurological symptoms including brain malformations, neurodevelopmental delay, muscle weakness and pathology and dysmorphic features. Electron and immunofluorescence microscopy on the muscle biopsy derived from the patient harbouring the homozygous p.[Pro1133Leu] missense variant revealed core-like zones of myofibrillar disintegration, autophagic vacuoles and accumulation of FLNc. Proteomic studies on the fibroblasts derived from the same patient showed dysregulation of a variety of proteins including FLNc and alpha-B-crystallin, a finding (confirmed by immunofluorescence) which is in line with the manifestation of symptoms associated with the syndromic phenotype of FILIP1opathy. The combined findings of this study show that the loss of functional FILIP1 leads to a recessive disorder characterized by neurological and muscular manifestations as well as dysmorphic features accompanied by perturbed proteostasis and myopathology.

Insane in the vembrane: filtering and transforming VCF/BCF files.

SUMMARY: We present vembrane as a command line variant call format (VCF)/binary call format (BCF) filtering tool that consolidates and extends the filtering functionality of previous software to meet any imaginable filtering use case. Vembrane exposes the VCF/BCF file type specification and its inofficial extensions by the annotation tools VEP and SnpEff as Python data structures. vembrane filter enables filtration by Python expressions, requiring only basic knowledge of the Python programming language. vembrane table allows users to generate tables from subsets of annotations or functions thereof. Finally, it is fast, by using pysam and relying on lazy evaluation. AVAILABILITY AND IMPLEMENTATION: Source code and installation instructions are available at github.com/vembrane/vembrane (doi: 10.5281/zenodo.7003981). SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Systematic analysis and prediction of genes associated with monogenic disorders on human chromosome X.

Disease gene discovery on chromosome (chr) X is challenging owing to its unique modes of inheritance. We undertook a systematic analysis of human chrX genes. We observe a higher proportion of disorder-associated genes and an enrichment of genes involved in cognition, language, and seizures on chrX compared to autosomes. We analyze gene constraints, exon and promoter conservation, expression, and paralogues, and report 127 genes sharing one or more attributes with known chrX disorder genes. Using machine learning classifiers trained to distinguish disease-associated from dispensable genes, we classify 247 genes, including 115 of the 127, as having high probability of being disease-associated. We provide evidence of an excess of variants in predicted genes in existing databases. Finally, we report damaging variants in CDK16 and TRPC5 in patients with intellectual disability or autism spectrum disorders. This study predicts large-scale gene-disease associations that could be used for prioritization of X-linked pathogenic variants.

No evidence for intervention-associated DNA methylation changes in monocytes of patients with posttraumatic stress disorder.

DNA methylation patterns can be responsive to environmental influences. This observation has sparked interest in the potential for psychological interventions to influence epigenetic processes. Recent studies have observed correlations between DNA methylation changes and therapy outcome. However, most did not control for changes in cell composition. This study had two aims: first, we sought to replicate therapy-associated changes in DNA methylation of commonly assessed candidate genes in isolated monocytes from 60 female patients with post-traumatic stress disorder (PTSD). Our second, exploratory goal was to identify novel genomic regions with substantial pre-to-post intervention DNA methylation changes by performing whole-genome bisulfite sequencing (WGBS) in two patients with PTSD. Equivalence testing and Bayesian analyses provided evidence against physiologically meaningful intervention-associated DNA methylation changes in monocytes of PTSD patients in commonly investigated target genes (NR3C1, FKBP5, SLC6A4, OXTR). Furthermore, WGBS yielded only a limited set of candidate regions with suggestive evidence of differential DNA methylation pre- to post-therapy. These differential DNA methylation patterns did not prove replicable when investigated in the entire cohort. We conclude that there is no evidence for major, recurrent intervention-associated DNA methylation changes in the investigated genes in monocytes of patients with PTSD.

Lasp1 regulates adherens junction dynamics and fibroblast transformation in destructive arthritis.

The LIM and SH3 domain protein 1 (Lasp1) was originally cloned from metastatic breast cancer and characterised as an adaptor molecule associated with tumourigenesis and cancer cell invasion. However, the regulation of Lasp1 and its function in the aggressive transformation of cells is unclear. Here we use integrative epigenomic profiling of invasive fibroblast-like synoviocytes (FLS) from patients with rheumatoid arthritis (RA) and from mouse models of the disease, to identify Lasp1 as an epigenomically co-modified region in chronic inflammatory arthritis and a functionally important binding partner of the Cadherin-11/ß-Catenin complex in zipper-like cell-to-cell contacts. In vitro, loss or blocking of Lasp1 alters pathological tissue formation, migratory behaviour and platelet-derived growth factor response of arthritic FLS. In arthritic human TNF transgenic mice, deletion of Lasp1 reduces arthritic joint destruction. Therefore, we show a function of Lasp1 in cellular junction formation and inflammatory tissue remodelling and identify Lasp1 as a potential target for treating inflammatory joint disorders associated with aggressive cellular transformation.

GC-rich repeat expansions: associated disorders and mechanisms.

Noncoding repeat expansions are a well-known cause of genetic disorders mainly affecting the central nervous system. Missed by most standard technologies used in routine diagnosis, pathogenic noncoding repeat expansions have to be searched for using specific techniques such as repeat-primed PCR or specific bioinformatics tools applied to genome data, such as ExpansionHunter. In this review, we focus on GC-rich repeat expansions, which represent at least one third of all noncoding repeat expansions described so far. GC-rich expansions are mainly located in regulatory regions (promoter, 5' untranslated region, first intron) of genes and can lead to either a toxic gain-of-function mediated by RNA toxicity and/or repeat-associated non-AUG (RAN) translation, or a loss-of-function of the associated gene, depending on their size and their methylation status. We herein review the clinical and molecular characteristics of disorders associated with these difficult-to-detect expansions.

A germ cell-specific ageing pattern in otherwise healthy men.

Life-long sperm production leads to the assumption that male fecundity remains unchanged throughout life. However, recently it was shown that paternal age has profound consequences for male fertility and offspring health. Paternal age effects are caused by an accumulation of germ cell mutations over time, causing severe congenital diseases. Apart from these well-described cases, molecular patterns of ageing in germ cells and their impact on DNA integrity have not been studied in detail. In this study, we aimed to assess the effects of 'pure' ageing on male reproductive health and germ cell quality. We assembled a cohort of 198 healthy men (18-84 years) for which end points such as semen and hormone profiles, sexual health and well-being, and sperm DNA parameters were evaluated. Sperm production and hormonal profiles were maintained at physiological levels over a period of six decades. In contrast, we identified a germ cell-specific ageing pattern characterized by a steady increase of telomere length in sperm and a sharp increase in sperm DNA instability, particularly after the sixth decade. Importantly, we found sperm DNA methylation changes in 236 regions, mostly nearby genes associated with neuronal development. By in silico analysis, we found that 10 of these regions are located in loci which can potentially escape the first wave of genome-wide demethylation after fertilization. In conclusion, human male germ cells present a unique germline-specific ageing process, which likely results in diminished fecundity in elderly men and poorer health prognosis for their offspring.

wg-blimp: an end-to-end analysis pipeline for whole genome bisulfite sequencing data.

BACKGROUND: Analysing whole genome bisulfite sequencing datasets is a data-intensive task that requires comprehensive and reproducible workflows to generate valid results. While many algorithms have been developed for tasks such as alignment, comprehensive end-to-end pipelines are still sparse. Furthermore, previous pipelines lack features or show technical deficiencies, thus impeding analyses. RESULTS: We developed wg-blimp (whole genome bisulfite sequencing methylation analysis pipeline) as an end-to-end pipeline to ease whole genome bisulfite sequencing data analysis. It integrates established algorithms for alignment, quality control, methylation calling, detection of differentially methylated regions, and methylome segmentation, requiring only a reference genome and raw sequencing data as input. Comparing wg-blimp to previous end-to-end pipelines reveals similar setups for common sequence processing tasks, but shows differences for post-alignment analyses. We improve on previous pipelines by providing a more comprehensive analysis workflow as well as an interactive user interface. To demonstrate wg-blimp's ability to produce correct results we used it to call differentially methylated regions for two publicly available datasets. We were able to replicate 112 of 114 previously published regions, and found results to be consistent with previous findings. We further applied wg-blimp to a publicly available sample of embryonic stem cells to showcase methylome segmentation. As expected, unmethylated regions were in close proximity of transcription start sites. Segmentation results were consistent with previous analyses, despite different reference genomes and sequencing techniques. CONCLUSIONS: wg-blimp provides a comprehensive analysis pipeline for whole genome bisulfite sequencing data as well as a user interface for simplified result inspection. We demonstrated its applicability by analysing multiple publicly available datasets. Thus, wg-blimp is a relevant alternative to previous analysis pipelines and may facilitate future epigenetic research.

Unstable TTTTA/TTTCA expansions in MARCH6 are associated with Familial Adult Myoclonic Epilepsy type 3.

Familial Adult Myoclonic Epilepsy (FAME) is a genetically heterogeneous disorder characterized by cortical tremor and seizures. Intronic TTTTA/TTTCA repeat expansions in SAMD12 (FAME1) are the main cause of FAME in Asia. Using genome sequencing and repeat-primed PCR, we identify another site of this repeat expansion, in MARCH6 (FAME3) in four European families. Analysis of single DNA molecules with nanopore sequencing and molecular combing show that expansions range from 3.3 to 14 kb on average. However, we observe considerable variability in expansion length and structure, supporting the existence of multiple expansion configurations in blood cells and fibroblasts of the same individual. Moreover, the largest expansions are associated with micro-rearrangements occurring near the expansion in 20% of cells. This study provides further evidence that FAME is caused by intronic TTTTA/TTTCA expansions in distinct genes and reveals that expansions exhibit an unexpectedly high somatic instability that can ultimately result in genomic rearrangements.

Genome-Wide Analysis of the Nucleosome Landscape in Individuals with Coffin-Siris Syndrome.

The switch/sucrose non-fermenting (SWI/SNF) complex is an ATP-dependent chromatin remodeller that regulates the spacing of nucleosomes and thereby controls gene expression. Heterozygous mutations in genes encoding subunits of the SWI/SNF complex have been reported in individuals with Coffin-Siris syndrome (CSS), with the majority of the mutations in ARID1B. CSS is a rare congenital disorder characterized by facial dysmorphisms, digital anomalies, and variable intellectual disability. We hypothesized that mutations in genes encoding subunits of the ubiquitously expressed SWI/SNF complex may lead to alterations of the nucleosome profiles in different cell types. We performed the first study on CSS-patient samples and investigated the nucleosome landscapes of cell-free DNA (cfDNA) isolated from blood plasma by whole-genome sequencing. In addition, we studied the nucleosome landscapes of CD14+ monocytes from CSS-affected individuals by nucleosome occupancy and methylome-sequencing (NOMe-seq) as well as their expression profiles. In cfDNA of CSS-affected individuals with heterozygous ARID1B mutations, we did not observe major changes in the nucleosome profile around transcription start sites. In CD14+ monocytes, we found few genomic regions with different nucleosome occupancy when compared to controls. RNA-seq analysis of CD14+ monocytes of these individuals detected only few differentially expressed genes, which were not in proximity to any of the identified differential nucleosome-depleted regions. In conclusion, we show that heterozygous mutations in the human SWI/SNF subunit ARID1B do not have a major impact on the nucleosome landscape or gene expression in blood cells. This might be due to functional redundancy, cell-type specificity, or alternative functions of ARID1B.

A hybrid parameter estimation algorithm for beta mixtures and applications to methylation state classification.

BACKGROUND: Mixtures of beta distributions are a flexible tool for modeling data with values on the unit interval, such as methylation levels. However, maximum likelihood parameter estimation with beta distributions suffers from problems because of singularities in the log-likelihood function if some observations take the values 0 or 1. METHODS: While ad-hoc corrections have been proposed to mitigate this problem, we propose a different approach to parameter estimation for beta mixtures where such problems do not arise in the first place. Our algorithm combines latent variables with the method of moments instead of maximum likelihood, which has computational advantages over the popular EM algorithm. RESULTS: As an application, we demonstrate that methylation state classification is more accurate when using adaptive thresholds from beta mixtures than non-adaptive thresholds on observed methylation levels. We also demonstrate that we can accurately infer the number of mixture components. CONCLUSIONS: The hybrid algorithm between likelihood-based component un-mixing and moment-based parameter estimation is a robust and efficient method for beta mixture estimation. We provide an implementation of the method ("betamix") as open source software under the MIT license.

Regions of common inter-individual DNA methylation differences in human monocytes: genetic basis and potential function.

BACKGROUND: There is increasing evidence for inter-individual methylation differences at CpG dinucleotides in the human genome, but the regional extent and function of these differences have not yet been studied in detail. For identifying regions of common methylation differences, we used whole genome bisulfite sequencing data of monocytes from five donors and a novel bioinformatic strategy. RESULTS: We identified 157 differentially methylated regions (DMRs) with four or more CpGs, almost none of which has been described before. The DMRs fall into different chromatin states, where methylation is inversely correlated with active, but not repressive histone marks. However, methylation is not correlated with the expression of associated genes. High-resolution single nucleotide polymorphism (SNP) genotyping of the five donors revealed evidence for a role of cis-acting genetic variation in establishing methylation patterns. To validate this finding in a larger cohort, we performed genome-wide association studies (GWAS) using SNP genotypes and 450k array methylation data from blood samples of 1128 individuals. Only 30/157 (19%) DMRs include at least one 450k CpG, which shows that these arrays miss a large proportion of DNA methylation variation. In most cases, the GWAS peak overlapped the CpG position, and these regions are enriched for CREB group, NF-1, Sp100 and CTCF binding motifs. In two cases, there was tentative evidence for a trans-effect by KRAB zinc finger proteins. CONCLUSIONS: Allele-specific DNA methylation occurs in discrete chromosomal regions and is driven by genetic variation in cis and trans, but in general has little effect on gene expression.

Epigenetic dynamics of monocyte-to-macrophage differentiation.

BACKGROUND: Monocyte-to-macrophage differentiation involves major biochemical and structural changes. In order to elucidate the role of gene regulatory changes during this process, we used high-throughput sequencing to analyze the complete transcriptome and epigenome of human monocytes that were differentiated in vitro by addition of colony-stimulating factor 1 in serum-free medium. RESULTS: Numerous mRNAs and miRNAs were significantly up- or down-regulated. More than 100 discrete DNA regions, most often far away from transcription start sites, were rapidly demethylated by the ten eleven translocation enzymes, became nucleosome-free and gained histone marks indicative of active enhancers. These regions were unique for macrophages and associated with genes involved in the regulation of the actin cytoskeleton, phagocytosis and innate immune response. CONCLUSIONS: In summary, we have discovered a phagocytic gene network that is repressed by DNA methylation in monocytes and rapidly de-repressed after the onset of macrophage differentiation.

Mutational dynamics between primary and relapse neuroblastomas.

Neuroblastoma is a malignancy of the developing sympathetic nervous system that is often lethal when relapse occurs. We here used whole-exome sequencing, mRNA expression profiling, array CGH and DNA methylation analysis to characterize 16 paired samples at diagnosis and relapse from individuals with neuroblastoma. The mutational burden significantly increased in relapsing tumors, accompanied by altered mutational signatures and reduced subclonal heterogeneity. Global allele frequencies at relapse indicated clonal mutation selection during disease progression. Promoter methylation patterns were consistent over disease course and were patient specific. Recurrent alterations at relapse included mutations in the putative CHD5 neuroblastoma tumor suppressor, chromosome 9p losses, DOCK8 mutations, inactivating mutations in PTPN14 and a relapse-specific activity pattern for the PTPN14 target YAP. Recurrent new mutations in HRAS, KRAS and genes mediating cell-cell interaction in 13 of 16 relapse tumors indicate disturbances in signaling pathways mediating mesenchymal transition. Our data shed light on genetic alteration frequency, identity and evolution in neuroblastoma.

Evolutionary origin and methylation status of human intronic CpG islands that are not present in mouse.

Imprinting of the human RB1 gene is due to the presence of a differentially methylated CpG island (CGI) in intron 2, which is part of a retrocopy derived from the PPP1R26 gene on chromosome 9. The murine Rb1 gene does not have this retrocopy and is not imprinted. We have investigated whether the RB1/Rb1 locus is unique with respect to these differences. For this, we have compared the CGIs from human and mouse by in silico analyses. We have found that the human genome does not only contain more CGIs than the mouse, but the proportion of intronic CGIs is also higher (7.7% vs. 3.5%). At least 2,033 human intronic CGIs are not present in the mouse. Among these CGIs, 104 show sequence similarities elsewhere in the human genome, which suggests that they arose from retrotransposition. We could narrow down the time points when most of these CGIs appeared during evolution. Their methylation status was analyzed in two monocyte methylome data sets from whole-genome bisulfite sequencing and in 18 published methylomes. Four CGIs, which are located in the RB1, ASRGL1, PARP11, and PDXDC1 genes, occur as methylated and unmethylated copies. In contrast to imprinted methylation at the RB1 locus, differential methylation of the ASRGL1 and PDXDC1 CGIs appears to be sequence dependent. Our study supports the notion that the epigenetic fate of the retrotransposed DNA depends on its sequence and selective forces at the integration site.

Mutation in the AP4B1 gene cause hereditary spastic paraplegia type 47 (SPG47) .

We recently identified a new locus for spastic paraplegia type 47 (SPG47) in a consanguineous Arabic family with two affected siblings with progressive spastic paraparesis,intellectual disability, seizures, periventricular white matter changes and thin corpus callosum. Using exome sequencing, we now identified a novel AP4B1 frameshift mutation (c.664delC) in this family. This mutation was homozygous in both affected siblings and heterozygous in both parents. The mutant allele was absent in 316 Caucasian and 200 ethnically matched control chromosomes. We propose that AP4B1 mutations cause SPG47 and should be considered in early onset spastic paraplegia with intellectual disability.