GLI1 confers resistance to PARP inhibitors by activating the DNA damage repair pathway.

Identifying the mechanisms of action of anticancer drugs is an important step in the development of new drugs. In this study, we established a comprehensive screening platform consisting of 68 oncogenes (MANO panel), encompassing 243 genetic variants, to identify predictive markers for drug efficacy. Validation was performed using drugs that targeted EGFR, BRAF, and MAP2K1, which confirmed the utility of this functional screening panel. Screening of a BRCA2-knockout DLD1 cell line (DLD1-KO) revealed that cells expressing SMO and GLI1 were resistant to olaparib. Gene set enrichment analysis identified genes associated with DNA damage repair that were enriched in cells overexpressing SMO and GLI1. The expression of genes associated with homologous recombination repair (HR), such as the FANC family and BRCA1/2, was significantly upregulated by GLI1 expression, which is indicative of PARP inhibitor resistance. Although not all representative genes of the nucleotide excision repair (NER) pathway were upregulated, NER activity was enhanced by GLI1. The GLI1 inhibitor was effective against DLD1-KO cells overexpressing GLI1 both in vitro and in vivo. Furthermore, the combination therapy of olaparib and GLI1 inhibitor exhibited a synergistic effect on DLD1-KO, suggesting the possible clinical application of GLI1 inhibitor targeting cancer with defective DNA damage repair. This platform enables the identification of biomarkers associated with drug sensitivity, and is a useful tool for drug development.

Genomic landscape of comprehensive genomic profiling in patients with malignant solid tumors in Japan.

BACKGROUND: Comprehensive genomic profiling (CGP) can aid the discovery of clinically useful, candidate antitumor agents; however, the variant annotations sometimes differ among the various types of CGP tests as well as the public database. The aim of this study is to clarify the genomic landscape of evaluating detected variants in patients with a malignant solid tumor. METHODS: The present, cross-sectional study used data from 57,084 patients with a malignant solid tumor who underwent CGP at the Center for Cancer Genomics and Advanced Therapeutics (C-CAT) between June 1, 2019 and August 18, 2023. The pathogenicity of the variants was annotated using public databases. RESULTS: As a result of re-annotation of the detected variants, 20.1% were pathogenic and 1.4% were benign. The mean number of pathogenic variants was 4.30 (95% confidence interval: 4.27-4.32) per patient. Of the entire cohort, 5.7% had no pathogenic variant. The co-occurrence of the genes depended on the tumor type. Germline findings were detected in 6.2%, 8.8%, and 15.8% of the patients using a tumor/normal panel, tumor-only panel, and liquid panel, respectively, with the most common gene being BRCA2 followed by TP53 and BRCA1. CONCLUSIONS: The detected variants should be re-annotated because several benign variants or variants of unknown significance were included in the CGP, and the genomic landscape derived from these results will help researchers and physicians interpret the results of CGP tests. The method of extracting presumptive, germline, pathogenic variants from patients using a tumor-only panel or circulating tumor DNA panel requires improvement.

A molecular glue RBM39-degrader induces synthetic lethality in cancer cells with homologous recombination repair deficiency.

E7820 and Indisulam (E7070) are sulfonamide molecular glues that modulate RNA splicing by degrading the splicing factor RBM39 via ternary complex formation with the E3 ligase adaptor DCAF15. To identify biomarkers of the antitumor efficacy of E7820, we treated patient-derived xenograft (PDX) mouse models established from 42 patients with solid tumors. The overall response rate was 38.1% (16 PDXs), and tumor regression was observed across various tumor types. Exome sequencing of the PDX genome revealed that loss-of-function mutations in genes of the homologous recombination repair (HRR) system, such as ATM, were significantly enriched in tumors that responded to E7820 (p = 4.5 × 103). Interestingly, E7820-mediated double-strand breaks in DNA were increased in tumors with BRCA2 dysfunction, and knockdown of BRCA1/2 transcripts or knockout of ATM, ATR, or BAP1 sensitized cancer cells to E7820. Transcriptomic analyses revealed that E7820 treatment resulted in the intron retention of mRNAs and decreased transcription, especially for HRR genes. This induced HRR malfunction probably leads to the synthetic lethality of tumor cells with homologous recombination deficiency (HRD). Furthermore, E7820, in combination with olaparib, exerted a synergistic effect, and E7820 was even effective in an olaparib-resistant cell line. In conclusion, HRD is a promising predictive biomarker of E7820 efficacy and has a high potential to improve the prognosis of patients with HRD-positive cancers.

