Genome-wide DNA methylation sequencing identifies epigenetic perturbations in the upper airways under long-term exposure to moderate levels of ambient air pollution.

While the link between exposure to high levels of ambient particulate matter (PM) and increased incidences of respiratory and cardiovascular diseases is widely recognized, recent epidemiological studies have shown that low PM concentrations are equally associated with adverse health effects. As DNA methylation is one of the main mechanisms by which cells regulate and stabilize gene expression, changes in the methylome could constitute early indicators of dysregulated signaling pathways. So far, little is known about PM-associated DNA methylation changes in the upper airways, the first point of contact between airborne pollutants and the human body. Here, we focused on cells of the upper respiratory tract and assessed their genome-wide DNA methylation pattern to explore exposure-associated early regulatory changes. Using a mobile epidemiological laboratory, nasal lavage samples were collected from a cohort of 60 adults that lived in districts with records of low (Simmerath) or moderate (Stuttgart) PM10 levels in Germany. PM10 concentrations were verified by particle measurements on the days of the sample collection and genome-wide DNA methylation was determined by enzymatic methyl sequencing at single-base resolution. We identified 231 differentially methylated regions (DMRs) between moderately and lowly PM10 exposed individuals. A high proportion of DMRs overlapped with regulatory elements, and DMR target genes were involved in pathways regulating cellular redox homeostasis and immune response. In addition, we found distinct changes in DNA methylation of the HOXA gene cluster whose methylation levels have previously been linked to air pollution exposure but also to carcinogenesis in several instances. The findings of this study suggest that regulatory changes in upper airway cells occur at PM10 levels below current European thresholds, some of which may be involved in the development of air pollution-related diseases.

Global hypomethylation in childhood asthma identified by genome-wide DNA-methylation sequencing preferentially affects enhancer regions.

BACKGROUND: Childhood asthma is a result of a complex interaction of genetic and environmental components causing epigenetic and immune dysregulation, airway inflammation and impaired lung function. Although different microarray based EWAS studies have been conducted, the impact of epigenetic regulation in asthma development is still widely unknown. We have therefore applied unbiased whole genome bisulfite sequencing (WGBS) to characterize global DNA-methylation profiles of asthmatic children compared to healthy controls. METHODS: Peripheral blood samples of 40 asthmatic and 42 control children aged 5-15 years from three birth cohorts were sequenced together with paired cord blood samples. Identified differentially methylated regions (DMRs) were categorized in genotype-associated, cell-type-dependent, or prenatally primed. Network analysis and subsequent natural language processing of DMR-associated genes was complemented by targeted analysis of functional translation of epigenetic regulation on the transcriptional and protein level. RESULTS: In total, 158 DMRs were identified in asthmatic children compared to controls of which 37% were related to the eosinophil content. A global hypomethylation was identified affecting predominantly enhancer regions and regulating key immune genes such as IL4, IL5RA, and EPX. These DMRs were confirmed in n = 267 samples and could be linked to aberrant gene expression. Out of the 158 DMRs identified in the established phenotype, 56 were perturbed already at birth and linked, at least in part, to prenatal influences such as tobacco smoke exposure or phthalate exposure. CONCLUSION: This is the first epigenetic study based on whole genome sequencing to identify marked dysregulation of enhancer regions as a hallmark of childhood asthma.

Differentially methylated regions within lung cancer risk loci are enriched in deregulated enhancers.

Genome-wide association studies (GWAS) have identified SNPs linked with lung cancer risk. Our aim was to discover the genes, non-coding RNAs, and regulatory elements within GWAS-identified risk regions that are deregulated in non-small cell lung carcinoma (NSCLC) to identify novel, clinically targetable genes and mechanisms in carcinogenesis. A targeted bisulphite-sequencing approach was used to comprehensively investigate DNA methylation changes occurring within lung cancer risk regions in 17 NSCLC and adjacent normal tissue pairs. We report differences in differentially methylated regions between adenocarcinoma and squamous cell carcinoma. Among the minimal regions found to be differentially methylated in at least 50% of the patients, 7 candidates were replicated in 2 independent cohorts (n = 27 and n = 87) and the potential of 6 as methylation-dependent regulatory elements was confirmed by functional assays. This study contributes to understanding the pathways implicated in lung cancer initiation and progression, and provides new potential targets for cancer treatment.

