Plasma Kidney Injury Molecule-1 for Preoperative Prediction of Renal Cell Carcinoma Versus Benign Renal Masses, and Association With Clinical Outcomes

PURPOSE

Both clear cell and papillary renal cell carcinomas (RCCs) overexpress kidney injury molecule-1 (KIM-1). We investigated whether plasma KIM-1 (pKIM-1) may be a useful risk stratification tool among patients with suspicious renal masses.

METHODS

Prenephrectomy pKIM-1 was measured in two independent cohorts of patients with renal masses. Cohort 1, from the prospective K2 trial, included 162 patients found to have clear cell RCC (cases) and 162 patients with benign renal masses (controls). Cohort 2 included 247 patients with small (cT1a) renal masses from an academic biorepository, of whom 184 had RCC. We assessed the relationship between pKIM-1, surgical pathology, and clinical outcomes.

RESULTS

In Cohort 1, pKIM-1 distinguished RCC versus benign masses with area under the receiver operating curve (AUC-ROC, 0.81 [95% CI, 0.76 to 0.86]). In Cohort 2 (cT1a only), pKIM-1 distinguished RCC versus benign masses (AUC-ROC, 0.74 [95% CI, 0.67 to 0.80]) and the addition of pKIM-1 to an established nomogram for predicting malignancy improved the model AUC-ROC (0.65 [95% CI, 0.57 to 0.74] v 0.78 [95% CI, 0.72 to 0.85]). A pKIM-1 cutpoint identified using Cohort 2 demonstrated sensitivity of 92.5% and specificity of 60% for identifying RCC in Cohort 1. In long-term follow-up of RCC cases (Cohort 1), higher prenephrectomy pKIM-1 was associated with worse metastasis-free survival (multivariable MFS hazard ratio [HR] 1.29 per unit increase in log pKIM-1, 95% CI, 1.10 to 1.53) and overall survival (multivariable OS HR 1.31 per unit increase in log pKIM-1, 95% CI, 1.10 to 1.54). In long-term follow-up of Cohort 2, no metastatic events occurred, consistent with the favorable prognosis of resected cT1a RCC.

CONCLUSION

Among patients with renal masses, pKIM-1 is associated with malignant pathology, worse MFS, and risk of death. pKIM-1 may be useful for selecting patients with renal masses for intervention versus surveillance.

Multi-ancestry genome-wide association study of kidney cancer identifies 63 susceptibility regions

Here, in a multi-ancestry genome-wide association study meta-analysis of kidney cancer (29,020 cases and 835,670 controls), we identified 63 susceptibility regions (50 novel) containing 108 independent risk loci. In analyses stratified by subtype, 52 regions (78 loci) were associated with clear cell renal cell carcinoma (RCC) and 6 regions (7 loci) with papillary RCC. Notably, we report a variant common in African ancestry individuals ( rs7629500 ) in the 3' untranslated region of VHL, nearly tripling clear cell RCC risk (odds ratio 2.72, 95% confidence interval 2.23-3.30). In cis-expression quantitative trait locus analyses, 48 variants from 34 regions point toward 83 candidate genes. Enrichment of hypoxia-inducible factor-binding sites underscores the importance of hypoxia-related mechanisms in kidney cancer. Our results advance understanding of the genetic architecture of kidney cancer, provide clues for functional investigation and enable generation of a validated polygenic risk score with an estimated area under the curve of 0.65 (0.74 including risk factors) among European ancestry individuals.

Increased chromatin accessibility drives transition to androgen receptor splice variant dependence in castration-resistant prostate cancer

Androgen receptor (AR) splice variants, of which ARv7 is the most common, are increased in prostate cancer (PC) that develops resistance to androgen signaling inhibitor drugs, but the extent to which these variants drive AR activity, and whether they have novel functions or dependencies, remain to be determined. We generated a subline of VCaP PC cells (VCaP16) that is resistant to the AR inhibitor enzalutamide (ENZ) and found that AR activity was independent of the full-length AR (ARfl), despite its continued high-level expression, and was instead driven by ARv7. The ARv7 cistrome and transcriptome in VCaP16 cells mirrored that of the ARfl in VCaP cells, although ARv7 chromatin binding was weaker, and strong ARv7 binding sites correlated with higher affinity ARfl binding sites across multiple models and clinical samples. Notably, although ARv7 expression in VCaP cells increased rapidly in response to ENZ, there was a long lag before it gained chromatin binding and transcriptional activity. This lag was associated with an increase in chromatin accessibility, with the AR and nuclear factor I (NFI) motifs being most enriched at these more accessible sites. Moreover, the transcriptional effects of combined NFIB and NFIX knockdown versus ARv7 knockdown were highly correlated. These findings indicate that ARv7 can drive the AR program, but that its activity is dependent on adaptations that increase chromatin accessibility to enhance its intrinsically weak chromatin binding.

The PENGUIN approach to reconstruct protein interactions at enhancer-promoter regions and its application to prostate cancer

We introduce Promoter-Enhancer-Guided Interaction Networks (PENGUIN), a method for studying protein-protein interaction (PPI) networks within enhancer-promoter interactions. PENGUIN integrates H3K27ac-HiChIP data with tissue-specific PPIs to define enhancer-promoter PPI networks (EPINs). We validated PENGUIN using cancer (LNCaP) and benign (LHSAR) prostate cell lines. Our analysis detected EPIN clusters enriched with the architectural protein CTCF, a regulator of enhancer-promoter interactions. CTCF presence was coupled with the prevalence of prostate cancer (PrCa) single nucleotide polymorphisms (SNPs) within the same EPIN clusters, suggesting functional implications in PrCa. Within the EPINs displaying enrichments in both CTCF and PrCa SNPs, we also show enrichment in oncogenes. We substantiated our identified SNPs through CRISPR/Cas9 knockout and RNAi screens experiments. Here we show that PENGUIN provides insights into the intricate interplay between enhancer-promoter interactions and PPI networks, which are crucial for identifying key genes and potential intervention targets. A dedicated server is available at https://penguin.life.bsc.es/ .

Liquid biopsy epigenomic profiling for cancer subtyping

Although circulating tumor DNA (ctDNA) assays are increasingly used to inform clinical decisions in cancer care, they have limited ability to identify the transcriptional programs that govern cancer phenotypes and their dynamic changes during the course of disease. To address these limitations, we developed a method for comprehensive epigenomic profiling of cancer from 1 ml of patient plasma. Using an immunoprecipitation-based approach targeting histone modifications and DNA methylation, we measured 1,268 epigenomic profiles in plasma from 433 individuals with one of 15 cancers. Our assay provided a robust proxy for transcriptional activity, allowing us to infer the expression levels of diagnostic markers and drug targets, measure the activity of therapeutically targetable transcription factors and detect epigenetic mechanisms of resistance. This proof-of-concept study in advanced cancers shows how plasma epigenomic profiling has the potential to unlock clinically actionable information that is currently accessible only via direct tissue sampling.

Grade Group 1 Prostate Cancers Exhibit Tumor-defining Androgen Receptor-driven Programs

Grade group 1 (GG1) primary prostate cancers with a pathologic Gleason score of 6 are considered indolent and generally not associated with fatal outcomes, so treatment is not indicated for most cases. These low-grade cancers have an overall negligible risk of locoregional progression and metastasis to distant organs, which is why there is an ongoing debate about whether these lesions should be reclassified as "noncancerous". However, the underlying molecular activity of key disease drivers, such as the androgen receptor (AR), have thus far not been thoroughly characterized in low-grade tumors. Therefore, we set out to delineate the AR chromatin-binding landscape in low-grade GG1 prostate cancers to gain insights into whether these AR-driven programs are actually tumor-specific or are normal prostate epithelium-like. These analyses showed that GG1 tumors do not harbor a distinct AR cistrome and, similar to higher-grade cancers, AR preferentially binds to tumor-defining cis-regulatory elements. Furthermore, the enhancer activity of these regions and the expression of their respective target genes were not significantly different in GG1 tumors. From an epigenetic perspective, this finding supports the cancer designation currently given to these low-grade tumors and clearly distinguishes them from noncancerous benign tissue. PATIENT SUMMARY: We characterized the molecular activity of the androgen receptor protein, which drives prostate cancer disease, in low-grade tumors. Our results show that these tumors are true cancers and are clearly separate from benign prostate tissue despite their low clinical aggressiveness.

A biallelic multiple nucleotide length polymorphism explains functional causality at 5p15.33 prostate cancer risk locus

