Clinical Implications of Monitoring ESR1 Mutations by Circulating Tumor DNA in Estrogen Receptor Positive Metastatic Breast Cancer: A Pilot Study

BACKGROUND: ESR1 mutations are frequently detected in ER+ MBC, and have been reported to be associated with endocrine therapy resistance. However, there are little researches to validate whether dynamic monitoring of ESR1 mutations could serve as a predictive plasma biomarker of acquired resistance to endocrine therapy. Therefore, in this study, we performed longitudinal circulating tumor DNA (ctDNA) detection to evaluate the clinical implications of monitoring ESR1 mutations. METHODS: We performed longitudinal dynamic mutation analyses of plasma samples from 45 patients with metastatic breast cancer (MBC) and sequencing paired biopsy tissues, using a targeted NGS panel of 425 genes. These patients were treated at the Second Affiliated Hospital of Dalian Medical University between January 2017 and February 2019 with written informed consent. RESULTS: Mutations profiles were highly concordant between plasma and paired tissue samples from 45 MBC patients (r = 0.96, P < 0.0001). ESR1 mutations were enriched in ER+ MBC patients after AI therapy (17.8%, 8/45). The median time from AI endocrine therapies to the initial detection of ESR1 mutation was 39 months (95% CI 21.32–57.57). Some hotspot mutations (Y537S (n = 5), Y537N (n = 1), D538G (n = 2), E380Q (n = 2)) and several rare mutations (L345SfsX7, 24fs, G344delinsGC) were identified in our cohort. In addition, we observed that two patients obtained multiple ESR1 mutations over the course of treatment (Y537N/Y537S/D538G, L345SfsX7/24fs/E380Q). Through dynamically monitoring ESR1 mutations by ctDNA, we demonstrated that the change of allele frequency of ESR1 mutations was an important biomarker, which could predict endocrine resistance of ER+ MBC in our study. We also observed that the combination of everolimus in four cases with acquired ESR1 mutations showed longer PFS than other therapies without everolimus. CONCLUSION: The dynamic monitoring of ESR1 mutations by ctDNA is a promising tool to predict endocrine therapy resistance in ER+ MBC patients.     

Detection of ctDNA with Personalized Molecular Barcode NGS and Its Clinical Significance in Patients with Early Breast Cancer

We attempted to detect circulating tumor DNA (ctDNA), taking advantage of molecular barcode next-generation sequencing (MB-NGS), which can be more easily customized to detect a variety of mutations with a high sensitivity than PCR-based methods. Sequencing with a gene panel consisting of the 13 most frequently mutated genes in breast tumors from stage I or II patients revealed 95 somatic mutations in the 12 genes in 62% (62/100) of tumors. Then, plasma DNA from each patient (n = 62) before surgery was analyzed via MB-NGS customized to each somatic mutation, resulting in the detection of ctDNA in 16.1% (10/62) of patients. ctDNA was significantly associated with biologically aggressive phenotypes, including large tumor size (P = .004), positive lymph node (P = .009), high histological grade (P < .001), negative ER (P = .018), negative PR (P = .017), and positive HER2 (P = .046). Furthermore, distant disease-free survival was significantly worse in patients with ctDNA (n = 10) than those without ctDNA (n = 52) (P < .001). Our results demonstrate that MB-NGS personalized to each mutation can detect ctDNA with a high sensitivity in early breast cancer patients at diagnosis, and it seems to have a potential to serve as a clinically useful tumor marker for predicting their prognosis.     

Exome Sequencing of Cell-Free DNA from Metastatic Cancer Patients Identifies Clinically Actionable Mutations Distinct from Primary Disease

