Spectrum of somatic mutational features of colorectal tumors in ancestrally diverse populations.

Ancestrally diverse and admixed populations, including the Hispanic/Latino/a/x/e community, are underrepresented in cancer genetic and genomic studies. Leveraging the Latino Colorectal Cancer Consortium, we analyzed whole exome sequencing data on tumor/normal pairs from 718 individuals with colorectal cancer (128 Latino, 469 non-Latino) to map somatic mutational features by ethnicity and genetic ancestry. Global proportions of African, East Asian, European, and Native American ancestries were estimated using ADMIXTURE. Associations between global genetic ancestry and somatic mutational features across genes were examined using logistic regression. TP53 , APC , and KRAS were the most recurrently mutated genes. Compared to non-Latino individuals, tumors from Latino individuals had fewer KRAS (OR=0.64, 95%CI=0.41-0.97, p=0.037) and PIK3CA mutations (OR=0.55, 95%CI=0.31-0.98, p=0.043). Genetic ancestry was associated with presence of somatic mutations in 39 genes (FDR-adjusted LRT p<0.05). Among these genes, a 10% increase in African ancestry was associated with significantly higher odds of mutation in KNCN (OR=1.34, 95%CI=1.09-1.66, p=5.74×10 -3 ) and TMEM184B (OR=1.53, 95%CI=1.10-2.12, p=0.011). Among RMGs, we found evidence of association between genetic ancestry and mutation status in CDC27 (LRT p=0.0084) and between SMAD2 mutation status and AFR ancestry (OR=1.14, 95%CI=1.00-1.30, p=0.046). Ancestry was not associated with tumor mutational burden. Individuals with above-average Native American ancestry had a lower frequency of microsatellite instable (MSI-H) vs microsatellite stable tumors (OR=0.45, 95%CI=0.21-0.99, p=0.048). Our findings provide new knowledge about the relationship between ancestral haplotypes and somatic mutational profiles that may be useful in developing precision medicine approaches and provide additional insight into genomic contributions to cancer disparities. Significance: Our data in ancestrally diverse populations adds essential information to characterize mutational features in the colorectal cancer genome. These results will help enhance equity in the development of precision medicine strategies.

Targeting CaMKK2 Inhibits Actin Cytoskeletal Assembly to Suppress Cancer Metastasis.

Triple-negative breast cancers (TNBC) tend to become invasive and metastatic at early stages in their development. Despite some treatment successes in early-stage localized TNBC, the rate of distant recurrence remains high, and long-term survival outcomes remain poor. In a search for new therapeutic targets for this disease, we observed that elevated expression of the serine/threonine kinase calcium/calmodulin (CaM)-dependent protein kinase kinase 2 (CaMKK2) is highly correlated with tumor invasiveness. In validation studies, genetic disruption of CaMKK2 expression or inhibition of its activity with small molecule inhibitors disrupted spontaneous metastatic outgrowth from primary tumors in murine xenograft models of TNBC. High-grade serous ovarian cancer (HGSOC), a high-risk, poor prognosis ovarian cancer subtype, shares many features with TNBC, and CaMKK2 inhibition effectively blocked metastatic progression in a validated xenograft model of this disease. Mechanistically, CaMKK2 increased the expression of the phosphodiesterase PDE1A, which hydrolyzed cyclic guanosine monophosphate (cGMP) to decrease the cGMP-dependent activity of protein kinase G1 (PKG1). Inhibition of PKG1 resulted in decreased phosphorylation of vasodilator-stimulated phosphoprotein (VASP), which in its hypophosphorylated state binds to and regulates F-actin assembly to facilitate cell movement. Together, these findings establish a targetable CaMKK2-PDE1A-PKG1-VASP signaling pathway that controls cancer cell motility and metastasis by impacting the actin cytoskeleton. Furthermore, it identifies CaMKK2 as a potential therapeutic target that can be exploited to restrict tumor invasiveness in patients diagnosed with early-stage TNBC or localized HGSOC. SIGNIFICANCE: CaMKK2 regulates actin cytoskeletal dynamics to promote tumor invasiveness and can be inhibited to suppress metastasis of breast and ovarian cancer, indicating CaMKK2 inhibition as a therapeutic strategy to arrest disease progression.

Ca 2+ /Calmodulin Dependent Protein Kinase Kinase-2 (CaMKK2) promotes Protein Kinase G (PKG)-dependent actin cytoskeletal assembly to increase tumor metastasis.

