Comparative genomics of primary prostate cancer and paired metastases: insights from 12 molecular case studies.

Primary prostate cancer (PCa) can show marked molecular heterogeneity. However, systematic analyses comparing primary PCa and matched metastases in individual patients are lacking. We aimed to address the molecular aspects of metastatic progression while accounting for the heterogeneity of primary PCa. In this pilot study, we collected 12 radical prostatectomy (RP) specimens from men who subsequently developed metastatic castration-resistant prostate cancer (mCRPC). We used histomorphology (Gleason grade, focus size, stage) and immunohistochemistry (IHC) (ERG and p53) to identify independent tumors and/or distinct subclones of primary PCa. We then compared molecular profiles of these primary PCa areas to matched metastatic samples using whole-exome sequencing (WES) and amplicon-based DNA and RNA sequencing. Based on combined pathology and molecular analysis, seven (58%) RP specimens harbored monoclonal and topographically continuous disease, albeit with some degree of intratumor heterogeneity; four (33%) specimens showed true multifocal disease; and one displayed monoclonal disease with discontinuous topography. Early (truncal) events in primary PCa included SPOP p.F133V (one patient), BRAF p.K601E (one patient), and TMPRSS2:ETS rearrangements (eight patients). Activating AR alterations were seen in nine (75%) mCRPC patients, but not in matched primary PCa. Hotspot TP53 mutations, found in metastases from three patients, were readily present in matched primary disease. Alterations in genes encoding epigenetic modifiers were observed in several patients (either shared between primary foci and metastases or in metastatic samples only). WES-based phylogenetic reconstruction and/or clonality scores were consistent with the index focus designated by pathology review in six out of nine (67%) cases. The three instances of discordance pertained to monoclonal, topographically continuous tumors, which would have been considered as unique disease in routine practice. Overall, our results emphasize pathologic and molecular heterogeneity of primary PCa, and suggest that comprehensive IHC-assisted pathology review and genomic analysis are highly concordant in nominating the 'index' primary PCa area. © 2022 The Authors. The Journal of Pathology published by John Wiley & Sons Ltd on behalf of The Pathological Society of Great Britain and Ireland.

Discordant and heterogeneous clinically relevant genomic alterations in circulating tumor cells vs plasma DNA from men with metastatic castration resistant prostate cancer.

Circulating tumor cell (CTC) and cell-free (cf) DNA-based genomic alterations are increasingly being used for clinical decision-making in oncology. However, the concordance and discordance between paired CTC and cfDNA genomic profiles remain largely unknown. We performed comparative genomic hybridization (CGH) on CTCs and cfDNA, and low-pass whole genome sequencing (lpWGS) on cfDNA to characterize genomic alterations (CNA) and tumor content in two independent prospective studies of 93 men with mCRPC treated with enzalutamide/abiraterone, or radium-223. Comprehensive analysis of 69 patient CTCs and 72 cfDNA samples from 93 men with mCRPC, including 64 paired samples, identified common concordant gains in FOXA1, AR, and MYC, and losses in BRCA1, PTEN, and RB1 between CTCs and cfDNA. Concordant PTEN loss and discordant BRCA2 gain were associated with significantly worse outcomes in Epic AR-V7 negative men with mCRPC treated with abiraterone/enzalutamide. We identified and externally validated CTC-specific genomic alternations that were discordant in paired cfDNA, even in samples with high tumor content. These CTC/cfDNA-discordant regions included key genomic regulators of lineage plasticity, osteomimicry, and cellular differentiation, including MYCN gain in CTCs (31%) that was rarely detected in cfDNA. CTC MYCN gain was associated with poor clinical outcomes in AR-V7 negative men and small cell transformation. In conclusion, we demonstrated concordance of multiple genomic alterations across CTC and cfDNA platforms; however, some genomic alterations displayed substantial discordance between CTC DNA and cfDNA despite the use of identical copy number analysis methods, suggesting tumor heterogeneity and divergent evolution associated with poor clinical outcomes.

Aldosterone-stimulating somatic gene mutations are common in normal adrenal glands.

