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.

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.

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.

MicroRNA-200c-141 and ∆Np63 are required for breast epithelial differentiation and branching morphogenesis.

The epithelial compartment of the breast contains two lineages, the luminal- and the myoepithelial cells. D492 is a breast epithelial cell line with stem cell properties that forms branching epithelial structures in 3D culture with both luminal- and myoepithelial differentiation. We have recently shown that D492 undergo epithelial to mesenchymal transition (EMT) when co-cultured with endothelial cells. This 3D co-culture model allows critical analysis of breast epithelial lineage development and EMT. In this study, we compared the microRNA (miR) expression profiles for D492 and its mesenchymal-derivative D492M. Suppression of the miR-200 family in D492M was among the most profound changes observed. Exogenous expression of miR-200c-141 in D492M reversed the EMT phenotype resulting in gain of luminal but not myoepithelial differentiation. In contrast, forced expression of ∆Np63 in D492M restored the myoepithelial phenotype only. Co-expression of miR-200c-141 and ∆Np63 in D492M restored the branching morphogenesis in 3D culture underlining the requirement for both luminal and myoepithelial elements for obtaining full branching morphogenesis in breast epithelium. Introduction of a miR-200c-141 construct in both D492 and D492M resulted in resistance to endothelial induced EMT. In conclusion, our data suggests that expression of miR-200c-141 and ∆Np63 in D492M can reverse EMT resulting in luminal- and myoepithelial differentiation, respectively, demonstrating the importance of these molecules in epithelial integrity in the human breast.

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.

Exosomes reflect the hypoxic status of glioma cells and mediate hypoxia-dependent activation of vascular cells during tumor development.

Hypoxia, or low oxygen tension, is a major regulator of tumor development and aggressiveness. However, how cancer cells adapt to hypoxia and communicate with their surrounding microenvironment during tumor development remain important questions. Here, we show that secreted vesicles with exosome characteristics mediate hypoxia-dependent intercellular signaling of the highly malignant brain tumor glioblastoma multiforme (GBM). In vitro hypoxia experiments with glioma cells and studies with patient materials reveal the enrichment in exosomes of hypoxia-regulated mRNAs and proteins (e.g., matrix metalloproteinases, IL-8, PDGFs, caveolin 1, and lysyl oxidase), several of which were associated with poor glioma patient prognosis. We show that exosomes derived from GBM cells grown at hypoxic compared with normoxic conditions are potent inducers of angiogenesis ex vivo and in vitro through phenotypic modulation of endothelial cells. Interestingly, endothelial cells were programmed by GBM cell-derived hypoxic exosomes to secrete several potent growth factors and cytokines and to stimulate pericyte PI3K/AKT signaling activation and migration. Moreover, exosomes derived from hypoxic compared with normoxic conditions showed increased autocrine, promigratory activation of GBM cells. These findings were correlated with significantly enhanced induction by hypoxic compared with normoxic exosomes of tumor vascularization, pericyte vessel coverage, GBM cell proliferation, as well as decreased tumor hypoxia in a mouse xenograft model. We conclude that the proteome and mRNA profiles of exosome vesicles closely reflect the oxygenation status of donor glioma cells and patient tumors, and that the exosomal pathway constitutes a potentially targetable driver of hypoxia-dependent intercellular signaling during tumor development.

DNA methylation and histone modifications regulate SOX11 expression in lymphoid and solid cancer cells.

