Defining the relationship between cellular and extracellular vesicle (EV) content in breast cancer via an integrative multi-omic analysis.

Much recent research has been dedicated to exploring the utility of extracellular vesicles (EVs) as circulating disease biomarkers. Underpinning this work is the assumption that the molecular cargo of EVs directly reflects the originating cell. Few attempts have been made, however, to empirically validate this on the -omic level. To this end, we have performed an integrative multi-omic analysis of a panel of breast cancer cell lines and corresponding EVs. Whole transcriptome analysis validated that the cellular transcriptome remained stable when cultured cells are transitioned to low serum or serum-free medium for EV collection. Transcriptomic profiling of the isolated EVs indicated a positive correlation between transcript levels in cells and EVs, including disease-associated transcripts. Analysis of the EV proteome verified that HER2 protein is present in EVs, however neither the estrogen (ER) nor progesterone (PR) receptor proteins are detected regardless of cellular expression. Using multivariate analysis, we derived an EV protein signature to infer cellular patterns of ER and HER2 expression, though the ER protein could not be directly detected. Integrative analyses affirmed that the EV proteome and transcriptome captured key phenotypic hallmarks of the originating cells, supporting the potential of EVs for non-invasive monitoring of breast cancers.

Transcriptional comparison of testicular adrenal rest tumors with fetal and adult tissues.

Background: Testicular adrenal rest tumors (TART) are a common complication of unknown cellular origin in patients with congenital adrenal hyperplasia (CAH). These benign tumors have both adrenal and testicular characteristics and are hypothesized to either derive from cells of adrenal origin from the fetal adrenogonadal primordium or by atypical differentiation of adult Leydig-progenitor cells. Objective: This study aims to unravel the identity and etiology of TART. Methods: Co-expression of adrenal-specific CYP11B1 and Leydig cell-specific HSD17B3 in TART was studied using immunohistochemistry. We studied the possibility of TART being derived from atypical differentiation of adult Leydig-progenitor cells by the quantification of adrenal-specific enzyme expression upon adrenocorticotrophic hormone (ACTH)-like stimulation of ex vivo cultured platelet-derived growth factor receptor alpha-positive cells. By comparing the transcriptome of TART (n = 16) with the transcriptome of fetal adrenal (n = 13), fetal testis (n = 5), adult adrenal (n = 11), and adult testis (n = 10) tissues, we explored the identity of TART. Results: We demonstrate co-expression of adrenal-specific CYP11B1 and testis-specific HSD17B3 in TART cells, indicating the existence of a distinct TART cell exhibiting both adrenal and testicular characteristics. Ex vivo cultured adult Leydig-progenitor cells did not express the ACTH-receptor MC2R but did express CYP11B1 upon stimulation. Unsupervised clustering of transcriptome data showed that TART was most similar to adult adrenal tissue, followed by adult testis tissue, and least similar to either fetal tissue. Conclusion: Our data suggest that TART is induced - most likely via activation of a cAMP/protein kinase A-dependent receptor - from a progenitor cell into a unique mature adrenal-like cell type, sometimes exhibiting both adrenal and testicular features.

Comprehensive methylome sequencing reveals prognostic epigenetic biomarkers for prostate cancer mortality.