The identification of the molecular drivers of cancer by sequencing is the backbone of precision medicine and the basis of personalized therapy; however, biopsies of primary tumors provide only a snapshot of the evolution of the disease and may miss potential therapeutic targets, especially in the metastatic setting. A liquid biopsy, in the form of cell-free DNA (cfDNA) sequencing, has the potential to capture the inter- and intra-tumoral heterogeneity present in metastatic disease, and, through serial blood draws, track the evolution of the tumor genome.In order to determine the clinical utility of cfDNA sequencing we performed whole-exome sequencing on cfDNA and tumor DNA from two patients with metastatic disease; only minor modifications to our sequencing and analysis pipelines were required for sequencing and mutation calling of cfDNA. The first patient had metastatic sarcoma and 47 of 48 mutations present in the primary tumor were also found in the cell-free DNA. The second patient had metastatic breast cancer and sequencing identified an ESR1 mutation in the cfDNA and metastatic site, but not in the primary tumor. This likely explains tumor progression on Anastrozole. Significant heterogeneity between the primary and metastatic tumors, with cfDNA reflecting the metastases, suggested separation from the primary lesion early in tumor evolution. This is best illustrated by an activating PIK3CA mutation (H1047R) which was clonal in the primary tumor, but completely absent from either the metastasis or cfDNA. Here we show that cfDNA sequencing supplies clinically actionable information with minimal risks compared to metastatic biopsies. This study demonstrates the utility of whole-exome sequencing of cell-free DNA from patients with metastatic disease. cfDNA sequencing identified an ESR1 mutation, potentially explaining a patient’s resistance to aromatase inhibition, and gave insight into how metastatic lesions differ from the primary tumor.     

Comparison of ESR1 Mutations in Tumor Tissue and Matched Plasma Samples from Metastatic Breast Cancer Patients

BACKGROUND: ESR1 mutation in circulating cell-free DNA (cfDNA) is emerging as a noninvasive biomarker of acquired resistance to endocrine therapy, but there is a paucity of data comparing the status of ESR1 gene in cfDNA with that in its corresponding tumor tissue. The objective of this study is to validate the degree of concordance of ESR1 mutations between plasma and tumor tissue. METHODS: ESR1 ligand-binding domain mutations Y537S, Y537N, Y537C, and D538G were analyzed using droplet digital PCR in 35 patients with metastatic breast cancer (MBC) (35 tumor tissue samples and 67 plasma samples). RESULTS: Of the 35 paired samples, 26 (74.3%) were concordant: one patient had detectable ESR1 mutations both plasma (ESR1 Y537S/Y537N) and tumor tissue (ESR1 Y537S/Y537C), and 25 had WT ESR1 alleles in both. Nine (25.7%) had discordance between the plasma and tissue results: five had mutations detected only in their tumor tissue (two Y537S, one Y537C, one D538G, and one Y537S/Y537N/D538G), and four had mutations detected only in their plasma (one Y537S, one Y537N, and two Y537S/Y537N/D538G). Furthermore, longitudinal plasma samples from 19 patients were used to assess changes in the presence of ESR1 mutations during treatment. Eleven patients had cfDNA ESR1 mutations over the course of treatment. A total of eight of 11 patients with MBC with cfDNA ESR1 mutations (72.7%) had the polyclonal mutations. CONCLUSION: We have shown the independent distribution of ESR1 mutations between plasma and tumor tissue in 35 patients with MBC.     

Extended RAS and BRAF Mutation Analysis Using Next-Generation Sequencing

Somatic mutations in KRAS, NRAS, and BRAF genes are related to resistance to anti-EGFR antibodies in colorectal cancer. We have established an extended RAS and BRAF mutation assay using a next-generation sequencer to analyze these mutations. Multiplexed deep sequencing was performed to detect somatic mutations within KRAS, NRAS, and BRAF, including minor mutated components. We first validated the technical performance of the multiplexed deep sequencing using 10 normal DNA and 20 formalin-fixed, paraffin-embedded (FFPE) tumor samples. To demonstrate the potential clinical utility of our assay, we profiled 100 FFPE tumor samples and 15 plasma samples obtained from colorectal cancer patients. We used a variant calling approach based on a Poisson distribution. The distribution of the mutation-positive population was hypothesized to follow a Poisson distribution, and a mutation-positive status was defined as a value greater than the significance level of the error rate (α = 2 x 10^-5). The cut-off value was determined to be the average error rate plus 7 standard deviations. Mutation analysis of 100 clinical FFPE tumor specimens was performed without any invalid cases. Mutations were detected at a frequency of 59% (59/100). KRAS mutation concordance between this assay and Scorpion-ARMS was 92% (92/100). DNA obtained from 15 plasma samples was also analyzed. KRAS and BRAF mutations were identified in both the plasma and tissue samples of 6 patients. The genetic screening assay using next-generation sequencer was validated for the detection of clinically relevant RAS and BRAF mutations using FFPE and liquid samples.     

