眼科

• 论著 • 上一篇    下一篇

糖尿病视网膜病变患者转录组学初步研究

刘颖1 何徐军2 陈建斌1 李志琛1 梅伟群1 张华北1 欧阳建1 钱佳丽1   

  1. 1解放军联勤保障部队第903医院内分泌科,杭州 310004; 2 浙江省人民医院省胃肠病学重点实验室,杭州 310013
  • 收稿日期:2019-08-31 出版日期:2020-03-25 发布日期:2020-03-25
  • 通讯作者: 刘颖,Email:liuying661188@163.com
  • 基金资助:
    南京军区医药卫生重点项目(11Z038);浙江省公益技术研究社会发展项目(2012C33120)

A preliminary study of transcriptome in patients with diabetic retinopathy

Liu Ying1, He Xujun2, Chen Jianbin1, Li Zhichen1, Mei Weiqun1, Zhang Huabei1, Ouyang Jian1, Qian Jiali1   

  1. 1Department of Endocrinology, the 903rd Hospital of Joint Logistics Support Force of Chinese PLA, Hangzhou 310004, China; 2 Zhejiang People's Hospital Key Laboratory of Gastroenterology, Hangzhou 310013, China
  • Received:2019-08-31 Online:2020-03-25 Published:2020-03-25
  • Contact: Liu Ying, Email: Liuying 661188@163.com
  • Supported by:
    the Key Medical and Health Projects of Nanjing Military Region (11Z038); the Social Development Projects of Public Welfare Technology Research in Zhejiang Province (2012C33120)

摘要: 目的 利用转录组学技术研究有无糖尿病视网膜病变(DR)患者差异表达基因,寻找DR相关基因。设计 前瞻性比较性病例系列。研究对象 14例2型糖尿病(T2DM)患者,其中 8例合并中度以上非增生性DR(NPDR)或增生性DR(PDR)的T2DM患者纳入DR组,6例不合并DR的T2DM患者纳入DM组。方法 采集每例患者外周血,分离白细胞,抽提总RNA,反转录为cDNA,构建转录组文库,Illumina HiSeqTM2000系统测序,对测序结果进行生物信息学分析。RT-qPCR检测外周血白细胞钙调蛋白1(CNN1)、溶血磷脂酸受体3(LPAR3)mRNA水平,验证转录组测序结果。以两组间差异基因表达量均数和中位数差异倍数上下调均≥2.0倍为标准筛选差异表达基因。主要指标 两组间差异表达基因,差异基因本体(GO)注释和富集、生化代谢与信号转导通路(Pathway)注释和富集、差异基因相互作用网络。两组间外周血白细胞CNN1、LPAR3的mRNA水平差异。结果 与DM组相比,DR组筛选出差异表达基因103个,其中上调41个,下调62个。下调最明显的为杀伤细胞免疫球蛋白样受体2DS1(KIR2DS1),上调最明显的为U105B小核仁RNA基因(RNU105B)。GO富集分析显示在生物学过程(GO-P)、分子功能(GO-F)中显著富集的GO条目为细胞过程、单一生物过程、结合功能。谷胱甘肽S转移酶mu1(GSTM1)、早期生长应答1(EGR-1)、LPAR3、CNN1在以上条目均有富集。Pathway显著性富集分析显示,差异基因富集最多的为抗原处理和呈现、自然杀伤细胞介导的细胞毒性通路,KIR2DS1均参与其中。差异基因相互作用网络分析示,差异表达基因中GAS1与多个基因之间存在相互作用。RT-qPCR结果显示,与DM组比较,DR组中的CNN1和 LPAR3的mRNA含量均升高(P均<0.05)。结论 采用转录组学技术发现,DR组患者外周血白细胞存在多种差异表达基因,下调的基因中KIR2DS1、GSTM1、EGR-1、GAS1,上调的基因中RNU105B、LPAR3、CNN1是否与DR的发生发展有关有待于进一步验证。(眼科, 2020, 29: 105-113)

关键词: 糖尿病视网膜病变, 转录组学, 基因

Abstract: Objective To study the differentially expressed genes and DR-related genes in patients with diabetic retinopathy (DR) by using transcriptome technology. Design Prospective comparative case series. Participants Fourteen patients with T2DM were recuited, among them, eight patients with moderate or higher non-proliferative DR (NPDR) or proliferative DR (PDR) were enrolled in the DR group and six patients without DR were enrolled in the DM group. Methods Peripheral blood of each patient was collected, white blood cells were isolated, total RNA was extracted, and retranscribed to DNA. Transcription group library was constructed. Illumina HiSeqTM2000 system was used to sequence the results. Bioinformatics analysis was used to analyze the sequencing results. RT-q PCR was used to detect the levels of calmodulin 1 (CNN1) and lysophosphatidic acid receptor 3 (LPAR3) in peripheral blood leukocytes. To validate transcriptome sequencing results. Main Outcome Measures Differentially expressed genes, GO annotation and enrichment, biochemical metabolism and pathway annotation and enrichment, and differential gene interaction network between the two groups. The levels of CNN1 and LPAR3 in peripheral leukocytes were different between the two groups. Results Compared with DM group, DR group screened 103 differentially expressed genes. Among them, 41 were up-regulated and 62 were down-regulated. The most significant down-regulation was killer cell immunoglobulin like receptor 2DS1 (KIR2DS1), and the most significant up-regulation was small nucleolar RNA gene (RNU105B). GO enrichment analysis showed that the GO items significantly enriched in biological process (GO-P) and molecular function (GO-F) were cell process, single biological process and binding function. Glutathione S-transferase mu1 (GSTM1), early growth response 1 (EGR-1), LPAR3 and CNN1 were all enriched in the above items. The analysis of significance enrichment of pathway showed that the most abundant differential genes were antigen treatment and presentation, natural killer cell-mediated cytotoxicity pathway, in which KIR2DS1 was involved. The analysis of gene interaction network showed that there was interaction between GAS1 and multiple genes. qPCR results indicated that the levels of CNN1 and LPAR3 in DR group were higher than those in DM group (all P<0.05). Conclusion Using transcriptome technology, we found that there were many differential expression genes in peripheral blood leukocytes of DR group. Among them, whether KIR2DS1, GSTM1, EGR-1, and GAS1 in the down regulated genes, and RNU105B, LPAR3, and CNN1 in the up regulated genes are related to the occurrence and development of DR remains to be further verified. (Ophthalmol CHN, 2020, 29: 105-113)

Key words: diabetic retinopathy, transcriptome, gene