Primary aldosteronism (PA) represents the most common cause of secondary hypertension, but little is known regarding its adrenal cellular origins. Recently, aldosterone-producing cell clusters (APCCs) with high expression of aldosterone synthase (CYP11B2) were found in both normal and PA adrenal tissue. PA-causing aldosterone-producing adenomas (APAs) harbor mutations in genes encoding ion channels/pumps that alter intracellular calcium homeostasis and cause renin-independent aldosterone production through increased CYP11B2 expression. Herein, we hypothesized that APCCs have APA-related aldosterone-stimulating somatic gene mutations. APCCs were studied in 42 normal adrenals from kidney donors. To clarify APCC molecular characteristics, we used microarrays to compare the APCC transcriptome with conventional adrenocortical zones [zona glomerulosa (ZG), zona fasciculata, and zona reticularis]. The APCC transcriptome was most similar to ZG but with an enhanced capacity to produce aldosterone. To determine if APCCs harbored APA-related mutations, we performed targeted next generation sequencing of DNA from 23 APCCs and adjacent normal adrenal tissue isolated from both formalin-fixed, paraffin-embedded, and frozen tissues. Known aldosterone driver mutations were identified in 8 of 23 (35%) APCCs, including mutations in calcium channel, voltage-dependent, L-type, α1D-subunit (CACNA1D; 6 of 23 APCCs) and ATPase, Na(+)/(K+) transporting, α1-polypeptide (ATP1A1; 2 of 23 APCCs), which were not observed in the adjacent normal adrenal tissue. Overall, we show three major findings: (i) APCCs are common in normal adrenals, (ii) APCCs harbor somatic mutations known to cause excess aldosterone production, and (iii) the mutation spectrum of aldosterone-driving mutations is different in APCCs from that seen in APA. These results provide molecular support for APCC as a precursor of PA.

Characterization of ADME gene variation in 21 populations by exome sequencing.

OBJECTIVE: Proteins involving absorption, distribution, metabolism, and excretion (ADME) play a critical role in drug pharmacokinetics. The type and frequency of genetic variation in the ADME genes differ among populations. The aim of this study was to systematically investigate common and rare ADME coding variation in diverse ethnic populations by exome sequencing. MATERIALS AND METHODS: Data derived from commercial exome capture arrays and next-generation sequencing were used to characterize coding variation in 298 ADME genes in 251 Northeast Asians and 1181 individuals from the 1000 Genomes Project. RESULTS: Approximately 75% of the ADME coding sequence was captured at high quality across the joint samples harboring more than 8000 variants, with 49% of individuals carrying at least one 'knockout' allele. ADME genes carried 50% more nonsynonymous variation than non-ADME genes (P=8.2×10) and showed significantly greater levels of population differentiation (P=7.6×10). Out of the 2135 variants identified that were predicted to be deleterious, 633 were not on commercially available ADME or general-purpose genotyping arrays. Forty deleterious variants within important ADME genes, with frequencies of at least 2% in at least one population, were identified as candidates for future pharmacogenetic studies. CONCLUSION: Exome sequencing was effective in accurately genotyping most ADME variants important for pharmacogenetic research, in addition to identifying rare or potentially de novo coding variants that may be clinically meaningful. Furthermore, as a class, ADME genes are more variable and less sensitive to purifying selection than non-ADME genes.

Intra-tumor genetic heterogeneity in rectal cancer.

Colorectal cancer arises in part from the cumulative effects of multiple gene lesions. Recent studies in selected cancer types have revealed significant intra-tumor genetic heterogeneity and highlighted its potential role in disease progression and resistance to therapy. We hypothesized the existence of significant intra-tumor genetic heterogeneity in rectal cancers involving variations in localized somatic mutations and copy number abnormalities. Two or three spatially disparate regions from each of six rectal tumors were dissected and subjected to the next-generation whole-exome DNA sequencing, Oncoscan SNP arrays, and targeted confirmatory sequencing and analysis. The resulting data were integrated to define subclones using SciClone. Mutant-allele tumor heterogeneity (MATH) scores, mutant allele frequency correlation, and mutation percent concordance were calculated, and copy number analysis including measurement of correlation between samples was performed. Somatic mutations profiles in individual cancers were similar to prior studies, with some variants found in previously reported significantly mutated genes and many patient-specific mutations in each tumor. Significant intra-tumor heterogeneity was identified in the spatially disparate regions of individual cancers. All tumors had some heterogeneity but the degree of heterogeneity was quite variable in the samples studied. We found that 67-97% of exonic somatic mutations were shared among all regions of an individual's tumor. The SciClone computational method identified 2-8 shared and unshared subclones in the spatially disparate areas in each tumor. MATH scores ranged from 7 to 41. Allele frequency correlation scores ranged from R(2)=0.69-0.96. Measurements of correlation between samples for copy number changes varied from R(2)=0.74-0.93. All tumors had some heterogeneity, but the degree was highly variable in the samples studied. The occurrence of significant intra-tumor heterogeneity may allow selected tumors to have a genetic reservoir to draw from in their evolutionary response to therapy and other challenges.

No large-effect low-frequency coding variation found for myocardial infarction.

