上颌尖牙阻生患者牙龈组织的转录组特征分析

口腔生物医学 ›› 2021, Vol. 12 ›› Issue (1): 34-37.

上颌尖牙阻生患者牙龈组织的转录组特征分析

陈彬¹，刘畅²，管晓燕¹，白国辉¹，王苗¹，钱声艳¹，刘建国¹

1. 遵义医科大学口腔疾病研究特色重点实验室
2. 贵州中医药大学药学院

收稿日期:2020-10-20 修回日期:2021-01-20 出版日期:2021-03-25 发布日期:2021-04-01
通讯作者: 刘建国 E-mail:a13985695604@163.com
基金资助:
贵州省科技厅科技合作项目;遵义市科技项目;遵义市红花岗区科计项目;遵义医科大学硕启

Transcriptome analysis of the gingival tissue from impacted maxillary permanent canine

1. The Special Key Laboratory of Oral Diseases Research, Institution of Higher Education in Guizhou Province
2. School of pharmacy, Guizhou University of traditional chinese medicine

Received:2020-10-20 Revised:2021-01-20 Online:2021-03-25 Published:2021-04-01
Contact: guo jianliu E-mail:a13985695604@163.com

摘要/Abstract

摘要： 目的：采集和分析上颌尖牙阻生患者牙龈组织的转录组信息，为其地域差异和病因基因理论奠定基础。方法：使用Resequencing Illumina HiSeq 2000测序平台高通量测序，采集遵义地区5例成人（25～35岁）上颌尖牙阻生患者牙龈组织的转录组信息，并采用FPKM方法计算基因在牙龈组织中的表达量。通过转录组比对，将基因映射蛋白质数据库进行功能注释。基于注释信息，挖掘上颌尖牙阻生患者的转录组特征。结果：高通量测序共获得223?530?890个高质量测序数据，FPKM方法鉴定出19?352个基因表达（FPKM＞0），其中99.63%的基因FPKM值小于1?000。转录组比对显示每个样本约15?412个基因发生了可变剪切，剪切次数约为16?167次。这些可变剪切类型中外显子选择性跳跃剪切所占比例最低（6.72%），外显子跳跃剪切比例最高（48.82%）。此外，利用这批转录组数据共检测出386?414个单核苷酸多态性位点（SNP），其中转换是最主要的突变类型（76.31%）。GO数据库注释显示，转录组SNP涉及到广泛的生物功能，且尤其富集于生物学过程和细胞组分。KEGG通路分析显示，SNP所在基因除了富集于代谢途径、细胞内吞、细胞凋亡、细胞衰老等相关通路，更多的基因则富集于癌症、信号通路以及免疫系统相关通路。结论：基于阻生尖牙牙龈组织的转录组研究显示牙龈和牙发育密切相关。

关键词:

单核苷酸多态性, 牙龈组织, 转录组, 牙周组织

Abstract: Objective:To collect and analyze the whole genome transcription information of gingival tissue of patients with impacted maxillary canine, and lay a foundation for the gene theory of the etiology. Methods:Using resequecing Illumina hiseq2000 sequencing platform, the transcriptome information of gingival tissue of 5 adult patients (25-35 years old) with impacted maxillary canines in Zunyi area was collected by high-throughput sequencing. The expression level and difference of gene in gingival tissue were analyzed and calculated by FPKM method. Results:223530890 high quality sequencing data were obtained, and 19352 genes were indentified to be expressed in the gingival tissue, among which 99.63% genes’ FPKM values were lower than 1000. 15412 genes had alternative splicing (AS), and the aplicing times were 16167. Mutually exclusive exons were found to be the rarest AS type (6.72%). While skipped exons found to be the major type (48.82%). Using these transcriptome data, 386414 SNPs were called, among which transitions were the main type (76.31%). According to GO database, transcriptome SNPs were involved in a wide range of biological functions, and are especially enriched in biological processes and cell components; KEGG pathway analysis shows that the genes of SNPs are enriched in metabolic pathways, endocytosis, apoptosis, cell aging and other related pathways, and more genes are enriched in cancer, signaling pathways and immune system related pathways. Conclusions:Genome transcription information of gingival tissue of patients with impacted maxillary canine showed gingiva were closely related to the development of teeth.

Key words:

single nucleotide polymorphism, ?gingival tissue, ?transcriptome, ?periodontal tissue

