PARKIN is not required to sustain OXPHOS function in adult mammalian tissues.

Loss-of-function variants in the PRKN gene encoding the ubiquitin E3 ligase PARKIN cause autosomal recessive early-onset Parkinson's disease (PD). Extensive in vitro and in vivo studies have reported that PARKIN is involved in multiple pathways of mitochondrial quality control, including mitochondrial degradation and biogenesis. However, these findings are surrounded by substantial controversy due to conflicting experimental data. In addition, the existing PARKIN-deficient mouse models have failed to faithfully recapitulate PD phenotypes. Therefore, we have investigated the mitochondrial role of PARKIN during ageing and in response to stress by employing a series of conditional Parkin knockout mice. We report that PARKIN loss does not affect oxidative phosphorylation (OXPHOS) capacity and mitochondrial DNA (mtDNA) levels in the brain, heart, and skeletal muscle of aged mice. We also demonstrate that PARKIN deficiency does not exacerbate the brain defects and the pro-inflammatory phenotype observed in mice carrying high levels of mtDNA mutations. To rule out compensatory mechanisms activated during embryonic development of Parkin-deficient mice, we generated a mouse model where loss of PARKIN was induced in adult dopaminergic (DA) neurons. Surprisingly, also these mice did not show motor impairment or neurodegeneration, and no major transcriptional changes were found in isolated midbrain DA neurons. Finally, we report a patient with compound heterozygous PRKN pathogenic variants that lacks PARKIN and has developed PD. The PARKIN deficiency did not impair OXPHOS activities or induce mitochondrial pathology in skeletal muscle from the patient. Altogether, our results argue that PARKIN is dispensable for OXPHOS function in adult mammalian tissues.

PRC2-mediated repression is essential to maintain identity and function of differentiated dopaminergic and serotonergic neurons.

How neurons can maintain cellular identity over an entire life span remains largely unknown. Here, we show that maintenance of identity in differentiated dopaminergic and serotonergic neurons is critically reliant on the Polycomb repressive complex 2 (PRC2). Deletion of the obligate PRC2 component, Eed, in these neurons resulted in global loss of H3K27me3, followed by a gradual activation of genes harboring both H3K27me3 and H3K9me3 modifications. Notably, H3K9me3 was lost at these PRC2 targets before gene activation. Neuronal survival was not compromised; instead, there was a reduction in subtype-specific gene expression and a progressive impairment of dopaminergic and serotonergic neuronal function, leading to behavioral deficits characteristic of Parkinson's disease and anxiety. Single-cell analysis revealed subtype-specific vulnerability to loss of PRC2 repression in dopamine neurons of the substantia nigra. Our study reveals that a PRC2-dependent nonpermissive chromatin state is essential to maintain the subtype identity and function of dopaminergic and serotonergic neurons.

Environmentally induced DNA methylation is inherited across generations in an aquatic keystone species.

Transgenerational inheritance of environmentally induced epigenetic marks can have significant impacts on eco-evolutionary dynamics, but the phenomenon remains controversial in ecological model systems. We used whole-genome bisulfite sequencing of individual water fleas (Daphnia magna) to assess whether environmentally induced DNA methylation is transgenerationally inherited. Genetically identical females were exposed to one of three natural stressors, or a de-methylating drug, and their offspring were propagated clonally for four generations under control conditions. We identified between 70 and 225 differentially methylated CpG positions (DMPs) in F1 individuals whose mothers were exposed to a natural stressor. Roughly half of these environmentally induced DMPs persisted until generation F4. In contrast, treatment with the drug demonstrated that pervasive hypomethylation upon exposure is reset almost completely after one generation. These results suggest that environmentally induced DNA methylation is non-random and stably inherited across generations in Daphnia, making epigenetic inheritance a putative factor in the eco-evolutionary dynamics of freshwater communities.

Mitochondrial dysfunction in adult midbrain dopamine neurons triggers an early immune response.

