Inhibiting DNA Methylation Causes an Interferon Response in Cancer via dsRNA Including Endogenous Retroviruses.

We show that DNA methyltransferase inhibitors (DNMTis) upregulate immune signaling in cancer through the viral defense pathway. In ovarian cancer (OC), DNMTis trigger cytosolic sensing of double-stranded RNA (dsRNA) causing a type I interferon response and apoptosis. Knocking down dsRNA sensors TLR3 and MAVS reduces this response 2-fold and blocking interferon beta or its receptor abrogates it. Upregulation of hypermethylated endogenous retrovirus (ERV) genes accompanies the response and ERV overexpression activates the response. Basal levels of ERV and viral defense gene expression significantly correlate in primary OC and the latter signature separates primary samples for multiple tumor types from The Cancer Genome Atlas into low versus high expression groups. In melanoma patients treated with an immune checkpoint therapy, high viral defense signature expression in tumors significantly associates with durable clinical response and DNMTi treatment sensitizes to anti-CTLA4 therapy in a pre-clinical melanoma model.

Young age at diagnosis is associated with worse prognosis in the Luminal A breast cancer subtype: a retrospective institutional cohort study.

PURPOSE: Although age is a recognized independent prognostic risk factor, its relative importance among molecular subtypes of Breast cancer (BCA) is not well documented. The aim of this study was to evaluate the prognostic role of age at diagnosis among different immunohistochemical subtypes of BCA. METHODS: We conducted a retrospective study of women with invasive BCA undergoing surgery at the Johns Hopkins Hospital, excluding patients presenting with stage IV breast cancer. Patients were stratified into three age groups: ≤ 40, 41-60, and > 60 years, and multivariable analysis was performed using Cox regression. We also identified differentially expressed genes (DEG) between age groups among BCA subtypes in the public TCGA dataset. Finally, we identified key driver genes within the DEGs using a weighted gene co-expression network analysis. RESULTS: Luminal A breast cancer patients had significantly lower 5 year disease-free survival (DFS) and distant metastasis-free survival (DMFS) in the ≤ 40 year age group compared to the 41-60 year age group, while the other molecular subtypes showed no significant association of DFS or DMFS with age. Age was a stronger outcome predictor than tumor grade or proliferative index in Luminal A BCA patients, but not other subtypes. BCA TCGA gene expression data were divided into two groups (≤ 40 years, > 40 years). We identified 374 DEGs in the Luminal A BCA subset, which were enriched in seven pathways and two modules of co-expressed genes. No age group-specific DEGs were identified in non-Luminal A subtypes. CONCLUSIONS: Age at diagnosis may be an important prognostic factor in Luminal A BCA.

Fecal Metabolome in Hmga1 Transgenic Mice with Polyposis: Evidence for Potential Screen for Early Detection of Precursor Lesions in Colorectal Cancer.

Because colorectal cancer (CRC) remains a leading cause of cancer mortality worldwide, more accessible screening tests are urgently needed to identify early stage lesions. We hypothesized that highly sensitive, metabolic profile analysis of stool samples will identify metabolites associated with early stage lesions and could serve as a noninvasive screening test. We therefore applied traveling wave ion mobility mass spectrometry (TWIMMS) coupled with ultraperformance liquid chromatography (UPLC) to investigate metabolic aberrations in stool samples in a transgenic model of premalignant polyposis aberrantly expressing the gene encoding the high mobility group A (Hmga1) chromatin remodeling protein. Here, we report for the first time that the fecal metabolome of Hmga1 mice is distinct from that of control mice and includes metabolites previously identified in human CRC. Significant alterations were observed in fatty acid metabolites and metabolites associated with bile acids (hypoxanthine xanthine, taurine) in Hmga1 mice compared to controls. Surprisingly, a marked increase in the levels of distinctive short, arginine-enriched, tetra-peptide fragments was observed in the transgenic mice. Together these findings suggest that specific metabolites are associated with Hmga1-induced polyposis and abnormal proliferation in intestinal epithelium. Although further studies are needed, these data provide a compelling rationale to develop fecal metabolomic analysis as a noninvasive screening tool to detect early precursor lesions to CRC in humans.

Exploring massive, genome scale datasets with the GenometriCorr package.

