Abstract
Diffuse glioma is a highly heterogeneous central nervous system tumor that is refractory to conventional therapy. Residual glioma cells escape from surgery and chemoradiotherapy, leading to lethal recurrence. Understanding the molecular mechanism of this recurrence process is critical to the development of successful therapies. Here, we analyzed whole-exome sequencing (WES) data of 97 paired primary and recurrent samples from 46 patients with glioma via a uniform pipeline. Clonality and phylogenetic analyses revealed that branching evolution was widespread in the recurrent process of gliomas. Recurrent tumors continued to evolve independently with chemoradiotherapy and harbored multiple recurrence-selected genetic alterations, such as amplification of PPFIBP1, PDE4DIP, and KRAS, deletion of TNFRSF14, DCC, CDKN2A, and MSH6, and mutations in ATRX, ARID1A, KEL, TP53, MSH6, and KMT2B. Meanwhile, truncal variants within partial driver genes were identified among primary and recurrent gliomas, suggesting that they might be ideal therapeutic targets. Intriguingly, the immunogenicity of recurrent gliomas did not increase significantly compared to the primary tumors. Genomic analysis of recurrent gliomas provided an opportunity to identify potentially clinically informative alterations not detected in clinically sampled primary tumors.