To date, single-nucleotide polymorphisms (SNPs) have been the most intensively investigated class of polymorphisms in genome wide associations studies (GWAS), however, other classes such as insertion-deletion or multiple nucleotide length polymorphism (MNLPs) may also confer disease risk. Multiple reports have shown that the 5p15.33 prostate cancer risk region is a particularly strong expression quantitative trait locus (eQTL) for Iroquois Homeobox 4 (IRX4) transcripts. Here, we demonstrate using epigenome and genome editing that a biallelic (21 and 47 base pairs (bp)) MNLP is the causal variant regulating IRX4 transcript levels. In LNCaP prostate cancer cells (homozygous for the 21 bp short allele), a single copy knock-in of the 47 bp long allele potently alters the chromatin state, enabling de novo functional binding of the androgen receptor (AR) associated with increased chromatin accessibility, Histone 3 lysine 27 acetylation (H3K27ac), and ~3-fold upregulation of IRX4 expression. We further show that an MNLP is amongst the strongest candidate susceptibility variants at two additional prostate cancer risk loci. We estimated that at least 5% of prostate cancer risk loci could be explained by functional non-SNP causal variants, which may have broader implications for other cancers GWAS. More generally, our results underscore the importance of investigating other classes of inherited variation as causal mediators of human traits.

Rewiring of master transcription factor cistromes during high-grade serous ovarian cancer development

The transcription factors MECOM, PAX8, SOX17 and WT1 are candidate master regulators of high-grade serous 'ovarian' cancer (HGSC), yet their cooperative role in the hypothesized tissue of origin, the fallopian tube secretory epithelium (FTSEC) is unknown. We generated 26 epigenome (CUT&TAG, CUT&RUN, ATAC-seq and HiC) data sets and 24 profiles of RNA-seq transcription factor knock-down followed by RNA sequencing in FTSEC and HGSC models to define binding sites and gene sets regulated by these factors in cis and trans. This revealed that MECOM, PAX8, SOX17 and WT1 are lineage-enriched, super-enhancer associated master regulators whose cooperative DNA-binding patterns and target genes are re-wired during tumor development. All four TFs were indispensable for HGSC clonogenicity and survival but only depletion of PAX8 and WT1 impaired FTSEC cell survival. These four TFs were pharmacologically inhibited by transcriptional inhibitors only in HGSCs but not in FTSECs. Collectively, our data highlights that tumor-specific epigenetic remodeling is tightly related to MECOM, PAX8, SOX17 and WT1 activity and these transcription factors are targetable in a tumor-specific manner through transcriptional inhibitors.

Analysis of evolutionary dynamics and clonal architecture in prostate cancer

The extent to which clinical and genomic characteristics associate with prostate cancer clonal architecture, tumor evolution, and therapeutic response remains unclear. Here, we reconstructed the clonal architecture and evolutionary trajectories of 845 prostate cancer tumors with harmonized clinical and molecular data. We observed that tumors from patients who self-reported as Black had more linear and monoclonal architectures, despite these men having higher rates of biochemical recurrence. This finding contrasts with prior observations relating polyclonal architecture to adverse clinical outcomes. Additionally, we utilized a novel approach to mutational signature analysis that leverages clonal architecture to uncover additional cases of homologous recombination and mismatch repair deficiency in primary and metastatic tumors and link the origin of mutational signatures to specific subclones. Broadly, prostate cancer clonal architecture analysis reveals novel biological insights that may be immediately clinically actionable and provide multiple opportunities for subsequent investigation.

Statement of significance

Tumors from patients who self-reported as Black demonstrate linear and monoclonal evolutionary trajectories yet experience higher rates of biochemical recurrence. In addition, analysis of clonal and subclonal mutational signatures identifies additional tumors with potentially actionable alterations such as deficiencies in mismatch repair and homologous recombination.

Nucleosome Patterns in Circulating Tumor DNA Reveal Transcriptional Regulation of Advanced Prostate Cancer Phenotypes

Advanced prostate cancers comprise distinct phenotypes, but tumor classification remains clinically challenging. Here, we harnessed circulating tumor DNA (ctDNA) to study tumor phenotypes by ascertaining nucleosome positioning patterns associated with transcription regulation. We sequenced plasma ctDNA whole genomes from patient-derived xenografts representing a spectrum of androgen receptor active (ARPC) and neuroendocrine (NEPC) prostate cancers. Nucleosome patterns associated with transcriptional activity were reflected in ctDNA at regions of genes, promoters, histone modifications, transcription factor binding, and accessible chromatin. We identified the activity of key phenotype-defining transcriptional regulators from ctDNA, including AR, ASCL1, HOXB13, HNF4G, and GATA2. To distinguish NEPC and ARPC in patient plasma samples, we developed prediction models that achieved accuracies of 97% for dominant phenotypes and 87% for mixed clinical phenotypes. Although phenotype classification is typically assessed by IHC or transcriptome profiling from tumor biopsies, we demonstrate that ctDNA provides comparable results with diagnostic advantages for precision oncology.

SIGNIFICANCE

This study provides insights into the dynamics of nucleosome positioning and gene regulation associated with cancer phenotypes that can be ascertained from ctDNA. New methods for classification in phenotype mixtures extend the utility of ctDNA beyond assessments of somatic DNA alterations with important implications for molecular classification and precision oncology. This article is highlighted in the In This Issue feature, p. 517.

Enhancer profiling identifies epigenetic markers of endocrine resistance and reveals therapeutic options for metastatic castration-resistant prostate cancer patients

Androgen Receptor (AR) signaling inhibitors, including enzalutamide, are treatment options for patients with metastatic castration-resistant prostate cancer (mCRPC), but resistance inevitably develops. Using metastatic samples from a prospective phase II clinical trial, we epigenetically profiled enhancer/promoter activities with H3K27ac chromatin immunoprecipitation followed by sequencing, before and after AR-targeted therapy. We identified a distinct subset of H3K27ac-differentially marked regions that associated with treatment responsiveness. These data were successfully validated in mCRPC patient-derived xenograft models (PDX). In silico analyses revealed HDAC3 as a critical factor that can drive resistance to hormonal interventions, which we validated in vitro . Using cell lines and mCRPC PDX tumors in vitro , we identified drug-drug synergy between enzalutamide and the pan-HDAC inhibitor vorinostat, providing therapeutic proof-of-concept. These findings demonstrate rationale for new therapeutic strategies using a combination of AR and HDAC inhibitors to improve patient outcome in advanced stages of mCRPC.

Epigenomic charting and functional annotation of risk loci in renal cell carcinoma

While the mutational and transcriptional landscapes of renal cell carcinoma (RCC) are well-known, the epigenome is poorly understood. We characterize the epigenome of clear cell (ccRCC), papillary (pRCC), and chromophobe RCC (chRCC) by using ChIP-seq, ATAC-Seq, RNA-seq, and SNP arrays. We integrate 153 individual data sets from 42 patients and nominate 50 histology-specific master transcription factors (MTF) to define RCC histologic subtypes, including EPAS1 and ETS-1 in ccRCC, HNF1B in pRCC, and FOXI1 in chRCC. We confirm histology-specific MTFs via immunohistochemistry including a ccRCC-specific TF, BHLHE41. FOXI1 overexpression with knock-down of EPAS1 in the 786-O ccRCC cell line induces transcriptional upregulation of chRCC-specific genes, TFCP2L1, ATP6V0D2, KIT, and INSRR, implicating FOXI1 as a MTF for chRCC. Integrating RCC GWAS risk SNPs with H3K27ac ChIP-seq and ATAC-seq data reveals that risk-variants are significantly enriched in allelically-imbalanced peaks. This epigenomic atlas in primary human samples provides a resource for future investigation.

Germline variants associated with toxicity to immune checkpoint blockade

Immune checkpoint inhibitors (ICIs) have yielded remarkable responses but often lead to immune-related adverse events (irAEs). Although germline causes for irAEs have been hypothesized, no individual variant associated with developing irAEs has been identified. We carried out a genome-wide association study of 1,751 patients on ICIs across 12 cancer types. We investigated two irAE phenotypes: (1) high-grade (3-5) and (2) all-grade events. We identified 3 genome-wide significant associations (P < 5 × 10-8) in the discovery cohort associated with all-grade irAEs: rs16906115 near IL7 (combined P = 3.6 × 10-11; hazard ratio (HR) = 2.1); rs75824728 near IL22RA1 (combined P = 3.5 × 10-8; HR = 1.8); and rs113861051 on 4p15 (combined P = 1.2 × 10-8, HR = 2.0); rs16906115 was replicated in 3 independent studies. The association near IL7 colocalized with the gain of a new cryptic exon for IL7, a critical regulator of lymphocyte homeostasis. Patients carrying the IL7 germline variant exhibited significantly increased lymphocyte stability after ICI initiation, which was itself predictive of downstream irAEs and improved survival.