UNLABELLED: We have created a statistically grounded tool for determining the correlation of genomewide data with other datasets or known biological features, intended to guide biological exploration of high-dimensional datasets, rather than providing immediate answers. The software enables several biologically motivated approaches to these data and here we describe the rationale and implementation for each approach. Our models and statistics are implemented in an R package that efficiently calculates the spatial correlation between two sets of genomic intervals (data and/or annotated features), for use as a metric of functional interaction. The software handles any type of pointwise or interval data and instead of running analyses with predefined metrics, it computes the significance and direction of several types of spatial association; this is intended to suggest potentially relevant relationships between the datasets. AVAILABILITY AND IMPLEMENTATION: The package, GenometriCorr, can be freely downloaded at http://genometricorr.sourceforge.net/. Installation guidelines and examples are available from the sourceforge repository. The package is pending submission to Bioconductor.

HMGA1 induces FGF19 to drive pancreatic carcinogenesis and stroma formation.

High mobility group A1 (HMGA1) chromatin regulators are upregulated in diverse tumors where they portend adverse outcomes, although how they function in cancer remains unclear. Pancreatic ductal adenocarcinomas (PDACs) are highly lethal tumors characterized by dense desmoplastic stroma composed predominantly of cancer-associated fibroblasts and fibrotic tissue. Here, we uncover an epigenetic program whereby HMGA1 upregulates FGF19 during tumor progression and stroma formation. HMGA1 deficiency disrupts oncogenic properties in vitro while impairing tumor inception and progression in KPC mice and subcutaneous or orthotopic models of PDAC. RNA sequencing revealed HMGA1 transcriptional networks governing proliferation and tumor-stroma interactions, including the FGF19 gene. HMGA1 directly induces FGF19 expression and increases its protein secretion by recruiting active histone marks (H3K4me3, H3K27Ac). Surprisingly, disrupting FGF19 via gene silencing or the FGFR4 inhibitor BLU9931 recapitulates most phenotypes observed with HMGA1 deficiency, decreasing tumor growth and formation of a desmoplastic stroma in mouse models of PDAC. In human PDAC, overexpression of HMGA1 and FGF19 defines a subset of tumors with extremely poor outcomes. Our results reveal what we believe is a new paradigm whereby HMGA1 and FGF19 drive tumor progression and stroma formation, thus illuminating FGF19 as a rational therapeutic target for a molecularly defined PDAC subtype.

Development of an automated liquid biopsy assay for methylated markers in advanced breast cancer.

Current molecular liquid biopsy assays to detect recurrence or monitor response to treatment require sophisticated technology, highly trained personnel, and a turnaround time of weeks. We describe the development and technical validation of an automated Liquid Biopsy for Breast Cancer Methylation (LBx-BCM) prototype, a DNA methylation detection cartridge assay that is simple to perform and quantitatively detects nine methylated markers within 4.5 h. LBx-BCM demonstrated high interassay reproducibility when analyzing exogenous methylated DNA (75-300 DNA copies) spiked into plasma (Coefficient of Variation, CV = 7.1 - 10.9%) and serum (CV = 19.1 - 36.1%). It also demonstrated high interuser reproducibility (Spearman r = 0.887, P < 0.0001) when samples of metastatic breast cancer (MBC, N = 11) and normal control (N = 4) were evaluated independently by two users. Analyses of interplatform reproducibility indicated very high concordance between LBx-BCM and the reference assay, cMethDNA, among 66 paired plasma samples (MBC N = 40, controls N = 26; Spearman r = 0.891; 95% CI = 0.825 - 0.933, P< 0.0001). LBx-BCM achieved a ROC AUC = 0.909 (95% CI = 0.836 - 0.982), 83% sensitivity and 92% specificity; cMethDNA achieved a ROC AUC = 0.896 (95% CI = 0.817 - 0.974), 83% sensitivity and 92% specificity in test set samples. The automated LBx-BCM cartridge prototype is fast, with performance levels equivalent to the highly sensitive, manual cMethDNA method. Future prospective clinical studies will evaluate LBx-BCM detection sensitivity and its ability to monitor therapeutic response during treatment for advanced breast cancer.

Modeling liquid association.

