Testing clonal relatedness of tumors using array comparative genomic hybridization: a statistical challenge

Chromosome 2q gain and epigenetic silencing of GATA3 in microglandular adenosis of the breast

Microglandular adenosis (MGA) represents a rare neoplasm of the mammary gland, which in a subset of cases may be associated with triple-negative breast cancer (BC). The biology of MGA is poorly understood. In this study, eight MGA cases (n = 4 with and n = 4 without associated BC) were subjected to a comprehensive characterization using immunohistochemistry, genome-wide DNA copy number (CN) profiling, fluorescence in situ hybridization (FISH), next-generation sequencing (NGS), and DNA methylation profiling using 850 K arrays and bisulfite pyrosequencing. Median patient age was 61 years (range 57-76 years). MGA lesions were estrogen receptor (ER)-negative, progesterone receptor-negative, HER2-negative, and S100-positive. DNA CN alterations (CNAs) were complex or limited to few gains and losses. CN gain on chromosome 2q was the most common CNA and was validated by FISH in five of eight cases. NGS demonstrated an average of two mutations per case (range 0-5) affecting 10 different genes (ARID1A, ATM, CTNNB1, FBXW7, FGFR2, MET, PIK3CA, PMS2, PTEN, and TP53). CNAs and mutations were similar in MGA and adjacent BC, indicating clonal relatedness. DNA methylation profiling identified aberrant hypermethylation of CpG sites within GATA3, a key transcription factor required for luminal differentiation. Immunohistochemistry showed regular GATA3 protein expression in the normal mammary epithelium and in ER-positive BC. Conversely, GATA3 was reduced or lost in all MGA cases tested (8/8). In conclusion, MGA is characterized by common CN gain on chromosome 2q and loss of GATA3. Epigenetic inactivation of GATA3 may provide a new clue to the peculiar biology of this rare neoplasia.

E-cadherin to P-cadherin switching in lobular breast cancer with tubular elements

Loss of E-cadherin expression due to mutation of the CDH1 gene is a characteristic feature of invasive lobular breast cancer (ILBC). Beta-catenin, which binds to the cytoplasmic domain of E-cadherin, is simultaneously downregulated, reflecting disassembly of adherens junctions (AJs) and loss of cell adhesion. E-cadherin to P-cadherin expression switching can rescue AJs and cell adhesion. However, P-cadherin has not been implicated in ILBC, so far. We aimed to characterize 13 ILBCs with exceptional histomorphology, which we termed ILBCs with tubular elements. The CDH1 mutational status was determined by next generation sequencing and whole-genome copy number (CN) profiling. Expression of cadherins was assessed by immunohistochemistry. ILBCs with tubular elements were ER-positive (13/13) and HER2-negative (13/13) and harbored deleterious CDH1 mutations (11/13) accompanied by loss of heterozygosity due to deletion of chromosome 16q22.1 (9/11). E-cadherin expression was lost or reduced in noncohesive tumor cells and in admixed tubular elements (13/13). Beta-catenin expression was lost in noncohesive tumor cells, but was retained in tubular elements (11/13), indicating focal rescue of AJ formation. N-cadherin and R-cadherin were always negative (0/13). Strikingly, P-cadherin was commonly positive (12/13) and immunoreactivity was accentuated in tubular elements. Adjacent lobular carcinoma in situ (LCIS) was always P-cadherin-negative (0/7). In a reference cohort of LCIS specimens, P-cadherin was constantly not expressed (0/25). In a reference cohort of invasive mammary carcinomas, P-cadherin-positive cases (36/268, 13%) were associated with triple-negative nonlobular breast cancer (P < 0.001). Compared with ILBCs from the reference cohort, P-cadherin expression was more common in ILBCs with tubular elements (12/13 versus 7/84, P < 0.001). In summary, E-cadherin to P-cadherin switching occurs in a subset of ILBCs. P-cadherin is the molecular determinant of a mixed-appearing histomorphology in ILBCs with tubular elements.

Atypical ductal hyperplasia is a multipotent precursor of breast carcinoma

The current model for breast cancer progression proposes independent 'low grade (LG)-like' and 'high grade (HG)-like' pathways but lacks a known precursor to HG cancer. We applied low-coverage whole-genome sequencing to atypical ductal hyperplasia (ADH) with and without carcinoma to shed light on breast cancer progression. Fourteen out of twenty isolated ADH cases harboured at least one copy number alteration (CNA), but had fewer aberrations than LG or HG ductal carcinoma in situ (DCIS). ADH carried more HG-like CNA than LG DCIS (e.g. 8q gain). Correspondingly, 64% (7/11) of ADH cases with synchronous HG carcinoma were clonally related, similar to LG carcinoma (67%, 6/9). This study represents a significant shift in our understanding of breast cancer progression, with ADH as a common precursor lesion to the independent 'low grade-like' and 'high grade-like' pathways. These data suggest that ADH can be a precursor of HG breast cancer and that LG and HG carcinomas can evolve from a similar ancestor lesion. We propose that although LG DCIS may be committed to a LG molecular pathway, ADH may remain multipotent, progressing to either LG or HG carcinoma. This multipotent nature suggests that some ADH cases could be more clinically significant than LG DCIS, requiring biomarkers for personalising management. Copyright © 2019 Pathological Society of Great Britain and Ireland. Published by John Wiley & Sons, Ltd.

