Search

Accurate identification and effective removal of unwanted variation is essential to derive meaningful biological results from RNA sequencing (RNA-seq) data, especially when the data come from large and complex studies. Using RNA-seq data from The Cancer Genome Atlas (TCGA), we examined several sources of unwanted variation and demonstrate here how these can significantly compromise various downstream analyses, including cancer subtype identification, association between gene expression and survival outcomes and gene co-expression analysis. We propose a strategy, called pseudo-replicates of pseudo-samples (PRPS), for deploying our recently developed normalization method, called removing unwanted variation III (RUV-III), to remove the variation caused by library size, tumor purity and batch effects in TCGA RNA-seq data. We illustrate the value of our approach by comparing it to the standard TCGA normalizations on several TCGA RNA-seq datasets. RUV-III with PRPS can be used to integrate and normalize other large transcriptomic datasets coming from multiple laboratories or platforms.

Women carrying germline mutations in BRCA1 are at a substantially elevated risk of breast cancer and their tumors typically have distinctive morphologic features. We hypothesized that constitutional methylation of the BRCA1 promoter region could give rise to such breast cancers in women. We selected 255 women diagnosed with breast cancer before the age of 40 years for whom BRCA1 germline mutations had not been identified. Of them, 52 had five or more of nine BRCA1 mutation-associated morphologic features (group 1), 39 had four (group 2), and 164 had three or less (group 3). The prevalence of detectable BRCA1 promoter methylation in peripheral blood DNA decreased from 31% to 10% to 5% across groups 1–3, respectively (P = 0.000002), and was significantly greater than the 4% frequency in unaffected controls (P = 0.004). Peripheral blood methylation was associated with a 3.5-fold (95% CI, 1.4–10.5) increased risk of having early onset breast cancer. Methylation was consistently mosaic in the peripheral blood where the estimated allelic frequency of BRCA1 promoter methylation ranged from 0.1% to 17%. Group 1 women, but not group 3 women, with detectable methylation of peripheral blood DNA had high levels of BRCA1 promoter methylation of their tumor DNA, indicating that constitutional BRCA1 methylation strongly predisposes toward the development of BRCA1 methylated tumors that then have features resembling BRCA1 mutated tumors. Screening peripheral blood for BRCA1 promoter methylation might thus predict early-onset breast cancers. This raises the possibility of chemoprevention or other intervention to diminish the risk of developing breast cancer in these women. Cancer Prev Res; 4(1); 23–33. ©2010 AACR. Cancer Prev Res; 4(1); 23–33. ©2010 AACR.

Supplementary Table S2 from Constitutional Methylation of the BRCA1 Promoter Is Specifically Associated with BRCA1 Mutation-Associated Pathology in Early-Onset Breast Cancer

Perspective on this Article from Detection of MGMT Promoter Methylation in Normal Individuals Is Strongly Associated with the T Allele of the rs16906252 MGMT Promoter Single Nucleotide Polymorphism

Perspective on This Article from Constitutional Methylation of the BRCA1 Promoter Is Specifically Associated with BRCA1 Mutation-Associated Pathology in Early-Onset Breast Cancer

PDF file - 28K

Supplementary Table 1 Sequencing Primers used for FcGR 2a ( H/R,rs1801274) and FcGR3a (F/V;rs396991) polymorphisms genotyping

HR pathway mutations