In 2002, Ker-Chau Li introduced the liquid association measure to characterize three-way interactions between genes, and developed a computationally efficient estimator that can be used to screen gene expression microarray data for such interactions. That study, and others published since then, have established the biological validity of the method, and clearly demonstrated it to be a useful tool for the analysis of genomic data sets. To build on this work, we have sought a parametric family of multivariate distributions with the flexibility to model the full range of trivariate dependencies encompassed by liquid association. Such a model could situate liquid association within a formal inferential theory. In this article, we describe such a family of distributions, a trivariate, conditional normal model having Gaussian univariate marginal distributions, and in fact including the trivariate Gaussian family as a special case. Perhaps the most interesting feature of the distribution is that the parameterization naturally parses the three-way dependence structure into a number of distinct, interpretable components. One of these components is very closely aligned to liquid association, and is developed as a measure we call modified liquid association. We develop two methods for estimating this quantity, and propose statistical tests for the existence of this type of dependence. We evaluate these inferential methods in a set of simulations and illustrate their use in the analysis of publicly available experimental data.

Retrospective analysis of cancer-specific gene expression panel for thyroid fine needle aspiration specimens.

BACKGROUND: While molecular testing is a promising strategy for preoperative assessment of cytologically indeterminate thyroid nodules, thyroid fine needle aspiration biopsy (FNA) presents unique challenges for molecular assays, including contaminating peripheral blood mononuclear cells (PBMC) and variable numbers of evaluable epithelial thyroid cells. Moreover, the newly recognized entity, noninvasive follicular thyroid neoplasm with papillary-like nuclear features (NIFTP), has added an additional challenge to the currently available molecular diagnostic platforms. New diagnostic tools are still needed to correctly distinguish benign and malignant thyroid nodules preoperatively. METHODS: Twenty-two transcript splice variants from 12 genes we previously identified as discriminating benign from malignant thyroid nodules were characterized in 80 frozen thyroid tumors from 8 histological subtypes. Isoforms detectable in PBMC were excluded, and the 5 most discriminating isoforms were further validated by real-time quantitative PCR (qPCR) on intraoperative FNA samples from 59 malignant tumors, 55 benign nodules, and 23 NIFTP samples. The qPCR threshold cycle values for each transcript were normalized to the thyrocyte-specific thyroid peroxidase isoform 1 (TPO1) and z-transformed. Receiver operating characteristic (ROC) analyses of the composite transcript scores were used to evaluate classification of thyroid FNAs by the 5-gene isoform expression panel. RESULTS: A molecular signature was developed by combining expression levels of specific isoforms of CDH3, FNDC4, HMGA2, KLK7, and PLAG1. FNAs containing at least 12-36 thyrocytes were sufficient for this assay. The 5-gene composite score achieved an area under the ROC curve (AUC) of 0.86 for distinguishing malignant from benign nodules, with a specificity of 91%, sensitivity of 75%, negative predictive value of 91%, and positive predictive value of 74%. CONCLUSION: Our newly developed 5-gene isoform expression panel distinguishes benign from malignant thyroid tumors and, may help distinguish benign from malignant thyroid nodules in the context of the new NIFTP subtype.

Methylated markers accurately distinguish primary central nervous system lymphomas (PCNSL) from other CNS tumors.

BACKGROUND: Definitive diagnosis of primary central nervous system lymphoma (PCNSL) requires invasive surgical brain biopsy, causing treatment delays. In this paper, we identified and validated tumor-specific markers that can distinguish PCNSL from other CNS tumors in tissues. In a pilot study, we tested these newly identified markers in plasma. RESULTS: The Methylation Outlier Detector program was used to identify markers in TCGA dataset of 48 diffuse large B-cell lymphoma (DLBCL) and 656 glioblastomas and lower-grade gliomas. Eight methylated markers clearly distinguished DLBCL from gliomas. Marker performance was verified (ROC-AUC of ≥ 0.989) in samples from several GEO datasets (95 PCNSL; 2112 other primary CNS tumors of 11 types). Next, we developed a novel, efficient assay called Tailed Amplicon Multiplexed-Methylation-Specific PCR (TAM-MSP), which uses two of the methylation markers, cg0504 and SCG3 triplexed with ACTB. FFPE tissue sections (25 cases each) of PCNSL and eight types of other primary CNS tumors were analyzed using TAM-MSP. TAM-MSP distinguished PCNSL from the other primary CNS tumors with 100% accuracy (AUC = 1.00, 95% CI 0.95-1.00, P < 0.001). The TAM-MSP assay also detected as few as 5 copies of fully methylated plasma DNA spiked into 0.5 ml of healthy plasma. In a pilot study of plasma from 15 PCNSL, 5 other CNS tumors and 6 healthy individuals, methylation in cg0504 and SCG3 was detectable in 3/15 PCNSL samples (20%). CONCLUSION: The Methylation Outlier Detector program identified methylated markers that distinguish PCNSL from other CNS tumors with accuracy. The high level of accuracy achieved by these markers was validated in tissues by a novel method, TAM-MSP. These studies lay a strong foundation for a liquid biopsy-based test to detect PCNSL-specific circulating tumor DNA.