Clonal relatedness in tumour pairs of breast cancer patients

Background

Molecular classification of tumour clonality is currently not evaluated in multiple invasive breast carcinomas, despite evidence suggesting common clonal origins. There is no consensus about which type of data (e.g. copy number, mutation, histology) and especially which statistical method is most suitable to distinguish clonal recurrences from independent primary tumours.

Methods

Thirty-seven invasive breast tumour pairs were stratified according to laterality and time interval between the diagnoses of the two tumours. In a multi-omics approach, tumour clonality was analysed by integrating clinical characteristics (n = 37), DNA copy number (n = 37), DNA methylation (n = 8), gene expression microarray (n = 7), RNA sequencing (n = 3), and SNP genotyping data (n = 3). Different statistical methods, e.g. the diagnostic similarity index (SI), were used to classify the tumours as clonally related recurrences or independent primary tumours.

Results

The SI and hierarchical clustering showed similar tendencies and the highest concordance with the other methods. Concordant evidence for tumour clonality was found in 46% (17/37) of patients. Notably, no association was found between the current clinical guidelines and molecular tumour features.

Conclusions

A more accurate classification of clonal relatedness between multiple breast tumours may help to mitigate treatment failure and relapse by integrating tumour-associated molecular features, clinical parameters, and statistical methods. Guidelines need to be defined with exact thresholds to standardise clonality testing in a routine diagnostic setting.

Electronic supplementary material

The online version of this article (10.1186/s13058-018-1022-y) contains supplementary material, which is available to authorized users.

Clonality assessment of multifocal lung adenocarcinoma by pathology evaluation and molecular analysis

The aim of this study was to explore morphologic and molecular features distinguishing between multifocal lung adenocarcinoma (MLA) and intrapulmonary metastases (IMs). Sixteen patients with MLAs, a total of 34 tumors, were reviewed. Four approaches were used: (1) array-comparative genomic hybridization (CGH) as a standard clonality assessment; (2) EGFR and KRAS mutational profiles as a supplementary method; (3) comprehensive histologic assessment (CHA) was method I in pathology evaluation; and (4) CHA combined with lepidic component analysis was method II. The lepidic component was divided into low grade and high grade according to extent of atypia; tumors with low-grade lepidic component were defined as primary. Eight patients were found to have IMs and 8 to have multiple primaries (MPs) by array-CGH; 7 had MPs and 9 had IMs by method I; 5 had MPs and 11 had IMs by method II. Compared with array-CGH, method I had a lower coincidence rate (65%) than method II (85%). Univariate analysis revealed that patients with MP had a better clinical outcome than those with IM only if the MPs were diagnosed by array-CGH (P = .034) or method II (P = .027) but not EGFR/KRAS mutation (P = .843) or method I (P = .493). Our results suggest that a low-grade lepidic component is a sign of a primary tumor. CHA combined with a low-grade lepidic component (method II) is more accurate clinically and more cost-effective in distinguishing MLAs from IMs. Also, EGFR mutation is not an appropriate molecular marker for clonality assessment.

Contralateral breast cancers: Independent cancers or metastases?

A cancer in the contralateral breast in a woman with a previous or synchronous breast cancer is typically considered to be an independent primary tumor. Emerging evidence suggests that in a small subset of these cases the second tumor represents a metastasis. We sought to investigate the issue using massively parallel sequencing targeting 254 genes recurrently mutated in breast cancer. We examined the tumor archives at Memorial Sloan Kettering Cancer Center for the period 1995-2006 to identify cases of contralateral breast cancer where surgery for both tumors was performed at the Center. We report results from 49 patients successfully analyzed by a targeted massively parallel sequencing assay. Somatic mutations and copy number alterations were defined by state-of-the-art algorithms. Clonal relatedness was evaluated by statistical tests specifically designed for this purpose. We found evidence that the tumors in contralateral breasts were clonally related in three cases (6%) on the basis of matching mutations at codons where somatic mutations are rare. Clinical data and the presence of similar patterns of gene copy number alterations were consistent with metastasis for all three cases. In three additional cases, there was a solitary matching mutation at a common PIK3CA locus. The results suggest that a subset of contralateral breast cancers represent metastases rather than independent primary tumors. Massively parallel sequencing analysis can provide important evidence to clarify the diagnosis. However, given the inter-tumor mutational heterogeneity in breast cancer, sufficiently large gene panels need to be employed to define clonality convincingly in all cases.