In high-grade serous ovarian carcinoma (HGSC), deleterious mutations in DNA repair gene RAD51C are established drivers of defective homologous recombination and are emerging biomarkers of PARP inhibitor (PARPi) sensitivity. RAD51C promoter methylation (meRAD51C) is detected at similar frequencies to mutations, yet its effects on PARPi responses remain unresolved.In this study, three HGSC patient-derived xenograft (PDX) models with methylation at most or all examined CpG sites in the RAD51C promoter show responses to PARPi. Both complete and heterogeneous methylation patterns were associated with RAD51C gene silencing and homologous recombination deficiency (HRD). PDX models lost meRAD51C following treatment with PARPi rucaparib or niraparib, where a single unmethylated copy of RAD51C was sufficient to drive PARPi resistance. Genomic copy number profiling of one of the PDX models using SNP arrays revealed that this resistance was acquired independently in two genetically distinct lineages.In a cohort of 12 patients with RAD51C-methylated HGSC, various patterns of meRAD51C were associated with genomic “scarring,” indicative of HRD history, but exhibited no clear correlations with clinical outcome. Differences in methylation stability under treatment pressure were also observed between patients, where one HGSC was found to maintain meRAD51C after six lines of therapy (four platinum-based), whereas another HGSC sample was found to have heterozygous meRAD51C and elevated RAD51C gene expression (relative to homozygous meRAD51C controls) after only neoadjuvant chemotherapy.As meRAD51C loss in a single gene copy was sufficient to cause PARPi resistance in PDX, methylation zygosity should be carefully assessed in previously treated patients when considering PARPi therapy.Significance:Homozygous RAD51C methylation is a positive predictive biomarker for sensitivity to PARP inhibitors, whereas a single unmethylated gene copy is sufficient to confer resistance.

Supplementary Figure 1 - PDF file 91K, Schematic diagram of experimental procedure and bioinformatic analysis

All Supplementary Methods for the manuscript

PDF file - 1609K

Supplementary Methods, Figures 1-5. Supplementary Figure 1. Outline of cohort selection. Supplementary Figure 2. Example of whole-genome single-nucleotide polymorphism (SNP)array profile and genome alteration print (GAP). Supplementary Figure 3. Pre-treatment serum CA125 levels. Supplementary Figure 4. Progression-free and overall survival in patients with SBT-EOC and EOC without adjacent SBT. Supplementary Figure 5. Progression-free and overall survival of patients with Grade 3 SBTEOC and Grade 3 EOC without adjacent SBT

Purpose: Women with epithelial ovarian cancer generally have a poor prognosis; however, a subset of patients has an unexpected dramatic and durable response to treatment. We sought to identify clinical, pathological, and molecular determinants of exceptional survival in women with high-grade serous cancer (HGSC), a disease associated with the majority of ovarian cancer deaths.Experimental Design: We evaluated the histories of 2,283 ovarian cancer patients and, after applying stringent clinical and pathological selection criteria, identified 96 with HGSC that represented significant outliers in terms of treatment response and overall survival. Patient samples were characterized immunohistochemically and by genome sequencing.Results: Different patterns of clinical response were seen: long progression-free survival (Long-PFS), multiple objective responses to chemotherapy (Multiple Responder), and/or greater than 10-year overall survival (Long-Term Survivors). Pathogenic germline and somatic mutations in genes involved in homologous recombination (HR) repair were enriched in all three groups relative to a population-based series. However, 29% of 10-year survivors lacked an identifiable HR pathway alteration, and tumors from these patients had increased Ki-67 staining. CD8+ tumor-infiltrating lymphocytes were more commonly present in Long-Term Survivors. RB1 loss was associated with long progression-free and overall survival. HR deficiency and RB1 loss were correlated, and co-occurrence was significantly associated with prolonged survival.Conclusions: There was diversity in the clinical trajectory of exceptional survivors associated with multiple molecular determinants of exceptional outcome in HGSC patients. Concurrent HR deficiency and RB1 loss were associated with favorable outcomes, suggesting that co-occurrence of specific mutations might mediate durable responses in such patients. Clin Cancer Res; 24(3); 569–80. ©2017 AACR.See related commentary by Peng and Mills, p. 508

