Using gene and gene-set association tests to identify lethal prostate cancer genes.

BACKGROUND: Recent advances in the detection and treatment of prostate cancer (PCa) have reduced morbidity and mortality from this common cancer. Despite these improvements, PCa remains the second leading cause of cancer death in men in the United States. Further understanding of the genetic underpinnings of lethal PCa is required to drive risk detection and prevention and ultimately reduce mortality. We therefore set out to identify germline variants associated with cases of lethal prostate cancer (LPCa). METHODS: Using a two-stage study design, we compared whole-exome sequencing data of 550 LPCa patients to 488 healthy male controls. Men were classified as having LPCa based on medical record review. Candidate genes were identified using gene- and gene-set-based rare truncating variant association tests. Case-control burden testing through Firth's penalized logistic regression and case-gnomAD allelic burden testing through a one-sided mid-p Fisher's exact test were conducted. Each gene's p-values from these tests were combined into an omnibus p-value for candidate gene selection. In the subsequent validation stage, genes were assessed using the UK Biobank and Firth's penalized logistic regression for each ancestry, combined through meta-analysis. RESULTS: Gene-based rare variant association tests identified 12 genes nominally associated with LPCa. Rare-variant association tests identified a gene set with a significantly higher burden of truncating germline mutations in LPCa patients than controls. Combining gene- and gene-set test results, four nominally significant genes (PPP1R3A, TG, PPFIBP2, and BTN3A3) were selected as candidates. Subsequent validation using the UK Biobank found that PPP1R3A was significantly associated with LPCa risk (odds ratio 2.34, CI 1.20-4.59). Specifically, pGln662ArgfsTer7 was identified as the predominant variant in PPP1R3A among LPCa patients in our dataset. CONCLUSIONS: Both individual gene and gene-set analyses identified candidates associated with LPCa. The novel association of PPP1R3A and LPCa risk merits further investigation.

Outcomes of 10 years of PSA screening for prostate cancer in Norwegian men with Lynch syndrome.

BACKGROUND: Pathogenic germline variants in the mismatch repair (MMR) genes are associated with an increased risk of prostate cancer (PCa). Since 2010 we have recommended MMR carriers annual PSA testing from the age of 40. Prospective studies of the outcome of long-term PSA screening are lacking. This study aimed to investigate the incidence and characteristics of PCa in Norwegian MMR carriers attending annual PSA screening (PSA threshold >3.0 ng/mL) to evaluate whether our recommendations should be continued. METHODS: This is a prospective observational study of 225 male MMR carriers who were recommended annual PSA screening by the Section of Inherited Cancer, Oslo University Hospital from 2010 and onwards. Incidence and tumor characteristics (age, PSA at diagnosis, Gleason score, TNM score) were described. IHC and MSI-analyses were done on available tumors. Standardized incidence ratio (SIR) was calculated based on data from the Cancer Registry of Norway. RESULTS: Twenty-two of 225 (9.8%) had been diagnosed with PCa, including 10/69 (14.5%) MSH2 carriers and 8/61 (13.1%) MSH6 carriers. Ten of 20 (50%) tumors had Gleason score ≥4 + 3 on biopsy and 6/11 (54.5%) had a pathological T3a/b stage. Eight of 17 (47.1%) tumors showed abnormal staining on IHC and 3/13 (23.1%) were MSI-high. SIR was 9.54 (95% CI 5.98-14.45) for all MMR genes, 13.0 (95% CI 6.23-23.9) for MSH2 and 13.74 for MSH6 (95% CI 5.93-27.08). CONCLUSIONS: Our results indicate that the MMR genes, and especially MSH2 and MSH6, are associated with a significant risk of PCa, and a high number of tumors show aggressive characteristics. While the impact of screening on patient outcomes remains to be more firmly established, the high SIR values we observe provide support for continued PSA screening of MSH2 and MSH6 carriers. Studies are needed to provide optimal recommendations for PSA-threshold and to evaluate whether MLH1 and PMS2 carriers should not be recommended screening.

De novo variants identified by trio whole exome sequencing of bladder exstrophy epispadias complex.

