Pan-cancer analysis of genomic scar patterns caused by homologous repair deficiency (HRD)

Prediction of homologous recombination deficiency from routine histology with attention-based multiple instance learning in nine different tumor types

BACKGROUND

Homologous recombination deficiency (HRD) is recognized as a pan-cancer predictive biomarker that potentially indicates who could benefit from treatment with PARP inhibitors (PARPi). Despite its clinical significance, HRD testing is highly complex. Here, we investigated in a proof-of-concept study whether Deep Learning (DL) can predict HRD status solely based on routine hematoxylin & eosin (H&E) histology images across nine different cancer types.

METHODS

We developed a deep learning pipeline with attention-weighted multiple instance learning (attMIL) to predict HRD status from histology images. As part of our approach, we calculated a genomic scar HRD score by combining loss of heterozygosity (LOH), telomeric allelic imbalance (TAI), and large-scale state transitions (LST) from whole genome sequencing (WGS) data of n = 5209 patients across two independent cohorts. The model's effectiveness was evaluated using the area under the receiver operating characteristic curve (AUROC), focusing on its accuracy in predicting genomic HRD against a clinically recognized cutoff value.

RESULTS

Our study demonstrated the predictability of genomic HRD status in endometrial, pancreatic, and lung cancers reaching cross-validated AUROCs of 0.79, 0.58, and 0.66, respectively. These predictions generalized well to an external cohort, with AUROCs of 0.93, 0.81, and 0.73. Moreover, a breast cancer-trained image-based HRD classifier yielded an AUROC of 0.78 in the internal validation cohort and was able to predict HRD in endometrial, prostate, and pancreatic cancer with AUROCs of 0.87, 0.84, and 0.67, indicating that a shared HRD-like phenotype occurs across these tumor entities.

CONCLUSIONS

This study establishes that HRD can be directly predicted from H&E slides using attMIL, demonstrating its applicability across nine different tumor types.

Genomic Insights Into High-Grade Infarct-Associated Bone Sarcomas

Sarcomas rarely develop in bones previously compromised by infarcts. These infarct-associated sarcomas often present as undifferentiated pleomorphic sarcomas (UPS), and their genetic characteristics are poorly understood. High-grade UPS of bone are typically treated with a combination of surgery and chemotherapy, similar to osteosarcoma. We conducted a detailed clinicopathologic and genomic analysis of 6 cases of intraosseous sarcomas arising from histologically and radiographically confirmed bone infarcts. We analyzed 523 genes for sequence-level mutations using next-generation sequencing with the TruSight Oncology 500 panel and utilized whole-genome single nucleotide polymorphism Microarray (OncoScan CNV) to detect copy number alterations and loss of heterozygosity (LOH). Genomic instability was assessed through homologous recombination deficiency (HRD) metrics, incorporating LOH, telomeric allelic imbalance, and large-scale state transitions. Fluorescence in situ hybridization and immunohistochemistry validated the findings. The cohort included 3 men and 3 women, with a median age of 70 years, and tumors located in the femur and tibia. Five of the 6 patients developed distant metastases. Treatment involved surgery and chemotherapy or immune checkpoint inhibitors. Genomic analysis revealed significant complexity and high HRD scores, ranging from 32 to 57 (with a cutoff of 32). Chromosome 12 alterations, including segmental amplification or chromothripsis, were observed in 4 cases. Notably, MDM2 amplification, confirmed by fluorescence in situ hybridization, was detected in 2 cases. Homozygous deletion of CDKN2A/B was observed in all six cases. Tumor mutational burden levels ranged from 2.4 to 7.9 mutations per megabase. Notable pathogenic mutations included H3-3A mutations (p.G35R and p.G35W), and mutations in HRAS, DNMT3A, NF2, PIK3CA, POLE, and TP53, each in one case. These results suggest that high-grade infarct-associated sarcomas of bone, whereas sharing high levels of structural variations with osteosarcoma, may exhibit potentially less frequent TP53 mutations and more common CDKN2A/B deletions. This points to the possibility that the mutation spectrum and disrupted pathways could be distinct from conventional osteosarcoma.

Targeting DNA Damage Response Deficiency in Thoracic Cancers

Thoracic cancers comprise non-small cell lung cancers (NSCLCs), small cell lung cancers (SCLCs) and malignant pleural mesotheliomas (MPM). Collectively, they account for the highest rate of death from malignancy worldwide. Genomic instability is a universal feature of cancer, which fuels mutations and tumour evolution. Deficiencies in DNA damage response (DDR) genes amplify genomic instability. Homologous recombination deficiency (HRD), resulting from BRCA1/BRCA2 inactivation, is exploited for therapeutic synthetic lethality with poly-ADP ribose polymerase (PARP) inhibitors in breast and ovarian cancers, as well as in prostate and pancreatic cancers. However, DDR deficiency and its therapeutic implications are less well established in thoracic cancers. Emerging evidence suggests that a subset of thoracic cancers may harbour DDR deficiency and may, thus, be effectively targeted with DDR agents. Here, we review the current evidence surrounding DDR in thoracic cancers and discuss the challenges and promise for achieving clinical benefit with such therapeutics.

Homologous recombination deficiency (HRD) is associated with better prognosis and possibly causes a non-inflamed tumour microenvironment in nasopharyngeal carcinoma

