Rare subclonal sequencing of breast cancers indicates putative metastatic driver mutations are predominately acquired after dissemination.

BACKGROUND: Evolutionary models of breast cancer progression differ on the extent to which metastatic potential is pre-encoded within primary tumors. Although metastatic recurrences often harbor putative driver mutations that are not detected in their antecedent primary tumor using standard sequencing technologies, whether these mutations were acquired before or after dissemination remains unclear. METHODS: To ascertain whether putative metastatic driver mutations initially deemed specific to the metastasis by whole exome sequencing were, in actuality, present within rare ancestral subclones of the primary tumors from which they arose, we employed error-controlled ultra-deep sequencing (UDS-UMI) coupled with FFPE artifact mitigation by uracil-DNA glycosylase (UDG) to assess the presence of 132 "metastasis-specific" mutations within antecedent primary tumors from 21 patients. Maximum mutation detection sensitivity was ~1% of primary tumor cells. A conceptual framework was developed to estimate relative likelihoods of alternative models of mutation acquisition. RESULTS: The ancestral primary tumor subclone responsible for seeding the metastasis was identified in 29% of patients, implicating several putative drivers in metastatic seeding including LRP5 A65V and PEAK1 K140Q. Despite this, 93% of metastasis-specific mutations in putative metastatic driver genes remained undetected within primary tumors, as did 96% of metastasis-specific mutations in known breast cancer drivers, including ERRB2 V777L, ESR1 D538G, and AKT1 D323H. Strikingly, even in those cases in which the rare ancestral subclone was identified, 87% of metastasis-specific putative driver mutations remained undetected. Modeling indicated that the sequential acquisition of multiple metastasis-specific driver or passenger mutations within the same rare subclonal lineage of the primary tumor was highly improbable. CONCLUSIONS: Our results strongly suggest that metastatic driver mutations are sequentially acquired and selected within the same clonal lineage both before, but more commonly after, dissemination from the primary tumor, and that these mutations are biologically consequential. Despite inherent limitations in sampling archival primary tumors, our findings indicate that tumor cells in most patients continue to undergo clinically relevant genomic evolution after their dissemination from the primary tumor. This provides further evidence that metastatic recurrence is a multi-step, mutation-driven process that extends beyond primary tumor dissemination and underscores the importance of longitudinal tumor assessment to help guide clinical decisions.

Prognostic correlations with the microbiome of breast cancer subtypes.

Alterations to the natural microbiome are linked to different diseases, and the presence or absence of specific microbes is directly related to disease outcomes. We performed a comprehensive analysis with unique cohorts of the four subtypes of breast cancer (BC) characterized by their microbial signatures, using a pan-pathogen microarray strategy. The signature (includes viruses, bacteria, fungi, and parasites) of each tumor subtype was correlated with clinical data to identify microbes with prognostic potential. The subtypes of BC had specific viromes and microbiomes, with ER+ and TN tumors showing the most and least diverse microbiome, respectively. The specific microbial signatures allowed discrimination between different BC subtypes. Furthermore, we demonstrated correlations between the presence and absence of specific microbes in BC subtypes with the clinical outcomes. This study provides a comprehensive map of the oncobiome of BC subtypes, with insights into disease prognosis that can be critical for precision therapeutic intervention strategies.

Computationally Derived Cribriform Area Index from Prostate Cancer Hematoxylin and Eosin Images Is Associated with Biochemical Recurrence Following Radical Prostatectomy and Is Most Prognostic in Gleason Grade Group 2.

