Single-cell transcriptome analysis reveals subtype-specific clonal evolution and microenvironmental changes in liver metastasis of pancreatic adenocarcinoma and their clinical implications.

BACKGROUND: Intratumoral heterogeneity (ITH) and tumor microenvironment (TME) of pancreatic ductal adenocarcinoma (PDAC) play important roles in tumor evolution and patient outcomes. However, the precise characterization of diverse cell populations and their crosstalk associated with PDAC progression and metastasis is still challenging. METHODS: We performed single-cell RNA sequencing (scRNA-seq) of treatment-naïve primary PDAC samples with and without paired liver metastasis samples to understand the interplay between ITH and TME in the PDAC evolution and its clinical associations. RESULTS: scRNA-seq analysis revealed that even a small proportion (22%) of basal-like malignant ductal cells could lead to poor chemotherapy response and patient survival and that epithelial-mesenchymal transition programs were largely subtype-specific. The clonal homogeneity significantly increased with more prevalent and pronounced copy number gains of oncogenes, such as KRAS and ETV1, and losses of tumor suppressor genes, such as SMAD2 and MAP2K4, along PDAC progression and metastasis. Moreover, diverse immune cell populations, including naïve SELLhi regulatory T cells (Tregs) and activated TIGIThi Tregs, contributed to shaping immunosuppressive TMEs of PDAC through cellular interactions with malignant ductal cells in PDAC evolution. Importantly, the proportion of basal-like ductal cells negatively correlated with that of immunoreactive cell populations, such as cytotoxic T cells, but positively correlated with that of immunosuppressive cell populations, such as Tregs. CONCLUSION: We uncover that the proportion of basal-like subtype is a key determinant for chemotherapy response and patient outcome, and that PDAC clonally evolves with subtype-specific dosage changes of cancer-associated genes by forming immunosuppressive microenvironments in its progression and metastasis.

Korea4K: whole genome sequences of 4,157 Koreans with 107 phenotypes derived from extensive health check-ups.

BACKGROUND: Phenome-wide association studies (PheWASs) have been conducted on Asian populations, including Koreans, but many were based on chip or exome genotyping data. Such studies have limitations regarding whole genome-wide association analysis, making it crucial to have genome-to-phenome association information with the largest possible whole genome and matched phenome data to conduct further population-genome studies and develop health care services based on population genomics. RESULTS: Here, we present 4,157 whole genome sequences (Korea4K) coupled with 107 health check-up parameters as the largest genomic resource of the Korean Genome Project. It encompasses most of the variants with allele frequency >0.001 in Koreans, indicating that it sufficiently covered most of the common and rare genetic variants with commonly measured phenotypes for Koreans. Korea4K provides 45,537,252 variants, and half of them were not present in Korea1K (1,094 samples). We also identified 1,356 new genotype-phenotype associations that were not found by the Korea1K dataset. Phenomics analyses further revealed 24 significant genetic correlations, 14 pleiotropic associations, and 127 causal relationships based on Mendelian randomization among 37 traits. In addition, the Korea4K imputation reference panel, the largest Korean variants reference to date, showed a superior imputation performance to Korea1K across all allele frequency categories. CONCLUSIONS: Collectively, Korea4K provides not only the largest Korean genome data but also corresponding health check-up parameters and novel genome-phenome associations. The large-scale pathological whole genome-wide omics data will become a powerful set for genome-phenome level association studies to discover causal markers for the prediction and diagnosis of health conditions in future studies.

Establishment of a patient-specific avatar organoid model derived from endoscopic ultrasonography-guided fine needle biopsy for timely clinical application in pancreatic ductal adenocarcinoma.