Mutational mechanisms shaping the coding and noncoding genome of germinal center derived B-cell lymphomas.

B cells have the unique property to somatically alter their immunoglobulin (IG) genes by V(D)J recombination, somatic hypermutation (SHM) and class-switch recombination (CSR). Aberrant targeting of these mechanisms is implicated in lymphomagenesis, but the mutational processes are poorly understood. By performing whole genome and transcriptome sequencing of 181 germinal center derived B-cell lymphomas (gcBCL) we identified distinct mutational signatures linked to SHM and CSR. We show that not only SHM, but presumably also CSR causes off-target mutations in non-IG genes. Kataegis clusters with high mutational density mainly affected early replicating regions and were enriched for SHM- and CSR-mediated off-target mutations. Moreover, they often co-occurred in loci physically interacting in the nucleus, suggesting that mutation hotspots promote increased mutation targeting of spatially co-localized loci (termed hypermutation by proxy). Only around 1% of somatic small variants were in protein coding sequences, but in about half of the driver genes, a contribution of B-cell specific mutational processes to their mutations was found. The B-cell-specific mutational processes contribute to both lymphoma initiation and intratumoral heterogeneity. Overall, we demonstrate that mutational processes involved in the development of gcBCL are more complex than previously appreciated, and that B cell-specific mutational processes contribute via diverse mechanisms to lymphomagenesis.

Lipomatous Solitary Fibrous Tumors Harbor Rare NAB2-STAT6 Fusion Variants and Show Up-Regulation of the Gene PPARG, Encoding for a Regulator of Adipocyte Differentiation.

Solitary fibrous tumors (SFTs) harbor activating NAB2-STAT6 gene fusions. Different variants of the NAB2-STAT6 gene fusion have been associated with distinct clinicopathologic features. Lipomatous SFTs are a morphologic variant of SFTs, characterized by a fat-forming tumor component. Our aim was to evaluate NAB2-STAT6 fusion variants and to further study the molecular genetic features in a cohort of lipomatous SFTs. A hybrid-capture-based next-generation sequencing panel was employed to detect NAB2-STAT6 gene fusions at the RNA level. In addition, the RNA expression levels of 507 genes were evaluated using this panel, and were compared with a control cohort of nonlipomatous SFTs. Notably, 5 of 11 (45%) of lipomatous SFTs in the current series harbored the uncommon NAB2 exon 4-STAT6 exon 4 gene fusion variant, which is observed in only 0.9% to 1.4% of nonlipomatous SFTs. Furthermore, lipomatous SFTs displayed significant differences in gene expression compared with their nonlipomatous counterparts, including up-regulation of the gene peroxisome proliferator activated receptor-γ (PPARG). Peroxisome proliferator activated receptor-γ is a nuclear receptor regulating adipocyte differentiation, providing a possible explanation for the fat-forming component in lipomatous SFTs. In summary, the current study provides a possible molecular genetic basis for the distinct morphologic features of lipomatous SFTs.

Gene Expression in Solitary Fibrous Tumors (SFTs) Correlates with Anatomic Localization and NAB2-STAT6 Gene Fusion Variants.