Genome-wide association studies have identified variants, primarily common, that are associated with coronary artery disease or myocardial infarction (MI), but have not tested the majority of the low frequency and rare variation in the genome. We explored the hypothesis that previously untested low frequency (1-5% minor allele frequency) and rare (<1% minor allele frequency) coding variants are associated with MI. We genotyped 2906 MI cases and 6738 non-MI controls from Norway using the Illumina HumanExome Beadchip, allowing for direct genotyping of 85 972 polymorphic coding variants as well as 48 known GWAS SNPs. We followed-up 34 coding variants in an additional 2350 MI cases and 2318 controls from Norway. We evaluated exome array coverage in a subset of these samples using whole exome sequencing (N = 151). The exome array provided successful genotyping for an estimated 72.5% of Norwegian loss-of-function or missense variants with frequency >1% and 66.2% of variants <1% frequency observed more than once. Despite 80% power in the two-stage study (N = 14 312) to detect association with low-frequency variants with high effect sizes [odds ratio (OR) >1.86 and >1.36 for 1 and 5% frequency, respectively], we did not identify any novel genes or single variants that reached significance. This suggests that low-frequency coding variants with large effect sizes (OR >2) may not exist for MI. Larger sample sizes may identify coding variants with more moderate effects.

Comprehensive genomic profiling of orbital and ocular adnexal lymphomas identifies frequent alterations in MYD88 and chromatin modifiers: new routes to targeted therapies.

Non-Hodgkin lymphoma of the orbit and ocular adnexa is the most common primary orbital malignancy. Treatments for low- (extra-nodal marginal zone and follicular lymphomas) and high-grade (diffuse large B-cell lymphoma) are associated with local and vision-threatening toxicities. High-grade lymphomas relapse frequently and exhibit poor survival rates. Despite advances in genomic profiling and precision medicine, orbital and ocular adnexal lymphomas remain poorly characterized molecularly. We performed targeted next-generation sequencing (NGS) profiling of 38 formalin-fixed, paraffin-embedded orbital and ocular adnexal lymphomas obtained from a single-center using a panel targeting near-term, clinically relevant genes. Potentially actionable mutations and copy number alterations were prioritized based on gain- and loss-of-function analyses, and catalogued, approved, and investigational therapies. Of 36 informative samples, including marginal zone lymphomas (n=20), follicular lymphomas (n=9), and diffuse large B-cell lymphomas (n=7), 53% harbored a prioritized alteration (median=1, range 0-5/sample). MYD88 was the most frequently altered gene in our cohort, with potentially clinically relevant hotspot gain-of-function mutations identified in 71% of diffuse large B-cell lymphomas and 25% of marginal zone lymphomas. Prioritized alterations in epigenetic modulators were common and included gain-of-function EZH2 and loss-of-function ARID1A mutations (14% of diffuse large B-cell lymphomas and 22% of follicular lymphomas contained alterations in each of these two genes). Single prioritized alterations were also identified in the histone methyltransferases KMT2B (follicular lymphoma) and KMT3B (diffuse large B-cell lymphoma). Loss-of-function mutations and copy number alterations in the tumor suppressors TP53 (diffuse large B-cell and follicular lymphoma), CDKN2A (diffuse large B-cell and marginal zone lymphoma), PTEN (diffuse large B-cell lymphoma), ATM (diffuse large B-cell lymphoma), and NF1 (diffuse large B-cell lymphoma), and gain-of-function mutations in the oncogenes HRAS (follicular lymphoma) and NRAS (diffuse large B-cell lymphoma) were also observed. Together, our study demonstrates that NGS can be used to profile routine formalin-fixed, paraffin-embedded orbital and ocular adnexal lymphomas for identification of somatic-driving alterations and nomination of potential therapeutic strategies.

Next generation sequencing of Cytokeratin 20-negative Merkel cell carcinoma reveals ultraviolet-signature mutations and recurrent TP53 and RB1 inactivation.

Merkel cell carcinoma is a rare but highly aggressive cutaneous neuroendocrine carcinoma. Cytokeratin 20 (CK20) is expressed in ~95% of Merkel cell carcinomas and is useful for distinction from morphologically similar entities including metastatic small-cell lung carcinoma. Lack of CK20 expression may make diagnosis of Merkel cell carcinoma more challenging, and has unknown biological significance. Approximately 80% of CK20-positive Merkel cell carcinomas are associated with the oncogenic Merkel cell polyomavirus. Merkel cell carcinomas lacking Merkel cell polyomavirus display distinct genetic changes from Merkel cell polyomavirus-positive Merkel cell carcinoma, including RB1 inactivating mutations. Unlike CK20-positive Merkel cell carcinoma, the majority of CK20-negative Merkel cell carcinomas are Merkel cell polyomavirus-negative, suggesting CK20-negative Merkel cell carcinomas predominantly arise through virus-independent pathway(s) and may harbor additional genetic differences from conventional Merkel cell carcinoma. Hence, we analyzed 15 CK20-negative Merkel cell carcinoma tumors (10 Merkel cell polyomavirus-negative, four Merkel cell polyomavirus-positive, and one undetermined) using the Ion Ampliseq Comprehensive Cancer Panel, which assesses copy number alterations and mutations in 409 cancer-relevant genes. Twelve tumors displayed prioritized high-level chromosomal gains or losses (average 1.9 per tumor). Non-synonymous high-confidence somatic mutations were detected in 14 tumors (average 11.9 per tumor). Assessing all somatic coding mutations, an ultraviolet-signature mutational profile was present, and more prevalent in Merkel cell polyomavirus-negative tumors. Recurrent deleterious tumor suppressor mutations affected TP53 (9/15, 60%), RB1 (3/15, 20%), and BAP1 (2/15, 13%). Oncogenic activating mutations included PIK3CA (3/15, 20%), AKT1 (1/15, 7%) and EZH2 (1/15, 7%). In conclusion, CK20-negative Merkel cell carcinoma display overlapping genetic changes with CK20-positive Merkel cell carcinoma, including RB1 mutations restricted to Merkel cell polyomavirus-negative tumors. However, some CK20-negative Merkel cell carcinomas harbor mutations not previously described in Merkel cell carcinoma. Hence, CK20-negative Merkel cell carcinomas harbor diverse oncogenic drivers which may represent therapeutic targets in individual tumors.