Embryonal and alveolar rhabdomyosarcoma in adolescents/young adults, adults and older adults: a population-based cohort study.

BACKGROUND: The clinical characteristics, outcomes, and prognostic factors of adult embryonal rhabdomyosarcomas (ERMS) and alveolar rhabdomyosarcomas (ARMS), particularly the differences among adolescents/young adults (AYA), adults, and older adults, remain unclear. We assessed the clinicopathological features and survival outcomes of adult patients with ERMS and ARMS in Japan and to compare these features among AYA, adult, and older adult patients. METHODS: We retrospectively analyzed data from the Bone and Soft Tissue Tumor Registry of Japan and enrolled patients aged ≥15 years with ERMS and ARMS. Disease-specific overall survival (DOS) was estimated using the Kaplan-Meier method, and a Cox regression model was used to identify prognostic factors. RESULTS: Among 184 patients with ERMS and ARMS (median age, 27 years; interquartile range, 18-49 years), a high rate of distant and regional nodal metastases was initially observed in 65 (35%) and 66 (36%) cases, respectively. Older age and distant metastasis at first presentation were statistically poor prognostic factors, and histological subtype and site of tumor origin were not associated with DOS. In patients with localized ERMS and ARMS, older age and nodal metastasis were poor prognostic factors; the 5-year DOS rates of patients with and without nodal metastasis were 23% and 72%, respectively. CONCLUSIONS: Older patients with rhabdomyosarcoma had a dismal prognosis, and distant metastasis was a poor prognostic factor. The prognostic factors differed between adult and pediatric patients with rhabdomyosarcoma; biological analyses, such as genome analysis of adult rhabdomyosarcoma and clinical trials with pediatric oncologists, are needed to improve the prognosis of adult rhabdomyosarcoma.

Prognostic factors and treatment outcomes in patients with pleomorphic rhabdomyosarcoma: a population-based cohort study.

BACKGROUND: Pleomorphic rhabdomyosarcoma is a rare sarcoma in adults. The clinical characteristics, outcomes and prognostic factors associated with pleomorphic rhabdomyosarcoma remain unclear. METHODS: We retrospectively analyzed data from the Bone and Soft Tissue Tumor Registry of Japan, and enrolled patients with pleomorphic rhabdomyosarcoma. Disease-specific overall survival, local recurrence-free survival and distant metastasis-free survival were estimated using the Kaplan-Meier method; Cox regression model was used to identify prognostic factors. RESULTS: In total, 182 patients with pleomorphic rhabdomyosarcoma were included. Median age was 63 (range 20-95) years. The lower extremity (48%) was the most frequent tumor origin site, while head and neck were rare (4%). A total of 43 patients (24%) had distant or regional nodal metastases at first presentation. In all cases, the 2-year and 5-year survival rates were 66.3% and 54.1%, respectively. Distant metastasis was a significant poor prognostic factor (Hazard ratio 6.65; 95% confidence intervals, 3.00-14.75, P < 0.0001), with median survival of such patients being 9.4 (95% confidence intervals: 5.3-12.2) months. In 134 localized cases, the 2-year and 5-year survival rates were 91.5% and 68.3%, respectively. Large tumor size and older age were associated with poorer prognosis. Through data from localized and locally curative cases extracted and adjusted by propensity score matching, we found that perioperative chemotherapy did not improve disease-specific overall survival, distant metastasis-free survival or local recurrence-free survival. CONCLUSIONS: Clinical characteristics and outcomes of pleomorphic rhabdomyosarcoma are similar to those of other high-grade soft tissue sarcomas. Pleomorphic rhabdomyosarcoma may be less chemosensitive, and a strategy other than the standard cytotoxic chemotherapy is required to improve its prognosis.