陈彬刘畅管晓燕白国辉王苗钱声艳刘建国. 上颌尖牙阻生患者牙龈组织的转录组特征分析[J]. 口腔生物医学, 2021, 12(1): 34-37.

参考文献 23

[1]	Venter IC, Adams MD, Myers EW, et al.The sequence of the human genome[J].Science, 2001, 291(5507):1304-1351
[2]	Lander ES, Linton LM, Birren B, et al.Initial sequencing and analysis ofthe human genome[J].Nature, 2001, 409(6822):860-921
[3]	陈柳, 华清泉.转录组在肿瘤中的应用及进展[J].中国医药导报, 2016, 13(6):58-61
[4]	Ding X, Zhu L, Ji T, et al.Long intergenic non-coding RNAs (LincRNAs) identified by RNA-Seq in breast cancer[J].PLoS One, 2014, 9(8):1-10
[5]	卢倩, 邹丽萍.转录组测序技术在癫痫中的应用[J].癫痫杂志, 2018, 4(2):112-116
[6]	Kukurba KR, Montgomery SB.RNA Sequencing and Analysis[J].Cold Spring Harb Protoc, 2015, 2015(11):951-969
[7]	刘安琪, 钱玉芬.上颌尖牙阻生原因及其对周围结构的影响[J].上海口腔医学, 2017, 26(6):684-687
[8]	赵海燕, 李姣, 丁一.牙被动萌出异常的研究进展[J].国际口腔医学杂志, 2013, 40(1):98-101
[9]	贺佳倩, 刘英群, 马海峰, 等.牙齿萌出过程中的表达与破骨细胞的关系[J].临床口腔医学杂志, 2008, 24(2):72-75
[10]	郭子杰, 全军华, 张铁梅, 等.大鼠牙龈组织表达护骨因子和护骨因子配体[J].现代口腔医学杂志, 2002, 16(4):289-291
[11]	金星爱, 祝莹颖, 薛欣, 等.在犬乳恒牙替换期牙龈上皮中的表达[J].哈尔滨医科大学学报, 2007, 41(4):333-334
[12]	Cock PJ, Fields CJ, Goto N.et alThe Sanger FASTQ file format for sequences with quality scores,and the SolexaIllumina FASTQ variants[J].Nucleic Acids Res, 2010, 38(6):1767-1771
[13]	Kim, D, Langmead, B.Salzberg,SL. .HISAT: a fast spliced aligner with low memory requirements[J].Nat. Methods, 2015, 12(4):357-360
[14]	Trapnell Cole, Roberts Adam, Goff Loyal et al.Differential gene and transcript expression analysis of RNA-seq experiments with TopHat and Cufflinks[J].Nature protocols, 2012, 7(3):562-578
[15]	McKenna A, Hanna M, Banks E et al.The Genome Analysis Toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data[J].Genome Reserch, 2010, 20(9):1297-1303
[16]	Shen S, Park J, Lu L et al.rMATS: Robust and flexible detection of differential alternative splicing from replicate RNA- Seq data[J].Proc. Natl Acad. Sci. USA, 2014, 111(51):5593-5601
[17]	王旭平, 杨胜林, 柳诚刚, 等.番鸭下丘脑组织的-特征分析[J].中国畜牧兽医, 2018, 45(3):604-611
[18]	冉茂良, 陈斌, 李智, 等.基于-测序数据鉴定和分析猪基因组可变剪接事件[J].中国科学:生命科学, 2016, 46(3):274-284
[19]	PAN Q, SHAI O, LEE J L, et al.Deep surveying of alterna tive splicing complexity in the human transcriptome by high throughput sequencing[J].Nature Genetics, 2008, 40(12):1413-1415
[20]	陈黎, 李国勤, 田勇, 等.北京鸭腹部脂肪组织的转录组特征分析[J].浙江农业学报, 2016, 28(5):743-747
[21]	GRAVELEY BR, BROOKS A N, CARLSON J W, et al.The developm -ental transcriptome of Drosophila melanogaster[J].Nature, 2011, 471(7339):473-479
[22]	KIME, MAGENA, ASTG.Different level sofalter native splicing among eukaryotes[J].Nucleic Acids Resech, 2007, 35(1):125-131
[23]	HITZEMANNR, BOTTOMLYD, DARAKJIANP, et al.Gene,behavior and next generation RNA Sequencing[J].Genes, Brain and Behavior, 2012, 12(1):1-12