BACKGROUND AND AIMS: Pancreatic ductal adenocarcinoma (PDAC) has the worst survival rate among tumors. At the time of diagnosis, more than 80% of PDACs are considered to be surgically unresectable, and there is an unmet need for treatment options in these inoperable PDACs. This study aimed to establish a patient-derived organoid (PDO) platform from endoscopic ultrasound-guided fine needle biopsy (EUS-FNB) collected at diagnosis and to determine its clinical applicability for the timely treatment of unresectable PDAC. METHODS: Patients with suspected PDAC were prospectively enrolled at the Samsung Medical Center from 2015 to 2019. PDAC tissues were acquired by means of EUS-FNB to establish PDAC PDOs, which were comprehensively analyzed for histology, genomic sequencing, and high-throughput screening (HTS) drug sensitivity test. RESULTS: PDAC PDOs were established with a success rate of 83.2% (94/113). It took approximately 3 weeks from acquiring minimal EUS-FNB specimens to generating sufficient PDAC PDOs for the simultaneous HTS drug sensitivity test and genomic sequencing. The high concordance between PDAC tissues and matched PDOs was confirmed, and whole-exome sequencing revealed the increased detection of genetic alterations in PDOs compared with EUS-FNB tissues. The HTS drug sensitivity test showed clinical correlation between the ex vivo PDO response and the actual chemotherapeutic response of the study patients in the real world (13 out of 15 cases). In addition, whole-transcriptome sequencing identified candidate genes associated with nab-paclitaxel resistance, such as ITGB7, ANPEP, and ST3GAL1. CONCLUSIONS: This PDAC PDO platform allows several therapeutic drugs to be tested within a short time window and opens the possibility for timely personalized medicine as a "patient avatar model" in clinical practice.

Impaired binding affinity of YTHDC1 with METTL3/METTL14 results in R-loop accumulation in myelodysplastic neoplasms with DDX41 mutation.

DEAD box helicase 41 (DDX41) mutations are the most prevalent predisposition to familial myelodysplastic syndrome (MDS). However, the precise roles of these variants in the pathogenesis of MDS have yet to be elucidated. Here, we discovered a novel mechanism by which DDX41 contributes to R-loop-induced DNA damage responses (DDR) in cooperation with the m6A-METTL complex (MAC) and YTHDC1 using DDX41 knockout (KO) and DDX41 knock-in (KI, R525H, Y259C) cell lines as well as primary samples from MDS patients. Compared to wild type (WT), DDX41 KO and KI led to increased levels of m6A RNA methylated R-loop. Interestingly, we found that DDX41 regulates m6A/R-loop levels by interacting with MAC components. Further, DDX41 promoted the recruitment of YTHDC1 to R-loops by promoting the binding between METTL3 and YTHDC1, which was dysregulated in DDX41-deficient cells, contributing to genomic instability. Collectively, we demonstrated that DDX41 plays a key role in the physiological control of R-loops in cooperation with MAC and YTHDC1. These findings provide novel insights into how defects in DDX41 influence MDS pathogenesis and suggest potential therapeutic targets for the treatment of MDS.

Enhancing Post-Surgical Rehabilitation Outcomes in Patients with Chronic Ankle Instability: Impact of Subtalar Joint Axis Balance Exercises Following Arthroscopic Modified Broström Operation.

Background and Objectives: This study aimed to evaluate the effects of subtalar joint axis-based balance exercises on the anterior talofibular ligament (ATFL) thickness, ankle strength, and ankle stability after an arthroscopic modified Broström operation (AMBO) for chronic ankle instability (CAI). Materials and Methods: The study included 47 patients diagnosed with CAI who underwent AMBO and were randomly divided into three groups: control (n = 11), general balance exercise (n = 17), and subtalar joint axis balance exercise (n = 19), regardless of the affected area. Participants in the exercise rehabilitation group performed exercises for 60 min twice a week for six weeks, starting six weeks after AMBO. ATFL thickness, ankle strength, and ankle dynamic stability were measured using musculoskeletal ultrasonography, Biodex, and Y-balance test, respectively, before and after treatment. Results: Compared with the remaining groups, the subtalar joint axis balance exercise group had reduced ATFL thickness (p = 0.000), improved ankle strength for eversion (p = 0.000) and inversion (p = 0.000), and enhanced ankle stability (p = 0.000). Conclusions: The study results suggest that subtalar joint axis-based balance exercises may contribute to the early recovery of the ankle joint after AMBO.

A role for mutations in AK9 and other genes affecting ependymal cells in idiopathic normal pressure hydrocephalus.

Idiopathic normal pressure hydrocephalus (iNPH) is an enigmatic neurological disorder that develops after age 60 and is characterized by gait difficulty, dementia, and incontinence. Recently, we reported that heterozygous CWH43 deletions may cause iNPH. Here, we identify mutations affecting nine additional genes (AK9, RXFP2, PRKD1, HAVCR1, OTOG, MYO7A, NOTCH1, SPG11, and MYH13) that are statistically enriched among iNPH patients. The encoded proteins are all highly expressed in choroid plexus and ependymal cells, and most have been associated with cilia. Damaging mutations in AK9, which encodes an adenylate kinase, were detected in 9.6% of iNPH patients. Mice homozygous for an iNPH-associated AK9 mutation displayed normal cilia structure and number, but decreased cilia motility and beat frequency, communicating hydrocephalus, and balance impairment. AK9+/- mice displayed normal brain development and behavior until early adulthood, but subsequently developed communicating hydrocephalus. Together, our findings suggest that heterozygous mutations that impair ventricular epithelial function may contribute to iNPH.