Automated and rapid detection of cancer in suspicious axillary lymph nodes in patients with breast cancer.

Preoperative staging of suspicious axillary lymph nodes (ALNs) allows patients to be triaged to ALN dissection or to sentinel lymph node biopsy (SLNB). Ultrasound-guided fine needle aspiration (FNA) and cytology of ALN is moderately sensitive but its clinical utility relies heavily on the cytologist's experience. We proposed that the 5-h automated GeneXpert system-based prototype breast cancer detection assay (BCDA) that quantitatively measures DNA methylation in ten tumor-specific gene markers could provide a facile, accurate test for detecting cancer in FNA of enlarged lymph nodes. We validated the assay in ALN-FNA samples from a prospective study of patients (N = 230) undergoing SLNB. In a blinded analysis of 218 evaluable LN-FNAs from 108 malignant and 110 benign LNs by histology, BCDA displayed a sensitivity of 90.7% and specificity of 99.1%, achieving an area under the ROC curve, AUC of 0.958 (95% CI: 0.928-0.989; P < 0.0001). Next, we conducted a study of archival FNAs of ipsilateral palpable LNs (malignant, N = 72, benign, N = 53 by cytology) collected in the outpatient setting prior to neoadjuvant chemotherapy (NAC). Using the ROC-threshold determined in the prospective study, compared to cytology, BCDA achieved a sensitivity of 94.4% and a specificity of 92.5% with a ROC-AUC = 0.977 (95% CI: 0.953-1.000; P < 0.0001). Our study shows that the automated assay detects cancer in suspicious lymph nodes with a high level of accuracy within 5 h. This cancer detection assay, scalable for analysis to scores of LN FNAs, could assist in determining eligibility of patients to different treatment regimens.

Measuring DNA Copy Number Variation Using High-Density Methylation Microarrays.

Genetic and epigenetic changes drive carcinogenesis, and their integrated analysis provides insights into mechanisms of cancer development. Computational methods have been developed to measure copy number variation (CNV) from methylation array data, including ChAMP-CNV, CN450K, and, introduced here, Epicopy. Using paired single nucleotide polymorphism (SNP) and methylation array data from the public The Cancer Genome Atlas repository, we optimized CNV calling and benchmarked the performance of these methods. We optimized the thresholds of all three methods and showed comparable performance across methods. Using Epicopy as a representative analysis of Illumina450K array, we show that Illumina450K-derived CNV methods achieve a sensitivity of 0.7 and a positive predictive value of 0.75 in identifying CNVs, which is similar to results achieved when comparing competing SNP microarray platforms with each other.

DNA Methylation Markers for Breast Cancer Detection in the Developing World.