BACKGROUND: Prostate cancer is a clinically heterogeneous disease with a subset of patients rapidly progressing to lethal-metastatic prostate cancer. Current clinicopathological measures are imperfect predictors of disease progression. Epigenetic changes are amongst the earliest molecular changes in tumourigenesis. To find new prognostic biomarkers to enable earlier intervention and improved outcomes, we performed methylome sequencing of DNA from patients with localised prostate cancer and long-term clinical follow-up. METHODS: We used whole-genome bisulphite sequencing (WGBS) to comprehensively map and compare DNA methylation of radical prostatectomy tissue between patients with lethal disease (n = 7) and non-lethal (n = 8) disease (median follow-up 19.5 years). Validation of differentially methylated regions (DMRs) was performed in an independent cohort (n = 185, median follow-up 15 years) using targeted multiplex bisulphite sequencing of candidate regions. Survival was assessed via univariable and multivariable analyses including clinicopathological measures (log-rank and Cox regression models). RESULTS: WGBS data analysis identified cancer-specific methylation patterns including CpG island hypermethylation, and hypomethylation of repetitive elements, with increasing disease risk. We identified 1420 DMRs associated with prostate cancer-specific mortality (PCSM), which showed enrichment for gene sets downregulated in prostate cancer and de novo methylated in cancer. Through comparison with public prostate cancer datasets, we refined the DMRs to develop an 18-gene prognostic panel. Applying this panel to an independent cohort, we found significant associations between PCSM and hypermethylation at EPHB3, PARP6, TBX1, MARCH6 and a regulatory element within CACNA2D4. Strikingly in a multivariable model, inclusion of CACNA2D4 methylation was a better predictor of PCSM versus grade alone (Harrell's C-index: 0.779 vs. 0.684). CONCLUSIONS: Our study provides detailed methylome maps of non-lethal and lethal prostate cancer and identifies novel genic regions that distinguish these patient groups. Inclusion of our DNA methylation biomarkers with existing clinicopathological measures improves prognostic models of prostate cancer mortality, and holds promise for clinical application.

Opportunities for Early Cancer Detection: The Rise of ctDNA Methylation-Based Pan-Cancer Screening Technologies.

The efficiency of conventional screening programs to identify early-stage malignancies can be limited by the low number of cancers recommended for screening as well as the high cumulative false-positive rate, and associated iatrogenic burden, resulting from repeated multimodal testing. The opportunity to use minimally invasive liquid biopsy testing to screen asymptomatic individuals at-risk for multiple cancers simultaneously could benefit from the aggregated diseases prevalence and a fixed specificity. Increasing both latter parameters is paramount to mediate high positive predictive value-a useful metric to evaluate a screening test accuracy and its potential harm-benefit. Thus, the use of a single test for multi-cancer early detection (stMCED) has emerged as an appealing strategy for increasing early cancer detection rate efficiency and benefit population health. A recent flurry of these stMCED technologies have been reported for clinical potential; however, their development is facing unique challenges to effectively improve clinical cost-benefit. One promising avenue is the analysis of circulating tumour DNA (ctDNA) for detecting DNA methylation biomarker fingerprints of malignancies-a hallmark of disease aetiology and progression holding the potential to be tissue- and cancer-type specific. Utilizing panels of epigenetic biomarkers could potentially help to detect earlier stages of malignancies as well as identify a tumour of origin from blood testing, useful information for follow-up clinical decision making and subsequent patient care improvement. Overall, this review collates the latest and most promising stMCED methodologies, summarizes their clinical performances, and discusses the specific requirements multi-cancer tests should meet to be successfully implemented into screening guidelines.

Dynamic Monitoring of EMT in CTCs as an Indicator of Cancer Metastasis.

Epithelial to mesenchymal transition (EMT) results in the genesis of circulating tumor cells (CTCs) from tumor sites and promotes the metastatic capability of CTCs in circulation. In this study, we develop a multiplex surface-enhanced Raman scattering nanotechnology for comprehensive characterization of EMT-associated phenotypes in CTCs, to monitor cancer metastasis. We observe the downregulation of the CTC marker (EpCAM) and the epithelial marker (E-cadherin), as well as the upregulation of a mesenchymal marker (N-cadherin) and a stem cell marker (ABCB5) during the transforming growth factor-ß-induced EMT process in breast cancer cell line models. Additionally, we also find changes in the heterogeneity levels of these selected markers in cells. With this method, we successfully detect the presence of disease in samples from breast cancer patients and characterize EMT-associated phenotypes in their CTCs. Overall, this approach and findings provide a new means for monitoring the EMT process in cancer, insights into the detailed mechanistic progress of the diseases, and have potential for detecting the early occurrence of cancer metastasis.

Network mapping of primary CD34+ cells by Ampliseq based whole transcriptome targeted resequencing identifies unexplored differentiation regulatory relationships.

With the exception of a few master transcription factors, regulators of neutrophil maturation are poorly annotated in the intermediate phenotypes between the granulocyte-macrophage progenitor (GMP) and the mature neutrophil phenotype. Additional challenges in identifying gene expression regulators in differentiation pathways relate to challenges wherein starting cell populations are heterogeneous in lineage potential and development, are spread across various states of quiescence, as well as sample quality and input limitations. These factors contribute to data variability make it difficult to draw simple regulatory inferences. In response we have applied a multi-omics approach using primary blood progenitor cells primed for homogeneous proliferation and granulocyte differentiation states which combines whole transcriptome resequencing (Ampliseq RNA) supported by droplet digital PCR (ddPCR) validation and mass spectrometry-based proteomics in a hypothesis-generation study of neutrophil differentiation pathways. Primary CD34+ cells isolated from human cord blood were first precultured in non-lineage driving medium to achieve an active, proliferating phenotype from which a neutrophil primed progenitor was isolated and cultured in neutrophil lineage supportive medium. Samples were then taken at 24-hour intervals over 9 days and analysed by Ampliseq RNA and mass spectrometry. The Ampliseq dataset depth, breadth and quality allowed for several unexplored transcriptional regulators and ncRNAs to be identified using a combinatorial approach of hierarchical clustering, enriched transcription factor binding motifs, and network mapping. Network mapping in particular increased comprehension of neutrophil differentiation regulatory relationships by implicating ARNT, NHLH1, PLAG1, and 6 non-coding RNAs associated with PU.1 regulation as cell-engineering targets with the potential to increase total neutrophil culture output. Overall, this study develops and demonstrates an effective new hypothesis generation methodology for transcriptome profiling during differentiation, thereby enabling identification of novel gene targets for editing interventions.

Regulation of Canonical Oncogenic Signaling Pathways in Cancer via DNA Methylation.

Disruption of signaling pathways that plays a role in the normal development and cellular homeostasis may lead to the dysregulation of cellular signaling and bring about the onset of different diseases, including cancer. In addition to genetic aberrations, DNA methylation also acts as an epigenetic modifier to drive the onset and progression of cancer by mediating the reversible transcription of related genes. Although the role of DNA methylation as an alternative driver of carcinogenesis has been well-established, the global effects of DNA methylation on oncogenic signaling pathways and the presentation of cancer is only emerging. In this article, we introduced a differential methylation parsing pipeline (MethylMine) which mined for epigenetic biomarkers based on feature selection. This pipeline was used to mine for biomarkers, which presented a substantial difference in methylation between the tumor and the matching normal tissue samples. Combined with the Data Integration Analysis for Biomarker discovery (DIABLO) framework for machine learning and multi-omic analysis, we revisited the TCGA DNA methylation and RNA-Seq datasets for breast, colorectal, lung, and prostate cancer, and identified differentially methylated genes within the NRF2-KEAP1/PI3K oncogenic pathway, which regulates the expression of cytoprotective genes, that serve as potential therapeutic targets to treat different cancers.

Comprehensive evaluation of targeted multiplex bisulphite PCR sequencing for validation of DNA methylation biomarker panels.

BACKGROUND: DNA methylation is a well-studied epigenetic mark that is frequently altered in diseases such as cancer, where specific changes are known to reflect the type and severity of the disease. Therefore, there is a growing interest in assessing the clinical utility of DNA methylation as a biomarker for diagnosing disease and guiding treatment. The development of an accurate loci-specific methylation assay, suitable for use on low-input clinical material, is crucial for advancing DNA methylation biomarkers into a clinical setting. A targeted multiplex bisulphite PCR sequencing approach meets these needs by allowing multiple DNA methylated regions to be interrogated simultaneously in one experiment on limited clinical material. RESULTS: Here, we provide an updated protocol and recommendations for multiplex bisulphite PCR sequencing (MBPS) assays for target DNA methylation analysis. We describe additional steps to improve performance and reliability: (1) pre-sequencing PCR optimisation which includes assessing the optimal PCR cycling temperature and primer concentration and (2) post-sequencing PCR optimisation to achieve uniform coverage of each amplicon. We use a gradient of methylated controls to demonstrate how PCR bias can be assessed and corrected. Methylated controls also allow assessment of the sensitivity of methylation detection for each amplicon. Here, we show that the MBPS assay can amplify as little as 0.625 ng starting DNA and can detect methylation differences of 1% with a sequencing coverage of 1000 reads. Furthermore, the multiplex bisulphite PCR assay can comprehensively interrogate multiple regions on 1-5 ng of formalin-fixed paraffin-embedded DNA or circulating cell-free DNA. CONCLUSIONS: The MBPS assay is a valuable approach for assessing methylated DNA regions in clinical samples with limited material. The optimisation and additional quality control steps described here improve the performance and reliability of this method, advancing it towards potential clinical applications in biomarker studies.

Optimizing Size Exclusion Chromatography for Extracellular Vesicle Enrichment and Proteomic Analysis from Clinically Relevant Samples.

The field of extracellular vesicle (EV) research has rapidly expanded in recent years, with particular interest in their potential as circulating biomarkers. Proteomic analysis of EVs from clinical samples is complicated by the low abundance of EV proteins relative to highly abundant circulating proteins such as albumin and apolipoproteins. To overcome this, size exclusion chromatography (SEC) has been proposed as a method to enrich EVs whilst depleting protein contaminants; however, the optimal SEC parameters for EV proteomics have not been thoroughly investigated. Here, quantitative evaluation and optimization of SEC are reported for separating EVs from contaminating proteins. Using a synthetic model system followed by cell line-derived EVs, it is found that a 10 mL Sepharose 4B column in PBS produces optimal resolution of EVs from background protein. By spiking-in cancer cell-derived EVs to healthy plasma, it is shown that some cancer EV-associated proteins are detectable by nano-LC-MS/MS when as little as 1% of the total plasma EV number are derived from a cancer cell line. These results suggest that an optimized SEC and nanoLC-MS/MS workflow may be sufficiently sensitive for disease EV protein biomarker discovery from patient-derived clinical samples.

A SERS microfluidic platform for targeting multiple soluble immune checkpoints.

Immune checkpoint blockade therapies are promising next generation immunotherapeutic treatments for cancer. Whilst sequential solid biopsies are an invaluable source of prognostic information, they are not feasible for monitoring therapeutic outcomes over time. Monitoring soluble immune checkpoint markers expression in body fluids could potentially be a better alternative. Current methods (e.g. ELISA) for detecting immune-checkpoint proteins mostly rely on the use of monoclonal antibodies which are expensive and time-consuming to manufacture and isolate. Herein, we report an integrated surface enhanced Raman scattering (SERS)-microfluidics device for the detection of immune checkpoint proteins which involves the use of i) nano yeast single chain variable fragment (scFv) as a promising alternative to monoclonal antibodies providing high stability at relative low-cost and simplicity for production, ii) graphene oxide functionalised surface to reduces the bio functionalization steps, thus avoiding the general paradigm of biotin-streptavidin chemistry and iii) a microfluidic platform enabling alternating current electrohydrodynamics (ac-EHD) induced nanomixing to enhance the target scFv binding and minimize the non-specific interactions. Specific and multiplex detection of immune checkpoint biomarkers is achieved by SERS based spectral encoding. Using this platform, we successfully demonstrated the detection of clinically relevant soluble immune checkpoints PD-1, PD-L1 and LAG-3 from as low as 100 fg/mL of analytes spiked in human serum.

