Durvalumab with or without tremelimumab plus chemotherapy in HRR non-mutated, platinum-resistant ovarian cancer (KGOG 3045): A phase II umbrella trial.

AIM: We investigated the efficacy and safety of durvalumab (D) with or without tremelimumab (T) in addition to single-agent chemotherapy (CT) in patients with platinum-resistant recurrent ovarian cancer (PROC) lacking homologous recombination repair (HRR) gene mutations. PATIENTS AND METHODS: KGOG 3045 was an open-label, investigator-initiated phase II umbrella trial. Patients with PROC without HRR gene mutations who had received ≥2 prior lines of therapy were enrolled. Patients with high PD-L1 expression (TPS ≥25%) were assigned to arm A (D + CT), whereas those with low PD-L1 expression were assigned to arm B (D + T75 + CT). After completing arm B recruitment, patients were sequentially assigned to arms C (D + T300 + CT) and D (D + CT). RESULTS: Overall, 58 patients were enrolled (5, 18, 17, and 18 patients in arms A, B, C, and D, respectively). The objective response rates were 20.0, 33.3, 29.4, and 22.2%, respectively. Grade 3-4 treatment-related adverse events were observed in 20.0, 66.7, 47.1, and 66.7 of patients, respectively, but were effectively managed. Multivariable analysis demonstrated that adding T to D + CT improved progression-free survival (adjusted HR, 0.435; 95% CI, 0.229-0.824; P = 0.011). Favorable response to chemoimmunotherapy was associated with MUC16 mutation (P = 0.0214), high EPCAM expression (P = 0.020), high matrix remodeling gene signature score (P = 0.017), and low FOXP3 expression (P = 0.047). Patients showing favorable responses to D + T + CT exhibited significantly higher EPCAM expression levels (P = 0.008) and matrix remodeling gene signature scores (P = 0.031) than those receiving D + CT. CONCLUSIONS: Dual immunotherapy with chemotherapy showed acceptable response rates and tolerable safety in HRR non-mutated PROC, warranting continued clinical investigation.

Randomized, two-arm, noncomparative phase 2 study of olaparib plus cediranib or durvalumab in HRR-mutated, platinum-resistant ovarian cancer: A substudy of KGOG 3045.

Choosing an optimal concomitant drug for combination with poly-ADP ribose polymerase (PARP) inhibitor based on patient-specific biomarker status may help increase to improve treatment efficacy in patients with ovarian cancer. However, the efficacy and safety of different PARP inhibitor-based combinations in patients with homologous recombination repair (HRR) mutations have not been evaluated in ovarian cancer. In this sub-study of Korean Gynecologic Oncology Group (KGOG) 3045, we compared the efficacy and safety of two olaparib-based combinations and biomarkers of patients with platinum-resistant ovarian cancer with HRR gene mutations. Patients were randomized to receive either olaparib (200 mg twice a day) + cediranib (30 mg daily) (Arm 1, n = 16) or olaparib (300 mg) + durvalumab (1,500 mg once every 4 weeks) (Arm 2, n = 14). The objective response rates for Arm 1 and Arm 2 were 50.0% and 42.9%, respectively. Most patients (83.3%) had BRCA mutations, which were similarly distributed between arms. Grade 3 or 4 treatment-related adverse events were observed in 37.5% and 35.7% of the patients, respectively, but all were managed properly. A high vascular endothelial growth factor signature was associated with favorable outcomes in Arm 1, whereas immune markers (PD-L1 expression [CPS ≥10], CD8, neutrophil-to-lymphocyte ratio and platelet-to-lymphocyte ratio) were associated with favorable outcomes in Arm 2. The activation of homologous recombination pathway upon disease progression was associated with poor response to subsequent therapy. Based on comprehensive biomarker profiling, including immunohistochemistry, whole-exome and RNA sequencing and whole blood-based analyses, we identified biomarkers that could help inform which of the two combination strategies is appropriate given a patient's biomarker status. Our findings have the potential to improve treatment outcome for patients with ovarian cancer in the PARP inhibitor era.