Triple-negative breast cancers (TNBCs) tend to become highly invasive early during cancer development. Despite some successes in the initial treatment of patients diagnosed with early-stage localized TNBC, the rate of metastatic recurrence remains high with poor long-term survival outcomes. Here we show that elevated expression of the serine/threonine-kinase, Calcium/Calmodulin (CaM)-dependent protein kinase kinase-2 (CaMKK2), is highly correlated with tumor invasiveness. We determined that genetic disruption of CaMKK2 expression, or inhibition of its activity, disrupted spontaneous metastatic outgrowth from primary tumors in murine xenograft models of TNBC. High-grade serous ovarian cancer (HGSOC), a high-risk, poor-prognosis ovarian cancer subtype, shares many genetic features with TNBC, and importantly, CaMKK2 inhibition effectively blocked metastatic progression in a validated xenograft model of this disease. Probing the mechanistic links between CaMKK2 and metastasis we defined the elements of a new signaling pathway that impacts actin cytoskeletal dynamics in a manner which increases cell migration/invasion and metastasis. Notably, CaMKK2 increases the expression of the phosphodiesterase PDE1A which decreases the cGMP-dependent activity of protein kinase G1 (PKG1). This inhibition of PKG1 results in decreased phosphorylation of Vasodilator-Stimulated Phosphoprotein (VASP), which in its hypophosphorylated state binds to and regulates F-actin assembly to facilitate contraction/cell movement. Together, these data establish a targetable CaMKK2-PDE1A-PKG1-VASP signaling pathway that controls cancer cell motility and metastasis. Further, it credentials CaMKK2 as a therapeutic target that can be exploited in the discovery of agents for use in the neoadjuvant/adjuvant setting to restrict tumor invasiveness in patients diagnosed with early-stage TNBC or localized HGSOC.

A widespread length-dependent splicing dysregulation in cancer.

Dysregulation of alternative splicing is a key molecular hallmark of cancer. However, the common features and underlying mechanisms remain unclear. Here, we report an intriguing length-dependent splicing regulation in cancers. By systematically analyzing the transcriptome of thousands of cancer patients, we found that short exons are more likely to be mis-spliced and preferentially excluded in cancers. Compared to other exons, cancer-associated short exons (CASEs) are more conserved and likely to encode in-frame low-complexity peptides, with functional enrichment in GTPase regulators and cell adhesion. We developed a CASE-based panel as reliable cancer stratification markers and strong predictors for survival, which is clinically useful because the detection of short exon splicing is practical. Mechanistically, mis-splicing of CASEs is regulated by elevated transcription and alteration of certain RNA binding proteins in cancers. Our findings uncover a common feature of cancer-specific splicing dysregulation with important clinical implications in cancer diagnosis and therapies.

Adaptive stress response genes associated with breast cancer subtypes and survival outcomes reveal race-related differences.

Aggressive breast cancer variants, like triple negative and inflammatory breast cancer, contribute to disparities in survival and clinical outcomes among African American (AA) patients compared to White (W) patients. We previously identified the dominant role of anti-apoptotic protein XIAP in regulating tumor cell adaptive stress response (ASR) that promotes a hyperproliferative, drug resistant phenotype. Using The Cancer Genome Atlas (TCGA), we identified 46-88 ASR genes that are differentially expressed (2-fold-change and adjusted p-value < 0.05) depending on PAM50 breast cancer subtype. On average, 20% of all 226 ASR genes exhibited race-related differential expression. These genes were functionally relevant in cell cycle, DNA damage response, signal transduction, and regulation of cell death-related processes. Moreover, 23% of the differentially expressed ASR genes were associated with AA and/or W breast cancer patient survival. These identified genes represent potential therapeutic targets to improve breast cancer outcomes and mitigate associated health disparities.

Deciphering associations between three RNA splicing-related genetic variants and lung cancer risk.