Homologous recombination deficiency (HRD) score is a reliable indicator of genomic instability. The significance of HRD in nasopharyngeal carcinoma (NPC), particularly its influence on prognosis and the immune microenvironment, has yet to be adequately explored. Understanding HRD status comprehensively can offer valuable insights for guiding precision treatment. We utilised three cohorts to investigate HRD status in NPC: the Zhujiang cohort from local collection and the Hong Kong (SRA288429) and Singapore (SRP035573) cohorts from public datasets. The GATK (genome analysis toolkit) best practice process was employed to investigate germline and somatic BRCA1/2 mutations and various bioinformatics tools and algorithms to examine the association between HRD status and clinical molecular characteristics. We found that individuals with a negative HRD status (no-HRD) exhibited a higher risk of recurrence [hazard ratio (HR), 1.43; 95% confidence interval (CI), 2.03-333.76; p = 0.012] in the Zhujiang cohort, whereas, in the Singapore cohort, they experienced a higher risk of mortality (HR, 26.04; 95% CI, 1.43-34.21; p = 0.016) compared with those in the HRD group. In vitro experiments demonstrated that NPC cells with BRCA1 knockdown exhibit heightened sensitivity to chemoradiotherapy. Furthermore, the HRD group showed significantly higher tumour mutational burden and tumour neoantigen burden levels than the no-HRD group. Immune infiltration analysis indicated that HRD tissues tend to have a non-inflamed tumour microenvironment. In conclusion, patients with HRD exhibit a comparatively favourable prognosis in NPC, possibly associated with a non-inflammatory immune microenvironment. These findings have positive implications for treatment stratification, enabling the selection of more precise and effective therapeutic approaches and aiding in the prediction of treatment response and prognosis to a certain extent.

Clinical and molecular significance of homologous recombination deficiency positive non-small cell lung cancer in Chinese population: An integrated genomic and transcriptional analysis

Objective

The clinical significance of homologous recombination deficiency (HRD) in breast cancer, ovarian cancer, and prostate cancer has been established, but the value of HRD in non-small cell lung cancer (NSCLC) has not been fully investigated. This study aimed to systematically analyze the HRD status of untreated NSCLC and its relationship with patient prognosis to further guide clinical care.

Methods

A total of 355 treatment-naïve NSCLC patients were retrospectively enrolled. HRD status was assessed using the AmoyDx Genomic Scar Score (GSS), with a score of ≥50 considered HRD-positive. Genomic, transcriptomic, tumor microenvironmental characteristics and prognosis between HRD-positive and HRD-negative patients were analyzed.

Results

Of the patients, 25.1% (89/355) were HRD-positive. Compared to HRD-negative patients, HRD-positive patients had more somatic pathogenic homologous recombination repair (HRR) mutations, higher tumor mutation burden (TMB) (P<0.001), and fewer driver gene mutations (P<0.001). Furthermore, HRD-positive NSCLC had more amplifications in PI3K pathway and cell cycle genes, MET and MYC in epidermal growth factor receptor (EGFR)/anaplastic lymphoma kinase (ALK) mutant NSCLC, and more PIK3CA and AURKA in EGFR/ALK wild-type NSCLC. HRD-positive NSCLC displayed higher tumor proliferation and immunosuppression activity. HRD-negative NSCLC showed activated signatures of major histocompatibility complex (MHC)-II, interferon (IFN)-γ and effector memory CD8+ T cells. HRD-positive patients had a worse prognosis and shorter progression-free survival (PFS) to targeted therapy (first- and third-generation EGFR-TKIs) (P=0.042). Additionally, HRD-positive, EGFR/ALK wild-type patients showed a numerically lower response to platinum-free immunotherapy regimens.

Conclusions

Unique genomic and transcriptional characteristics were found in HRD-positive NSCLC. Poor prognosis and poor response to EGFR-TKIs and immunotherapy were observed in HRD-positive NSCLC. This study highlights potential actionable alterations in HRD-positive NSCLC, suggesting possible combinational therapeutic strategies for these patients.

Identification of CXCL13 as a Promising Biomarker for Immune Checkpoint Blockade Therapy and PARP Inhibitor Therapy in Ovarian Cancer

Ovarian cancer has poor response rates to immune checkpoint blockade (ICB) therapy, despite the use of genomic sequencing to identify molecular targets. Homologous recombination deficiency (HRD) is a conventional indicator of genomic instability (GI) and has been used as a marker for targeted therapies. Indicators reflecting HRD status have shown potential in predicting the efficacy of ICB treatment. Public databases, including TCGA, ICGC, and GEO, were used to obtain data. HRD scores, neoantigen load, and TMB were obtained from the TCGA cohort. Candidate biomarkers were validated in multiple databases, such as the Imvigor210 immunotherapy cohort and the open-source single-cell sequencing database. Immunohistochemistry was performed to further validate the results in independent cohorts. CXCL10, CXCL11, and CXCL13 were found to be significantly upregulated in HRD tumors and exhibited prognostic value. A comprehensive analysis of the tumor immune microenvironment (TIME) revealed that CXCL13 expression positively correlated with neoantigen load and immune cell infiltration. In addition, single-cell sequencing data and clinical trial results supported the utility of CXCL13 as a biomarker for ICB therapy. Not only does CXCL13 serve as a biomarker reflecting HRD status, but it also introduces a potentially novel perspective on prognostic biomarkers for ICB in ovarian cancer.

Leveraging Off-Target Reads in Panel Sequencing for Homologous Recombination Repair Deficiency Screening in Tumor

Targeted tumor only sequencing has become a standard practice in cancer diagnostics. This study aims to develop an approach for robust copy number variant calling in tumor samples using only off-target region (OTR) reads. We also established a clinical use case for homologous recombination deficiency (HRD) score estimation (HRDest) using the sum of telomeric-allelic imbalance and large-scale state transition scores without the need for loss of heterozygosity information. A strong correlation was found between HRD score and the sum of telomeric-allelic imbalance + large-scale state transition in The Cancer Genome Atlas cohort (ρ = 0.99, P < 2.2 × 10-16) and in a clinical in-house cohort of 34 tumors (ρ = 0.9, P = 5.1 × 10-13) comparing whole-exome sequencing and targeted sequencing data. HRDest scores from 1086 clinical cases were compared with The Cancer Genome Atlas data set. There were no significant differences in HRD score distribution within the analyzed tumor types. As a control, commercially available HRD standards were also sequenced, and the HRDest scores obtained from the OTR reads were well within the HRD reference range provided by the manufacturer. In conclusion, OTR reads of tumor-only panel sequencing can be used to determine genome-wide copy number variant profiles and to approximate HRD scores.

High-throughput molecular assays for inclusion in personalised oncology trials - State-of-the-art and beyond