BACKGROUND: The neural transcription factor SOX11 is present at specific stages during embryo development with a very restricted expression in adult tissue, indicating precise regulation of transcription. SOX11 is strongly up-regulated in some malignancies and have a functional role in tumorgenesis. With the aim to explore differences in epigenetic regulation of SOX11 expression in normal versus neoplastic cells, we investigated methylation and histone modifications related to the SOX11 promoter and the possibility to induce re-expression using histone deacetylase (HDAC) or EZH2 inhibitors. METHODS: The epigenetic regulation of SOX11 was investigated in distinct non-malignant cell populations (n = 7) and neoplastic cell-lines (n = 42) of different cellular origins. DNA methylation was assessed using bisulfite sequencing, methylation-specific melting curve analysis, MethyLight and pyrosequencing. The presence of H3K27me3 was assessed using ChIP-qPCR. The HDAC inhibitors Vorinostat and trichostatin A were used to induce SOX11 in cell lines with no endogenous expression. RESULTS: The SOX11 promoter shows a low degree of methylation and strong enrichment of H3K27me3 in non-malignant differentiated cells, independent of cellular origin. Cancers of the B-cell lineage are strongly marked by de novo methylation at the SOX11 promoter in SOX11 non-expressing cells, while solid cancer entities display a more varying degree of SOX11 promoter methylation. The silencing mark H3K27me3 was generally present at the SOX11 promoter in non-expressing cells, and an increased enrichment was observed in cancer cells with a low degree of SOX11 methylation compared to cells with dense methylation. Finally, we demonstrate that the HDAC inhibitors (vorinostat and trichostatin A) induce SOX11 expression in cancer cells with low levels of SOX11 methylation. CONCLUSIONS: We show that SOX11 is strongly marked by repressive histone marks in non-malignant cells. In contrast, SOX11 regulation in neoplastic tissues is more complex involving both DNA methylation and histone modifications. The possibility to re-express SOX11 in non-methylated tissue is of clinical relevance, and was successfully achieved in cell lines with low levels of SOX11 methylation. In breast cancer patients, methylation of the SOX11 promoter was shown to correlate with estrogen receptor status, suggesting that SOX11 may be functionally re-expressed during treatment with HDAC inhibitors in specific patient subgroups.

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.

The gene expression landscape of breast cancer is shaped by tumor protein p53 status and epithelial-mesenchymal transition.

INTRODUCTION: Gene expression data derived from clinical cancer specimens provide an opportunity to characterize cancer-specific transcriptional programs. Here, we present an analysis delineating a correlation-based gene expression landscape of breast cancer that identifies modules with strong associations to breast cancer-specific and general tumor biology. METHODS: Modules of highly connected genes were extracted from a gene co-expression network that was constructed based on Pearson correlation, and module activities were then calculated using a pathway activity score. Functional annotations of modules were experimentally validated with an siRNA cell spot microarray system using the KPL-4 breast cancer cell line, and by using gene expression data from functional studies. Modules were derived using gene expression data representing 1,608 breast cancer samples and validated in data sets representing 971 independent breast cancer samples as well as 1,231 samples from other cancer forms. RESULTS: The initial co-expression network analysis resulted in the characterization of eight tightly regulated gene modules. Cell cycle genes were divided into two transcriptional programs, and experimental validation using an siRNA screen showed different functional roles for these programs during proliferation. The division of the two programs was found to act as a marker for tumor protein p53 (TP53) gene status in luminal breast cancer, with the two programs being separated only in luminal tumors with functional p53 (encoded by TP53). Moreover, a module containing fibroblast and stroma-related genes was highly expressed in fibroblasts, but was also up-regulated by overexpression of epithelial-mesenchymal transition factors such as transforming growth factor beta 1 (TGF-beta1) and Snail in immortalized human mammary epithelial cells. Strikingly, the stroma transcriptional program related to less malignant tumors for luminal disease and aggressive lymph node positive disease among basal-like tumors. CONCLUSIONS: We have derived a robust gene expression landscape of breast cancer that reflects known subtypes as well as heterogeneity within these subtypes. By applying the modules to TP53-mutated samples we shed light on the biological consequences of non-functional p53 in otherwise low-proliferating luminal breast cancer. Furthermore, as in the case of the stroma module, we show that the biological and clinical interpretation of a set of co-regulated genes is subtype-dependent.

Gene expression profiling of primary male breast cancers reveals two unique subgroups and identifies N-acetyltransferase-1 (NAT1) as a novel prognostic biomarker.