Limited efforts have been made in assessing the effect of genome-wide profiling of RNA splicing-related variation on lung cancer risk. In the present study, we first identified RNA splicing-related genetic variants linked to lung cancer in a genome-wide profiling analysis and then conducted a two-stage (discovery and replication) association study in populations of European ancestry. Discovery and validation were conducted sequentially with a total of 29,266 cases and 56,450 controls from both the Transdisciplinary Research in Cancer of the Lung and the International Lung Cancer Consortium as well as the OncoArray database. For those variants identified as significant in the two datasets, we further performed stratified analyses by smoking status and histological type and investigated their effects on gene expression and potential regulatory mechanisms. We identified three genetic variants significantly associated with lung cancer risk: rs329118 in JADE2 (P = 8.80E-09), rs2285521 in GGA2 (P = 4.43E-08), and rs198459 in MYRF (P = 1.60E-06). The combined effects of all three SNPs were more evident in lung squamous cell carcinomas (P = 1.81E-08, P = 6.21E-08, and P = 7.93E-04, respectively) than in lung adenocarcinomas and in ever smokers (P = 9.80E-05, P = 2.70E-04, and P = 2.90E-05, respectively) than in never smokers. Gene expression quantitative trait analysis suggested a role for the SNPs in regulating transcriptional expression of the corresponding target genes. In conclusion, we report that three RNA splicing-related genetic variants contribute to lung cancer susceptibility in European populations. However, additional validation is needed, and specific splicing mechanisms of the target genes underlying the observed associations also warrants further exploration.

Characterization of a castrate-resistant prostate cancer xenograft derived from a patient of West African ancestry.

BACKGROUND: Prostate cancer is a clinically and molecularly heterogeneous disease, with highest incidence and mortality among men of African ancestry. To date, prostate cancer patient-derived xenograft (PCPDX) models to study this disease have been difficult to establish because of limited specimen availability and poor uptake rates in immunodeficient mice. Ancestrally diverse PCPDXs are even more rare, and only six PCPDXs from self-identified African American patients from one institution were recently made available. METHODS: In the present study, we established a PCPDX from prostate cancer tissue from a patient of estimated 90% West African ancestry with metastatic castration resistant disease, and characterized this model's pathology, karyotype, hotspot mutations, copy number, gene fusions, gene expression, growth rate in normal and castrated mice, therapeutic response, and experimental metastasis. RESULTS: This PCPDX has a mutation in TP53 and loss of PTEN and RB1. We have documented a 100% take rate in mice after thawing the PCPDX tumor from frozen stock. The PCPDX is castrate- and docetaxel-resistant and cisplatin-sensitive, and has gene expression patterns associated with such drug responses. After tail vein injection, the PCPDX tumor cells can colonize the lungs of mice. CONCLUSION: This PCPDX, along with others that are established and characterized, will be useful pre-clinically for studying the heterogeneity of prostate cancer biology and testing new therapeutics in models expected to be reflective of the clinical setting.

A Zebrafish Model of Metastatic Colonization Pinpoints Cellular Mechanisms of Circulating Tumor Cell Extravasation.

Metastasis is a multistep process in which cells must detach, migrate/invade local structures, intravasate, circulate, extravasate, and colonize. A full understanding of the complexity of this process has been limited by the lack of ability to study these steps in isolation with detailed molecular analyses. Leveraging a comparative oncology approach, we injected canine osteosarcoma cells into the circulation of transgenic zebrafish with fluorescent blood vessels in a biologically dynamic metastasis extravasation model. Circulating tumor cell clusters that successfully extravasated the vasculature as multicellular units were isolated under intravital imaging (n = 6). These extravasation-positive tumor cell clusters sublines were then molecularly profiled by RNA-Seq. Using a systems-level analysis, we pinpointed the downregulation of KRAS signaling, immune pathways, and extracellular matrix (ECM) organization as enriched in extravasated cells (p < 0.05). Within the extracellular matrix remodeling pathway, we identified versican (VCAN) as consistently upregulated and central to the ECM gene regulatory network (p < 0.05). Versican expression is prognostic for a poorer metastasis-free and overall survival in patients with osteosarcoma. Together, our results provide a novel experimental framework to study discrete steps in the metastatic process. Using this system, we identify the versican/ECM network dysregulation as a potential contributor to osteosarcoma circulating tumor cell metastasis.

A prospective trial of abiraterone acetate plus prednisone in Black and White men with metastatic castrate-resistant prostate cancer.