Factors predictive of second-line chemotherapy in soft tissue sarcoma: An analysis of the National Genomic Profiling Database.

Of the drugs used in second-line chemotherapy for soft tissue sarcoma (STS), trabectedin is effective for liposarcoma and leiomyosarcoma (L-sarcoma), eribulin for liposarcoma, and pazopanib for non-liposarcoma. The indications for these drugs in STS other than L-sarcoma have not been established. Here we explored the prognosis, mutation profiles, and drug-response factors in STS using real-world big data. Clinicogenomic data on 1761 patients with sarcoma who underwent FoundationOne CDx were obtained from a national database in Japan. Patients with TP53 and KDM2D mutations had a significantly shorter survival period of 253 (95% CI, 99-404) and 330 (95% CI, 20-552) days, respectively, than those without mutations. Non-supervised clustering based on mutation profiles generated 13 tumor clusters. The response rate (RR) to trabectedin was highest in an MDM2-amplification cluster (odds ratio [OR]: 2.2; p = 0.2). The RR was lowest for eribulin in an MDM2-amplification cluster (OR: 0.4; p = 0.03) and highest in a TERT-mutation cluster (OR: 3.0; p = 0.03). The RR was highest for pazopanib in a PIK3CA/PTEN-wild type cluster (OR: 2.1; p = 0.03). In particular, patients harboring mutations in genes regulating the PI3K/Akt/mTOR pathway had a lower RR than patients without mutations (OR: 0.3; p = 0.04). In STS, mutation profiles were more useful in predicting the drug response than histology. The present study demonstrated the potential of tailored therapy guided by mutation profiles established by comprehensive genomic profiling testing in optimizing second-line chemotherapy for STS. The findings of this study will hopefully contribute some valuable insights into enhancing STS treatment strategies and outcomes.

Calcified deep somatic soft tissue leiomyoma of the leg: A case report.

Deep somatic soft tissue leiomyoma is a rare benign smooth muscle tumor, with around 60 cases reported thus far in the literature. Herein, we report a huge leiomyoma in deep somatic soft tissue with massive calcification in the peroneus longus muscle. A 29-year-old male patient presented with a hard, spindle-shaped, 15-cm-diameter mass of the left lower leg. Plain radiography showed marked calcification within the lesion, and magnetic resonance imaging revealed a well-circumscribed lesion with strong contrast enhancement within the peroneus longus. A core needle biopsy was performed to rule out a soft tissue sarcoma, and deep somatic soft tissue leiomyoma was preoperatively diagnosed. A marginal excision was performed, and the patient recovered without event and has been recurrence-free for 16 months. As far as we know, the present soft tissue leiomyoma was one of the largest to be documented, the previously reported cases being up to 14 cm in maximum diameter. Leiomyomas of deep somatic soft tissue should be included in the differential diagnosis of a huge, deep-seated tumor with marked calcification as well as sarcomas, especially in younger patients.

KDM2B-Rearranged Soft Tissue Sarcomas Expand the Concept of BCOR-Associated Sarcoma.