Combined Model-Based Prediction for Non-Invasive Prenatal Screening.

The risk of chromosomal abnormalities in the child increases with increasing maternal age. Although non-invasive prenatal testing (NIPT) is a safe and effective prenatal screening method, the accuracy of the test results needs to be improved owing to various testing conditions. We attempted to achieve a more accurate and robust prediction of chromosomal abnormalities by combining multiple methods. Here, three different methods, namely standard Z-score, normalized chromosome value, and within-sample reference bin, were used for 1698 reference and 109 test samples of whole-genome sequencing. The logistic regression model combining the three methods achieved a higher accuracy than any single method. In conclusion, the proposed method offers a promising approach for increasing the reliability of NIPT.

Renal Cell Carcinoma-Infiltrating CD3low Vγ9Vδ1 T Cells Represent Potentially Novel Anti-Tumor Immune Players.

Due to the highly immunogenic nature of renal cell carcinoma (RCC), the tumor microenvironment (TME) is enriched with various innate and adaptive immune subsets. In particular, gamma-delta (γδ) T cells can act as potent attractive mediators of adoptive cell transfer immunotherapy because of their unique properties such as non-reliance on major histocompatibility complex expression, their ability to infiltrate human tumors and recognize tumor antigens, relative insensitivity to immune checkpoint molecules, and broad tumor cytotoxicity. Therefore, it is now critical to better characterize human γδ T-cell subsets and their mechanisms in RCCs, especially the stage of differentiation. In this study, we aimed to identify γδ T cells that might have adaptive responses against RCC progression. We characterized γδ T cells in peripheral blood and tumor-infiltrating lymphocytes (TILs) in freshly resected tumor specimens from 20 RCC patients. Furthermore, we performed a gene set enrichment analysis on RNA-sequencing data from The Cancer Genome Atlas (TCGA) derived from normal kidneys and RCC tumors to ascertain the association between γδ T-cell infiltration and anti-cancer immune activity. Notably, RCC-infiltrating CD3low Vγ9Vδ1 T cells with a terminally differentiated effector memory phenotype with up-regulated activation/exhaustion molecules were newly detected as predominant TILs, and the cytotoxic activity of these cells against RCC was confirmed in vitro. In an additional analysis of the TCGA RCC dataset, γδ T-cell enrichment scores correlated strongly with those for CTLs, Th1 cells, "exhausted" T cells, and M1 macrophages, suggesting active involvement of γδ T cells in anti-tumor rather than pro-tumor activity, and Vδ1 cells were more abundant than Vδ2 or Vδ3 cells in RCC tumor samples. Thus, we posit that Vγ9Vδ1 T cells may represent an excellent candidate for adoptive immunotherapy in RCC patients with a high risk of relapse after surgery.

SIDR: simultaneous isolation and parallel sequencing of genomic DNA and total RNA from single cells.

Simultaneous sequencing of the genome and transcriptome at the single-cell level is a powerful tool for characterizing genomic and transcriptomic variation and revealing correlative relationships. However, it remains technically challenging to analyze both the genome and transcriptome in the same cell. Here, we report a novel method for simultaneous isolation of genomic DNA and total RNA (SIDR) from single cells, achieving high recovery rates with minimal cross-contamination, as is crucial for accurate description and integration of the single-cell genome and transcriptome. For reliable and efficient separation of genomic DNA and total RNA from single cells, the method uses hypotonic lysis to preserve nuclear lamina integrity and subsequently captures the cell lysate using antibody-conjugated magnetic microbeads. Evaluating the performance of this method using real-time PCR demonstrated that it efficiently recovered genomic DNA and total RNA. Thorough data quality assessments showed that DNA and RNA simultaneously fractionated by the SIDR method were suitable for genome and transcriptome sequencing analysis at the single-cell level. The integration of single-cell genome and transcriptome sequencing by SIDR (SIDR-seq) showed that genetic alterations, such as copy-number and single-nucleotide variations, were more accurately captured by single-cell SIDR-seq compared with conventional single-cell RNA-seq, although copy-number variations positively correlated with the corresponding gene expression levels. These results suggest that SIDR-seq is potentially a powerful tool to reveal genetic heterogeneity and phenotypic information inferred from gene expression patterns at the single-cell level.

