Breast and prostate cancers harbor common somatic copy number alterations that consistently differ by race and are associated with survival

Pan-cancer studies of somatic copy number alterations (SCNAs) have demonstrated common SCNA patterns across cancer types, but despite demonstrable differences in aggressiveness of some cancers by race, pan-cancer SCNA variation by race has not been explored. This study investigated a) racial differences in SCNAs in both breast and prostate cancer, b) the degree to which they are shared across cancers, and c) the impact of these shared, race-differentiated SCNAs on cancer survival. Utilizing data from The Cancer Genome Atlas (TCGA), SCNAs were identified using GISTIC 2.0, and in each tumor type, differences in SCNA magnitude between African Americans (AA) and European Americans (EA) were tested using linear regression. Unsupervised hierarchical clustering of the copy number of genes residing in race-differentiated SCNAs shared between tumor types was used to identify SCNA-defined patient groups, and Cox proportional hazards regression was used to test for association between those groups and overall/progression-free survival (PFS). We identified SCNAs that differed by race in breast (n = 58 SCNAs; permutation p < 10− 4) and prostate tumors (n = 78 SCNAs; permutation p = 0.006). Six race-differentiated SCNAs common to breast and prostate found at chromosomes 5q11.2-q14.1, 5q15-q21.1, 8q21.11-q21.13, 8q21.3-q24.3, 11q22.3, and 13q12.3-q21.3 had consistent differences by race across both tumor types, and all six were of higher magnitude in AAs, with the chromosome 8q regions being the only amplifications. Higher magnitude copy number differences in AAs were also identified at two of these race-differentiated SCNAs in two additional hormonally-driven tumor types: endometrial (8q21.3-q24.3 and 13q12.3-q21.3) and ovarian (13q12.3-q21.3) cancers. Race differentiated SCNA-defined patient groups were significantly associated with survival differences in both cancer types, and these groups also differentiated within triple negative breast cancers based on PFS. While the frequency of the SCNA-defined patient groups differed by race, their effects on survival did not. This study identified race-differentiated SCNAs shared by two related cancers. The association of SCNA-defined patient groups with survival demonstrates the clinical significance of combinations of these race-differentiated genomic aberrations, and the higher frequency of these alterations in AA relative to EA patients may explain racial disparities in risk of aggressive breast and prostate cancer.