Purpose: High-grade serous carcinoma (HGSC) accounts for the majority of epithelial ovarian cancer deaths. Genomic and functional data suggest that approximately half of unselected HGSC have disruption of the BRCA pathway and defects in homologous recombination repair (HRR). Pathway disruption is regarded as imparting a BRCAness phenotype. We explored the molecular changes in HGSC arising in association with specific BRCA1/BRCA2 somatic or germline mutations and in those with BRCA1 DNA promoter methylation.Experimental Design: We describe gene expression and copy number analysis of two large cohorts of HGSC in which both germline and somatic inactivation of HRR has been measured.Results:BRCA1 disruptions were associated with the C2 (immunoreactive) molecular subtype of HGSC, characterized by intense intratumoral T-cell infiltration. We derived and validated a predictor of BRCA1 mutation or methylation status, but could not distinguish BRCA2 from wild-type tumors. DNA copy number analysis showed that cases with BRCA1 mutation were significantly associated with amplification both at 8q24 (frequencies: BRCA1 tumors 50%, BRCA2 tumors 32%, and wild-type tumors 9%) and regions of the X-chromosome specifically dysregulated in basal-like breast cancer (BLBC; BRCA1 62%, BRCA2 34%, and wild-type 35%). Tumors associated with BRCA1/BRCA2 mutations shared a negative association with amplification at 19p13 (BRCA1 0%, BRCA2 3%, and wild-type 20%) and 19q12 (BRCA1 6%, BRCA2 3%, and wild-type 29%).Conclusion: The molecular differences between tumors associated with BRCA1 compared with BRCA2 mutations are in accord with emerging clinical and pathologic data and support a growing appreciation of the relationship between HGSC and BLBC. Clin Cancer Res; 19(13); 3474–84. ©2013 AACR.

Supplementary tables 1, 3 and 4 - XLSX file 31K, 1. Performance summary of exome sequencing and capture kit used for each sample 3. List of genes for each melanoma subtype 4. List of potential actionable mutations in BRAF/NRAS wild-type tumors

Revised Supplementary Data, containing Supplementary Methods, Supplementary Figures, Supplementary Tables Supplementary Figure S1. Outline of cohort selection and analyses. Supplementary Figure S2. Clinical response and therapy course of 96 patients with exceptional responses to chemotherapy. Supplementary Figure S3. Distribution and type of TP53 mutations. Supplementary Figure S4. RB1 protein expression altered by genomic inactivation. Supplementary Figure S5. Characterization of CD8 and Ki-67 in tumors according to homologous recombination mutation status. Supplementary Table S2 Immunohistochemical analysis: primary antibodies and staining conditions Supplementary Table S3 Homologous recombination and DNA repair panel Supplementary Table S6 Comparison of molecular alteration prevalence between clinical subgroups Supplementary Table S7 Patient characteristics of tissue microarray cohort

All supplementary tables for the manuscript

PDF file - 338K

Supplementary Highlighted Methods. This document contains the revised supplementary methods and results with changes highlighted

Supplementary Table 2 - PDF file 16632K, Single nucleotide variants identified for each sample

Purpose: Two germline Fc-γ receptor (FCGR) polymorphisms, rs1801274 [FCGR2A;His(H)131Arg(R)] and rs396991 [FCGR3A;Phe(F)158Val(V)] produce altered proteins through amino acid substitutions; both are reported to be associated with cetuximab-related outcomes. We performed a validation of these polymorphisms in NCIC CTG CO.17, a randomized trial of cetuximab monotherapy in refractory, metastatic colorectal cancer expressing EGFR.Experimental Design: DNA extracted from formalin-fixed paraffin-embedded tissue was genotyped. In addition to log-rank tests, Cox proportional hazard models assessed their relationships with overall (OS) and progression-free survival (PFS), adjusting for clinically important prognostic factors, along with a polymorphism–treatment arm interaction term.Results: Somatic KRAS status was wild-type for exon 2 in 153 (52%) of 293 patients, from whom tumor DNA was available. For FCGR2A H/H, a genotype–treatment interaction for KRAS wild-type patients was observed for OS (P = 0.03). In KRAS wild-type patients carrying FCGR2A H/H, cetuximab (vs. no cetuximab) improved survival substantially, with adjusted HRs (aHR) of 0.36 (OS) and 0.19 (PFS) and absolute benefits of 5.5 months (OS; P = 0.003) and 3.7 months (PFS; P = 0.02). In contrast, patients carrying FCGR2A R alleles (H/R or R/R) had aHRs of only 0.78 (OS; 2.8-month benefit) and 0.53 (PFS; 1.6-month benefit). No relationships were found for rs396991 (FCGR3A).Conclusions: In the CO.17 trial, cetuximab worked best for patients with KRAS wild-type colorectal cancers carrying FCGR2A H/H genotypes. Significantly lower benefits were observed in patients carrying germline FCGR2A R alleles. Clin Cancer Res; 22(10); 2435–44. ©2016 AACR.