Multi-Purpose Utility of Circulating Plasma DNA Testing in Patients with Advanced Cancers

Tumor genomic instability and selective treatment pressures result in clonal disease evolution; molecular stratification for molecularly targeted drug administration requires repeated access to tumor DNA. We hypothesized that circulating plasma DNA (cpDNA) in advanced cancer patients is largely derived from tumor, has prognostic utility, and can be utilized for multiplex tumor mutation sequencing when repeat biopsy is not feasible. We utilized the Sequenom MassArray System and OncoCarta panel for somatic mutation profiling. Matched samples, acquired from the same patient but at different time points were evaluated; these comprised formalin-fixed paraffin-embedded (FFPE) archival tumor tissue (primary and/or metastatic) and cpDNA. The feasibility, sensitivity, and specificity of this high-throughput, multiplex mutation detection approach was tested utilizing specimens acquired from 105 patients with solid tumors referred for participation in Phase I trials of molecularly targeted drugs. The median cpDNA concentration was 17 ng/ml (range: 0.5–1600); this was 3-fold higher than in healthy volunteers. Moreover, higher cpDNA concentrations associated with worse overall survival; there was an overall survival (OS) hazard ratio of 2.4 (95% CI 1.4, 4.2) for each 10-fold increase in cpDNA concentration and in multivariate analyses, cpDNA concentration, albumin, and performance status remained independent predictors of OS. These data suggest that plasma DNA in these cancer patients is largely derived from tumor. We also observed high detection concordance for critical ‘hot-spot’ mutations (KRAS, BRAF, PIK3CA) in matched cpDNA and archival tumor tissue, and important differences between archival tumor and cpDNA. This multiplex sequencing assay can be utilized to detect somatic mutations from plasma in advanced cancer patients, when safe repeat tumor biopsy is not feasible and genomic analysis of archival tumor is deemed insufficient. Overall, circulating nucleic acid biomarker studies have clinically important multi-purpose utility in advanced cancer patients and further studies to pursue their incorporation into the standard of care are warranted.     

Oligonucleotide gap-fill ligation for mutation detection and sequencing in situ

In clinical diagnostics a great need exists for targeted in situ multiplex nucleic acid analysis as the mutational status can offer guidance for effective treatment. One well-established method uses padlock probes for mutation detection and multiplex expression analysis directly in cells and tissues. Here, we use oligonucleotide gap-fill ligation to further increase specificity and to capture molecular substrates for in situ sequencing. Short oligonucleotides are joined at both ends of a padlock gap probe by two ligation events and are then locally amplified by target-primed rolling circle amplification (RCA) preserving spatial information. We demonstrate the specific detection of the A3243G mutation of mitochondrial DNA and we successfully characterize a single nucleotide variant in the ACTB mRNA in cells by in situ sequencing of RCA products generated by padlock gap-fill ligation. To demonstrate the clinical applicability of our assay, we show specific detection of a point mutation in the EGFR gene in fresh frozen and formalin-fixed, paraffin-embedded (FFPE) lung cancer samples and confirm the detected mutation by in situ sequencing. This approach presents several advantages over conventional padlock probes allowing simpler assay design for multiplexed mutation detection to screen for the presence of mutations in clinically relevant mutational hotspots directly in situ.     

Development and Validation of a Next-Generation Sequencing Assay for BRCA1 and BRCA2 Variants for the Clinical Laboratory