INTRODUCTION: Male breast cancer (MBC) is a rare and inadequately characterized disease. The aim of the present study was to characterize MBC tumors transcriptionally, to classify them into comprehensive subgroups, and to compare them with female breast cancer (FBC). METHODS: A total of 66 clinicopathologically well-annotated fresh frozen MBC tumors were analyzed using Illumina Human HT-12 bead arrays, and a tissue microarray with 220 MBC tumors was constructed for validation using immunohistochemistry. Two external gene expression datasets were used for comparison purposes: 37 MBCs and 359 FBCs. RESULTS: Using an unsupervised approach, we classified the MBC tumors into two subgroups, luminal M1 and luminal M2, respectively, with differences in tumor biological features and outcome, and which differed from the intrinsic subgroups described in FBC. The two subgroups were recapitulated in the external MBC dataset. Luminal M2 tumors were characterized by high expression of immune response genes and genes associated with estrogen receptor (ER) signaling. Luminal M1 tumors, on the other hand, despite being ER positive by immunohistochemistry showed a lower correlation to genes associated with ER signaling and displayed a more aggressive phenotype and worse prognosis. Validation of two of the most differentially expressed genes, class 1 human leukocyte antigen (HLA) and the metabolizing gene N-acetyltransferase-1 (NAT1), respectively, revealed significantly better survival associated with high expression of both markers (HLA, hazard ratio (HR) 3.6, P = 0.002; NAT1, HR 2.5, P = 0.033). Importantly, NAT1 remained significant in a multivariate analysis (HR 2.8, P = 0.040) and may thus be a novel prognostic marker in MBC. CONCLUSIONS: We have detected two unique and stable subgroups of MBC with differences in tumor biological features and outcome. They differ from the widely acknowledged intrinsic subgroups of FBC. As such, they may constitute two novel subgroups of breast cancer, occurring exclusively in men, and which may consequently require novel treatment approaches. Finally, we identified NAT1 as a possible prognostic biomarker for MBC, as suggested by NAT1 positivity corresponding to better outcome.

Characterisation of amplification patterns and target genes at chromosome 11q13 in CCND1-amplified sporadic and familial breast tumours.

Amplification of chromosomal region 11q13, containing the cell cycle regulatory gene CCND1, is frequently found in breast cancer and other malignancies. It is associated with the favourable oestrogen receptor (ER)-positive breast tumour phenotype, but also with poor prognosis and treatment failure. 11q13 spans almost 14 Mb and contains more than 200 genes and is affected by various patterns of copy number gains, suggesting complex mechanisms and selective pressure during tumour progression. In this study, we used 32 k tiling BAC array CGH to analyse 94 CCND1-amplified breast tumours from sporadic, hereditary, and familial breast cancers to fine map chromosome 11q13. A set containing 281 CCND1-non-amplified breast tumours was used for comparisons. We used gene expression data to further validate the functional effect of gene amplification. We identified six core regions covering 11q13.1-q14.1 that were amplified in different combinations. The major core contained CCND1, whereas two cores were found proximal of CCND1 and three distal. The majority of the CCND1-amplified tumours were ER-positive and classified as luminal B. Furthermore, we found that CCND1 amplification is associated with a more aggressive phenotype within histological grade 2 tumours and luminal A subtype tumours. Amplification was equally prevalent in familial and sporadic tumours, but strikingly rare in BRCA1- and BRCA2-mutated tumours. We conclude that 11q13 includes many potential target genes in addition to CCND1.

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.

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.

Landscape of somatic allelic imbalances and copy number alterations in HER2-amplified breast cancer.