Supplementary Tables 1-7. Supplementary Table 1. SBT-EOC cases used for molecular analyses Supplementary Table 2. Mutations Screened on the Sequenom/OncoMap3 Platform Supplementary Table 3. Primer sequences for the indicated exons of TP53, ERBB2 and NRAS Supplementary Table 4. Clinico-pathologic characteristics of the subset of SBT-EOC used for molecular analysis compared to the entire SBT-EOC cohort Supplementary Table 5A. Somatic mutations found in borderline and invasive ovarian tumors Supplementary Table 5B. Genes in which somatic mutations were found according to grade and co-existing serous borderline regions Supplementary Table 6. Genes differentially expressed between paired borderline and invasive regions of SBTEOC (p

Purpose: Low-grade serous ovarian carcinomas (LGSC) are Ras pathway-mutated, TP53 wild-type, and frequently associated with borderline tumors. Patients with LGSCs respond poorly to platinum-based chemotherapy and may benefit from pathway-targeted agents. High-grade serous carcinomas (HGSC) are TP53-mutated and are thought to be rarely associated with borderline tumors. We sought to determine whether borderline histology associated with grade 2 or 3 carcinoma was an indicator of Ras mutation, and we explored the molecular relationship between coexisting invasive and borderline histologies.Experimental Design: We reviewed >1,200 patients and identified 102 serous carcinomas with adjacent borderline regions for analyses, including candidate mutation screening, copy number, and gene expression profiling.Results: We found a similar frequency of low, moderate, and high-grade carcinomas with coexisting borderline histology. BRAF/KRAS alterations were common in LGSC; however, we also found recurrent NRAS mutations. Whereas borderline tumors harbored BRAF/KRAS mutations, NRAS mutations were restricted to carcinomas, representing the first example of a Ras oncogene with an obligatory association with invasive serous cancer. Coexisting borderline and invasive components showed nearly identical genomic profiles. Grade 2 cases with coexisting borderline included tumors with molecular features of LGSC, whereas others were typical of HGSC. However, all grade 3 carcinomas with coexisting borderline histology were molecularly indistinguishable from typical HGSC.Conclusion: Our findings suggest that NRAS is an oncogenic driver in serous ovarian tumors. We demonstrate that borderline histology is an unreliable predictor of Ras pathway aberration and underscore an important role for molecular classification in identifying patients that may benefit from targeted agents. Clin Cancer Res; 20(24); 6618–30. ©2014 AACR.

Purpose: The mutation load in melanoma is generally high compared with other tumor types due to extensive UV damage. Translation of exome sequencing data into clinically relevant information is therefore challenging. This study sought to characterize mutations identified in primary cutaneous melanomas and correlate these with clinicopathologic features.Experimental Design: DNA was extracted from 34 fresh-frozen primary cutaneous melanomas and matched peripheral blood. Tumor histopathology was reviewed by two dermatopathologists. Exome sequencing was conducted and mutation rates were correlated with age, sex, tumor site, and histopathologic variables. Differences in mutations between categories of solar elastosis, pigmentation, and BRAF/NRAS mutational status were investigated.Results: The average mutation rate was 12 per megabase, similar to published results in metastases. The average mutation rate in severely sun damaged (SSD) skin was 21 per Mb compared with 3.8 per Mb in non-SSD skin (P = 0.001). BRAF/NRAS wild-type (WT) tumors had a higher average mutation rate compared with BRAF/NRAS–mutant tumors (27 vs. 5.6 mutations per Mb; P = 0.0001). Tandem CC>TT/GG>AA mutations comprised 70% of all dinucleotide substitutions and were more common in tumors arising in SSD skin (P = 0.0008) and in BRAF/NRAS WT tumors (P = 0.0007). Targetable and potentially targetable mutations in WT tumors, including NF1, KIT, and NOTCH1, were spread over various signaling pathways.Conclusion: Melanomas arising in SSD skin have higher mutation loads and contain a spectrum of molecular subtypes compared with BRAF- and NRAS-mutant tumors indicating multigene screening approaches and combination therapies may be required for management of these patients. Clin Cancer Res; 19(17); 4589–98. ©2013 AACR.

All supplementary figures for the manuscript

Graft-derived cell-free DNA (gdcfDNA) quantification is a promising, minimally invasive tool for detecting acute T cell-mediated rejection (ATCMR) following liver transplantation (LT). We investigated the utility of measuring hepatocyte-specific methylation in cfDNA (HS-cfDNA) to quantify gdcfDNA, examining its accuracy in detecting ATCMR in a prospective, cross-sectional study. Blood was collected from LT recipients immediately prior to graft biopsy for suspected rejection. HS-cfDNA was quantified using droplet-digital polymerase chain reaction. Prebiopsy liver function tests (LFTs) and HS-cfDNA levels were correlated with biopsy results and the primary outcome of treated biopsy-proven acute rejection (tBPAR). A total of 51 patients were recruited; 37 had evidence of rejection on biopsy and 20 required treatment. As much as 11 patients needed inpatient treatment for rejection. HS-cfDNA significantly outperformed LFTs in identifying patients with tBPAR, particularly those needing inpatient treatment (area under the curve, 73.0%; 95% confidence interval, 55.4%-90.6%; P = 0.01). At a threshold of33.5% of the total cfDNA fraction, HS-cfDNA had a specificity of 97%, correctly excluding tBPAR in 30/31 patients. Quantifying graft-specific methylation in cfDNA has a major advantage over previous gdcfDNA techniques: it does not require genotyping/sequencing, lending it greater feasibility for translation into transplantation care. Low levels of HS-cfDNA were a strong negative predictor for tBPAR (negative predictive value, 86%) and may have a future role in triaging patients prior to invasive graft biopsies.

(1) Background: Epithelial&ndash, mesenchymal plasticity (EMP) is a dynamic process whereby epithelial carcinoma cells reversibly acquire morphological and invasive characteristics typical of mesenchymal cells. Identifying the methylation differences between epithelial and mesenchymal states may assist in the identification of optimal DNA methylation biomarkers for the blood-based monitoring of cancer. (2) Methods: Methylation-sensitive high-resolution melting (MS-HRM) was used to examine the promoter methylation status of a panel of established and novel markers in a range of breast cancer cell lines spanning the epithelial&ndash, mesenchymal spectrum. Pyrosequencing was used to validate the MS-HRM results. (3) Results: VIM, DKK3, and CRABP1 were methylated in the majority of epithelial breast cancer cell lines, while methylation of GRHL2, MIR200C, and CDH1 was restricted to mesenchymal cell lines. Some markers that have been used to assess minimal residual disease such as AKR1B1 and APC methylation proved to be specific for epithelial breast cell lines. However, RASSF1A, RAR&beta, TWIST1, and SFRP2 methylation was seen in both epithelial and mesenchymal cell lines, supporting their suitability for a multimarker panel. (4) Conclusions: Profiling DNA methylation shows a distinction between epithelial and mesenchymal phenotypes. Understanding how DNA methylation varies between epithelial and mesenchymal phenotypes may lead to more rational selection of methylation-based biomarkers for circulating tumour DNA analysis.

BACKGROUND AND AIM: The role of circulating mitochondrial DNA (cmtDNA) in transplantation remains to be elucidated. cmtDNA may be released into the circulation as a consequence of liver injury; yet recent work also suggests a causative role for cmtDNA leading to hepatocellular injury. We hypothesized that elevated cmtDNA would be associated with adverse events after liver transplantation (LT) and conducted an observational cohort study. METHODS: Twenty-one patients were enrolled prospectively prior to LT. RESULTS: Postoperative complications were observed in 47.6% (n = 10). Seven patients (33.3%) had early allograft dysfunction (EAD), and six patients (28.5%) experienced acute cellular rejection within 6 months of LT. cmtDNA levels were significantly elevated in all recipients after LT compared with healthy controls and preoperative samples (1 361 937 copies/mL [IQR 586 781-3 399 687] after LT; 545 531 copies/mL [IQR 238 562-1 381 015] before LT; and 194 562 copies/mL [IQR 182 359-231 515] in healthy controls) and returned to normal levels by 5 days after transplantation. cmtDNA levels were particularly elevated in those who developed EAD in the early postoperative period (P

In high-grade serous ovarian carcinoma (HGSC), deleterious mutations in DNA repair gene RAD51C are established drivers of defective homologous recombination and are emerging biomarkers of PARP inhibitor (PARPi) sensitivity. RAD51C promoter methylation (meRAD51C) is detected at similar frequencies to mutations, yet its effects on PARPi responses remain unresolved. In this study, three HGSC patient-derived xenograft (PDX) models with methylation at most or all examined CpG sites in the RAD51C promoter show responses to PARPi. Both complete and heterogeneous methylation patterns were associated with RAD51C gene silencing and homologous recombination deficiency (HRD). PDX models lost meRAD51C following treatment with PARPi rucaparib or niraparib, where a single unmethylated copy of RAD51C was sufficient to drive PARPi resistance. Genomic copy number profiling of one of the PDX models using SNP arrays revealed that this resistance was acquired independently in two genetically distinct lineages. In a cohort of 12 patients with RAD51C-methylated HGSC, various patterns of meRAD51C were associated with genomic “scarring,” indicative of HRD history, but exhibited no clear correlations with clinical outcome. Differences in methylation stability under treatment pressure were also observed between patients, where one HGSC was found to maintain meRAD51C after six lines of therapy (four platinum-based), whereas another HGSC sample was found to have heterozygous meRAD51C and elevated RAD51C gene expression (relative to homozygous meRAD51C controls) after only neoadjuvant chemotherapy. As meRAD51C loss in a single gene copy was sufficient to cause PARPi resistance in PDX, methylation zygosity should be carefully assessed in previously treated patients when considering PARPi therapy. Significance: Homozygous RAD51C methylation is a positive predictive biomarker for sensitivity to PARP inhibitors, whereas a single unmethylated gene copy is sufficient to confer resistance.

Due to advances in modern medicine, liver transplantation has revolutionised the prognosis of many previously incurable liver diseases. This progress has largely been due to advances in immunosuppressant therapy. However, despite the judicious use of immunosuppression, many liver transplant recipients still experience complications such as rejection, which necessitates diagnosis via invasive liver biopsy. There is a clear need for novel, minimally-invasive tests to optimise immunosuppression and improve patient outcomes. An emerging biomarker in this ‘‘precision medicine’‘ liver transplantation field is that of donor-specific cell free DNA. In this review, we detail the background and methods of detecting this biomarker, examine its utility in liver transplantation and discuss future research directions that may be most impactful.