基于高通量测序分析骨肉瘤发生和转移的分子特征

许明贤; 李瑞希; 郭伟堂; 程鸣哲; 邹昌业

doi:10.13471/j.cnki.j.sun.yat-sen.univ(med.sci).2022.0614

您当前的位置：

首页 >

文章列表页 >

基于高通量测序分析骨肉瘤发生和转移的分子特征

临床研究 | 更新时间：2024-02-19

- 基于高通量测序分析骨肉瘤发生和转移的分子特征
- Molecular Characteristics of Osteosarcoma Occurrence and Metastasis Based on High Throughput Sequencing Technology
- 中山大学学报（医学科学版） 2022年43卷第6期页码：985-994
- 作者机构：
  
  1.中山大学附属第一医院骨肿瘤科，广东广州510080
  2.中山大学附属第一医院影像科，广东广州510080
- 作者简介：
  
  许明贤，硕士，研究方向：骨肿瘤，E-mail：hsumingxian@foxmail.com
- 基金信息：
  
  国家自然科学基金(81001194)
- DOI：10.13471/j.cnki.j.sun.yat-sen.univ(med.sci).2022.0614
  中图分类号： R738.1
- 纸质出版日期：2022-11-20，
  
  收稿日期：2022-05-06，
扫描看全文
许明贤,李瑞希,郭伟堂等.基于高通量测序分析骨肉瘤发生和转移的分子特征[J].中山大学学报（医学科学版）,2022,43(06):985-994.

XU Ming-xian,LI Rui-xi,GUO Wei-tang,et al.Molecular Characteristics of Osteosarcoma Occurrence and Metastasis Based on High Throughput Sequencing Technology[J].Journal of Sun Yat-sen University(Medical Sciences),2022,43(06):985-994.
许明贤,李瑞希,郭伟堂等.基于高通量测序分析骨肉瘤发生和转移的分子特征[J].中山大学学报（医学科学版）,2022,43(06):985-994. DOI： 10.13471/j.cnki.j.sun.yat-sen.univ(med.sci).2022.0614.

XU Ming-xian,LI Rui-xi,GUO Wei-tang,et al.Molecular Characteristics of Osteosarcoma Occurrence and Metastasis Based on High Throughput Sequencing Technology[J].Journal of Sun Yat-sen University(Medical Sciences),2022,43(06):985-994. DOI： 10.13471/j.cnki.j.sun.yat-sen.univ(med.sci).2022.0614.

摘要

目的

本研究通过第二代高通量测序（NGS）检测配对骨肉瘤转移灶和原发灶标本，分析骨肉瘤原发灶与转移灶的基因图谱差异，以期发现促进骨肉瘤发生和转移的相关分子及可能机制。

方法

12例转移患者的原发灶与转移病灶肿瘤组织样本，使用第二代高通量测序（NGS）进行检测，其中9对panel检测（678个基因）以及3对全外显子检测（WES），分析比较骨肉瘤原发灶与转移灶的基因图谱差异：基因拷贝数变异通过EXCAVATOR检测；DNA层面的融合突变通过Lumpy检测；RNA层面的融合突变通过Defuse + STAR-Fusion，并通过Circos图显示突变的分布；Metascape 用于分析差异基因的 GO 和 KEGG 信号通路富集；使用Pyclone+Citup / LICHEE进行克隆进化分析。

结果

骨肉瘤基因总体突变模式主要以基因扩增为主，主要包括FLCN（37.5%），GID4（37.5%），TP53（33.3%），ATRX（25%），CALR（25%），CCND3（25%），CCNE1（25%），NOCR1（25%），TFEB（25%），VEGFA（25%）等基因。GO聚类结果提示细胞周期通路突变频率最高。KT1、PLCG2、EGFR这3个基因在转移灶显著富集的信号通路中多次出现，可能与骨肉瘤的转移密切相关。转移灶肿瘤负荷突变（TMB）频率显著高于原发灶（

=0.013）。3例进行WES检测的转移患者均呈现线性克隆进化，提示骨肉瘤转移基因的突变可能呈次序性累积。同时我们确定了可能在骨肉瘤进展中发挥作用的新候选基因，包括PLCG2、MYO15A与PEX6。

结论

骨肉瘤发生与转移的基因突变模式主要以基因扩增为主，骨肉瘤患者肿瘤负荷突变频率在转移灶中显著高于原发灶，转移患者存在相互关联的线性基因克隆进化。本研究发现了可能在骨肉瘤进展中发挥作用的3个新基因PLCG2、MYO15A和PEX6。