The role of PDGFRA as a therapeutic target in young colorectal cancer patients.

BACKGROUND: Young patients with colorectal cancer (CRC) exhibit poor prognoses compared to older patients due to the difficulty in early diagnosis and treatment. However, the underlying molecular characteristics are still unclear. METHODS: We conducted a comprehensive analysis of 49 CRC patients without hereditary CRC using the whole-exome and RNA sequencing with tumor and matched normal samples. A total of 594 TCGA samples and 4 patient-derived cells were utilized for validation. RESULTS: Consensus molecular subtype 4 (CMS4) (53.85%) and CMS2 (38.46%) were enriched in the young (≤ 40 years) and old (> 60 years) age groups, respectively. A CMS4-associated gene, platelet-derived growth factor receptor α (PDGFRA), was significantly upregulated in young patients with CRC (FC = 3.21, p = 0.0001) and was negatively correlated with age (p = 0.0001, R = - 0.526). Moreover, PDGFRA showed a positive co-expression with metastasis-related genes in young CRC patients. In vitro validation confirmed that young patient-derived cells (PDCs) showed an enriched expression of PDGFRA compared to old PDCs and a reduced proliferation rate by knockdown of PDGFRA. Furthermore, young CRC patients were more sensitive to regorafenib, a PDGFRA-targeting drug, than old CRC patients. CONCLUSIONS: Our study suggests that CRC in young patients is associated with CMS4 and PDGFRA. In addition, PDGFRA may serve potential of novel therapeutic strategies and represent a predictive biomarker of response to regorafenib for young CRC patients.

Correlations between metabolic texture features, genetic heterogeneity, and mutation burden in patients with lung cancer.

PURPOSE: This study investigated the correlations between parameters of 18F-fluorodeoxyglucose (FDG) uptake on positron emission tomography (PET) scan and indices of genetic properties, heterogeneity index (HI), and tumor mutation burden (TMB), in patients with lung cancer. METHODS: We produced 106 PET indices for each tumor site that underwent genomic analysis in a total of 176 study subjects (age, 62.0 ± 10.0 y; males, 68.2%), comprising 101 adenocarcinoma (ADC), 29 squamous cell carcinoma (SQCC), and 46 small cell lung cancer (SCLC) patients. We then examined the correlations of the PET parameters with genetic properties of HI and TMB, according to pathology and tumor site. RESULTS: Comparisons between PET parameters and the genetic properties with false discovery rate (FDR) correction revealed that the surface standard uptake value (SUV) entropy of SUV statistics had a significant correlation with HI only in patients with SCLC who underwent a genetic test in lymph nodes (r = 0.592, p = 0.028), whereas PET parameters did not show a significant correlation with HI or TMB in patients with SCLC who underwent a genetic test in lung tissue. In patients with ADC and SQCC, there was no significant correlation between PET parameters and the genetic properties. Although SUVmax showed raw p values less than 0.05 in correlation with HI (r = 0.315, raw p = 0.048) and TMB (r = 0.206, raw p = 0.043) in ADC, and SUVpeak had a raw p value less than 0.05 in correlation with HI (r = 0.394, raw p = 0.046) in SQCC, these parameters were not significant when corrected by FDR. CONCLUSIONS: In this study, surface SUV entropy had a significant correlation with HI in SCLC. Regarding other PET parameters and tumors, no significant correlation with genetic parameters existed.

Evaluation of somatic copy number estimation tools for whole-exome sequencing data.