The objective of this study was to design and validate a next-generation sequencing assay (NGS) to detect BRCA1 and BRCA2 mutations. We developed an assay using random shearing of genomic DNA followed by RNA bait tile hybridization and NGS sequencing on both the Illumina MiSeq and Ion Personal Gene Machine (PGM). We determined that the MiSeq Reporter software supplied with the instrument could not detect deletions greater than 9 base pairs. Therefore, we developed an alternative alignment and variant calling software, Quest Sequencing Analysis Pipeline (QSAP), that was capable of detecting large deletions and insertions. In validation studies, we used DNA from 27 stem cell lines, all with known deleterious BRCA1 or BRCA2 mutations, and DNA from 67 consented control individuals who had a total of 352 benign variants. Both the MiSeq/QSAP combination and PGM/Torrent Suite combination had 100% sensitivity for the 379 known variants in the validation series. However, the PGM/Torrent Suite combination had a lower intra- and inter-assay precision of 96.2% and 96.7%, respectively when compared to the MiSeq/QSAP combination of 100% and 99.4%, respectively. All PGM/Torrent Suite inconsistencies were false-positive variant assignments. We began commercial testing using both platforms and in the first 521 clinical samples MiSeq/QSAP had 100% sensitivity for BRCA1/2 variants, including a 64-bp deletion and a 10-bp insertion not identified by PGM/Torrent Suite, which also suffered from a high false-positive rate. Neither the MiSeq nor PGM platform with their supplied alignment and variant calling software are appropriate for a clinical laboratory BRCA sequencing test. We have developed an NGS BRCA1/2 sequencing assay, MiSeq/QSAP, with 100% analytic sensitivity and specificity in the validation set consisting of 379 variants. The MiSeq/QSAP combination has sufficient performance for use in a clinical laboratory.     

Screening of Duchenne Muscular Dystrophy (DMD) Mutations and Investigating Its Mutational Mechanism in Chinese Patients

Duchenne muscular dystrophy (DMD) is a common X-linked recessive disease of muscle degeneration and death. In order to provide accurate and reliable genetic counseling and prenatal diagnosis, we screened DMD mutations in a cohort of 119 Chinese patients using multiplex ligation-dependent probe amplification (MLPA) and denaturing high performance liquid chromatography (DHPLC) followed by Sanger sequencing. In these unrelated DMD patients, we identified 11 patients with DMD small mutations (9.2%) and 81 patients with DMD deletions/duplications (del/dup) (68.1%), of which 64 (79.0%) were deletions, 16 (19.8%) were duplications, and one (1.2%) was both deletion and duplication. Furthermore, we analyzed the frequency of DMD breakpoint in the 64 deletion cases by calculating exon-deletion events of certain exon interval that revealed a novel mutation hotspot boundary. To explore why DMD rearrangement breakpoints were predisposed to specific regions (hotspot), we precisely characterized junction sequences of breakpoints at the nucleotide level in 21 patients with exon deleted/duplicated in DMD with a high-resolution SNP microarray assay. There were no exactly recurrent breakpoints and there was also no significant difference between single-exon del/dup and multiple-exon del/dup cases. The data from the current study provided a comprehensive strategy to detect DMD mutations for clinical practice, and identified two deletion hotspots at exon 43–55 and exon 10–23 by calculating exon-deletion events of certain exon interval. Furthermore, this is the first study to characterize DMD breakpoint at the nucleotide level in a Chinese population. Our observations provide better understanding of the mechanism for DMD gene rearrangements.     

Microsatellites in the estrogen receptor (ESR1, ESR2) and androgen receptor (AR) genes and breast cancer risk in African American and Nigerian women