BACKGROUND: The presence of invasive cribriform adenocarcinoma (ICC), an expanse of cells containing punched-out lumina uninterrupted by stroma, in radical prostatectomy (RP) specimens has been associated with biochemical recurrence (BCR). However, ICC identification has only moderate inter-reviewer agreement. OBJECTIVE: To investigate quantitative machine-based assessment of the extent and prognostic utility of ICC, especially within individual Gleason grade groups. DESIGN, SETTING, AND PARTICIPANTS: A machine learning approach was developed for ICC segmentation using 70 RP patients and validated in a cohort of 749 patients from four sites whose median year of surgery was 2007 and with median follow-up of 28 mo. ICC was segmented on one representative hematoxylin and eosin RP slide per patient and the fraction of tumor area composed of ICC, the cribriform area index (CAI), was measured. OUTCOME MEASUREMENTS AND STATISTICAL ANALYSIS: The association between CAI and BCR was measured in terms of the concordance index (c index) and hazard ratio (HR). RESULTS AND LIMITATIONS: CAI was correlated with BCR (c index 0.62) in the validation set of 411 patients with ICC morphology, especially those with Gleason grade group 2 cancer (n = 192; c index 0.66), and was less prognostic when patients without ICC were included (c index 0.54). A doubling of CAI in the group with ICC morphology was prognostic after controlling for Gleason grade, surgical margin positivity, preoperative prostate-specific antigen level, pathological T stage, and age (HR 1.19, 95% confidence interval 1.03-1.38; p = 0.018). CONCLUSIONS: Automated image analysis and machine learning could provide an objective, quantitative, reproducible, and high-throughput method of quantifying ICC area. The performance of CAI for grade group 2 cancer suggests that for patients with little Gleason 4 pattern, the ICC fraction has a strong prognostic role. PATIENT SUMMARY: Machine-based measurement of a specific cell pattern (cribriform; sieve-like, with lots of spaces) in images of prostate specimens could improve risk stratification for patients with prostate cancer. In the future, this could help in expanding the criteria for active surveillance.

Stable and discriminating features are predictive of cancer presence and Gleason grade in radical prostatectomy specimens: a multi-site study.

Site variation in fixation, staining, and scanning can confound automated tissue based image classifiers for disease characterization. In this study we incorporated stability into four feature selection methods for identifying the most robust and discriminating features for two prostate histopathology classification tasks. We evaluated 242 morphology features from N = 212 prostatectomy specimens from four sites for automated cancer detection and grading. We quantified instability as the rate of significant cross-site feature differences. We mapped feature stability and discriminability using 188 non-cancerous and 210 cancerous regions via 3-fold cross validation, then held one site out, creating independent training and testing sets. In training, one feature set was selected only for discriminability, another for discriminability and stability. We trained a classifier with each feature set, testing on the hold out site. Experiments were repeated with 117 Gleason grade 3 and 112 grade 4 regions. Stability was calculated across non-cancerous regions. Gland shape features yielded the best stability and area under the receiver operating curve (AUC) trade-off while co-occurrence texture features were generally unstable. Our stability-informed method produced a cancer detection AUC of 0.98 ± 0.05 and increased average Gleason grading AUC by 4.38%. Color normalization of the images tended to exacerbate feature instability.

Accurate and reproducible invasive breast cancer detection in whole-slide images: A Deep Learning approach for quantifying tumor extent.

With the increasing ability to routinely and rapidly digitize whole slide images with slide scanners, there has been interest in developing computerized image analysis algorithms for automated detection of disease extent from digital pathology images. The manual identification of presence and extent of breast cancer by a pathologist is critical for patient management for tumor staging and assessing treatment response. However, this process is tedious and subject to inter- and intra-reader variability. For computerized methods to be useful as decision support tools, they need to be resilient to data acquired from different sources, different staining and cutting protocols and different scanners. The objective of this study was to evaluate the accuracy and robustness of a deep learning-based method to automatically identify the extent of invasive tumor on digitized images. Here, we present a new method that employs a convolutional neural network for detecting presence of invasive tumor on whole slide images. Our approach involves training the classifier on nearly 400 exemplars from multiple different sites, and scanners, and then independently validating on almost 200 cases from The Cancer Genome Atlas. Our approach yielded a Dice coefficient of 75.86%, a positive predictive value of 71.62% and a negative predictive value of 96.77% in terms of pixel-by-pixel evaluation compared to manually annotated regions of invasive ductal carcinoma.

Evaluating stability of histomorphometric features across scanner and staining variations: prostate cancer diagnosis from whole slide images.