ctDNA transiting into urine is ultrashort and facilitates noninvasive liquid biopsy of HPV+ oropharyngeal cancer.

BACKGROUNDTransrenal cell-free tumor DNA (TR-ctDNA), which transits from the bloodstream into urine, has the potential to enable noninvasive cancer detection for a wide variety of nonurologic cancer types.MethodsUsing whole-genome sequencing, we discovered that urine TR-ctDNA fragments across multiple cancer types are predominantly ultrashort (<50 bp) and, therefore, likely to be missed by conventional ctDNA assays. We developed an ultrashort droplet digital PCR assay to detect TR-ctDNA originating from HPV-associated oropharyngeal squamous cell carcinoma (HPV+ OPSCC) and confirmed that assaying ultrashort DNA is critical for sensitive cancer detection from urine samples.ResultsTR-ctDNA was concordant with plasma ctDNA for cancer detection in patients with HPV+ OPSCC. As proof of concept for using urine TR-ctDNA for posttreatment surveillance, in a small longitudinal case series, TR-ctDNA showed promise for noninvasive detection of recurrence of HPV+ OPSCC.ConclusionOur data indicate that focusing on ultrashort fragments of TR-ctDNA will be important for realizing the full potential of urine-based cancer diagnostics. This has implications for urine-based detection of a wide variety of cancer types and for facilitating access to care through at-home specimen collections.FundingNIH grants R33 CA229023, R21 CA225493; NIH/National Cancer Institute grants U01 CA183848, R01 CA184153, and P30CA046592; American Cancer Society RSG-18-062-01-TBG; American Cancer Society Mission Boost grant MBGI-22-056-01-MBG; and the A. Alfred Taubman Medical Research Institute.

Development and validation of an integrative pan-solid tumor predictor of PD-1/PD-L1 blockade benefit.

BACKGROUND: Anti-PD-1 and PD-L1 (collectively PD-[L]1) therapies are approved for many advanced solid tumors. Biomarkers beyond PD-L1 immunohistochemistry, microsatellite instability, and tumor mutation burden (TMB) may improve benefit prediction. METHODS: Using treatment data and genomic and transcriptomic tumor tissue profiling from an observational trial (NCT03061305), we developed Immunotherapy Response Score (IRS), a pan-tumor predictive model of PD-(L)1 benefit. IRS real-world progression free survival (rwPFS) and overall survival (OS) prediction was validated in an independent cohort of trial patients. RESULTS: Here, by Cox modeling, we develop IRS-which combines TMB with CD274, PDCD1, ADAM12 and TOP2A quantitative expression-to predict pembrolizumab rwPFS (648 patients; 26 tumor types; IRS-High or -Low groups). In the 248 patient validation cohort (248 patients; 24 tumor types; non-pembrolizumab PD-[L]1 monotherapy treatment), median rwPFS and OS are significantly longer in IRS-High vs. IRS-Low patients (rwPFS adjusted hazard ratio [aHR] 0.52, p = 0.003; OS aHR 0.49, p = 0.005); TMB alone does not significantly predict PD-(L)1 rwPFS nor OS. In 146 patients treated with systemic therapy prior to pembrolizumab monotherapy, pembrolizumab rwPFS is only significantly longer than immediately preceding therapy rwPFS in IRS-High patients (interaction test p = 0.001). In propensity matched lung cancer patients treated with first-line pembrolizumab monotherapy or pembrolizumab+chemotherapy, monotherapy rwPFS is significantly shorter in IRS-Low patients, but is not significantly different in IRS-High patients. Across 24,463 molecularly-evaluable trial patients, 7.6% of patients outside of monotherapy PD-(L)1 approved tumor types are IRS-High/TMB-Low. CONCLUSIONS: The validated, predictive, pan-tumor IRS model can expand PD-(L)1 monotherapy benefit outside currently approved indications.

