Research ArticleCancer

Sex differences in GBM revealed by analysis of patient imaging, transcriptome, and survival data

See allHide authors and affiliations

Science Translational Medicine  02 Jan 2019:
Vol. 11, Issue 473, eaao5253
DOI: 10.1126/scitranslmed.aao5253

Brain tumors’ battle of the sexes

In recent years, there has been increasing interest in personalized medicine for cancer, considering the unique biology of each tumor and patient to optimize therapeutic approaches. Sex differences play a role in patient outcomes, and Yang et al. determined that in the case of the brain tumor glioblastoma, these go beyond hormonal influences and appear to be intrinsic to the tumor cells themselves. The authors found that the sex of the patient correlates not only with prognosis but also with responses to different treatments, suggesting that it may be an important factor to consider when optimizing the therapeutic regimen for each patient.

Abstract

Sex differences in the incidence and outcome of human disease are broadly recognized but, in most cases, not sufficiently understood to enable sex-specific approaches to treatment. Glioblastoma (GBM), the most common malignant brain tumor, provides a case in point. Despite well-established differences in incidence and emerging indications of differences in outcome, there are few insights that distinguish male and female GBM at the molecular level or allow specific targeting of these biological differences. Here, using a quantitative imaging–based measure of response, we found that standard therapy is more effective in female compared with male patients with GBM. We then applied a computational algorithm to linked GBM transcriptome and outcome data and identified sex-specific molecular subtypes of GBM in which cell cycle and integrin signaling are the critical determinants of survival for male and female patients, respectively. The clinical relevance of cell cycle and integrin signaling pathway signatures was further established through correlations between gene expression and in vitro chemotherapy sensitivity in a panel of male and female patient-derived GBM cell lines. Together, these results suggest that greater precision in GBM molecular subtyping can be achieved through sex-specific analyses and that improved outcomes for all patients might be accomplished by tailoring treatment to sex differences in molecular mechanisms.

INTRODUCTION

Current epidemiological data indicate that sex differences exist in the incidence of cardiovascular disease, disorders of the immune system, depression, addiction, asthma, and cancers (14), including glioblastoma (GBM) (5). Although sex differences in disease incidence and severity may parallel variation in circulating sex hormone concentrations, in many cases, sex differences exist across all stages of life, indicating some independence from acute hormone action (3, 6). Sex differences in GBM are evident in all age groups and therefore cannot be solely the consequence of activational effects of sex hormones (5, 711). Enumerating the molecular bases for sex differences in GBM is likely to reveal fundamental modulators of cancer risk and outcome as well as guide sex-specific components of precision medicine approaches to cancer treatment.

Identifying the basis for sex differences in cancer biology cannot be accomplished by analysis of merged male and female datasets. Instead, it requires comparison of results from parallel analyses of male and female data. The importance of this was recently highlighted in a study of asthma, a disease driven by both genetic and environmental factors, which occurs in twice as many boys as girls. Mersha et al. (2) examined the influence of genetic variants on asthma, including an analysis of shared and sex-specific variant effects. Of 47 variants that correlated with asthma risk in the sex-specific analyses, only 21 were detected in the combined analysis, suggesting that biologically important mechanisms of disease were obscured by a “net canceling effect” that arose from opposing effects of genetic variation in sexes. A similar effect was observed in neurofibromatosis 1 (NF1)–associated low-grade gliomas. Despite equal tumor incidence in males and females, polymorphisms in AC8 in patients with NF1 increased the risk of low-grade glioma in female patients but reduced the risk in male patients (12). The effect of AC8 polymorphisms, which may be related to the mechanistic role of cAMP (cyclic adenosine 3′,5′-monophosphate) in NF1-assoiated glioma (13, 14), was unapparent without a sex-specific analysis because of the net canceling effect.

Whereas low-grade glioma incidence is nearly identical in males and females, malignant brain tumors in general occur more commonly in males, regardless of patient age or geographical location (5, 11, 15, 16). As shown in recent reports, GBM occurs with a male-to-female ratio of 1.6:1 (5, 810). In particular, although the understanding of molecular subtypes of GBM is still evolving (17), three of the four originally described transcriptional subtypes of GBM—mesenchymal, proneural, and neural GBM—exhibit a 2:1 male-to-female incidence ratio, whereas classical GBM occurs with equal incidence (15, 18). To date, analyses of the transcriptome data from which these molecular subtypes were derived have been performed with merged male-female data and have not yielded insights into the molecular basis for sex differences in GBM incidence.

In addition to sex differences in incidence, emerging analyses suggest that patient outcomes may also differ between males and females in the pediatric (19) and adult patient populations with GBM (20). In a study analyzing more than 27,000 patients, Trifiletti et al. (10) found that female sex was associated with longer survival, as did Ostrom et al. (20) in an analysis of 5372 GBM cases from the National Cancer Institute’s Surveillance, Epidemiology, and End Results (SEER) program and an additional 228 GBM cases from the Ohio Brain Tumor Study. Similarly, female patients exhibited longer survival from gliosarcoma (8), and being female was associated with better outcome in a nomogram for predicting survival in patients with GBM (9). Thus, the elucidation of sex-specific mechanisms in GBM has the potential to improve outcome for all patients by refining our understanding of disease causation and treatment response.

Here, we performed quantitative analyses of therapeutic responses in male and female patients with GBM using a validated magnetic resonance imaging (MRI)–based method for calculating tumor growth velocities. We also applied a computational algorithm to male and female GBM transcriptome data to gain insights into the relevance and biological basis of sex differences in GBM. Our studies indicate that standard treatment is more effective for females than for males with GBM and that, for the current standard of care—surgery, radiation, and temozolomide (TMZ) chemotherapy—survival in males is correlated with the expression of cell cycle regulators, whereas in females it is correlated with the expression of integrin signaling pathway components. These studies provide a coherent view of sex differences in GBM biology and their clinical ramifications. They support the development of diagnostics and treatments that incorporate sex differences in GBM biology.

RESULTS

Standard treatment is more effective in female compared with male patients with GBM