In the last decades, the development of high-throughput molecular assays has revolutionised cancer diagnostics, paving the way for the concept of personalised cancer medicine. This progress has been driven by the introduction of such technologies through biomarker-driven oncology trials. In this review, strengths and limitations of various state-of-the-art sequencing technologies, including gene panel sequencing (DNA and RNA), whole-exome/whole-genome sequencing and whole-transcriptome sequencing, are explored, focusing on their ability to identify clinically relevant biomarkers with diagnostic, prognostic and/or predictive impact. This includes the need to assess complex biomarkers, for example microsatellite instability, tumour mutation burden and homologous recombination deficiency, to identify patients suitable for specific therapies, including immunotherapy. Furthermore, the crucial role of biomarker analysis and multidisciplinary molecular tumour boards in selecting patients for trial inclusion is discussed in relation to various trial concepts, including drug repurposing. Recognising that today's exploratory techniques will evolve into tomorrow's routine diagnostics and clinical study inclusion assays, the importance of emerging technologies for multimodal diagnostics, such as proteomics and in vivo drug sensitivity testing, is also discussed. In addition, key regulatory aspects and the importance of patient engagement in all phases of a clinical trial are described. Finally, we propose a set of recommendations for consideration when planning a new precision cancer medicine trial.

A Systematic Review of the Prevalence of Germline BRCA mutations in North Asia Breast Cancer Patients

OBJECTIVE

The BRCA1/2 mutation status testing is the global standard of care for breast cancer patients with a family history of cancer. BRCA1/2 mutations are known to be ethno-specific. For some ethnic groups of the Northern Asia (Buryats, Yakuts, Altaians, Tuvans, Khakasses, etc.) the founder mutations in the BRCA1/2 genes have not been revealed. This systematic review was conducted to assess the prevalence of BRCA1/2 mutation in breast cancer patients inhabiting Eastern Europe and Northern Asia (or Siberia).

METHODS

A total of 23,561 studies published between 2014 and 2024 were analyzed, of which 55 were included in the review. The literature search was conducted using RusMed, Cyberleninka, Google Scholar, eLibrary, NCBI databases (n=5) and conference papers.

RESULTS

The founder mutations (c.5266dupC and/or c.181T>G) of BRCA1 gene that were frequently observed in the Slav peoples were also identified in Chechens, Armenians, Bashkirs, Ukrainians, Mordovians, Mari, Kabardians, Tatars, Uzbeks, Kyrgyz, Ossetians, Khanty indigenous peoples and Adygs. For Chechens, Kabardians, Ingush, Buryats, Khakasses, Sakha, Tuvans and Armenians, rare pathogenic variants of the BRCA1/2, ATM, СНЕК2, BRIP1, NBN, PTEN, TP53, PMS1, XPA, LGR4, BRWD1 and PALB2 genes were found. No data are available about the frequency of pathogenic BRCA1/2 mutations for ethnic groups, such as the Udmurts, Komi, Tajiks, Tabasarans, and Nogais indigenous people.

CONCLUSION

This is the first systematic review that provides the spectrum of BRCA mutations in ethnic groups of breast cancer patients inhabiting Eastern Europe and Northern Asia. It has been shown that the mutations are ethnospecific (varied widely within groups) and not all groups are equally well studied. Further studies on the ethnic specificity of BRCA gene mutations are required.

High Concordance of Different Assays in the Determination of Homologous Recombination Deficiency-Associated Genomic Instability in Ovarian Cancer

PURPOSE

Poly(ADP-ribose) polymerase inhibitors (PARPi) have shown promising clinical results in the treatment of ovarian cancer. Analysis of biomarker subgroups consistently revealed higher benefits for patients with homologous recombination deficiency (HRD). The test that is most often used for the detection of HRD in clinical studies is the Myriad myChoice assay. However, other assays can also be used to assess biomarkers, which are indicative of HRD, genomic instability (GI), and BRCA1/2 mutation status. Many of these assays have high potential to be broadly applied in clinical routine diagnostics in a time-effective decentralized manner. Here, we compare the performance of a multitude of alternative assays in comparison with Myriad myChoice in high-grade serous ovarian cancer (HGSOC).

METHODS

DNA from HGSOC samples was extracted from formalin-fixed paraffin-embedded tissue blocks of cases previously run with the Myriad myChoice assay, and GI was measured by multiple molecular assays (CytoSNP, AmoyDx, Illumina TSO500 HRD, OncoScan, NOGGO GISv1, QIAseq HRD Panel and whole genome sequencing), applying different bioinformatics algorithms.

RESULTS

Application of different assays to assess GI, including Myriad myChoice, revealed high concordance of the generated scores ranging from very substantial to nearly perfect fit, depending on the assay and bioinformatics pipelines applied. Interlaboratory comparison of assays also showed high concordance of GI scores.

CONCLUSION

Assays for GI assessment not only show a high concordance with each other but also in correlation with Myriad myChoice. Thus, almost all of the assays included here can be used effectively to assess HRD-associated GI in the clinical setting. This is important as PARPi treatment on the basis of these tests is compliant with European Medicines Agency approvals, which are methodologically not test-bound.

A deep-learning-based genomic status estimating framework for homologous recombination deficiency detection from low-pass whole genome sequencing

Genome-wide sequencing allows for prediction of clinical treatment responses and outcomes by estimating genomic status. Here, we developed Genomic Status scan (GSscan), a long short-term memory (LSTM)-based deep-learning framework, which utilizes low-pass whole genome sequencing (WGS) data to capture genomic instability-related features. In this study, GSscan directly surveys homologous recombination deficiency (HRD) status independent of other existing biomarkers. In breast cancer, GSscan achieved an AUC of 0.980 in simulated low-pass WGS data, and obtained a higher HRD risk score in clinical BRCA-deficient breast cancer samples (p = 1.3 × 10-4, compared with BRCA-intact samples). In ovarian cancer, GSscan obtained higher HRD risk scores in BRCA-deficient samples in both simulated data and clinical samples (p = 2.3 × 10-5 and p = 0.039, respectively, compared with BRCA-intact samples). Moreover, HRD-positive patients predicted by GSscan showed longer progression-free intervals in TCGA datasets (p = 0.0011) treated with platinum-based adjuvant chemotherapy, outperforming existing low-pass WGS-based methods. Furthermore, GSscan can accurately predict HRD status using only 1 ng of input DNA and a minimum sequencing coverage of 0.02 × , providing a reliable, accessible, and cost-effective approach. In summary, GSscan effectively and accurately detected HRD status, and provide a broadly applicable framework for disease diagnosis and selecting appropriate disease treatment.