BACKGROUND: Retrospective analyses of randomized trials suggest that Black men with metastatic castration-resistant prostate cancer (mCRPC) have longer survival than White men. The authors conducted a prospective study of abiraterone acetate plus prednisone to explore outcomes by race. METHODS: This race-stratified, multicenter study estimated radiographic progression-free survival (rPFS) in Black and White men with mCRPC. Secondary end points included prostate-specific antigen (PSA) kinetics, overall survival (OS), and safety. Exploratory analysis included genome-wide genotyping to identify single nucleotide polymorphisms associated with progression in a model incorporating genetic ancestry. One hundred patients self-identified as White (n = 50) or Black (n = 50) were enrolled. Eligibility criteria were modified to facilitate the enrollment of individual Black patients. RESULTS: The median rPFS for Black and White patients was 16.6 and 16.8 months, respectively; their times to PSA progression (TTP) were 16.6 and 11.5 months, respectively; and their OS was 35.9 and 35.7 months, respectively. Estimated rates of PSA decline by ≥50% in Black and White patients were 74% and 66%, respectively; and PSA declines to <0.2 ng/mL were 26% and 10%, respectively. Rates of grade 3 and 4 hypertension, hypokalemia, and hyperglycemia were higher in Black men. CONCLUSIONS: Multicenter prospective studies by race are feasible in men with mCRPC but require less restrictive eligibility. Despite higher comorbidity rates, Black patients demonstrated rPFS and OS similar to those of White patients and trended toward greater TTP and PSA declines, consistent with retrospective reports. Importantly, Black men may have higher side-effect rates than White men. This exploratory genome-wide analysis of TTP identified a possible candidate marker of ancestry-dependent treatment outcomes.

Dysregulation of the ESRP2-NF2-YAP/TAZ axis promotes hepatobiliary carcinogenesis in non-alcoholic fatty liver disease.

BACKGROUND & AIMS: Non-alcoholic fatty liver disease (NAFLD), the hepatic correlate of the metabolic syndrome, is a major risk factor for hepatobiliary cancer (HBC). Although chronic inflammation is thought to be the root cause of all these diseases, the mechanism whereby it promotes HBC in NAFLD remains poorly understood. Herein, we aim to evaluate the hypothesis that inflammation-related dysregulation of the ESRP2-NF2-YAP/TAZ axis promotes HB carcinogenesis. METHODS: We use murine NAFLD models, liver biopsies from patients with NAFLD, human liver cancer registry data, and studies in liver cancer cell lines. RESULTS: Our results confirm the hypothesis that inflammation-related dysregulation of the ESRP2-NF2-YAP/TAZ axis promotes HB carcinogenesis, supporting a model whereby chronic inflammation suppresses hepatocyte expression of ESRP2, an RNA splicing factor that directly targets and activates NF2, a tumor suppressor that is necessary to constrain YAP/TAZ activation. The resultant loss of NF2 function permits sustained YAP/TAZ activity that drives hepatocyte proliferation and de-differentiation. CONCLUSION: Herein, we report on a novel mechanism by which chronic inflammation leads to sustained activation of YAP/TAZ activity; this imposes a selection pressure that favors liver cells with mutations enabling survival during chronic oncogenic stress. LAY SUMMARY: Non-alcoholic fatty liver disease (NAFLD) increases the risk of hepatobiliary carcinogenesis. However, the underlying mechanism remains unknown. Our study demonstrates that chronic inflammation suppresses hepatocyte expression of ESRP2, an adult RNA splicing factor that activates NF2. Thus, inactive (fetal) NF2 loses the ability to activate Hippo kinases, leading to the increased activity of downstream YAP/TAZ and promoting hepatobiliary carcinogenesis in chronically injured livers.

Differential alternative RNA splicing and transcription events between tumors from African American and White patients in The Cancer Genome Atlas.

Individuals of African ancestry suffer disproportionally from higher incidence, aggressiveness, and mortality for particular cancers. This disparity likely results from an interplay among differences in multiple determinants of health, including differences in tumor biology. We used The Cancer Genome Atlas (TCGA) SpliceSeq and TCGA aggregate expression datasets and identified differential alternative RNA splicing and transcription events (ARS/T) in cancers between self-identified African American (AA) and White (W) patients. We found that retained intron events were enriched among race-related ARS/T. In addition, on average, 12% of the most highly ranked race-related ARS/T overlapped between any two analyzed cancers. Moreover, the genes undergoing race-related ARS/T functioned in cancer-promoting pathways, and a number of race-related ARS/T were associated with patient survival. We built a web-application, CanSplice, to mine genomic datasets by self-identified race. The race-related targets have the potential to aid in the development of new biomarkers and therapeutics to mitigate cancer disparity.

Biological Aspects of Cancer Health Disparities.

Racial and ethnic disparities span the continuum of cancer care and are driven by a complex interplay among social, psychosocial, lifestyle, environmental, health system, and biological determinants of health. Research is needed to identify these determinants of cancer health disparities and to develop interventions to achieve cancer health equity. Herein, we focus on the overall burden of ancestry-related molecular alterations, the functional significance of the alterations in hallmarks of cancer, and the implications of the alterations for precision oncology and immuno-oncology. In conclusion, we reflect on the importance of estimating ancestry, improving diverse racial and ethnic participation in cancer clinical trials, and examining the intersection among determinants of cancer health disparities.