Genetic variants in hormone receptor genes may be crucial predisposing factors for breast cancer, and microsatellites in the estrogen receptor (ESR1, ESR2) and androgen receptor (AR) genes have been suggested to play a role. We studied 258 African-American (AA) women with breast cancer and 259 hospital-based controls, as well as 349 Nigerian (NG) female breast cancer patients and 296 community controls. Three microsatellites, ESR1_TA, ESR2_CA and AR_CAG, in the ESR1, ESR2 and AR genes, respectively, were genotyped. Their repeat lengths were then analyzed as continuous and dichotomous variables. Analyses of continuous variables showed no association with breast cancer risk in either AA or NG at ESR1_TA; AA cases had shorter repeats in the long allele of ESR2_CA than AA controls (Mann-Whitney P  = 0.036; logistic regression P  = 0.04, OR  = 0.91, 95% CI 0.83–1.00), whereas NG patients had longer repeats in the short allele than NG controls (Mann-Whitney P  = 0.0018; logistic regression P  = 0.04, OR  = 1.06, 95% CI 1.00–1.11); and AA cases carried longer repeats in the short allele of AR_CAG than AA controls (Mann-Whitney P  = 0.038; logistic regression P  = 0.03, OR  = 1.08, 95% CI 1.01–1.15). When allele sizes were categorized as dichotomous variables, we discovered that women with two long alleles of ESR2_CA had increased risk of breast cancer (OR  = 1.38, 95% CI 1.10–1.74; P  = 0.006). This is the first study to investigate these three microsatellites in hormonal receptor genes in relation to breast cancer risk in an indigenous African population. After adjusting for multiple-testing, our findings suggest that ESR2_CA is associated with breast cancer risk in Nigerian women, whereas ESR1_TA and AR_CAG seem to have no association with the disease among African American or Nigerian women.     

Efficient Genotyping of KRAS Mutant Non-Small Cell Lung Cancer Using a Multiplexed Droplet Digital PCR Approach

Droplet digital PCR (ddPCR) can be used to detect low frequency mutations in oncogene-driven lung cancer. The range of KRAS point mutations observed in NSCLC necessitates a multiplex approach to efficient mutation detection in circulating DNA. Here we report the design and optimisation of three discriminatory ddPCR multiplex assays investigating nine different KRAS mutations using PrimePCR™ ddPCR™ Mutation Assays and the Bio-Rad QX100 system. Together these mutations account for 95% of the nucleotide changes found in KRAS in human cancer. Multiplex reactions were optimised on genomic DNA extracted from KRAS mutant cell lines and tested on DNA extracted from fixed tumour tissue from a cohort of lung cancer patients without prior knowledge of the specific KRAS genotype. The multiplex ddPCR assays had a limit of detection of better than 1 mutant KRAS molecule in 2,000 wild-type KRAS molecules, which compared favourably with a limit of detection of 1 in 50 for next generation sequencing and 1 in 10 for Sanger sequencing. Multiplex ddPCR assays thus provide a highly efficient methodology to identify KRAS mutations in lung adenocarcinoma.     

Downsloping High-Frequency Hearing Loss Due to Inner Ear Tricellular Tight Junction Disruption by a Novel ILDR1 Mutation in the Ig-Like Domain

The immunoglobulin (Ig)-like domain containing receptor 1 (ILDR1) gene encodes angulin-2/ILDR1, a recently discovered tight junction protein, which forms tricellular tight junction (tTJ) structures with tricellulin and lipolysis-stimulated lipoprotein receptor (LSR) at tricellular contacts (TCs) in the inner ear. Previously reported recessive mutations within ILDR1 have been shown to cause severe to profound nonsyndromic sensorineural hearing loss (SNHL), DFNB42. Whole-exome sequencing of a Korean multiplex family segregating partial deafness identified a novel homozygous ILDR1 variant (p.P69H) within the Ig-like domain. To address the pathogenicity of p.P69H, the angulin-2/ILDR1 p.P69H variant protein, along with the previously reported pathogenic ILDR1 mutations, was expressed in angulin-1/LSR knockdown epithelial cells. Interestingly, partial mislocalization of the p.P69H variant protein and tricellulin at TCs was observed, in contrast to a severe mislocalization and complete failure of tricellulin recruitment of the other reported ILDR1 mutations. Additionally, three-dimensional protein modeling revealed that angulin-2/ILDR1 contributed to tTJ by forming a homo-trimer structure through its Ig-like domain, and the p.P69H variant was predicted to disturb homo-trimer formation. In this study, we propose a possible role of angulin-2/ILDR1 in tTJ formation in the inner ear and a wider audiologic phenotypic spectrum of DFNB42 caused by mutations within ILDR1.     

RASSF1A Promoter Methylation Levels Positively Correlate with Estrogen Receptor Expression in Breast Cancer Patients