Genomic Features of Homologous Recombination Deficiency in Breast Cancer: Impact on Testing and Immunotherapy

Genomic instability is one of the well-established hallmarks of cancer. The homologous recombination repair (HRR) pathway plays a critical role in correcting the double-stranded breaks (DSB) due to DNA damage in human cells. Traditionally, the BRCA1/2 genes in the HRR pathway have been tested for their association with breast cancer. However, defects in the HRR pathway (HRD, also termed 'BRCAness'), which has up to 50 genes, have been shown to be involved in tumorigenesis and treatment susceptibility to poly-ADP ribose polymerase inhibitors (PARPis), platinum-based chemotherapy, and immune checkpoint inhibitors (ICIs). A reliable consensus on HRD scores is yet to be established. Emerging evidence suggests that only a subset of breast cancer patients benefit from ICI-based immunotherapy. Currently, albeit with limitations, the expression of programmed death-ligand 1 (PDL1) and tumor mutational burden (TMB) are utilized as biomarkers to predict the favorable outcomes of ICI therapy in breast cancer patients. Preclinical studies demonstrate an interplay between the HRR pathway and PDL1 expression. In this review, we outline the current understanding of the role of HRD in genomic instability leading to breast tumorigenesis and delineate outcomes from various clinical trials. Furthermore, we discuss potential strategies for combining HRD-targeted therapy with immunotherapy to achieve the best healthcare outcomes in breast cancer patients.

cGAS-STING pathway expression correlates with genomic instability and immune cell infiltration in breast cancer

Genomic instability, as caused by oncogene-induced replication stress, can lead to the activation of inflammatory signaling, involving the cGAS-STING and JAK-STAT pathways. Inflammatory signaling has been associated with pro-tumorigenic features, but also with favorable response to treatment, including to immune checkpoint inhibition. In this study, we aim to explore relations between inflammatory signaling, markers of replication stress, and immune cell infiltration in breast cancer. Expression levels of cGAS-STING signaling components (STING, phospho-TBK1, and phospho-STAT1), replication stress markers (γH2AX and pRPA), replication stress-related proto-oncogenes (Cyclin E1 and c-Myc) and immune cell markers (CD20, CD4, and CD57) are determined immunohistochemically on primary breast cancer samples (n = 380). RNA-sequencing data from TCGA (n = 1082) and METABRIC (n = 1904) are used to calculate cGAS-STING scores. pTBK1, pSTAT1 expression and cGAS-STING pathway scores are all increased in triple-negative breast cancers compared to other subtypes. Expression of γH2AX, pRPA, Cyclin E1, c-Myc, and immune cell infiltration positively correlate with p-STAT1 expression (P < 0.001). Additionally, we observe significant positive associations between expression of pTBK1 and γH2AX, pRPA, c-Myc, and number of CD4+ cells and CD20+ cells. Also, cGAS-STING scores are correlated with genomic instability metrics, such as homologous recombination deficiency (P < 0.001) and tumor mutational burden (P < 0.01). Moreover, data from the I-SPY2 clinical trial (n = 71) confirms that higher cGAS-STING scores are observed in breast cancer patients who responded to immunotherapy combined with chemotherapy. In conclusion, the cGAS-STING pathway is highly expressed in TNBCs and is correlated with genomic instability and immune cell infiltration.

Challenges and Future Directions in the Management of Tumor Mutational Burden-High (TMB-H) Advanced Solid Malignancies

A standardized assessment of Tumor Mutational Burden (TMB) poses challenges across diverse tumor histologies, treatment modalities, and testing platforms, requiring careful consideration to ensure consistency and reproducibility. Despite clinical trials demonstrating favorable responses to immune checkpoint inhibitors (ICIs), not all patients with elevated TMB exhibit benefits, and certain tumors with a normal TMB may respond to ICIs. Therefore, a comprehensive understanding of the intricate interplay between TMB and the tumor microenvironment, as well as genomic features, is crucial to refine its predictive value. Bioinformatics advancements hold potential to improve the precision and cost-effectiveness of TMB assessments, addressing existing challenges. Similarly, integrating TMB with other biomarkers and employing comprehensive, multiomics approaches could further enhance its predictive value. Ongoing collaborative endeavors in research, standardization, and clinical validation are pivotal in harnessing the full potential of TMB as a biomarker in the clinic settings.

Low microsatellite instability: A distinct instability type in gastric cancer?

PURPOSE

We recently showed that low microsatellite instability (MSI-L) is associated with a good response to platinum/5-fluorouracil (5-FU) neoadjuvant chemotherapy (CTx) in gastric cancer. The purpose of this study was to characterize the instability pattern and to investigate an association of MSI-L tumors with mutations in genes of DNA repair pathways and with total tumor mutation burden (TMB).

METHODS

MSI patterns were compared between 67 MSI high (-H) and 35 MSI-L tumors. Whole-exome sequencing was performed in 34 microsatellite stable (MSS) and 20 MSI-L tumors after or without neoadjuvant CTx.

RESULTS

Of the 35 MSI-L tumors, 33 tumors had instability at a dinucleotide repeat marker. In the homologous recombination (HR) pathway, 10 of the 34 (29%) MSS and 10 of the 20 (50%) MSI-L tumors showed variants (p = 0.154). In the DNA damage tolerance pathway, 6 of the 34 (18%) MSS and 7 of the 20 (35%) MSI-L tumors had variants (p = 0.194). The HR deficiency score was similar in both tumor groups. TMB was significantly higher in MSI-L compared to MSS tumors after CTx (p = 0.046). In the MSS and MSI-L tumors without CTx no difference was observed (p = 1.00).

CONCLUSION

MSI-L due to instability at dinucleotide repeat markers was associated with increased TMB after neoadjuvant CTx treatment, indicating sensitivity to platinum/5-FU CTx. If confirmed in further studies, this could contribute to refined chemotherapeutic options including immune-based strategies for GC patients with MSI-L tumors.