IL7 genetic variation and toxicity to immune checkpoint blockade in patients with melanoma

Treatment with immune checkpoint blockade (ICB) frequently triggers immune-related adverse events (irAEs), causing considerable morbidity. In 214 patients receiving ICB for melanoma, we observed increased severe irAE risk in minor allele carriers of rs16906115, intronic to IL7. We found that rs16906115 forms a B cell-specific expression quantitative trait locus (eQTL) to IL7 in patients. Patients carrying the risk allele demonstrate increased pre-treatment B cell IL7 expression, which independently associates with irAE risk, divergent immunoglobulin expression and more B cell receptor mutations. Consistent with the role of IL-7 in T cell development, risk allele carriers have distinct ICB-induced CD8+ T cell subset responses, skewing of T cell clonality and greater proportional repertoire occupancy by large clones. Finally, analysis of TCGA data suggests that risk allele carriers independently have improved melanoma survival. These observations highlight key roles for B cells and IL-7 in both ICB response and toxicity and clinical outcomes in melanoma.

Extensive androgen receptor enhancer heterogeneity in primary prostate cancers underlies transcriptional diversity and metastatic potential

Androgen receptor (AR) drives prostate cancer (PCa) development and progression. AR chromatin binding profiles are highly plastic and form recurrent programmatic changes that differentiate disease stages, subtypes and patient outcomes. While prior studies focused on concordance between patient subgroups, inter-tumor heterogeneity of AR enhancer selectivity remains unexplored. Here we report high levels of AR chromatin binding heterogeneity in human primary prostate tumors, that overlap with heterogeneity observed in healthy prostate epithelium. Such heterogeneity has functional consequences, as somatic mutations converge on commonly-shared AR sites in primary over metastatic tissues. In contrast, less-frequently shared AR sites associate strongly with AR-driven gene expression, while such heterogeneous AR enhancer usage also distinguishes patients' outcome. These findings indicate that epigenetic heterogeneity in primary disease is directly informative for risk of biochemical relapse. Cumulatively, our results illustrate a high level of AR enhancer heterogeneity in primary PCa driving differential expression and clinical impact.

Copy Number Variants Are Ovarian Cancer Risk Alleles at Known and Novel Risk Loci

BACKGROUND

Known risk alleles for epithelial ovarian cancer (EOC) account for approximately 40% of the heritability for EOC. Copy number variants (CNVs) have not been investigated as EOC risk alleles in a large population cohort.

METHODS

Single nucleotide polymorphism array data from 13 071 EOC cases and 17 306 controls of White European ancestry were used to identify CNVs associated with EOC risk using a rare admixture maximum likelihood test for gene burden and a by-probe ratio test. We performed enrichment analysis of CNVs at known EOC risk loci and functional biofeatures in ovarian cancer-related cell types.

RESULTS

We identified statistically significant risk associations with CNVs at known EOC risk genes; BRCA1 (PEOC = 1.60E-21; OREOC = 8.24), RAD51C (Phigh-grade serous ovarian cancer [HGSOC] = 5.5E-4; odds ratio [OR]HGSOC = 5.74 del), and BRCA2 (PHGSOC = 7.0E-4; ORHGSOC = 3.31 deletion). Four suggestive associations (P < .001) were identified for rare CNVs. Risk-associated CNVs were enriched (P < .05) at known EOC risk loci identified by genome-wide association study. Noncoding CNVs were enriched in active promoters and insulators in EOC-related cell types.

CONCLUSIONS

CNVs in BRCA1 have been previously reported in smaller studies, but their observed frequency in this large population-based cohort, along with the CNVs observed at BRCA2 and RAD51C gene loci in EOC cases, suggests that these CNVs are potentially pathogenic and may contribute to the spectrum of disease-causing mutations in these genes. CNVs are likely to occur in a wider set of susceptibility regions, with potential implications for clinical genetic testing and disease prevention.

A Distinct Chromatin State Drives Therapeutic Resistance in Invasive Lobular Breast Cancer

Most invasive lobular breast cancers (ILC) are of the luminal A subtype and are strongly hormone receptor-positive. Yet, ILC is relatively resistant to tamoxifen and associated with inferior long-term outcomes compared with invasive ductal cancers (IDC). In this study, we sought to gain mechanistic insights into these clinical findings that are not explained by the genetic landscape of ILC and to identify strategies to improve patient outcomes. A comprehensive analysis of the epigenome of ILC in preclinical models and clinical samples showed that, compared with IDC, ILC harbored a distinct chromatin state linked to gained recruitment of FOXA1, a lineage-defining pioneer transcription factor. This resulted in an ILC-unique FOXA1-estrogen receptor (ER) axis that promoted the transcription of genes associated with tumor progression and poor outcomes. The ILC-unique FOXA1-ER axis led to retained ER chromatin binding after tamoxifen treatment, which facilitated tamoxifen resistance while remaining strongly dependent on ER signaling. Mechanistically, gained FOXA1 binding was associated with the autoinduction of FOXA1 in ILC through an ILC-unique FOXA1 binding site. Targeted silencing of this regulatory site resulted in the disruption of the feed-forward loop and growth inhibition in ILC. In summary, ILC is characterized by a unique chromatin state and FOXA1-ER axis that is associated with tumor progression, offering a novel mechanism of tamoxifen resistance. These results underscore the importance of conducting clinical trials dedicated to patients with ILC in order to optimize treatments in this breast cancer subtype.

SIGNIFICANCE

A unique FOXA1-ER axis in invasive lobular breast cancer promotes disease progression and tamoxifen resistance, highlighting a potential therapeutic avenue for clinical investigations dedicated to this disease. See related commentary by Blawski and Toska, p. 3668.

Genetic determinants of chromatin reveal prostate cancer risk mediated by context-dependent gene regulation

Many genetic variants affect disease risk by altering context-dependent gene regulation. Such variants are difficult to study mechanistically using current methods that link genetic variation to steady-state gene expression levels, such as expression quantitative trait loci (eQTLs). To address this challenge, we developed the cistrome-wide association study (CWAS), a framework for identifying genotypic and allele-specific effects on chromatin that are also associated with disease. In prostate cancer, CWAS identified regulatory elements and androgen receptor-binding sites that explained the association at 52 of 98 known prostate cancer risk loci and discovered 17 additional risk loci. CWAS implicated key developmental transcription factors in prostate cancer risk that are overlooked by eQTL-based approaches due to context-dependent gene regulation. We experimentally validated associations and demonstrated the extensibility of CWAS to additional epigenomic datasets and phenotypes, including response to prostate cancer treatment. CWAS is a powerful and biologically interpretable paradigm for studying variants that influence traits by affecting transcriptional regulation.

The Prostate Cancer Androgen Receptor Cistrome in African American Men Associates with Upregulation of Lipid Metabolism and Immune Response

African-American (AA) men are more likely to be diagnosed with and die from prostate cancer than European American (EA) men. Despite the central role of the androgen receptor (AR) transcription factor in prostate cancer, little is known about the contribution of epigenetics to observed racial disparities. We performed AR chromatin immunoprecipitation sequencing on primary prostate tumors from AA and EA men, finding that sites with greater AR binding intensity in AA relative to EA prostate cancer are enriched for lipid metabolism and immune response genes. Integration with transcriptomic and metabolomic data demonstrated coinciding upregulation of lipid metabolism gene expression and increased lipid levels in AA prostate cancer. In a metastatic prostate cancer cohort, upregulated lipid metabolism associated with poor prognosis. These findings offer the first insights into ancestry-specific differences in the prostate cancer AR cistrome. The data suggest a model whereby increased androgen signaling may contribute to higher levels of lipid metabolism, immune response, and cytokine signaling in AA prostate tumors. Given the association of upregulated lipogenesis with prostate cancer progression, our study provides a plausible biological explanation for the higher incidence and aggressiveness of prostate cancer observed in AA men.

SIGNIFICANCE

With immunotherapies and inhibitors of metabolic enzymes in clinical development, the altered lipid metabolism and immune response in African-American men provides potential therapeutic opportunities to attenuate racial disparities in prostate cancer.

Addition of Germline Testing to Tumor-Only Sequencing Improves Detection of Pathogenic Germline Variants in Men With Advanced Prostate Cancer

PURPOSE