Quantitative histomorphometry (QH) is the process of computerized feature extraction from digitized tissue slide images to predict disease presence, behavior, and outcome. Feature stability between sites may be compromised by laboratory-specific variables including dye batch, slice thickness, and the whole slide scanner used. We present two new measures, preparation-induced instability score and latent instability score, to quantify feature instability across and within datasets. In a use case involving prostate cancer, we examined QH features which may detect cancer on whole slide images. Using our method, we found that five feature families (graph, shape, co-occurring gland tensor, sub-graph, and texture) were different between datasets in 19.7% to 48.6% of comparisons while the values expected without site variation were 4.2% to 4.6%. Color normalizing all images to a template did not reduce instability. Scanning the same 34 slides on three scanners demonstrated that Haralick features were most substantively affected by scanner variation, being unstable in 62% of comparisons. We found that unstable feature families performed significantly worse in inter- than intrasite classification. Our results appear to suggest QH features should be evaluated across sites to assess robustness, and class discriminability alone should not represent the benchmark for digital pathology feature selection.

Self-renewal of CD133(hi) cells by IL6/Notch3 signalling regulates endocrine resistance in metastatic breast cancer.

The mechanisms of metastatic progression from hormonal therapy (HT) are largely unknown in luminal breast cancer. Here we demonstrate the enrichment of CD133(hi)/ER(lo) cancer cells in clinical specimens following neoadjuvant endocrine therapy and in HT refractory metastatic disease. We develop experimental models of metastatic luminal breast cancer and demonstrate that HT can promote the generation of HT-resistant, self-renewing CD133(hi)/ER(lo)/IL6(hi) cancer stem cells (CSCs). HT initially abrogates oxidative phosphorylation (OXPHOS) generating self-renewal-deficient cancer cells, CD133(hi)/ER(lo)/OXPHOS(lo). These cells exit metabolic dormancy via an IL6-driven feed-forward ER(lo)-IL6(hi)-Notch(hi) loop, activating OXPHOS, in the absence of ER activity. The inhibition of IL6R/IL6-Notch pathways switches the self-renewal of CD133(hi) CSCs, from an IL6/Notch-dependent one to an ER-dependent one, through the re-expression of ER. Thus, HT induces an OXPHOS metabolic editing of luminal breast cancers, paradoxically establishing HT-driven self-renewal of dormant CD133(hi)/ER(lo) cells mediating metastatic progression, which is sensitive to dual targeted therapy.

Computer-extracted Features Can Distinguish Noncancerous Confounding Disease from Prostatic Adenocarcinoma at Multiparametric MR Imaging.

PURPOSE: To determine the best features to discriminate prostate cancer from benign disease and its relationship to benign disease class and cancer grade. MATERIALS AND METHODS: The institutional review board approved this study and waived the need for informed consent. A retrospective cohort of 70 patients (age range, 48-70 years; median, 62 years), all of whom were scheduled to undergo radical prostatectomy and underwent preoperative 3-T multiparametric magnetic resonance (MR) imaging, including T2-weighted, diffusion-weighted, and dynamic contrast material-enhanced imaging, were included. The digitized prostatectomy slides were annotated for cancer and noncancerous disease and coregistered to MR imaging with an interactive deformable coregistration scheme. Computer-identified features for each of the noncancerous disease categories (eg, benign prostatic hyperplasia [BPH], prostatic intraepithelial neoplasia [PIN], inflammation, and atrophy) and prostate cancer were extracted. Feature selection was performed to identify the features with the highest discriminatory power. The performance of these five features was evaluated by using the area under the receiver operating characteristic curve (AUC). RESULTS: High-b-value diffusion-weighted images were more discriminative in distinguishing BPH from prostate cancer than apparent diffusion coefficient, which was most suitable for distinguishing PIN from prostate cancer. The focal appearance of lesions on dynamic contrast-enhanced images may help discriminate atrophy and inflammation from cancer. Which imaging features are discriminative for different benign lesions is influenced by cancer grade. The apparent diffusion coefficient appeared to be the most discriminative feature in identifying high-grade cancer. Classification results showed increased performance by taking into account specific benign types (AUC = 0.70) compared with grouping all noncancerous findings together (AUC = 0.62). CONCLUSION: The best features with which to discriminate prostate cancer from noncancerous benign disease depend on the type of benign disease and cancer grade. Use of the best features may result in better diagnostic performance.