Combination of genomic instability score and TP53 status for prognosis prediction in lung adenocarcinoma

The genomic instability (GI) /homologous recombination deficiency (HRD) score, calculated as the sum of the events of loss of heterozygosity (LOH), large-scale state transition (LST) and telomere allele imbalance (TAI), is used to guide the choice of treatment in several cancers, but its relationship with genomic features, clinicopathological characteristics and prognosis in lung cancer is poorly understood, which could lead to population bias in prospective studies. We retrospectively analyzed 1011 lung cancer patients whose tumor samples were successfully profiled by high-throughput sequencing panel including GI/HRD score. Alterations of many cancer suppressor genes were associated with higher GI/HRD scores, biallelic inactivation of TP53 was correlated with a high GI/HRD score. A combination of two gene alterations exhibited a higher GI/HRD scores than single gene alterations. The GI/HRD score was associated with advanced stages in lung adenocarcinoma but not in lung squamous cell carcinoma. Furthermore, patients with higher GI/HRD scores had significantly shorter overall survival and progression-free survival than patients with lower GI/HRD scores. Finally, patients with a combination of a higher GI/HRD scores and TP53 alteration exhibited an extremely poor prognosis compared with patients with a lower GI/HRD scores and wild-type TP53 (overall survival, training cohort, hazard ratio (HR) = 8.56, P < 0.001; validation cohort, HR = 6.47, P < 0.001; progression-free survival, HR = 4.76, P < 0.001). Our study revealed the prognostic value of the GI/HRD score in lung adenocarcinoma, but not for all lung cancer. Moreover, the combination of the GI/HRD score and TP53 status could be a promising strategy to predict the prognosis of patients with lung adenocarcinoma.

DNA repair biomarkers to guide usage of combined PARP inhibitors and chemotherapy: A meta-analysis and systematic review

PURPOSE

The addition of PARP inhibitors to chemotherapy has been assessed in > 80 clinical trials across multiple malignancies, on the premise that PARP inhibitors will increase chemotherapy effectiveness regardless of whether cancers have underlying disruption of DNA repair pathways. Consequently, the majority of combination therapy trials have been performed on patients without biomarker selection, despite the use of homologous recombination deficiency to dictate use of PARP inhibitors in the maintenance setting. An unresolved question is whether biomarkers are needed to identify patients who respond to combination PARP inhibitors and chemotherapy.

METHODS

A systematic literature review identified studies using PARP inhibitors in combination with chemotherapy versus chemotherapy alone, where the study included a biomarker of DNA repair function (BRCA1, BRCA2, homologous recombination deficiency test, ATM, ERCC1, SLFN11). Hazard ratios (HR) were pooled in a meta-analysis using generic inverse-variance, and fixed or random effects modelling. Subgroup analyses were conducted on biomarker selection and type of malignancy.

RESULTS

Nine studies comprising 2547 patients met the inclusion criteria. Progression-free survival (PFS) was significantly better in patients with a DNA repair biomarker (HR: 0.57, 95% CI: 0.48-0.68, p < 0.00001), but there was no benefit in patients who lacked a biomarker (HR: 0.94, 95% CI: 0.82-1.08, p = 0.38). Subgroup analysis showed that BRCA status and SLFN11 biomarkers could predict benefit, and biomarker-driven benefit occurred in ovarian, breast and small cell lung cancers. The addition of PARP inhibitors to chemotherapy was associated with increased grade 3/4 side effects, and particularly neutropenia.

CONCLUSIONS

Combination therapy only improves PFS in patients with identifiable DNA repair biomarkers. This indicates that PARP inhibitors do not sensitise patients to chemotherapy treatment, except where their cancer has a homologous recombination defect, or an alternative biomarker of altered DNA repair. While effective in patients with DNA repair biomarkers, there is a risk of high-grade haematological side-effects with the use of combination therapy. Thus, the benefit in PFS from combination therapy must be weighed against potential adverse effects, as individual arms of treatment can also confer benefit.

Integrative bioinformatics analysis of WDHD1: a potential biomarker for pan-cancer prognosis, diagnosis, and immunotherapy

BACKGROUND

Although WD repeat and high-mobility group box DNA binding protein 1 (WDHD1) played an essential role in DNA replication, chromosome stability, and DNA damage repair, the panoramic picture of WDHD1 in human tumors remains unclear. Hence, this study aims to comprehensively characterize WDHD1 across 33 human cancers.

METHODS

Based on publicly available databases such as TCGA, GTEx, and HPA, we used a bioinformatics approach to systematically explore the genomic features and biological functions of WDHD1 in pan-cancer.

RESULTS

WDHD1 mRNA levels were significantly increased in more than 20 types of tumor tissues. Elevated WDHD1 expression was associated with significantly shorter overall survival (OS) in 10 tumors. Furthermore, in uterine corpus endometrial carcinoma (UCEC) and liver hepatocellular carcinoma (LIHC), WDHD1 expression was significantly associated with higher histological grades and pathological stages. In addition, WDHD1 had a high diagnostic value among 16 tumors (area under the ROC curve [AUC] > 0.9). Functional enrichment analyses suggested that WDHD1 probably participated in many oncogenic pathways such as E2F and MYC targets (false discovery rate [FDR] < 0.05), and it was involved in the processes of DNA replication and DNA damage repair (p.adjust < 0.05). WDHD1 expression also correlated with the half-maximal inhibitory concentrations (IC50) of rapamycin (4 out of 10 cancers) and paclitaxel (10 out of 10 cancers). Overall, WDHD1 was negatively associated with immune cell infiltration and might promote tumor immune escape. Our analysis of genomic alterations suggested that WDHD1 was altered in 1.5% of pan-cancer cohorts and the "mutation" was the predominant type of alteration. Finally, through correlation analysis, we found that WDHD1 might be closely associated with tumor heterogeneity, tumor stemness, mismatch repair (MMR), and RNA methylation modification, which were all processes associated with the tumor progression.

CONCLUSIONS