Whole-exome sequencing (WES) has become a standard method for detecting genetic variants in human diseases. Although the primary use of WES data has been the identification of single nucleotide variations and indels, these data also offer a possibility of detecting copy number variations (CNVs) at high resolution. However, WES data have uneven read coverage along the genome owing to the target capture step, and the development of a robust WES-based CNV tool is challenging. Here, we evaluate six WES somatic CNV detection tools: ADTEx, CONTRA, Control-FREEC, EXCAVATOR, ExomeCNV and Varscan2. Using WES data from 50 kidney chromophobe, 50 bladder urothelial carcinoma, and 50 stomach adenocarcinoma patients from The Cancer Genome Atlas, we compared the CNV calls from the six tools with a reference CNV set that was identified by both single nucleotide polymorphism array 6.0 and whole-genome sequencing data. We found that these algorithms gave highly variable results: visual inspection reveals significant differences between the WES-based segmentation profiles and the reference profile, as well as among the WES-based profiles. Using a 50% overlap criterion, 13-77% of WES CNV calls were covered by CNVs from the reference set, up to 21% of the copy gains were called as losses or vice versa, and dramatic differences in CNV sizes and CNV numbers were observed. Overall, ADTEx and EXCAVATOR had the best performance with relatively high precision and sensitivity. We suggest that the current algorithms for somatic CNV detection from WES data are limited in their performance and that more robust algorithms are needed.

Neurofibromatosis 2 (NF2) controls the invasiveness of glioblastoma through YAP-dependent expression of CYR61/CCN1 and miR-296-3p.

Glioblastoma multiforme (GBM) is the most common and aggressive type of primary brain tumor derived from non-neuronal glial cells. Neurofibromatosis 2 (NF2) protein, also termed as merlin, is a well-known tumor suppressor; however, the molecular mechanism underlying this effect has not yet been fully defined. To investigate the role of NF2 in the invasiveness of GBM, we used two GBM cell lines: NF2-expressing T98G cells and NF2-deficient A172 cells. Knockdown of NF2 increased the invasiveness of T98G cells, whereas NF2-overexpressing A172 cells showed decreased invasive activity. Moreover, re-expression of NF2 reversed the high invasiveness of NF2-silenced T98G cells, indicating that NF2 negatively regulates GBM invasiveness. We further found that the NF2-mediated regulation of invasiveness was dependent on YAP and TEAD2 expression levels. NF2 also controlled the expression of YAP targets, including cysteine-rich angiogenic inducer 61 (CYR61/CCN1), by regulating the nuclear localization of YAP. Silencing of CYR61/CCN1 blocked the increased invasiveness of T98G cells, suggesting that CYR61/CCN1 is required for NF2-mediated invasiveness. Through microRNA microarray analysis, we found that NF2 negatively regulates the expression of miR-296-3p. Overexpression of miR-296-3p suppressed the expression of STAT5A, induced the phosphorylation of STAT3 by downregulating SOCS2, and increased the invasiveness of T98G cells. Taken together, we demonstrate that NF2 negatively controls the invasiveness of GBM through YAP-dependent induction of CYR61/CCN1 and miR-296-3p.

Genome-wide screen identifies novel machineries required for both ciliogenesis and cell cycle arrest upon serum starvation.

Biogenesis of the primary cilium, a cellular organelle mediating various signaling pathways, is generally coordinated with cell cycle exit/re-entry. Although the dynamic cell cycle-associated profile of the primary cilium has been largely accepted, the mechanism governing the link between ciliogenesis and cell cycle progression has been poorly understood. Using a human genome-wide RNAi screen, we identify genes encoding subunits of the spliceosome and proteasome as novel regulators of ciliogenesis. We demonstrate that 1) the mRNA processing-related hits are essential for RNA expression of molecules acting in cilia disassembly, such as AURKA and PLK1, and 2) the ubiquitin-proteasome systems (UPS)-involved hits are necessary for proteolysis of molecules acting in cilia assembly, such as IFT88 and CPAP. In particular, we show that these screen hit-associated mechanisms are crucial for both cilia assembly and cell cycle arrest in response to serum withdrawal. Finally, our data suggest that the mRNA processing mechanism may modulate the UPS-dependent decay of cilia assembly regulators to control ciliary resorption-coupled cell cycle re-entry.

Deciphering intratumor heterogeneity using cancer genome analysis.

