Abstract
Black men die from prostate cancer twice as often as White men, a disparity likely due to inherited genetics, modifiable cancer risk factors, and healthcare access. It is incompletely understood how and why tumor genomes differ by self-reported race and genetic ancestry.
Among 2,069 men with prostate cancer (1,841 self-reported White, 63 Asian, 165 Black) with access to clinical-grade sequencing at the same cancer center, prevalence of tumor and germline alterations was assessed in cancer driver genes reported to have different alteration prevalence by race.
Clinical characteristics such as prostate-specific antigen and age at diagnosis as well as cancer stage at sample procurement differed by self-reported race. However, most genomic differences persisted when adjusting for clinical characteristics. Tumors from Black men harbored fewer PTEN mutations and more AR alterations than those from White men. Tumors from Asian men had more FOXA1 mutations and more ZFHX3 alterations than White men. Despite fewer TP53 mutations, tumors from Black men had more aneuploidy, particularly chromosome arm 8q gains, an adverse prognostic factor. Genetic ancestry was associated with similar tumor alterations as self-reported race, but also with modifiable cancer risk factors. Community-level average income was associated with chr8q gains after adjusting for race and ancestry.
Tumor genomics differed by race even after accounting for clinical characteristics. Equalizing access to care may not fully eliminate such differences. Therapies for alterations more common in racial minorities are needed. Tumor genomic differences should not be assumed to be entirely due to germline genetics.
Prostate cancer has the greatest racial disparities in death rates in the United States of all cancers. As one contributor, tumor genomes have been suggested to differ by self-reported race. This study demonstrates that the same differences in tumor genomics as by self-reported race are observed if comparing men by genetic ancestry. There were notable differences in clinical factors by self-reported race among the men who underwent tumor-normal sequencing, but differences in tumor genomics persisted after adjusting for clinical characteristics. Tumors from Black men had more poor-prognosis chr8q gains, which beyond race and ancestry were also associated with a crude, area-level measure of socioeconomic status. These results underscore that therapies for genomic alterations more common in racial minorities are needed. Differences in tumor genomes should neither be assumed to result solely from access to care nor solely from germline genetics.
Introduction
Each year in the United States, there are 38 deaths from prostate cancer per 100,000 Black men, compared with 18 per 100,000 White men (1), a disparity that is greater than for any other cancer. It is critical to understand the contributions of the potential underlying causes, including differences in inherited (germline) genetics, modifiable risk factors, and healthcare access. Prostate cancer is highly heritable (2), and germline risk loci for prostate cancer differ by ancestry (3). In addition, the burden of modifiable risk factors for prostate cancer–related death, such as smoking (4), is higher among Black men (5). Racial and ethnic groups have unequal healthcare access to cancer screening, diagnostics, and treatment. Consequently, prostate cancer among Black men may be diagnosed and treated at later stages with shorter survival after diagnosis.
Another potential contributor to racial differences is underlying differences in molecular tumor features that may indicate different cancer etiologies as well as higher or lower sensitivity of the tumor to targeted therapies. Prostate tumors differ by race in key somatic DNA alterations (6–10). For example, in the largest study to date (8), metastatic prostate tumors from Black men, compared with those from White men, were reported to have more DNA mutations in the androgen receptor AR, in BRAF, in DNA repair genes, and in genes considered “actionable,” and fewer TP53 alterations than tumors from Asian men. These differences suggest that for some therapies, treatment benefits may be different by race. However, previous genomic studies had limited or no data on differences in healthcare access or cancer risk factors, and they typically relied on either race or genetic ancestry.
A first question is whether differences in tumor genomics are a consequence of disparities in healthcare access. For example, if Black men had tumor sequencing later in their disease course than White men, finding more AR alterations in their tumors could merely reflect barriers to healthcare access. Conversely, if genomic differences persisted once consequences of healthcare access were accounted for, then race may need to be considered for precision prevention and treatment selection (11, 12).
A second question is the interpretation of “race” and to what extent genetic ancestry is associated with genomic tumor differences associated with self-reported race. Race is a self-ascribed or socially-ascribed identity, whereas genetic ancestry describes components of an individual's genetic origin (13). Yet both concepts are highly correlated and must be carefully disentangled.
Finally, a third question is to what extent differences in tumor genomics are driven by germline genetics or could be reduced by acting on racial disparities in cancer risk factors. We addressed these questions in one of the largest, clinical-grade sequencing cohorts of prostate cancer.
Materials and Methods
Study population
This study included men with histologically confirmed prostate cancer who had access to care at a comprehensive cancer center and who reported race on intake forms (Supplementary Fig. S1). Patients gave written informed consent for this institutional review board–approved study that was conducted in accordance with the U.S. Common Rule. Freely available definitions of data elements for a clinical-genomic research database, their validation, and the analysis pipeline were described elsewhere (14). Area-level median household income per 5-digit ZIP code was from the U.S. Census 2015–2019 American Community Survey.
Genomics
Tumor profiling used archival or fresh formalin-fixed, paraffin-embedded tumor samples. Germline DNA from peripheral blood was only analyzed with specific consent. An FDA-cleared targeted hybridization capture-based gene panel, MSK-IMPACT, including 341 to 468 genes depending on its version, was used to capture single-nucleotide variants, small insertions and deletions, copy-number alterations, and structural rearrangements (15).
TMPRSS2 and ERG alterations, more common in men of European ancestry (16), were not included in analyses given the suboptimal sensitivity of DNA panel sequencing for gene fusions like TMPRSS2:ERG and the analytical focus on absolute differences. Oncogenic or likely oncogenic function was annotated using OncoKB 2.8 (September 17, 2020; ref. 17).
For copy-number gains and deletions, absolute log2 copy-number deviations >0.2 were required. The fraction of copy number–altered genome was defined as the proportion of the genome with such a copy-number change relative to the size of the genome with copy number profiled. Estimates were corrected for pathologist-assessed tumor cellularity; there were no differences in tumor cellularity by race (Supplementary Fig. S7). CNTools v1.46, copy number v1.30, and, for calling arm-level events, ASCETS (18) with breadth of coverage >0.5 were used.
To infer genetic ancestry, we used the 1000 Genomes Project (1KGP) with potentially admixed individuals excluded as reference. We selected 5,072 bi-allelic autosomal single-nucleotide polymorphism (SNP) markers from 1KGP with minor allele frequency >1% and within MSK-IMPACT468 bait intervals. For each patient, we genotyped these markers on the matched normal BAM files using GATK v4.0 Pileup, merged with 1KGP reference data and pruned with PLINK v1.9 (parameters −indep-pairwise 500kb 50 0.2) to keep only those markers that were in linkage equilibrium. We then performed supervised ADMIXTURE v1.3 (19) to estimate proportions of African (AFR), European (EUR), East Asian (EAS), NAM (Native American; not included in analyses), and South Asian (SAS) ancestries.
Statistical analyses
Prevalence differences for genomic alterations by race (categorical) or ancestry (each ancestry modeled as a continuous variable) were estimated using marginal standardization of logistic regression, with bias-corrected accelerated (BCa) 95% confidence intervals (CI) from parametric bootstrapping with 10,000 repeats (20). Regression models were used for differences in age at diagnosis and time to sequencing (quantile regression for the median), counts of actionable alterations (Poisson model), arithmetic mean ratios of fraction genome altered (generalized linear model with Gaussian distribution and log link), prevalence differences in chr8q gain by continuously modeled income (binomial model with identity link), and survival from sequencing to death or last contact (Cox regression with separate baseline hazards by disease extent at sequencing).
Data availability
Clinical and genomic data from this study are available via cBioPortal (https://www.cbioportal.org/study/summary?id=prad_msk_stopsack_2021).
Results
Clinical characteristics and cancer risk factors
We first compared clinical features of 2,069 men with prostate cancer by self-reported race including 63 Asian men, 165 Black men, and 1,841 White men (Fig. 1A; Table 1). Despite being 3.0 years younger at diagnosis (95% CI, 0.6 to 5.4), Black men had median serum prostate-specific antigen (PSA) levels at diagnosis of 17 ng/mL, higher than White men (9 ng/mL; difference, 8.2 ng/mL; 95% CI, 0.1 to 16.3). Asian men had a more advanced cancer stage at sample procurement (32% castration-resistant) than White men (19%; difference, 13 percentage points; 95% CI, 1 to 24). Time to tumor sequencing was similar by race in our study (Table 1).
. | . | By self-reported race . | ||
---|---|---|---|---|
. | Overall . | White . | Asian . | Black . |
N | 2,069 | 1,841 | 63 | 165 |
Age at diagnosis (years) | 63 (56–68) | 63 (56–68) | 63 (57–68) | 60 (53–64) |
Age at sample (years) | 66 (59–72) | 66 (59–72) | 66 (60–70) | 62 (56–69) |
Diagnosis to sample (months) | 3 (0–43) | 3 (0–42) | 4 (0–29) | 3 (0–50) |
Diagnosis to sequencing (months) | 26 (6–83) | 27 (6–83) | 18 (7–60) | 24 (6–84) |
Smoking status | ||||
Never | 1,099 (55%) | 981 (55%) | 38 (60%) | 80 (49%) |
Former | 783 (39%) | 698 (39%) | 24 (38%) | 61 (38%) |
Current | 132 (7%) | 110 (6%) | 1 (2%) | 21 (13%) |
Unknown | 55 | 52 | 0 | 3 |
Area-level household income | 115,844 | 119,884 | 119,780 | 72,568 |
($/year) | (84,643–156,558) | (90,293–158,420) | (89,288–160,935) | (54,886–96,086) |
Unknown | 40 | 39 | 0 | 1 |
PSA at diagnosis (ng/mL) | 9.2 (5.3–26.8) | 8.8 (5.2–23.3) | 17.9 (6.8–97.0) | 17.0 (7.0–94.5) |
Unknown | 116 | 104 | 2 | 10 |
Gleason score | ||||
<7 | 171 (9%) | 159 (9%) | 3 (5%) | 9 (6%) |
3+4 | 287 (15%) | 249 (15%) | 9 (16%) | 29 (19%) |
4+3 | 336 (18%) | 294 (17%) | 12 (21%) | 30 (19%) |
8 | 381 (20%) | 338 (20%) | 11 (19%) | 32 (21%) |
9–10 | 726 (38%) | 648 (38%) | 23 (40%) | 55 (35%) |
Unknown | 168 | 153 | 5 | 10 |
Stage (M) | ||||
M0 | 1,501 (73%) | 1,341 (73%) | 39 (62%) | 121 (73%) |
M1 | 568 (27%) | 500 (27%) | 24 (38%) | 44 (27%) |
Sample tissue | ||||
Prostate | 1,352 (65%) | 1,217 (66%) | 34 (54%) | 101 (61%) |
Lymph node | 331 (16%) | 284 (15%) | 15 (24%) | 32 (19%) |
Bone | 191 (9%) | 166 (9%) | 8 (13%) | 17 (10%) |
Visceral | 125 (6%) | 110 (6%) | 3 (5%) | 12 (7%) |
Other soft tissue | 70 (3%) | 64 (3%) | 3 (5%) | 3 (2%) |
Extent of disease at sample | ||||
Localized | 700 (34%) | 630 (34%) | 16 (25%) | 54 (33%) |
Regional nodes | 362 (17%) | 324 (18%) | 9 (14%) | 29 (18%) |
Metastatic hormone-sensitive | 578 (28%) | 512 (28%) | 18 (29%) | 48 (29%) |
Nonmetastatic castration-resistant | 17 (1%) | 16 (1%) | 0 (0%) | 1 (1%) |
Metastatic castration-resistant | 405 (20%) | 353 (19%) | 20 (32%) | 32 (19%) |
Metastatic, variant histology | 7 (0%) | 6 (0%) | 0 (0%) | 1 (1%) |
Tumor cellularity by pathologist (%) | 40 (20–50) | 40 (20–50) | 40 (20–50) | 40 (20–50) |
Unknown | 11 | 9 | 1 | 1 |
Mean sequencing coverage (X) | 611 (461–754) | 612 (461–753) | 648 (539–793) | 581 (453–737) |
Consent for germline analyses | 1,781 (86%) | 1,583 (86%) | 54 (86%) | 144 (87%) |
. | . | By self-reported race . | ||
---|---|---|---|---|
. | Overall . | White . | Asian . | Black . |
N | 2,069 | 1,841 | 63 | 165 |
Age at diagnosis (years) | 63 (56–68) | 63 (56–68) | 63 (57–68) | 60 (53–64) |
Age at sample (years) | 66 (59–72) | 66 (59–72) | 66 (60–70) | 62 (56–69) |
Diagnosis to sample (months) | 3 (0–43) | 3 (0–42) | 4 (0–29) | 3 (0–50) |
Diagnosis to sequencing (months) | 26 (6–83) | 27 (6–83) | 18 (7–60) | 24 (6–84) |
Smoking status | ||||
Never | 1,099 (55%) | 981 (55%) | 38 (60%) | 80 (49%) |
Former | 783 (39%) | 698 (39%) | 24 (38%) | 61 (38%) |
Current | 132 (7%) | 110 (6%) | 1 (2%) | 21 (13%) |
Unknown | 55 | 52 | 0 | 3 |
Area-level household income | 115,844 | 119,884 | 119,780 | 72,568 |
($/year) | (84,643–156,558) | (90,293–158,420) | (89,288–160,935) | (54,886–96,086) |
Unknown | 40 | 39 | 0 | 1 |
PSA at diagnosis (ng/mL) | 9.2 (5.3–26.8) | 8.8 (5.2–23.3) | 17.9 (6.8–97.0) | 17.0 (7.0–94.5) |
Unknown | 116 | 104 | 2 | 10 |
Gleason score | ||||
<7 | 171 (9%) | 159 (9%) | 3 (5%) | 9 (6%) |
3+4 | 287 (15%) | 249 (15%) | 9 (16%) | 29 (19%) |
4+3 | 336 (18%) | 294 (17%) | 12 (21%) | 30 (19%) |
8 | 381 (20%) | 338 (20%) | 11 (19%) | 32 (21%) |
9–10 | 726 (38%) | 648 (38%) | 23 (40%) | 55 (35%) |
Unknown | 168 | 153 | 5 | 10 |
Stage (M) | ||||
M0 | 1,501 (73%) | 1,341 (73%) | 39 (62%) | 121 (73%) |
M1 | 568 (27%) | 500 (27%) | 24 (38%) | 44 (27%) |
Sample tissue | ||||
Prostate | 1,352 (65%) | 1,217 (66%) | 34 (54%) | 101 (61%) |
Lymph node | 331 (16%) | 284 (15%) | 15 (24%) | 32 (19%) |
Bone | 191 (9%) | 166 (9%) | 8 (13%) | 17 (10%) |
Visceral | 125 (6%) | 110 (6%) | 3 (5%) | 12 (7%) |
Other soft tissue | 70 (3%) | 64 (3%) | 3 (5%) | 3 (2%) |
Extent of disease at sample | ||||
Localized | 700 (34%) | 630 (34%) | 16 (25%) | 54 (33%) |
Regional nodes | 362 (17%) | 324 (18%) | 9 (14%) | 29 (18%) |
Metastatic hormone-sensitive | 578 (28%) | 512 (28%) | 18 (29%) | 48 (29%) |
Nonmetastatic castration-resistant | 17 (1%) | 16 (1%) | 0 (0%) | 1 (1%) |
Metastatic castration-resistant | 405 (20%) | 353 (19%) | 20 (32%) | 32 (19%) |
Metastatic, variant histology | 7 (0%) | 6 (0%) | 0 (0%) | 1 (1%) |
Tumor cellularity by pathologist (%) | 40 (20–50) | 40 (20–50) | 40 (20–50) | 40 (20–50) |
Unknown | 11 | 9 | 1 | 1 |
Mean sequencing coverage (X) | 611 (461–754) | 612 (461–753) | 648 (539–793) | 581 (453–737) |
Consent for germline analyses | 1,781 (86%) | 1,583 (86%) | 54 (86%) | 144 (87%) |
aValues are count (percent) or median (quartile 1–quartile 3).
These differences in clinical cancer features prompted us to investigate two examples of potential upstream factors (Fig. 1B). Socioeconomic status has profound influences on cancer risk and survival (21) through cancer risk factors and healthcare use and access. Black men lived in ZIP code areas with an average household income 36% lower than that of White men (95% CI, 30% to 42%). There were no noticeable differences between Asian and White men. Twice as many Black men (13%) as White men (6%) were current smokers at the time of cancer diagnosis (difference, 7 percentage points; 95% CI, 2 to 12), a strong risk factor for prostate cancer–related death (4).
Tumor genomics by self-reported race before and after accounting for clinical features
We quantified alteration prevalence in oncogenes and tumor suppressors in tumor DNA (Fig. 2; Supplementary Fig. S2). Our focus was on genes that putatively differed by race in previous studies (7–10). Unlike previous studies, we considered all types of DNA alterations (mutations, copy-number alterations, and fusions) but limited analyses to those classified as oncogenic or likely oncogenic based on a curated precision oncology database (17). We first compared expected tumor genomics by self-reported race because genetic ancestry is typically unknown in clinical practice until after a sequencing test is performed.
Alterations in the androgen receptor gene (AR) occurred in 19% of tumors from Asian men and 18% among Black men, as compared with 12% among White men. Potential differences in AR alteration prevalence between tumors from Asian and White men were completely attenuated (prevalence difference, 0 percentage points; 95% CI, −5 to 6) after adjusting for clinical features (age, disease extent, PSA, stage, sample type) and thus accounting for the more advanced disease state when samples had been procured from Asian men. In contrast, AR alterations were slightly more common in tumors from Black men than White men even after adjusting for clinical features (difference, 5 percentage points; 95% CI, 1 to 8).
Other differences in tumor genomics by race persisted when adjusting for clinical features (Fig. 2), suggesting that the observed differences in tumor genomics were not driven exclusively by differences in disease extent. The strongest differences were in FOXA1, ZFHX3, TP53, and PTEN, with a smaller difference in MYC (copy-number gains) and BRAF alterations (Fig. 2).
In FOXA1, 27% of tumors from Asian men had alterations compared to 13% in White men (adjusted difference, 15 percentage points; 95% CI, 6 to 26), a difference driven by class-1 FOXA1 mutations in the forkhead domain (22), which occurred in 21% of tumors from Asian men compared with 8% from White men (Supplementary Fig. S3A). ZFHX3 alterations, both mutations and copy-number alterations, were enriched in tumors from Asian men (17%) compared with White men (7%; adjusted difference, 8 percentage points; 95% CI, 1 to 17; Supplementary Fig. S2).
TP53 alterations were less common in tumors from Black men (19%) than White men (30%; adjusted difference, −10 percentage points; 95% CI, −3 to −17). This difference was likely driven by fewer TP53 missense mutations among tumors from Black men (11%) than White men (18%; adjusted difference, −4 percentage points; 95% CI, −9 to 1; Supplementary Fig. S3B). Differences were not explained by frameshift or truncating mutations or single hotspots, such as the most common mutation at amino acid R273 (prevalence among White men, 2%).
Alterations in the tumor suppressor PTEN were less common in tumors from Asian and Black men (14%) compared with White men (20%; adjusted difference Black vs. White, −7 percentage points; 95% CI, −12 to −1), particularly due to fewer PTEN mutations (8% in Black vs. 3% in White; adjusted difference, −4 percentage points; 95% CI, −7 to 0; Supplementary Fig. S3C).
In contrast to a previous study (8), we did not observe differences in the prevalence of DNA repair alterations (Fig. 2) or other potentially actionable alterations, nor did we find the (low) prevalence of microsatellite instability and the tumor mutation burden to differ by race (Supplementary Table S1). When combining gene-level alterations in oncogenic signaling pathways (23), differences in PI3K, RAS/RAF/MAPK, and MYC signaling pathways (Supplementary Fig. S4) were consistent with the observed differences in PTEN, BRAF, and MYC, respectively. In addition to these confirmatory analyses, we conducted an exploratory analysis of frequently altered genes (Supplementary Fig. S5), without additional candidates.
Tumor genomics and cancer risk factors by genetic ancestry
SNPs captured by the sequencing panels allowed us to infer genetic ancestry, which was strongly associated with but, as expected, not perfectly predictive of self-reported race (Fig. 3A). We then assessed how higher proportions of specific genetic ancestries were associated with tumor alterations that differed by self-reported race. We qualitatively observed the same associations between genetic ancestry and tumor alterations as between self-reported race and tumor alterations (Fig. 3B). For example, TP53 alteration prevalence was −13 percentage points lower (95% CI, −19 to −5) in tumors from men with 100% AFR genetic ancestry markers than with 0% AFR ancestry.
We then probed to what extent ancestry can be interpreted causally as the influence of germline genetics or whether ancestry was also associated with nongenetic factors. We again encountered differences in smoking prevalence and in community-wide median household income, a unidimensional and imprecise measure of socioeconomic status, when comparing men by genetic ancestry (Fig. 3C). Among men with 100% AFR ancestry, 15% were predicted to live in a community with a median household income of at least $100,000 annually, compared with 68% among men with 0% AFR ancestry (difference, 52 percentage points; 95% CI, 46 to 58).
We assessed to what extent differences in tumor alteration prevalence could directly result from alterations present in the germline, which were included in all analyses. Consent was given by 1,781 (86%) men for germline genetic testing for alterations in genes associated with heritable cancer risk, without differences in consent by self-reported race. While germline alterations were included in all analyses, none of the genes with prevalence differences by race and ancestry had more than a single germline mutation among the 1,781 patients (Supplementary Table S2).
TP53 alterations and genomic instability
Tumor aneuploidy is strongly tied to TP53 alterations and genomic instability (24). Despite fewer TP53 alterations, tumors from Black men did not have a lower proportion of the tumor genome affected by copy-number alterations (Fig. 4A). When adjusting for TP53 alteration status, tumors from Black men, but not Asian men, had slightly more copy-number alterations than those from White men.
Patterns of copy-number alterations appeared generally similar by race, except for chromosome arm 8q gains (Fig. 4B; Supplementary Fig. S6). Forty-nine percent of tumors from Black men and 47% of tumors from Asian men had chr8q gains, as compared with 37% of tumors from White men. Adjusting for clinical factors attenuated the difference in chr8q gain prevalence for Asian men (difference to White, 5 percentage points; 95% CI, −6 to 17) but not for Black men (difference to White, 11 percentage points; 95% CI, 4 to 18; Fig. 4C).
As expected, chr8q gains were an adverse prognostic factor. Over a median of 30 months, 441 deaths occurred among 1,690 White men with follow-up data after sequencing, 15 among 58 Asian men (precluding survival analyses), and 45 among 156 Black men. These data were not indicative of a generally worse prognosis among Black men in our study population, but we cannot rule out a difference (HR compared with White men, 1.15; 95% CI, 0.85 to 1.57). Tumors with chr8q gain, compared with tumors without chr8q gain but the same disease extent at sequencing, had a worse prognosis both among White men (HR, 1.61; 95% CI, 1.33 to 1.96) and among Black men (HR, 2.00; 95% CI, 1.00 to 4.01).
The difference in chr8q gain prevalence by race was mirrored by genetic ancestry. Tumors from men with 100% AFR ancestry had a 13 percentage-point higher prevalence of chr8q gains compared with 0% AFR ancestry (95% CI, 4 to 23). Yet even after adjusting for both self-reported race and genetic ancestry on an individual level, area-level median household income was associated with chr8q gain (Fig. 4D). For each $50,000/year increase, men of the same self-reported race and genetic ancestry were 2 percentage points less likely to have tumors with chr8q gain (95% CI, −4 to 0). Conversely, adjusting for community-level income reduced the above-noted difference in chr8q gains from 13 to 10 percentage points (95% CI, 0 to 20).
Discussion
Even within a group of men that ultimately had access to healthcare and tumor sequencing at the same cancer center, clinical factors and cancer risk factors differed noticeably by self-reported race and ancestry. These features are an essential part of the perspective on differences in tumor genomics. Tumor genomics still differed between men with similar prostate cancer clinical features but different self-reported race or genetic ancestry, in particular, prevalence of AR and BRAF alterations, class-1 FOXA1 mutations, ZFHX3 alterations, mutations of PTEN, missense mutations in TP53, and gains of chromosome arm 8q. Despite fewer TP53 alterations, tumors from Black men and those with greater AFR ancestry had more aneuploidy at chr8q, which may contribute to their aggressiveness. Our results suggest that self-reported race carries some information on what tumor genomic profiles to expect in clinic before a sequencing test is obtained. They also highlight the need for targeted treatments for genomic alterations more common in non-White men.
Our findings on differences in alteration prevalence of commonly altered genes in prostate cancer are generally consistent with previous studies (8–10, 25). Some men in this study were included in a previous study (8) through Project GENIE (26). Our data are thus not an independent validation that alteration prevalence differs. Instead, the importance of our findings is the demonstration that most differences persist when restricting to oncogenic driver alterations, when accounting for different clinical features between men with different race, and when considering genetic ancestry.
The higher prevalence of AR and BRAF alterations in tumors from Black men is directly relevant in light of currently available targeted therapies (8). FOXA1, in which class-1 mutations in the DNA-binding forkhead domain were more common among tumors from Asian men, has been described functionally (22) but is not yet therapeutically actionable. In contrast to the class-1 FOXA1 mutations that occur early during carcinogenesis (22), a previously noted higher prevalence of (generally late-occurring) AR alterations in tumors from Asian men (8) was explained by more advanced disease among Asian men and thus should not influence clinical decision making. A lower prevalence of PTEN alterations in non-White men is relevant for PI3K inhibitors, which have reached phase III testing in prostate cancer (27).
Ancestry differences in aneuploidy, albeit in varying directions, have been observed in other cancer types (28). In our study, tumors with fewer TP53 alterations in Black men, specifically missense mutations, would be expected to have fewer chromosome arm losses and gains, given the strong association between TP53 alterations and aneuploidy burden (24). Surprisingly, tumors from Black men and those with higher AFR ancestry were no less aneuploid and instead harbored more chr8q gains, which are a poor prognostic indicator in primary prostate cancer (29) as well as in Black and White men in our study. chr8q gains explain differences in MYC alterations. Therapeutic targets on chr8q beyond MYC might be particularly relevant for Black men.
Many germline risk loci for prostate cancer are located at chr8q24, and some differ by ancestry (30). In our study, genetic ancestry had similar associations with chr8q gain and the other somatic alterations as did self-reported race. Germline genetics are clearly a contributor to tumor genomics. However, chr8q gains were also associated with a crude measure of socioeconomic status of neighborhoods, which differed strongly between Black and White men and by genetic ancestry. This crude, ecological measure does not allow pinpointing causes, which may include limited education opportunities and social services, neighborhood violence, environmental pollution, food insecurity, and many others (12), in addition to individual-level risk factors. Yet that modifiable risk factors have a strong impact is evidenced by prostate cancer rates among Asian-American men, which are lower than among Whites but substantially higher than for men of Asian ancestry living in Asia (31)—without changes in germline genetics upon immigration to the U.S. Our findings illustrate that studies on genetic ancestry as a potential cause of differences in tumor genomics need to consider confounding by nongenetic factors.
The underrepresentation of Asian and Black men in our study partially reflects disparities in access to healthcare. Increasing representation of racial and ethnic minorities will be important to allow for assessing differences by Hispanic ethnicity, between EAS and SAS ancestries, and by ancestry within each of the ancestry superpopulations, particularly between Northern, Eastern, and Western African ancestry that have been subsumed under AFR (11, 32). Other facets of tumor genomics, and a richer set of cancer risk factors, need to be assessed. Even more important are population-based, rather than hospital-based study designs, without which we cannot know the denominator of men at risk of developing cancers with certain genomic features—a prerequisite for causal attributions and thus interventions that address causes, rather than just correlates of racial differences.
Taken together, tumor genomes differed by self-reported race even after accounting for differences in clinical features between men of different race. Cautiously considering race and tumor genomics jointly has potential to contribute to reducing the significant disparities in prostate cancer death rates that exist in the United States today (13). Differences in tumor genomics by race and ancestry should not be assumed to be completely genetically determined and thus immutable in the future.
Authors' Disclosures
S.M. McBride reports grants from Janssen, AstraZeneca, and Roche outside the submitted work. M.J. Morris reports personal fees from Oric Pharmaceuticals, Curium, Athenex, Exelexis, and AstraZeneca outside the submitted work. D.E. Rathkopf reports uncompensated consulting or advisory roles for Genentech/Roche, Janssen Oncology, AstraZeneca, and Myovant, as well as institutional research funding from AstraZeneca, BMS/Celgene, Ferring, Genentech/Roche, Janssen Oncology, Medivation, Millennium, Novartis, Taiho Pharmaceutical, Takeda, TRACON Pharma, and Phosplatin. D.C. Danila reports research support from US Department of Defense, American Society of Clinical Oncology, Prostate Cancer Foundation, Stand Up 2 Cancer, Janssen Research & Development, Astellas, Medivation, Agensys, Genentech, and CreaTV, as well as consultancy for Angle LLT, Axiom LLT, Janssen Research & Development, AstraZeneca, Clovis, Astellas, Medivation, Pfizer, and Agensys. K.A. Autio reports grants from AstraZeneca, Pfizer, Amgen, and Parker Institute for Cancer Immunotherapy outside the submitted work. H.I. Scher reports grants from NCI during the conduct of the study. H.I. Scher also reports personal fees from Asterias Biotherapeutics, Bayer, Ambry Genetics Corp/Konica Minolta Inc, Pfizer Inc, Sun Pharmaceuticals Industries Inc, WCG Oncology, and Elsevier Ltd; grants from Epic Sciences, Illumina Inc, Janssen, Menarini Silicon Biosystems, and Thermo Fisher; and non-financial support from Amgen, ESSA Pharma, Janssen, Menarini Silicon Biosystems, Amgen, Asterias Biotherapeutics, Bayer, ESSA Pharma Inc, Phosplatin, Pfizer Inc, and WCG Oncology outside the submitted work. L.A. Mucci reports grants from AstraZeneca and Janssen, as well as personal fees from Bayer outside the submitted work. D.B. Solit reports personal fees from Pfizer, Loxo/Lilly Oncology, Vividion Therapeutics, Scorpion Therapeutics, FORE Therapeutics, and BridgeBio outside the submitted work. Y. Chen reports grants from NCI during the conduct of the study, as well as other support from Oric Pharmaceuticals outside the submitted work. M.F. Berger reports personal fees from Eli Lilly and personal fees from PetDx outside the submitted work. W. Abida reports personal fees and other support from Clovis Oncology and ORIC Pharmaceuticals; personal fees from Janssen, Daiichi Sankyo, OncLive, Roche, Medscape, and Aptitute Health; other support from AstraZeneca, Zenith Epigenetics, and Epizyme outside the submitted work. No disclosures were reported by the other authors.
Authors' Contributions
K.H. Stopsack: Conceptualization, data curation, software, formal analysis, funding acquisition, visualization, writing–original draft, writing–review and editing. S. Nandakumar: Data curation, investigation, visualization, writing–review and editing. K. Arora: Data curation, formal analysis, methodology, writing–review and editing. B. Nguyen: Data curation, investigation, writing–review and editing. S.E. Vasselman: Data curation, writing–review and editing. B. Nweji: Data curation, writing–review and editing. S.M. McBride: Resources, writing–review and editing. M.J. Morris: Resources, writing–review and editing. D.E. Rathkopf: Resources, writing–review and editing. S.F. Slovin: Resources, writing–review and editing. D.C. Danila: Resources, writing–review and editing. K.A. Autio: Resources, writing–review and editing. H.I. Scher: Resources, funding acquisition, writing–review and editing. L.A. Mucci: Investigation, writing–review and editing. D.B. Solit: Resources, funding acquisition, writing–review and editing. M. Gönen: Validation, writing–review and editing. Y. Chen: Resources, writing–review and editing. M.F. Berger: Methodology, writing–review and editing. N. Schultz: Resources, supervision, investigation, writing–review and editing. W. Abida: Resources, data curation, supervision, funding acquisition, writing–review and editing. P.W. Kantoff: Resources, supervision, funding acquisition, writing–review and editing.
Acknowledgments
This work was funded in part by the NCI (1P01CA228696, to P.W. Kantoff; P30CA008748, Cancer Center Support grant; P50CA092629, Prostate Cancer SPORE), the Department of Defense (Early Investigator Research Award W81XWH-18–1-0330, to K.H. Stopsack; Physician Research Award W81XWH-17–1-0124, to W. Abida), and the Marie-Josée and Henry R. Kravis Center for Molecular Oncology. K.H. Stopsack, S.M. McBride, D.E. Rathkopf, L.A. Mucci, N. Schultz, and W. Abida are Prostate Cancer Foundation Young Investigators. We are grateful for contributions of the Molecular Diagnostics Service and Marc Ladanyi, MD, to the MSK-IMPACT platform.
The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked advertisement in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.
Note: Supplementary data for this article are available at Clinical Cancer Research Online (http://clincancerres.aacrjournals.org/).