The aim of this study was to investigate the relationship between the promoter methylation in five cancer-associated genes and clinicopathologic features for identification of molecular markers of tumor metastatic potential and hormone therapy response efficiency in breast cancer. The methylation levels in paraffin-embedded tumor tissues, plasma, and blood cells from 151 sporadic breast cancer patients and blood samples of 50 controls were evaluated by quantitative multiplex methylation-specific polymerase chain reaction. DNA methylation of RAS-association domain family member 1 (RASSF1A), estrogen receptor 1 (ESR1), cadherin 1, type 1, E-cadherin (CDH1), TIMP metallopeptidase inhibitor 3 (TIMP3) and spleen tyrosine kinase (SYK) genes was detected in the tumors of 124, 19, 15, 15, and 6 patients with mean levels of 48.45%, 3.81%, 2.36%, 27.55%, and 10.81%, respectively. Plasma samples exhibited methylation in the same genes in 25, 10, 15, 17, and 3 patients with levels of 22.54%, 17.20%, 22.87%, 31.93%, and 27.42%, respectively. Cumulative methylation results confirmed different spectra in tumor and plasma samples. Simultaneous methylation in tumors and plasma were shown in less than 17% of patients. RASSF1A methylation levels in tumor samples statistically differ according to tumor size (P = .029), estrogen receptor (ER) and progesterone receptor (PR) status (P = .000 and P = .004), and immunohistochemical subtype (P = .000). Moreover, the positive correlation was found between RASSF1A methylation levels and percentage of cancer cells expressing ER and PR. The direct relationship between RASSF1A promoter methylation and expression of ER could aid the prognosis of hormonal therapy response.     

ESR1 Amplification in Breast Cancer by Optimized RNase FISH: Frequent but Low-Level and Heterogeneous

Prevalence of ESR1 amplification in breast cancer is highly disputed and discrepancies have been related to different technical protocols and different scoring approaches. In addition, pre-mRNA artifacts have been proposed to influence outcome of ESR1 FISH analysis. We analyzed ESR1 gene copy number status combining an improved RNase FISH protocol with multiplex ligation-dependent probe amplification (MLPA) after laser microdissection. FISH showed a high prevalence of ESR1 gains and amplifications despite RNase treatment but MLPA did not confirm ESR1 copy number increases detected by FISH in more than half of cases. We suggest that the combination of the ESR1-specific intra-tumor heterogeneity and low-level copy number increase accounts for these discrepancies.     

Analysis of PALB2 Gene in BRCA1/BRCA2 Negative Spanish Hereditary Breast/Ovarian Cancer Families with Pancreatic Cancer Cases

Background

The PALB2 gene, also known as FANCN, forms a bond and co-localizes with BRCA2 in DNA repair. Germline mutations in PALB2 have been identified in approximately 1% of familial breast cancer and 3–4% of familial pancreatic cancer. The goal of this study was to determine the prevalence of PALB2 mutations in a population of BRCA1/BRCA2 negative breast cancer patients selected from either a personal or family history of pancreatic cancer.

Methods

132 non-BRCA1/BRCA2 breast/ovarian cancer families with at least one pancreatic cancer case were included in the study. PALB2 mutational analysis was performed by direct sequencing of all coding exons and intron/exon boundaries, as well as multiplex ligation-dependent probe amplification.

Results

Two PALB2 truncating mutations, the c.1653T>A (p.Tyr551Stop) previously reported, and c.3362del (p.Gly1121ValfsX3) which is a novel frameshift mutation, were identified. Moreover, several PALB2 variants were detected; some of them were predicted as pathological by bioinformatic analysis. Considering truncating mutations, the prevalence rate of our population of BRCA1/2-negative breast cancer patients with pancreatic cancer is 1.5%.

Conclusions

The prevalence rate of PALB2 mutations in non-BRCA1/BRCA2 breast/ovarian cancer families, selected from either a personal or family pancreatic cancer history, is similar to that previously described for unselected breast/ovarian cancer families. Future research directed towards identifying other gene(s) involved in the development of breast/pancreatic cancer families is required.