Solitary fibrous tumors (SFTs) harbor recurrent NAB2-STAT6 gene fusions, promoting constitutional up-regulation of oncogenic early growth response 1 (EGR1)-dependent gene expression. SFTs with the most common canonical NAB2 exon 4-STAT6 exon 2 fusion variant are often located in the thorax (pleuropulmonary) and are less cellular with abundant collagen. In contrast, SFTs with NAB2 exon 6-STAT6 exon 16/17 fusion variants typically display a cellular round to ovoid cell morphology and are often located in the deep soft tissue of the retroperitoneum and intra-abdominal pelvic region or in the meninges. Here, we employed next-generation sequencing-based gene expression profiling to identify significant differences in gene expression associated with anatomic localization and NAB2-STAT6 gene fusion variants. SFTs with the NAB2 exon 4-STAT6 exon 2 fusion variant showed a transcriptional signature enriched for genes involved in DNA binding, gene transcription, and nuclear localization, whereas SFTs with the NAB2 exon 6-STAT6 exon 16/17 fusion variants were enriched for genes involved in tyrosine kinase signaling, cell proliferation, and cytoplasmic localization. Specific transcription factor binding motifs were enriched among differentially expressed genes in SFTs with different fusion variants, implicating co-transcription factors in the modification of chimeric NGFI-A binding protein 2 (NAB2)-STAT6-dependent deregulation of EGR1-dependent gene expression. In summary, this study establishes a potential molecular biologic basis for clinicopathologic differences in SFTs with distinct NAB2-STAT6 gene fusion variants.

Unraveling most abundant mutational signatures in head and neck cancer.

Genomic alterations are a driving force in the multistep process of head and neck cancer (HNC) and result from the interaction of exogenous environmental exposures and endogenous cellular processes. Each of these processes leaves a characteristic pattern of mutations on the tumor genome providing the unique opportunity to decipher specific signatures of mutational processes operative during HNC pathogenesis and to address their prognostic value. Computational analysis of whole exome sequencing data of the HIPO-HNC (Heidelberg Center for Personalized Oncology-head and neck cancer) (n = 83) and TCGA-HNSC (The Cancer Genome Atlas-Head and Neck Squamous Cell Carcinoma) (n = 506) cohorts revealed five common mutational signatures (Catalogue of Somatic Mutations in Cancer [COSMIC] Signatures 1, 2, 3, 13 and 16) and demonstrated their significant association with etiological risk factors (tobacco, alcohol and HPV16). Unsupervised hierarchical clustering identified four clusters (A, B, C1 and C2) of which Subcluster C2 was enriched for cases with a higher frequency of signature 16 mutations. Tumors of Subcluster C2 had significantly lower p16INK4A expression accompanied by homozygous CDKN2A deletion in almost one half of cases. Survival analysis revealed an unfavorable prognosis for patients with tumors characterized by a higher mutation burden attributed to signature 16 as well as cases in Subcluster C2. Finally, a LASSO-Cox regression model was applied to prioritize clinically relevant signatures and to establish a prognostic risk score for head and neck squamous cell carcinoma patients. In conclusion, our study provides a proof of concept that computational analysis of somatic mutational signatures is not only a powerful tool to decipher environmental and intrinsic processes in the pathogenesis of HNC, but could also pave the way to establish reliable prognostic patterns.

Genome-Wide DNA Methylation Profiling in Early Stage I Lung Adenocarcinoma Reveals Predictive Aberrant Methylation in the Promoter Region of the Long Noncoding RNA PLUT: An Exploratory Study.

INTRODUCTION: Surgical procedure is the treatment of choice in early stage I lung adenocarcinoma. However, a considerable number of patients experience recurrence within the first 2 years after complete resection. Suitable prognostic biomarkers that identify patients at high risk of recurrence (who may probably benefit from adjuvant treatment) are still not available. This study aimed at identifying methylation markers for early recurrence that may become important tools for the development of new treatment modalities. METHODS: Genome-wide DNA methylation profiling was performed on 30 stage I lung adenocarcinomas, comparing 14 patients with early metastatic recurrence with 16 patients with a long-term relapse-free survival period using methylated-CpG-immunoprecipitation followed by high-throughput next-generation sequencing. The differentially methylated regions between the two subgroups were validated for their prognostic value in two independent cohorts using the MassCLEAVE assay, a high-resolution quantitative methylation analysis. RESULTS: Unsupervised clustering of patients in the discovery cohort on the basis of differentially methylated regions identified patients with shorter relapse-free survival (hazard ratio: 2.23; 95% confidence interval: 0.66-7.53; p = 0.03). In two validation cohorts, promoter hypermethylation of the long noncoding RNA PLUT was significantly associated with shorter relapse-free survival (hazard ratio: 0.54; 95% confidence interval: 0.31-0.93; p < 0.026) and could be reported as an independent prognostic factor in the multivariate Cox regression analysis. CONCLUSIONS: Promoter hypermethylation of the long noncoding RNA PLUT is predictive in patients with early stage I adenocarcinoma at high risk for early recurrence. Further studies are needed to validate its role in carcinogenesis and its use as a biomarker to facilitate patient selection and risk stratification.