Background
Both breast and prostate cancer exhibit racial disparities in incidence and outcomes that could be tied to underlying differences in genetics and disease biology. Relative to European American (EA) men, African American (AA) men are at a 1.7-fold increased risk of being diagnosed with prostate cancer and a 2.3-fold increased risk of dying from the disease [1]. While breast cancer incidence is slightly higher among EA women compared with AAs, mortality rates from this disease are more than 40% higher among AA women [1]. This mortality gap for AAs with breast cancer may be explained in part by the differential incidence of biologically-aggressive tumors that are negative for the estrogen receptor, the progesterone receptor, and the ERBB2 gene (commonly called triple negative breast cancer, TBNC), a phenotype that is twice as likely to occur in AA compared to EA women [2]. While molecular differences between these two racial groups have been investigated separately in both breast and prostate cancer, the degree to which somatic alterations that differ by race are shared across these two hormonally-driven cancers has not.
As breast and prostate cancer share some of the same biologic pathways, such as the steroid-hormone metabolism and insulin-like growth factor signaling, which when altered can lead to either cancer, a biologic rationale exists to investigate shared somatic genetic mechanisms between these two hormonally-driven cancers [3]. This concept of a shared somatic genetic etiology between breast and prostate cancer is supported by a recent study which used a gene expression panel, the PAM50, developed to molecularly subtype breast cancer, to classify prostate tumors into luminal-and basal-like subtypes [4]. Recent pan-cancer analyses have also shown breast and prostate cancer to have similar levels of insertion/deletion mutations [5] and intra-tumor heterogeneity [6], suggesting a similar mutational landscape may exist between the two cancers.
At a somatic genomic level, the mutational landscape of breast and prostate cancers are characterized by a large array of point mutations, but most tumors also demonstrate significant structural rearrangements of chromosomes and large segments of DNA gain and loss [7,8]. These somatic copy number alterations (SCNAs) are associated with disease recurrence and survival in both breast [9][10][11] and prostate cancer [12][13][14]; however, these SCNA studies have included subjects principally of European ancestry. We have previously shown that SCNA biomarkers for aggressive prostate cancer established in EA populations may have utility in AAs [15], but the observed racial differences in SCNA frequencies in prostate cancer [16] suggest that larger somatic molecular studies of AAs are needed to better characterize these differences in frequency that may be related to disparities in disease outcome. Race-differentiated SCNAs have been shown to lead to gene expression differences in tumor immune response in prostate cancer [16] and associate with previously defined molecular subtypes in breast cancer [17].
While ethnic differences in SCNAs have been investigated before in both breast and prostate cancer, the datasets are typically small in size and with limited validation of the findings. Employing a novel approach that quantifies the magnitude of SCNA alteration, we used data from AAs and EAs in The Cancer Genome Atlas (TCGA) to identify SCNAs that 1) differ by race in breast and prostate cancer separately and 2) validate across these two hormonally-driven cancers. We further validated the race-differentiated SCNAs shared between breast and prostate cancer in two other hormonally driven cancersendometrial and ovarian. We also tested the hypothesis that germline genetic African ancestry at an SCNA in TCGA admixed AAs (i.e. those whose genomes are comprised of ancestry from both Africa and Europe) could explain the observed differences in magnitude of SCNAs between AAs and EAs in both breast and prostate cancer. Finally, using an unsupervised approach, we evaluated whether underlying combinations of race-differentiated SCNAs defined patient groups that were associated with survival, and using gene expression data, we identified genes where expression was associated with both race and copy number within race-differentiated SCNAs.