Bladder exstrophy epispadias complex (BEEC) encompasses a spectrum of conditions ranging from mild epispadias to the most severe form: omphalocele-bladder exstrophy-imperforate anus-spinal defects (OEIS). BEEC involves abnormalities related to anatomical structures that are proposed to have a similar underlying etiology and pathogenesis. In general, BEEC, is considered to arise from a sequence of events in embryonic development and is believed to be a multi-etiological disease with contributions from genetic and environmental factors. Several genes have been implicated and mouse models have been generated, including a knockout model of p63, which is involved in the synthesis of stratified epithelium. Mice lacking p63 have undifferentiated ventral urothelium. MNX1 has also been implicated. In addition, cigarette smoking, diazepam and clomid have been implied as environmental factors due to their relative association. By in large, the etiology and pathogenesis of human BEEC is unknown. We performed de novo analysis of whole exome sequencing (WES) of germline samples from 31 unrelated trios where the probands have a diagnosis of BEEC syndrome. We also evaluated the DECIPHER database to identify copy number variants (CNVs) in genes in individuals with the search terms "bladder exstrophy" in an attempt to identify additional candidate genes within these regions. Several de novo variants were identified; however, a candidate gene is still unclear. This data further supports the multi-etiological nature of BEEC.

FASN Gene Methylation is Associated with Fatty Acid Synthase Expression and Clinical-genomic Features of Prostate Cancer.

Fatty acid synthase (FASN) catalyzes the synthesis of long-chain saturated fatty acids and is overexpressed during prostatic tumorigenesis, where it is the therapeutic target in several ongoing trials. However, the mechanism of FASN upregulation in prostate cancer remains unclear. Here, we examine FASN gene CpG methylation pattern by InfiniumEPIC profiling and whole-genome bisulfite sequencing across multiple racially diverse primary and metastatic prostate cancer cohorts, comparing with FASN protein expression as measured by digitally quantified IHC assay and reverse phase protein array analysis or FASN gene expression. We demonstrate that the FASN gene body is hypomethylated and overexpressed in primary prostate tumors compared with benign tissue, and FASN gene methylation is significantly inversely correlated with FASN protein or gene expression in both primary and metastatic prostate cancer. Primary prostate tumors with ERG gene rearrangement have increased FASN expression and we find evidence of FASN hypomethylation in this context. FASN expression is also significantly increased in prostate tumors from carriers of the germline HOXB13 G84E mutation compared with matched controls, consistent with a report that HOXB13 may contribute to epigenetic regulation of FASN in vitro. However, in contrast to previous studies, we find no significant association of FASN expression or methylation with self-identified race in models that include ERG status across two independent primary tumor cohorts. Taken together, these data support a potential epigenetic mechanism for FASN regulation in the prostate which may be relevant for selecting patients responsive to FASN inhibitors. SIGNIFICANCE: Here, we leverage multiple independent primary and metastatic prostate cancer cohorts to demonstrate that FASN gene body methylation is highly inversely correlated with FASN gene and protein expression. This finding may shed light on epigenetic mechanisms of FASN regulation in prostate cancer and provides a potentially useful biomarker for selecting patients in future trials of FASN inhibitors.

Characterization of HOXB13 expression patterns in localized and metastatic castration-resistant prostate cancer.

HOXB13 is a key lineage homeobox transcription factor that plays a critical role in the differentiation of the prostate gland. Several studies have suggested that HOXB13 alterations may be involved in prostate cancer development and progression. Despite its potential biological relevance, little is known about the expression of HOXB13 across the disease spectrum of prostate cancer. To this end, we validated a HOXB13 antibody using genetic controls and investigated HOXB13 protein expression in murine and human developing prostates, localized prostate cancers, and metastatic castration-resistant prostate cancers. We observed that HOXB13 expression increases during later stages of murine prostate development. All localized prostate cancers showed HOXB13 protein expression. Interestingly, lower HOXB13 expression levels were observed in higher-grade tumors, although no significant association between HOXB13 expression and recurrence or disease-specific survival was found. In advanced metastatic prostate cancers, HOXB13 expression was retained in the majority of tumors. While we observed lower levels of HOXB13 protein and mRNA levels in tumors with evidence of lineage plasticity, 84% of androgen receptor-negative castration-resistant prostate cancers and neuroendocrine prostate cancers (NEPCs) retained detectable levels of HOXB13. Notably, the reduced expression observed in NEPCs was associated with a gain of HOXB13 gene body CpG methylation. In comparison to the commonly used prostate lineage marker NKX3.1, HOXB13 showed greater sensitivity in detecting advanced metastatic prostate cancers. Additionally, in a cohort of 837 patients, 383 with prostatic and 454 with non-prostatic tumors, we found that HOXB13 immunohistochemistry had a 97% sensitivity and 99% specificity for prostatic origin. Taken together, our studies provide valuable insight into the expression pattern of HOXB13 during prostate development and cancer progression. Furthermore, our findings support the utility of HOXB13 as a diagnostic biomarker for prostate cancer, particularly to confirm the prostatic origin of advanced metastatic castration-resistant tumors. © 2023 The Pathological Society of Great Britain and Ireland.