Guidelines recommend somatic and germline testing for men with advanced prostate cancer (PCa). Barriers to widespread implementation result in underutilization of germline testing. Somatic testing alone risks missing pathogenic germline variants (PGVs). We sought to determine whether the addition of germline testing to tumor-only sequencing improves detection of PGVs in men with advanced PCa. Secondarily, we sought to define the added value of combining somatic and germline testing to optimize detection of clinically actionable alterations.

PATIENTS AND METHODS

We analyzed results of independent germline testing and tumor-only sequencing from 100 men with advanced PCa from a prospective clinical trial (ClinicalTrials.gov identifier: NCT03328091). The primary outcome was the proportion of PGVs not reported with tumor-only sequencing. The secondary outcome was the association of locus-specific loss of heterozygosity for PGVs in homologous recombination genes with clinical-genomic features.

RESULTS

In the 100 men who underwent germline testing and tumor-only sequencing, 24 PGVs were identified, 17 of which were clinically actionable, in 23 patients. Tumor-only sequencing failed to report four (17%) of the PGVs. One additional PGV (4.2%) had variant allele frequency on tumor-sequencing below the threshold for follow-up germline testing. When integrating tumor-only sequencing with germling testing results, 33% of patients harbored clinically actionable alterations. Rates of locus-specific loss of heterozygosity were higher for BRCA2 PGVs in castration-resistant PCa than PGVs in other homologous recombination genes in hormone-sensitive PCa (P = .029).

CONCLUSION

Tumor-only sequencing failed to report more than 20% of PGVs in men with advanced PCa. These findings strongly support guideline recommendations for universal germline and somatic testing in this population. Combining tumor and germline sequencing doubled the chance of detecting a clinically actionable alteration.

Drug-induced epigenomic plasticity reprograms circadian rhythm regulation to drive prostate cancer towards androgen-independence

In prostate cancer, androgen receptor (AR)-targeting agents are very effective in various disease stages. However, therapy resistance inevitably occurs and little is known about how tumor cells adapt to bypass AR suppression. Here, we performed integrative multi-omics analyses on tissues isolated before and after 3 months of AR-targeting enzalutamide monotherapy from high-risk prostate cancer patients enrolled in a neoadjuvant clinical trial. Transcriptomic analyses demonstrated that AR inhibition drove tumors towards a neuroendocrine-like disease state. Additionally, epigenomic profiling revealed massive enzalutamide-induced reprogramming of pioneer factor FOXA1 - from inactive chromatin sites towards active cis-regulatory elements that dictate pro-survival signals. Notably, treatment-induced FOXA1 sites were enriched for circadian clock component ARNTL. Post-treatment ARNTL levels associated with poor outcome, and ARNTL knockout strongly decreased prostate cancer cell growth. Our data highlight a remarkable cistromic plasticity of FOXA1 following AR-targeted therapy, and revealed an acquired dependency on circadian regulator ARNTL, a novel candidate therapeutic target.

Response to supraphysiological testosterone is predicted by a distinct androgen receptor cistrome

The androgen receptor (AR) is a master transcription factor that regulates prostate cancer (PC) development and progression. Inhibition of AR signaling by androgen deprivation is the first-line therapy with initial efficacy for advanced and recurrent PC. Paradoxically, supraphysiological levels of testosterone (SPT) also inhibit PC progression. However, as with any therapy, not all patients show a therapeutic benefit, and responses differ widely in magnitude and duration. In this study, we evaluated whether differences in the AR cistrome before treatment can distinguish between SPT-responding (R) and -nonresponding (NR) tumors. We provide the first preclinical evidence to our knowledge that SPT-R tumors exhibit a distinct AR cistrome when compared with SPT-NR tumors, indicating a differential biological role of the AR. We applied an integrated analysis of ChIP-Seq and RNA-Seq to the pretreatment tumors and identified an SPT-R signature that distinguishes R and NR tumors. Because transcriptomes of SPT-treated clinical specimens are not available, we interrogated available castration-resistant PC (CRPC) transcriptomes and showed that the SPT-R signature is associated with improved survival and has the potential to identify patients who would respond to SPT. These findings provide an opportunity to identify the subset of patients with CRPC who would benefit from SPT therapy.

MYC drives aggressive prostate cancer by disrupting transcriptional pause release at androgen receptor targets

c-MYC (MYC) is a major driver of prostate cancer tumorigenesis and progression. Although MYC is overexpressed in both early and metastatic disease and associated with poor survival, its impact on prostate transcriptional reprogramming remains elusive. We demonstrate that MYC overexpression significantly diminishes the androgen receptor (AR) transcriptional program (the set of genes directly targeted by the AR protein) in luminal prostate cells without altering AR expression. Analyses of clinical specimens reveal that concurrent low AR and high MYC transcriptional programs accelerate prostate cancer progression toward a metastatic, castration-resistant disease. Data integration of single-cell transcriptomics together with ChIP-seq uncover an increase in RNA polymerase II (Pol II) promoter-proximal pausing at AR-dependent genes following MYC overexpression without an accompanying deactivation of AR-bound enhancers. Altogether, our findings suggest that MYC overexpression antagonizes the canonical AR transcriptional program and contributes to prostate tumor initiation and progression by disrupting transcriptional pause release at AR-regulated genes.

CREB5 reprograms FOXA1 nuclear interactions to promote resistance to androgen receptor targeting therapies

Metastatic castration resistant prostate cancers (mCRPC) are treated with therapies that antagonize the androgen receptor (AR). Nearly all patients develop resistance to AR-targeted therapies (ART). Our previous work identified CREB5 as an upregulated target gene in human mCRPC that promoted resistance to all clinically-approved ART. The mechanisms by which CREB5 promotes progression of mCRPC or other cancers remains elusive. Integrating ChIP-seq and rapid immunoprecipitation and mass spectroscopy of endogenous proteins (RIME), we report that cells overexpressing CREB5 demonstrate extensive reprogramming of nuclear protein-protein interactions in response to the ART agent enzalutamide. Specifically, CREB5 physically interacts with AR, the pioneering actor FOXA1, and other known co-factors of AR and FOXA1 at transcription regulatory elements recently found to be active in mCRPC patients. We identified a subset of CREB5/FOXA1 co-interacting nuclear factors that have critical functions for AR transcription (GRHL2, HOXB13) while others (TBX3, NFIC) regulated cell viability and ART resistance and were amplified or overexpressed in mCRPC. Upon examining the nuclear protein interactions and the impact of CREB5 expression on the mCRPC patient transcriptome, we found CREB5 was associated with Wnt signaling and epithelial to mesenchymal transitions, implicating these pathways in CREB5/FOXA1-mediated ART resistance. Overall, these observations define the molecular interactions among CREB5, FOXA1, and pathways that promote ART resistance.

Androgen receptor reprogramming demarcates prognostic, context-dependent gene sets in primary and metastatic prostate cancer

The androgen receptor (AR) is a prostate master transcription factor. It binds to genetic enhancers, where it regulates gene activity and plays a fundamental role in prostate pathophysiology. Previous work has demonstrated that AR-DNA binding is systematically and consistently reprogrammed during prostate tumorigenesis and disease progression. We charted these reprogrammed AR sites and identified genes proximal to them. We were able to devise gene lists based on AR status within specific histological contexts: normal prostate epithelium, primary prostate tumor, and metastatic prostate cancer. We evaluated expression of the genes in these gene sets in subjects from two distinct clinical cohorts-men treated with surgery for localized prostate cancer and men with metastatic prostate cancer. Among men with localized prostate cancer, expression of genes proximal to AR sites lost in the transition from normal prostate to prostate tumor was associated with clinical outcome. Among men with metastatic disease, expression of genes proximal to AR sites gained in metastatic tumors was associated with clinical outcome. These results are consistent with the notion that AR is fundamental to both maintaining differentiation in normal prostate tissue and driving de-differentiation in advanced prostate cancer. More broadly, the study demonstrates the power of incorporating context-dependent epigenetic data into genetic analyses.

HOXB13 suppresses de novo lipogenesis through HDAC3-mediated epigenetic reprogramming in prostate cancer

HOXB13, a homeodomain transcription factor, critically regulates androgen receptor (AR) activities and androgen-dependent prostate cancer (PCa) growth. However, its functions in AR-independent contexts remain elusive. Here we report HOXB13 interaction with histone deacetylase HDAC3, which is disrupted by the HOXB13 G84E mutation that has been associated with early-onset PCa. Independently of AR, HOXB13 recruits HDAC3 to lipogenic enhancers to catalyze histone deacetylation and suppress lipogenic regulators such as fatty acid synthase. Analysis of human tissues reveals that the HOXB13 gene is hypermethylated and downregulated in approximately 30% of metastatic castration-resistant PCa. HOXB13 loss or G84E mutation leads to lipid accumulation in PCa cells, thereby promoting cell motility and xenograft tumor metastasis, which is mitigated by pharmaceutical inhibition of fatty acid synthase. In summary, we present evidence that HOXB13 recruits HDAC3 to suppress de novo lipogenesis and inhibit tumor metastasis and that lipogenic pathway inhibitors may be useful to treat HOXB13-low PCa.