Integrative Analysis of Multi-omics Data Identified EGFR and PTGS2 as Key Nodes in a Gene Regulatory Network Related to Immune Phenotypes in Head and Neck Cancer.

PURPOSE: Malignant progression exhibits a tightly orchestrated balance between immune effector response and tolerance. However, underlying molecular principles that drive the establishment and maintenance of the tumor immune phenotype remain to be elucidated. EXPERIMENTAL DESIGN: We trained a novel molecular classifier based on immune cell subsets related to programmed death-ligand 1 (PD-L1) and interferon γ (IFNγ) expression, which revealed distinct subgroups with higher (cluster A) or lower (subcluster B3) cytotoxic immune phenotypes. Integrative analysis of multi-omics data was conducted to identify differences in genetic and epigenetic landscapes as well as their impact on differentially expressed genes (DEG) among immune phenotypes. A prognostic gene signature for immune checkpoint inhibition (ICI) was established by a least absolute shrinkage and selection operator (LASSO)-Cox regression model. RESULTS: Mutational landscape analyses unraveled a higher frequency of CASP8 somatic mutations in subcluster A1, while subcluster B3 exhibited a characteristic pattern of copy-number alterations affecting chemokine signaling and immune effector response. The integrative multi-omics approach identified EGFR and PTGS2 as key nodes in a gene regulatory network related to the immune phenotype, and several DEGs related to the immune phenotypes were affected by EGFR inhibition in tumor cell lines. Finally, we established a prognostic gene signature by a LASSO-Cox regression model based on DEGs between nonprogressive disease and progressive disease subgroups for ICI. CONCLUSIONS: Our data highlight a complex interplay between genetic and epigenetic events in the establishment of the tumor immune phenotype and provide compelling experimental evidence that a patient with squamous cell carcinoma of the head and neck at higher risk for ICI treatment failure might benefit from a combination with EGFR inhibition.

NOTCH target gene HES5 mediates oncogenic and tumor suppressive functions in hepatocarcinogenesis.

NOTCH receptor signaling plays a pivotal role in liver homeostasis and hepatocarcinogenesis. However, the role of NOTCH pathway mutations and the NOTCH target gene HES5 in liver tumorigenesis are poorly understood. Here we performed whole-exome sequencing of 54 human HCC specimens and compared the prevalence of NOTCH pathway component mutations with the TCGA-LIHC cohort (N = 364). In addition, we functionally characterized the NOTCH target HES5 and the patient-derived HES5-R31G mutation in vitro and in an orthotopic mouse model applying different oncogenic backgrounds, to dissect the role of HES5 in different tumor subgroups in vivo. We identified nonsynonymous mutations in 14 immediate NOTCH pathway genes affecting 24.1% and 16.8% of HCC patients in the two independent cohorts, respectively. Among these, the HES5-R31G mutation was predicted in silico to have high biological relevance. Functional analyses in cell culture showed that HES5 reduced cell migration and clonogenicity. Further analyses revealed that the patient-derived HES5-R31G mutant protein was non-functional due to loss of DNA binding and greatly reduced nuclear localization. Furthermore, HES5 exhibited a negative feedback loop by directly inhibiting the NOTCH target HES1 and downregulated the pro-proliferative MYC targets ODC1 and LDHA. Interestingly, HES5 inhibited MYC-dependent hepatocarcinogenesis, whereas it promoted AKT-dependent liver tumor formation and stem cell features in a murine model. Thus, NOTCH pathway component mutations are commonly observed in HCC. Furthermore, the NOTCH target gene HES5 has both pro- and anti-tumorigenic functions in liver cancer proposing a driver gene dependency and it promotes tumorigenesis with its interaction partner AKT.