RNA splicing and aggregate gene expression differences in lung squamous cell carcinoma between patients of West African and European ancestry.

OBJECTIVES: Despite disparities in lung cancer incidence and mortality, the molecular landscape of lung cancer in patients of African ancestry remains underexplored, and race-related differences in RNA splicing remain unexplored. MATERIALS AND METHODS: We identified differentially spliced genes (DSGs) and differentially expressed genes (DEGs) in biobanked lung squamous cell carcinoma (LUSC) between patients of West African and European ancestry, using ancestral genotyping and Affymetrix Clariom D array. DSGs and DEGs were validated independently using the National Cancer Institute Genomic Data Commons. Associated biological processes, overlapping canonical pathways, enriched gene sets, and cancer relevance were identified using Gene Ontology Consortium, Ingenuity Pathway Analysis, Gene Set Enrichment Analysis, and CancerMine, respectively. Association with LUSC survival was conducted using The Cancer Genome Atlas. RESULTS: 4,829 DSGs and 267 DEGs were identified, including novel targets in NSCLC as well as genes identified previously to have relevance to NSCLC. RNA splicing events within 3 DSGs as well as 1 DEG were validated in the independent cohort. 853 DSGs and 29 DEGs have been implicated as potential drivers, oncogenes and/or tumor suppressor genes. Biological processes enriched among DSGs and DEGs included metabolic process, biological regulation, and multicellular organismal process and, among DSGs, ion transport. Overlapping canonical pathways among DSGs included neuronal signaling pathways and, among DEGs, cell metabolism involving biosynthesis. Gene sets enriched among DSGs included KRAS Signaling, UV Response, E2 F Targets, Glycolysis, and Coagulation. 355 RNA splicing events within DSGs and 18 DEGs show potential association with LUSC patient survival. CONCLUSION: These DSGs and DEGs, which show potential biological and clinical relevance, could have the ability to drive novel biomarker and therapeutic development to mitigate LUSC disparities.

Changing the landscape of non-small cell lung cancer disparities.

In the United States, lung and bronchus cancers are the second most common types of cancer and are responsible for the largest number of deaths from cancer, with African Americans suffering disproportionately from lung and bronchus cancers. This disparity likely results from a complex interplay among social, psycho-social, lifestyle, environmental, health system, and biological determinants of health. Toward improving outcomes for lung cancer patients of all races and ethnicities and mitigating lung cancer disparities, in this commentary, we bring forward biological factors that contribute to lung cancer disparities, efforts to identify, functionally characterize, and modulate novel ancestry-related RNA splicing-related targets in lung cancer for precision intervention, and translational and clinical research needs to improve outcomes for lung cancer patients of all races and ethnicities and mitigate lung cancer disparities.

Alternative RNA Splicing as a Potential Major Source of Untapped Molecular Targets in Precision Oncology and Cancer Disparities.

Studies of alternative RNA splicing (ARS) have the potential to provide an abundance of novel targets for development of new biomarkers and therapeutics in oncology, which will be necessary to improve outcomes for patients with cancer and mitigate cancer disparities. ARS, a key step in gene expression enabling individual genes to encode multiple proteins, is emerging as a major driver of abnormal phenotypic heterogeneity. Recent studies have begun to identify RNA splicing-related genetic and genomic variation in tumors, oncogenes dysregulated by ARS, RNA splice variants driving race-related cancer aggressiveness and drug response, spliceosome-dependent transformation, and RNA splicing-related immunogenic epitopes in cancer. In addition, recent studies have begun to identify and test, preclinically and clinically, approaches to modulate and exploit ARS for therapeutic application, including splice-switching oligonucleotides, small molecules targeting RNA splicing or RNA splice variants, and combination regimens with immunotherapies. Although ARS data hold such promise for precision oncology, inclusion of studies of ARS in translational and clinical cancer research remains limited. Technologic developments in sequencing and bioinformatics are being routinely incorporated into clinical oncology that permit investigation of clinically relevant ARS events, yet ARS remains largely overlooked either because of a lack of awareness within the clinical oncology community or perceived barriers to the technical complexity of analyzing ARS. This perspective aims to increase such awareness, propose immediate opportunities to improve identification and analysis of ARS, and call for bioinformaticians and cancer researchers to work together to address the urgent need to incorporate ARS into cancer biology and precision oncology.