Dopamine (DA) neurons of the midbrain are at risk to become affected by mitochondrial damage over time and mitochondrial defects have been frequently reported in Parkinson's disease (PD) patients. However, the causal contribution of adult-onset mitochondrial dysfunction to PD remains uncertain. Here, we developed a mouse model lacking Mitofusin 2 (MFN2), a key regulator of mitochondrial network homeostasis, in adult midbrain DA neurons. The knockout mice develop severe and progressive DA neuron-specific mitochondrial dysfunction resulting in neurodegeneration and parkinsonism. To gain further insights into pathophysiological events, we performed transcriptomic analyses of isolated DA neurons and found that mitochondrial dysfunction triggers an early onset immune response, which precedes mitochondrial swelling, mtDNA depletion, respiratory chain deficiency and cell death. Our experiments show that the immune response is an early pathological event when mitochondrial dysfunction is induced in adult midbrain DA neurons and that neuronal death may be promoted non-cell autonomously by the cross-talk and activation of surrounding glial cells.

Mutational patterns and clonal evolution from diagnosis to relapse in pediatric acute lymphoblastic leukemia.

The mechanisms driving clonal heterogeneity and evolution in relapsed pediatric acute lymphoblastic leukemia (ALL) are not fully understood. We performed whole genome sequencing of samples collected at diagnosis, relapse(s) and remission from 29 Nordic patients. Somatic point mutations and large-scale structural variants were called using individually matched remission samples as controls, and allelic expression of the mutations was assessed in ALL cells using RNA-sequencing. We observed an increased burden of somatic mutations at relapse, compared to diagnosis, and at second relapse compared to first relapse. In addition to 29 known ALL driver genes, of which nine genes carried recurrent protein-coding mutations in our sample set, we identified putative non-protein coding mutations in regulatory regions of seven additional genes that have not previously been described in ALL. Cluster analysis of hundreds of somatic mutations per sample revealed three distinct evolutionary trajectories during ALL progression from diagnosis to relapse. The evolutionary trajectories provide insight into the mutational mechanisms leading relapse in ALL and could offer biomarkers for improved risk prediction in individual patients.

Tissue-specific transcriptional imprinting and heterogeneity in human innate lymphoid cells revealed by full-length single-cell RNA-sequencing.

The impact of the microenvironment on innate lymphoid cell (ILC)-mediated immunity in humans remains largely unknown. Here we used full-length Smart-seq2 single-cell RNA-sequencing to unravel tissue-specific transcriptional profiles and heterogeneity of CD127+ ILCs across four human tissues. Correlation analysis identified gene modules characterizing the migratory properties of tonsil and blood ILCs, and signatures of tissue-residency, activation and modified metabolism in colon and lung ILCs. Trajectory analysis revealed potential differentiation pathways from circulating and tissue-resident naïve ILCs to a spectrum of mature ILC subsets. In the lung we identified both CRTH2+ and CRTH2- ILC2 with lung-specific signatures, which could be recapitulated by alarmin-exposure of circulating ILC2. Finally, we describe unique TCR-V(D)J-rearrangement patterns of blood ILC1-like cells, revealing a subset of potentially immature ILCs with TCR-δ rearrangement. Our study provides a useful resource for in-depth understanding of ILC-mediated immunity in humans, with implications for disease.

Author Correction: Mutational and putative neoantigen load predict clinical benefit of adoptive T cell therapy in melanoma.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Analysis of DNA methylation patterns in the tumor immune microenvironment of metastatic melanoma.

The presence of immune cells in the tumor microenvironment has been associated with response to immunotherapies across several cancer types, including melanoma. Despite its therapeutic relevance, characterization of the melanoma immune microenvironments remains insufficiently explored. To distinguish the immune microenvironment in a cohort of 180 metastatic melanoma clinical specimens, we developed a method using promoter CpG methylation of immune cell type-specific genes extracted from genome-wide methylation arrays. Unsupervised clustering identified three immune methylation clusters with varying levels of immune CpG methylation that are related to patient survival. Matching protein and gene expression data further corroborated the identified epigenetic characterization. Exploration of the possible immune exclusion mechanisms at play revealed likely dependency on MITF protein level and PTEN loss-of-function events for melanomas unresponsive to immunotherapies (immune-low). To understand whether melanoma tumors resemble other solid tumors in terms of immune methylation characteristics, we explored 15 different solid tumor cohorts from TCGA. Low-dimensional projection based on immune cell type-specific methylation revealed grouping of the solid tumors in line with melanoma immune methylation clusters rather than tumor types. Association of survival outcome with immune cell type-specific methylation differed across tumor and cell types. However, in melanomas immune cell type-specific methylation was associated with inferior patient survival. Exploration of the immune methylation patterns in a pan-cancer context suggested that specific immune microenvironments might occur across the cancer spectrum. Together, our findings underscore the existence of diverse immune microenvironments, which may be informative for future immunotherapeutic applications.