Sex differences in GBM incidence have been repeatedly reported (5, 711), and recent studies have suggested that being female is associated with better outcome from GBM in both adults and children (810, 19, 20). The introduction of TMZ as a component of trimodal care for adults with GBM has improved outcomes and highlighted factors, such as O6-methylguanine DNA methyltransferase (MGMT) promoter methylation, that affect response and survival (21, 22). Thus, we wondered whether sex differences in GBM survival are a consequence of differential treatment effects on male versus female patients. To answer this question, we used an MRI data–based analysis, with which the velocity of radial tumor expansion can be determined (2326). Growth velocity, which correlates with outcome (27, 28), was measured approximately every 2 months in a cohort of 63 patients with GBM (40 males and 23 females) treated with standard-of-care surgery, focal irradiation (XRT), and systemic TMZ chemotherapy (fig. S1) (21, 22). Analysis of serial magnetic resonance (MR) images obtained during postradiation maintenance TMZ treatment indicated that female patients exhibited a greater response to treatment than male patients. Although the initial tumor growth velocities were similarly distributed in male and female patients (P = 0.3985, Wilcoxon rank sum test; Fig. 1A), a steady and significant decline in growth velocity during TMZ treatment was evident only for female patients [P = 0.02569 (female trend test), P = 0.1186 (male trend test); Fig. 1A]. To determine whether the initial TMZ velocity correlated equivalently with survival across sexes, we fit a Cox proportional hazard model with the main effect of velocity (in continuous scale) for the male and female populations separately and found that the velocity had a sex-specific impact on survival (male, P = 0.302; female, P = 0.0161). To visualize the sex-specific effect, we performed an iterative Kaplan-Meier (KM) analysis in male and female GBMs separately to divide male and female patients into high versus low velocity groups and tested the survival difference using the log rank test. For female patients, lower first TMZ velocity was associated with a significantly longer survival compared with higher velocity [median survival, 3090 days versus 681 days (P = 0.00817); Fig. 1B]. In contrast, male patients exhibited no statistically significant correlation between survival and velocity [median survival, 1111 days versus 533 days (P = 0.263); Fig. 1B]. To determine whether first TMZ velocities correlated with Revised Assessment in Neuro-Oncology (RANO) clinical response criteria, we compared survival, growth velocity, and RANO measures for this patient cohort. No significant correlations were detected (fig. S2). Although these data cannot distinguish between the therapeutic effects of radiation versus TMZ, they do suggest that females with GBM may benefit more from standard treatment than males with GBM and that this difference in response, which is detectable using tumor growth velocity measures, may contribute to their survival advantage.

Fig. 1 Sex differences in MRI-based metrics of therapeutic responses and their correlation with survival.

(A) Tumor growth velocities calculated from serial T1 gadolinium (Gd)–enhanced MR images exhibit progressive decline for female (n = 23) but not for male (n = 40) patients with GBM. (B) Velocity of tumor growth (low velocity, green line; high velocity, purple line) over the first TMZ imaging interval (1 to 3, 28-day cycles of TMZ) stratifies female survival (P = 0.00817, log rank) but not male survival (P = 0.263). (C) Histograms of pretreatment D and rho values in all available MRI cases (independent 53 and 318 GBM case series) for male (n = 227) and female (n = 144) patients. (D) Pretreatment D significantly stratifies survival among females (n = 144; P = 0.0071, log rank) and not among males (n = 227, P = 0.61). (E) High pretreatment rho is associated with worse survival outcomes for both females (n = 144; P = 0.032, log rank) and males (n = 227, P = 0.0037).

To validate and further investigate the basis for this difference in response, we applied an established mathematical model of glioma proliferation and invasion (25, 26, 29). We examined the presurgical MRIs (T1 gadolinium and T2 sequences) of 53 patients from the original growth velocity cohort combined with an additional independent cohort of 318 patients, for a total of 371 newly diagnosed patients with GBM (227 males and 144 females). We found that the distribution of estimated net infiltration rates (D, in square millimeters per year) and net proliferation rates [rho (ρ), 1/year] did not differ between males and females before surgery (Fig. 1C). Next, we sought to determine which of the two components (24, 26, 30) was predictive of survival in male and female patients, separately. The sex-specific median of each component was used to unbiasedly dichotomize the patients into low and high groups. Female patients with low D (≤23.03 mm2/year) had significantly longer survival [median overall survival (OS), 589 days versus 390 days (P = 0.0071)] compared with female patients with high D (>23.03 mm2/year; Fig. 1D, left). This was in contrast to male patients for whom survival did not differ [median OS, 540 days versus 450 days (P = 0.614)] as a function of D (≤28.993 mm2/year versus >28.993 mm2/year; Fig. 1D, right). Rho, in contrast to D, stratified survival for both males and females. Females with low rho (≤18.25 per year) exhibited median OS of 542 days as compared with those with high rho (>18.25 per year), who exhibited median OS of 415 days (P = 0.032; Fig. 1E, left). Males with low rho (≤18.25 per year) exhibited median OS of 596 days as compared with those with high rho (>18.25 per year), who exhibited median OS of 410 days (P = 0.0037; Fig. 1E, right). In the independent analysis of the expansion cohort of 318 patients (195 males and 123 females), D and rho had effects on survival in female patients that were similar to those described in the discovery cohort (227 males and 144 females), but neither D nor rho stratified survival for male patients (fig. S3). Together with the established sex differences in incidence, these data suggest that the biology of male and female GBM may be distinct and that outcomes for all patients might be improved if therapies were better tailored to patient sex.

Sex differences in GBM biology are revealed by JIVE decomposition

To gain insight into potential sex differences in GBM biology, we examined the transcriptome data available through The Cancer Genome Atlas (TCGA). Using Joint and Individual Variation Explained (JIVE) to integratively decompose the male and the female transcriptome data of the TCGA dataset into three orthogonal components, we identified the joint structure that was common to both sexes, the individual structure that was specific to each sex, and, additionally, the residuals (fig. S4). The heat maps of the male joint structure across the male patients with GBM and the female joint structure across the female patients with GBM indicated that the joint structures extracted by JIVE closely captured the dominant molecular signatures defining the TCGA GBM subtypes (fig. S5 and table S1). However, the joint component only explained ~45% of the total variance in the transcriptomes for each sex, whereas the sex-specific components, independent of the joint components, explained a large proportion of the remaining variability. Specifically, the male-specific component accounted for 38.5% of the total variability in the male transcriptome, and the female-specific component explained 33.6% of the total variability in the female transcriptome (fig. S6). The extracted male and female individual components exhibited distinct patterns compared with their counterpart joint structure, and more importantly, the male-specific component showed distinct patterns compared with the female-specific component (Fig. 2, A and B). We hypothesized that focused analyses of the extracted sex-specific components would reveal which gliomagenic mechanisms are most characteristic of male versus female GBM.