Therapies activating the immune system (checkpoint inhibitors) have revolutionized the treatment of patients with advanced cancer, however new molecular tests may better identify patients who could benefit. Using treatment data and clinical molecular test results, we report the development and validation of Immunotherapy Response Score (IRS) to predict checkpoint inhibitor benefit. Across patients with more than 20 advanced cancer types, IRS better predicted checkpoint inhibitor benefit than currently available tests. Data from >20,000 patients showed that IRS identifies ~8% of patients with advanced cancer who may dramatically benefit from checkpoint inhibitors but would not receive them today based on currently available tests. Our approach may help clinicians to decide which patients should receive checkpoint inhibitors to treat their disease.

Validation of Immunotherapy Response Score as Predictive of Pan-solid Tumor Anti-PD-1/PD-L1 Benefit.

Immunotherapy response score (IRS) integrates tumor mutation burden (TMB) and quantitative expression biomarkers to predict anti-PD-1/PD-L1 [PD-(L)1] monotherapy benefit. Here, we evaluated IRS in additional cohorts. Patients from an observational trial (NCT03061305) treated with anti-PD-(L)1 monotherapy were included and assigned to IRS-High (-H) versus -Low (-L) groups. Associations with real-world progression-free survival (rwPFS) and overall survival (OS) were determined by Cox proportional hazards (CPH) modeling. Those with available PD-L1 IHC treated with anti-PD-(L)1 with or without chemotherapy were separately assessed. Patients treated with PD-(L)1 and/or chemotherapy (five relevant tumor types) were assigned to three IRS groups [IRS-L divided into IRS-Ultra-Low (-UL) and Intermediate-Low (-IL), and similarly assessed]. In the 352 patient anti-PD-(L)1 monotherapy validation cohort (31 tumor types), IRS-H versus IRS-L patients had significantly longer rwPFS and OS. IRS significantly improved CPH associations with rwPFS and OS beyond microsatellite instability (MSI)/TMB alone. In a 189 patient (10 tumor types) PD-L1 IHC comparison cohort, IRS, but not PD-L1 IHC nor TMB, was significantly associated with anti-PD-L1 rwPFS. In a 1,103-patient cohort (from five relevant tumor types), rwPFS did not significantly differ in IRS-UL patients treated with chemotherapy versus chemotherapy plus anti-PD-(L)1, nor in IRS-H patients treated with anti-PD-(L)1 versus anti-PD-(L)1 + chemotherapy. IRS associations were consistent across subgroups, including both Europeans and non-Europeans. These results confirm the utility of IRS utility for predicting pan-solid tumor PD-(L)1 monotherapy benefit beyond available biomarkers and demonstrate utility for informing on anti-PD-(L)1 and/or chemotherapy treatment. Significance: This study confirms the utility of the integrative IRS biomarker for predicting anti-PD-L1/PD-1 benefit. IRS significantly improved upon currently available biomarkers, including PD-L1 IHC, TMB, and MSI status. Additional utility for informing on chemotherapy, anti-PD-L1/PD-1, and anti-PD-L1/PD-1 plus chemotherapy treatments decisions is shown.

Development of a Whole-urine, Multiplexed, Next-generation RNA-sequencing Assay for Early Detection of Aggressive Prostate Cancer.

BACKGROUND: Despite biomarker development advances, early detection of aggressive prostate cancer (PCa) remains challenging. We previously developed a clinical-grade urine test (Michigan Prostate Score [MiPS]) for individualized aggressive PCa risk prediction. MiPS combines serum prostate-specific antigen (PSA), the TMPRSS2:ERG (T2:ERG) gene fusion, and PCA3 lncRNA in whole urine after digital rectal examination (DRE). OBJECTIVE: To improve on MiPS with a novel next-generation sequencing (NGS) multibiomarker urine assay for early detection of aggressive PCa. DESIGN, SETTING, AND PARTICIPANTS: Preclinical development and validation of a post-DRE urine RNA NGS assay (Urine Prostate Seq [UPSeq]) assessing 84 PCa transcriptomic biomarkers, including T2:ERG, PCA3, additional PCa fusions/isoforms, mRNAs, lncRNAs, and expressed mutations. Our UPSeq model was trained on 73 patients and validated on a held-out set of 36 patients representing the spectrum of disease (benign to grade group [GG] 5 PCa). OUTCOME MEASUREMENTS AND STATISTICAL ANALYSIS: The area under the receiver operating characteristic curve (AUC) of UPSeq was compared with PSA, MiPS, and other existing models/biomarkers for predicting GG ≥3 PCa. RESULTS AND LIMITATIONS: UPSeq demonstrated high analytical accuracy and concordance with MiPS, and was able to detect expressed germline HOXB13 and somatic SPOP mutations. In an extreme design cohort (n = 109; benign/GG 1 vs GG ≥3 PCa, stratified to exclude GG 2 cancer in order to capture signal difference between extreme ends of disease), UPSeq showed differential expression for T2:ERG.T1E4 (1.2 vs 78.8 median normalized reads, p < 0.00001) and PCA3 (1024 vs 2521, p = 0.02), additional T2:ERG splice isoforms, and other candidate biomarkers. Using machine learning, we developed a 15-transcript model on the training set (n = 73) that outperformed serum PSA and sequencing-derived MiPS in predicting GG ≥3 PCa in the held-out validation set (n = 36; AUC 0.82 vs 0.69 and 0.69, respectively). CONCLUSIONS: These results support the potential utility of our novel urine-based RNA NGS assay to supplement PSA for improved early detection of aggressive PCa. PATIENT SUMMARY: We have developed a new urine-based test for the detection of aggressive prostate cancer, which promises improvement upon current biomarker tests.