Clinical and Functional Analyses of an African-ancestry Gain-of-function HOXB13 Variant Implicated in Aggressive Prostate Cancer.

BACKGROUND: Recent reports have uncovered a HOXB13 variant (X285K) predisposing to prostate cancer in men of West African ancestry. The clinical relevance and protein function associated with this inherited variant are unknown. OBJECTIVE: To determine the clinical relevance of HOXB13 (X285K) in comparison with HOXB13 (G84E) and BRCA2 pathogenic/likely pathogenic (P/LP) variants, and to elucidate the oncogenic mechanisms of the X285K protein. DESIGN, SETTING, AND PARTICIPANTS: Real-world data were collected from 21,393 men with prostate cancer undergoing genetic testing from 2019 to 2022, and in vitro cell-line models were established for the evaluation of oncogenic functions associated with the X285K protein. OUTCOME MEASUREMENTS AND STATISTICAL ANALYSIS: Genetic testing results were compared among patient groups according to self-reported race/ethnicity, Gleason scores, and American Joint Committee on Cancer stages using the exact test. Oncogenic functions of X285K were evaluated by RNA sequencing, chromatin immunoprecipitation sequencing, and Western blot analyses. RESULTS AND LIMITATIONS: HOXB13 (X285K) was significantly enriched in self-reported Black (1.01%) versus White (0.01%) patients. We observed a trend of more aggressive disease in the HOXB13 (X285K) and BRCA2 P/LP carriers than in the HOXB13 (G84E) carriers. Replacement of the wild-type HOXB13 protein with the X285K protein resulted in a gain of an E2F/MYC signature, validated by the elevated expression of cyclin B1 and c-Myc, without affecting the androgen response signature. Elevated expression of cyclin B1 and c-Myc was explained by enhanced binding of the X285K protein to the promoters and enhancers of these genes. The limitations of the study are the lack of complete clinical outcome data for all patients studied and the use of a single cell line in the functional analysis. CONCLUSIONS: HOXB13 (X285K) is significantly enriched in self-reported Black patients, and X285K carriers detected in the real-world clinical setting have aggressive prostate cancer features similar to the BRCA2 carriers. Functional studies revealed a unique gain-of-function oncogenic mechanism of X285K protein in regulating E2F/MYC signatures. PATIENT SUMMARY: The HOXB13 (X285K) variant is clinically and functionally linked to aggressive prostate cancer, supporting genetic testing for X285K in Black men and early disease screening of carriers of this variant.

Identification of long noncoding RNAs with aberrant expression in prostate cancer metastases.