Fig. 2 Heat maps of joint and sex-specific expression components of TCGA GBM transcriptome data revealed by JIVE.

The heat maps visualize each expression component. Each row represents a gene, and each column a patient sample. For each patient, there are two color codes presented above the heat map. These identify their assignment to sex-specific clusters and to TCGA molecular subtypes (gray indicates unassigned samples). Samples were ordered by sex-specific clusters. The original female (A) and male (B) expression data were decomposed into the shared expression component common to both sexes (“Joint”) and the expression component individual to each sex (“Female specific” and “Male specific”) and residuals as indicated. The female-relevant heat maps (A) show 283 signature genes that define the five female-specific clusters, and the male-relevant heat maps (B) show 293 signature genes that define the five male-specific clusters. (C) The Venn diagram of male and female signature genes indicates that 116 genes are in common.

Sex-specific clusters are identified using the TCGA sex-specific transcriptome expression

To identify sex-specific patient subgroups, we performed independent hierarchical clustering on the male- and female-specific components from the JIVE decomposition. Weighted and unweighted consensus clustering was applied to the sex-specific expression to evaluate the robustness of sex-specific clustering (fig. S7, A and B). To determine the optimal number of sex-specific clusters, we varied the total number of clusters from two to six for each sex (fig. S7, C and D) and examined the cumulative density function (CDF) curves for the consensus matrices (fig. S7, E and F). We compared the resultant increase in area under the CDF curves (fig. S7, G and H) when the total number of clusters increased by one. The similarity among samples from each sex-specific cluster was examined to remove samples with great dissimilarity to the majority of samples in the cluster based on the Silhouette scores [fig. S7, A and B, right panels]. Five male (mc1–5) and five female (fc1–5) clusters were thus identified as optimally capturing the transcriptomic subtypes within male and female TCGA data. The five male and five female clusters were defined by sets of 293 and 283 genes (Fig. 2C), respectively, with 116 in common but 177 unique to the male clusters and 167 unique to the female clusters (table S2).

Cases from multiple TCGA molecular subtypes (18) were distributed to each of the five male or five female clusters (Fig. 2), indicating successful separation of the individual components from the joint structure components and increasing the likelihood that this approach could reveal sex effects on gliomagenic mechanisms. The one exception was fc3, 70% of which were proneural subtype tumors with IDH1 mutations (seven IDH1 mutants and three wild types (WTs)]. In contrast, male proneural subtype tumors with IDH1 mutations were distributed across three of the male clusters, suggesting that IDH1 mutations may have sex-specific effects in GBM.

Sex-specific clusters are robust to excluding IDH1 mutant cases

Current diagnostic criteria indicate that IDH1 mutant and WT GBM are two separate diseases (31). Thus, we examined whether the definition of the sex-specific clusters was robust to excluding IDH1 mutant cases. We removed the IDH1 mutant and glioma CpG island methylator phenotype (G-CIMP) cases from the TCGA, GSE13041, GSE16011, and Rembrandt datasets. We followed the same procedure (independent JIVE analysis, consensus clustering, and determination of optimal total number of sex-specific clusters) to identify sex-specific clusters in IDH1 WT cases. Most of the samples (65 to 96.2%) were in agreement with their cluster assignments from the initial analysis, and mc5 was rediscovered in the IDH1 WT cases (Fig. 3). Because fc3 was predominantly composed of IDH1 mutant cases, it was substantially diminished in this analysis.

Fig. 3 Sex-specific survival effects of IDH mutation.

(A) OS benefit of fc3 and mc5 is demonstrated in the combined TCGA, GSE13041, GSE16011, and REMBRANDT datasets. See table S3 for P values and hazard ratios (HRs). (B) OS for IDH1 WT cases indicates that both fc3 and mc5 exert effects on survival in the absence of IDH1 mutation. (C) OS in IDH1 mutant cases indicates that male-specific clusters are still associated with an effect on survival. The numbers of female IDH1 mutant cases not assigned to fc3 are n = 3, 2, 3, and 6 in fc1, fc2, fc4, and fc5, respectively, using TCGA and GSE16011 samples in combination (see table S7). (D) IDH1 mutation confers a similar survival benefit in males and females with GBM. (E) The survival benefit of fc3 is independent of IDH1 status. In contrast, IDH1 status exerts a significant effect on survival in mc5 cases. P = 4.3 × 10−4 for the comparison between mc5 cases with and without IDH1 mutation. Overall log rank test P value is shown comparing across all the groups presented in each panel (table S8 shows the P values and HRs for all pairwise comparisons).

Survival differences exist among sex-specific clusters in the TCGA GBM cohort