Development and Validation of StrataNGS, a Multiplex PCR, Semiconductor Sequencing-Based Comprehensive Genomic Profiling Test.

Despite widespread use in targeted tumor testing, multiplex PCR/semiconductor (Ion Torrent) sequencing-based assessment of all comprehensive genomic profiling (CGP) variant classes has been limited. Herein, we describe the development and validation of StrataNGS, a 429-gene, multiplex PCR/semiconductor sequencing-based CGP laboratory-developed test performed on co-isolated DNA and RNA from formalin-fixed, paraffin-embedded tumor specimens with ≥2 mm2 tumor surface area. Validation was performed in accordance with MolDX CGP validation guidelines using 1986 clinical formalin-fixed, paraffin-embedded samples and an in-house developed optimized bioinformatics pipeline. Across CGP variant classes, accuracy ranged from 0.945 for tumor mutational burden (TMB) status to >0.999 for mutations and gene fusions, positive predictive value ranged from 0.915 for TMB status to 1.00 for gene fusions, and reproducibility ranged from 0.998 for copy number alterations to 1.00 for splice variants and insertions/deletions. StrataNGS TMB estimates were highly correlated to those from whole exome- or FoundationOne CDx-determined TMB (Pearson r = 0.998 and 0.960, respectively); TMB reproducibility was 0.996 (concordance correlation coefficient). Limit of detection for all variant classes was <20% tumor content. Together, we demonstrate that multiplex PCR/semiconductor sequencing-based tumor tissue CGP is feasible using optimized bioinformatic approaches described herein.

Real-World Performance of a Comprehensive Genomic Profiling Test Optimized for Small Tumor Samples.

PURPOSE: Tissue-based comprehensive genomic profiling (CGP) is increasingly used for treatment selection in patients with advanced cancer; however, tissue availability may limit widespread implementation. Here, we established real-world CGP tissue availability and assessed CGP performance on consecutively received samples. MATERIALS AND METHODS: We conducted a post hoc, nonprespecified analysis of 32,048 consecutive tumor tissue samples received for StrataNGS, a multiplex polymerase chain reaction (PCR)-based comprehensive genomic profiling (PCR-CGP) test, as part of an ongoing observational trial (NCT03061305). Sample characteristics and PCR-CGP performance were assessed across all tested samples, including exception samples not meeting minimum input quality control (QC) requirements (< 20% tumor content [TC], < 2 mm2 tumor surface area [TSA], DNA or RNA yield < 1 ng/µL, or specimen age > 5 years). Tests reporting ≥ 1 prioritized alteration or meeting TC and sequencing QC were considered successful. For prostate carcinoma and lung adenocarcinoma, tests reporting ≥ 1 actionable or informative alteration or meeting TC and sequencing QC were considered actionable. RESULTS: Among 31,165 (97.2%) samples where PCR-CGP was attempted, 10.7% had < 20% TC and 59.2% were small (< 25 mm2 tumor surface area). Of 31,101 samples evaluable for input requirements, 8,089 (26.0%) were exceptions not meeting requirements. However, 94.2% of the 31,101 tested samples were successfully reported, including 80.5% of exception samples. Positive predictive value of PCR-CGP for ERBB2 amplification in exceptions and/or sequencing QC-failure breast cancer samples was 96.7%. Importantly, 84.0% of tested prostate carcinomas and 87.9% of lung adenocarcinomas yielded results informing treatment selection. CONCLUSION: Most real-world tissue samples from patients with advanced cancer desiring CGP are limited, requiring optimized CGP approaches to produce meaningful results. An optimized PCR-CGP test, coupled with an inclusive exception testing policy, delivered reportable results for > 94% of samples, potentially expanding the proportion of CGP-testable patients and impact of biomarker-guided therapies.

CDKN1B Deletions are Associated with Metastasis in African American Men with Clinically Localized, Surgically Treated Prostate Cancer.