Subclonal analysis in a lobular breast cancer with classical and solid growth pattern mimicking a solid-papillary carcinoma

Recently, a new variant of invasive lobular breast cancer (ILBC) with solid-papillary-like growth pattern has been described. We present a case of ILBC with solid-papillary-like growth pattern in the main tumour mass and classical invasive lobular growth pattern in adjacent satellite foci. The two tumour components were subjected to comprehensive molecular analyses. Both components were ER/PR-positive, HER2-negative, and showed a complete loss of E-cadherin and beta-catenin protein expression, as determined by immunohistochemistry. Gene expression profiling classified the main tumour and a satellite focus as luminal-B and luminal-A subtypes, respectively. Whole-genome copy number profiles were highly similar in both tumour components. Shared copy number alterations (CNAs) included gains of chromosome 1q21.1-q43 and losses of chromosome 16q11.2-q24.3, the locus of the CDH1/E-cadherin tumour suppressor gene. CNAs detected only in the main tumour included a gain of chromosome 20q12-q13.33 and a loss of chromosome 1p36.33-p34.3, which has recently been associated with the solid variant of ILBC. Next generation sequencing revealed an identical, truncating CDH1 mutation (p.G169fs*5) in both tumour components confirming a common clonal ancestry. In conclusion, we confirm the recently described variant of ILBC with solid-papillary-like growth pattern and provide evidence that it evolves from classical ILBC by subclonal evolution.

Estimating the probability of clonal relatedness of pairs of tumors in cancer patients

Next generation sequencing panels are being used increasingly in cancer research to study tumor evolution. A specific statistical challenge is to compare the mutational profiles in different tumors from a patient to determine the strength of evidence that the tumors are clonally related, that is, derived from a single, founder clonal cell. The presence of identical mutations in each tumor provides evidence of clonal relatedness, although the strength of evidence from a match is related to how commonly the mutation is seen in the tumor type under investigation. This evidence must be weighed against the evidence in favor of independent tumors from non-matching mutations. In this article, we frame this challenge in the context of diagnosis using a novel random effects model. In this way, by analyzing a set of tumor pairs, we can estimate the proportion of cases that are clonally related in the sample as well as the individual diagnostic probabilities for each case. The method is illustrated using data from a study to determine the clonal relationship of lobular carcinoma in situ with subsequent invasive breast cancers, where each tumor in the pair was subjected to whole exome sequencing. The statistical properties of the method are evaluated using simulations, demonstrating that the key model parameters are estimated with only modest bias in small samples in most configurations.

Pathologic Staging and Survival of Patients With Synchronous Bilateral Lung Carcinomas

OBJECTIVES

To compare survival data in patients with resected bilateral synchronous pulmonary carcinomas with survival data from patients with lung cancer in pStages I through IV and to evaluate the usefulness of comprehensive histologic evaluation (CHE) of tumor histologic patterns to distinguish between synchronous primaries and intrapulmonary metastases.

METHODS

Ten-year overall survival (OS) data from 18 patients with 44 resected synchronous bilateral lung cancers, classified as "synchronous primaries" or "metastases" using CHE, were compared with survival data of 2,879 patients with lung cancer in pStages I through IV.

RESULTS

Forty and four tumors from 16 and two patients, respectively, were classified as synchronous primaries and metastases. There were no significant differences in survival between these 18 patients and pStage I controls or between the synchronous primaries and the metastases patient groups. However, there were significant differences in survival between the patients with resected synchronous bilateral tumors and each of the pStage II through IV control groups (P < .05).

CONCLUSIONS

Patients with resected synchronous bilateral lung cancers had similar 10-year OS to patients with stage I disease, regardless of CHE data. Most resected tumors were synchronous primaries by CHE.

Mutational Profiling Can Establish Clonal or Independent Origin in Synchronous Bilateral Breast and Other Tumors

Background