Epigenetically reprogrammed methylation landscape drives the DNA self-assembly and serves as a universal cancer biomarker.

Epigenetic reprogramming in cancer genomes creates a distinct methylation landscape encompassing clustered methylation at regulatory regions separated by large intergenic tracks of hypomethylated regions. This methylation landscape that we referred to as Methylscape is displayed by most cancer types, thus may serve as a universal cancer biomarker. To-date most research has focused on the biological consequences of DNA Methylscape changes whereas its impact on DNA physicochemical properties remains unexplored. Herein, we examine the effect of levels and genomic distribution of methylcytosines on the physicochemical properties of DNA to detect the Methylscape biomarker. We find that DNA polymeric behaviour is strongly affected by differential patterning of methylcytosine, leading to fundamental differences in DNA solvation and DNA-gold affinity between cancerous and normal genomes. We exploit these Methylscape differences to develop simple, highly sensitive and selective electrochemical or colorimetric one-step assays for the detection of cancer. These assays are quick, i.e., analysis time ≤10 minutes, and require minimal sample preparation and small DNA input.

MiR-29b-1-5p is altered in BRCA1 mutant tumours and is a biomarker in basal-like breast cancer.

Depletion of BRCA1 protein in mouse mammary glands results in defects in lactational development and increased susceptibility to mammary cancer. Extensive work has focussed on the role of BRCA1 in the normal breast and in the development of breast cancer, the cell of origin for BRCA1 tumours and the protein-coding genes altered in BRCA1 deficient cells. However, the role of non-coding RNAs in BRCA1-deficient cells is poorly understood. To evaluate miRNA expression in BRCA1 deficient mammary cells, RNA sequencing was performed on the mammary glands of Brca1 knockout mice. We identified 140 differentially expressed miRNAs, 9 of which were also differentially expressed in human BRCA1 breast tumours or familial non-BRCA1 patients and during normal gland development. We show that BRCA1 binds to putative cis-elements in promoter regions of the miRNAs with the potential to regulate their expression, and that four miRNAs (miR-29b-1-5p, miR-664, miR-16-2 and miR-744) significantly stratified the overall survival of basal-like tumours. Importantly the prognostic value of miR-29b-1-5p was higher in significance than several commonly used clinical biomarkers. These results emphasise the role of Brca1 in modulating expression of miRNAs and highlights the potential for BRCA1 regulated miRNAs to be informative biomarkers associated with BRCA1 loss and survival in breast cancer.