PURPOSE: The potential biological determinants of aggressive prostate cancer in African American (AA) men are unknown. Here we characterize prostate cancer genomic alterations in the largest cohort to date of AA men with clinical follow-up for metastasis, with the aim to elucidate the key molecular drivers associated with poor prognosis in this population. EXPERIMENTAL DESIGN: Targeted sequencing was retrospectively performed on 205 prostate tumors from AA men treated with radical prostatectomy (RP) to examine somatic genomic alterations and percent of the genome with copy-number alterations (PGA). Cox proportional hazards analyses assessed the association of genomic alterations with risk of metastasis. RESULTS: At RP, 71% (145/205) of patients had grade group ≥3 disease, and 49% (99/202) were non-organ confined. The median PGA was 3.7% (IQR = 0.9%-9.4%) and differed by pathologic grade (P < 0.001) and stage (P = 0.02). Median follow-up was 5 years. AA men with the highest quartile of PGA had increased risks of metastasis (multivariable: HR = 13.45; 95% CI, 2.55-70.86; P = 0.002). The most common somatic mutations were SPOP (11.2%), FOXA1 (8.3%), and TP53 (3.9%). The most common loci altered at the copy number level were CDKN1B (6.3%), CHD1 (4.4%), and PTEN (3.4%). TP53 mutations and deep deletions in CDKN1B were associated with increased risks of metastasis on multivariable analyses (TP53: HR = 9.5; 95% CI, 2.2-40.6; P = 0.002; CDKN1B: HR = 6.7; 95% CI, 1.3-35.2; P = 0.026). CONCLUSIONS: Overall, PGA, somatic TP53 mutations, and a novel finding of deep deletions in CDKN1B were associated with poor prognosis in AA men. These findings require confirmation in additional AA cohorts.

A conserved intratumoral regulatory T cell signature identifies 4-1BB as a pan-cancer target.

Despite advancements in targeting the immune checkpoints program cell death protein 1 (PD-1), programmed death ligand 1 (PD-L1), and cytotoxic T lymphocyte-associated protein 4 (CTLA-4) for cancer immunotherapy, a large number of patients and cancer types remain unresponsive. Current immunotherapies focus on modulating an antitumor immune response by directly or indirectly expanding antitumor CD8 T cells. A complementary strategy might involve inhibition of Tregs that otherwise suppress antitumor immune responses. Here, we sought to identify functional immune molecules preferentially expressed on tumor-infiltrating Tregs. Using genome-wide RNA-Seq analysis of purified Tregs sorted from multiple human cancer types, we identified a conserved Treg immune checkpoint signature. Using immunocompetent murine tumor models, we found that antibody-mediated depletion of 4-1BB-expressing cells (4-1BB is also known as TNFRSF9 or CD137) decreased tumor growth without negatively affecting CD8 T cell function. Furthermore, we found that the immune checkpoint 4-1BB had a high selectivity for human tumor Tregs and was associated with worse survival outcomes in patients with multiple tumor types. Thus, antibody-mediated depletion of 4-1BB-expressing Tregs represents a strategy with potential activity across cancer types.

Multiclonality and Marked Branched Evolution of Low-Grade Endometrioid Endometrial Carcinoma.

The molecular events driving low-grade endometrioid endometrial carcinoma (LGEC) development-like in many cancers-are incompletely understood. Hence, here we performed multiregion, comprehensive somatic molecular profiling of routinely processed formalin-fixed, paraffin-embedded (FFPE) material from 13 cases of LGEC totaling 64 minute, spatially defined cell populations ranging from presumed precursor lesions through invasive LGEC. Shared driving PTEN, PIK3R1, or PIK3CA mutations support clonal origin of the samples in each case, except for two cases with two clonally distinct neoplastic populations, consistent with unexpected multiclonality in LGEC development. Although substantial heterogeneity in driving somatic alterations was present across populations in nearly all cases, these alterations were usually clonal in a given population, supporting continued selection and clonal sweeping of driving alterations in populations with both precursor and LGEC histology. Importantly, CTNNB1 mutational status, which has been proposed as both prognostic and predictive in LGEC, was frequently heterogeneous and subclonal, occurring both exclusively in precursor or cancer populations in different cases. Whole-transcriptome profiling of coisolated RNA from 12 lesions (from 5 cases) was robust and confirmed histologic and molecular heterogeneity, including activated Wnt signaling in CTNNB1-mutant versus wild-type populations. Taken together, we demonstrate clinically relevant multiclonality and intratumoral heterogeneity during LGEC development with important implications for diagnosis, prognosis, and therapeutic prediction. More broadly, our methodology is broadly scalable to enable high-throughput genomic and transcriptomic characterization of precursor and invasive cancer populations from routine FFPE specimens. IMPLICATIONS: Multiregion profiling of LGEC populations using a highly scalable approach demonstrates clinically relevant multiclonality and intratumoral heterogeneity.

Circulating Tumor Cell-Based Molecular Classifier for Predicting Resistance to Abiraterone and Enzalutamide in Metastatic Castration-Resistant Prostate Cancer.