Prostate cancer (PCa) is the second-most common cause of male cancer-related death in western industrialized countries, and the emergence of metastases is a key challenge in the treatment of PCa. Accumulating studies have shown that long noncoding RNAs (lncRNAs) play an important role in the regulation of diverse cellular and molecular processes during the development and progression of cancer. Here, we utilized a unique cohort of castration-resistant prostate cancer metastases (mCRPC) and corresponding localized tumors and RNA sequencing (RNA-seq). First, we showed that patient-to-patient variability accounted for most of the variance in lncRNA expression between the samples, suggesting that genomic alterations in the samples are the main drivers of lncRNA expression in PCa metastasis. Subsequently, we identified 27 lncRNAs with differential expression (DE-lncRNAs) between metastases and corresponding primary tumors, suggesting that they are mCRPC-specific lncRNAs. Analyses of potential regulation by transcription factors (TFs) revealed that approximately half of the DE-lncRNAs have at least one binding site for the androgen receptor in their regulatory regions. In addition, TF enrichment analysis revealed the enrichment of binding sites for PCa-associated TFs, such as FOXA1 and HOXB13, in the regulatory regions of the DE-lncRNAs. In a cohort of prostatectomy-treated prostate tumors, four of the DE-lncRNAs showed association with progression-free time and two of them (lnc-SCFD2-2 and lnc-R3HCC1L-8) were independent prognostic markers. Our study highlights several mCRPC-specific lncRNAs that might be important in the progression of the disease to the metastatic stage and may also serve as potential biomarkers for aggressive PCa.

Evidence of Novel Susceptibility Variants for Prostate Cancer and a Multiancestry Polygenic Risk Score Associated with Aggressive Disease in Men of African Ancestry.

BACKGROUND: Genetic factors play an important role in prostate cancer (PCa) susceptibility. OBJECTIVE: To discover common genetic variants contributing to the risk of PCa in men of African ancestry. DESIGN, SETTING, AND PARTICIPANTS: We conducted a meta-analysis of ten genome-wide association studies consisting of 19378 cases and 61620 controls of African ancestry. OUTCOME MEASUREMENTS AND STATISTICAL ANALYSIS: Common genotyped and imputed variants were tested for their association with PCa risk. Novel susceptibility loci were identified and incorporated into a multiancestry polygenic risk score (PRS). The PRS was evaluated for associations with PCa risk and disease aggressiveness. RESULTS AND LIMITATIONS: Nine novel susceptibility loci for PCa were identified, of which seven were only found or substantially more common in men of African ancestry, including an African-specific stop-gain variant in the prostate-specific gene anoctamin 7 (ANO7). A multiancestry PRS of 278 risk variants conferred strong associations with PCa risk in African ancestry studies (odds ratios [ORs] >3 and >5 for men in the top PRS decile and percentile, respectively). More importantly, compared with men in the 40-60% PRS category, men in the top PRS decile had a significantly higher risk of aggressive PCa (OR = 1.23, 95% confidence interval = 1.10-1.38, p = 4.4 × 10-4). CONCLUSIONS: This study demonstrates the importance of large-scale genetic studies in men of African ancestry for a better understanding of PCa susceptibility in this high-risk population and suggests a potential clinical utility of PRS in differentiating between the risks of developing aggressive and nonaggressive disease in men of African ancestry. PATIENT SUMMARY: In this large genetic study in men of African ancestry, we discovered nine novel prostate cancer (PCa) risk variants. We also showed that a multiancestry polygenic risk score was effective in stratifying PCa risk, and was able to differentiate risk of aggressive and nonaggressive disease.

Association of self-identified race and genetic ancestry with the immunogenomic landscape of primary prostate cancer.

The genomic and immune landscapes of prostate cancer differ by self-identified race. However, few studies have examined the genome-wide copy number landscape and immune content of matched cohorts with genetic ancestry data and clinical outcomes. Here, we assessed prostate cancer somatic copy number alterations (sCNA) and tumor immune content of a grade-matched, surgically treated cohort of 145 self-identified Black (BL) and 145 self-identified White (WH) patients with genetic ancestry estimation. A generalized linear model adjusted with age, preoperative prostate-specific antigen (PSA), and Gleason Grade Group and filtered for germline copy number variations (gCNV) identified 143 loci where copy number varied significantly by percent African ancestry, clustering on chromosomes 6p, 10q, 11p, 12p, and 17p. Multivariable Cox regression models adjusted for age, preoperative PSA levels, and Gleason Grade Group revealed that chromosome 8q gains (including MYC) were significantly associated with biochemical recurrence and metastasis, independent of genetic ancestry. Finally, Treg density in BL and WH patients was significantly correlated with percent genome altered, and these findings were validated in the TCGA cohort. Taken together, our findings identify specific sCNA linked to genetic ancestry and outcome in primary prostate cancer and demonstrate that Treg infiltration varies by global sCNA burden in primary disease.