Parallel profiling of cancer cells and proteins using a graphene oxide functionalized ac-EHD SERS immunoassay.

Circulating biomarkers have emerged as promising non-invasive, real-time surrogates for cancer diagnosis, prognosis and monitoring of the therapeutic response. Current bio-sensing techniques mostly involve detection of either circulating cells or proteins which are inadequate in unfolding complex pathologic transformations. Herein, we report parallel detection of cellular and molecular markers (protein) for cancer using a multiplex platform featuring (i) graphene oxide (GO) functionalization that increases the active surface area and more importantly reduces the functionalization steps for rapid detection, (ii) alternating-current electrohydrodynamic (ac-EHD) fluid flow that provides delicate micro-mixing to enhance target-sensor interactions thereby minimizing non-specific binding and (iii) surface enhanced Raman scattering (SERS) for multiplex detection. We find that our platform possesses high sensitivity for detecting both proteins and cells. More importantly, this platform not only detects the cancer cells but also can simultaneously monitor the heterogeneous expression of cell surface proteins which could be clinically useful to determine effective patient therapy. We demonstrate the specific and sensitive detection of breast cancer cells from a mixture of non-target cells and report the heterogeneous expression of human epidermal growth factor receptor 2 (HER2) proteins on the individual cancer cell surface. Concurrently, we detect as low as 100 fg mL-1 HER2 and Mucin 16 proteins spiked in blood serum.