Sarcomas with BCOR genetic alterations (BCOR-associated sarcomas) represent a recently recognized family of soft tissue and bone tumors characterized by BCOR fusion, BCOR internal tandem duplication, or YWHAE::NUTM2B fusion. Histologically, the tumors demonstrate oval to spindle cell proliferation in a variably vascular stroma and overexpression of BCOR and SATB2. Herein, we describe 3 soft tissue sarcomas with KDM2B fusions that phenotypically and epigenetically match BCOR-associated sarcomas. The cases included 1 infant, 1 adolescent, and 1 older patient. All tumors showed histologic findings indistinguishable from those of BCOR-associated sarcomas and were originally diagnosed as such based on the phenotype. However, none of the tumors had BCOR or YWHAE genetic alterations. Instead, targeted RNA sequencing identified in-frame KDM2B::NUTM2B, KDM2B::CREBBP, and KDM2B::DUX4 fusions. KDM2B fusions were validated using reverse-transcription PCR, Sanger sequencing, and in situ hybridization assays. Genome-wide DNA methylation analysis matched all 3 tumors with BCOR-associated sarcomas using the Deutsches Krebsforschungszentrum (DKFZ) classifier and t-distributed stochastic neighbor embedding analysis. One localized tumor showed a flat genome-wide copy number profile, and the patient remained disease-free after treatment. The other tumors showed multiple copy number alterations, including MDM2/CDK4 amplification and/or CDKN2A/B loss, and both tumors metastasized, leading to the patient's death in one of the cases. When tested using KDM2B immunohistochemistry, all 3 KDM2B-rearranged sarcomas showed diffuse strong staining, and all 13 sarcomas with BCOR genetic alterations also demonstrated diffuse, strong, or weak staining. By contrast, among 72 mimicking tumors, only a subset of synovial sarcomas showed focal or diffuse weak KDM2B expression. In conclusion, our study suggests that KDM2B-rearranged soft tissue sarcomas belong to the BCOR-associated sarcoma family and expand its molecular spectrum. This may be related to the known molecular relationship between KDM2B and BCOR in the polycomb repressive complex 1.1. Immunohistochemical analysis of KDM2B is a potentially valuable diagnostic tool for BCOR-associated sarcomas, including those with KDM2B rearrangement.

Extensive analysis of 59 sarcoma-related fusion genes identified pazopanib as a potential inhibitor to COL1A1-PDGFB fusion gene.

Sarcomas are malignant mesenchymal tumors that are extremely rare and divergent. Fusion genes are involved in approximately 30% of sarcomas as driver oncogenes; however, their detailed functions are not fully understood. In this study, we determined the functional significance of 59 sarcoma-related fusion genes. The transforming potential and drug sensitivities of these fusion genes were evaluated using a focus formation assay (FFA) and the mixed-all-nominated-in-one (MANO) method, respectively. The transcriptome was also examined using RNA sequencing of 3T3 cells transduced with each fusion gene. Approximately half (28/59, 47%) of the fusion genes exhibited transformation in the FFA assay, which was classified into five types based on the resulting phenotype. The sensitivity to 12 drugs including multityrosine kinase inhibitors was assessed using the MANO method and pazopanib was found to be more effective against cells expressing the COL1A1-PDGFB fusion gene compared with the others. The downstream MAPK/AKT pathway was suppressed at the protein level following pazopanib treatment. The fusion genes were classified into four subgroups by cluster analysis of the gene expression data and gene set enrichment analysis. In summary, the oncogenicity and drug sensitivity of 59 fusion genes were simultaneously evaluated using a high-throughput strategy. Pazopanib was selected as a candidate drug for sarcomas harboring the COL1A1-PDGFB fusion gene. This assessment could be useful as a screening platform and provides a database to evaluate customized therapy for fusion gene-associated sarcomas.

PTPN11 variant may be a prognostic indicator of IDH-wildtype glioblastoma in a comprehensive genomic profiling cohort.

