Context-dependent DNA methylation signatures in animal livestock.

DNA methylation is an important epigenetic modification that is widely conserved across animal genomes. It is widely accepted that DNA methylation patterns can change in a context-dependent manner, including in response to changing environmental parameters. However, this phenomenon has not been analyzed in animal livestock yet, where it holds major potential for biomarker development. Building on the previous identification of population-specific DNA methylation in clonal marbled crayfish, we have now generated numerous base-resolution methylomes to analyze location-specific DNA methylation patterns. We also describe the time-dependent conversion of epigenetic signatures upon transfer from one environment to another. We further demonstrate production system-specific methylation signatures in shrimp, river-specific signatures in salmon and farm-specific signatures in chicken. Together, our findings provide a detailed resource for epigenetic variation in animal livestock and suggest the possibility for origin tracing of animal products by epigenetic fingerprinting.

Morphologic and genetic heterogeneity in breast fibroepithelial lesions-a comprehensive mapping study.

Breast fibroepithelial lesions (FELs) encompass the common fibroadenoma (FA) and relatively rare phyllodes tumour (PT); the latter entity is usually classified as benign, borderline or malignant. Intratumoural heterogeneity is frequently present in these tumours, making accurate histologic evaluation challenging. Despite their rarity, PTs are an important clinical problem due to their propensity for recurrence and, in the case of malignant PT, metastasis. Surgical excision is the mainstay of management. Recent work has uncovered myriad genetic alterations in breast FELs. In this study, exome sequencing was performed on seven cases of morphologically heterogeneous breast FELs, including FAs, PTs of all grades, and a case of metaplastic spindle cell carcinoma arising in PT, in order to elucidate their intratumoural genetic repertoire. Gene mutations identified encompassed cell signalling, tumour suppressor, DNA repair and cell cycle regulating pathways. Mutations common to multiple tumour regions generally showed higher variant allele frequency. Frequent mutations included MED12, TP53, RARA and PIK3CA. Histological observations of increased cellular density and pleomorphism correlated with mutational burden. Phylogenetic analyses revealed disparate pathways of possible tumour progression. In summary, histological heterogeneity correlated with genetic changes in breast FELs.

A novel genomic panel as an adjunctive diagnostic tool for the characterization and profiling of breast Fibroepithelial lesions.

BACKGROUND: Known collectively as breast fibroepithelial lesions (FELs), the common fibroadenomas (FAs) and the rarer phyllodes tumors (PTs) are a heterogenous group of biphasic neoplasms. Owing to limited tissue availability, inter-observer variability, overlapping histological features and heterogeneity of these lesions, diagnosing them accurately on core biopsies is challenging. As the choice management option depends on the histological diagnosis; a novel 16-gene panel assay was developed to improve the accuracy of preoperative diagnosis on core biopsy specimens. METHODS: Using this 16-gene panel, targeted amplicon-based sequencing was performed on 275 formalin-fixed, paraffin-embedded (FFPE) breast FEL specimens, archived at the Singapore General Hospital, from 2008 to 2012. RESULTS: In total, 167 FAs, 24 benign, 14 borderline and 6 malignant PTs, were profiled. Compared to FAs, PTs had significantly higher mutation rates in the TERT promoter (p <  0.001), RARA (p <  0.001), FLNA, RB1 and TP53 (p = 0.002, 0.020 and 0.018, respectively). In addition to a higher mutational count (p <  0.001), TERT promoter (p <  0.001), frameshift, nonsense and splice site (p = 0.001, < 0.001 and 0.043, respectively) mutations were also frequently observed in PTs. A multivariate logistic regression model was built using these as variables and a predictive scoring system was developed. It classifies a FEL at low or high risk (score <  1 and ≥ 1, respectively) of being a PT. This scoring system has good discrimination (ROC area = 0.773, 95% CI: 0.70 to 0.85), calibration (p = 0.945) and is significant in predicting PTs (p <  0.001). CONCLUSION: This novel study demonstrates the ability to extract DNA of sufficient quality and quantity for targeted sequencing from FFPE breast core biopsy specimens, along with their successful characterization and profiling using our customized 16-gene panel. Prospective work includes validating the utility of this promising 16-gene panel assay as an adjunctive diagnostic tool in clinical practice.

