Abstract 7007: 5hmc-sequencing identifies candidate genes in Puerto Rican Hispanic/Latino men with aggressive prostate cancer
Authors: Patel Manishkumar S. , Almubarak Muosa , Matta Jaime , Ortiz-Sánchez Carmen , Encarnacion Jarline , Ruiz-Deya Gilberto , Dutil Julie , Dhillon Jasreman , Yamoah Kosj , Berglund Anders , Park Hyun , Wang Liang , Park Jong Y. 
Published: March 2024

Introduction: Puerto Rican (PR) Hispanic/Latino (H/L) men are an understudied population which has the highest prostate cancer (PCa) specific mortality when compared with other Hispanic populations. According to the recent PR Cancer Registry data, PCa is the leading cancer type in terms of incidence (38% of all cancer cases) and mortality (16% of all cancer deaths) in PR H/L men. PR H/L men also have significantly higher PCa-specific mortality than non-Hispanic Whites (NHW) and non-Hispanic Blacks (NHB). So far, little is known about higher mortality in PR H/L men. It is believed that epigenetic changes of key genes play a critical role in aggressive tumors. Therefore, our study was aimed to identify key 5-hydroxymethylcytosine (5hmC) changes in PR H/L men with aggressive PCa.

Methods: We prepared sequencing libraries using the enriched 5hmC DNA from 42 prostate tumor and 44 adjacent normal FFPE samples. We performed sequencing on an Illumina NextSeq 500 with 75bp single end. We used DESeq2 package to identify differentially methylated genes (DMGs). We also downloaded additional 5hmC data (Database: EGAS00001004942) of 51 localized PCa and 7 adjacent normal for confirmation analysis. We used Gleason score (GS) as a cut-off to identify DMGs in aggressive PCa patients (GS- 7 (3+4) or less, Non-aggressive; GS- 7 (4+3) and above, Aggressive).

Results: We identified 808 DMGs in tumors compared to adjacent normal tissues (FDR&lt;0.05, log2FC&gt;|0.4|). Pathway analysis of DMGs demonstrated that cell-cycle, cell division, and DNA repair related pathways were most upregulated in tumors. Since 5hmC abundance is highly correlated with gene expression levels, we further investigated 808 DMGs with TCGA prostate cancer gene expression data. This analysis identified 59/73 DMGs (80.1%, (FDR&lt;0.05, delta&gt;|1|) with significant 5hmC and gene expression changes in the same direction. We also performed differential methylation analysis in another dataset and observed that 129 DMGs were common with our results. Finally, we demonstrated 111 DMGs (p&lt;0.05) in aggressive tumors (n=7) compared to non-aggressive tumors (n=15). These 111 DMGs include important 5hmC hypomethylated genes (BCCIP and KLK10) and hypermethylated genes (COX6C, ARMC2, PVT1) in aggressive PCa patients.

Conclusions: Our study identifies 59 genes having coordinated epigenetic and transcriptomic changes in PR H/L men and 111 DMGs in aggressive tumors. The coordinated 59 genes include tumor suppressor genes such as DKK3 and PRDM8 with downregulated 5hmC and gene expression levels. Also, previous studies reported downregulation of BCCIP and KLK10 and upregulation of COX6C, ARMC2, PVT1 expression in PCa. The coordinated changes in 5hmC could be critical in aggressive PCa biology, hence in future, these genes will be validated as candidate biomarkers and targets for aggressive PCa among PR H/L men.

Citation Format: Manishkumar S. Patel, Muosa Almubarak, Jaime Matta, Carmen Ortiz-Sánchez, Jarline Encarnacion, Gilberto Ruiz-Deya, Julie Dutil, Jasreman Dhillon, Kosj Yamoah, Anders Berglund, Hyun Park, Liang Wang, Jong Y. Park. 5hmc-sequencing identifies candidate genes in Puerto Rican Hispanic/Latino men with aggressive prostate cancer [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2024; Part 1 (Regular Abstracts); 2024 Apr 5-10; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2024;84(6_Suppl):Abstract nr 7007.