Partially methylated domains are hypervariable in breast cancer and fuel widespread CpG island hypermethylation.

Global loss of DNA methylation and CpG island (CGI) hypermethylation are key epigenomic aberrations in cancer. Global loss manifests itself in partially methylated domains (PMDs) which extend up to megabases. However, the distribution of PMDs within and between tumor types, and their effects on key functional genomic elements including CGIs are poorly defined. We comprehensively show that loss of methylation in PMDs occurs in a large fraction of the genome and represents the prime source of DNA methylation variation. PMDs are hypervariable in methylation level, size and distribution, and display elevated mutation rates. They impose intermediate DNA methylation levels incognizant of functional genomic elements including CGIs, underpinning a CGI methylator phenotype (CIMP). Repression effects on tumor suppressor genes are negligible as they are generally excluded from PMDs. The genomic distribution of PMDs reports tissue-of-origin and may represent tissue-specific silent regions which tolerate instability at the epigenetic, transcriptomic and genetic level.

Author Correction: Landscape of somatic mutations in 560 breast cancer whole-genome sequences.

In the Methods section of this Article, 'greater than' should have been 'less than' in the sentence 'Putative regions of clustered rearrangements were identified as having an average inter-rearrangement distance that was at least 10 times greater than the whole-genome average for the individual sample. '. The Article has not been corrected.

Spatially and functionally distinct subclasses of breast cancer-associated fibroblasts revealed by single cell RNA sequencing.

Cancer-associated fibroblasts (CAFs) are a major constituent of the tumor microenvironment, although their origin and roles in shaping disease initiation, progression and treatment response remain unclear due to significant heterogeneity. Here, following a negative selection strategy combined with single-cell RNA sequencing of 768 transcriptomes of mesenchymal cells from a genetically engineered mouse model of breast cancer, we define three distinct subpopulations of CAFs. Validation at the transcriptional and protein level in several experimental models of cancer and human tumors reveal spatial separation of the CAF subclasses attributable to different origins, including the peri-vascular niche, the mammary fat pad and the transformed epithelium. Gene profiles for each CAF subtype correlate to distinctive functional programs and hold independent prognostic capability in clinical cohorts by association to metastatic disease. In conclusion, the improved resolution of the widely defined CAF population opens the possibility for biomarker-driven development of drugs for precision targeting of CAFs.

Author Correction: A comprehensive map coupling histone modifications with gene regulation in adult dopaminergic and serotonergic neurons.

In the originally published version of this Article, financial support was not fully acknowledged. The PDF and HTML versions of the Article have now been corrected to include support to Thomas Perlmann provided by Knut and Alice Wallenberg Foundation (grant 2013.0075) and Swedish Research Council (VR; grant 2016-02506).

A comprehensive map coupling histone modifications with gene regulation in adult dopaminergic and serotonergic neurons.

The brain is composed of hundreds of different neuronal subtypes, which largely retain their identity throughout the lifespan of the organism. The mechanisms governing this stability are not fully understood, partly due to the diversity and limited size of clinically relevant neuronal populations, which constitute a technical challenge for analysis. Here, using a strategy that allows for ChIP-seq combined with RNA-seq in small neuronal populations in vivo, we present a comparative analysis of permissive and repressive histone modifications in adult midbrain dopaminergic neurons, raphe nuclei serotonergic neurons, and embryonic neural progenitors. Furthermore, we utilize the map generated by our analysis to show that the transcriptional response of midbrain dopaminergic neurons following 6-OHDA or methamphetamine injection is characterized by increased expression of genes with promoters dually marked by H3K4me3/H3K27me3. Our study provides an in vivo genome-wide analysis of permissive/repressive histone modifications coupled to gene expression in these rare neuronal subtypes.

Mutational and putative neoantigen load predict clinical benefit of adoptive T cell therapy in melanoma.

Adoptive T-cell therapy (ACT) is a highly intensive immunotherapy regime that has yielded remarkable response rates and many durable responses in clinical trials in melanoma; however, 50-60% of the patients have no clinical benefit. Here, we searched for predictive biomarkers to ACT in melanoma. Whole exome- and transcriptome sequencing and neoantigen prediction were applied to pre-treatment samples from 27 patients recruited to a clinical phase I/II trial of ACT in stage IV melanoma. All patients had previously progressed on other immunotherapies. We report that clinical benefit is associated with significantly higher predicted neoantigen load. High mutation and predicted neoantigen load are significantly associated with improved progression-free and overall survival. Further, clinical benefit is associated with the expression of immune activation signatures including a high MHC-I antigen processing and presentation score. These results improve our understanding of mechanisms behind clinical benefit of ACT in melanoma.