PURPOSE: An unmet need in low-resource countries is an automated breast cancer detection assay to prioritize women who should undergo core breast biopsy and pathologic review. Therefore, we sought to identify and validate a panel of methylated DNA markers to discriminate between cancer and benign breast lesions using cells obtained by fine-needle aspiration (FNA).Experimental Design: Two case-control studies were conducted comparing cancer and benign breast tissue identified from clinical repositories in the United States, China, and South Africa for marker selection/training (N = 226) and testing (N = 246). Twenty-five methylated markers were assayed by Quantitative Multiplex-Methylation-Specific PCR (QM-MSP) to select and test a cancer-specific panel. Next, a pilot study was conducted on archival FNAs (49 benign, 24 invasive) from women with mammographically suspicious lesions using a newly developed, 5-hour, quantitative, automated cartridge system. We calculated sensitivity, specificity, and area under the receiver-operating characteristic curve (AUC) compared with histopathology for the marker panel. RESULTS: In the discovery cohort, 10 of 25 markers were selected that were highly methylated in breast cancer compared with benign tissues by QM-MSP. In the independent test cohort, this panel yielded an AUC of 0.937 (95% CI = 0.900-0.970). In the FNA pilot, we achieved an AUC of 0.960 (95% CI = 0.883-1.0) using the automated cartridge system. CONCLUSIONS: We developed and piloted a fast and accurate methylation marker-based automated cartridge system to detect breast cancer in FNA samples. This quick ancillary test has the potential to prioritize cancer over benign tissues for expedited pathologic evaluation in poorly resourced countries.

Methylomic Analysis of Ovarian Cancers Identifies Tumor-Specific Alterations Readily Detectable in Early Precursor Lesions.

PURPOSE: High-grade serous ovarian carcinoma (HGSOC) typically remains undiagnosed until advanced stages when peritoneal dissemination has already occurred. Here, we sought to identify HGSOC-specific alterations in DNA methylation and assess their potential to provide sensitive and specific detection of HGSOC at its earliest stages. EXPERIMENTAL DESIGN: MethylationEPIC genome-wide methylation analysis was performed on a discovery cohort comprising 23 HGSOC, 37 non-HGSOC malignant, and 36 histologically unremarkable gynecologic tissue samples. The resulting data were processed using selective bioinformatic criteria to identify regions of high-confidence HGSOC-specific differential methylation. Quantitative methylation-specific real-time PCR (qMSP) assays were then developed for 8 of the top-performing regions and analytically validated in a cohort of 90 tissue samples. Lastly, qMSP assays were used to assess and compare methylation in 30 laser-capture microdissected (LCM) fallopian tube epithelia samples obtained from cancer-free and serous tubal intraepithelial carcinoma (STIC) positive women. RESULTS: Bioinformatic selection identified 91 regions of robust, HGSOC-specific hypermethylation, 23 of which exhibited an area under the receiver-operator curve (AUC) value ≥ 0.9 in the discovery cohort. Seven of 8 top-performing regions demonstrated AUC values between 0.838 and 0.968 when analytically validated by qMSP in a 90-patient cohort. A panel of the 3 top-performing genes (c17orf64, IRX2, and TUBB6) was able to perfectly discriminate HGSOC (AUC 1.0). Hypermethylation within these loci was found exclusively in LCM fallopian tube epithelia from women with STIC lesions, but not in cancer-free fallopian tubes. CONCLUSIONS: A panel of methylation biomarkers can be used to accurately identify HGSOC, even at precursor stages of the disease.

Monitoring of Serum DNA Methylation as an Early Independent Marker of Response and Survival in Metastatic Breast Cancer: TBCRC 005 Prospective Biomarker Study.

Purpose Epigenetic alterations measured in blood may help guide breast cancer treatment. The multisite prospective study TBCRC 005 was conducted to examine the ability of a novel panel of cell-free DNA methylation markers to predict survival outcomes in metastatic breast cancer (MBC) using a new quantitative multiplex assay (cMethDNA). Patients and Methods Ten genes were tested in duplicate serum samples from 141 women at baseline, at week 4, and at first restaging. A cumulative methylation index (CMI) was generated on the basis of six of the 10 genes tested. Methylation cut points were selected to maximize the log-rank statistic, and cross-validation was used to obtain unbiased point estimates. Logistic regression or Cox proportional hazard models were used to test associations between the CMI and progression-free survival (PFS), overall survival (OS), and disease status at first restaging. The added value of the CMI in predicting survival outcomes was evaluated and compared with circulating tumor cells (CellSearch). Results Median PFS and OS were significantly shorter in women with a high CMI (PFS, 2.1 months; OS, 12.3 months) versus a low CMI (PFS, 5.8 months; OS, 21.7 months). In multivariable models, among women with MBC, a high versus low CMI at week 4 was independently associated with worse PFS (hazard ratio, 1.79; 95% CI, 1.23 to 2.60; P = .002) and OS (hazard ratio, 1.75; 95% CI, 1.21 to 2.54; P = .003). An increase in the CMI from baseline to week 4 was associated with worse PFS ( P < .001) and progressive disease at first restaging ( P < .001). Week 4 CMI was a strong predictor of PFS, even in the presence of circulating tumor cells ( P = .004). Conclusion Methylation of this gene panel is a strong predictor of survival outcomes in MBC and may have clinical usefulness in risk stratification and disease monitoring.