Our pan-cancer analysis of WDHD1 provides valuable insights into the genomic characterization and biological functions of WDHD1 in human cancers and offers some theoretical support for the future use of WDHD1-targeted therapies, immunotherapies, and chemotherapeutic combinations for the management of tumors.

A novel marker integrating multiple genetic alterations better predicts platinum sensitivity in ovarian cancer than HRD score

Introduction: Platinum-based chemotherapy is the first-line treatment strategy for ovarian cancer patients. The dismal prognosis of ovarian cancer was shown to be stringently associated with the heterogeneity of tumor cells in response to this therapy, therefore understanding platinum sensitivity in ovarian cancer would be helpful for improving patients' quality of life and clinical outcomes. HRDetect, utilized to characterize patients' homologous recombination repair deficiency, was used to predict patients' response to platinum-based chemotherapy. However, whether each of the single features contributing to HRD score is associated with platinum sensitivity remains elusive. Methods: We analyzed the whole-exome sequencing data of 196 patients who received platinum-based chemotherapy from the TCGA database. Genetic features were determined individually to see if they could indicate patients' response to platinum-based chemotherapy and prognosis, then integrated into a Pt-score employing LASSO regression model to assess its predictive performance. Results and discussion: Multiple genetic features, including bi-allelic inactivation of BRCA1/2 genes and genes involved in HR pathway, multiple somatic mutations in genes involved in DNA damage repair (DDR), and previously reported HRD-related features, were found to be stringently associated with platinum sensitivity and improved prognosis. Higher contributions of mutational signature SBS39 or ID6 predicted improved overall survival. Besides, arm-level loss of heterozygosity (LOH) of either chr4p or chr5q predicted significantly better disease-free survival. Notably, some of these features were found independent of HRD. And SBS3, an HRD-related feature, was found irrelevant to platinum sensitivity. Integrated all candidate markers using the LASSO model to yield a Pt-score, which showed better predictive ability compared to HRDetect in determining platinum sensitivity and predicting patients' prognosis, and this performance was validated in an independent cohort. The outcomes of our study will be instrumental in devising effective strategies for treating ovarian cancer with platinum-based chemotherapy.

High EVI1 and PARP1 expression as favourable prognostic markers in high-grade serous ovarian carcinoma

BACKGROUND

Mechanisms of development and progression of high-grade serous ovarian cancer (HGSOC) are poorly understood. EVI1 and PARP1, part of TGF-ß pathway, are upregulated in cancers with DNA repair deficiencies with DNA repair deficiencies and may influce disease progression and survival. Therefore we questioned the prognostic significance of protein expression of EVI1 alone and in combination with PARP1 and analyzed them in a cohort of patients with HGSOC.

METHODS

For 562 HGSOC patients, we evaluated EVI1 and PARP1 expression by immunohistochemical staining on tissue microarrays with QuPath digital semi-automatic positive cell detection.

RESULTS

High EVI1 expressing (> 30% positive tumor cells) HGSOC were associated with improved progression-free survival (PFS) (HR = 0.66, 95% CI: 0.504-0.852, p = 0.002) and overall survival (OS) (HR = 0.45, 95% CI: 0.352-0.563, p < 0.001), including multivariate analysis. Most interestingly, mutual high expression of both proteins identifies a group with particularly good prognosis. Our findings were proven technically and clinically using bioinformatical data sets for single-cell sequencing, copy number variation and gene as well as protein expression.

CONCLUSIONS

EVI1 and PARP1 are robust prognostic biomarkers for favorable prognosis in HGSOC and imply further research with respect to their reciprocity.

HRD Testing of Ovarian Cancer in Routine Practice: What Are We Dealing With?

Assessment of homologous recombination deficiency (HRD) status is now essential for ovarian cancer patient management. The aim of our study was to analyze the influence of ethnic variations, tumor purity, and neoadjuvant chemotherapy (CT) on the determination of HRD scores as well as to evaluate feasibility of HRD testing with the Amoy HRD Focus Assay in routine clinical practice. The HRD status, including the BRCA status and genomic scar score (GSS), was analyzed in 452 ovarian cancer specimens. The successful rate of HRD testing was 86% (388/452). The BRCA mutational rate was 29% (114/388); 252 samples (65%) were classified as HRD-positive. Our data demonstrate the feasibility of internal HRD testing by the AmoyDx HRD Focus Panel for high-grade serous ovarian cancer (HGSOC), showing results similar to other methods. The HRD rate in the Russian population is very similar to those of other European populations, as is the BRCA mutation frequency. The most substantial contribution to HRD level diversity is testing criteria depending on intrahospital arrangements. The analysis shows that biallelic BRCA alterations had higher GSS compared with those with monoallelic inactivation, consistent with positive HRD status. The study indicates that grades 1-2 of the pathological response caused by chemotherapy affect HRD scores and suggests controlling for tumor purity of 40% or more as a critical factor for GSS measurement.

Accurate tumor purity determination is critical for the analysis of homologous recombination deficiency (HRD)

Homologous recombination deficiency (HRD) is a predictive marker for response to poly (ADP-ribose) polymerase inhibitors (PARPi) in ovarian carcinoma. HRD scores have entered routine diagnostics, but the influence of algorithms, parameters and confounders has not been analyzed comprehensively. A series of 100 poorly differentiated ovarian carcinoma samples was analyzed using whole exome sequencing (WES) and genotyping. Tumor purity was determined using conventional pathology, digital pathology, and two bioinformatic methods. HRD scores were calculated from copy number profiles determined by Sequenza and by Sclust either with or without fixed tumor purity. Tumor purity determination by digital pathology combined with a tumory purity informed variant of Sequenza served as reference method for HRD scoring. Seven tumors had deleterious mutations in BRCA1/2, 12 tumors had deleterious mutations in other homologous recombination repair (HRR) genes, 18 tumors had variants of unknown significance (VUS) in BRCA1/2 or other HRR genes, while the remaining 63 tumors had no relevant alterations. Using the reference method for HRD scoring, 68 tumors were HRD-positive. HRDsum determined by WES correlated strongly with HRDsum determined by single nucleotide polymorphism (SNP) arrays (R = 0.85). Conventional pathology systematically overestimated tumor purity by 8% compared to digital pathology. All investigated methods agreed on classifying the deleterious BRCA1/2-mutated tumors as HRD-positive, but discrepancies were observed for some of the remaining tumors. Discordant HRD classification of 11% of the tumors was observed comparing the tumor purity uninformed default of Sequenza and the reference method. In conclusion, tumor purity is a critical factor for the determination of HRD scores. Assistance by digital pathology helps to improve accuracy and imprecision of its estimation.