Oncogenic activation of the STAT3 pathway drives PD-L1 expression in natural killer/T-cell lymphoma.

Mature T-cell lymphomas, including peripheral T-cell lymphoma (PTCL) and extranodal NK/T-cell lymphoma (NKTL), represent a heterogeneous group of non-Hodgkin lymphomas with dismal outcomes and limited treatment options. To determine the extent of involvement of the JAK/STAT pathway in this malignancy, we performed targeted capture sequencing of 188 genes in this pathway in 171 PTCL and NKTL cases. A total of 272 nonsynonymous somatic mutations in 101 genes were identified in 73% of the samples, including 258 single-nucleotide variants and 14 insertions or deletions. Recurrent mutations were most frequently located in STAT3 and TP53 (15%), followed by JAK3 and JAK1 (6%) and SOCS1 (4%). A high prevalence of STAT3 mutation (21%) was observed specifically in NKTL. Novel STAT3 mutations (p.D427H, E616G, p.E616K, and p.E696K) were shown to increase STAT3 phosphorylation and transcriptional activity of STAT3 in the absence of cytokine, in which p.E616K induced programmed cell death-ligand 1 (PD-L1) expression by robust binding of activated STAT3 to the PD-L1 gene promoter. Consistent with these findings, PD-L1 was overexpressed in NKTL cell lines harboring hotspot STAT3 mutations, and similar findings were observed by the overexpression of p.E616K and p.E616G in the STAT3 wild-type NKTL cell line. Conversely, STAT3 silencing and inhibition decreased PD-L1 expression in STAT3 mutant NKTL cell lines. In NKTL tumors, STAT3 activation correlated significantly with PD-L1 expression. We demonstrated that STAT3 activation confers high PD-L1 expression, which may promote tumor immune evasion. The combination of PD-1/PD-L1 antibodies and STAT3 inhibitors might be a promising therapeutic approach for NKTL, and possibly PTCL.

Whole-Genome and Epigenomic Landscapes of Etiologically Distinct Subtypes of Cholangiocarcinoma.

Cholangiocarcinoma (CCA) is a hepatobiliary malignancy exhibiting high incidence in countries with endemic liver-fluke infection. We analyzed 489 CCAs from 10 countries, combining whole-genome (71 cases), targeted/exome, copy-number, gene expression, and DNA methylation information. Integrative clustering defined 4 CCA clusters-fluke-positive CCAs (clusters 1/2) are enriched in ERBB2 amplifications and TP53 mutations; conversely, fluke-negative CCAs (clusters 3/4) exhibit high copy-number alterations and PD-1/PD-L2 expression, or epigenetic mutations (IDH1/2, BAP1) and FGFR/PRKA-related gene rearrangements. Whole-genome analysis highlighted FGFR2 3' untranslated region deletion as a mechanism of FGFR2 upregulation. Integration of noncoding promoter mutations with protein-DNA binding profiles demonstrates pervasive modulation of H3K27me3-associated sites in CCA. Clusters 1 and 4 exhibit distinct DNA hypermethylation patterns targeting either CpG islands or shores-mutation signature and subclonality analysis suggests that these reflect different mutational pathways. Our results exemplify how genetics, epigenetics, and environmental carcinogens can interplay across different geographies to generate distinct molecular subtypes of cancer.Significance: Integrated whole-genome and epigenomic analysis of CCA on an international scale identifies new CCA driver genes, noncoding promoter mutations, and structural variants. CCA molecular landscapes differ radically by etiology, underscoring how distinct cancer subtypes in the same organ may arise through different extrinsic and intrinsic carcinogenic processes. Cancer Discov; 7(10); 1116-35. ©2017 AACR.This article is highlighted in the In This Issue feature, p. 1047.

Loss of tumor suppressor KDM6A amplifies PRC2-regulated transcriptional repression in bladder cancer and can be targeted through inhibition of EZH2.