INTRODUCTION: Human epidermal growth factor receptor 2 (HER2)-amplified breast cancer represents a clinically well-defined subgroup due to availability of targeted treatment. However, HER2-amplified tumors have been shown to be heterogeneous at the genomic level by genome-wide microarray analyses, pointing towards a need of further investigations for identification of recurrent copy number alterations and delineation of patterns of allelic imbalance. METHODS: High-density whole genome array-based comparative genomic hybridization (aCGH) and single nucleotide polymorphism (SNP) array data from 260 HER2-amplified breast tumors or cell lines, and 346 HER2-negative breast cancers with molecular subtype information were assembled from different repositories. Copy number alteration (CNA), loss-of-heterozygosity (LOH), copy number neutral allelic imbalance (CNN-AI), subclonal CNA and patterns of tumor DNA ploidy were analyzed using bioinformatical methods such as genomic identification of significant targets in cancer (GISTIC) and genome alteration print (GAP). The patterns of tumor ploidy were confirmed in 338 unrelated breast cancers analyzed by DNA flow cytometry with concurrent BAC aCGH and gene expression data. RESULTS: A core set of 36 genomic regions commonly affected by copy number gain or loss was identified by integrating results with a previous study, together comprising > 400 HER2-amplified tumors. While CNN-AI frequency appeared evenly distributed over chromosomes in HER2-amplified tumors, not targeting specific regions and often < 20% in frequency, the occurrence of LOH was strongly associated with regions of copy number loss. HER2-amplified and HER2-negative tumors stratified by molecular subtypes displayed different patterns of LOH and CNN-AI, with basal-like tumors showing highest frequencies followed by HER2-amplified and luminal B cases. Tumor aneuploidy was strongly associated with increasing levels of LOH, CNN-AI, CNAs and occurrence of subclonal copy number events, irrespective of subtype. Finally, SNP data from individual tumors indicated that genomic amplification in general appears as monoallelic, that is, it preferentially targets one parental chromosome in HER2-amplified tumors. CONCLUSIONS: We have delineated the genomic landscape of CNAs, amplifications, LOH, and CNN-AI in HER2-amplified breast cancer, but also demonstrated a strong association between different types of genomic aberrations and tumor aneuploidy irrespective of molecular subtype.

High-resolution genomic profiling of male breast cancer reveals differences hidden behind the similarities with female breast cancer.

Male breast cancer (MBC) is extremely rare and poorly characterized on the molecular level. Using high-resolution genomic data, we aimed to characterize MBC by genomic imbalances and to compare it with female breast cancer (FBC), and further to investigate whether the genomic profiles hold any prognostic information. Fifty-six fresh frozen MBC tumors were analyzed using high-resolution tiling BAC arrays. Significant regions in common between cases were assessed using Genomic Identification of Significant Targets in Cancer (GISTIC) analysis. A publicly available genomic data set of 359 FBC tumors was used for reference purposes. The data revealed a broad pattern of aberrations, confirming that MBC is a heterogeneous tumor type. Genomic gains were more common in MBC than in FBC and often involved whole chromosome arms, while losses of genomic material were less frequent. The most common aberrations were similar between the genders, but high-level amplifications were more common in FBC. We identified two genomic subgroups among MBCs; male-complex and male-simple. The male-complex subgroup displayed striking similarities with the previously reported luminal-complex FBC subgroup, while the male-simple subgroup seems to represent a new subgroup of breast cancer occurring only in men. There are many similarities between FBC and MBC with respect to genomic imbalances, but there are also distinct differences as revealed by high-resolution genomic profiling. MBC can be divided into two comprehensive genomic subgroups, which may be of prognostic value. The male-simple subgroup appears notably different from any genomic subgroup so far defined in FBC.

CD44 isoforms are heterogeneously expressed in breast cancer and correlate with tumor subtypes and cancer stem cell markers.