Abstract

Objective

The purpose of this study was to explore the difference in expression of gene between the osteosarcoma metastasis and primary specimens through NGS detection， in order to find the molecular mechanism of inducing osteosarcoma occurrence and metastasis.

Methods

Next generation sequencing （NGS） was used to detect the tumor tissue samples of paired primary and metastatic lesions in 12 patients with metastasis， including 9 pairs of Panel detection （678 genes） and 3 pairs of whole exon detection （WES）. Copy number variation was detected by Excavator. DNA mutations were detected by Lumpy fusion， and RNA mutations were detected with Defuse + Star fusion. The distribution of mutations was shown by Circos. Metascape was used to analyze the enrichment of Go and KEGG signal pathways of differential genes. Pyclone+Citup/LICHEE was used for clonal evolution analysis.

Results

The most frequent gene mutation mode of osteosarcoma was gene amplification， high-frequency mutation genes were FLCN（37.5%），GID4（37.5%），TP53（33.3%），ATRX（25%），CALR（25%），CCND3（25%），CCNE1（25%），NOCR1（25%），TFEB（25%），VEGFA（25%）， and the mutation frequency of cell cycle pathway was the highest. KT1， PLCG2 and EGFR were closely related to the metastasis of osteosarcoma. The frequency of tumor mutation burden （TMB） in metastatic lesions was significantly higher than primary lesions （P = 0.013）. Three patients detected by WES showed linear clonal evolution， suggesting that the mutations of osteosarcoma metastasis genes may be accumulated sequentially. We identified new candidate genes that may play an important role in the progression of osteosarcoma， including PLCG2， MYO15A and PEX6.

Conclusions

The most frequent gene mutation pattern of osteosarcoma occurrence and metastasis is gene amplification， and the frequency of tumor mutation burden （TMB） in metastatic lesions is significantly higher than that in primary lesions. There is an interrelated linear gene clonal evolution in patients with metastasis. This study suggests that PLCG2， MYO15A and PEX6 may play an important role in the progression of osteosarcoma.

关键词

骨肉瘤肺转移第二代高通量测序基因图谱基因突变

Keywords

osteosarcomalung metastasisNGSgene mapgene mutation

references

Rainusso N， Wang LL， Yustein JT. The adolescent and young adult with cancer： state of the art-bone tumors［J］. Curr Oncol Rep， 2013， 15（4）： 296-307.

Sayles LC， Breese MR， Koehne AL， et al. Genome-informed targeted therapy for osteosarcoma［J］. Cancer Discov， 2019， 9（1）： 46-63.

Chen X， Bahrami A， Pappo A， et al. Recurrent somatic structural variations contribute to tumorigenesis in pediatric osteosarcoma［J］. Cell Rep， 2014， 7（1）： 104-112.

Perry JA， Kiezun A， Tonzi P， et al. Complementary genomic approaches highlight the PI3K/mTOR pathway as a common vulnerability in osteosarcoma［J］. Proc Natl Acad Sci U S A， 2014， 111（51）： E5564-5573.

Wu CC， Beird HC， Andrew Livingston J， et al. Immuno-genomic landscape of osteosarcoma［J］. Nat Commun， 2020， 11（1）： 1008.

International Cancer Genome C， Hudson TJ， Anderson W， et al. International network of cancer genome projects［J］. Nature， 2010， 464（7291）： 993-998.

Forbes SA， Beare D， Gunasekaran P， et al. COSMIC： exploring the world's knowledge of somatic mutations in human cancer［J］. Nucleic Acids Res， 2015， 43（Database issue）： D805-811.

Hutter C， Zenklusen JC. The cancer genome atlas： creating lasting value beyond its data［J］. Cell， 2018， 173（2）： 283-285.

Pareek CS， Smoczynski R， Tretyn A. Sequencing technologies and genome sequencing［J］. J Appl Genet， 2011， 52（4）： 413-435.

Andersson C， Fagman H， Hansson M， et al. Profiling of potential driver mutations in sarcomas by targeted next generation sequencing［J］. Cancer Genet， 2016， 209（4）： 154-160.

Parsons DW， Roy A， Yang Y， et al. Diagnostic yield of clinical tumor and germline whole-exome sequencing for children with solid tumors［J］. JAMA Oncol， 2016， 2（5）： 616-624.