Synchronous tumors can be independent primary tumors or a primary-metastatic (clonal) pair, which may have clinical implications. Mutational profiling of tumor DNA is increasingly common in the clinic. We investigated whether mutational profiling can distinguish independent from clonal tumors in breast and other cancers, using a carefully defined test based on the Clonal Likelihood Score (CLS = 100 x # shared high confidence (HC) mutations/ # total HC mutations).

Methods

Statistical properties of a formal test using the CLS were investigated. A high CLS is evidence in favor of clonality; the test is implemented as a one-sided binomial test of proportions. Test parameters were empirically determined using 16,422 independent breast tumor pairs and 15 primary-metastatic tumor pairs from 10 cancer types using The Cancer Genome Atlas.

Results

We validated performance of the test with its established parameters, using five published data sets comprising 15,758 known independent tumor pairs (maximum CLS = 4.1%, minimum p-value = 0.48) and 283 known tumor clonal pairs (minimum CLS 13%, maximum p-value <0.01), across renal cell, testicular, and colorectal cancer. The CLS test correctly classified all validation samples but one, which it appears may have been incorrectly classified in the published data. As proof-of-concept we then applied the CLS test to two new cases of invasive synchronous bilateral breast cancer at our institution, each with one hormone receptor positive (ER+/PR+/HER2-) lobular and one triple negative ductal carcinoma. High confidence mutations were identified by exome sequencing and results were validated using deep targeted sequencing. The first tumor pair had CLS of 81% (p-value < 10–15), supporting clonality. In the second pair, no common mutations of 184 variants were validated (p-value >0.99), supporting independence. A plausible molecular mechanism for the shift from hormone receptor positive to triple negative was identified in the clonal pair.

Conclusion

We have developed the statistical properties of a carefully defined Clonal Likelihood Score test from mutational profiling of tumor DNA. Under identified conditions, the test appears to reliably distinguish between synchronous tumors of clonal and of independent origin in several cancer types. This approach may have scientific and clinical utility.

Integrative genomics with mediation analysis in a survival context

DNA copy number aberrations (DCNA) and subsequent altered gene expression profiles may have a major impact on tumor initiation, on development, and eventually on recurrence and cancer-specific mortality. However, most methods employed in integrative genomic analysis of the two biological levels, DNA and RNA, do not consider survival time. In the present note, we propose the adoption of a survival analysis-based framework for the integrative analysis of DCNA and mRNA levels to reveal their implication on patient clinical outcome with the prerequisite that the effect of DCNA on survival is mediated by mRNA levels. The specific aim of the paper is to offer a feasible framework to test the DCNA-mRNA-survival pathway. We provide statistical inference algorithms for mediation based on asymptotic results. Furthermore, we illustrate the applicability of the method in an integrative genomic analysis setting by using a breast cancer data set consisting of 141 invasive breast tumors. In addition, we provide implementation in R.

Genomic aberrations in normal tissue adjacent to HER2-amplified breast cancers: field cancerization or contaminating tumor cells?

Field cancerization effects as well as isolated tumor cell foci extending well beyond the invasive tumor margin have been described previously to account for local recurrence rates following breast conserving surgery despite adequate surgical margins and breast radiotherapy. To look for evidence of possible tumor cell contamination or field cancerization by genetic effects, a pilot study (Study 1: 12 sample pairs) followed by a verification study (Study 2: 20 sample pairs) were performed on DNA extracted from HER2-positive breast tumors and matching normal adjacent mammary tissue samples excised 1-3 cm beyond the invasive tumor margin. High-resolution molecular inversion probe (MIP) arrays were used to compare genomic copy number variations, including increased HER2 gene copies, between the paired samples; as well, a detailed histologic and immunohistochemical (IHC) re-evaluation of all Study 2 samples was performed blinded to the genomic results to characterize the adjacent normal tissue composition bracketing the DNA-extracted samples. Overall, 14/32 (44 %) sample pairs from both studies produced genome-wide evidence of genetic aberrations including HER2 copy number gains within the adjacent normal tissue samples. The observed single-parental origin of monoallelic HER2 amplicon haplotypes shared by informative tumor-normal pairs, as well as commonly gained loci elsewhere on 17q, suggested the presence of contaminating tumor cells in the genomically aberrant normal samples. Histologic and IHC analyses identified occult 25-200 μm tumor cell clusters overexpressing HER2 scattered in more than half, but not all, of the genomically aberrant normal samples re-evaluated, but in none of the genomically normal samples. These genomic and microscopic findings support the conclusion that tumor cell contamination rather than genetic field cancerization represents the likeliest cause of local clinical recurrence rates following breast conserving surgery, and mandate caution in assuming the genomic normalcy of histologically benign appearing peritumor breast tissue.

Clonality: an R package for testing clonal relatedness of two tumors from the same patient based on their genomic profiles

SUMMARY

If a cancer patient develops multiple tumors, it is sometimes impossible to determine whether these tumors are independent or clonal based solely on pathological characteristics. Investigators have studied how to improve this diagnostic challenge by comparing the presence of loss of heterozygosity (LOH) at selected genetic locations of tumor samples, or by comparing genomewide copy number array profiles. We have previously developed statistical methodology to compare such genomic profiles for an evidence of clonality. We assembled the software for these tests in a new R package called 'Clonality'. For LOH profiles, the package contains significance tests. The analysis of copy number profiles includes a likelihood ratio statistic and reference distribution, as well as an option to produce various plots that summarize the results.

AVAILABILITY

Bioconductor (http://bioconductor.org/packages/release/bioc/html/Clonality.html) and http://www.mskcc.org/mskcc/html/13287.cfm.