Intratumor heterogeneity within individual cancer tissues underlies the numerous phenotypes of cancer. Tumor subclones ultimately affect therapeutic outcomes due to their distinct molecular features. Drug-resistant subclones are present at a low frequency in tissues at the time of biopsy, but can also arise as a result of acquired somatic mutations. A number of different approaches have been utilized to understand the nature of intratumor heterogeneity. Clonal analysis using whole exome or genome sequencing data can help monitor subclones in the context of tumor progression. Multiregional biopsies permit the molecular characterization of subclones within tumors. Deep sequencing has also provided researchers with the ability to measure the low allele fraction variant within a small number of cells. Ultimately, single-cell sequencing will enable the identification of every minor population within a tumor microenvironment. In the clinical context, the ability to identify and monitor the subclonal architecture of a tumor is valuable for the development of precise cancer therapeutic methods.

Deregulation of Retinaldehyde Dehydrogenase 2 Leads to Defective Angiogenic Function of Endothelial Colony-Forming Cells in Pediatric Moyamoya Disease.

OBJECTIVE--: Moyamoya disease (MMD) is a common cause of childhood stroke, in which the abnormal function of the endothelial colony-forming cell (ECFC) plays a key role in the pathogenesis of the disease. This study was designed to identify genes involved in MMD pathogenesis using gene expression profiling and to understand the defective function of MMD ECFCs. APPROACH AND RESULTS--: We compared gene expression profiles of ECFCs isolated from patients with MMD and normal controls. Among the differentially expressed genes, we selected a gene with the most downregulated expression, retinaldehyde dehydrogenase 2 (RALDH2). The activity of RALDH2 in MMD ECFCs was assessed by in vitro tube formation assay and in vivo Matrigel plug assay in the presence of all-trans retinoic acid. The transcriptional control of RALDH2 was tested using ChIP assays on acetyl-histone H3. In the results, MMD ECFCs inefficiently formed capillary tubes in vitro and capillaries in vivo, a defect restored by all-trans retinoic acid treatment. Knockdown of RALDH2 mRNA in normal ECFCs also induced decreased activity of capillary formation in vitro. The decreased level of RALDH2 mRNA in MMD ECFCs was attributed to defective acetyl-histone H3 binding to the promoter region. CONCLUSIONS--: From these results, we conclude that the expression of RALDH2 was epigenetically suppressed in ECFCs from patients with MMD, which may play a key role in their functional impairment.

The diagnostic application of targeted re-sequencing in Korean patients with retinitis pigmentosa.

BACKGROUND: Identification of the causative genes of retinitis pigmentosa (RP) is important for the clinical care of patients with RP. However, a comprehensive genetic study has not been performed in Korean RP patients. Moreover, the genetic heterogeneity found in sensorineural genetic disorders makes identification of pathogenic mutations challenging. Therefore, high throughput genetic testing using massively parallel sequencing is needed. RESULTS: Sixty-two Korean patients with nonsyndromic RP (46 patients from 18 families and 16 simplex cases) who consented to molecular genetic testing were recruited in this study and targeted exome sequencing was applied on 53 RP-related genes. Causal variants were characterised by selecting exonic and splicing variants, selecting variants with low allele frequency (below 1 %), and discarding the remaining variants with quality below 20. The variants were additionally confirmed by an inheritance pattern and cosegregation test of the families, and the rest of the variants were prioritised using in-silico prediction tools. Finally, causal variants were detected from 10 of 18 familial cases (55.5 %) and 7 of 16 simplex cases (43.7 %) in total. Novel variants were detected in 13 of 20 (65 %) candidate variants. Compound heterozygous variants were found in four of 7 simplex cases. CONCLUSION: Panel-based targeted re-sequencing can be used as an effective molecular diagnostic tool for RP.

RAD51 as an immunohistochemistry-based marker of poly(ADP-ribose) polymerase inhibitor resistance in ovarian cancer.