Zehir A， Benayed R， Shah RH， et al. Mutational landscape of metastatic cancer revealed from prospective clinical sequencing of 10，000 patients［J］. Nat Med， 2017， 23（6）： 703-713.

Chakravarty D， Gao J， Phillips SM， et al. OncoKB： a precision oncology knowledge base［J］. JCO Precis Oncol， 2017：PO.17.00011.

Metzker ML. Sequencing technologies - the next generation［J］. Nat Rev Genet， 2010， 11（1）： 31-46.

Luo Z， Xiao L， Li J， et al. Integrative analysis reveals driver long non-coding RNAs in osteosarcoma［J］. Medicine （Baltimore）， 2019， 98（6）： e14302.

Czarnecka AM， Synoradzki K， Firlej W， et al. Molecular biology of osteosarcoma［J］. Cancers （Basel）， 2020， 12（8）：2130.

Stricker HM， Rommerswinkel N， Keil S， et al. The phospholipase D inhibitor FIPI potently blocks EGF-induced calcium signaling in human breast cancer cells［J］. Cell Commun Signal， 2021， 19（1）： 43.

Semenza GL. Targeting HIF-1 for cancer therapy［J］. Nat Rev Cancer， 2003， 3（10）： 721-732.

Wang J， Ni J， Song D， et al. MAT1 facilitates the lung metastasis of osteosarcoma through upregulation of AKT1 expression［J］. Life Sci， 2019， 234： 116771.

Sun L， Li J， Yan B. Gene expression profiling analysis of osteosarcoma cell lines［J］. Mol Med Rep， 2015， 12（3）： 4266-4272.

Freeman SS， Allen SW， Ganti R， et al. Copy number gains in EGFR and copy number losses in PTEN are common events in osteosarcoma tumors［J］. Cancer， 2008， 113（6）： 1453-1461.

Do SI， Jung WW， Kim HS， et al. The expression of epidermal growth factor receptor and its downstream signaling molecules in osteosarcoma［J］. Int J Oncol， 2009， 34（3）： 797-803.

Bridge JA， Nelson M， Mccomb E， et al. Cytogenetic findings in 73 osteosarcoma specimens and a review of the literature［J］. Cancer Genet Cytogenet， 1997， 95（1）： 74-87.

Tarkkanen M， Karhu R， Kallioniemi A， et al. Gains and losses of DNA sequences in osteosarcomas by comparative genomic hybridization［J］. Cancer Res， 1995， 55（6）： 1334-1338.

Tarkkanen M， Elomaa I， Blomqvist C， et al. DNA sequence copy number increase at 8q： a potential new prognostic marker in high-grade osteosarcoma［J］. Int J Cancer， 1999， 84（2）： 114-121.

Chalmers ZR， Connelly CF， Fabrizio D， et al. Analysis of 100，000 human cancer genomes reveals the landscape of tumor mutational burden［J］. Genome Med， 2017， 9（1）： 34.

Johnson DB， Frampton GM， Rioth MJ， et al. Targeted next generation sequencing identifies markers of response to PD-1 blockade［J］. Cancer Immunol Res， 2016， 4（11）： 959-967.

Ramos KJ， Somayaji R， Nichols DP， et al. Comparative effectiveness research in pediatric respiratory disease： promise and pitfalls［J］. Paediatr Drugs， 2018， 20（1）： 1-7.

Chen C， Xie L， Ren T， et al. Immunotherapy for osteosarcoma： fundamental mechanism， rationale， and recent breakthroughs［J］. Cancer Lett， 2021， 500： 1-10.

Pich O， Muinos F， Lolkema MP， et al. The mutational footprints of cancer therapies［J］. Nat Genet， 2019， 51（12）： 1732-1740.

Rihova K， Ducka M， Zambo IS， et al. Transcription factor c-Myb： novel prognostic factor in osteosarcoma［J］. Clin Exp Metastasis， 2022， 39（2）： 375-390.

浏览量

下载量

CSCD

文章被引用时，请邮件提醒。

提交

工具集

关联资源

RNA结合蛋白PTBP1与磷酸化-AXL共表达在骨肉瘤中的临床意义

骨肉瘤细胞外泌体激活肺组织成纤维细胞炎性信号通路

聚乙二醇修饰脂质体阿霉素联合顺铂单次给药治疗晚期骨肉瘤的剂量递增试验

槲皮素诱导耐甲氨蝶呤骨肉瘤细胞系U2-OS/MTX300凋亡

膜反向斑点杂交法检测广东人G6PD基因6种常见突变