Breast cancer genome and transcriptome integration implicates specific mutational signatures with immune cell infiltration.

A recent comprehensive whole genome analysis of a large breast cancer cohort was used to link known and novel drivers and substitution signatures to the transcriptome of 266 cases. Here, we validate that subtype-specific aberrations show concordant expression changes for, for example, TP53, PIK3CA, PTEN, CCND1 and CDH1. We find that CCND3 expression levels do not correlate with amplification, while increased GATA3 expression in mutant GATA3 cancers suggests GATA3 is an oncogene. In luminal cases the total number of substitutions, irrespective of type, associates with cell cycle gene expression and adverse outcome, whereas the number of mutations of signatures 3 and 13 associates with immune-response specific gene expression, increased numbers of tumour-infiltrating lymphocytes and better outcome. Thus, while earlier reports imply that the sheer number of somatic aberrations could trigger an immune-response, our data suggests that substitutions of a particular type are more effective in doing so than others.

Consensus of gene expression phenotypes and prognostic risk predictors in primary lung adenocarcinoma.

Transcriptional profiling of lung adenocarcinomas has identified numerous gene expression phenotype (GEP) and risk prediction (RP) signatures associated with patient outcome. However, classification agreement between signatures, underlying transcriptional programs, and independent signature validation are less studied. We classified 2395 transcriptional adenocarcinoma profiles, assembled from 17 public cohorts, using 11 GEP and seven RP signatures, finding that 16 signatures were associated with patient survival in the total cohort and in multiple individual cohorts. For significant signatures, total cohort hazard ratios were ~2 in univariate analyses (mean=1.95, range=1.4-2.6). Strong classification agreement between signatures was observed, especially for predicted low-risk patients by adenocarcinoma-derived signatures. Expression of proliferation-related genes correlated strongly with GEP subtype classifications and RP scores, driving the gene signature association with prognosis. A three-group consensus definition of samples across 10 GEP classifiers demonstrated aggregation of samples with specific smoking patterns, gender, and EGFR/KRAS mutations, while survival differences were only significant when patients were divided into low- or high-risk. In summary, our study demonstrates a consensus between GEPs and RPs in lung adenocarcinoma through a common underlying transcriptional program. This consensus generalizes reported problems with current signatures in a clinical context, stressing development of new adenocarcinoma-specific single sample predictors for clinical use.

Landscape of somatic mutations in 560 breast cancer whole-genome sequences.

We analysed whole-genome sequences of 560 breast cancers to advance understanding of the driver mutations conferring clonal advantage and the mutational processes generating somatic mutations. We found that 93 protein-coding cancer genes carried probable driver mutations. Some non-coding regions exhibited high mutation frequencies, but most have distinctive structural features probably causing elevated mutation rates and do not contain driver mutations. Mutational signature analysis was extended to genome rearrangements and revealed twelve base substitution and six rearrangement signatures. Three rearrangement signatures, characterized by tandem duplications or deletions, appear associated with defective homologous-recombination-based DNA repair: one with deficient BRCA1 function, another with deficient BRCA1 or BRCA2 function, the cause of the third is unknown. This analysis of all classes of somatic mutation across exons, introns and intergenic regions highlights the repertoire of cancer genes and mutational processes operating, and progresses towards a comprehensive account of the somatic genetic basis of breast cancer.

An integrated genomics analysis of epigenetic subtypes in human breast tumors links DNA methylation patterns to chromatin states in normal mammary cells.