Whole-genome fingerprint of the DNA methylome during chemically induced differentiation of the human AML cell line HL-60/S4.

Epigenomic regulation plays a vital role in cell differentiation. The leukemic HL-60/S4 [human myeloid leukemic cell line HL-60/S4 (ATCC CRL-3306)] promyelocytic cell can be easily differentiated from its undifferentiated promyelocyte state into neutrophil- and macrophage-like cell states. In this study, we present the underlying genome and epigenome architecture of HL-60/S4 through its differentiation. We performed whole-genome bisulphite sequencing of HL-60/S4 cells and their differentiated counterparts. With the support of karyotyping, we show that HL-60/S4 maintains a stable genome throughout differentiation. Analysis of differential Cytosine-phosphate-Guanine dinucleotide methylation reveals that most methylation changes occur in the macrophage-like state. Differential methylation of promoters was associated with immune-related terms. Key immune genes, CEBPA, GFI1, MAFB and GATA1 showed differential expression and methylation. However, we observed the strongest enrichment of methylation changes in enhancers and CTCF binding sites, implying that methylation plays a major role in large-scale transcriptional reprogramming and chromatin reorganisation during differentiation. Correlation of differential expression and distal methylation with support from chromatin capture experiments allowed us to identify putative proximal and long-range enhancers for a number of immune cell differentiation genes, including CEBPA and CCNF Integrating expression data, we present a model of HL-60/S4 differentiation in relation to the wider scope of myeloid differentiation.

Impact of cancer mutational signatures on transcription factor motifs in the human genome.

BACKGROUND: Somatic mutations in cancer genomes occur through a variety of molecular mechanisms, which contribute to different mutational patterns. To summarize these, mutational signatures have been defined using a large number of cancer genomes, and related to distinct mutagenic processes. Each cancer genome can be compared to this reference dataset and its exposure to one or the other signature be determined. Given the very different mutational patterns of these signatures, we anticipate that they will have distinct impact on genomic elements, in particular motifs for transcription factor binding sites (TFBS). METHODS: We used the 30 mutational signatures from the COSMIC database, and derived a theoretical framework to infer the impact of these signatures on the alteration of transcription factor (TF) binding motifs from the JASPAR database. Hence, we translated the trinucleotide mutation frequencies of the signatures into alteration frequencies of specific TF binding motifs, leading either to creation or disruption of these motifs. RESULTS: Motif families show different susceptibility to alterations induced by the mutational signatures. For certain motifs, a high correlation is observed between the TFBS motif creation and disruption events related to the information content of the motif. Moreover, we observe striking patterns regarding for example the Ets-motif family, for which a high impact of UV induced signatures is observed. Our model also confirms the susceptibility of specific transcription factor motifs to deamination processes. CONCLUSION: Our results show that the mutational signatures have different impact on the binding motifs of transcription factors and that for certain high complexity motifs there is a strong correlation between creation and disruption, related to the information content of the motif. This study represents a background estimation of the alterations due purely to mutational signatures in the absence of additional contributions, e.g. from evolutionary processes.

Nuclear NR4A3 Immunostaining Is a Specific and Sensitive Novel Marker for Acinic Cell Carcinoma of the Salivary Glands.