Objective: Effective functional biomarkers that can be readily used in clinical practice to predict poly(ADP-ribose) polymerase inhibitor (PARPi) sensitivity are lacking. With the widespread adoption of PARPi maintenance therapy in ovarian cancer, particularly in patients with BRCA mutation or HR deficiencies, accurately identifying de novo or acquired resistance to PARPi has become critical in clinical practice. We investigated RAD51 immunohistochemistry (IHC) as a functional biomarker for predicting PARPi sensitivity in ovarian cancer. Methods: Ovarian cancer patients who had received PARPi and had archival tissue samples prior to PARPi exposure ("pre-PARPi") and/or after progression on PARPi ("post-PARPi") were selected. RAD51 IHC expression was semi-quantitatively evaluated using the H-score in geminin (a G2/S phase marker)- and γH2AX (a DNA damage marker)-positive tissues. A RAD51 H-score of 20 was used as the cutoff value. Results: In total, 72 samples from 56 patients were analyzed. The median RAD51 H-score was 20 (range: 0-90) overall, 10 (0-190) in pre-PARPi samples (n = 34), and 25 (1-170) in post-PARPi samples (n = 19). Among patients with BRCA mutations, RAD51-low patients had better progression-free survival (PFS) after PARPi treatment than RAD51-high patients (P = 0.029). No difference was found in PFS with respect to the genomic scar score (P = 0.930). Analysis of matched pre- and post-PARPi samples collected from 15 patients indicated an increase in the RAD51 H-score upon progression on PARPi, particularly among pre-PARPi low-RAD51-expressing patients. Conclusion: RAD51 is a potential functional IHC biomarker of de novo and acquired PARPi resistance in BRCA-mutated ovarian cancer and can be used to fine-tune ovarian cancer treatment.

Combined transcriptome, genetic diversity and metabolite profiling in tomato fruit reveals that the ethylene response factor SlERF6 plays an important role in ripening and carotenoid accumulation.

Solanum lycopersicum (tomato) and its wild relatives harbor genetic diversity that yields heritable variation in fruit chemistry that could be exploited to identify genes regulating their synthesis and accumulation. Carotenoids, for example, are essential in plant and animal nutrition, and are the visual indicators of ripening for many fruits, including tomato. Whereas carotenoid synthesis is well characterized, factors regulating flux through the pathway are poorly understood at the molecular level. To exploit the impact of tomato genetic diversity on carotenoids, Solanum pennellii introgression lines were used as a source of defined natural variation and as a resource for the identification of candidate regulatory genes. Ripe fruits were analyzed for numerous fruit metabolites and transcriptome profiles generated using a 12,000 unigene oligoarray. Correlation analysis between carotenoid content and gene expression profiles revealed 953 carotenoid-correlated genes. To narrow the pool, subnetwork analysis of carotenoid-correlated transcription revealed 38 candidates. One candidate for impact on trans-lycopene and ß-carotene accumulation was functionally charaterized, SlERF6, revealing that it indeed influences carotenoid biosynthesis and additional ripening phenotypes. Reduced expression of SlERF6 by RNAi enhanced both carotenoid and ethylene levels during fruit ripening, demonstrating an important role for SlERF6 in ripening, integrating the ethylene and carotenoid synthesis pathways.

Tomato Functional Genomics Database: a comprehensive resource and analysis package for tomato functional genomics.

Tomato Functional Genomics Database (TFGD) provides a comprehensive resource to store, query, mine, analyze, visualize and integrate large-scale tomato functional genomics data sets. The database is functionally expanded from the previously described Tomato Expression Database by including metabolite profiles as well as large-scale tomato small RNA (sRNA) data sets. Computational pipelines have been developed to process microarray, metabolite and sRNA data sets archived in the database, respectively, and TFGD provides downloads of all the analyzed results. TFGD is also designed to enable users to easily retrieve biologically important information through a set of efficient query interfaces and analysis tools, including improved array probe annotations as well as tools to identify co-expressed genes, significantly affected biological processes and biochemical pathways from gene expression data sets and miRNA targets, and to integrate transcript and metabolite profiles, and sRNA and mRNA sequences. The suite of tools and interfaces in TFGD allow intelligent data mining of recently released and continually expanding large-scale tomato functional genomics data sets. TFGD is available at http://ted.bti.cornell.edu.