Direct prediction of Homologous Recombination Deficiency from routine histology in ten different tumor types with attention-based Multiple Instance Learning: a development and validation study

Background

Homologous Recombination Deficiency (HRD) is a pan-cancer predictive biomarker that identifies patients who benefit from therapy with PARP inhibitors (PARPi). However, testing for HRD is highly complex. Here, we investigated whether Deep Learning can predict HRD status solely based on routine Hematoxylin & Eosin (H&E) histology images in ten cancer types.

Methods

We developed a fully automated deep learning pipeline with attention-weighted multiple instance learning (attMIL) to predict HRD status from histology images. A combined genomic scar HRD score, which integrated loss of heterozygosity (LOH), telomeric allelic imbalance (TAI) and large-scale state transitions (LST) was calculated from whole genome sequencing data for n=4,565 patients from two independent cohorts. The primary statistical endpoint was the Area Under the Receiver Operating Characteristic curve (AUROC) for the prediction of genomic scar HRD with a clinically used cutoff value.

Results

We found that HRD status is predictable in tumors of the endometrium, pancreas and lung, reaching cross-validated AUROCs of 0.79, 0.58 and 0.66. Predictions generalized well to an external cohort with AUROCs of 0.93, 0.81 and 0.73 respectively. Additionally, an HRD classifier trained on breast cancer yielded an AUROC of 0.78 in internal validation and was able to predict HRD in endometrial, prostate and pancreatic cancer with AUROCs of 0.87, 0.84 and 0.67 indicating a shared HRD-like phenotype is across tumor entities.

Conclusion

In this study, we show that HRD is directly predictable from H&E slides using attMIL within and across ten different tumor types.

Germline variants associated with breast cancer in Khakass women of North Asia

INTRODUCTION

Variants in the BRCA1/2 genes are responsible for familial breast cancer. Numerous studies showed a different spectrum of BRCA variants among breast cancer patients of different Ethnicity origin. In the available literature, no previous research has focused on breast cancer-associated variants among the Khakass people (the indigenous people of the Russian Federation).

METHODS

Twenty-six Khakass breast cancer patients were enrolled in the study. Genomic DNA was isolated from blood samples and used to prepare libraries using a Hereditary Cancer Solution kit. Next-generation sequencing (NGS) was performed using the MiSeq System (Illumina, USA).

RESULTS

In our study, 12% of patients (3/26) carried a single pathogenic variant; 54% of patients (14/26) carried variants of uncertain significance (VUS) or conflicting variants; and 35% of patients (9/26) did not carry any clinically significant variants. Germline pathogenic variant in the ATM gene (rs780619951, NC_000011.10:g.108259022C > T) was identified in two unrelated patients with a family history of cancer (7.6%, 2/26). The pathogenic truncating variant in the ATM gene (p. R805* or c.2413C > T) leads to the nonfunctional version of the protein. This variant has been earlier reported in individuals with a family history of breast cancer.

CONCLUSIONS

Our pilot study describes the germline variant in the ATM gene associated with breast cancer in Khakass women of North Asia.

Homologous recombination deficiency in triple-negative breast cancer: Multi-scale transcriptomics reveals distinct tumor microenvironments and limitations in predicting immunotherapy response

BACKGROUND

Triple-negative breast cancer (TNBC) is the most aggressive subtype of breast cancer and has the highest proportion of homologous recombination deficiency (HRD). HRD has been considered a biomarker of response to immune checkpoint inhibitors (ICIs), but the reality is more complicated. A comprehensive comparison of the tumor microenvironment (TME) in HRD and non-HRD TNBC samples may be helpful.

METHODS

Datasets from single-cell, spatial, and bulk RNA-sequencing were collected to explore the role of HRD in the development of TME at multiple scales. Based on the findings in the TME, machine learning algorithms were used to construct a response prediction model in eleven ICI therapy cohorts.

RESULTS

A more exhausted phenotype of T cells and a more tolerogenic phenotype of dendritic cells were found in the non-HRD group. HRD reprograms the predominant phenotype of cancer-associated fibroblasts (CAFs) from myofibroblastic CAFs to inflammatory-like CAFs. As interactions between myofibroblastic CAFs and other cells, DPP4-chemokines associated with reduced immune cell recruitment were unique in the non-HRD group. The prediction model based on DPP4-related genes had acceptable performance in predicting response, prognosis, and immune cell content. Higher HRD scores in bulk RNA-sequencing samples indicated more activated immune cell function, but not higher immune cell content, which may be affected by factors such as antigen-presenting capacity.

CONCLUSIONS

Based on multi-scale transcriptomics, our findings comprehensively reveal differences in the TME between HRD and non-HRD samples. Combining HRD with the prediction model or other methods for assessing immune cell content, may better predict response to ICIs in TNBC.

Homologous Recombination Repair Gene Mutations to Predict Olaparib Plus Bevacizumab Efficacy in the First-Line Ovarian Cancer PAOLA-1/ENGOT-ov25 Trial

PURPOSE

The PAOLA-1/ENGOT-ov25 trial of maintenance olaparib plus bevacizumab for newly diagnosed advanced high-grade ovarian cancer demonstrated a significant progression-free survival (PFS) benefit over placebo plus bevacizumab, particularly in patients with homologous recombination deficiency (HRD)-positive tumors. We explored whether mutations in non-BRCA1 or BRCA2 homologous recombination repair (non-BRCA HRRm) genes predicted benefit from olaparib plus bevacizumab in PAOLA-1.

METHODS