Study subjects and molecular data
Data from The Cancer Genome Atlas (TCGA) were used for this study. TCGA somatic copy number (Affymetrix SNP 6 array) and germline genotype (Affymetrix SNP 6 array) data from the 22 autosomal chromosomes were downloaded from TCGA and are available in dbGaP, accession number phs000178.v11.p8 (https://ftp.ncbi.nlm.nih.gov/dbgap/stud ies/phs000178/phs000178.v11.p8/). These molecular data types were measured on 893 female breast cancer cases (719 EAs and 174 AAs) and 313 prostate cancer cases (270 EAs and 43 AAs). Self-reported AAs with a genome-wide African ancestry proportion (estimated as described in the Supplementary Methods) less than 20% and self-reported EAs with a genome-wide African ancestry proportion greater than 20% were removed from analysis. Samples with missing age-at-diagnosis and tumor severity data were also removed. After these quality control steps, 886 breast cancers and 309 prostate cancers remained, and the cancer specific characteristics are shown in Table 1.

Statistical methods
Regional chromosomal ancestry was inferred using RFMix [18] (see Supplementary Methods). Genome-wide SCNA burden for each tumor was calculated as the percentage of the genome where the log 2 (copy number/2) differed (either positive or negative) from the tumor median by more than 0.1. Overall survival (OS) and progression-free survival Table 1 Characteristics of TCGA African American (AA) and European American (EA) breast and prostate tumors (PFS) were downloaded from the pan-cancer TCGA publication [19] and events were censored at 10 years.
To identify commonly altered SCNAs, Gistic2 was utilized [20]. Briefly, Gistic2 models log 2 (copy number/2) at each SNP location across tumor genomes and identifies frequent SCNAs through the use of both the frequency and quantitative degree of a given alteration in the sample of tumors being analyzed. For each tumor type, SCNAs were identified using Gistic2 in AAs and EAs separately to account for the differences in samples size (i.e. > 3 times more EAs in both tumor types). The resulting non-identical SCNA boundaries between AAs and EAs were resolved such that a partially overlapping region could be divided into as many as three subregions based on the overlapping status: identified only in AA, identified only in EA, or identified in both groups. The approach to the reconciliation of the degree of overlap between SCNAs between the two groups is detailed in the Supplementary Methods and Figure S1.
To identify SCNAs that differ between AA and EA individuals, we quantified the copy number magnitude across each SCNA by calculating the area under the copy number curve (cnAUC) (See Supplementary Methods and Figure  S2). For this measure, positive area (above the null value of two copies) corresponds to amplification and negative area (below the null value of two copies) corresponds to deletion.
To assess the association between cnAUC and race for each SCNA, the following linear model was used: For the dichotomous race variable, EA race was treated as the referent group. As such, a positive race difference in amplified SCNAs and a negative race difference in loss SCNAs represent higher magnitude copy number amplifications and deletions, respectively, in AA relative to EA tumors. These models were also adjusted for age at diagnosis and tumor pathology (pathological stage for the breast, endometrial, ovarian, lung and kidney cancer analyses; and Gleason grade for the prostate cancer analyses). To be inclusive and account for the lower number of AA prostate tumors (n = 42) relative to AA breast tumors (n = 172), a p-value threshold of 0.1 was used to identify potential SCNAs with differences between races. To allow the comparison of race coefficients of SCNAs with different lengths, the coefficients were standardized by each SCNA's length. To assess the statistical significance of the total number of racedifferentiated SCNAs observed both within each tumor type and shared between tumor types, permutation tests were performed, with details presented in Supplementary Methods. A similar linear model was used to quantify the relationship between SCNAs and regional ancestry (see Supplementary Methods) To identify underlying patient groups based on the race-differentiated SCNAs shared by both breast and prostate cancer, unsupervised hierarchical clustering, using the complete linkage algorithm and Euclidean distance, was applied to the log 2 (copy number/2) of genes residing in the six consistently race-differentiated SCNAs. The gap statistic [21] was used to determine the optimal number of clusters, where the optimal number of clusters, k, is the first k where k + 1 clusters does not lead to an increase of the gap statistic. Three patient groups in both breast and prostate cancer were determined to be optimal and were visualized using a copy number-based heatmap. Chi-square tests were used to assess the association between patient group and race or TNBC status, and log-rank tests were used to test for differences in Kaplan-Meier survival curves among different groups. Cox proportional hazards models were used to test the association of patient group with 10-year OS and PFS while adjusting for race, age-at-diagnosis, tumor pathology, and SCNA burden.
To investigate differences in transcriptional gene expression within the six consistently race-differentiated SCNAs, the raw, gene-level mapped Illumina RNASeqV2 sequencing reads were analyzed using the R limma package [22]. Briefly, the raw read counts were normalized by total number of mapped reads in each sample, yielding counts per million (cpm),and further log transformed (log-cpm). The log-cpm values allow for the use of standard linear regression modeling to identify differentially expressed genes by race. Age-at-diagnosis, tumor pathology (pathological stage for breast cancer and Gleason grade for prostate cancer) were also included in the model. The Benjamini-Hochberg false discovery rate (FDR) was used to adjust p-values for multiple comparisons. For those genes with significant differential expression by race, we also assessed whether copy number at the gene-level was associated with gene expression using linear regression.

SCNAs in TCGA breast tumors and differences by race
Given the higher frequency of tumor data for EAs in comparison to AAs in TCGA breast and prostate tumor sets, we elected to identify SCNAs separately by race group in order not to obscure SCNAs that might be specific to AAs. From the 714 EA TCGA breast tumors (Table 1), Gistic2 identified 55 SCNAs ( Fig. 1a and b). Of these, 24 (43.6%) were amplifications, and 31 (54.6%) were deletions. In the 172 AA TCGA breast tumors (Table 1), 57 SCNAs were identified ( Fig. 1a and b). Of these, 27 (47.4%) were amplifications, and 30 (52.6%) were deletions. To integrate SCNAs identified separately in EA and AA breast tumors, non-overlapping SCNAs between the two groups were retained as is, while overlapping SCNAs were sub-divided based on the degree of overlap (see Supplementary Methods and Figure S1). This integration resulted in a total of 132 unique SCNAs.
Overall, AA breast tumors had a higher percentage of their genomes affected by SCNAs (termed "SCNA burden"; average 47.5%) in comparison to EAs (39.6%; p = 3*10 − 4 ; Figure S3). To assess specific regional SCNA differences between AA and EA tumors, the degree of alteration at each SCNA was quantified as the area under the copy number ( log 2 copy number 2 ) curve (cnAUC) for each patient (see Methods, Supplementary Methods, and Figure S2). For these analyses, positive race differences in amplified SCNAs (i.e. AA cnAUCgain -EA cnAUCgain > 0) and negative race differences in deletion SCNAs (i.e. AA cnAUCloss -EA cnAUCloss < 0) represent higher magnitude copy number gains and losses in AA relative to EA tumors, respectively. In breast tumors, there were 58 SCNAs (44% of 132) that showed significant differences by race (individual SCNA p < 0.1; Table S1), and this total number of breast cancer race-differentiated SCNAs was greater than what would be expected by chance alone (permutation p < 10 − 4 ; see Supplementary Methods and Figure S4). Of these 58 SCNAs, 25 were amplifications and 33 were deletions. For the majority of both amplifications and deletions, AA breast tumors had higher magnitude SCNAs relative to their EA counterparts. Specifically, of the 25 amplifications, 18 (72%) had positive differences (i.e. more gain in AA tumors), and in the 33 deletions, 22 (66.7%) had negative differences (i.e. more loss in AA tumors). In particular, for the top five most significant race-differentiated SCNAs, all had higher magnitude copy number changes in AAs, and these SCNAs fell into the following two loci: two adjacent deletions on chromosome 5 (5q11-5q14, p = 3.48*10 − 6 ; 5q14-5q22, p = 1.57*10 − 5 ) and three contiguous amplifications on chromosome 10 (10p15, p = 7.24*10 − 5 ; 10p15, p = 9.92*10 − 5 ; 10p15-10p13, p = 1.87*10 − 5 ).
A total of 21 SCNAs (28% of 74; 5 amplifications and 16 deletions) significantly differed by race (p < 0.1, Table S2), and this total number of prostate cancer race-differentiated SCNAs was greater than what would be expected by chance alone (permutation p = 0.0064; see Supplementary Methods and Figure S4). Similar to the breast tumors, AA prostate tumors displayed a greater percentage of higher magnitude SCNAs at both amplifications (80%) and deletions (81%).
To further validate the six race-differentiated SCNAs, we tested for racial differences in the following four additional tumor types with sufficient data for AAs in TCGA: endometrial, ovarian, lung, and kidney tumors. The results of the SCNA differences by race for these four tumor types are presented in Table S3. From the six alterations and consistent with the breast and prostate results, the amplification SCNA on 8q21.3-q24.3 and deletion SCNA on 13q12.3-21.3 were more extreme in AAs relative to EAs in endometrial cancer; this was also the case with the deletion SCNA on 13q12.3-21.3 in ovarian cancer (Table S3).
As AAs are admixed (i.e. have genomes composed of DNA from more than one ancestral population, African and European ancestral populations for AAs), we investigated whether germline chromosomal ancestry (either African or European) differences at each SCNA can explain the greater magnitude of SCNAs observed in admixed AA tumors. A consistent association with local ancestry in these SCNAs would be indicated by: 1) a positive difference in amplification SCNAs reflecting increasing gain associated with increasing regional chromosomal African ancestry or; 2) a negative difference in deletion SCNAs reflecting more loss associated with increasing regional chromosomal African ancestry. By these criteria, the amplification at 8q21.3-8q24.3 and the deletion at 13q12.3-13q21.3 in both breast and prostate tumors showed African ancestry differences consistent in direction with the race association analyses, but these did not reach statistical significance in either tumor type (Table S4).
Unsupervised hierarchical clustering using copy number of genes residing within the six race-differentiated SCNAs identified patient groups with differences in survival To assess the clinical relevance of the six racedifferentiated SCNAs consistent across both breast and prostate cancer, unsupervised hierarchical clustering was first used to classify samples based on the copy number of genes residing within these alteration regions in breast and prostate cancer separately. In breast cancer, three patient groups (labelled as BRG1-3) were identified (Fig. 3a) based on the gap statistic [21] (Figure S9A), and BRG2 and 3 contained tumors with more copy number changes. Consistently, AAs were more likely classified as BRG2 and BRG3 (Table S5, p = 0.01). The three patient groups also exhibited differences in both 10-year overall survival (OS) and 10-year PFS, with BRG2 and 3 group patients having worse OS ( Fig. 3b; hazard ratio [HR] HR BRG3vs1 = 1.69, p = 0.007; HR BRG2vs1 = 1.65, p = 0.389) and PFS (Fig. 3c; HR BRG3vs1 = 1.36, p = 0.075; HR BRG2vs1 = 1.76, p = 0.214) relative to BRG1. In a multivariate model, Fig. 2 Race-differentiated SCNAs identified in TCGA breast and prostate tumors. a Circos plot of all race-differentiated SCNAs (p < 0.1) identified in TCGA breast (outer circle) and prostate (inner circle) tumors. The red color indicates SCNAs with higher magnitude amplification in African Americans (AAs), and the orange color indicates SCNAs with higher magnitude amplification in European Americans (EAs). The dark blue color indicates SCNAs with higher magnitude deletions in AAs, while the light blue shows SCNAs with higher magnitude deletions in EAs. The nine shaded regions indicate alterations that exists in both tumor types. Among those nine, the six SCNAs with consistent differences by race across tumor types are indicated by stars. As an example of one consistent region, average log 2 (copy number/2) profiles on chromosome 5 were plotted for each race in b breast and c prostate tumors. For each tumor type, the top panel corresponds to the complete chromosome 5 profile, and the lower panel corresponds to chromosome 5q15-21.1, where the actual boundaries are from the prostate tumor SCNA adjusting for race, age, pathologic stage, and SCNA burden, BRG3 had poorer survival relative to BRG1 (Table 2; OS HR BRG3vs1 = 1.72, p = 0.01; and PFS HR BRG3vs1 = 1.34, p = 0.117). Interestingly, no significant differences were observed between AA and EA breast tumors by BRG for 10-year OS (race interaction p-value = 0.393) or PFS (race interaction p-value = 0.347; Figure S10).
For prostate cancer, three patient groups were identified ( Figure S9B) and labelled as PRG1-3 (Fig. 4a). AAs appeared more frequently in PRG2 (Table S5; p = 0.069). Due to the limited number of OS events, analyses were restricted to10-year PFS. PRG2 and PRG3 patients demonstrated worse 10-year PFS compared with PRG1 (Fig. 4b, HR PRG2vs1 = 2.14, p-value = 0.111; HR PRG3vs1 = 2.72, p-value = 0.001). In a multivariate model, adjusting for race, age, Gleason grade, and SCNA burden, PRG3 patients had suggestive worse survival relative to PRG1 patients ( Table 2, HR PRG3vs1 = 2.04, p = 0.056). Consistent with the breast cancer results, there was no evidence that the association between PRGs and PFS differed between AAs and EAs (race interaction p-value = 0.977; Figure S11). Differentially expressed genes in the six race-differentiated SCNAs common to breast and prostate cancer To identify race-differences in gene expression consistent with the differences in copy number magnitude within the six race-differentiated SCNAs in both breast and prostate tumors, we performed differential gene expression analyses in both tumor types using TCGA RNA-Seq expression data. Of 522 unique mRNA transcripts in theses six race-differentiated SCNAs, after adjusting for multiple comparisons using a false discovery rate (FDR) of 5%, a total of 18 mRNAs (Table S6) were consistently differentially expressed by race in both tumors. Among these 18 genes, we further assessed whether copy number was associated with gene expression, and the results are summarized in Table S7. We identified 15 genes in breast cancer and 16 genes in prostate cancer where gene expression was significantly (FDR < 0.05) associated with copy number, and 15 of these were shared between breast and prostate cancer.