Summaries of Affymetrix GeneChip probe level data.

High density oligonucleotide array technology is widely used in many areas of biomedical research for quantitative and highly parallel measurements of gene expression. Affymetrix GeneChip arrays are the most popular. In this technology each gene is typically represented by a set of 11-20 pairs of probes. In order to obtain expression measures it is necessary to summarize the probe level data. Using two extensive spike-in studies and a dilution study, we developed a set of tools for assessing the effectiveness of expression measures. We found that the performance of the current version of the default expression measure provided by Affymetrix Microarray Suite can be significantly improved by the use of probe level summaries derived from empirically motivated statistical models. In particular, improvements in the ability to detect differentially expressed genes are demonstrated.

Comprehensive Pan-Genomic Characterization of Adrenocortical Carcinoma.

We describe a comprehensive genomic characterization of adrenocortical carcinoma (ACC). Using this dataset, we expand the catalogue of known ACC driver genes to include PRKAR1A, RPL22, TERF2, CCNE1, and NF1. Genome wide DNA copy-number analysis revealed frequent occurrence of massive DNA loss followed by whole-genome doubling (WGD), which was associated with aggressive clinical course, suggesting WGD is a hallmark of disease progression. Corroborating this hypothesis were increased TERT expression, decreased telomere length, and activation of cell-cycle programs. Integrated subtype analysis identified three ACC subtypes with distinct clinical outcome and molecular alterations which could be captured by a 68-CpG probe DNA-methylation signature, proposing a strategy for clinical stratification of patients based on molecular markers.

Modular network construction using eQTL data: an analysis of computational costs and benefits.

BACKGROUND: In this paper, we consider analytic methods for the integrated analysis of genomic DNA variation and mRNA expression (also named as eQTL data), to discover genetic networks that are associated with a complex trait of interest. Our focus is the systematic evaluation of the trade-off between network size and network search efficiency in the construction of these networks. RESULTS: We developed a modular approach to network construction, building from smaller networks to larger ones, thereby reducing the search space while including more variables in the analysis. The goal is achieving a lower computational cost while maintaining high confidence in the resulting networks. As demonstrated in our simulation results, networks built in this way have low node/edge false discovery rate (FDR) and high edge sensitivity comparing to greedy search. We further demonstrate our method in a data set of cellular responses to two chemotherapeutic agents: docetaxel and 5-fluorouracil (5-FU), and identify biologically plausible networks that might describe resistances to these drugs. CONCLUSION: In this study, we suggest that guided comprehensive searches for parsimonious networks should be considered as an alternative to greedy network searches.

Do breast cancer cell lines provide a relevant model of the patient tumor methylome?

It is well documented that tumor cells undergo dramatic genetic and epigenetic changes during initial establishment as cell lines and in subsequent serial passaging, and that the resultant cell lines may have evolved significantly from the primary tumors from which they were derived. This has potential implications due to their widespread use in drug response experiments and studies of genomic function. One approach to optimizing the design of such cell line studies is to identify and use the cell lines that faithfully recapitulate critical features of primary tumors. To evaluate the epigenetic fidelity of breast cancer cell lines in the context of primary tumors, we performed methylation profiling of 55 well-characterized breast cancer cell lines on the Illumina HumanMethylation27 BeadChip platform, and compared them to publicly available methylation profiles of primary breast tumors. We found that the DNA methylation profiles of breast cancer cell lines largely retain the features that characterize primary tumors, although there are crucial differences as well. We describe these similarities and differences between primary tumors and breast cancer cell lines in detail, and develop a quantitative measure of similarity that is used to score each cell line with respect to how faithfully its methylation profile mirrors that of primary tumors.