PURPOSE: Glioblastoma (GBM) is the most common type of primary malignant brain tumor and has a poor prognosis. Identifying novel targets and stratification strategies is urgently needed to improve patient survival. The present study aimed to identify clinically relevant genomic alterations in IDH-wildtype GBM using data from comprehensive genomic profiling (CGP) assays performed nationwide in Japan. METHODS: The CGP assay results of 392 IDH-wildtype GBM cases performed between October 2019 and February 2023 obtained from the Center for Cancer Genomics and Advanced Therapeutics were retrospectively analyzed. RESULTS: The median patient age was 52.5 years, and 207 patients (53%) were male. In the 286 patients for whom survival information was available, a protein-tyrosine phosphatase non-receptor type 11 (PTPN11) variant detected in 20 patients (6.8%) was extracted as the gene associated with significantly shorter overall survival (p = 0.002). Multivariate analysis demonstrated that the PTPN11 variant and poor performance status were independent prognostic indicators. In contrast, no prognostic impact was observed in the cohort in The Cancer Genome Atlas data. The discrepancy in the prognostic impact of the PTPN11 variant from these two pools might have resulted from differences in the biases affecting the survival of patients who underwent a CGP assay, including left-truncation and right-censored bias. However, survival simulation done to adjust for these biases showed that the prognostic impact of the PTPN11 variant was also significant. CONCLUSIONS: The PTPN11 variant was a negative prognostic indicator of IDH-wildtype GBM in the patient cohort with the CGP assay.

BRAF mutations and concurrent alterations in patients with soft tissue sarcoma.

BRAF alterations, including V600E and non-V600E mutations and fusions, in soft tissue sarcoma (STS) have been identified in a limited case series. Here, we aimed to evaluate the frequency of BRAF mutations and concurrent alterations in STS to understand their therapeutic action. In this retrospective analysis, we included data from 1964 patients with advanced STS who underwent comprehensive genomic profiling tests at hospitals in Japan between June 2019 and March 2023. The prevalence of BRAF and recurrent concurrent gene alterations were also investigated. BRAF mutations were detected in 24 (1.2%) of 1964 STS patients, with a median age of 47 (range 1-69) years. BRAF V600E was detected in 11 (0.6%) of the 1964 patients with STS, BRAF non-V600E mutations in 9 (4.6%), and BRAF fusions were detected in 4 (0.2%). BRAF V600E was identified in 4 (0.2%) cases of malignant peripheral nerve sheath tumors. The most common concurrent alteration was CDKN2A (11 cases, 45.8%), and the frequency was equivalent to that of the BRAF V600E (5/11 cases, 45.5%) and non-V600E (5/9 cases, 55.6%) groups. Recurrent concurrent alterations, such as TERT promoter mutations (7 cases, 29.2%), were detected at the same frequency in the V600E and non-V600E groups. In contrast, TP53 alterations (4/9 cases, 44.4%) and mitogen-activated protein kinase (MAPK)-activating genes, including NF1, GNAQ, and GNA11 (3/9 cases, 33.3%), were identified as relatively higher in the non-V600E group than in the V600E group (each 1/11 case, 9.1%). We identified BRAF alterations at a rate of 1.2% in all patients with advanced STS. Among them, BRAF V600E and BRAF fusions account for 45.8% and 16.7%, respectively. Collectively, our findings support the clinical characteristics and therapeutic strategies for patients with BRAF-altered advanced STS.

Prognostic benefit of comprehensive genomic profiling in clinical practice remains uncertain.

The overall survival of patients who received genomically matched therapy was not significantly longer than that of patients receiving treatment only other than genomically matched therapy in the breast invasive ductal carcinoma (A), colorectal adenocarcinoma (B), and pancreatic adenocarcinoma (C) cohorts.

The pathogenic role of the BRCA2 c.7847C>T (p.Ser2616Phe) variant in breast and ovarian cancer predisposition.