Recently, we discovered the recurrent genomic rearrangement [t(4;9)(q13;q31)] enabling upregulation of the transcription factor Nuclear Receptor Subfamily 4 Group A Member 3 (NR4A3) through enhancer hijacking as the oncogenic driver event in acinic cell carcinoma (AciCC) of the salivary glands. In the current study, we evaluated the usefulness of NR4A3 immunostaining and NR4A3 fluorescence in situ hybridization (FISH) in the differential diagnosis of AciCC, comparing a total of 64 AciCCs including 17% cases with high-grade transformation, 29 secretory (mammary analog) carcinomas (MASC), and 70 other salivary gland carcinomas. Nuclear NR4A3 immunostaining was a highly specific (100%) and sensitive (98%) marker for AciCC with only 1 negative case, whereas NR4A3 FISH was less sensitive (84%). None of the MASCs or other salivary gland carcinomas displayed any nuclear NR4A3 immunostaining. The recently described HTN3-MSANTD3 gene fusion was observed in 4 of 49 (8%) evaluable AciCCs, all with nuclear NR4A3 immunostaining. In summary, NR4A3 immunostaining is a highly specific and sensitive marker for AciCC, which may be especially valuable in cases with high-grade transformation and in "zymogen granule"-poor examples within the differential diagnostic spectrum of AciCC and MASC.

Somatic mutations and promotor methylation of the ryanodine receptor 2 is a common event in the pathogenesis of head and neck cancer.

Genomic sequencing projects unraveled the mutational landscape of head and neck squamous cell carcinoma (HNSCC) and provided a comprehensive catalog of somatic mutations. However, the limited number of significant cancer-related genes obtained so far only partially explains the biological complexity of HNSCC and hampers the development of novel diagnostic biomarkers and therapeutic targets. We pursued a multiscale omics approach based on whole-exome sequencing, global DNA methylation and gene expression profiling data derived from tumor samples of the HIPO-HNC cohort (n = 87), and confirmed new findings with datasets from The Cancer Genome Atlas (TCGA). Promoter methylation was confirmed by MassARRAY analysis and protein expression was assessed by immunohistochemistry and immunofluorescence staining. We discovered a set of cancer-related genes with frequent somatic mutations and high frequency of promoter methylation. This included the ryanodine receptor 2 (RYR2), which showed variable promoter methylation and expression in both tumor samples and cell lines. Immunohistochemical staining of tissue sections unraveled a gradual loss of RYR2 expression from normal mucosa via dysplastic lesion to invasive cancer and indicated that reduced RYR2 expression in adjacent tissue and precancerous lesions might serve as risk factor for unfavorable prognosis and upcoming malignant conversion. In summary, our data indicate that impaired RYR2 function by either somatic mutation or epigenetic silencing is a common event in HNSCC pathogenesis. Detection of RYR2 expression and/or promoter methylation might enable risk assessment for malignant conversion of dysplastic lesions.

Enhancer hijacking activates oncogenic transcription factor NR4A3 in acinic cell carcinomas of the salivary glands.

The molecular pathogenesis of salivary gland acinic cell carcinoma (AciCC) is poorly understood. The secretory Ca-binding phosphoprotein (SCPP) gene cluster at 4q13 encodes structurally related phosphoproteins of which some are specifically expressed at high levels in the salivary glands and constitute major components of saliva. Here we report on recurrent rearrangements [t(4;9)(q13;q31)] in AciCC that translocate active enhancer regions from the SCPP gene cluster to the region upstream of Nuclear Receptor Subfamily 4 Group A Member 3 (NR4A3) at 9q31. We show that NR4A3 is specifically upregulated in AciCCs, and that active chromatin regions and gene expression signatures in AciCCs are highly correlated with the NR4A3 transcription factor binding motif. Overexpression of NR4A3 in mouse salivary gland cells increases expression of known NR4A3 target genes and has a stimulatory functional effect on cell proliferation. We conclude that NR4A3 is upregulated through enhancer hijacking and has important oncogenic functions in AciCC.

The whole-genome landscape of medulloblastoma subtypes.