Discussion
Breast and prostate cancer are two hormonally-driven cancers with established racial disparities in both incidence and severity in women and men, respectively. SCNA events can be considered a primary driver of carcinogenesis across these two cancers, and in our analysis of TCGA breast and prostate cancer copy number alteration data, we demonstrated that a majority of the race-differentiated SCNAs for either tumor type was of greater magnitude in AAs. Further, inter-tumor comparisons demonstrated six SCNAs common to both cancers that had higher magnitude SCNAs in AAs. Two of these race-differentiated SCNA were validated in endometrial and ovarian cancer, two other hormonally-driven cancers. Using unsupervised hierarchical clustering on copy number of genes residing in these six SCNAs, we demonstrated three patient groups in breast cancer and three patient groups in prostate cancer with distinct survival profiles, and these patient groups further differentiated within TNBC based on PFS. Finally, within these six racedifferentiated SCNAs shared by both breast and prostate tumors, 15 mRNA transcripts were differentially expressed by race and copy number in both tumors types, consistent with observed SCNA differences, and these included known cancer genes.
AA women with breast cancer may have greater intratumor genetic heterogeneity and increased basal gene expression even within the triple-negative phenotype, suggesting a more aggressive tumor biology among AA patients [23]. However, a recent analysis of TCGA data on 178 TNBCs found no difference in somatic mutation counts by race [24], which suggests that structural DNA changes may play a greater role in driving racial tumor biology differences. In our current analysis of TCGA breast tumors, 69% of the SCNAs identified were of greater magnitude in AAs. A previous study [17] of 53 AA and 206 EA breast tumors using array comparative genomic hybridization (aCGH) with 4320 probes showed nine race-differentiated SCNAs in triple-negative cases on chromosomes 5q, 8q, 9q, 10q, 14q and 15q, where AAs had greater changes in all nine SCNAs. Four (5q31, 8q21-24, 14q32, and 15q12) of these overlapped with regions identified in our analyses.
In a recent analysis of TCGA breast cancer data by Hou et al. [25], AAs had more loss in four regions (8p23, 13q14, 11q23, and 21q21) and more gains in six regions (8q24, 8p12, 8q11, 19q12, 5p15, and 12q15) based on analysis of nominal copy number deletion and amplification calls. Using our magnitude-based approach (i.e. cnAUC), which additionally adjusted for age at diagnosis and pathological stage, loss on 8p23, 13q14, 11q23 and gains on 8q24, 8q11 and 19q12 were greater in AAs similar to Hou et al. We did not identify racedifferentiated SNCAs in the other four regions, but we did identify additional race-differentiated SCNAs (19 deletions and 11 amplification; Fig. 1a and b and Table S1).
While our breast cancer analyses utilized a similar TCGA tumor dataset as Huo et al. [25], there are distinctions between the approaches that may have led to differences in our findings beyond the slight difference in the number of breast tumors available at the time of analysis. As mentioned above, Huo et al. discretized SCNAs into nominal gains and losses, and their analysis compared the differences in the frequency of SCNAs between the two races. In comparison, we used a novel quantitative measure across each SCNA for each tumor, which not only captures differences in SCNA frequency but also incorporates the alteration magnitude. As the use of this additional information identified more race-differentiated SCNAs relative to Huo et al. [25], our approach may be more sensitive to the detection of such differences. Another important distinction is subject inclusion. While both studies utilized the germline genetic data available to estimate genetic ancestry and the agreement between ancestry and self-reported race to select subjects, we applied a more stringent threshold for study entry. For example, Huo et al. included 86 subjects with missing self-reported race data by assigning them an ethnicity most consistent with their genetic ancestry. In contrast, we chose to only include subjects whose selfreported race agreed with their genetic ancestry. Our approach is likely more conservative, which again may impact differences in findings. Finally, Hou et al. did not adjust for pathologic stage in their analysis, which is correlated with the SCNAs. Not accounting for pathologic stage may confound the relationship between SCNAs and race, potentially producing a false positive effect that is simply due to differences in pathologic stage between AAs and EAs.
In addition to breast cancer, growing evidence also suggests that prostate tumors of AA men undergo somatic changes that are differ from EA men. For example, Khani et al. [26] showed that tumors of AA men have less ERG rearrangement and PTEN deletion but more SPINK1 overexpression than clinically similar tumors of EA men. In prostate tumors, we identified five amplification and 16 deletion SCNAs that differed by race. Among these SCNAs, loss at 6q13-22, 13q13-14, and 16q11-24, and gains of 8q24 also showed higher frequency in 20 AA prostate tumors from a previous genome-wide SCNA scan [27]. In another study [16] where the SCNAs of 28 AA and 180 EA prostate tumors were profiled by aCGH, 23 regions were identified as race-differentiated, and in 21 (91.3%) regions, AA tumors had more copy number changes. The greater frequency of deletions in AA tumors observed on chromosomes 1q31 and 16q22 were consistent with our findings. Combined with the literature on race-differentiated SCNAs in breast tumors, these studies collectively suggest that the genomes of AA breast and prostate tumors are less stable and harbor greater numbers of SCNAs compared with their EA counterparts.
While the mechanism(s) giving rise to this difference in stability are unknown, one possible explanation is differential expression of common fragile sites (CFS). CFS expression is known to give rise to SCNAs, including deletions that lead to tumor suppressor gene deactivation in multiple cancer types [26], and currently, there are 125 known CFS spread throughout the human genome. Using CFS genomic coordinates provided in the HumCFS database [28], we found that four of the six race-differentiated SCNAs common to both breast and prostate cancer overlapped with CFS (5q11.2-5q14.1 with FRA5H, 5q15-5q21.1 with FRA5B and FRA5F, 8q21.3-8q24.3 with FRA8D, 13q12.3-13q21.3 with FRA13B and FRA13C), including those that were also observed in endometrial (8q21.3-q24.3 and 13q12.3-q21.3) and ovarian (13q12.3-q21.3) cancers. This observation raises the provocative question of whether CFS expression in tumors differs by race of the individual. While evaluating CFS expression as a possible mechanism underlying race-differentiated SCNAs is outside of the scope of this paper, our observation suggests that it is worthy of further investigation.
Using the gene-level copy number information within the six race-differentiated SCNAs that have consistent directional differences by race across both tumor types, we identified three patient groups in both breast and prostate tumors through unsupervised hierarchical clustering. Groups containing samples with more frequent copy number changes in these six regions are associated with poorer OS and PFS in breast cancer and worse PFS in prostate cancer. These differences remain after adjusting for self-reported race, age-at-diagnosis, tumor pathology, and genome-wide SCNA burden, highlighting the clinical significance of identified SCNA patient groups. Poorer breast and prostate cancer outcomes in AAs compared to EAs is at least partly explained by the higher incidence rates of biologically aggressive disease in AAs. In breast cancer, AAs had a two-fold higher incidence rates of TNBC compared with EAs [29], and evidence suggests that prostate cancer progresses more rapidly for AA men, who experience a higher prostate cancer volume at diagnosis and a four-times higher incidence of metastatic disease relative to EA men [30]. A previous analysis of TCGA data reveled that AA TNBC patients had a shorter time to disease progression compared with EAs, but no disparity with hormone receptor-positive or HER2/neu-positive patients [23]. Other studies have demonstrated similar rates of disease progression for AA and EA TNBC patients but poorer survival for AA cases with hormone receptor positive tumors [23,31,32]. Our analyses examining the effect of race within each identified SCNA group demonstrated a lack of statistically significant differences in survival between AAs and EAs. This observation favors the hypothesis that breast and prostate cancer outcome disparity between AA and EA is due to a higher frequency of deleterious SCNAs in AA relative to EA tumors. Furthermore, we demonstrated that in TNBC and hormone receptor-positive breast cancer, survival differences still exists between SCNA-defined groups, suggesting that these patient groups may provide additional prognostic information to currently defined tumor subtypes. Independent studies will be needed to confirm the association between the SCNA patient groups we have identified and survival.
Within the six SCNAs that were found to have consistent differences in magnitude by race across breast and prostate tumors, we identified 15 genes that were differentially expressed by race and copy number, many of which have known cancer related functions. For example, TOP1MT on 8q24.3, a mitochondrial DNA topoisomerase, was the gene with the most significant difference in expression by race in both tumor types. TOP1MT has c-myc binding sites in its promoter and follows the expression of the master regulator oncogene, c-myc [33]. The well-characterized tumor suppressor gene RB1 (retinoblastoma susceptibility gene) on 13q14.2 displayed significantly lower expression in AA breast and prostate tumors, which correlates with the greater loss of DNA we observed in AAs relative to EAs. RB1 was the first molecularly defined tumor suppressor by suppressing cell cycle progression and controlling chromatin remodeling and cell death [34,35]. RB1 loss is associated with tumor progression in both breast and prostate cancers [36][37][38]. The long non-coding RNA PVT1 within the amplification SCNA identified on 8q24.21 is overexpressed in both AA breast and prostate tumors, and it was found that the overexpression of the MYC oncogene in tumors depends on PVT1 expression [39]. In summary, these findings suggest that genes known to be involved in tumorigenesis and cancer progression, are similarly differentially expressed by race and copy number and may, therefore, underlie racial differences in tumor biology linked to the six race-differentiated SCNAs between AA and EA breast and prostate tumors.
There are multiple limitations in this study that should be mentioned. First and foremost is the underrepresentation of AAs (and other minorities) in TCGA relative to EAs. An important consequence of this reality is the limitation in statistical power to detect underlying racial differences in SCNAs both within and across tumor types. While we were able to identify statistically significant aggregate racial differences in SCNAs using a permutation approach, there are likely additional, potentially more subtle racial differences in SCNAs that were missed in this analysis. Further, our ability to explore how these differences were shared across different tumor types was similarly curtailed by this limitation. Also, while we were able to adjust our analyses for important variables that could confound both the relationship between race and SCNAs and the relationship between SCNAs and survival, this list of potential confounders was limited in TCGA. Together, coupled with the positive findings from this study, these limitations suggest that larger genomic studies with racial balance and comprehensive and consistent clinical/pathologic data collection are necessary to make further progress towards understanding how race influences cancer biology and outcomes.

Conclusions
In summary, using TCGA data, we demonstrated that a majority of race-differentiated SCNAs in breast and prostate tumors have greater magnitude alterations in AAs relative to EAs and that this was the case for all of six race-differentiated SCNAs that were shared consistently across these tumor types. Combined, these findings suggest that breast and prostate tumor genomes of AAs may be more susceptible to SCNAs that underlie a more aggressive phenotype. This information may be useful for deciding whether adjuvant treatment is needed after surgery in the subset of patients with SCNAs associated with worse outcomes. Finally, while this study focused on breast and prostate tumors, our additional findings from endometrial and ovarian tumors that are consistent with a subset of the racially differentiated SCNAs in breast and prostate cancers suggest that this genomic phenomenon may be representative of other tumor types with racially disparate outcomes, which warrants further studies with appropriate race representation.