Single-cell RNA Sequencing Reveals Novel Cellular Factors for Response to Immunosuppressive Therapy in Aplastic Anemia.

Aplastic anemia (AA) is a lethal hematological disorder; however, its pathogenesis is not fully understood. Although immunosuppressive therapy (IST) is a major treatment option for AA, one-third of patients do not respond to IST and its resistance mechanism remains elusive. To understand AA pathogenesis and IST resistance, we performed single-cell RNA sequencing (scRNA-seq) of bone marrow (BM) from healthy controls and patients with AA at diagnosis. We found that CD34+ early-stage erythroid precursor cells and PROM1+ hematopoietic stem cells were significantly depleted in AA, which suggests that the depletion of CD34+ early-stage erythroid precursor cells and PROM1+ hematopoietic stem cells might be one of the major mechanisms for AA pathogenesis related with BM-cell hypoplasia. More importantly, we observed the significant enrichment of CD8+ T cells and T cell-activating intercellular interactions in IST responders, indicating the association between the expansion and activation of T cells and the positive response of IST in AA. Taken together, our findings represent a valuable resource offering novel insights into the cellular heterogeneity in the BM of AA and reveal potential biomarkers for IST, building the foundation for future precision therapies in AA.

Predicting preterm birth using machine learning techniques in oral microbiome.

Preterm birth prediction is essential for improving neonatal outcomes. While many machine learning techniques have been applied to predict preterm birth using health records, inflammatory markers, and vaginal microbiome data, the role of prenatal oral microbiome remains unclear. This study aimed to compare oral microbiome compositions between a preterm and a full-term birth group, identify oral microbiome associated with preterm birth, and develop a preterm birth prediction model using machine learning of oral microbiome compositions. Participants included singleton pregnant women admitted to Jeonbuk National University Hospital between 2019 and 2021. Subjects were divided into a preterm and a full-term birth group based on pregnancy outcomes. Oral microbiome samples were collected using mouthwash within 24 h before delivery and 16S ribosomal RNA sequencing was performed to analyze taxonomy. Differentially abundant taxa were identified using DESeq2. A random forest classifier was applied to predict preterm birth based on the oral microbiome. A total of 59 women participated in this study, with 30 in the preterm birth group and 29 in the full-term birth group. There was no significant difference in maternal clinical characteristics between the preterm and the full-birth group. Twenty-five differentially abundant taxa were identified, including 22 full-term birth-enriched taxa and 3 preterm birth-enriched taxa. The random forest classifier achieved high balanced accuracies (0.765 ± 0.071) using the 9 most important taxa. Our study identified 25 differentially abundant taxa that could differentiate preterm and full-term birth groups. A preterm birth prediction model was developed using machine learning of oral microbiome compositions in mouthwash samples. Findings of this study suggest the potential of using oral microbiome for predicting preterm birth. Further multi-center and larger studies are required to validate our results before clinical applications.

Genomic characteristics of triple negative apocrine carcinoma: a comparison to triple negative breast cancer.