A multi-level investigation of the genetic relationship between endometriosis and ovarian cancer histotypes

Endometriosis is associated with increased risk of epithelial ovarian cancers (EOCs). Using data from large endometriosis and EOC genome-wide association meta-analyses, we estimate the genetic correlation and evaluate the causal relationship between genetic liability to endometriosis and EOC histotypes, and identify shared susceptibility loci. We estimate a significant genetic correlation (r_g) between endometriosis and clear cell (r_g = 0.71), endometrioid (r_g = 0.48), and high-grade serous (r_g = 0.19) ovarian cancer, associations supported by Mendelian randomization analyses. Bivariate meta-analysis identified 28 loci associated with both endometriosis and EOC, including 19 with evidence for a shared underlying association signal. Differences in the shared risk suggest different underlying pathways may contribute to the relationship between endometriosis and the different histotypes. Functional annotation using transcriptomic and epigenomic profiles of relevant tissues/cells highlights several target genes. This comprehensive analysis reveals profound genetic overlap between endometriosis and EOC histotypes with valuable genomic targets for understanding the biological mechanisms linking the diseases.

•

Endometriosis is genetically correlated with CCOC, ENOC, and HGSOC

•

Genetic liability to endometriosis confers risk of these EOC histotypes

•

Profound colocalization of genetic associations at endometriosis and EOC risk loci

•

Functional annotation highlights shared target genes elucidating the genetic link

A genome-scale CRISPR screen reveals PRMT1 as a critical regulator of androgen receptor signaling in prostate cancer

Androgen receptor (AR) signaling is the central driver of prostate cancer across disease states. While androgen deprivation therapy (ADT) is effective in the initial treatment of prostate cancer, resistance to ADT or to next-generation androgen pathway inhibitors invariably arises, most commonly through the re-activation of the AR axis. Thus, orthogonal approaches to inhibit AR signaling in advanced prostate cancer are essential. Here, via genome-scale CRISPR-Cas9 screening, we identify protein arginine methyltransferase 1 (PRMT1) as a critical mediator of AR expression and signaling. PRMT1 regulates the recruitment of AR to genomic target sites and the inhibition of PRMT1 impairs AR binding at lineage-specific enhancers, leading to decreased expression of key oncogenes, including AR itself. In addition, AR-driven prostate cancer cells are uniquely susceptible to combined AR and PRMT1 inhibition. Our findings implicate PRMT1 as a key regulator of AR output and provide a preclinical framework for co-targeting of AR and PRMT1 in advanced prostate cancer.

Detecting neuroendocrine prostate cancer through tissue-informed cell-free DNA methylation analysis

PURPOSE

Neuroendocrine prostate cancer (NEPC) is a resistance phenotype that emerges in men with metastatic castration-resistant prostate adenocarcinoma (CR-PRAD) and has important clinical implications, but is challenging to detect in practice. Herein, we report a novel tissue-informed epigenetic approach to non-invasively detect NEPC.

EXPERIMENTAL DESIGN

We first performed methylated immunoprecipitation and high-throughput sequencing (MeDIP-seq) on a training set of tumors, identified differentially methylated regions between NEPC and CR-PRAD, and built a model to predict the presence of NEPC (termed NEPC Risk Score). We then performed MeDIP-seq on cell-free DNA (cfDNA) from two independent cohorts of men with NEPC or CR-PRAD and assessed the accuracy of the model to predict the presence NEPC.

RESULTS

The test cohort comprised cfDNA samples from 48 men, 9 with NEPC and 39 with CR-PRAD. NEPC Risk Scores were significantly higher in men with NEPC than CR-PRAD (P=4.3x10^-7) and discriminated between NEPC and CR-PRAD with high accuracy (AUROC 0.96). The optimal NEPC Risk Score cut-off demonstrated 100% sensitivity and 90% specificity for detecting NEPC. The independent, multi-institutional validation cohort included cfDNA from 53 men, including 12 with NEPC and 41 with CR-PRAD. NEPC Risk Scores were significantly higher in men with NEPC than CR-PRAD (P=7.5x10^-12) and perfectly discriminated NEPC from CR-PRAD (AUROC 1.0). Applying the pre-defined NEPC Risk Score cut-off to the validation cohort resulted in 100% sensitivity and 95% specificity for detecting NEPC.

CONCLUSIONS

Tissue-informed cfDNA methylation analysis is a promising approach for non-invasive detection of NEPC in men with advanced prostate cancer.

H3K27ac HiChIP in prostate cell lines identifies risk genes for prostate cancer susceptibility

Genome-wide association studies (GWASs) have identified more than 200 prostate cancer (PrCa) risk regions, which provide potential insights into causal mechanisms. Multiple lines of evidence show that a significant proportion of PrCa risk can be explained by germline causal variants that dysregulate nearby target genes in prostate-relevant tissues, thus altering disease risk. The traditional approach to explore this hypothesis has been correlating GWAS variants with steady-state transcript levels, referred to as expression quantitative trait loci (eQTLs). In this work, we assess the utility of chromosome conformation capture (3C) coupled with immunoprecipitation (HiChIP) to identify target genes for PrCa GWAS risk loci. We find that interactome data confirm previously reported PrCa target genes identified through GWAS/eQTL overlap (e.g., MLPH). Interestingly, HiChIP identifies links between PrCa GWAS variants and genes well-known to play a role in prostate cancer biology (e.g., AR) that are not detected by eQTL-based methods. HiChIP predicted enhancer elements at the AR and NKX3-1 prostate cancer risk loci, and both were experimentally confirmed to regulate expression of the corresponding genes through CRISPR interference (CRISPRi) perturbation in LNCaP cells. Our results demonstrate that looping data harbor additional information beyond eQTLs and expand the number of PrCa GWAS loci that can be linked to candidate susceptibility genes.

Predicting master transcription factors from pan-cancer expression data

Critical developmental “master transcription factors” (MTFs) can be subverted during tumorigenesis to control oncogenic transcriptional programs. Current approaches to identifying MTFs rely on ChIP-seq data, which is unavailable for many cancers. We developed the CaCTS (Cancer Core Transcription factor Specificity) algorithm to prioritize candidate MTFs using pan-cancer RNA sequencing data. CaCTS identified candidate MTFs across 34 tumor types and 140 subtypes including predictions for cancer types/subtypes for which MTFs are unknown, including e.g. PAX8, SOX17, and MECOM as candidates in ovarian cancer (OvCa). In OvCa cells, consistent with known MTF properties, these factors are required for viability, lie proximal to superenhancers, co-occupy regulatory elements globally, co-bind loci encoding OvCa biomarkers, and are sensitive to pharmacologic inhibition of transcription. Our predictions of MTFs, especially for tumor types with limited understanding of transcriptional drivers, pave the way to therapeutic targeting of MTFs in a broad spectrum of cancers.

Allele-specific epigenetic activity in prostate cancer and normal prostate tissue implicates prostate cancer risk mechanisms

Genome-wide association studies (GWASs) of prostate cancer have identified >250 significant risk loci, but the causal variants and mechanisms for these loci remain largely unknown. Here, we sought to identify and characterize risk-harboring regulatory elements by integrating epigenomes from primary prostate tumor and normal tissues of 27 individuals across the H3K27ac, H3K4me3, and H3K4me2 histone marks and FOXA1 and HOXB13 transcription factors. We identified 7,371 peaks with significant allele specificity (allele-specific quantitative trait locus [asQTL] peaks). Showcasing their relevance to prostate cancer risk, H3K27ac T-asQTL peaks were the single annotation most enriched for prostate cancer GWAS heritability (40×), significantly higher than corresponding non-asQTL H3K27ac peaks (14×) or coding regions (14×). Surprisingly, fine-mapped GWAS risk variants were most significantly enriched for asQTL peaks observed in tumors, including asQTL peaks that were differentially imbalanced with respect to tumor-normal states. These data pinpointed putative causal regulatory elements at 20 GWAS loci, of which 11 were detected only in the tumor samples. More broadly, tumor-specific asQTLs were enriched for expression QTLs in benign tissues as well as accessible regions found in stem cells, supporting a hypothesis where some germline variants become reactivated during or after transformation and can be captured by epigenomic profiling of the tumor. Our study demonstrates the power of allele specificity in chromatin signals to uncover GWAS mechanisms, highlights the relevance of tumor-specific regulation in the context of cancer risk, and prioritizes multiple loci for experimental follow-up.

