Temporal and spatial characteristics of tumor evolution in a mouse model of oral squamous cell carcinoma.

OBJECTIVES: We aimed to elucidate the temporal and spatial characteristics of tumor evolution in an oral squamous cell carcinoma (OSCC) mouse model with higher burden of lymphatic metastasis through high-throughput sequencing. METHODS: The OSCC model was established in 9 mice. DNA was extracted from the tumors of primary tongue lesions and disseminated tumor cells (DTCs) of submandibular gland lymph nodes and bone marrow, and then whole genome sequencing was performed. After bioinformatics analysis, somatic single-nucleotide variants (SSNVs) and copy number variations (CNVs) data were obtained. Based on SSNVs, clonal architecture and ancestor-descendant relationships among tumor cell subclones were elucidated. RESULTS: A total of 238 tumor-related SSNVs with 120 high-frequency mutated genes were obtained from 36 samples of 9 mice by whole-genome sequencing. The number of unique SSNVs in the primary lesion, submandibular lymph node and bone marrow was greater than the number of shared SSNVs. Furthermore, the primary lesion-originated subclones, which were identified by SSNVs, were also detected in submandibular lymph nodes in the early stage of oral carcinogenesis. Moreover, at different histopathological stages, unique subclones were also identified in DTCs isolated from lymph nodes. CONCLUSION: Tumor heterogeneity is significant in primary tumor cells and disseminated tumor cells. OSCC cells probably disseminate to lymph nodes in the early stage of oral carcinogenesis. OSCC is characterized by polyclonal dissemination, and the evolutionary trajectory of DTCs is potentially dominated by the tumor microenvironment.