Abstract
Like BRCA2, MAGEC3 is an ovarian cancer predisposition gene that has been shown to have prognostic significance in ovarian cancer patients. Despite the clinical significance of each individual gene, no studies have been done to assess the clinical significance of their combined expression. We therefore sought to determine the relationship between MAGEC3 and BRCA2 expression in ovarian cancer and their association with patient characteristics and outcomes. Immuno-histochemical staining was quantitated on tumor microarrays of human tumor samples obtained from 357 patients with epithelial ovarian cancer to ascertain BRCA2 expression levels. In conjunction with our previously published MAGEC3 expression data, we observed a weak inverse correlation of MAGEC3 with BRCA2 expression (r= -0.15; p<0.05) in cases with full-length BRCA2. Patients with optimal cytoreduction, loss of MAGEC3, and detectable BRCA2 expression had better overall (median OS: 127.9 vs 65.3 months, p=0.035) and progression-free (median PFS: 85.3 vs 18.8 months, p=0.002) survival compared to patients that were BRCA2 expressors with MAGEC3 normal levels. Our results suggest that combined expression of MAGEC3 and BRCA2 serves as a better predictor of prognosis than each marker alone.
Simple Summary Early detection of ovarian cancer is a significant clinical challenge, with most women presenting with advanced stages of disease at initial diagnosis. The aim of this study was to evaluate the role of MAGEC3 and BRCA2 in epithelial ovarian cancer progression. We evaluated the effect of MAGEC3 and BRCA2 on the other’s expression. We tested this in humans using immunohistochemical staining of human tumor samples obtained from patients with epithelial ovarian cancer (n = 357). We found a weak inverse correlation between MAGEC3 and BRCA2 expression in epithelial ovarian cancers. Further, our data suggests that the combined expression of MAGEC3 and BRCA2 may be a better predictor of outcomes in patients than the individual markers alone.
1. Introduction
Ovarian cancer is a heterogenous disease and it remains the leading cause of mortality from gynecologic malignancies, with more than 19,000 new cases and 12,810 deaths expected in the United States in 2022 [1]. Epithelial ovarian cancer (EOC), which accounts for over 90% of all ovarian cancers, is comprised of five distinct subtypes, amongst which high-grade serous ovarian carcinoma (HGSOC) is the most commonly diagnosed [2]. HGSOC is also the most lethal subtype, accounting for about 70% of all epithelial ovarian cancer cases [3, 4] and 70%–80% of all ovarian cancer deaths [5] due to it being diagnosed at an advanced stage [6].
Since most women present with advanced stages of disease at initial diagnosis [4, 7, 8], early detection of ovarian cancer remains an increasingly important clinical challenge [2, 8, 9]. There is a compelling argument for novel ovarian cancer screening methodologies, given that the current strategy results in early detection of ovarian cancer in only ∼20% of patients [10]. Advances in ovarian cancer diagnosis include the use of genetic testing to identify mutated genes associated with epithelial ovarian cancer in patients with a known family history of breast and ovarian cancer, such as TP53, PIK3CA, KRAS, and BRCA1/2 [2, 8].
BRCA1/2 genes encode proteins involved in the repair of DNA double-strand breaks via the homologous recombination pathway [11-14]. Germline and somatic mutations, as well as promoter methylation, are the common causes of BRCA1/2 dysfunction [15, 16] which impairs DNA repair leading to increased risk of tumorigenesis [17]. BRCA2 mutations specifically are associated with a 15-30% lifetime risk of ovarian cancer [18]. Most of the deleterious mutations that have been identified in ovarian cancer families are frameshift (insertions or deletions) mutations or nonsense mutation which result in premature termination of the protein-coding sequences [19-22]. Resultant BRCA2 mutant proteins retain their N-terminal transactivation domain which is enables its recombinational repair and checkpoint roles [23-25] but the C-terminal, which harbors the nuclear localization signals (NLS), is lost [25]. The inability to translocate into the nucleus explains why these mutants are nonfunctional and unable to facilitate the DNA repair process in most cases [20, 25, 26].
BRCA1/2 mutations serve as a marker of vulnerability to PARP inhibitors (olaparib, rucaparib, and niraparib). These have become viable treatment options for patients with advanced epithelial ovarian cancer due to enhanced clinical effectiveness and tolerable toxicity [27-31]. While women with BRCA1/2 mutations retain the best response to PARP inhibitors [32], tumors with genomic instability still show clinically significant response to therapy [33], prompting researchers to consider whether other markers like ATM, FANC A/F, CHK2, RAD51B/C, and CDK12 could be used as de facto surrogate markers of deficient homologous repair [34, 35] to expand the pool of women who could benefit from PARP inhibitor therapy.
One such potential marker is MAGEC3. MAGEC3 is a member of the “Type II” melanoma antigen gene (MAGE) family by DNA sequence homology and has low but ubiquitous protein expression in most tissues with high expression in testes and tonsil [36]. In the context of cancer, members of the family have shown varying functions in DNA repair and cell cycle regulation, [37, 38] as well as potential application as markers of stemness in cancer cells [39]. Biochemically, MAGE proteins bind specifically to E3 RING ubiquitin ligases through their MAGE homology domains (MHDs) and regulate protein ubiquitination by enhancing E3 ligase activity [37, 38, 40, 41]. Previously, we described a linkage between early-onset ovarian cancers [42] and MAGEC3 [43, 44] and more recently found that MAGEC3 is a potential prognostic biomarker in ovarian cancer, as loss levels were associated with a favorable progression-free survival (PFS) in ovarian cancer patients [36].
The relationship between both BRCA2 and MAGEC3 protein expression and their clinical significance in ovarian carcinoma is unknown. Immunohistochemistry (IHC) is a low-cost and easily accessible technique that allows for the query of both N- and C-terminal domains of the BRCA2 protein, including protein truncation of the C-terminal of the BRCA2 protein which severely diminishes protein function [45]. Similarly, IHC has also been used to detect expression of the MAGEC3 protein [36, 46]. Herein, we describe the relationship between MAGEC3 and BRCA2 in patient samples using IHC and tested the hypothesis that their co-expression is correlated with clinical outcome.
2. Materials and Methods
2.1. Ovarian Cancer Patients
Two cohorts of patients were used: (1) Patients with advanced high-grade serous ovarian carcinoma treated at Roswell Park Cancer Institute (RPCI) who had previously had BRCA2 genetic testing. This cohort was used for the pilot experiment and classified into three categories based on the results of the genetic testing: BRCA2-Wildtype, BRCA2-Mutant, and BRCA2-Loss. (2) Patients with a primary diagnosis of ovarian cancer treated at RPCI by maximal debulking surgery and first-line platinum-based chemotherapy between 2000 and 2012 were obtained from a curated database under an institutional review board-approved protocol. Tumor blocks from these cohorts were taken during primary debulking surgery and prior to the initiation of systemic therapy.
In total, there were n = 3 ovarian cancer patients from cohort 1 with sequenced tumor samples processed as whole tissue sections for the pilot experiment and n = 357 ovarian cancer patients from cohort 2 with unsequenced tumor samples for the tissue microarray analysis. The detailed handling protocol for the tissue specimens has been previously described [47].
2.2. Whole Tissue Sections
The whole tissue sections derived from formalin-fixed and paraffin-embedded tumor tissues obtained from three sets of sequenced ovarian tumor cases in cohort 1 to qualify antibody reactivity. IHC H-scores were generated using n=30 photographs of randomly selected fields on each slide.
2.3. Tissue Microarray (TMA) Construction
The tissue microarrays (TMA) were constructed using formalin-fixed and paraffin-embedded tumor tissues punched from each donor block from cohort 2. Triplicate TMAs (containing three representative cores from each tumor) were prepared and stained to assess tumor heterogeneity. These TMAs included all epithelial ovarian cancer subtypes, non-tumor control tissues, and reference cores for normalization across different blocks.
2.4. Antibody Validation
The Human Protein Atlas [48] (HPA052067) was used to vet the BRCA2 antibodies, in accordance with the International Working Group for Antibody Validation [49] and were evaluated and authorized for IHC and supported for antigen specificity. The BRCA2 antibodies are directed against the BRCA2 C-terminal amino acids (2587–2601 amino acids, Abcam ab53887, Cambridge, MA, USA) and BRCA2 N-terminal amino acids (100–150 amino acids, Proteintech 19791-1-AP, Rosemont, IL, USA).
The MAGEC3 antibody has been previously described [36].
2.5. Immnuhistochemistry
Formalin-fixed paraffin sections were cut at a thickness of 4μm, placed on glass slides and heated at 60°C for 1hour. The TMA sections were then deparaffinized with xylene and ethanol and rinsed in water. For antigen retrieval, slides were pretreated in an antigen retrieval high pH solution (catalog #GV804, Dako, Santa Clara, CA, USA) for 20 min and allowed to cool for 20 min at room temperature. The slides were incubated with a solution of 3% H2O2 for 10 minutes. “BRCA2 N-terminal” antibody (catalog #19791-1-AP, Proteintech, Rosemont, IL, USA) was added to one set of slides at 1:500 dilution while “BRCA2 C-terminal” antibody (catalog# ab53887, Abcam, Cambridge, MA, USA) was added to the other set of slides at 1:500 dilution and incubated for 60 minutes at room temperature. This was followed by incubation with biotinylated goat anti-rabbit secondary antibody (catalog # 14708, Cell Signaling Technology, Boston, MA, USA), 1:1000 dilution for 15 minutes at room temperature. ABC reagent (catalog #PK 6100, Vector Labs, Burlingame, CA) was used for signal enhancement and applied for 30 minutes. Slides were incubated for 5 minutes with 3,3’-diaminobenzidine (DAB) substrate (catalog #K3467, Dako, Carpinteria, CA) and then counterstained for 20 seconds with DAKO Hematoxylin (catalog #CS700, Dako, Carpinteria, CA). Slides were dehydrated in several baths of graded alcohols and xylenes prior to being cover slipped.
TMAs were previously stained for MAGEC3 as described [36].
2.6. Digital Image Analysis
Digital images of both IHC-stained whole tissue sections and IHC-stained TMA slides were captured using a whole slide scanner at 10X magnification (Leica ScanScope, Aperio Technologies, UK), and the digital images were stored in SVS format (Aperio). The digital images were retrieved using a file management web interface (eSlide Manager, Aperio) and reviewed with the server software (ImageScope, Aperio, version 12.3.3.5048). Digital images of representative whole tissue section fields and TMA cores were extracted through the server software (ImageScope, Aperio, version 12.3.3.5048) at 10× magnification and quantified into H-scores using an automated IHC profiler built in a digital image analysis software (ImageJ version 1.52a, National Institute of Health, USA) [50].
2.7. Statistical Methods
After the IHC scoring was complete, the optimal cutpoint of the IHC staining score for BRCA2 expression was determined using H-scores obtained for the C-terminal of the sequenced ovarian tumor cases with the aid of the R package, cutpointsr [179]. We stratified the BRCA2 C-terminal H-scores into two different groups (BRCA2 expressors and BRCA2 non-expressors).
All statistical analyses were performed using the R 3.1.2 statistical computing language. Associations between BRCA2/MAGEC3 expression and categorical variables were tested by chi-square tests and continuous variables were tested by one way ANOVA tests. Survival probabilities were calculated by Kaplan-Meier analysis using log-rank testing. A nominal significance threshold of 0.05 was used unless otherwise specified. The multivariate analysis included stage, categorized as early (I, II, IIIA, or IIIB) or late (IIIC or IV), grade (1 vs. 2/3), debulking status (optimal vs. suboptimal), and platinum-sensitive versus refractory disease.
3. Results
3.1. Validation of IHC Measurements of BRCA2 Protein Expression in Cancer Tissue
We first tested whether IHC BRCA2 protein expression confirmed the predicted expression by molecularly characterized and sequenced ovarian tumor samples. The adopted two-antibody IHC strategy, which utilizes antibodies specific towards the N-terminus and C-terminus, allows for discrimination between the presence of normal protein, protein truncation, and loss of protein expression in the tumor samples (Figure 1A). These antibodies were applied to three cases confirmed to have BRCA2 wild type (WT), BRCA2 truncating mutation in exon 11 (Mutant), and total loss of BRCA2 (Loss) through sequencing. Visual analysis of IHC staining shows that, as expected, wildtype samples have adequate staining observed with both the N-terminal and C-terminal antibodies, mutant samples have staining only at the N-terminal and not the C-terminal, and BRCA2 loss has poor staining intensity seen with for both antibodies (Figure 1B).
To further validate this strategy, the staining intensity and percentage of positive tumor cells was quantified using an automated IHC H-scoring algorithm built in ImageJ. The resulting H-scores show that N-terminal expression levels are significantly different between BRCA2 wild type, BRCA2 mutant, and BRCA2 loss in sequenced cases (WT vs Mutant: t-test, p<0.001; WT vs Loss: p<0.001; Mutant vs Loss: p<0.001) (Figure S1). However, given that a truncating mutation renders the protein non-functional, samples that demonstrate lack of C-terminal staining should functionally mimic BRCA2 loss. The difference in C-terminal expression levels when BRCA2 wildtype was compared to BRCA2 mutant and BRCA2 loss cases was observed to be statistically significant (WT vs Mutant: t-test, p<0.001; WT vs Loss: p<0.001), while the difference in C-terminal expression levels for the BRCA2 mutant and BRCA2 loss cases was not statistically significant (Mutant vs Loss: t-test, p=0.9) (Figure 2A). Given its ability to differentiate wild-type cases from the functionally impaired mutant and loss cases, C-terminal scoring was used for further analyses. C-terminal scores for the BRCA2 mutant and BRCA2 loss cases were combined and termed “BRCA2 non-expressors” while the BRCA2 wild-type cases were termed “BRCA2 expressors.” When the same IHC strategy was applied to 357 unsequenced patient samples with unknown BRCA2 status, we used the C-terminal expression levels of the BRCA2 expressors and non-expressors to determine an optimal cutoff point for stratifying the unsequenced cases. The optimal cutoff point utilized to stratify the larger unsequenced cohort of patients into “BRCA2 expressors” and “BRCA2 non-expressors” was a H-score of 49.5 (Figure 2B).
MAGEC3 protein expression quantitation by IHC for patients in this cohort was previously reported in a recent publication from our lab [36]. Expression levels were dichotomized at the median into either “MAGEC3 Normal” or “MAGEC3 Loss” based on the findings that higher levels of MAGEC3 were not significantly different from expression in normal ovary tissue (t-test, p = 0.368), while lower levels of MAGEC3 were significantly lower (t.test, p <0.001). Given our interest in the combined effect of MAGEC3 and BRCA2 expression levels on patient prognosis, we used MAGEC3 expression data from the prior publication as well as the BRCA2 levels ascertained in this study for subsequent analyses.
3.2. Clinical Characteristics of Ovarian Cancer Microarray Patient Population
A total of 357 tumor tissues from patients with ovarian and primary peritoneal cancer were analyzed using IHC. From the demographic information of de-identified patients in this study (Table 1), the mean age of diagnosis was 63 years (range 21-93 years) and 52.1% of patients were diagnosed after 2006. Most of the patients were white (334; 94.9%), non-Hispanic (350; 99.4%), and presented with advanced/late-stage disease (81.5%), poorly differentiated tumors (73.5%) and serous histology (79.6%). Platinum-sensitive disease was demonstrated in 156 of the 357 patients (55.5%) with 125 patients classified as platinum resistant or treatment refractory (42.5%). For survival data, the patients were followed for an average of 65.2 months (maximum 274.2 months).
Using the optimal cutoff point calculated above, 299/257 (83.8%) of the patients were classified as BRCA2 expressors while 58/357 (16.2%) were BRCA2 non-expressors. Using the previously reported stratification of MAGEC3, 173/357 (48.5%) of patients in this cohort are MAGEC3 loss while 184/357 (51.5%) are MAGEC3 normal. When considering the combination of BRCA2 and MAGEC3 levels, 34 (9.5%) were BRCA2 non-expressors and showed MAGEC3 loss levels, 139 (40%) were BRCA2 expressors with MAGEC3 loss levels, 24 (6.7%) were BRCA2 non-expressors with MAGEC3 normal levels, and 160 (44.8%) were BRCA2 expressors with MAGEC3 normal levels. Through extensive statistical analysis, BRCA2 expression and normal MAGEC3 levels were associated with cases ascertained after 2006 (Chi-sq test, p<0.001), ovary as primary site of origin (Chi-sq test, p<0.001), poorly/undifferentiated tumors (Chi-sq test, p=0.03) and serous histology (Chi-sq test, p<0.001).
3.3. BRCA2 Protein Shows a Weak Negative Correlation with MAGEC3 Expression in Ovarian Cancer Tumor Samples
The H-scores from the semi-quantitative analysis of MAGEC3 TMAs in our previous study [36] were correlated with the H-scores of BRCA2 TMAs in this current study using a two-sided Pearson’s correlation test. The initial assessment of the relationship between MAGEC3 and BRCA2 expression in the tumor samples revealed no correlation (Figure 3A). However, when the cases were stratified based on BRCA2 expression, the MAGEC3 H-scores were shown to be inversely correlated with the H-scores of BRCA2 expressors in the ovarian tumor samples (r = -0.15; p<0.05) (Figure 3B), while a positive correlation was observed in the BRCA2 non-expressors (r=0.32; p<0.05) (Figure 3C).
3.4. MAGEC3 and BRCA2 Association with Prognosis in Epithelial Ovarian Cancer
The impact of age at diagnosis, disease stage, grade, histology, cytoreduction status, and MAGEC3/BRCA2 expression on clinical outcome was analyzed using both univariate and multivariate Cox proportional hazards regression analyses (Table 2). The age at diagnosis, grade, and R0 cytoreduction were univariately associated with prognosis in epithelial ovarian cancers. Multivariate regression analysis indicated that patients with advanced stage disease (HR = 3.28, 95% CI 2.05-5.26, p<0.001) and poorly/undifferentiated tumors (HR = 1.34, 95% CI 1.02-1.76, p=0.04) were statistically significant predictors of survival.
Independent evaluation of the relationship between BRCA2 expression and patient outcomes in our cohort of advanced ovarian cancer cases revealed that BRCA2 expression levels did not impact overall survival (log-rank p = 0.198) and progression-free survival (log-rank p = 0.15) in patients with optimal cytoreduction (Figure S2), whereas independent evaluation of MAGEC3 expression in our previous study revealed that women with MAGEC3 loss had better progression-free survival (log-rank p=0.002) [36]. In patients with optimal cytoreduction, subgroup analyses of combined expression of MAGEC3 and BRCA2 showed that patients that were BRCA2 expressors with MAGEC3 loss levels had better overall survival compared to patients that were BRCA2 expressors with MAGEC3 normal levels (median OS: 127.9 vs 65.3 months; log-rank p=0.035) (Figure 4A). This was also observed for progression free survival (median PFS: 85.3 vs 18.8 months; log-rank p=0.002) (Figure 4B). Conversely, in BRCA2 non-expressor cases, having MAGEC3 loss levels conferred no survival benefit over MAGEC3 normal levels in overall survival (log-rank p=0.677) nor progression-free survival (log-rank p=0.753) (Figure 4C+D).
4. Discussion
The aims of this study were to evaluate BRCA2 and MAGEC3 for their influence in epithelial ovarian cancer progression and to assess the clinical significance of their combined expression. Like BRCA2 [51], evaluation of MAGEC3 protein expression has been identified as a novel prognostic biomarker in ovarian cancer [36]. However, the relationship between both BRCA2 and MAGEC3 protein expression and their clinical significance in ovarian carcinoma is unknown.
The growing importance of IHC as a tool for detecting protein dysfunction via protein expression profiles, combined with its ease of use and practicality in identifying the protein of interest within the tumor’s cellular compartment, makes it a highly appealing technique to measure validated biomarkers. It proved to be a viable and cost-effective approach for determining BRCA2 protein expression in tumor samples for this study, as similarly reported by other studies [52, 53], as well as for determining MAGEC3 protein expression [36]. We used a dual antibody strategy to discriminate differential expression of BRCA2 at both the N-terminal and C-terminal. Because of its ability to differentiate the wild-type case from mutant and loss cases, C-terminal expression was used to define BRCA2 expression in a large unsequenced cohort. Prior studies also support the idea that the C-terminal antibody could be useful in screening cancers for BRCA mutations as well as BRCA2 protein expression in patients with unknown mutation status [45, 52, 54]. Using this approach, we identified 83.8% (299/357) of patients as BRCA2 expressors and 16.2% (58/357) as BRCA2 non-expressors.
In the assessment of the relationship between MAGEC3 and BRCA2 proteins in the tumor samples from the unsequenced cohort, there was no significant correlation of protein expression when looking at all cases. However, when we looked at only BRCA2 expressors, there was an inverse correlation with MAGEC3 expression. An interesting conjecture emerges that, like other proteins in the MAGE family, MAGEC3 is interacting with an E3 ligase which leads to the ubiquitination and subsequent degradation of a protein, in this case BRCA2. Because this downregulation was only observed in BRCA2 expressors, it is possible that this downregulation will only occur in cells expressing a full-length BRCA2 protein. Preliminary in vitro data from our lab supports this hypothesis, as over-expression of MAGEC3 results in a reduction of BRCA2 protein level in both fibrosarcoma (HT1080) and ovarian cancer (SKOV3) cell lines [55]. It is important to note that both the HT1080 and SKOV3 were reported to have a wild-type BRCA2 gene [56].
In the unsequenced cohort of patients, survival analyses for overall and progression-free survival showed no significant relationship between BRCA2 expression and patient outcomes; whereas our previous study of the same cohort showed that women with MAGEC3 loss had better progression-free survival [36]. For ovarian cancer patients with optimal cytoreduction, combined expression of MAGEC3 and BRCA2 showed that patients that were BRCA2 expressors with MAGEC3 loss levels had better overall and progression-free survival compared with patients who were BRCA2 expressors with normal MAGEC3 levels. Interestingly, in patients that were BRCA2 non-expressors, MAGEC3 loss levels conferred no survival advantage over MAGEC3 normal levels for both overall survival and progression-free survival. Given our results showing that MAGEC3 down-regulates BRCA2, we anticipated that higher levels of MAGEC3 would mimic BRCA2 loss status resulting in susceptibility to platinum-based therapies and prolonged survival. However, we note that the correlation observed in these patients was weak and that the in vitro experiments which showed downregulation of BRCA2 were performed at supraphysiological levels of MAGEC3 which would enhance its downregulating effects. A possible explanation for the trends observed is likely related to MAGEC3’s additional role in augmenting DNA repair [36]. Additional unpublished data from our lab revealed that ovarian cancer cells expressing MAGEC3 fare better under cisplatin (CDDP) insult by clearing cisplatin adducts [55]. These results led us to conclude that MAGEC3 is a prognostic marker in cases that are BRCA2 expressors and that its role in DNA repair, as opposed to its weaker role in downregulating BRCA2, is responsible for the trends observed. Because MAGEC3 only retains its prognostic ability in BRCA2 expressors, the combined expression remains a better predictor of outcomes in patients than the individual markers alone.
Limitations of our study include the use of only three sequenced ovarian tumor cases with known genomic status to develop our IHC assay, as well as the use of retrospective cohorts for the ovarian cancer TMAs. To reduce batch-processing effects and resulting slide-to-slide variation, the scores were normalized using normal control cores. Slide staining was done by the pathology core facility at Roswell Park Comprehensive Cancer Center using appropriate antibodies for IHC.
Despite these limitations, our study has considerable strengths including the analysis of a large cohort of patients to provide insight into the relationship between MAGEC3 and BRCA2 protein expression in ovarian cancer cases. Additionally, the integration and utilization of the automated semi-quantitative IHC profiler built into Image J for assessment of the expression levels of the proteins eliminated any visual bias and high-level observer variation that might have occurred with conventional and manual pathological TMA slide assessment. This was reported as effective in similar studies [57, 58].
5. Conclusions
The combined expression of MAGEC3 and BRCA2 serves as a better predictor of prognosis than either marker alone, as MAGEC3 only retains its prognostic significance in BRCA2 expressors. Ovarian cancer cases expressing BRCA2 with low levels of MAGEC3 fare better than those with normal levels. This is contrary to expectations, given MAGEC3’s role in downregulating BRCA2. Given the role of MAGEC3 in augmenting DNA repair, we propose that loss of this augmentation explains the favorable trend of patients with low MAGEC3 levels. Further work is needed to confirm this hypothesis.
Data Availability
All data produced in the present study are available upon reasonable request to the authors
Supplementary Materials
The following supporting information can be downloaded at: www.mdpi.com/xxx/s1, Figure S1: BRCA2 N-terminal expression in the sequenced ovarian tumor cases.; Figure S2: Kaplan–Meier survival curves for independent evaluation of the relationship between BRCA2 expression and patient outcomes in the larger cohort of patients with unsequenced advanced ovarian cancer cases.
Author Contributions
Conceptualization, K.H.E., E.B.O. and I.A.; methodology, E.B.O., J.E.J., I.A., M.M., S.R.R. and K.H.E.; software, K.H.E.; formal analysis, E.B.O., J.E.J., E.I., S.R.R., E.M.D. and K.H.E.; investigation, E.B.O., J.E.J., I.A., M.M., K.H. and K.H.E.; resources, I.A.,M.M. and K.H.E; data curation, E.B.O., J.E.J., I.A. and K.H.E.; writing—original draft preparation, E.B.O. and K.H.E.; writing—review and editing, E.B.O., J.E.J., I.A., M.M.D., E.I., E.M.D., S.R.R., K.L.N, J.E.O. and K.H.E.; visualization, E.B.O, J.E.J., I.A., M.M.D., S.R.R., K.H.E., M.M. and K.H.; supervision, S.R.R., K.L.N., J.E.O. and K.H.E.; project administration, E.B.O.; funding acquisition, K.H.E. All authors have read and agreed to the published version of the manuscript.
Funding
This research was funded by Department of Defense grants OC170368 and PC180449 (to K.H.E.), the Roswell Park Cancer Center Support Grant P30CA016056, and grants from the Roswell Park Alliance Foundation.
Institutional Review Board Statement
The study was conducted in accordance with the Declaration of Helsinki and approved by the Institutional Review Board of Roswell Park Comprehensive Cancer Center (protocol code BDR109119, 2 December 2019).
Informed Consent Statement
Our research was approved as a non-human subject research study; thus, patient consent was waived. The research was classified as non-human subject research because it utilized de-identified tissue specimens obtained for other objectives, such as under another study design, as part of an ongoing biorepository, or as residual material following a clinical intervention. If the specimens were coded, the holder of the re-identification key was prohibited by institutional regulations or a direct agreement from providing the key to the investigator or other research employees under any circumstances. Furthermore, the investigator agreed not to seek for Identifiers or re-identification. Roswell Park’s biorepositories use an honest broker system to link specimens with associated data and provide it to researchers in a HIPAA-compliant, de-identified format. The study proposal was submitted to the Institutional Assessment Board office for review of human subject protections such as privacy and confidentiality, as well as the safeguarding of identifiable health information, in accordance with Roswell Park policy.
Data Availability Statement
The data and code used in this manuscript are available upon request.
Conflicts of Interest
The authors declare no conflict of interest.
Acknowledgments
We would like to acknowledge the Roswell Park Comprehensive Cancer Center pathology core and experimental tumor model resource for their assistance in preparing and reviewing the tissue microarrays for our analysis.