Genetic ancestry and clinical outcomes to immune checkpoint inhibitors among seven common cancers

10536

Background: Prior studies and clinical trials report associations between self-reported race and clinical outcomes to Immune Checkpoint Inhibitors (ICIs). However, comprehensive studies of ancestry-associated differences in clinical outcomes have not been performed. We derived genetic ancestry scores and assessed clinical outcomes in 1341 patients with cancer treated with ICIs. Methods: Patients at the Dana-Farber Cancer Institute treated with ICIs only and with relevant cancer types and targeted exome sequencing data (Oncopanel) were included. Relevant cancer types included colorectal adenocarcinoma (CRC), esophagogastric adenocarcinoma (EGC), head and neck squamous cell carcinoma (HNSCC), melanoma, non-small cell lung cancer (NSCLC), renal cell carcinoma (RCC), and urothelial carcinoma (UC). We developed a bioinformatics pipeline to infer fine-scale genetic ancestry for each patient (n=1341) directly from tumor sequencing data by leveraging off and on-target sequenced reads and external ancestry reference panels. Three ancestry scores were determined (African, East Asian, European). Overall survival (OS) and time-to-treatment failure (TTF) were compared by Cox logistic regression between ancestral populations. Hazard ratio (HR) was derived using multivariable analysis, adjusted for single versus combination therapy, prior lines of therapy, and tumor mutational burden (TMB, as percentiles). Results: Median follow-up was 37.8 months (m; interquartile range: 35.7-39.5m). Common cancer types included CRC (n=52), EGC (n=114), HNSCC (n=88), melanoma (n=274), NSCLC (n=571), RCC (n=99), and UC (n=143). A higher East Asian ancestry (EAS) was significantly associated with worse OS ( p=0.03) and TTF ( p=0.002) in patients with RCC, independent of the histologic subtype (Table). There was no significant association between any of the three ancestral populations and clinical outcomes in the other 6 cancer types. Conclusions: We described clinical outcomes to ICIs across three global populations in 7 cancers. As the medical field re-evaluates the use of self-reported race in clinical decision-making, we utilize a novel ancestry pipeline that can be readily applied to tumor-only sequencing panels and better characterize non-white populations. We find no ancestry differences in clinical outcomes except in patients with RCC treated with ICIs which will require future validation. We plan to analyze genomic correlates of response by ancestry in each of the cancer types to better understand these diverge clinical behaviors.[Table: see text]

Impact of Pathogenic Germline DNA Damage Repair alterations on Response to Intense Neoadjuvant Androgen Deprivation Therapy in High-risk Localized Prostate Cancer

BACKGROUND

Intense neoadjuvant androgen deprivation therapy (ADT) before radical prostatectomy (RP) is an investigational approach to reduce recurrence rates in men with high-risk localized prostate cancer (PCa). The impact of germline DNA damage repair (gDDR) gene alterations on response to intense neoadjuvant ADT is not known.

OBJECTIVE

To evaluate the prevalence of gDDR alterations among men with localized PCa at high risk of recurrence and evaluate their impact on response to intense neoadjuvant ADT.

DESIGN, SETTING, AND PARTICIPANTS

We performed germline panel sequencing for 201 men with intermediate- and high-risk localized PCa from five randomized multicenter clinical trials of intense neoadjuvant ADT before RP.

INTERVENTION

Intense neoadjuvant ADT followed by RP.

OUTCOME MEASUREMENTS AND STATISTICAL ANALYSIS

The prevalence of pathogenic gDDR alterations and their association with exceptional pathologic response (complete response or minimal residual disease, defined as residual tumor with the largest cross-section dimension ≤5 mm) to intense neoadjuvant ADT and rates of post-RP biochemical recurrence.

RESULTS AND LIMITATIONS

Pathogenic gDDR alterations were detected in 19 (9.5%) of the 201 PCa patients. The most frequently altered genes were BRCA2 (n = 6; 3.0%) and ATM (n = 4; 2.0%). Patients with gDDR alterations exhibited similar rates of exceptional pathologic response (26% vs 22%), pT3 disease (42% vs 53%), lymph node involvement (5.3% vs 10%), extraprostatic extension (35% vs 54%), and positive margins (5.3% vs 13%) to patients without gDDR alterations (all p > 0.05). The 3-yr biochemical recurrence-free survival was also similar at 45% (95% confidence interval 7.9-78%) for men with gDDR alterations and 55% (95% confidence interval 44-64%) for men without gDDR alterations.

CONCLUSIONS

gDDR alterations are common among men with intermediate- and high-risk localized PCa. Men with gDDR alterations appear to have a comparable response to intense neoadjuvant ADT to that among men without gDDR alterations and should not be excluded from consideration for this treatment approach.

PATIENT SUMMARY

Intense therapy to inhibit the production of androgen hormones (eg, testosterone) before surgery may minimize the risk of cancer recurrence for men with high-risk localized prostate cancer. Inherited mutations in certain DNA repair genes are associated with particularly high rates of recurrence. We found that men with these mutations respond equally well to this intense androgen inhibition before surgery as men without the mutations.

Reprogramming of the FOXA1 cistrome in treatment-emergent neuroendocrine prostate cancer

Lineage plasticity, the ability of a cell to alter its identity, is an increasingly common mechanism of adaptive resistance to targeted therapy in cancer. An archetypal example is the development of neuroendocrine prostate cancer (NEPC) after treatment of prostate adenocarcinoma (PRAD) with inhibitors of androgen signaling. NEPC is an aggressive variant of prostate cancer that aberrantly expresses genes characteristic of neuroendocrine (NE) tissues and no longer depends on androgens. Here, we investigate the epigenomic basis of this resistance mechanism by profiling histone modifications in NEPC and PRAD patient-derived xenografts (PDXs) using chromatin immunoprecipitation and sequencing (ChIP-seq). We identify a vast network of cis-regulatory elements (N~15,000) that are recurrently activated in NEPC. The FOXA1 transcription factor (TF), which pioneers androgen receptor (AR) chromatin binding in the prostate epithelium, is reprogrammed to NE-specific regulatory elements in NEPC. Despite loss of dependence upon AR, NEPC maintains FOXA1 expression and requires FOXA1 for proliferation and expression of NE lineage-defining genes. Ectopic expression of the NE lineage TFs ASCL1 and NKX2-1 in PRAD cells reprograms FOXA1 to bind to NE regulatory elements and induces enhancer activity as evidenced by histone modifications at these sites. Our data establish the importance of FOXA1 in NEPC and provide a principled approach to identifying cancer dependencies through epigenomic profiling.

Epigenetic and transcriptional analysis reveals a core transcriptional program conserved in clonal prostate cancer metastases

The epigenomic regulation of transcriptional programs in metastatic prostate cancer is poorly understood. We studied the epigenomic landscape of prostate cancer drivers using transcriptional profiling and ChIP‐seq in four clonal metastatic tumors derived from a single prostate cancer patient. Our epigenomic analyses focused on androgen receptor (AR), which is a key oncogenic driver in prostate cancer, the AR pioneer factor FOXA1, chromatin insulator CCCTC‐Binding Factor, as well as for modified histones H3K27ac and H3K27me3. The vast majority of AR binding sites were shared among healthy prostate, primary prostate cancer, and metastatic tumor samples, signifying core AR‐driven transcriptional regulation within the prostate cell lineage. Genes associated with core AR‐binding events were significantly enriched for essential genes in prostate cancer cell proliferation. Remarkably, the metastasis‐specific active AR binding sites showed no differential transcriptional output, indicating a robust transcriptional program across metastatic samples. Combined, our data reveal a core transcriptional program in clonal metastatic prostate cancer, despite epigenomic differences in the AR cistrome.

Trans-ethnic variation in germline variants of patients with renal cell carcinoma

Prior studies of the renal cell carcinoma (RCC) germline landscape investigated predominantly patients of European ancestry. We examine the frequency of germline pathogenic and likely pathogenic (P/LP) variants in 1,829 patients with RCC from various ancestries. Overall, P/LP variants are found in 17% of patients, among whom 10.3% harbor one or more clinically actionable variants with potential preventive or therapeutic utility. Patients of African ancestry with RCC harbor significantly more P/LP variants in FH compared to patients of non-African ancestry with RCC and African controls from the Genome Aggregation Database (gnomAD). Patients of non-African ancestry have significantly more P/LP variants in CHEK2 compared to patients of African ancestry with RCC and non-Finnish Europeans controls. Non-Africans with RCC have more actionable variants compared to Africans with RCC. This work helps understand the underlying biological differences in RCC between Africans and non-Africans and paves the way to more comprehensive genomic characterization of underrepresented populations.