To establish the importance of the sex-specific clusters, we next determined whether the sex-specific clusters in the TCGA data were associated with differences in survival outcomes. KM analyses of male and female clusters confirmed that survival differences exist among both male and female clusters (Fig. 4, A and B). Not surprisingly, fc3, in which 70% of the cases are IDH1 mutant, exhibited significantly better disease-free survival (DFS) with a median time to progression (TTP) of 1758 days compared with each of the other four female clusters [fc1, 259 days (P = 3.3 × 10−5); fc2, 289 days (P = 5 × 10−4); fc4, 182 days (P = 1.64 × 10−4; fc5, 350 days (P = 9.6 × 10−5); Fig. 4A and table S3]. In contrast, although IDH1 mutant cases segregated nearly equally to mc2, mc3, and mc5, only mc3 (median TTP, 408 days) and mc5 (median TTP, 262 days) were associated with prolonged DFS compared with other male clusters [mc1, 240 days (P = 1.2 × 10−2 versus mc3); mc2, 186 days (P = 7.1 × 10−3 versus mc3, P = 2.8 × 10−2 versus mc5); mc4, 158 days (P = 7.3 × 10−3 versus mc3, P = 1.6 × 10−2 versus mc5; Fig. 4B and table S3]. This finding suggested that an interaction may exist between IDH1 mutation and sex-specific cluster features in males but not in females in the determination of survival. To further evaluate this possibility, we separated the IDH1 mutant patients from the sex-specific clusters (Fig. 4, C and D). Only three cases in fc3 were IDH WT, and each of them was alive at 5 years [median DFS for fc3 was not calculable; fc1, 256 days (P = 1.4 × 10−2); fc2, 274 days (P = 1.7 × 10−2); fc4, 182 days (P = 3.4 × 10−2); fc5, 350 days (P = 1.2 × 10−2); Fig. 4C and table S4]. Similarly, the DFS benefit of mc3 and mc5 remained intact after removal of the IDH1 mutant cases [mc3, 408 days; mc5, 262 days; mc1, 240 days (P = 4.2 × 10−3 versus mc3); mc2, 176 days (P = 1.5 × 10−3 versus mc3, P = 1.5 × 10−2 versus mc5); mc4, 158 days (P = 3.1 × 10−3 versus mc3, P = 2.1 × 10−2 versus mc5); Fig. 4D and table S5]. These results suggest that the survival effects of fc3, mc3, and mc5 may be independent of IDH1 mutational status.

Fig. 4 DFS of sex-specific clusters in TCGA GBM dataset and OS of sex-specific clusters in three independent datasets combined.

(A) DFS in TCGA-derived female clusters (fc1–5). (B) DFS in TCGA-derived male clusters (mc1–5). (C) DFS in TCGA-derived female clusters (fc1–5) in which IDH1 mutant cases are plotted as an independent cluster. (D) DFS in TCGA-derived male clusters (mc1–5) in which IDH1 mutant cases are plotted as an independent cluster. Independent samples combining the GSE13041, GSE16011, and REMBRANDT datasets were assigned to sex-specific clusters, and the superiority of OS of fc3 (E) and mc5 (F) was validated in the independent samples. Overall log rank test P value is shown comparing across all the groups presented in each panel (see tables S3 and S4 for the P values and HRs for all pairwise comparisons).

Survival patterns of sex-specific clusters were independently validated

The transcriptome data of GSE13041, GSE16011, and REMBRANDT were decomposed with the JIVE principal components from the TCGA data analysis (fig. S4), and the independent samples were assigned to the TCGA-derived sex-specific clusters on the basis of the nearest-neighbor algorithm. We then sought to validate the male and female cluster-specific survival profiles using all the independent samples. We were limited to an analysis of OS by data availability. The OS benefit of fc3 and mc5 was validated in these datasets (Fig. 4, E and F, and table S3). Using all the samples of the datasets under analyses (TCGA, GSE13041, GSE16011, and REMBRANDT), median OS for fc3 was 1172 days, as compared with 416 days for fc1 (P = 5.6 × 10−5), 378 days for fc2 (P = 1.2 × 10−7), 423 days for fc4 (P = 8.3 × 10−8), and 359 days for fc5 (P = 4.2 × 10−7) (Fig. 3A, left, and table S3). Median OS for mc5 was 620 days, as compared with 422 days for mc1 (P = 1.4 × 10−6), 360 days for mc2 (P = 8.8 × 10−9), 398 days for mc3 (P = 7.9 × 10−5), and 387 days for mc4 (P = 6.2 × 10−6) (Fig. 3A, right). Of the three validation datasets, only GSE16011 specified IDH mutational status. In this dataset, IDH1 mutant tumors were disproportionately distributed to fc3 but more broadly to multiple male clusters (tables S6 and S7), similar to IDH1 distribution in the TCGA samples.

IDH1 mutation status interacts with sex-specific clusters

IDH1 mutation confers a better prognosis in GBM (32). The survival advantage of mc5 and fc3 was observed, irrespective of IDH1 status, and for males, IDH1 mutation distributed more equally across the clusters without consistent survival benefits (Fig. 3, B and C). In IDH1 WT cases from the combined data (TCGA, GSE16011, and REMBRANDT), mc5 was still correlated with the longest survival among the male clusters (HRs, 0.61 to 0.65 and P = 0.0039 to 0.022; Fig. 3B and table S5), and fc3 had HRs ranging from 0.26 to 0.29 compared with the other female clusters (P = 0.0032 to 0.0093; Fig. 3B and table S4). In IDH1 mutant cases from the combined data (TCGA and GSE16011), the sample size was too small and lacked sufficient power to render statistical significance on survival comparisons between mc5 and fc3 versus all the other male (mc1–4) and female (fc1, fc2, fc4, and fc5) clusters, respectively, but the estimated HRs for mc5 compared with the other male clusters, and fc3 compared with the other female clusters were always below 1 [for mc5: HRs, 0.15 to 0.54, P = 0.017 versus mc1, P = 6.3 × 10−5 versus mc4, P = 0.12 versus mc2, and P = 0.16 versus mc3; for fc3: HRs, 0.19 to 0.91, P = 0.013 versus fc2; tables S6 and S7]. Thus, IDH1 mutation was validated as a good prognostic feature for both males and females (Fig. 3D). However, IDH1 mutation interacts with fc3 and mc5 cluster features differently (interaction P = 0.07; Fig. 3E and table S8). Fc3 conferred longer survival regardless of the IDH1 mutation status. In contrast, IDH1 mutation further stratified survival differences among mc5 cases (Fig. 3E), such that IDH1 mutant mc5 GBM showed comparable or even slightly better survival than IDH1 mutant fc3 GBM [HR = 0.79, confidence interval = 0.34 to 1.9 (P = 0.59); table S8], although statistically not significant. Together with the broader distribution of IDH1 mutation cases across all male sex–specific clusters, these findings indicate that IDH1 mutation interacts with sex in the determination of survival.

Sex-specific clusters show differing survival patterns by TCGA molecular subtype

To gain further insights into cluster-specific effects on survival, we compared the survival differences of the male- and female-specific clusters when stratified by the original Verhaak subtypes (18). We found a consistent cluster effect in which neural, mesenchymal, and proneural specimens in mc5 and fc3 exhibited better survival than tumors of these same Verhaak subtypes that had clustered to mc1–4 or fc1, fc2, fc4, and fc5 (fig. S8). Neither male nor female cluster effects were evident for the classical subtype tumors, the only subtype for which there is no sex difference in incidence (15). These data suggest that for those molecular subtypes of GBM in which sex affects tumor incidence, sex also affects patient survival. In addition, these findings indicate that sex can modulate the impact of specific gliomagenic mechanisms on survival but that not all mechanisms, such as those underlying classical subtype tumors, will be sensitive to the effects of sex.

Pathway analysis indicates that survival in males and females with GBM may be dependent on different mechanisms

The unequal effects of sex on survival for tumors of different molecular subtypes suggest that the effects of sex are not mediated solely by factors such as sex hormones, whose actions would distribute equivalently across patients of a given sex regardless of their molecular subtype. Instead, these findings indicate that either tumor cell–intrinsic sex differences or an interaction between tumor cell–intrinsic and microenvironmental sex differences determines responsiveness to treatment and patient survival. To gain insight into possible mechanisms underlying sex-specific survival benefits, we compared the survival and transcriptome expression of fc3 and mc5.

Median survival for fc3 was 1172 days compared with 620 days for mc5 (Fig. 5A). To test whether similar or distinct mechanisms accounted for these sex differences in survival, we asked what distinguished fc3 and mc5 from the other female and male clusters, respectively. One hundred ninety-seven transcripts distinguished mc5 from the other male clusters, and 123 transcripts distinguished fc3 from the other female clusters (table S2). Using the Genomatix Suite for pathway analysis, we found that 13 transcripts belonging to calcium/calmodulin signaling, synaptic, and other neuronal function pathways were shared between mc5 and fc3 (Fig. 5B). Examination of the female-specific transcripts revealed the integrin signaling pathway as the most significant pathway (adjusted P < 0.001) that distinguished fc3 from other female clusters (Fig. 5C and table S9), with nine transcripts in the pathway (labeled fc3.9). Six of the nine transcripts from this pathway [PLAT (33), CHL1 (34, 35), FERMT1 (36), PCDH8 (37), IGFBP2 (38, 39), and POSTN (40)] have known roles in glioma, and three (PLAT, IGFBP2, and POSTN) can distinguish proneural from classical high-grade glioma gene signatures (41). Six of the nine genes (AK5, AMIGO2, PLAT, CHL1, PCDH8, and IGFBP2) were down-regulated in fc3 compared with other female clusters, suggesting that better survival in fc3 patients is favored by tumors with reduced integrin signaling (fig. S9).

Fig. 5 Analysis of genes and pathways that mediate better survival.

(A) In the combined dataset, the survival of females assigned to fc3 (median survival, 1172 days) was compared with the survival of males assigned to mc5 (median survival, 620 days). (B to D) Genes that distinguished fc3 and mc5 from other female and male clusters, respectively, were compared (see table S2). Pathways in all analyses were prioritized by the combination of the numbers of genes from the pathway involved and the corrected P value for the relevance of the pathway. (B) Calcium/calmodulin signaling was the most significantly involved shared pathway between fc3 and mc5 (adjusted P < 0.001). (C) The integrin signaling pathway was the most significant female-specific pathway (adjusted P < 0.001; table S9). Genes that were up- and down-regulated in fc3 compared with the other female clusters are in red and blue boxes, respectively. (D) Cell cycle regulation was the most significant male-specific pathway (adjusted P < 0.001; table S9). Genes that were up- and down-regulated in mc5 compared with the other male clusters are in red and blue boxes, respectively. See table S2 for the complete gene lists and statistics for each analysis.

Better outcome in mc5 was significantly (adjusted P < 0.001) associated with cell cycle regulation pathways (Fig. 5D). Seventeen transcripts (labeled mc5.17) were components of this pathway, and they included known critical regulators of mitosis such as CDC20 (37, 38), CKS2 (39), PRC1 (40), NUSAP1 (41), PBK (42), cyclin B1 and B2 (43), and KIF20A (44). Fifteen of the 17 transcripts were significantly down-regulated in mc5 compared with the other male clusters (P ≤ 0.0061 for differences in the original expression data; P ≤ 1.7 × 10−6 for difference in male-specific expression data) and approached the expression observed in fc3 (Fig. 6, A and B, and fig. S10, A, A’, B, and B’). NEFH and NEFM were the exceptions, with each exhibiting greater expression in mc5 compared with each of the other male clusters. This suggests that treatment response and survival in males are determined by lower activity in factors that promote cell cycle progression.

Fig. 6 mc5-defining genes and OS in the merged TCGA, GSE16011, and GSE13041 datasets.

(A) Density plots for sex-specific expression of male (in blue) and female (in red) GBM specimens of three mc5-defining genes (BIRC5, KIF20A, and CCNB2). The overlay in male and female plots indicates near-identical expression in the populations. (B) Expression of each gene by sex and sex-specific clusters is presented as boxplots. (C) High and low expression groups for each gene were defined relative to the level of expression that distinguished mc5 from the other male clusters (see the “Overall sex-specific survival effects” section in the Supplementary Materials). The survival effects of differences in expression were determined for males and females. Each gene exerted a greater effect on survival in males compared with females. P values from the Cox regression model are labeled in red for comparisons between survival curves of female patients with GBM with low versus high expression of each and labeled in blue for the same survival analysis of male patients with GBM. The P value labeled in green refers to the interaction of sex and expression of a gene in the Cox regression models. Parallel analyses of the fc3-defining genes and the other mc5-defining genes are presented in figs. S9 and S10, respectively.

Each of the 9 genes that distinguished fc3 (fc3.9) and each of the 17 genes that distinguished mc5 (mc5.17) from other female and male clusters were similarly expressed in male and female patients with GBM overall (Fig. 6A and figs. S9 and S10, A and A’). Thus, we wondered whether these genes might exert sex-specific effects on survival. For each transcript, we separated all male and female cases into low and high expression groups based on the amount of expression that distinguished mc5 or fc3 from the other male or female clusters. We then determined the effect on OS for each transcript in each sex separately. Last, we compared the effect of the whole gene set on OS between males and females in the combined dataset. None of the distinguishing genes of fc3 exhibited a differential effect on survival in males compared with females (fig. S9). In contrast, although each of the down-regulated cell cycle pathway genes in mc5 affected OS in both males and females, they exhibited a greater effect, as evidenced by smaller HRs, in males compared with females (Fig. 6C and fig. S10, C and C’). Comparing the survival effects of the gene set in males and females, the HRs of the 17 genes were significantly higher in males than in their female counterparts (P = 4.6 × 10−5, Wilcoxon signed-rank test), indicating that the gene set as a whole exerted a greater effect in males than in females, despite almost overlapping expression density of each gene in males and females (Fig. 6A and figs. S9 and S10).

Expression of sex-specific cluster-defining genes correlates with chemotherapy sensitivity

Sex differences in GBM survival could result from many different cellular, tissue, or organismal factors. To further evaluate the potential prognostic value of the mc5.17 and fc3.9 gene signatures, we performed dose-response analyses for TMZ, etoposide, lomustine (CCNU), and vincristine (VCR) in five male and four female primary GBM cell lines (fig. S11) to determine how the expression of mc5.17 and fc3.9 specific genes correlated with median inhibitory concentration (IC50) values.

Only one cell line (male B66) demonstrated appreciable MGMT expression, as measured by Western blot analysis (fig. S12). The TMZ IC50 of this line was less than that of two other male cell lines with no MGMT expression (figs. S11 and S12), indicating that MGMT expression was not a dominant determinant of TMZ resistance in these assays. Absolute IC50 values were calculated from each dose-response curve and correlated with gene expression, as determined by the Illumina HumanHT-12 v4 expression microarray for each cell line. Overall, male cell lines did not exhibit significantly higher absolute IC50 values than female cell lines (Fig. 7A). To determine whether the expression of mc5.17 and fc3.9 genes stratified responses for male and female cell lines, respectively, we calculated Spearman rank correlation coefficients between IC50 values and gene expression. Spearman rank correlation coefficients comparing the expression of the 17 genes in mc5.17 and IC50s were, on average, positive for male cell lines, indicating that low expression of mc5.17 genes correlated with low IC50 values (high treatment efficacy) for each of the four agents (TMZ, etoposide, CCNU, and VCR) (Fig. 7B). In contrast, in female cell lines, low expression of the mc5.17 genes predicted high IC50 values (low treatment efficacy). As a negative control, the distribution of the averaged correlation coefficient of 1000 random gene sets of the same size as mc5.17 with 17 randomly selected genes centered around 0, indicating no correlation, as expected. When the relationship between each fc3.9 gene and IC50 values for these drugs in male or female cell lines was analyzed, treatment efficacy in female but not in male cell lines was predicted by fc3.9 genes. Again, 1000 random gene sets of the same size (9 randomly selected genes) were not correlated with IC50 for any drug (Fig. 7C).

Fig. 7 Expression of cluster-defining genes and response to common chemotherapeutics in vitro.

(A) Absolute IC50 values for TMZ, etoposide, CCNU, and VCR for five male and four female patient-derived GBM cell lines were calculated from six-point dose-response curves for each cell line. Boxplots of IC50 across cell lines by sex are presented (horizontal bar indicates median). Median male and female IC50 values were not significantly different based on two-sample t test. Spearman correlation coefficients of IC50 values for each drug with expression of mc5.17 genes (B), fc3.9 genes (C), or random gene sets are shown in male and female cell lines. For mc5.17 and fc3.9 genes, box plots represent the distribution of the 17 or 9 cluster-defining genes, respectively, and for random gene sets, the box plots represent the distribution of the Olkin-averaged Spearman correlation coefficient (54) of 17 or 9 randomly selected genes per random gene set for 1000 random gene sets. *P < 0.01 compared with random gene sets for each sex. (D) Quantification of the percentage of phospho-histone H3 (pHH3)–positive nuclei in male GBM cells implanted in male (black bars) or female (white bars) nude mice. Tumor-bearing mice were treated with vehicle [dimethyl sulfoxide (DMSO)], TMZ (21 mg/kg per day × 5 doses), or etoposide (20 mg/kg every other day × 3 doses), and pHH3 positivity was determined in a blinded fashion.

These results indicated that sex-specific expression of these genes is predictive of survival in patients with GBM and the in vitro efficacy of common chemotherapies. Together, these findings suggested that survival in patients with GBM may be related to cell-intrinsic sex differences in mechanisms that broadly affect treatment response. To further evaluate the possibility of cell intrinsic sex-specific determinants of treatment response, we generated flank xenografts using a well-described male murine GBM model with complete loss of neurofibromin and p53 function, and activation of the EGF receptor (15, 42). One million male GBM cells were implanted into male (n = 14) or female (n = 14) nude mice. We focused on male cells alone because we have previously characterized their more consistent in vivo tumor-forming potential as compared with their female counterparts. After establishing steady tumor growth, we treated tumor-bearing mice with TMZ (21 mg/kg per day × 5 doses), etoposide (20 mg/kg every other day × 3 doses), or vehicle (DMSO). We chose these two drugs as representative of agents with higher and lower IC50 values when tested in male GBM cells, respectively. We evaluated acute treatment responses by measuring the drugs’ effects on proliferation using quantification of the percentage of nuclei that were positive for the mitotic marker pHH3. Consistent with our earlier results, we found that male cells formed tumors in recipient mice regardless of their sex. However, proliferation was significantly greater (P = 0.0092) in tumors growing in male compared with female mice (Fig. 7D). We further found that, consistent with their in vitro IC50 values, etoposide, but not TMZ, significantly (P = 0.0082) reduced tumor cell proliferation. We did not detect an interaction between recipient mouse sex and drug effects (P = 0.5092), indicating that in this model, although the sex of the microenvironment can influence tumor growth rates, cell-intrinsic effects determine chemotherapy response (table S10).

DISCUSSION

Sex differences are increasingly recognized as important determinants of human health and disease. Although sex differences in incidence, disease phenotype, and outcome are well described and broadly recognized, the molecular bases for sex differences beyond acute hormone actions are poorly understood. Among the obstacles to improved understanding of sex differences is the inconsistent application of methodologies into laboratory-based and clinical research design that can adequately detect and quantify sex differences. As an example, current epidemiological data indicate that in the United States, the male-to-female incidence ratio for GBM is 1.6:1 (5). Although substantial sex differences in the incidence of GBM and other brain cancers have been recognized for decades, large-scale analyses continue to most commonly merge data from both sexes, obscuring discovery of valuable information contained in the sex differences.

Recent exceptions illustrate the value of using sex differences to highlight important elements of cancer biology and clinical response. We found that sex-specific, cell-intrinsic responses to loss of p53 function render male astrocytes more vulnerable to malignant transformation compared with female astrocytes (15). These findings may well relate to the sex differences in glioma incidence and are consistent with other data describing sexual dimorphism in the p53 pathway, including radiographic sex differences in men and women with GBM as a function of their p53 mutational status (43). Understanding the molecular basis for sexual dimorphism in the p53 pathway and what it means with regard to cancer biology and clinical oncology remains an important area of research.

These studies emphasize that analyses without consideration of sex can obscure critical elements of biology and, in aggregate, highlight the importance of parallel but separate analyses of male and female cells, male and female animals, and male and female patients. Here, we applied the JIVE algorithm to decompose male and female GBM transcriptome datasets of TCGA into joint and sex-specific components. We found that male and female patients with GBM cluster into five distinct male and female subtypes that are distinguished by gene expression and survival. These clusters, which were identified using the TCGA transcriptome dataset, were subsequently validated in three independent datasets. Although GBM has recently been identified as a “low sex effect” cancer at the transcriptome level (44), our analyses indicate that even genes with similar expression in males and females can impart substantial sex-specific effects on survival and yield mechanistically important information. Together with the sex-specific effects of p53 loss (15) and Arlm1 variants (45), these data suggest that the cellular and organismal sex context of gene expression affects the consequences of oncogenic events. A similar mechanism was invoked to explain the sex-specific effects of AC8 polymorphisms, which increased the risk of low-grade glioma in females but decreased the risk in males with NF1 (12).

Most compelling in this regard are the molecular features of the longer surviving subtypes of male and female patients. IDH1 mutation is a molecular marker of a distinct form of GBM that is associated with better outcome (46, 47). Here, we demonstrate that IDH1 mutation exhibits sex-specific survival benefit. In the combined dataset, almost all IDH1 mutant female tumors were assigned to fc3. This was the only female cluster with distinctly better outcome. In contrast, IDH1 mutations were distributed across all male clusters. Moreover, although the numbers were small, fc3 conferred a survival advantage regardless of IDH1 mutational status, whereas IDH1 mutation further stratified survival in mc5 cases. Thus, the predictive value of IDH1 mutation can be better defined in a sex-specific context. This finding is in contrast to a recent immunohistochemical analysis of IDH1 mutation in a single cohort of 105 patients (48). In this cohort, there were a total of nine IDH1 mutant tumors: four in males and five in females. The difference in survival for male patients (n = 61) with and without IDH1 mutations reached statistical significance. This was not true for the female patients (n = 44), but the sample size was small and the results were not validated in an independent cohort. The sex-specific impact of IDH1 mutation on survival will require additional evaluation.

Fc3 and mc5 shared a distinguishing signature in calcium/calmodulin signaling, with a particular representation of genes essential for synaptic function. They diverged in other molecular features, with mc5 exhibiting down-regulation of mitotic spindle and cell cycle regulatory genes and fc3 exhibiting a down-regulation of integrin signaling pathway components. Most compelling was the sex-specific effect on survival of genes within the cell cycle regulatory pathway, despite the fact that the component transcripts were similarly expressed in male and female tumors. These observations are consistent with the hypothesis that sex effects in cancer cannot simply be defined by gene expression but rather need to include the potential sex differences in gene effect. A similar observation regarding sex differences in MGMT promoter methylation was recently published (48).

Among the striking results of this study is the potential harmonization between the sex differences in gene expression, in vitro drug sensitivity, MRI measures of tumor dispersion and proliferation, MRI measures of treatment response, and patient survival. The gene expression analysis identified down-regulation of cell cycle progression and down-regulation of integrin signaling as correlated with best survival in male and female patients, respectively. Expression of the 17 and 9 gene signatures that distinguished the longest surviving male and female cohorts, respectively, also correlated, in a sex-specific manner, with in vitro drug sensitivity as measured by IC50 values for a panel of primary GBM cell lines. Moreover, we found evidence that MRI-based predictors of survival may differ for males and females with GBM. These predictors are based on measures of rates of proliferation and invasion. It would be premature to overemphasize the relationship between the concordance in these measures and patient outcomes, but such multiscale correlations are the goal of projects such as The Human Tumor Atlas. In this regard, these findings may provide an example of how sex differences in cancer can be productively incorporated into these efforts.

The current study has several limitations that should be considered. First, the MR image analysis was performed retrospectively on a cohort collected from multiple medical institutions over the course of many years. Although the inclusion criteria were designed to mitigate interpatient variability, imaging and treatments may have varied between institutions and patients. In addition, the MR images were segmented and validated by trained individuals, which could introduce some interoperator variability into calculated growth velocities. Second, although all genomic data, including DNA sequence and copy number variation, mRNA expression, as well as protein expression and posttranslational modification data, should be collectively analyzed for sex-specific features, such a comprehensive data source of a large enough sample size for each component has yet to be established for sex-specific modeling. TCGA, the largest repository of GBM profiling results, contains gene expression data for over 500 cases but more limited data of other types. Thus, we defined sex-specific clusters primarily using mRNA gene expression data. We did, however, perform independent analyses after excluding IDH1 mutant cases to better reflect current thinking regarding IDH1 mutant GBM as a distinct disease entity from IDH1 WT GBM. In some cases, IDH1 mutation information was not available and G-CIMP status was used when appropriate as a surrogate. This analysis indicated that the sex-specific clusters were evident before and after excluding IDH1 mutant cases, suggesting that they are robust across GBM disease types. With regard to clinical outcomes, we analyzed both DFS and OS in the TCGA dataset. Although sex-specific effects were greatest for DFS, this parameter was not reported in all the validation datasets, and thus, we were limited to OS for the merged analyses. Missing IDH1 and G-CIMP information in most of the validation samples resulted in relatively small sample sizes when IDH1 mutation status was considered for survival analysis. Third, although evaluating chemosensitivity in the panel of patient-derived cell lines yielded important data, the relatively small number of cell lines limits their interpretation. Additional studies, with more cellular isolates and in vivo treatments, particularly those that might yield insights into the mechanisms underlying sex-specific effects of chemotherapy, will be necessary before we can rationally apply these results to clinical trial design. Last, although it is not yet possible to ascribe a specific fraction of the survival differences between male and female patients with GBM to any of the sex differences we describe, the current study does suggest that greater precision in the stratification of patients with GBM may be achieved through sex-specific molecular subtyping and that improvements in GBM outcome might be possible with sex-specific approaches to treatment, including blocking cell cycle progression in male patients and targeting integrin signaling in female patients.

MATERIALS AND METHODS

Study design

This study was designed to investigate sex differences in GBM incidence and outcome. We performed three kinds of analyses to achieve this goal. First, we applied a previously validated MR image analysis method to calculate tumor growth velocities, and an established mathematical model to estimate tumor proliferation and invasion rates. These parameters were derived from patient image and clinical data retrospectively collected from multiple institutions and sourced through the clinical research database at the Mayo Clinic (Phoenix). We also evaluated sex-specific correlations between these growth parameters and survival. Second, we derived sex-specific molecular subtypes of GBM through discovery using TCGA data (49) and validation in three additional datasets: GSE16011 (50), GSE13041 (51), and the REMBRANDT study (52, 53). Third, we measured sex differences in the in vitro cytotoxic effects of four common chemotherapeutics in a panel of nine (five male and four female) patient-derived GBM cell isolates. In addition, we evaluated the relative contributions of cellular sex and the sex of the microenvironment to therapeutic responses of two different chemotherapeutics by parallel implantation of male murine GBM cells into equal numbers of male and female mice. All pathology analyses were performed in a blinded fashion. All animal and human studies were approved by the appropriate Animal and Human Studies committees at the Mayo Clinic (Phoenix) and the Washington University School of Medicine.

Detailed Materials and Methods can be found in the Supplementary Materials.

SUPPLEMENTARY MATERIALS

www.sciencetranslationalmedicine.org/cgi/content/full/11/473/eaao5253/DC1

Detailed Materials and Methods

Fig. S1. Sample timeline for MRI analysis of patient data.

Fig. S2. Relation between RANO criteria, first tumor growth velocity, and survival.

Fig. S3. Validation of sex differences in correlations between presurgical D and rho and survival.

Fig. S4. Geometric illustration of JIVE decomposition.

Fig. S5. Joint structure of male and female GBM samples relative to TCGA subtype signatures.

Fig. S6. Expression variation in TCGA data explained by the JIVE components.

Fig. S7. Identification of the sex-specific clusters by consensus clustering for female and male GBMs.

Fig. S8. Sex-specific cluster survival benefits in merged TCGA, GSE16011, and GSE13041 datasets.

Fig. S9. fc3-defining genes and survival.

Fig. S10. mc5-defining genes and survival.

Fig. S11. TMZ, etoposide, VCR, and CCNU dose-response curves for five male and four female primary GBM cell lines.

Fig. S12. Western blot analysis of MGMT expression in female and male patient-derived GBM cell lines.

Fig. S13. Distribution of the coefficients of variation in the TCGA data.

Fig. S14. Analysis flowchart.

Table S1. Sex and TCGA subtype assignments in the four datasets.

Table S2. Gene expression signatures for sex-specific clusters.

Table S3. Median survival time, P values, and HRs associated with survival differences in sex-specific clusters.

Table S4. Survival by IDH1 status—Female IDH1 WT.

Table S5. Survival by IDH1 status—Male IDH1 WT.

Table S6. Survival by IDH1 status—Male IDH1 mutant.

Table S7. Survival by IDH1 status—Female IDH1 mutant.

Table S8. Survival analysis of sex—IDH1 interaction.

Table S9. Pathway analysis of fc3 and mc5.

Table S10. Statistical analysis of the interaction between sex and chemotherapy treatment in vivo.

Table S11. Patient-specific proliferation, invasion, hypoxia, necrosis, angiogenesis (PINHA) values.

Table S12. Parameters for calculation of tumor growth velocities.

REFERENCES AND NOTES

Acknowledgments: The content is solely the responsibility of the authors and does not necessarily represent the official views of the NIH or government. Funding: This work was supported by grants from the NIH [R01 CA174737 (to J.B.R.); R01 NS060752, R01 CA164371, U54 CA210180, U54 CA143970, and U54 CA193489 (to K.R.S.); and K08 NS081105 and R01 NS094670 (to A.H.K.)]; the Children’s Discovery Institute of Washington University and Joshua’s Great Things (to J.B.R.); the James S. McDonnell Foundation, the Ivy Foundation, and the Mayo Clinic (to K.R.S.); and NIH U01 CA168397 and Ben & Catherine Ivy Foundation (to M.E.B.). Author contributions: J.B.R. and K.R.S. conceived the experiments. W.Y., P.W., J.L., K.R.S., and J.B.R. designed the experiments. W.Y., S.J.T., N.M.W., E.C., K.W.S., N.W., J.D.L., M.E.B., J.S.B.-S., K.R.S., J.L., and J.B.R. analyzed the data. A.H.K. generated the GBM cell lines. S.J.T., N.M.W., and J.B.R. generated the histopathology data. All authors contributed to the writing and editing of the manuscript. Competing interests: A.H.K. served as a paid consultant for Monterris and received a Stryker Research Grant to evaluate outcomes after neurosurgical use of their dural substitute. K.R.S. served as a scientific advisor for the James S. McDonnell Foundation. All the other authors declare that they have no competing interests. Data and materials availability: In accordance with Mayo Clinic institutional review board (IRB) restrictions, all data used for imaging-based analyses in this manuscript are provided as tables in the Supplementary Materials. As is typical of all protected patient data, we have a clinical research repository that can only be shared per specific IRB requirements. By reasonable request, by contacting K. R. Swanson, we may be able to create a data sharing agreement to share other data with interested institutions, following our institutional guidelines as well as those of the recipient. All other data associated with this study are present in the paper or Supplementary Materials.
View Abstract

Navigate This Article