Detection of aberrant protein phosphorylation in cancer using direct gold-protein affinity interactions.

Protein phosphorylation is one of the most prominent post-translational mechanisms for protein regulation, which is frequently impaired in cancer. Through the covalent addition of phosphate groups to certain amino-acids, the interactions of former residues with nearby amino-acids are drastically altered, resulting in major changes of protein conformation that impacts its biological function. Herein, we report that these conformational changes can also disturb the protein's ability to interact with and adsorb onto bare gold surfaces. We exploited this feature to develop a simple electrochemical method for detecting the aberrant phosphorylation of EGFR protein in several lung cancer cell lines. This method, which required as low as 10ng/µL (i.e., 50ng) of purified EGFR protein, also enabled monitoring cell sensitivity to tyrosine kinase inhibitors (TKI) - a common drug used for restoring the function of aberrantly phosphorylated proteins in lung cancer. The reported strategy based on direct gold-protein affinity interactions avoids the conventional paradigm of requiring a phospho-specific antibody for detection and could be a potential alternative of widely used mass spectrometry.

Rapid Molecular Profiling of Myeloproliferative Neoplasms Using Targeted Exon Resequencing of 86 Genes Involved in JAK-STAT Signaling and Epigenetic Regulation.

Myeloproliferative neoplasms (MPNs) are a heterogeneous group of blood disorders characterized by excess production of mature blood cells and an increased risk of late transformation to acute myeloid leukemia or primary myelofibrosis. Approximately 15% of MPN cases do not carry mutations in JAK2, CALR, or MPL and are thus often referred to as triple-negative cases. These are caused by a diverse set of rare mutations in cytokine receptors, JAK-STAT signaling pathway components, or epigenetic modifiers. In addition, some cases diagnosed as MPN are reactive rather than clonal disorders, so a negative result from a genetic screen can be informative. To obtain a comprehensive rapid molecular diagnosis for most MPNs, we developed an assay to detect genetic mutations (single nucleotide variants and/or small insertions/deletions) in 86 genes using targeted exon resequencing (AmpliSeq) and a bench-top semiconductor machine (Ion Torrent Personal Genome Machine). Our assay reliably detects well characterized mutations in JAK2, CALR, and MPL, but also rarer mutations in ASXL1, TET2, SH2B3, and other genes. Some of these mutations are novel. We find multiple mutations in advanced cases, suggesting co-operation between Janus kinase-STAT pathway mutations and epigenetic mutations in disease progression. This assay can be used to follow molecular progression, clonal heterogeneity, and drug resistance in MPNs.