Integrative molecular characterization of sarcomatoid and rhabdoid renal cell carcinoma

Sarcomatoid and rhabdoid (S/R) renal cell carcinoma (RCC) are highly aggressive tumors with limited molecular and clinical characterization. Emerging evidence suggests immune checkpoint inhibitors (ICI) are particularly effective for these tumors, although the biological basis for this property is largely unknown. Here, we evaluate multiple clinical trial and real-world cohorts of S/R RCC to characterize their molecular features, clinical outcomes, and immunologic characteristics. We find that S/R RCC tumors harbor distinctive molecular features that may account for their aggressive behavior, including BAP1 mutations, CDKN2A deletions, and increased expression of MYC transcriptional programs. We show that these tumors are highly responsive to ICI and that they exhibit an immune-inflamed phenotype characterized by immune activation, increased cytotoxic immune infiltration, upregulation of antigen presentation machinery genes, and PD-L1 expression. Our findings build on prior work and shed light on the molecular drivers of aggressivity and responsiveness to ICI of S/R RCC.

Sarcomatoid and rhabdoid tumours are highly aggressive forms of renal cell carcinoma that are also responsive to immunotherapy. In this study, the authors perform a comprehensive molecular characterization of these tumours discovering an enrichment of specific alterations and an inflamed phenotype.

CDK4/6 inhibition reprograms the breast cancer enhancer landscape by stimulating AP-1 transcriptional activity

Pharmacologic inhibitors of cyclin-dependent kinases 4 and 6 (CDK4/6) were designed to induce cancer cell cycle arrest. Recent studies have suggested that these agents also exert other effects, influencing cancer cell immunogenicity, apoptotic responses, and differentiation. Using cell-based and mouse models of breast cancer together with clinical specimens, we show that CDK4/6 inhibitors induce remodeling of cancer cell chromatin characterized by widespread enhancer activation, and that this explains many of these effects. The newly activated enhancers include classical super-enhancers that drive luminal differentiation and apoptotic evasion, as well as a set of enhancers overlying endogenous retroviral elements that is enriched for proximity to interferon-driven genes. Mechanistically, CDK4/6 inhibition increases the level of several Activator Protein-1 (AP-1) transcription factor proteins, which are in turn implicated in the activity of many of the new enhancers. Our findings offer insights into CDK4/6 pathway biology and should inform the future development of CDK4/6 inhibitors.

Ovarian Cancer Risk Variants Are Enriched in Histotype-Specific Enhancers and Disrupt Transcription Factor Binding Sites

Quantifying the functional effects of complex disease risk variants can provide insights into mechanisms underlying disease biology. Genome-wide association studies have identified 39 regions associated with risk of epithelial ovarian cancer (EOC). The vast majority of these variants lie in the non-coding genome, where they likely function through interaction with gene regulatory elements. In this study we first estimated the heritability explained by known common low penetrance risk alleles for EOC. The narrow sense heritability (hg2) of EOC overall and high-grade serous ovarian cancer (HGSOCs) were estimated to be 5%-6%. Partitioned SNP heritability across broad functional categories indicated a significant contribution of regulatory elements to EOC heritability. We collated epigenomic profiling data for 77 cell and tissue types from Roadmap Epigenomics and ENCODE, and from H3K27Ac ChIP-seq data generated in 26 ovarian cancer and precursor-related cell and tissue types. We identified significant enrichment of risk single-nucleotide polymorphisms (SNPs) in active regulatory elements marked by H3K27Ac in HGSOCs. To further investigate how risk SNPs in active regulatory elements influence predisposition to ovarian cancer, we used motifbreakR to predict the disruption of transcription factor binding sites. We identified 469 candidate causal risk variants in H3K27Ac peaks that are predicted to significantly break transcription factor (TF) motifs. The most frequently broken motif was REST (p value = 0.0028), which has been reported as both a tumor suppressor and an oncogene. Overall, these systematic functional annotations with epigenomic data improve interpretation of EOC risk variants and shed light on likely cells of origin.

Implementation of Germline Testing for Prostate Cancer: Philadelphia Prostate Cancer Consensus Conference 2019

PURPOSE

Germline testing (GT) is a central feature of prostate cancer (PCA) treatment, management, and hereditary cancer assessment. Critical needs include optimized multigene testing strategies that incorporate evolving genetic data, consistency in GT indications and management, and alternate genetic evaluation models that address the rising demand for genetic services.

METHODS

A multidisciplinary consensus conference that included experts, stakeholders, and national organization leaders was convened in response to current practice challenges and to develop a genetic implementation framework. Evidence review informed questions using the modified Delphi model. The final framework included criteria with strong (> 75%) agreement (Recommend) or moderate (50% to 74%) agreement (Consider).

RESULTS

Large germline panels and somatic testing were recommended for metastatic PCA. Reflex testing-initial testing of priority genes followed by expanded testing-was suggested for multiple scenarios. Metastatic disease or family history suggestive of hereditary PCA was recommended for GT. Additional family history and pathologic criteria garnered moderate consensus. Priority genes to test for metastatic disease treatment included BRCA2, BRCA1, and mismatch repair genes, with broader testing, such as ATM, for clinical trial eligibility. BRCA2 was recommended for active surveillance discussions. Screening starting at age 40 years or 10 years before the youngest PCA diagnosis in a family was recommended for BRCA2 carriers, with consideration in HOXB13, BRCA1, ATM, and mismatch repair carriers. Collaborative (point-of-care) evaluation models between health care and genetic providers was endorsed to address the genetic counseling shortage. The genetic evaluation framework included optimal pretest informed consent, post-test discussion, cascade testing, and technology-based approaches.

CONCLUSION

This multidisciplinary, consensus-driven PCA genetic implementation framework provides novel guidance to clinicians and patients tailored to the precision era. Multiple research, education, and policy needs remain of importance.

Abstract A30: Super-enhancer-associated long noncoding RNA UCA1 interacts directly with AMOT to inhibit Hippo signaling pathway in epithelial ovarian cancer

Long noncoding RNAs (lncRNAs) have emerged as critical regulators of tumorigenesis and yet their mechanistic role remains challenging to characterize. Here, we present LncRNA Interpreter, a strategy integrating functional proteomics with interacting proteins to characterize lncRNAs. Its robustness is exemplified by lncRNA Urothelial Cancer Associated 1 (UCA1), a driver of ovarian cancer development. Reverse phase protein array (RPPA) analysis indicates that UCA1 inhibits Hippo signaling pathway and in vivo RNA antisense purification (iRAP) of UCA1 interacting proteins identified angiomotin (AMOT), a known YAP regulator, as direct binding partner. Loss-of-function experiments show that AMOT mediates Hippo signaling pathway inhibition by UCA1. UCA1 enhances the AMOT-YAP interaction to prevent YAP phosphorylation and facilitate its nuclear translocation. Together, our LncRNA Interpreter pipeline identified UCA1 as an lncRNA regulator of the Hippo signaling pathway and highlighted the UCA1-AMOT-YAP signaling axis in ovarian cancer development. LncRNA Interpreter can readily be applied to other lncRNAs implicated in complex diseases.

Citation Format: Xianzhi Lin, Tassja J. Spindler, Marcos A de Souza Fonseca, Rosario I. Corona, Ji-Heui Seo, Felipe S Dezem, Lewyn Li, Janet M. Lee, Henry W. Long, Thomas A. Sellers, Beth Y. Karlan, Houtan Noushmehr, Matthew L. Freedman, Simon A. Gayther, Kate Lawrenson. Super-enhancer-associated long noncoding RNA UCA1 interacts directly with AMOT to inhibit Hippo signaling pathway in epithelial ovarian cancer [abstract]. In: Proceedings of the AACR Special Conference on the Hippo Pathway: Signaling, Cancer, and Beyond; 2019 May 8-11; San Diego, CA. Philadelphia (PA): AACR; Mol Cancer Res 2020;18(8_Suppl):Abstract nr A30.

A rare variant of African ancestry activates 8q24 lncRNA hub by modulating cancer associated enhancer