Eight hundred and six patients were randomly assigned (2:1). Tumors were analyzed using the Myriad MyChoice HRD Plus assay to assess non-BRCA HRRm and HRD status; HRD was based on a genomic instability score (GIS) of ≥ 42. In this exploratory analysis, PFS was assessed in patients harboring deleterious mutations using six non-BRCA HRR gene panels, three devised for this analysis and three previously published.

RESULTS

The non-BRCA HRRm prevalence ranged from 30 of 806 (3.7%) to 79 of 806 (9.8%) depending on the gene panel used, whereas 152 of 806 (18.9%) had non-BRCA1 or BRCA2 mutation HRD-positive tumors. The majority of tumors harboring non-BRCA HRRm had a low median GIS; however, a GIS of > 42 was observed for tumors with mutations in five HRR genes (BLM, BRIP1, RAD51C, PALB2, and RAD51D). Rates of gene-specific biallelic loss were variable (0% to 100%) in non-BRCA HRRm tumors relative to BRCA1-mutated (99%) or BRCA2-mutated (86%) tumors. Across all gene panels tested, hazard ratios for PFS (95% CI) ranged from 0.92 (0.51 to 1.73) to 1.83 (0.76 to 5.43).

CONCLUSION

Acknowledging limitations of small subgroup sizes, non-BRCA HRRm gene panels were not predictive of PFS benefit with maintenance olaparib plus bevacizumab versus placebo plus bevacizumab in PAOLA-1, irrespective of the gene panel tested. Current gene panels exploring HRRm should not be considered a substitute for HRD determined by BRCA mutation status and genomic instability testing in first-line high-grade ovarian cancer.

Increased HRD score in cisplatin resistant penile cancer cells

BACKGROUND/INTRODUCTION

Penile cancer is a rare disease in demand for new therapeutic options. Frequently used combination chemotherapy with 5 fluorouracil (5-FU) and cisplatin (CDDP) in patients with metastatic penile cancer mostly results in the development of acquired drug resistance. Availability of cell culture models with acquired resistance against standard therapy could help to understand molecular mechanisms underlying chemotherapy resistance and to identify candidate treatments for an efficient second line therapy.

METHODS

We generated a cell line from a humanpapilloma virus (HPV) negative penile squamous cell carcinoma (UKF-PEC-1). This cell line was subject to chronic exposure to chemotherapy with CDDP and / or 5-FU to induce acquired resistance in the newly established chemo-resistant sublines (PEC-1^rCDDP²⁵⁰⁰, adapted to 2500 ng/ml CDDP; UKF-PEC-1^r5-FU⁵⁰⁰, adapted to 500 ng/ml 5- FU; UKF-PEC1^rCDDP²⁵⁰⁰/^r5-FU⁵⁰⁰, adapted to 2500 ng/ml CDDP and 500 ng/ml 5 -FU). Afterwards cell line pellets were formalin-fixed, paraffin embedded and subject to sequencing as well as testing for homologous recombination deficiency (HRD). Additionally, exemplary immunohistochemical stainings for p53 and gammaH2AX were applied for verification purposes. Finally, UKF-PEC-1^rCDDP²⁵⁰⁰, UKF-PEC-1^r5-FU⁵⁰⁰, UKF-PEC1^rCDDP²⁵⁰⁰/^r5-FU⁵⁰⁰, and UKF-PEC-3 (an alternative penis cancer cell line) were tested for sensitivity to paclitaxel, docetaxel, olaparib, and rucaparib.

RESULTS AND CONCLUSIONS

The chemo-resistant sublines differed in their mutational landscapes. UKF-PEC-1^rCDDP²⁵⁰⁰ was characterized by an increased HRD score, which is supposed to be associated with increased PARP inhibitor and immune checkpoint inhibitor sensitivity in cancer. However, UKF-PEC-1^rCDDP²⁵⁰⁰ did not display sensitivity to PARP inhibitors.

Prediction of homologous recombination deficiency from cancer gene expression data

Objective

Homologous recombination deficiency (HRD) is the main mechanism of tumorigenesis in some cancers. HRD causes abnormal double-strand break repair, resulting in genomic scars. Some scoring HRD tests have been approved as companion diagnostics of polyadenosine diphosphate-ribose polymerase (PARP) inhibitor treatment. This study aimed to build an HRD prediction model using gene expression data from various cancer types.

Methods

The cancer genome atlas data were used for HRD prediction modeling. A total of 10,567 cases of 33 cancer types were included, and expression data from 5128 out of 20,502 genes were included as predictors. A penalized logistic regression model was chosen as a modeling technique.

Results

The area under the curve of the receiver operating characteristic curve of HRD status prediction was 0.98 for the training set and 0.93 for the test set. The accuracy of HRD status prediction was 0.93 for the training set and 0.88 for the test set.

Conclusions

Our study suggests that the HRD prediction model based on penalized logistic regression using gene expression data can be used to select patients for treatment with PARP inhibitors.

Integrative Genomic Tests in Clinical Oncology

Many clinical decisions in oncology practice rely on the presence or absence of an alteration in a single genetic locus, be it a pathogenic variant in a hereditary cancer gene or activating mutation in a drug target. In addition, there are integrative tests that produce continuous variables and evaluate complex characteristics of the entire tumor genome. Microsatellite instability (MSI) analysis identifies tumors with the accumulation of mutations in short repetitive nucleotide sequences. This procedure is utilized in Lynch syndrome diagnostic pipelines and for the selection of patients for immunotherapy. MSI analysis is well-established for colorectal malignancies, but its applications in other cancer types lack standardization and require additional research. Homologous repair deficiency (HRD) indicates tumor sensitivity to PARP inhibitors and some cytotoxic drugs. HRD-related "genomic scars" are manifested by a characteristic pattern of allelic imbalances, accumulation of deletions with flanking homology, and specific mutation signatures. The detection of the genetic consequences of HRD is particularly sophisticated and expensive, as it involves either whole genome sequencing (WGS) or the utilization of large next-generation sequencing (NGS) panels. Tumor mutation burden (TMB) can be determined by whole exome sequencing (WES) or middle-throughput NGS multigene testing. Although TMB is regarded as an agnostic indicator of tumor sensitivity to immunotherapy, the clinical utility of this test is proven only for a few cancer types.