While circulating tumor cell (CTC)-based detection of AR-V7 has been demonstrated to predict patient response to second-generation androgen receptor therapies, the rarity of AR-V7 expression in metastatic castrate-resistant prostate cancer (mCRPC) suggests that other drivers of resistance exist. We sought to use a multiplex gene expression platform to interrogate CTCs and identify potential markers of resistance to abiraterone and enzalutamide. 37 patients with mCRPC initiating treatment with enzalutamide (n = 16) or abiraterone (n = 21) were prospectively enrolled for CTC collection and gene expression analysis using a panel of 89 prostate cancer-related genes. Gene expression from CTCs was correlated with PSA response and radioclinical progression-free survival (PFS) using Kaplan-Meier and Cox regression analyses. Twenty patients (54%) had detectable CTCs. At a median follow-up of 11.3 months, increased expression of the following genes was significantly associated with shorter PSA PFS and radioclinical PFS: AR, AR-V7, PSA, PSCA, TSPAN8, NKX3.1, and WNT5B. Additionally, high SPINK1 expression was associated with increased PFS. A predictive model including all eight genes gave an area under the curve (AUC) of 0.84 for PSA PFS and 0.86 for radioclinical PFS. In comparison, the AR-V7 only model resulted in AUC values of 0.65 and 0.64.These data demonstrate that clinically relevant information regarding gene expression can be obtained from whole blood using a CTC-based approach. Multigene classifiers in this setting may allow for the development of noninvasive predictive biomarkers to guide clinical management.

Somatic mutations in adrenocortical carcinoma with primary aldosteronism or hyperreninemic hyperaldosteronism.

Several somatic mutations specific to aldosterone-producing adenomas (APAs) have been described. A small proportion of adrenocortical carcinomas (ACCs) are associated with hyperaldosteronism, either primary aldosteronism or hyperreninemic hyperaldosteronism. However, it is unknown whether they harbor mutations of the same spectrum as APAs. The objective of this study is to describe the clinical phenotype and molecular genotype of ACCs with hyperaldosteronism, particularly the analysis for common APA-associated genetic changes. Patients were identified by retrospective chart review at a specialized referral center and by positive staining for CYP11B2 of tissue microarrays. Twenty-five patients with ACC and hyperaldosteronism were initially identified by retrospective chart review, and tissue for further analysis was available on 13 tumors. Seven patients were identified by positive staining for CYP11B2 in a tissue microarray, of which two were already identified in the initial chart review. Therefore, a total number of 18 patients with a diagnosis of ACC and features of either primary aldosteronism or hyperreninemic hyperaldosteronism were therefore included in the final study. Mutational status for a select list of oncogenes, tumor suppressor genes and genes known to carry mutations in APAs were analyzed by next-generation sequencing. Review of clinical data suggested autonomous aldosterone production in the majority of cases, while for some cases, hyperreninemic hyperaldosteronism was the more likely mechanism. The mutational landscape of ACCs associated with hyperaldosteronism was not different from ACCs with a different hormonal phenotype. None of the ACCs harbored mutations of known APA-associated genes, suggesting an alternative mechanism conferring aldosterone production.

Biologic Significance of Magnetic Resonance Imaging Invisibility in Localized Prostate Cancer.

PURPOSE: Multiparametric magnetic resonance imaging (mpMRI) is used widely for prostate cancer (PCa) evaluation. Approximately 35% of aggressive tumors, however, are not visible on mpMRI. We sought to identify the molecular alterations associated with mpMRI-invisible tumors and determine whether mpMRI visibility is associated with PCa prognosis. METHODS: Discovery and validation cohorts included patients who underwent mpMRI before radical prostatectomy and were found to harbor both mpMRI-visible (Prostate Imaging and Reporting Data System 3 to 5) and -invisible (Prostate Imaging and Reporting Data System 1 or 2) foci on surgical pathology. Next-generation sequencing was performed to determine differential gene expression between mpMRI-visible and -invisible foci. A genetic signature for tumor mpMRI visibility was derived in the discovery cohort and assessed in an independent validation cohort. Its association with long-term oncologic outcomes was evaluated in a separate testing cohort. RESULTS: The discovery cohort included 10 patients with 26 distinct PCa foci on surgical pathology, of which 12 (46%) were visible and 14 (54%) were invisible on preoperative mpMRI. Next-generation sequencing detected prioritized genetic mutations in 14 (54%) tumor foci (n = 8 mpMRI visible, n = 6 mpMRI invisible). A nine-gene signature (composed largely of cell organization/structure genes) associated with mpMRI visibility was derived (area under the curve = 0.89), and the signature predicted MRI visibility with 75% sensitivity and 100% specificity (area under the curve = 0.88) in the validation cohort. In the testing cohort (n = 375, median follow-up 8 years) there was no significant difference in biochemical recurrence, distant metastasis, or cancer-specific mortality in patients with predicted mpMRI-visible versus -invisible tumors (all P > .05). CONCLUSION: Compared with mpMRI-invisible disease, mpMRI-visible tumors are associated with underexpression of cellular organization genes. mpMRI visibility does not seem to be predictive of long-term cancer outcomes, highlighting the need for biopsy strategies that detect mpMRI-invisible tumors.