Substantial numbers of variants of unknown significance (VUSs) have been identified in BRCA1/2 through genetic testing, which poses a significant clinical challenge because the contribution of these VUSs to cancer predisposition has not yet been determined. Here, we report 10 Japanese patients from seven families with breast or ovarian cancer harboring the BRCA2 c.7847C>T (p.Ser2616Phe) variant that was interpreted as a VUS. This variant recurs only in families from Japan and has not been reported in the global general population databases. A Japanese patient with Fanconi anemia with compound heterozygous variants c.7847C>T (p.Ser2616Phe) and c.475+1G>A in BRCA2 was reported. In silico predictions and quantitative cosegregation analysis suggest a high probability of pathogenicity. The clinical features of the variant carriers were not specific to, but were consistent with, those of patients with hereditary breast and ovarian cancer. A validated functional assay, called the mixed-all-nominated-in-one-BRCA (MANO-B) method and the accurate BRCA companion diagnostic (ABCD) test, demonstrated the deleterious effects of the variant. Altogether, following the American College of Medical Genetics and Genomics and the Association for Molecular Pathology (ACMG/AMP) guidelines, this variant satisfied the "PS3," "PM2," "PM3," and "PP3" criteria. We thus conclude that the BRCA2 c.7847C>T (p.Ser2616Phe) variant is a "likely pathogenic" variant that is specifically observed in the Japanese population, leading to a breast and ovarian cancer predisposition.

High-Throughput Functional Evaluation of MAP2K1 Variants in Cancer.

Activating mutations in mitogen-activated protein kinase kinase 1 (MAP2K1) are involved in a variety of cancers and may be classified according to their RAF dependence. Sensitivity to combined BRAF and MEK treatments is associated with co-mutations of MAP2K1 and BRAF; however, the significance of less frequent MAP2K1 mutations is largely unknown. The transforming potential and drug sensitivity of 100 MAP2K1 variants were evaluated using individual assays and the mixed-all-nominated-in-one method. In addition, A375, a melanoma cell line harboring the BRAF V600E mutation, was used to evaluate the function of the MAP2K1 variants in combination with active RAF signaling. Among a total of 67 variants of unknown significance, 16 were evaluated as oncogenic or likely oncogenic. The drug sensitivity of the individual variants did not vary with respect to BRAF inhibitors, MEK inhibitors (MEKi), or their combination. Sensitivity to BRAF inhibitors was associated with the RAF dependency of the MAP2K1 variants, whereas resistance was higher in RAF-regulated or independent variants compared with RAF-dependent variants. Thus, the synergistic effect of BRAF and MEKis may be observed in RAF-regulated and RAF-dependent variants. MAP2K1 variants exhibit differential sensitivity to BRAF and MEKis, suggesting the importance of individual functional analysis for the selection of optimal treatments for each patient. This comprehensive evaluation reveals precise functional information and provides optimal combination treatment for individual MAP2K1 variants.

Role of perioperative chemotherapy and radiotherapy for localized high-grade malignant peripheral nerve sheath tumor at the extremities and trunk wall: a population-based cohort study.

BACKGROUND: Primary tumor resection is the mainstay of treatment for malignant peripheral nerve sheath tumors. However, the efficacy of perioperative chemotherapy and radiotherapy for malignant peripheral nerve sheath tumors has not been elucidated. METHODS: This retrospective analysis based on a Japanese registry included patients with localized malignant peripheral nerve sheath tumors arising at the extremities and trunk wall. Disease-specific overall survival and local recurrence-free survival were estimated using the Kaplan-Meier method. A Cox regression model was used to identify prognostic factors. Characteristics of groups with or without chemotherapy were adjusted using propensity score matching. RESULTS: In total, 291 patients were included. The 5-year disease-specific overall survival rate was 70.6%. Multivariate analysis of disease-specific overall survival revealed that deep-seated tumors were a poor prognostic factor, but perioperative chemotherapy was not associated with disease-specific overall survival (hazard ratio, 0.81; 95% confidence interval, 0.45-1.43, P = 0.46). Local recurrence was observed in 55 patients (19.0%), and surgical margins (R1 and R2) were significant risk factors. Overall, perioperative chemotherapy did not prolong disease-specific overall survival (5-year disease-specific overall survival: 74.1% vs. 69.3%, P = 0.75) and had limited efficacy in the group with tumor size ≥ 5 cm, although the difference was not statistically significant (5-year disease-specific overall survival: 77.2% vs. 68.6%, respectively, P = 0.13). After adjustment by propensity score matching, perioperative chemotherapy significantly prolonged disease-specific overall survival (5-year disease-specific overall survival: 74.9% vs. 57.1%, P = 0.03), but this effect was not observed in local recurrence-free survival. In all patients, perioperative radiotherapy did not correlate with local recurrence-free survival (hazard ratio, 1.43; 95% confidence interval 0.78-2.62, P = 0.25). CONCLUSIONS: Perioperative chemotherapy had limited efficacy for disease-specific overall survival in patients with localized malignant peripheral nerve sheath tumors.