Apocrine carcinoma is a rare breast cancer subtype. As such, the genomic characteristics of apocrine carcinoma with triple negative immunohistochemical results (TNAC), which has been treated as triple negative breast cancer (TNBC), have not been revealed. In this study, we evaluated the genomic characteristics of TNAC compared to TNBC with low Ki-67 (LK-TNBC). In the genetic analysis of 73 TNACs and 32 LK-TNBCs, the most frequently mutated driver gene in TNAC was TP53 (16/56, 28.6%), followed by PIK3CA (9/56, 16.1%), ZNF717 (8/56, 14.3%), and PIK3R1 (6/56, 10.71%). Mutational signature analysis showed enrichment of defective DNA mismatch repair (MMR)-related signatures (SBS6 and SBS21) and the SBS5 signature in TNAC, whereas an APOBEC activity-associated mutational signature (SBS13) was more prominent in LK-TNBC (Student's t test, p < 0.05). In intrinsic subtyping, 38.4% of TNACs were classified as luminal A, 27.4% as luminal B, 26.0% as HER2-enriched (HER2-E), 2.7% as basal, and 5.5% as normal-like. The basal subtype was the most dominant subtype (43.8%) in LK-TNBC (p < 0.001), followed by luminal B (21.9%), HER2-E (21.9%), and luminal A (12.5%). In the survival analysis, TNAC had a five-year disease-free survival (DFS) rate of 92.2% compared to 59.1% for LK-TNBC (P = 0.001) and a five-year overall survival (OS) rate of 95.3% compared to 74.6% for LK-TNBC (P = 0.0099). TNAC has different genetic characteristics and better survival outcomes than LK-TNBC. In particular, normal-like and luminal A subtypes in TNAC have much better DFS and OS than other intrinsic subtypes. Our findings are expected to impact medical practice for patients diagnosed with TNAC.

Experimental systems for the analysis of mutational signatures: no 'one-size-fits-all' solution.

Cells constantly accumulate mutations, which are caused by replication errors, as well as through the action of endogenous and exogenous DNA-damaging agents. Mutational patterns reflect the status of DNA repair machinery and the history of genotoxin exposure of a given cellular clone. Computationally derived mutational signatures can shed light on the origins of cancer. However, to understand the etiology of cancer signatures, they need to be compared with experimental signatures, which are obtained from the isogenic cell lines or organisms under controlled conditions. Experimental mutational patterns were instrumental in understanding the nature of signatures caused by mismatch repair and BRCA deficiencies. Here, we describe how different cell lines and model organisms were used in recent years to decipher mutational signatures observed in cancer genomes and provide examples of how data from different experimental systems complement and support each other.

Uniportal versus biportal endoscopic spine surgery: a comprehensive review.

INTRODUCTION: Endoscopic spine surgery techniques are minimally invasive alternatives to conventional open surgery for degenerative spinal diseases. Clinical studies and meta-analyses have proven the usefulness of uniportal full-endoscopic spine procedures. However, a steep learning curve is a critical barrier for endoscopic procedures. Recently, biportal endoscopic spine surgeries have been developed to make it easier for spine surgeons to learn and perform. Consequently, the biportal approach has gained popularity among aspiring endoscopic spine surgeons. This review compared the characteristics of uniportal and biportal surgeries to help spine surgeons perform endoscopic procedures more effectively. AREA COVERED: The review analyzed English-language clinical literature in Core databases and compared uniportal and biportal endoscopic spine surgery techniques. Clinical studies have compared the technical principles of both techniques, and the authors suggested appropriate strategies for learning and practicing endoscopic procedures. EXPERT OPINION: Uniportal, full-endoscopic spine surgery is a minimally invasive procedure that preserves muscles and uses a keyhole approach under local anesthesia. In contrast, biportal surgery is more familiar to a surgeon and can be performed more widely, although the keyhole approach is limited. Aspiring endoscopic spine surgeons can learn either method according to the surgeon's preference and the clinical situation.

A heterozygous mutation in UBE2H in a patient with developmental delay leads to an aberrant brain development in zebrafish.

BACKGROUND: Ubiquitin-related rare diseases are generally characterized by developmental delays and mental retardation, but the exact incidence or prevalence is not yet fully understood. The clinical application of next-generation sequencing for pediatric seizures and developmental delay of unknown causes has become common in studies aimed at identification of a causal gene in patients with ubiquitin-related rare diseases that cannot be diagnosed using conventional fluorescence in situ hybridization or chromosome microarray tests. Our study aimed to investigate the effects of ubiquitin-proteasome system on ultra-rare neurodevelopmental diseases, through functional identification of candidate genes and variants. METHODS: In our present work, we carried out genome analysis of a patient with clinical phenotypes of developmental delay and intractable convulsion, to identify causal mutations. Further characterization of the candidate gene was performed using zebrafish, through gene knockdown approaches. Transcriptomic analysis using whole embryos of zebrafish knockdown morphants and additional functional studies identified downstream pathways of the candidate gene affecting neurogenesis. RESULTS: Through trio-based whole-genome sequencing analysis, we identified a de novo missense variant of the ubiquitin system-related gene UBE2H (c.449C>T; p.Thr150Met) in the proband. Using zebrafish, we found that Ube2h is required for normal brain development. Differential gene expression analysis revealed activation of the ATM-p53 signaling pathway in the absence of Ube2h. Moreover, depletion of ube2h led to induction of apoptosis, specifically in the differentiated neural cells. Finally, we found that a missense mutation in zebrafish, ube2h (c.449C>T; p.Thr150Met), which mimics a variant identified in a patient with neurodevelopmental defects, causes aberrant Ube2h function in zebrafish embryos. CONCLUSION: A de novo heterozygous variant in the UBE2H c.449C>T (p.Thr150Met) has been identified in a pediatric patient with global developmental delay and UBE2H is essential for normal neurogenesis in the brain.

Anion receptors with nitrone C-H hydrogen bond donors.

We report the use of nitrone C-H groups as hydrogen bond donors for binding anions. Acyclic anion receptors with two nitrone moieties were synthesized by condensation reactions between aryl-aldehydes and m-phenylenedihydroxylamines. The solid-state structure of dinitrones revealed C-H⋯O hydrogen bonds. Anion binding properties were investigated using NMR titration with halide and phosphate salts.

GPR143 controls ESCRT-dependent exosome biogenesis and promotes cancer metastasis.

Exosomes transport a variety of macromolecules and modulate intercellular communication in physiology and disease. However, the regulation mechanisms that determine exosome contents during exosome biogenesis remain poorly understood. Here, we find that GPR143, an atypical GPCR, controls the endosomal sorting complex required for the transport (ESCRT)-dependent exosome biogenesis pathway. GPR143 interacts with HRS (an ESCRT-0 Subunit) and promotes its association to cargo proteins, such as EGFR, which subsequently enables selective protein sorting into intraluminal vesicles (ILVs) in multivesicular bodies (MVBs). GPR143 is elevated in multiple cancers, and quantitative proteomic and RNA profiling of exosomes in human cancer cell lines showed that the GPR143-ESCRT pathway promotes secretion of exosomes that carry unique cargo, including integrins signaling proteins. Through gain- and loss-of-function studies in mice, we show that GPR143 promotes metastasis by secreting exosomes and increasing cancer cell motility/invasion through the integrin/FAK/Src pathway. These findings provide a mechanism for regulating the exosomal proteome and demonstrate its ability to promote cancer cell motility.

RAMP: response-aware multi-task learning with contrastive regularization for cancer drug response prediction.

The accurate prediction of cancer drug sensitivity according to the multiomics profiles of individual patients is crucial for precision cancer medicine. However, the development of prediction models has been challenged by the complex crosstalk of input features and the resistance-dominant drug response information contained in public databases. In this study, we propose a novel multidrug response prediction framework, response-aware multitask prediction (RAMP), via a Bayesian neural network and restrict it by soft-supervised contrastive regularization. To utilize network embedding vectors as representation learning features for heterogeneous networks, we harness response-aware negative sampling, which applies cell line-drug response information to the training of network embeddings. RAMP overcomes the prediction accuracy limitation induced by the imbalance of trained response data based on the comprehensive selection and utilization of drug response features. When trained on the Genomics of Drug Sensitivity in Cancer dataset, RAMP achieved an area under the receiver operating characteristic curve > 89%, an area under the precision-recall curve > 59% and an $\textrm{F}_1$ score > 52% and outperformed previously developed methods on both balanced and imbalanced datasets. Furthermore, RAMP predicted many missing drug responses that were not included in the public databases. Our results showed that RAMP will be suitable for the high-throughput prediction of cancer drug sensitivity and will be useful for guiding cancer drug selection processes. The Python implementation for RAMP is available at https://github.com/hvcl/RAMP.

Cellular plasticity and immune microenvironment of malignant pleural effusion are associated with EGFR-TKI resistance in non-small-cell lung carcinoma.

Malignant pleural effusion (MPE) is a complication of lung cancer that can be used as an alternative method for tissue sampling because it is generally simple and minimally invasive. Our study evaluated the diagnostic potential of non-small-cell lung carcinoma (NSCLC)-associated MPE in terms of understanding tumor heterogeneity and identifying response factors for EGFR tyrosine kinase inhibitor (TKI) therapy. We performed a single-cell RNA sequencing analysis of 31,743 cells isolated from the MPEs of 9 patients with NSCLC (5 resistant and 4 sensitive to EGFR TKI) with EGFR mutations. Interestingly, lung epithelial precursor-like cells with upregulated GNB2L1 and CAV1 expression were enriched in the EGFR TKI-resistant group. Moreover, GZMK upregulated transitional effector T cells, and plasmacytoid dendritic cells were significantly enriched in the EGFR TKI-resistant patients. Our results suggest that cellular plasticity and immunosuppressive microenvironment in MPEs are potentially associated with the TKI response of patients with EGFR-mutated NSCLC.

Randomized multicenter phase II trial of prophylactic irradiation of para-aortic lymph nodes in advanced cervical cancer according to tumor hypoxia: Korean Radiation Oncology Group (KROG 07-01) study.

We conducted a prospective phase II study on whether extended-field irradiation (EFI) confers survival benefits depending on hypoxic markers in locally advanced uterine cervical cancer (LAUCC). RNA-seq was performed to identify immune and hypoxic gene signatures. A total of 288 patients were randomized to either EFI or pelvic radiotherapy (PRT). All patients completed chemoradiotherapy. Overall, significantly higher 5-year para-aortic recurrence free survival (PARFS) rate occurred in EFI (97.6%) than in PRT group (87.2%), with marginal tendency to improve disease-free survival (DFS; 78% vs 70%, P = .066). Subgroup analyses were performed based on carbonic anhydrase 9 (CA9)-only positive, CA9/hypoxia-inducible factor (HIF) double positive and CA9 negative. In the CA9-only positive, EFI successfully increased 5-year PARFS (100% vs 76.4%, P = .010), resulting in significantly improved long-term DFS (85.7% vs 54.7%, P = .023) compared to the PRT, while there was no such benefit of EFI in the CA9/HIFs double positive. RNA-seq analysis identified distinct immunehigh subgroup with negative correlation with hypoxia gene signatures (R = -.37, P < .01), which showed a higher 5-year DFS than the immunelow (P = .032). Hypoxia-related genes were upregulated in the CA9/HIFs double positive compared to CA9 negative (P < .05). Only 17.4% of patients in CA9-negative group showed immunelow signatures, while 40.0% of patients in the double-positive group exhibited immunelow signatures. In conclusion, EFI improved PARFS significantly in all patients, but therapeutic efficacy of EFI in terms of improved DFS was solely observed in CA9-only positive LAUCC, and not in CA9/HIFs double-positive subgroup. RNA-seq analysis suggested that hypoxia-induced immunosuppression may be related to treatment resistance in LAUCC.

PWWP2B promotes DNA end resection and homologous recombination.

Genome instability is one of the leading causes of gastric cancers. However, the mutational landscape of driver genes in gastric cancer is poorly understood. Here, we investigate somatic mutations in 25 Korean gastric adenocarcinoma patients using whole-exome sequencing and show that PWWP2B is one of the most frequently mutated genes. PWWP2B mutation correlates with lower cancer patient survival. We find that PWWP2B has a role in DNA double-strand break repair. As a nuclear protein, PWWP2B moves to sites of DNA damage through its interaction with UHRF1. Depletion of PWWP2B enhances cellular sensitivity to ionizing radiation (IR) and impairs IR-induced foci formation of RAD51. PWWP2B interacts with MRE11 and participates in homologous recombination via promoting DNA end-resection. Taken together, our data show that PWWP2B facilitates the recruitment of DNA repair machinery to sites of DNA damage and promotes HR-mediated DNA double-strand break repair. Impaired PWWP2B function might thus cause genome instability and promote gastric cancer development.

A functional genomic approach to actionable gene fusions for precision oncology.

Fusion genes represent a class of attractive therapeutic targets. Thousands of fusion genes have been identified in patients with cancer, but the functional consequences and therapeutic implications of most of these remain largely unknown. Here, we develop a functional genomic approach that consists of efficient fusion reconstruction and sensitive cell viability and drug response assays. Applying this approach, we characterize ~100 fusion genes detected in patient samples of The Cancer Genome Atlas, revealing a notable fraction of low-frequency fusions with activating effects on tumor growth. Focusing on those in the RTK-RAS pathway, we identify a number of activating fusions that can markedly affect sensitivity to relevant drugs. Last, we propose an integrated, level-of-evidence classification system to prioritize gene fusions systematically. Our study reiterates the urgent clinical need to incorporate similar functional genomic approaches to characterize gene fusions, thereby maximizing the utility of gene fusions for precision oncology.