Identifying in vivo DCE MRI markers associated with microvessel architecture and gleason grades of prostate cancer.

BACKGROUND: To identify computer extracted in vivo dynamic contrast enhanced (DCE) MRI markers associated with quantitative histomorphometric (QH) characteristics of microvessels and Gleason scores (GS) in prostate cancer. METHODS: This study considered retrospective data from 23 biopsy confirmed prostate cancer patients who underwent 3 Tesla multiparametric MRI before radical prostatectomy (RP). Representative slices from RP specimens were stained with vascular marker CD31. Tumor extent was mapped from RP sections onto DCE MRI using nonlinear registration methods. Seventy-seven microvessel QH features and 18 DCE MRI kinetic features were extracted and evaluated for their ability to distinguish low from intermediate and high GS. The effect of temporal sampling on kinetic features was assessed and correlations between those robust to temporal resolution and microvessel features discriminative of GS were examined. RESULTS: A total of 12 microvessel architectural features were discriminative of low and intermediate/high grade tumors with area under the receiver operating characteristic curve (AUC) > 0.7. These features were most highly correlated with mean washout gradient (WG) (max rho = -0.62). Independent analysis revealed WG to be moderately robust to temporal resolution (intraclass correlation coefficient [ICC] = 0.63) and WG variance, which was poorly correlated with microvessel features, to be predictive of low grade tumors (AUC = 0.77). Enhancement ratio was the most robust (ICC = 0.96) and discriminative (AUC = 0.78) kinetic feature but was moderately correlated with microvessel features (max rho = -0.52). CONCLUSION: Computer extracted features of prostate DCE MRI appear to be correlated with microvessel architecture and may be discriminative of low versus intermediate and high GS.

Risk factors for gyrA and parC mutations in Pseudomonas aeruginosa.

OBJECTIVE: The major mechanism of fluoroquinolone (FQ) resistance in Pseudomonas aeruginosa (PSA) is modification of target proteins in DNA gyrase and topoisomerase IV, most commonly the gyrA and parC subunits. The objective of this study was to determine risk factors for PSA with and without gyrA or parC mutations. DESIGN: Case-case-control study SETTING: Two adult academic acute-care hospitals PATIENTS: Case 1 study participants had a PSA isolate on hospital day 3 or later with any gyrA or parC mutation; case 2 study participants had a PSA isolate on hospital day 3 or later without these mutations. Controls were a random sample of all inpatients with a stay of 3 days or more. METHODS: Each case group was compared to the control group in separate multivariate models on the basis of demographics and inpatient antibiotic exposure, and risk factors were qualitatively compared. RESULTS: Of 298 PSA isolates, 172 (57.7%) had at least 1 mutation. Exposure to vancomycin and other agents with extended Gram-positive activity was a risk factor for both cases (case 1 odds ratio [OR], 1.09; 95% confidence interval [CI], 1.04-1.13; OR, 1.14; 95% CI, 1.03-1.26; case 2 OR, 1.09; 95% CI, 1.03-1.14; OR, 1.13; 95% CI, 1.01-1.25, respectively). CONCLUSIONS: Exposure to agents with extended Gram-positive activity is a risk factor for isolation of PSA overall but not for gyrA/parC mutations. FQ exposure is not associated with isolation of PSA with mutations.

Supervised multi-view canonical correlation analysis (sMVCCA): integrating histologic and proteomic features for predicting recurrent prostate cancer.