Single-nucleotide polymorphisms of stemness genes predicted to regulate RNA splicing, microRNA and oncogenic signaling are associated with prostate cancer survival.

Prostate cancer (PCa) is a clinically and molecularly heterogeneous disease, with variation in outcomes only partially predicted by grade and stage. Additional tools to distinguish indolent from aggressive disease are needed. Phenotypic characteristics of stemness correlate with poor cancer prognosis. Given this correlation, we identified single-nucleotide polymorphisms (SNPs) of stemness-related genes and examined their associations with PCa survival. SNPs within stemness-related genes were analyzed for association with overall survival of PCa in the Prostate, Lung, Colorectal and Ovarian Cancer Screening Trial. Significant SNPs predicted to be functional were selected for linkage disequilibrium analysis and combined and stratified analyses. Identified SNPs were evaluated for association with gene expression. SNPs of CD44 (rs9666607), ABCC1 (rs35605 and rs212091) and GDF15 (rs1058587) were associated with PCa survival and predicted to be functional. A role for rs9666607 of CD44 and rs35605 of ABCC1 in RNA splicing regulation, rs212091 of ABCC1 in miRNA binding site activity and rs1058587 of GDF15 in causing an amino acid change was predicted. These SNPs represent potential novel prognostic markers for overall survival of PCa and support a contribution of the stemness pathway to PCa patient outcome.

Corrigendum: Alternative splicing promotes tumour aggressiveness and drug resistance in African American prostate cancer.

This corrects the article DOI: 10.1038/ncomms15921.

Alternative splicing promotes tumour aggressiveness and drug resistance in African American prostate cancer.

Clinical challenges exist in reducing prostate cancer (PCa) disparities. The RNA splicing landscape of PCa across racial populations has not been fully explored as a potential molecular mechanism contributing to race-related tumour aggressiveness. Here, we identify novel genome-wide, race-specific RNA splicing events as critical drivers of PCa aggressiveness and therapeutic resistance in African American (AA) men. AA-enriched splice variants of PIK3CD, FGFR3, TSC2 and RASGRP2 contribute to greater oncogenic potential compared with corresponding European American (EA)-expressing variants. Ectopic overexpression of the newly cloned AA-enriched variant, PIK3CD-S, in EA PCa cell lines enhances AKT/mTOR signalling and increases proliferative and invasive capacity in vitro and confers resistance to selective PI3Kδ inhibitor, CAL-101 (idelalisib), in mouse xenograft models. High PIK3CD-S expression in PCa specimens associates with poor survival. These results highlight the potential of RNA splice variants to serve as novel biomarkers and molecular targets for developmental therapeutics in aggressive PCa.

Associations between RNA splicing regulatory variants of stemness-related genes and racial disparities in susceptibility to prostate cancer.

Evidence suggests that cells with a stemness phenotype play a pivotal role in oncogenesis, and prostate cells exhibiting this phenotype have been identified. We used two genome-wide association study (GWAS) datasets of African descendants, from the Multiethnic/Minority Cohort Study of Diet and Cancer (MEC) and the Ghana Prostate Study, and two GWAS datasets of non-Hispanic whites, from the Prostate, Lung, Colorectal and Ovarian (PLCO) Cancer Screening Trial and the Breast and Prostate Cancer Cohort Consortium (BPC3), to analyze the associations between genetic variants of stemness-related genes and racial disparities in susceptibility to prostate cancer. We evaluated associations of single-nucleotide polymorphisms (SNPs) in 25 stemness-related genes with prostate cancer risk in 1,609 cases and 2,550 controls of non-Hispanic whites (4,934 SNPs) and 1,144 cases and 1,116 controls of African descendants (5,448 SNPs) with correction by false discovery rate ≤0.2. We identified 32 SNPs in five genes (TP63, ALDH1A1, WNT1, MET and EGFR) that were significantly associated with prostate cancer risk, of which six SNPs in three genes (TP63, ALDH1A1 and WNT1) and eight EGFR SNPs showed heterogeneity in susceptibility between these two racial groups. In addition, 13 SNPs in MET and one in ALDH1A1 were found only in African descendants. The in silico bioinformatics analyses revealed that EGFR rs2072454 and SNPs in linkage with the identified SNPs in MET and ALDH1A1 (r2 > 0.6) were predicted to regulate RNA splicing. These variants may serve as novel biomarkers for racial disparities in prostate cancer risk.