Trithorax-like group complex containing KDM6A acts antagonistically to Polycomb-repressive complex 2 (PRC2) containing EZH2 in maintaining the dynamics of the repression and activation of gene expression through H3K27 methylation. In urothelial bladder carcinoma, KDM6A (a H3K27 demethylase) is frequently mutated, but its functional consequences and therapeutic targetability remain unknown. About 70% of KDM6A mutations resulted in a total loss of expression and a consequent loss of demethylase function in this cancer type. Further transcriptome analysis found multiple deregulated pathways, especially PRC2/EZH2, in KDM6A-mutated urothelial bladder carcinoma. Chromatin immunoprecipitation sequencing analysis revealed enrichment of H3K27me3 at specific loci in KDM6A-null cells, including PRC2/EZH2 and their downstream targets. Consequently, we targeted EZH2 (an H3K27 methylase) and demonstrated that KDM6A-null urothelial bladder carcinoma cell lines were sensitive to EZH2 inhibition. Loss- and gain-of-function assays confirmed that cells with loss of KDM6A are vulnerable to EZH2. IGFBP3, a direct KDM6A/EZH2/H3K27me3 target, was up-regulated by EZH2 inhibition and contributed to the observed EZH2-dependent growth suppression in KDM6A-null cell lines. EZH2 inhibition delayed tumor onset in KDM6A-null cells and caused regression of KDM6A-null bladder tumors in both patient-derived and cell line xenograft models. In summary, our study demonstrates that inactivating mutations of KDM6A, which are common in urothelial bladder carcinoma, are potentially targetable by inhibiting EZH2.

MED12 protein expression in breast fibroepithelial lesions: correlation with mutation status and oestrogen receptor expression.

AIMS: Recent reports have identified recurrent MED12 somatic mutations in fibroadenomas and phyllodes tumours. The frequency and type of somatic mutations were noted to be similar to those of uterine leiomyomas. We aimed to investigate protein expression of MED12, correlating it to MED12 mutational status and expression of oestrogen receptors (ER). METHODS: Immunohistochemistry was performed on a total of 232 fibroepithelial lesions (100 fibroadenomas, 132 phyllodes tumours) diagnosed at the Department of Pathology, Singapore General Hospital using MED12, ERα and ERß antibodies. Expressions were evaluated in both stroma and epithelium, and correlated with MED12 mutational status. RESULTS: MED12 mutation was significantly associated with high MED12 protein expression (H-score >150) in the stroma (p=0.029), but not in the epithelium. It was not associated with ERα and ERß protein expression in both stroma and epithelium. MED12 protein expression was significantly correlated with ERα in epithelial (p=0.007) and ERß in stromal (p=0.049) components. MED12 was not significantly different between fibroadenomas and phyllodes tumours. Epithelial expression of ERα was significantly higher in fibroadenomas (p<0.001) than in phyllodes tumours. Conversely, both epithelial and stromal expression of ERß was significantly higher in phyllodes tumours (p<0.001). CONCLUSIONS: Positive associations observed between MED12 and ERα, ERß immunohistochemical expression suggest a biological interplay between the proteins. The lack of significant association of MED12 mutation with ER protein expression indicates a need to further explore the functional impact of MED12 mutations on the ER signalling pathway in breast fibroepithelial lesions.

Genomic landscapes of breast fibroepithelial tumors.

Breast fibroepithelial tumors comprise a heterogeneous spectrum of pathological entities, from benign fibroadenomas to malignant phyllodes tumors. Although MED12 mutations have been frequently found in fibroadenomas and phyllodes tumors, the landscapes of genetic alterations across the fibroepithelial tumor spectrum remain unclear. Here, by performing exome sequencing of 22 phyllodes tumors followed by targeted sequencing of 100 breast fibroepithelial tumors, we observed three distinct somatic mutation patterns. First, we frequently observed MED12 and RARA mutations in both fibroadenomas and phyllodes tumors, emphasizing the importance of these mutations in fibroepithelial tumorigenesis. Second, phyllodes tumors exhibited mutations in FLNA, SETD2 and KMT2D, suggesting a role in driving phyllodes tumor development. Third, borderline and malignant phyllodes tumors harbored additional mutations in cancer-associated genes. RARA mutations exhibited clustering in the portion of the gene encoding the ligand-binding domain, functionally suppressed RARA-mediated transcriptional activation and enhanced RARA interactions with transcriptional co-repressors. This study provides insights into the molecular pathogenesis of breast fibroepithelial tumors, with potential clinical implications.