Phase II Randomized Preoperative Window-of-Opportunity Study of the PI3K Inhibitor Pictilisib Plus Anastrozole Compared With Anastrozole Alone in Patients With Estrogen Receptor-Positive Breast Cancer.

PURPOSE: Preclinical data support a key role for the PI3K pathway in estrogen receptor-positive breast cancer and suggest that combining PI3K inhibitors with endocrine therapy may overcome resistance. This preoperative window study assessed whether adding the PI3K inhibitor pictilisib (GDC-0941) can increase the antitumor effects of anastrozole in primary breast cancer and aimed to identify the most appropriate patient population for combination therapy. PATIENTS AND METHODS: In this randomized, open-label phase II trial, postmenopausal women with newly diagnosed operable estrogen receptor-positive, human epidermal growth factor receptor 2 (HER2)-negative breast cancers were recruited. Participants were randomly allocated (2:1, favoring the combination) to 2 weeks of preoperative treatment with anastrozole 1 mg once per day (n = 26) or the combination of anastrozole 1 mg with pictilisib 260 mg once per day (n = 49). The primary end point was inhibition of tumor cell proliferation as measured by change in Ki-67 protein expression between tumor samples taken before and at the end of treatment. RESULTS: There was significantly greater geometric mean Ki-67 suppression of 83.8% (one-sided 95% CI, ≥ 79.0%) for the combination and 66.0% (95% CI, ≤ 75.4%) for anastrozole (geometric mean ratio [combination:anastrozole], 0.48; 95% CI, ≤ 0.72; P = .004). PIK3CA mutations were not predictive of response to pictilisib, but there was significant interaction between response to treatment and molecular subtype (P = .03); for patients with luminal B tumors, the combination:anastrozole geometric mean ratio of Ki-67 suppression was 0.37 (95% CI, ≤ 0.67; P = .008), whereas no significant Ki-67 response was observed for pictilisib in luminal A tumors (1.01; P = .98). Multivariable analysis confirmed Ki-67 response to the combination treatment of patients with luminal B tumors irrespective of progesterone receptor status or baseline Ki-67 expression. CONCLUSION: Adding pictilisib to anastrozole significantly increases suppression of tumor cell proliferation in luminal B primary breast cancer.

Germline polymorphisms in an enhancer of PSIP1 are associated with progression-free survival in epithelial ovarian cancer.

Women with epithelial ovarian cancer (EOC) are usually treated with platinum/taxane therapy after cytoreductive surgery but there is considerable inter-individual variation in response. To identify germline single-nucleotide polymorphisms (SNPs) that contribute to variations in individual responses to chemotherapy, we carried out a multi-phase genome-wide association study (GWAS) in 1,244 women diagnosed with serous EOC who were treated with the same first-line chemotherapy, carboplatin and paclitaxel. We identified two SNPs (rs7874043 and rs72700653) in TTC39B (best P=7x10-5, HR=1.90, for rs7874043) associated with progression-free survival (PFS). Functional analyses show that both SNPs lie in a putative regulatory element (PRE) that physically interacts with the promoters of PSIP1, CCDC171 and an alternative promoter of TTC39B. The C allele of rs7874043 is associated with poor PFS and showed increased binding of the Sp1 transcription factor, which is critical for chromatin interactions with PSIP1. Silencing of PSIP1 significantly impaired DNA damage-induced Rad51 nuclear foci and reduced cell viability in ovarian cancer lines. PSIP1 (PC4 and SFRS1 Interacting Protein 1) is known to protect cells from stress-induced apoptosis, and high expression is associated with poor PFS in EOC patients. We therefore suggest that the minor allele of rs7874043 confers poor PFS by increasing PSIP1 expression.

Exemplary multiplex bisulfite amplicon data used to demonstrate the utility of Methpat.