T-Cell Receptor Repertoire Characteristics Associated with Prognostic Significance in High-Grade Serous Ovarian Carcinoma.

High-grade serous ovarian carcinoma (HGSOC) is a fatal gynecological malignancy. Somatic recombination occurring during T-cell receptor (TCR) development results in TCR diversity, and the TCR repertoire, thus produced, is associated with immune response. This study analyzed the difference in the TCR repertoire and their prognostic significance in 51 patients with HGSOC. The patient's clinical characteristics, gene expression pattern, TCR clonotypes, and degree of tumor-infiltrating leukocytes (TILs) were analyzed, and the patients were divided into groups depending on their recurrence pattern, tumor-infiltrating leukocyte (TIL) score, and homologous recombinant repair pathway deficiency (HRD)-associated mutations. The TCR repertoire was low in patients with recurrence and showed the expansion of eight TCR segments. Interestingly, a few genes correlated with the TCRs also showed a difference in expression according to the prognosis. Among them, seven genes were related to immune responses and KIAA1199 was up-regulated in ovarian cancer. Our study shows that the differences in the TCR repertoire in patients with ovarian cancer and their associated immune pathways could affect the prognosis of HGSOC.

Predictive biomarkers for the responsiveness of recurrent glioblastomas to activated killer cell immunotherapy.

BACKGROUND: Recurrent glioblastoma multiforme (GBM) is a highly aggressive primary malignant brain tumor that is resistant to existing treatments. Recently, we reported that activated autologous natural killer (NK) cell therapeutics induced a marked increase in survival of some patients with recurrent GBM. METHODS: To identify biomarkers that predict responsiveness to NK cell therapeutics, we examined immune profiles in tumor tissues using NanoString nCounter analysis and compared the profiles between 5 responders and 7 non-responders. Through a three-step data analysis, we identified three candidate biomarkers (TNFRSF18, TNFSF4, and IL12RB2) and performed validation with qRT-PCR. We also performed immunohistochemistry and a NK cell migration assay to assess the function of these genes. RESULTS: Responders had higher expression of many immune-signaling genes compared with non-responders, which suggests an immune-active tumor microenvironment in responders. The random forest model that identified TNFRSF18, TNFSF4, and IL12RB2 showed a 100% accuracy (95% CI 73.5-100%) for predicting the response to NK cell therapeutics. The expression levels of these three genes by qRT-PCR were highly correlated with the NanoString levels, with high Pearson's correlation coefficients (0.419 (TNFRSF18), 0.700 (TNFSF4), and 0.502 (IL12RB2)); their prediction performance also showed 100% accuracy (95% CI 73.54-100%) by logistic regression modeling. We also demonstrated that these genes were related to cytotoxic T cell infiltration and NK cell migration in the tumor microenvironment. CONCLUSION: We identified TNFRSF18, TNFSF4, and IL12RB2 as biomarkers that predict response to NK cell therapeutics in recurrent GBM, which might provide a new treatment strategy for this highly aggressive tumor.

Ultrafast prediction of somatic structural variations by filtering out reads matched to pan-genome k-mer sets.

Variant callers typically produce massive numbers of false positives for structural variations, such as cancer-relevant copy-number alterations and fusion genes resulting from genome rearrangements. Here we describe an ultrafast and accurate detector of somatic structural variations that reduces read-mapping costs by filtering out reads matched to pan-genome k-mer sets. The detector, which we named ETCHING (for efficient detection of chromosomal rearrangements and fusion genes), reduces the number of false positives by leveraging machine-learning classifiers trained with six breakend-related features (clipped-read count, split-reads count, supporting paired-end read count, average mapping quality, depth difference and total length of clipped bases). When benchmarked against six callers on reference cell-free DNA, validated biomarkers of structural variants, matched tumour and normal whole genomes, and tumour-only targeted sequencing datasets, ETCHING was 11-fold faster than the second-fastest structural-variant caller at comparable performance and memory use. The speed and accuracy of ETCHING may aid large-scale genome projects and facilitate practical implementations in precision medicine.