Dramatic response to entrectinib in a patient with malignant peripheral nerve sheath tumor harboring novel SNRNP70-NTRK3 fusion gene.

Neurotropic tropomyosin receptor kinase (NTRK) gene rearrangements have been reported in limited cases of sarcomas; however, to date, there has been only one report of such rearrangements in malignant peripheral nerve sheath tumors (MPNSTs). Herein, we describe a 51-year-old male patient with a buttock tumor arising from the sciatic nerve, which was diagnosed as MPNST with positive S-100 staining, negative SOX10 staining, and loss of trimethylation at lysine 27 of histone H3 (H3K27me3) confirmed by immunohistochemistry. Soon after the resection of the primary tumor, the patient was found to have pulmonary and lymph node metastases. Chemotherapy with eribulin and trabectedin showed limited effects. However, the patient responded rapidly to pazopanib, but severe side effects caused discontinuation of the treatment. RNA panel testing revealed a novel fusion gene between Small Nuclear Ribonucleoprotein U1 Subunit 70 (SNRNP70) gene and NTRK3 gene. Furthermore, loss of NF1, SUZ12, and CDKN2A genes was confirmed by DNA panel testing, which is compatible with a histological diagnosis of MPNST. SNRNP70 possesses a coiled-coiled domain and seems to induce constitutive activation of NTRK3 through dimerization. In fact, immunohistochemistry revealed diffuse staining of pan-TRK within tumor cells. Treatment with entrectinib, which is an NTRK inhibitor, showed a quick and durable response for 10 months. Although NTRK rearrangements are very rare in MPNST, this case highlights the importance of genetic testing in MPNST, especially using an RNA panel for the detection of rare fusion genes.

Selection bias due to delayed comprehensive genomic profiling in Japan.

Patients with advanced cancer undergo comprehensive genomic profiling in Japan only after treatment options have been exhausted. Patients with a very poor prognosis were not able to undergo profiling tests, resulting in a selection bias called length bias, which makes accurate survival analysis impossible. The actual impact of length bias on the overall survival of patients who have undergone profiling tests is unclear, yet appropriate methods for adjusting for length bias have not been developed. To assess the length bias in overall survival, we established a simulation-based model for length bias adjustment. This study utilized clinicogenomic data of 8813 patients with advanced cancer who underwent profiling tests at hospitals throughout Japan between June 2019 and April 2022. Length bias was estimated by the conditional Kendall τ statistics and was significantly positive for 13 of the 15 cancer subtypes, suggesting a worse prognosis for patients who underwent profiling tests in early timing. The median overall survival time in colorectal, breast, and pancreatic cancer from the initial survival-prolonging chemotherapy with adjustment for length bias was 937 (886-991), 1225 (1152-1368), and 585 (553-617) days, respectively (median; 95% credible interval). Adjusting for length bias made it possible to analyze the prognostic relevance of oncogenic mutations and treatments. In total, 12 tumor-specific oncogenic mutations correlating with poor survival were detected after adjustment. There was no difference in survival between FOLFIRINOX (leucovorin, fluorouracil, irinotecan, and oxaliplatin) or gemcitabine with nab-paclitaxel-treated groups as first-line chemotherapy for pancreatic cancer. Adjusting for length bias is an essential part of utilizing real-world clinicogenomic data.