In this work, we present a new methodology to facilitate prediction of recurrent prostate cancer (CaP) following radical prostatectomy (RP) via the integration of quantitative image features and protein expression in the excised prostate. Creating a fused predictor from high-dimensional data streams is challenging because the classifier must 1) account for the "curse of dimensionality" problem, which hinders classifier performance when the number of features exceeds the number of patient studies and 2) balance potential mismatches in the number of features across different channels to avoid classifier bias towards channels with more features. Our new data integration methodology, supervised Multi-view Canonical Correlation Analysis (sMVCCA), aims to integrate infinite views of highdimensional data to provide more amenable data representations for disease classification. Additionally, we demonstrate sMVCCA using Spearman's rank correlation which, unlike Pearson's correlation, can account for nonlinear correlations and outliers. Forty CaP patients with pathological Gleason scores 6-8 were considered for this study. 21 of these men revealed biochemical recurrence (BCR) following RP, while 19 did not. For each patient, 189 quantitative histomorphometric attributes and 650 protein expression levels were extracted from the primary tumor nodule. The fused histomorphometric/proteomic representation via sMVCCA combined with a random forest classifier predicted BCR with a mean AUC of 0.74 and a maximum AUC of 0.9286. We found sMVCCA to perform statistically significantly (p < 0.05) better than comparative state-of-the-art data fusion strategies for predicting BCR. Furthermore, Kaplan-Meier analysis demonstrated improved BCR-free survival prediction for the sMVCCA-fused classifier as compared to histology or proteomic features alone.

Co-occurring gland angularity in localized subgraphs: predicting biochemical recurrence in intermediate-risk prostate cancer patients.

Quantitative histomorphometry (QH) refers to the application of advanced computational image analysis to reproducibly describe disease appearance on digitized histopathology images. QH thus could serve as an important complementary tool for pathologists in interrogating and interpreting cancer morphology and malignancy. In the US, annually, over 60,000 prostate cancer patients undergo radical prostatectomy treatment. Around 10,000 of these men experience biochemical recurrence within 5 years of surgery, a marker for local or distant disease recurrence. The ability to predict the risk of biochemical recurrence soon after surgery could allow for adjuvant therapies to be prescribed as necessary to improve long term treatment outcomes. The underlying hypothesis with our approach, co-occurring gland angularity (CGA), is that in benign or less aggressive prostate cancer, gland orientations within local neighborhoods are similar to each other but are more chaotically arranged in aggressive disease. By modeling the extent of the disorder, we can differentiate surgically removed prostate tissue sections from (a) benign and malignant regions and (b) more and less aggressive prostate cancer. For a cohort of 40 intermediate-risk (mostly Gleason sum 7) surgically cured prostate cancer patients where half suffered biochemical recurrence, the CGA features were able to predict biochemical recurrence with 73% accuracy. Additionally, for 80 regions of interest chosen from the 40 studies, corresponding to both normal and cancerous cases, the CGA features yielded a 99% accuracy. CGAs were shown to be statistically signicantly ([Formula: see text]) better at predicting BCR compared to state-of-the-art QH methods and postoperative prostate cancer nomograms.

A unique pancreatic tumor with exclusive hepatocytic differentiation.

Only 7 cases of pancreatic tumor with hepatocytic differentiation have been reported in the literature, including 6 cases of hepatoid carcinoma and one case of hepatoid adenoma. Diagnosis of hepatoid carcinoma depends on recognition of characteristic histological features, supported by other evidence linked to hepatic lineage including alpha-fetoprotein production, positive immunoreactivity to liver synthesized proteins, and in situ hybridization detection of albumin mRNA. In addition, a synchronous focus of carcinoma arising in pancreatic ducts, islet cells, or acinar cells is essential. We report a unique case of pancreatic tumor with exclusive hepatocytic differentiation. In this tumor, we were unable to find a synchronous focus of carcinoma arising in pancreatic ducts, islet cells, or acinar cells, ruling out the possibility of its being hepatoid carcinoma. Long term follow-up can help to determine whether this tumor is benign or malignant. The patterns of reticulin staining and immunohistochemical staining are suggestive of malignancy, but mitotic activity is low and nuclear pleomorphism is minimal.