BACKGROUND: DNA methylation is a complex epigenetic marker that can be analyzed using a wide variety of methods. Interpretation and visualization of DNA methylation data can mask complexity in terms of methylation status at each CpG site, cellular heterogeneity of samples and allelic DNA methylation patterns within a given DNA strand. Bisulfite sequencing is considered the gold standard, but visualization of massively parallel sequencing results remains a significant challenge. FINDINGS: We created a program called Methpat that facilitates visualization and interpretation of bisulfite sequencing data generated by massively parallel sequencing. To demonstrate this, we performed multiplex PCR that targeted 48 regions of interest across 86 human samples. The regions selected included known gene promoters associated with cancer, repetitive elements, known imprinted regions and mitochondrial genomic sequences. We interrogated a range of samples including human cell lines, primary tumours and primary tissue samples. Methpat generates two forms of output: a tab-delimited text file for each sample that summarizes DNA methylation patterns and their read counts for each amplicon, and a HTML file that summarizes this data visually. Methpat can be used with publicly available whole genome bisulfite sequencing and reduced representation bisulfite sequencing datasets with sufficient read depths. CONCLUSIONS: Using Methpat, complex DNA methylation data derived from massively parallel sequencing can be summarized and visualized for biological interpretation. By accounting for allelic DNA methylation states and their abundance in a sample, Methpat can unmask the complexity of DNA methylation and yield further biological insight in existing datasets.

DNA methylation of oestrogen-regulated enhancers defines endocrine sensitivity in breast cancer.

Expression of oestrogen receptor (ESR1) determines whether a breast cancer patient receives endocrine therapy, but does not guarantee patient response. The molecular factors that define endocrine response in ESR1-positive breast cancer patients remain poorly understood. Here we characterize the DNA methylome of endocrine sensitivity and demonstrate the potential impact of differential DNA methylation on endocrine response in breast cancer. We show that DNA hypermethylation occurs predominantly at oestrogen-responsive enhancers and is associated with reduced ESR1 binding and decreased gene expression of key regulators of ESR1 activity, thus providing a novel mechanism by which endocrine response is abated in ESR1-positive breast cancers. Conversely, we delineate that ESR1-responsive enhancer hypomethylation is critical in transition from normal mammary epithelial cells to endocrine-responsive ESR1-positive cancer. Cumulatively, these novel insights highlight the potential of ESR1-responsive enhancer methylation to both predict ESR1-positive disease and stratify ESR1-positive breast cancer patients as responders to endocrine therapy.

Multiplex bisulfite PCR resequencing of clinical FFPE DNA.

BACKGROUND: The clinical utility of DNA methylation as a predictive or prognostic biomarker requires scalable resequencing protocols for bisulfite-converted DNA. Key features of any validation method should be adaptability for low- or high-throughput needs and high reproducibility, and should only require minimal amounts of precious clinical sample as input material. Critically, this method should also deliver robust results when working with bisulfite-converted DNA extracted from formalin-fixed, paraffin-embedded (FFPE) blocks. RESULTS: We report here for the first time on comparison studies between the Fluidigm Access Array system and multiplex assays for multiplex bisulfite PCR resequencing. The requirement of the Fluidigm Access Array system for high template amounts and its sensitivity to variations in template quality rendered it unsuitable for bisulfite PCR applications utilizing FFPE DNA. In response to this limitation, we established a multiplex bisulfite PCR assay capable of delivering robust methylation data using minimal amounts of FFPE clinical DNA. To evaluate the parameters and reproducibility of this assay, 57 amplicons were used to prepare sequencing libraries in triplicate for 13 FFPE tumour samples, as well as a series of 5 methylated controls (0%, 25%, 50%, 75%, and 100%). Analysis of this data demonstrated that this multiplex assay had high reproducibility (mean standard deviation of 1.4% for methylation values), was low cost, required low sample input (50 ng of DNA or less), and could be scaled for both low- and high-throughput needs. Notably, ExoSAP-IT (exonuclease I) treatment to remove residual primers in bisulfite resequencing libraries appeared to degrade the library and generate a high-molecular weight smear which may impact on the degree of methylation assessed. CONCLUSIONS: Multiplex bisulfite PCR assays represent a convenient and scalable method for validation and screening of methylated DNA regions from archival FFPE DNA. Moreover, the library construction process detailed here can be rapidly optimized and implemented with a minimal amount of work, can be performed using the standard equipment found in any molecular biology laboratory, and can be easily adapted for use on both genomic DNA and bisulfite DNA applications. However, in preparing bisulfite libraries for sequencing, the use of ExoSAP-IT is not recommended due to potential off-target nuclease effects which may impact downstream methylation analysis.