BACKGROUND: Aberrant DNA methylation is frequently observed in breast cancer. However, the relationship between methylation patterns and the heterogeneity of breast cancer has not been comprehensively characterized. METHODS: Whole-genome DNA methylation analysis using Illumina Infinium HumanMethylation450 BeadChip arrays was performed on 188 human breast tumors. Unsupervised bootstrap consensus clustering was performed to identify DNA methylation epigenetic subgroups (epitypes). The Cancer Genome Atlas data, including methylation profiles of 669 human breast tumors, was used for validation. The identified epitypes were characterized by integration with publicly available genome-wide data, including gene expression levels, DNA copy numbers, whole-exome sequencing data, and chromatin states. RESULTS: We identified seven breast cancer epitypes. One epitype was distinctly associated with basal-like tumors and with BRCA1 mutations, one epitype contained a subset of ERBB2-amplified tumors characterized by multiple additional amplifications and the most complex genomes, and one epitype displayed a methylation profile similar to normal epithelial cells. Luminal tumors were stratified into the remaining four epitypes, with differences in promoter hypermethylation, global hypomethylation, proliferative rates, and genomic instability. Specific hyper- and hypomethylation across the basal-like epitype was rare. However, we observed that the candidate genomic instability drivers BRCA1 and HORMAD1 displayed aberrant methylation linked to gene expression levels in some basal-like tumors. Hypomethylation in luminal tumors was associated with DNA repeats and subtelomeric regions. We observed two dominant patterns of aberrant methylation in breast cancer. One pattern, constitutively methylated in both basal-like and luminal breast cancer, was linked to genes with promoters in a Polycomb-repressed state in normal epithelial cells and displayed no correlation with gene expression levels. The second pattern correlated with gene expression levels and was associated with methylation in luminal tumors and genes with active promoters in normal epithelial cells. CONCLUSIONS: Our results suggest that hypermethylation patterns across basal-like breast cancer may have limited influence on tumor progression and instead reflect the repressed chromatin state of the tissue of origin. On the contrary, hypermethylation patterns specific to luminal breast cancer influence gene expression, may contribute to tumor progression, and may present an actionable epigenetic alteration in a subset of luminal breast cancers.

Prognostic and Chemotherapy Predictive Value of Gene-Expression Phenotypes in Primary Lung Adenocarcinoma.

PURPOSE: Primary lung adenocarcinoma remains a deadly disease. Gene-expression phenotypes (GEPs) in adenocarcinoma have potential to provide clinically relevant disease stratification for improved prognosis and treatment prediction, given appropriate clinical and methodologic validation. EXPERIMENTAL DESIGN: 2,395 transcriptional adenocarcinoma profiles were assembled from 17 public cohorts and classified by a nearest centroid GEP classifier into three subtypes: terminal respiratory unit (TRU), proximal-proliferative, and proximal-inflammatory, and additionally scored by five transcriptional metagenes representing different biologic processes, including proliferation. Prognostic- and chemotherapy-predictive associations of the subtypes were analyzed by univariate and multivariate analysis using overall survival or distant metastasis-free survival as endpoints. RESULTS: Overall, GEPs were associated with patient outcome in both univariate and multivariate analyses, although not in all individual cohorts. The prognostically relevant division was between TRU- and non-TRU-classified cases, with expression of proliferation-associated genes as a key prognostic component. In contrast, GEP classification was not predictive of adjuvant chemotherapy response. GEP classification showed stability to random perturbations of genes or samples and alterations to classification procedures (typically <10% of cases/cohort switching subtype). High classification variability (>20% of cases switching subtype) was observed when removing larger or entire fractions of a single subtype, due to gene-centering shifts not addressable by the classifier. CONCLUSIONS: In a large-scale evaluation, we show that GEPs add prognostic value to standard clinicopathologic variables in lung adenocarcinoma. Subject to classifier refinement and confirmation in prospective cohorts, GEPs have potential to affect the prognostication of adenocarcinoma patients through a molecularly driven disease stratification.

DNA methylation subgroups in melanoma are associated with proliferative and immunological processes.

BACKGROUND: DNA methylation at CpG dinucleotides is modified in tumorigenesis with potential impact on transcriptional activity. METHODS: We used the Illumina 450 K platform to evaluate DNA methylation patterns of 50 metastatic melanoma tumors, with matched gene expression data. RESULTS: We identified three different methylation groups and validated the groups in independent data from The Cancer Genome Atlas. One group displayed hypermethylation of a developmental promoter set, genome-wide demethylation, increased proliferation and activity of the SWI/SNF complex. A second group had a methylation pattern resembling stromal and leukocyte cells, over-expressed an immune signature and had improved survival rates in metastatic tumors (p < 0.05). A third group had intermediate methylation levels and expressed both proliferative and immune signatures. The methylation groups corresponded to some degree with previously identified gene expression phenotypes. CONCLUSIONS: Melanoma consists of divergent methylation groups that are distinguished by promoter methylation, proliferation and content of immunological cells.