Exome sequencing identifies highly recurrent MED12 somatic mutations in breast fibroadenoma.

Fibroadenomas are the most common breast tumors in women under 30 (refs. 1,2). Exome sequencing of eight fibroadenomas with matching whole-blood samples revealed recurrent somatic mutations solely in MED12, which encodes a Mediator complex subunit. Targeted sequencing of an additional 90 fibroadenomas confirmed highly frequent MED12 exon 2 mutations (58/98, 59%) that are probably somatic, with 71% of mutations occurring in codon 44. Using laser capture microdissection, we show that MED12 fibroadenoma mutations are present in stromal but not epithelial mammary cells. Expression profiling of MED12-mutated and wild-type fibroadenomas revealed that MED12 mutations are associated with dysregulated estrogen signaling and extracellular matrix organization. The fibroadenoma MED12 mutation spectrum is nearly identical to that of previously reported MED12 lesions in uterine leiomyoma but not those of other tumors. Benign tumors of the breast and uterus, both of which are key target tissues of estrogen, may thus share a common genetic basis underpinned by highly frequent and specific MED12 mutations.

Exome sequencing identifies distinct mutational patterns in liver fluke-related and non-infection-related bile duct cancers.

The impact of different carcinogenic exposures on the specific patterns of somatic mutation in human tumors remains unclear. To address this issue, we profiled 209 cholangiocarcinomas (CCAs) from Asia and Europe, including 108 cases caused by infection with the liver fluke Opisthorchis viverrini and 101 cases caused by non-O. viverrini-related etiologies. Whole-exome sequencing (n = 15) and prevalence screening (n = 194) identified recurrent somatic mutations in BAP1 and ARID1A, neither of which, to our knowledge, has previously been reported to be mutated in CCA. Comparisons between intrahepatic O. viverrini-related and non-O. viverrini-related CCAs demonstrated statistically significant differences in mutation patterns: BAP1, IDH1 and IDH2 were more frequently mutated in non-O. viverrini CCAs, whereas TP53 mutations showed the reciprocal pattern. Functional studies demonstrated tumor suppressive functions for BAP1 and ARID1A, establishing the role of chromatin modulators in CCA pathogenesis. These findings indicate that different causative etiologies may induce distinct somatic alterations, even within the same tumor type.

Regulation of the boundaries of accessible chromatin.

Regulatory regions maintain nucleosome-depleted, open chromatin status but simultaneously require the presence of nucleosomes for specific histone modifications. It remains unclear how these can be achieved for proper regulatory function. Here we demonstrate that nucleosomes positioned within accessible chromatin regions near the boundaries provide platforms for histone modifications while preventing the occlusion of regulatory elements. These boundary nucleosomes were particularly enriched for active or poised regulatory marks in human, such as histone acetylations, H3K4 methylations, H3K9me3, H3K79me2, and H4K20me1. Additionally, we found that based on a genome-wide profiling of ~100 recombinant yeast strains, the location of open chromatin borders tends to vary mostly within 150 bp upon genetic perturbation whereas this positional variation increases in proportion to the sequence preferences of the underlying DNA for nucleosome formation. More than 40% of the local boundary shifts were associated with genetic variation in cis- or trans-acting factors. A sizeable fraction of the identified genetic factors was also associated with nearby gene expression, which was correlated with the distance between the transcription start site (tss) and the boundary that faces the tss. Taken together, the variation in the width of accessible chromatin regions may arise in conjunction with the modulation of the boundary nucleosomes by post-translational modifications or by chromatin regulators and in association with the activity of nearby gene transcription.