Current therapies for medulloblastoma, a highly malignant childhood brain tumour, impose debilitating effects on the developing child, and highlight the need for molecularly targeted treatments with reduced toxicity. Previous studies have been unable to identify the full spectrum of driver genes and molecular processes that operate in medulloblastoma subgroups. Here we analyse the somatic landscape across 491 sequenced medulloblastoma samples and the molecular heterogeneity among 1,256 epigenetically analysed cases, and identify subgroup-specific driver alterations that include previously undiscovered actionable targets. Driver mutations were confidently assigned to most patients belonging to Group 3 and Group 4 medulloblastoma subgroups, greatly enhancing previous knowledge. New molecular subtypes were differentially enriched for specific driver events, including hotspot in-frame insertions that target KBTBD4 and 'enhancer hijacking' events that activate PRDM6. Thus, the application of integrative genomics to an extensive cohort of clinical samples derived from a single childhood cancer entity revealed a series of cancer genes and biologically relevant subtype diversity that represent attractive therapeutic targets for the treatment of patients with medulloblastoma.

Recurrent Somatic PDGFRB Mutations in Sporadic Infantile/Solitary Adult Myofibromas But Not in Angioleiomyomas and Myopericytomas.

Infantile myofibroma (MF) is an uncommon benign myofibroblastic tumor of infancy and childhood. Solitary adult MF shares similar features with infantile MF. The lesions occur in 3 clinicopathologic settings: solitary, multicentric, and generalized and can be either sporadic or familial. Traditionally, infantile MF has been included in the spectrum of infantile hemangiopericytoma. The recent World Health Organization classification listed MF, angioleiomyoma, and myopericytoma under the general heading of perivascular tumors in the sense of a morphologic spectrum of perivascular myoid cell neoplasms. Although activating germline PDGFRB mutations have recently been linked to familial infantile MF, the molecular pathogenesis of sporadic infantile and adult solitary MF remained unclear. In this study, we analyzed 25 solitary MFs without evidence of familial disease (9 infantile and 16 adult MFs) to address the question whether somatic PDGFRB mutations might be responsible for the sporadic form of the disease. Given the presumed histogenetic link of MF to myopericytoma and angioleiomyoma, we additionally analyzed a control group of 6 myopericytomas and 9 angioleiomyomas for PDGFRB mutations. We detected PDGFRB mutations in 6/8 (75%) analyzable infantile and in 11/16 (69%) adult MFs but in none of the angioleiomyomas or myopericytomas. In 2 infantile MFs, additional sequencing of the germline confirmed the somatic nature of PDGFRB mutations. To our knowledge, this is the first study reporting apparently somatic recurrent PDGFRB mutations as molecular driver events in the majority of sporadic infantile and adult solitary MFs. Our results suggest molecular distinctness of MF as compared with angioleiomyoma/myopericytoma. Investigation of more cases including those with atypical and worrisome features, as well as other mimickers in the heterogenous morphologic spectrum of MF, is mandatory for validating the potential diagnostic value of PDGFRB mutation testing as a possible surrogate in difficult-to-classify lesions.

Prenatal maternal stress and wheeze in children: novel insights into epigenetic regulation.

Psychological stress during pregnancy increases the risk of childhood wheeze and asthma. However, the transmitting mechanisms remain largely unknown. Since epigenetic alterations have emerged as a link between perturbations in the prenatal environment and an increased disease risk we used whole genome bisulfite sequencing (WGBS) to analyze changes in DNA methylation in mothers and their children related to prenatal psychosocial stress and assessed its role in the development of wheeze in the child. We evaluated genomic regions altered in their methylation level due to maternal stress based of WGBS data of 10 mother-child-pairs. These data were complemented by longitudinal targeted methylation and transcriptional analyses in children from our prospective mother-child cohort LINA for whom maternal stress and wheezing information was available (n = 443). High maternal stress was associated with an increased risk for persistent wheezing in the child until the age of 5. Both mothers and children showed genome-wide alterations in DNA-methylation specifically in enhancer elements. Deregulated neuroendocrine and neurotransmitter receptor interactions were observed in stressed mothers and their children. In children but not in mothers, calcium- and Wnt-signaling required for lung maturation in the prenatal period were epigenetically deregulated and could be linked with wheezing later in children's life.