BACKGROUND: The CD44 cell adhesion molecule is aberrantly expressed in many breast tumors and has been implicated in the metastatic process as well as in the putative cancer stem cell (CSC) compartment. We aimed to investigate potential associations between alternatively spliced isoforms of CD44 and CSCs as well as to various breast cancer biomarkers and molecular subtypes. METHODS: We used q-RT-PCR and exon-exon spanning assays to analyze the expression of four alternatively spliced CD44 isoforms as well as the total expression of CD44 in 187 breast tumors and 13 cell lines. ALDH1 protein expression was determined by IHC on TMA. RESULTS: Breast cancer cell lines showed a heterogeneous expression pattern of the CD44 isoforms, which shifted considerably when cells were grown as mammospheres. Tumors characterized as positive for the CD44+/CD24- phenotype by immunohistochemistry were associated to all isoforms except the CD44 standard (CD44S) isoform, which lacks all variant exons. Conversely, tumors with strong expression of the CSC marker ALDH1 had elevated expression of CD44S. A high expression of the CD44v2-v10 isoform, which retain all variant exons, was correlated to positive steroid receptor status, low proliferation and luminal A subtype. The CD44v3-v10 isoform showed similar correlations, while high expression of CD44v8-v10 was correlated to positive EGFR, negative/low HER2 status and basal-like subtype. High expression of CD44S was associated with strong HER2 staining and also a subgroup of basal-like tumors. Unsupervised hierarchical cluster analysis of CD44 isoform expression data divided tumors into four main clusters, which showed significant correlations to molecular subtypes and differences in 10-year overall survival. CONCLUSIONS: We demonstrate that individual CD44 isoforms can be associated to different breast cancer subtypes and clinical markers such as HER2, ER and PgR, which suggests involvement of CD44 splice variants in specific oncogenic signaling pathways. Efforts to link CD44 to CSCs and tumor progression should consider the expression of various CD44 isoforms.

High Myc pathway activity and low stage of neuronal differentiation associate with poor outcome in neuroblastoma.

The childhood cancer neuroblastoma arises in the developing sympathetic nervous system and is a genotypically and phenotypically heterogeneous disease. Prognostic markers of poor survival probability include amplification of the MYCN oncogene and an undifferentiated morphology. Whereas these features discriminate high- from low-risk patients with precision, identification of poor outcome low- and intermediate-risk patients is more challenging. In this study, we analyze two large neuroblastoma microarray datasets using a priori-defined gene expression signatures. We show that differential overexpression of Myc transcriptional targets and low expression of genes involved in sympathetic neuronal differentiation predicts relapse and death from disease. This was evident not only for high-risk patients but was also robust in identifying groups of poor prognosis patients who were otherwise judged to be at low- or intermediate-risk for adverse outcome. These data suggest that pathway-specific gene expression profiling might be useful in the clinic to adjust treatment strategies for children with neuroblastoma.

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.

Genome-wide transcription factor binding site/promoter databases for the analysis of gene sets and co-occurrence of transcription factor binding motifs.

BACKGROUND: The use of global gene expression profiling is a well established approach to understand biological processes. One of the major goals of these investigations is to identify sets of genes with similar expression patterns. Such gene signatures may be very informative and reveal new aspects of particular biological processes. A logical and systematic next step is to reduce the identified gene signatures to the regulatory components that induce the relevant gene expression changes. A central issue in this context is to identify transcription factors, or transcription factor binding sites (TFBS), likely to be of importance for the expression of the gene signatures. RESULTS: We develop a strategy that efficiently produces TFBS/promoter databases based on user-defined criteria. The resulting databases constitute all genes in the Santa Cruz database and the positions for all TFBS provided by the user as position weight matrices. These databases are then used for two purposes, to identify significant TFBS in the promoters in sets of genes and to identify clusters of co-occurring TFBS. We use two criteria for significance, significantly enriched TFBS in terms of total number of binding sites for the promoters, and significantly present TFBS in terms of the fraction of promoters with binding sites. Significant TFBS are identified by a re-sampling procedure in which the query gene set is compared with typically 10(5) gene lists of similar size randomly drawn from the TFBS/promoter database. We apply this strategy to a large number of published ChIP-Chip data sets and show that the proposed approach faithfully reproduces ChIP-Chip results. The strategy also identifies relevant TFBS when analyzing gene signatures obtained from the MSigDB database. In addition, we show that several TFBS are highly correlated and that co-occurring TFBS define functionally related sets of genes. CONCLUSIONS: The presented approach of promoter analysis faithfully reproduces the results from several ChIP-Chip and MigDB derived gene sets and hence may prove to be an important method in the analysis of gene signatures obtained through ChIP-Chip or global gene expression experiments. We show that TFBS are organized in clusters of co-occurring TFBS that together define highly coherent sets of genes.