Genetic variation at the 8q24 locus is linked with the greater susceptibility to prostate cancer in men of African ancestry. One such African ancestry specific rare variant, rs72725854 (A>G/T) (~6% allele frequency) has been associated with a ~2-fold increase in prostate cancer risk. However, the functional relevance of this variant is unknown. Here we show that the variant rs72725854 is present in a prostate cancer-specific enhancer at 8q24 locus. Chromatin-conformation capture and dCas9 mediated enhancer blocking establish a direct regulatory link between this enhancer and lncRNAs PCAT1, PRNCR1 and PVT1. The risk allele ('T') is associated with higher expression of PCAT1, PVT1 and c-myc in prostate tumors. Further, enhancer with the risk allele gains response to androgen stimulation by recruiting the transcription factor SPDEF whereas, non-risk alleles remain non-responsive. Elevated expression of these lncRNAs and c-myc in risk allele carriers may explain their greater susceptibility to prostate cancer.

Detection of renal cell carcinoma using plasma and urine cell-free DNA methylomes

Improving early cancer detection has the potential to substantially reduce cancer-related mortality. Cell-free methylated DNA immunoprecipitation and high-throughput sequencing (cfMeDIP-seq) is a highly sensitive assay capable of detecting early-stage tumors. We report accurate classification of patients across all stages of renal cell carcinoma (RCC) in plasma (area under the receiver operating characteristic (AUROC) curve of 0.99) and demonstrate the validity of this assay to identify patients with RCC using urine cell-free DNA (cfDNA; AUROC of 0.86).

Detection of urothelial carcinoma using plasma cell-free methylated DNA

5046

Background: Methylation profiling of circulating cell-free DNA (cfDNA) is a promising approach for non-invasive tumor detection due to the presence of tissue-specific epigenetic signatures that are detectable in cfDNA. Cell-free methylated DNA immunoprecipitation and high-throughput sequencing (cfMedDIP-seq) is a sensitive, low-input, cost-effective, bisulfite-free approach to profiling cfDNA methylomes, capable of detecting and classifying various tumor types. We tested the feasibility of cfMeDIP-seq to detect urothelial carcinoma (UC) in plasma samples. Methods: We performed cfMeDIP-seq on plasma samples from 43 patients (pts): 18 metastatic UC (UC) pts, 12 pre-cystectomy non-metastatic UC pts, and 13 cancer-free controls. Six (50%) of pre-cystectomy cases were non-muscle invasive UC. cfDNA was immunoprecipitated and enriched using an antibody targeting 5-methylcytosine and PCR-amplified to create a sequence-ready library. The top differentially methylated regions (DMRs) between UC and control samples were used to train a regularized binomial generalized linear model using 80% of the samples as a training set. The 20% of withheld test samples were then assigned a probability of being UC or control. This process was repeated 100 times. Results: The average amount (standard deviation) of cfDNA isolated from 1 ml of UC plasma samples was 29.2 (27.4) ng/µL and 8.02 (3.58) ng/µL in cancer-free controls. We identified 9,826 DMRs in plasma samples at an adjusted p-value of < 0.01, which partitioned UC and control samples. Iterative training and classification of held out samples using the top 300 DMRs resulted in a mean AUROC of 0.987. Conclusions: cfMeDIP-seq is an interesting new approach for non-invasive detection of UC. cfMeDIP-seq demonstrates high sensitivity to detect UC across all stages of UC, including non-muscle invasive disease.

Non-coding somatic mutations converge on the PAX8 pathway in ovarian cancer

The functional consequences of somatic non-coding mutations in ovarian cancer (OC) are unknown. To identify regulatory elements (RE) and genes perturbed by acquired non-coding variants, here we establish epigenomic and transcriptomic landscapes of primary OCs using H3K27ac ChIP-seq and RNA-seq, and then integrate these with whole genome sequencing data from 232 OCs. We identify 25 frequently mutated regulatory elements, including an enhancer at 6p22.1 which associates with differential expression of ZSCAN16 (P = 6.6 × 10-4) and ZSCAN12 (P = 0.02). CRISPR/Cas9 knockout of this enhancer induces downregulation of both genes. Globally, there is an enrichment of single nucleotide variants in active binding sites for TEAD4 (P = 6 × 10-11) and its binding partner PAX8 (P = 2×10-10), a known lineage-specific transcription factor in OC. In addition, the collection of cis REs associated with PAX8 comprise the most frequently mutated set of enhancers in OC (P = 0.003). These data indicate that non-coding somatic mutations disrupt the PAX8 transcriptional network during OC development.

Mammalian SWI/SNF Complex Genomic Alterations and Immune Checkpoint Blockade in Solid Tumors

Prior data have variably implicated the inactivation of the mammalian SWItch/Sucrose Non-Fermentable (mSWI/SNF) complex with increased tumor sensitivity to immune checkpoint inhibitors (ICI). Herein, we examined the association between mSWI/SNF variants and clinical outcomes to ICIs. We correlated somatic loss-of-function (LOF) variants in a predefined set of mSWI/SNF genes (ARID1A, ARID1B, SMARCA4, SMARCB1, PBRM1, and ARID2) with clinical outcomes in patients with cancer treated with systemic ICIs. We identified 676 patients from Dana-Farber Cancer Institute (DFCI, Boston, MA) and 848 patients from a publicly available database from Memorial Sloan Kettering Cancer Center (MSKCC, New York, NY) who met the inclusion criteria. Multivariable analyses were conducted and adjusted for available baseline factors and tumor mutational burden. Median follow-up was 19.6 (17.6-22.0) months and 28.0 (25.0-29.0) months for the DFCI and MSKCC cohorts, respectively. Seven solid tumor subtypes were examined. In the DFCI cohort, LOF variants of mSWI/SNF did not predict improved overall survival (OS), time-to-treatment failure (TTF), or disease control rate. Only patients with renal cell carcinoma with mSWI/SNF LOF showed significantly improved OS and TTF with adjusted HRs (95% confidence interval) of 0.33 (0.16-0.7) and 0.49 (0.27-0.88), respectively, and this was mostly driven by PRBM1 In the MSKCC cohort, where only OS was captured, LOF mSWI/SNF did not correlate with improved outcomes across any tumor subtype. We did not find a consistent association between mSWI/SNF LOF variants and improved clinical outcomes to ICIs, suggesting that mSWI/SNF variants should not be considered as biomarkers of response to ICIs.

Response to olaparib or carboplatin in a real-world cohort of men with DNA damage repair (DDR) deficient metastatic castration-resistant prostate cancer (mCRPC)

43

Background: PARP inhibitors (PARPi) and platinum chemotherapy cause synthetic lethality and result in clinical benefit in men with mCRPC whose tumors harbor deleterious alterations in DDR genes. BRCA2 is the strongest predictor, but data supports a role for BRCA1, ATM, and other genes. PARPi are more expensive than platinum chemotherapy; whether they are more effective is not known. We evaluated outcomes in men with DDR-deficient mCRPC treated with PARPi or platinum chemotherapy. Methods: We identified men treated at Dana-Farber Cancer Institute who received olaparib or carboplatin for mCRPC and underwent tumor mutation profiling. Progression-free survival (PFS) was determined using Prostate Cancer Working Group 3 guidelines. Pathogenic germline or somatic single-nucleotide variants, insertions/deletions, and copy number loss in BRCA2, BRCA1, and ATM, were identified. Log-rank test was utilized to compare differences in PFS. Results: We identified 27 men who received olaparib: 19/27 (70.4%) received prior taxane chemotherapy and 3/27 (11.1%) prior carboplatin; 16 had mono- or bi-allelic alterations in BRCA2 and 2 had alterations in BRCA1 or ATM. We identified 41 men treated with carboplatin: 39/41 (95.1%) received carboplatin in combination with a taxane; 36/41 (87.8%) received prior taxane chemotherapy and 7/41 (17.1%) prior PARPi; 16 had BRCA2 alterations and 5 had alterations in BRCA1 or ATM. Among men with BRCA2 alterations, there was no significant difference in PFS (HR 0.71, 95% CI 0.45-1.11; p = 0.13) between those treated with olaparib (4.9 months [m]) versus carboplatin (5.4 m). Similarly, no difference in PFS (HR 0.80, 95% CI 0.54-1.16; p = 0.24) was observed among men with BRCA2, BRCA1, or ATM alterations treated with olaparib (3.8 m) versus carboplatin (3.6 m). Conclusions: In a small real-world cohort of men with mCRPC harboring deleterious alterations in BRCA2, BRCA1, and ATM, there was no difference in PFS between those treated with olaparib versus carboplatin. These findings support further investigation into the comparative effectiveness and cost-effectiveness of PARPi and platinum chemotherapy in men with DDR-deficient mCRPC.