Biomarker Assessment of Homologous Recombination Deficiency in Epithelial Ovarian Cancer: Association With Progression-Free Survival After Surgery
Huan Yi1† 
Linhong Li3† 
Jimiao Huang1 
Zhiming Ma3 Hongping Li3 Jian Chen2 Xiangqin Zheng1 Jingjing Chen4 Haixin He2 
Jianrong Song1*- 1Department of Gynecology Oncology, Fujian Key Laboratory of Women and Children’s Critical Diseases Research [Fujian Maternity and Child Health Hospital(Fujian Women and Children’s Hospital)], Fujian Maternity and Child Health Hospital College of Clinical Medicine for Obstetrics and Gynecology and Pediatrics, Fujian Medical University, Fuzhou, China
- 2Department of Gynecology Oncology, Fujian Provincial Cancer Hospital, Fuzhou, China
- 3Research and Development Division, Oriomics Biotech Inc, Hangzhou, China
- 4Department of Obstetrics and Gynecology, Fuding General Hospital, Fuding, China
Identifying BRCA mutations and homologous recombination deficiency (HRD) is the key to choosing patients for poly (ADP-ribose) polymerase inhibitor (PARPi) therapy. At present, a large amount of research focuses on the application of HRD detection in ovarian cancer. However, few studies have discussed the relationship between HRD detection and postoperative survival in patients with epithelial ovarian cancer (EOC). This study included 38 consecutive patients with EOC who underwent cytoreduction surgery. Owing to tissue availability, only 29 patients underwent molecular profiling and survival analysis. Overall, 21 (72.4%) tumors had HRD scores of ≥42. Mutations in BRCA were observed in 5/29 (17.2%) patients. In this cohort, an HRD score of ≥42 was more common in serous ovarian tumors. We found no statistically significant association between homologous recombination repair (HRR) genes and HRD scores except for tumor protein P53 (TP53) mutation. We also found a strong positive association between HRD scores and chromosomal instability (CIN). In the survival analysis, an HRD score of >23 was correlated with better postoperative progression-free survival (pPFS). With increased depth of research, an appropriate HRD score threshold may serve as a prognostic tool and should be assessed in future studies to predict the clinical value of PARPi.
Introduction
Ovarian cancer is the third most frequent gynecologic cancer and the leading cause of gynecologic cancer death in China (Association, 2021; Sung et al., 2021). Due to the lack of effective early screening methods, most patients are diagnosed at a late stage. The 5-years survival rate is lower than 40% in China (Jiang et al., 2018). The World Health Organization (WHO) classifies ovarian cancers into those involving epithelial cells, germ cells, and mesenchyme (Kaku et al., 2003). Epithelial ovarian cancer (EOC) is the major type of ovarian malignancy, and its histologic classifications is as follows: high-grade serous (HGS) (up to 70%), low-grade serous, endometrioid, clear cell, and mucinous (Ledermann et al., 2016). The primary standard treatment approaches for EOC, such as debulking surgery and adjuvant platinum chemotherapy, can be used to achieve good initial response rates, but most ovarian cancer patients will eventually relapse (Deraco et al., 2012; Hiss, 2012; Zheng and Gao, 2012).
The homologous recombination (HR) pathway is an important DNA damage response mechanism that safeguards genome stability by repairing double-strand DNA breaks with high fidelity. Among the numerous albuminoidal factors of the HR pathway, BRCA proteins are essential for the complete HR repair response (Wang and Lindahl, 2016; Pastink et al., 2001; del Rivero and Kohn, 2017). Germline or somatic BRCA mutations and BRCA gene promoter methylation are the main causes of homologous recombination deficiency (HRD), but other genetic abnormalities in the homologous recombination repair (HRR) pathway may also contribute to HRD (Xu et al., 2010; Pilié et al., 2019; Kalachand et al., 2020). Furthermore, homologous recombination proficiency (HRP) can be restored by reverse mutation of various HRR pathway genes, such as BRCA1/2, PALB2, and RAD51C/D, indicating that HRD status is a complex and dynamic phenotype (Norquist et al., 2011; Goodall et al., 2017; Kondrashova et al., 2017; Weigelt et al., 2017; Lin et al., 2019). Due to the abnormal function of genes in the HR pathway, cells with HRD are more sensitive to platinum drugs and poly (ADP- ribose) polymerase inhibitors (PARPi) (Telli et al., 2016; Stronach et al., 2018; Stover et al., 2020). In recent years, multiple clinical studies have shown that PARPi treatment produces better progression-free survival (PFS) in patients with BRCA1/2 mutation (Mirza et al., 2016; Coleman et al., 2017; Pujade-Lauraine et al., 2017). Expanding the use of PARPi by identifying other HRD biomarkers is increasingly the focus of research. Currently, HRD tests can be divided into three categories: 1) HRR pathway-related genes, 2) genomic “scars” or signatures of mutations that represent patterns of genome instability and 3) functional assays that provide a real-time read out of HRD or HRP (Miller et al., 2020; Chiang et al., 2021).
An HRD score ≥42 is correlated with longer survival in EOC patients cured with platinum drugs or PARPi (Stronach et al., 2018; Chen et al., 2021). However, it is uncertain whether HRD status is associated with postoperative PFS (pPFS) in EOC. Therefore, we investigated the relationship between HRD status and prognosis in postoperative EOC patients, and we found that postoperative EOC can be divided into high and low-risk groups based on HRD status.
Methods
Study Populations
A consecutive series of 38 EOC patients who underwent surgical treatment from the Fujian Maternal and Child Health Hospital (from July 2017, to July 2021) were prospectively enrolled for HRD analyses. We included biopsy-proven epithelial ovarian cancer patients with sufficient tissue sample collection (five to ten 15 micron formalin-fixed paraffin embedded (FFPE) tissue sections) for molecular testing. Patients without sufficient tumor specimens and clinical information were excluded from this study. Chemotherapy regimens for patients receiving neoadjuvant chemotherapy (NACT) or primary cytoreductive surgery (PCS) included cisplatin, carboplatin, Nida platinum and paclitaxel. The study was permitted by the local Ethics Committees of the Fujian Maternal and Child Health Hospital, and informed consent was acquired from each patient. Further details of the enrolled patients are shown in Table 1.
TABLE 1. Demographic and clinical characteristics (n = 29).
Clinical Data Collection
Patients were followed up from July 2013 to December 2021 every 6 months. Clinical information was acquired from the internal database of the hospital. PFS was assessed based on Response Evaluation Criteria in Solid Tumours (RECIST) version 1.1 criteria, measured as the time from primary surgery to disease progression or death as a result of any cause, whichever was earlier, and ultrasonography, computed tomography (CT), and magnetic resonance imaging (MRI) were performed at baseline and after six to eight cycles of chemotherapy (Bogaerts et al., 2009).
Capture-Based Targeted DNA Sequencing
DNA from FFPE specimens was purified with the TIANamp FFPE DNA Kit (DP331, TIANGEN Biotechnology, China) according to the manufacturer’s instructions. The extracted DNA specimens were quantified with a Qubit dsDNA HS Assay Kit (Invitrogen) and cleaved into 150–200-bp fragments with Covaris M220 (Covaris, Brighton, United Kingdom). Qsep100 (Bioptic, Taiwan, China) was used for fragment quality control. The damaged DNA was repaired by a VAHTS® DNA Damage Repair Kit (Vazyme, Nanjing, China). A DNA fragment library was created by a nano DNA library preparation kit (Nanodigmbio, Nanjing, China), and then PCR amplification and purification were performed. Afterwards, the library was captured by a x-Gen® Hybridization and Wash Kit (IDT). The designed HRD panel contained 112 genes and 32,501 single nucleotide polymorphisms (HRDcare®, Oriomics Biotech Inc, Zhejiang, China). NadPrep NanoBlocker (Nanodigmbio, Nanjing, China) was then used to block the barcoded adapters. A Qubit 3.0 fluorimeter (Invitrogen) and Qsep100 were used for the final library quality control. Sequencing was performed on an Illumina NextSeq 550 Dx platform, and the average coverage depth of the captured area was 150×.
Analysis of Sequencing Data
Fastp 0.20.0 software was used to evaluate the quality of the raw sequencing data (Chen et al., 2018). The sequencing data were aligned to HG19 utilizing Burrows–Wheeler Aligner v.0.7.17. AMtools v.1.7 was utilized to transform SAM files to BAM files and to dispose the BAM files by chromosomal coordinates (Li et al., 2009). Copy number variants (CNVs) and single nucleotide variation (SNV) identification was performed using facets v.0.16.0 and VarDict v.1.8.2. Mutations were then filtered according to the following criteria (Association, 2021): variant allele frequency (VAF) greater than 1% with at least three supporting reads; and (Sung et al., 2021) ExAC, 1,000 Genomes or dbSNP databases of <0.05. The remaining variants were annotated with ANNOVAR(Ng and Henikoff, 2003). The HRD score, including assessments of loss of heterozygosity, telomeric allelic imbalance, and large-scale state transition, was calculated utilizing the scarHRD software, as previously described (Telli et al., 2016; Sztupinszki et al., 2018; (Timms et al., 2014).
Analysis of the Differentially Expressed Genes and Functional Enrichment Analyses
The Cancer Genome Atlas (TCGA) EOC data set was used to examine the differences in mRNA levels between the high (≥ 42) and low (< 42) HRD score groups using a cutoff padj value of < 0.05 and a |log2FoldChange| value of > 1. DESeq2 (version 1.24.0) was used to acquire differentially expressed genes (DEGs). Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses of the DEGs were conducted by clusterProfiler (version 3.12.0).
Statistical Methods
All statistical analyses were performed using the R software package (version 3.6.1). Kaplan-Meier analyses were used to evaluate the relationship between gene mutation status and PFS based on the survival package (version 3.1–11). The survminer package (version 0.4.6) was used to determine the optimal cutoff point for the HRD score. All p values were two-sided unless otherwise noted, and values <0.05 were considered statistically significant.
Results
Clinical Characteristics of the Patients
Thirty-eight patients with EOC were analyzed for HRD score and BRCA1 and BRCA2 mutation status. All patients underwent NACT or PCS, and some patients subsequently received interval cytoreductive surgery. The research predated PARPi approvals for first-line maintenance treatment; therefore, few patients chose this therapy. The median age of the EOC patients was 52 years (range 34–70 years). Most patients had stage III and IV disease (75.9%) and serous histology (75.8%). The median PFS after the primary surgery was 13 months (range 2–79 months). All patients had a Karnofsky performance scale (KPS) score of 90–100 and a Zubrod performance status (ZPS) score of 0–1 at baseline.
Characterization of HRD
The HRD score analysis was successful in 29/38 samples (76.3%). Nine patients were excluded because of their low tumor DNA content. Patient demographics and clinical data are shown in Table 1. For the 29 tumor tissue samples that met the experimental requirements, we calculated the HRD score for each sample, referring to the three chromosomal abnormalities previously approved by the FDA in Myriad’s myChoice HRD, which included LOH, TAI, and LST. HRD positivity was determined by either a test value of ≥42 or the presence of a BRCA pathogenic or likely pathogenic mutation. Overall, 21 (72.4%) tumor samples had HRD scores of ≥42, which is similar to a previous report (Chen et al., 2021). The proportion of HRD scores of ≥42 in our study was higher than that reported in Western countries (Patel et al., 2018). Mutations in BRCA were observed in 5/29 (17.2%) patients. An HRD score of ≥42 was more common in serous ovarian tumors in our cohort (p = 0.0037; Figure 1A). Patients with advanced ovarian tumors had a higher proportion of HRD scores of ≥42, although this result was not statistically significant (Figure 1B).
FIGURE 1. Distribution of HRD scores with different histologies (A) and tumor stages (B). ** indicates p < 0.01; ns indicates no significant difference.
HRR Gene Mutation and Copy Number Variation
Pathogenic mutations in BRCA were detected in five (5/29, 17.2%) samples and included four germline mutations and one somatic mutation (Table 2). We then evaluated the association between the HRD score and mutations in all HRR pathway genes, not just BRCA mutations. These ovarian tumors had universal Tumor Protein P53 (TP53) mutations and many other genetic mutations, most of which were HRR pathway gene mutations (Figure 2). We found no statistically significant associations between the whole mutation statuses of the genes of the HRR pathway and HRD scores. Nevertheless, individual analyses of these HRR mutations showed that TP53 mutation was associated with high HRD scores (Figure 3, p = 0.031), which has also been reported in colorectal and prostate cancer ((Smeby et al., 2020; Lotan et al., 2021).
TABLE 2. HRD testing characteristics in patients.
FIGURE 2. HRR gene mutation landscape and HRD scores.
FIGURE 3. Association of HRD scores (A), loss of heterozygosity (LOH) (B), telomeric allelic imbalance (TAI) (C) and large-scale state transition (LST) (D) with TP53 gene mutations.
Next, we analyzed the CNVs of HRR-related genes and identified the relationship between CNV and HRD scores. The higher the HRD score was, the more abnormal the copy number. (Supplementary Figures S1–S29). We also found a significant correlation between HRD scores and chromosomal instability (CIN) (Figure 4, p = 4.7 e-09, R = 0.86).
FIGURE 4. Relationship between HRD scores and chromosomal instability (CIN).
Survival
Of all patients who completed the HRD score assessment (n = 29), 22 patients had an HRD score of>23. Among the 24 patients without somatic or germline BRCA mutations who achieved HRD score detection, 17 had an HRD score greater than 23 (Figure 5).
FIGURE 5. HRD score values in ovarian cancer. Values of HRD scores above 23 (dashed line) were regarded as HRD-high, and values below 23 were regarded as HRD-low.
The result of a postoperative progression-free survival (pPFS) univariate analysis utilizing different HRD score thresholds is shown in Table 3. The median pPFS was 13 months (range 2–79 months). Although patients with an HRD score of ≥42 had improved pPFS, there were no significant differences between groups. Notably, we found that an HRD score of >23 was statistically correlated with better pPFS (HR 0.21, 95% CI 0.048–0.97; p = 0.046) (Table 3). Multivariate analysis had similar results. An HRD score of >23 (HR 0.019, 95% CI 0.00095–0.36; p = 0.008) was significantly associated with improved pPFS when adjusting for histology and TP53 mutation (Table 4).
TABLE 3. Postoperative progression-free survival in patients.
TABLE 4. Multivariate analysis for postoperative progression-free survival.
In addition, our results showed that patients with low HRD scores had significantly worse pPFS than those with high HRD scores, which was statistically significant by adjusting the HRD score threshold to 23 (p = 0.18, Figure 6A, p = 0.028, Figure 6B). In relapsed patients, a higher HRD score was associated with a shorter platinum-free interval (PFI) (Figure 7).
FIGURE 6. The association between diverse HRD score thresholds and postoperative progression-free survival: (A) 42; (B) 23.
FIGURE 7. Relationship between HRD scores and platinum-free interval (PFI).
Functional Enrichment Analysis of DEGs
A total of 474 DEGs were identified, 252 of which were upregulated and 222 of which were downregulated. To further explore the association of the HRD score and survival, GO and KEGG enrichment analyses were conducted. GO enrichment analysis showed that the DEGs between the high and low HRD score groups were remarkably enriched for substrate-specific channel activity, ligand-gated ion channel activity (molecular function [MF]), icosanoid metabolic process, unsaturated fatty acid metabolic process (biological process [BP]), and synaptic membrane and postsynaptic membrane components (cellular component [CC]). KEGG pathway analysis revealed that the DEGs were dominant in neuroactive ligand-receptor interaction and chemical carcinogenesis-receptor activation (Figure 8).
FIGURE 8. GO and pathway enrichment analysis of the DEGs: (A) Molecular Function; (B) Biological Process; (C) Cellular Component. (D) KEGG pathway enrichment analysis of the DEGs.
Discussion
By enabling an assessment of genomic instability and BRCA mutations, HRD scores are biomarkers of HR pathway defects, and former studies have shown that these biomarkers predict beneficial responses to DNA-damaging factors in certain tumor types (Telli et al., 2016; Pujade-Lauraine et al., 2017; Mirza et al., 2016). The frontline therapy that is standard for advanced ovarian cancer patients is paclitaxel/carboplatin combined treatment, and the FDA has already approved PARPi for frontline maintenance therapy. With the increasing clinical demand and interest in the assessment of tumor molecular features to treat EOC, we tried to determine the characteristics of HRD biomarkers and their correlation with survival outcomes.
First, the quality of the clinical samples is very important for HRD score detection. In a previous study (Miller et al., 2020), (How et al., 2021), approximately 80%–90% of the samples were completely tested for HRD. Failed detection was frequently due to the poor quality of tissue samples. The percentage of completed testing was lower (29/38, 76.3%) in this study. This result may be due to our relatively poor sample quality. At the same time, we found that high-grade serous ovarian cancer had a higher HRD score than those of other histological ovarian cancers. It has been reported that more than 45% of patients with HGS ovarian cancer may have HRD abnormalities, and the clinical and molecular characteristics of HGS ovarian cancer are associated with genome instability (Moschetta et al., 2016; Hodgson et al., 2018). Our study revealed that the Chinese population had a higher proportion of HRD scores of ≥42 than those of other populations. We suspected that this finding may be due to differences in ethnicity and the different subtypes of ovarian cancer in the study cohort. Furthermore, TP53 gene mutations have been detected in more than 80% of patients. An unexpected finding of our study was the positive association between the TP53 mutation status and the HRD score. RAD51 may be the connection between HRD and wild-type TP53 in EOC. RAD51 plays an important role in homologous recombination by interactingwith DNA, mediating homologous pairing and strand exchange and assisting in double-strand break repair (Arias ‐ Lopez et al., 2006). TP53 downregulates the expression of the crucial HR proteins, including RAD51, thereby inhibiting inapposite DNA repair (Bindra et al., 2004). Nevertheless, prevalent variations in TP53 in tumors lead to overexpression of RAD51 and increased resistance to DNA damage. This finding is consistent with our conclusion, which may interpret the high HRD score in the previously mentioned patients with TP53 variations.
Finally, we analyzed the relationship between the HRD score and PFS after surgery. Some researchers have found that adjusting the HRD threshold may increase the accuracy and the predictive ability of the score for assessing beneficial responses to PARPi or platinum drugs (Stronach et al., 2018; How et al., 2021; Li et al., 2020). In univariate analysis, HRD score thresholds of 42 and 23 were both correlated with better pPFS in EOC patients. Nevertheless, only the result obtained from patients with an HRD score of 23 approached statistical significance (p = 0.046). Multivariate analysis showed the same results (p = 0.008). Our exploratory research of diverse HRD score thresholds in postoperative EOC tumors showed that HRD status was significantly correlated with pPFS (p = 0.028) upon altering the threshold from 42 to 23. This finding indicated that EOC patients with low HRD scores had significantly worse pPFS than those with high HRD scores. Interestingly, we found that the HRD score was negatively correlated with the PFI in relapsed patients. Furthermore, we analyzed the DEGs between EOC patients with high and low HRD scores in TCGA and found that these genes were mainly enriched in GO analysis pathways related to ion channel activity and unsaturated fatty acid metabolic processes. These particular gene expression signatures may be responsible for the positive correlation between the HRD score and the postoperative survival in EOC patients. By performing HRD detection in tissue samples from patients with postoperative EOC, we effectively assessed the postoperative survival of patients. Patients with low HRD scores might have an earlier postoperative recurrence; therefore, postoperative monitoring should be emphasized.
Our study has several limitations. First, our study was conducted with a small specimen size and short follow-up time. However, patient sample collection and follow-up are ongoing. Soon, we will conduct a larger cohort study with a longer follow-up time to fully demonstrate the practical value of HRD detection in Chinese ovarian cancer patients. Second, the adjustment of threshold from 42 to 23 might increase the false positive rate of PARPi-sensitive patients, which may lead to over-treatment and medical resources waste. Finally, further research is required to evaluate the forecasting power of these two thresholds for PARPi treatment benefits.
Conclusion
In summary, EOC therapy is becoming increasingly personalized, with the increasing knowledge of the molecular characteristics guiding therapy stratification. The HRD scoring test offers a worthy adjunct to gathering evidence of HRD that may have been missed by only detecting somatic BRCA mutation. An HRD score threshold of 23 strongly correlated with an improvement in pPFS unlike the current threshold of 42. This new threshold may help better predict postoperative patient survival.
Data Availability Statement
The data supporting the results of this research can be obtained from the corresponding authors upon reasonable request.
Ethics Statement
The study was approved by local Ethics Committees of Fujian Maternal and Child Health Hospital and written informed consent was obtained from all patients.
Author Contributions
HY and LL made equal contributions to this work. JS, LL, HH, and HY conceived the study. JH, JC, XZ, and JC were responsible for specimen collection and storage. ZM and LL analysed the data. HY and LL composed the first draft of the manuscript. HL, HH, JS, and HY edited the manuscript. All authors approved the final manuscript.
Funding
This study was funded by Special Health Subsidy of Fujian Provincial Finance Department (2100206 50502 302), Health Care Youth and Middle-aged Training Project of Fujian Province (2021GGB014), National Science Foundation of Fujian Province (2020J01336, 2019J01512) and Youth National Science Foundation of China (NSFC, No. 81802604).
Conflict of Interest
LL, HL, and ZM were employed by Oriomics Biotech Inc.
The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The reviewer QX declared a shared affiliation with the author(s) JC and HH to the handling editor at the time of review.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Acknowledgments
The authors thank all patients and their families who participated in this research.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmolb.2022.906922/full#supplementary-material
References
Arias‐Lopez, C., Lazaro‐Trueba, I., Kerr, P., Lord, C. J., Dexter, T., Iravani, M., et al. (2006). p53 Modulates Homologous Recombination by Transcriptional Regulation of the RAD51 Gene. EMBO Rep. 7 (2), 219–224. doi:10.1038/sj.embor.7400587
Association, G. O. C. O. C. A. C. (2021). Guidelines for Diagnosis and Treatment of Ovarian Malignant Tumors (2021 Edition). China Oncol. 31 (06), 490–500. doi:10.19401/j.cnki.1007-3639.2021.06.07
Bindra, R. S., Schaffer, P. J., Meng, A., Woo, J., Måseide, K., Roth, M. E., et al. (2004). Down-regulation of Rad51 and Decreased Homologous Recombination in Hypoxic Cancer Cells. Mol. Cell Biol. 24 (19), 8504–8518. doi:10.1128/MCB.24.19.8504-8518.2004
Bogaerts, J., Ford, R., Sargent, D., Schwartz, L. H., Rubinstein, L., Lacombe, D., et al. (2009). Individual Patient Data Analysis to Assess Modifications to the RECIST Criteria. Eur. J. Cancer 45 (2), 248–260. doi:10.1016/j.ejca.2008.10.027
Chen, S., Zhou, Y., Chen, Y., and Gu, J. (2018). Fastp: an Ultra-fast All-In-One FASTQ Preprocessor. Bioinformatics 34 (17), i884–i890. doi:10.1093/bioinformatics/bty560
Chen, X., Ni, J., Xu, X., Guo, W., Cheng, X., Zhao, Q., et al. (2021). Homologous Recombination Deficiency Associated with Response to Poly (ADP-Ribose) Polymerase Inhibitors in Ovarian Cancer Patients: The First Real-word Data from China. Front. Oncol. 11, 746571. doi:10.3389/fonc.2021.746571
Chiang, Y.-C., Lin, P.-H., and Cheng, W.-F. (2021). Homologous Recombination Deficiency Assays in Epithelial Ovarian Cancer: Current Status and Future Direction. Front. Oncol. 11, 675972. doi:10.3389/fonc.2021.675972
Coleman, R. L., Oza, A. M., Lorusso, D., Aghajanian, C., Oaknin, A., Dean, A., et al. (2017). Rucaparib Maintenance Treatment for Recurrent Ovarian Carcinoma after Response to Platinum Therapy (ARIEL3): a Randomised, Double-Blind, Placebo-Controlled, Phase 3 Trial. Lancet 390 (10106), 1949–1961. doi:10.1016/S0140-6736(17)32440-6
del Rivero, J., and Kohn, E. C. (2017). PARP Inhibitors: The Cornerstone of DNA Repair-Targeted Therapies. Oncol. Willist. Park) 31 (4), 265–273.
Deraco, M., Virzì, S., Iusco, D., Puccio, F., Macrì, A., Famulari, C., et al. (2012). Secondary Cytoreductive Surgery and Hyperthermic Intraperitoneal Chemotherapy for Recurrent Epithelial Ovarian Cancer: a Multi-Institutional Study. BJOG Int. J. Obstetrics Gynaecol. 119 (7), 800–809. doi:10.1111/j.1471-0528.2011.03207.x
Goodall, J., Mateo, J., Yuan, W., Mossop, H., Porta, N., Miranda, S., et al. (2017). Circulating Cell-free DNA to Guide Prostate Cancer Treatment with PARP Inhibition. Cancer Discov. 7 (9), 1006–1017. doi:10.1158/2159-8290.CD-17-0261
Hiss, D. (2012). Optimizing Molecular-Targeted Therapies in Ovarian Cancer: the Renewed Surge of Interest in Ovarian Cancer Biomarkers and Cell Signaling Pathways. J. Oncol. 2012, 737981. doi:10.1155/2012/737981
Hodgson, D. R., Dougherty, B. A., Lai, Z., Fielding, A., Grinsted, L., Spencer, S., et al. (2018). Candidate Biomarkers of PARP Inhibitor Sensitivity in Ovarian Cancer beyond the BRCA Genes. Br. J. Cancer 119 (11), 1401–1409. doi:10.1038/s41416-018-0274-8
How, J. A., Jazaeri, A. A., Fellman, B., Daniels, M. S., Penn, S., Solimeno, C., et al. (2021). Modification of Homologous Recombination Deficiency Score Threshold and Association with Long-Term Survival in Epithelial Ovarian Cancer. Cancers 13 (5), 946. doi:10.3390/cancers13050946
Jiang, X., Tang, H., and Chen, T. (2018). Epidemiology of Gynecologic Cancers in China. J. Gynecol. Oncol. 29 (1), e7. doi:10.3802/jgo.2018.29.e7
Kaku, T., Ogawa, S., Kawano, Y., Ohishi, Y., Kobayashi, H., Hirakawa, T., et al. (2003). Histological Classification of Ovarian Cancer. Med. electron Microsc. 36 (1), 9–17. doi:10.1007/s007950300002
Kalachand, R. D., Stordal, B., Madden, S., Chandler, B., Cunningham, J., Goode, E. L., et al. (2020). BRCA1 Promoter Methylation and Clinical Outcomes in Ovarian Cancer: An Individual Patient Data Meta-Analysis. J. Natl. Cancer Inst. 112 (12), 1190–1203. doi:10.1093/jnci/djaa070
Kondrashova, O., Nguyen, M., Shield-Artin, K., Tinker, A. V., Teng, N. N. H., Harrell, M. I., et al. (2017). Secondary Somatic Mutations Restoring RAD51C and RAD51D Associated with Acquired Resistance to the PARP Inhibitor Rucaparib in High-Grade Ovarian Carcinoma. Cancer Discov. 7 (9), 984–998. doi:10.1158/2159-8290.CD-17-0419
Ledermann, J. A., Drew, Y., and Kristeleit, R. S. (2016). Homologous Recombination Deficiency and Ovarian Cancer. Eur. J. cancer 60, 49–58. doi:10.1016/j.ejca.2016.03.005
Li, H., Handsaker, B., Wysoker, A., Fennell, T., Ruan, J., Homer, N., et al. (2009). The Sequence Alignment/Map Format and SAMtools. Bioinformatics 25 (16), 2078–2079. doi:10.1093/bioinformatics/btp352
Li, H., Tu, J., Zhao, Z., Chen, L., Qu, Y., Li, H., et al. (2020). Molecular Signatures of BRCAness Analysis Identifies PARP Inhibitor Niraparib as a Novel Targeted Therapeutic Strategy for Soft Tissue Sarcomas. Theranostics 10 (21), 9477–9494. doi:10.7150/thno.45763
Lin, K. K., Harrell, M. I., Oza, A. M., Oaknin, A., Ray-Coquard, I., Tinker, A. V., et al. (2019). BRCA Reversion Mutations in Circulating Tumor DNA Predict Primary and Acquired Resistance to the PARP Inhibitor Rucaparib in High-Grade Ovarian Carcinoma. Cancer Discov. 9 (2), 210–219. doi:10.1158/2159-8290.CD-18-0715
Lotan, T. L., Kaur, H. B., Salles, D. C., Murali, S., Schaeffer, E. M., Lanchbury, J. S., et al. (2021). Homologous Recombination Deficiency (HRD) Score in Germline BRCA2- versus ATM-Altered Prostate Cancer. Mod. Pathol. 34 (6), 1185–1193. doi:10.1038/s41379-020-00731-4
Miller, R. E., Leary, A., Scott, C. L., Serra, V., Lord, C. J., Bowtell, D., et al. (2020). ESMO Recommendations on Predictive Biomarker Testing for Homologous Recombination Deficiency and PARP Inhibitor Benefit in Ovarian Cancer. Ann. Oncol. 31 (12), 1606–1622. doi:10.1016/j.annonc.2020.08.2102
Mirza, M. R., Monk, B. J., Herrstedt, J., Oza, A. M., Mahner, S., Redondo, A., et al. (2016). Niraparib Maintenance Therapy in Platinum-Sensitive, Recurrent Ovarian Cancer. N. Engl. J. Med. 375 (22), 2154–2164. doi:10.1056/NEJMoa1611310
Moschetta, M., George, A., Kaye, S. B., and Banerjee, S. (2016). BRCA Somatic Mutations and Epigenetic BRCA Modifications in Serous Ovarian Cancer. Ann. Oncol. 27 (8), 1449–1455. doi:10.1093/annonc/mdw142
Ng, P. C., and Henikoff, S. (2003). SIFT: Predicting Amino Acid Changes that Affect Protein Function. Nucleic Acids Res. 31 (13), 3812–3814. doi:10.1093/nar/gkg509
Norquist, B., Wurz, K. A., Pennil, C. C., Garcia, R., Gross, J., Sakai, W., et al. (2011). Secondary Somatic Mutations Restoring BRCA1/2 Predict Chemotherapy Resistance in Hereditary Ovarian Carcinomas. Jco 29 (22), 3008–3015. doi:10.1200/JCO.2010.34.2980
Pastink, A., Eeken, J. C. J., and Lohman, P. H. M. (2001). Genomic Integrity and the Repair of Double-Strand DNA Breaks. Mutat. Research/Fundamental Mol. Mech. Mutagen. 480-481, 37–50. doi:10.1016/s0027-5107(01)00167-1
Patel, J. N., Braicu, I., Timms, K. M., Solimeno, C., Tshiaba, P., Reid, J., et al. (2018). Characterisation of Homologous Recombination Deficiency in Paired Primary and Recurrent High-Grade Serous Ovarian Cancer. Br. J. Cancer 119 (9), 1060–1066. doi:10.1038/s41416-018-0268-6
Pilié, P. G., Gay, C. M., Byers, L. A., O'Connor, M. J., and Yap, T. A. (2019). PARP Inhibitors: Extending Benefit beyond BRCA-Mutant Cancers. Clin. Cancer Res. 25 (13), 3759–3771. doi:10.1158/1078-0432.CCR-18-0968
Pujade-Lauraine, E., Ledermann, J. A., Selle, F., Gebski, V., Penson, R. T., Oza, A. M., et al. (2017). Olaparib Tablets as Maintenance Therapy in Patients with Platinum-Sensitive, Relapsed Ovarian Cancer and a BRCA1/2 Mutation (SOLO2/ENGOT-Ov21): a Double-Blind, Randomised, Placebo-Controlled, Phase 3 Trial. Lancet Oncol. 18 (9), 1274–1284. doi:10.1016/S1470-2045(17)30469-2
Smeby, J., Kryeziu, K., Berg, K. C. G., Eilertsen, I. A., Eide, P. W., Johannessen, B., et al. (2020). Molecular Correlates of Sensitivity to PARP Inhibition beyond Homologous Recombination Deficiency in Pre-clinical Models of Colorectal Cancer Point to Wild-type TP53 Activity. EBioMedicine 59, 102923. doi:10.1016/j.ebiom.2020.102923
Stover, E. H., Fuh, K., Konstantinopoulos, P. A., Matulonis, U. A., and Liu, J. F. (2020). Clinical Assays for Assessment of Homologous Recombination DNA Repair Deficiency. Gynecol. Oncol. 159 (3), 887–898. doi:10.1016/j.ygyno.2020.09.029
Stronach, E. A., Paul, J., Timms, K. M., Hughes, E., Brown, K., Neff, C., et al. (2018). Biomarker Assessment of HR Deficiency, Tumor BRCA1/2 Mutations, and CCNE1 Copy Number in Ovarian Cancer: Associations with Clinical Outcome Following Platinum Monotherapy. Mol. Cancer Res. 16 (7), 1103–1111. doi:10.1158/1541-7786.MCR-18-0034
Sung, H., Ferlay, J., Siegel, R. L., Laversanne, M., Soerjomataram, I., Jemal, A., et al. (2021). Global Cancer Statistics 2020: GLOBOCAN Estimates of Incidence and Mortality Worldwide for 36 Cancers in 185 Countries. CA A Cancer J. Clin. 71 (3), 209–249. doi:10.3322/caac.21660
Sztupinszki, Z., Diossy, M., Krzystanek, M., Reiniger, L., Csabai, I., Favero, F., et al. (2018). Migrating the SNP Array-Based Homologous Recombination Deficiency Measures to Next Generation Sequencing Data of Breast Cancer. NPJ Breast Cancer 4, 16. doi:10.1038/s41523-018-0066-6
Telli, M. L., Timms, K. M., Reid, J., Hennessy, B., Mills, G. B., Jensen, K. C., et al. (2016). Homologous Recombination Deficiency (HRD) Score Predicts Response to Platinum-Containing Neoadjuvant Chemotherapy in Patients with Triple-Negative Breast Cancer. Clin. Cancer Res. 22 (15), 3764–3773. doi:10.1158/1078-0432.CCR-15-2477
Timms, K. M., Abkevich, V., Hughes, E., Neff, C., Reid, J., Morris, B., et al. (2014). Association of BRCA1/2defects with Genomic Scores Predictive of DNA Damage Repair Deficiency Among Breast Cancer Subtypes. Breast Cancer Res. 16 (6), 475. doi:10.1186/s13058-014-0475-x
Wang, J., and Lindahl, T. (2016). Maintenance of Genome Stability. Genomics, Proteomics Bioinforma. 14 (3), 119–121. doi:10.1016/j.gpb.2016.06.001
Weigelt, B., Comino-Méndez, I., de Bruijn, I., Tian, L., Meisel, J. L., García-Murillas, I., et al. (2017). Diverse BRCA1 and BRCA2 Reversion Mutations in Circulating Cell-free DNA of Therapy-Resistant Breast or Ovarian Cancer. Clin. Cancer Res. 23 (21), 6708–6720. doi:10.1158/1078-0432.CCR-17-0544
Xu, X., Gammon, M. D., Zhang, Y., Cho, Y. H., Wetmur, J. G., Bradshaw, P. T., et al. (2010). Gene Promoter Methylation Is Associated with Increased Mortality Among Women with Breast Cancer. Breast Cancer Res. Treat. 121 (3), 685–692. doi:10.1007/s10549-009-0628-2
Keywords: epithelial ovarian cancer, homologous recombination deficiency, homologous recombination deficiency score, chromosomal instability, postoperative progression-free survival
Citation: Yi H, Li L, Huang J, Ma Z, Li H, Chen J, Zheng X, Chen J, He H and Song J (2022) Biomarker Assessment of Homologous Recombination Deficiency in Epithelial Ovarian Cancer: Association With Progression-Free Survival After Surgery. Front. Mol. Biosci. 9:906922. doi: 10.3389/fmolb.2022.906922
Received: 29 March 2022; Accepted: 17 May 2022;
Published: 13 June 2022.
Edited by:
Teng Ma, Capital Medical University, ChinaReviewed by:
Sravya Palavalasa, University of Michigan, United StatesYadong Yang, Biochain (Beijing) Science & Technology Inc., China
Qin Xu, Fujian Provincial Cancer Hospital, China
Copyright © 2022 Yi, Li, Huang, Ma, Li, Chen, Zheng, Chen, He and Song. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Jianrong Song, sjr0318@163.com
†These authors have contributed equally to this work and share first authorship