中国寄生虫学与寄生虫病杂志 ›› 2023, Vol. 41 ›› Issue (3): 270-278.doi: 10.12140/j.issn.1000-7423.2023.03.002

• 论著 • 上一篇    下一篇

弓形虫慢性感染小鼠脑转录组分析及与抑郁相关的犬尿氨酸通路的验证

张驰(), 陈嘉婷, 辛紫萱, 杨莉莉, 杨梓瀚, 彭鸿娟*()   

  1. 南方医科大学公共卫生学院病原生物学系,广东省热带病研究重点实验室,广州 510515
  • 收稿日期:2022-11-08 修回日期:2023-02-01 出版日期:2023-06-30 发布日期:2023-06-15
  • 通讯作者: *彭鸿娟(1971-),女,博士,教授,从事病原体-宿主相互作用机制研究。E-mail:hongjuan@smu.edu.cn
  • 作者简介:张驰(1995-),男,硕士研究生,从事病原生物学研究。E-mail:zhangchismu@126.com
  • 基金资助:
    国家自然科学基金面上项目(82272364);国家自然科学基金面上项目(81971954);广州市重点实验室基础研究计划(202102100001)

Transcriptome analysis of mice brain chronically infected with Toxoplasma gondii and validation of the kynurenine pathway associated with depression

ZHANG Chi(), CHEN Jiating, XIN Zixuan, YANG Lili, YANG Zihan, PENG Hongjuan*()   

  1. Department of Pathogen Biology, School of Public Health, Guangdong Provincial Key Laboratory of Tropical Medicine, Southern Medical University, Guangzhou 510515, China
  • Received:2022-11-08 Revised:2023-02-01 Online:2023-06-30 Published:2023-06-15
  • Contact: *E-mail: hongjuan@smu.edu.cn
  • Supported by:
    General Program of National Natural Science Foundation of China(82272364);General Program of National Natural Science Foundation of China(81971954);Basic Research Program of Guangzhou Key Laboratory(202102100001)

摘要:

目的 筛选并分析刚地弓形虫慢性感染小鼠脑转录组差异表达基因(DEG),分析与抑郁相关的犬尿氨酸(KYN)通路DEG的相对转录水平,为探究弓形虫慢性感染导致小鼠抑郁样症状的机制提供理论依据。方法 18只SV129雄性小鼠随机平均分为感染组和对照组。感染组小鼠腹腔注射ME49株速殖子120个(200 μl),对照组注射等体积的磷酸缓冲液,感染后3个月收集感染组和对照组小鼠脑组织,提取小鼠脑组织总RNA进行转录组测序,筛选DEG,对DEG进行聚类分析、基因本体(GO)功能注释分析、京都基因与基因组百科全书(KEGG)功能注释和富集分析。选取与抑郁症相关的KYN通路的8个DEG,分别为γ干扰素(IFN-γ)、吲哚胺2,3-双加氧酶1(IDO1)、IDO2、犬尿氨酸酶(KYNU)、犬尿氨酸-3-单氧化酶(KMO)、3-羟基邻氨基苯甲酸3,4-二加氧酶(3-HAO)、波形蛋白(Vim)和脑源性神经营养因子(BDNF),以甘油醛-3-磷酸脱氢酶基因为内参,实时荧光定量PCR(qRT-PCR)检测各个基因的相对转录水平。结果 感染组和对照组小鼠脑转录组的DEG共2 295个,其中上调2 016个,下调279个。GO分析结果显示,生物过程中富集最显著的是定位,共257个DEG;细胞组分中富集最显著的是蛋白复合物,共425个DEG;分子功能中富集最显著的是分子转导活性,共177个DEG。生物过程、细胞组分和分子功能富集DEG数量最多的分别是细胞过程、细胞组分和结合,分别有1 039、1 240、1 088个DEG。KEGG分析结果显示,功能注释分析上调居前3位的代谢通路分别为免疫系统、信号转导、病毒性传染病,下调居前3位的分别为信号转导、信号分子和相互作用、免疫系统;功能富集分析结果显示77条通路富集显著。与抑郁症相关的信号通路有肿瘤坏死因子、神经活性配体-受体相互作用、NF-kappa B、JAK-STAT、坏死性凋亡、细胞凋亡、趋化因子、KYN通路等。qRT-PCR结果显示,以对照组小鼠相对转录水平为100%,感染组小鼠IFN-γ、IDO1、IDO2、KYNU、KMO、3-HAO和Vim等7个基因的相对转录水平分别为3 023.08%、355.52%、190.17%、496.55%、339.92%、212.74%、507.34%,较对照组的转录水平显著上调(t = 3.782、3.749、3.226、2.908、2.533、5.656、2.948,均P < 0.05或0.01);BDNF的相对转录水平为63.32%,转录水平显著下调(t = 2.398,P < 0.05)。IFN-γ、IDO1、IDO2、KYNU、KMO、3-HAO、BDNF、Vim等8个基因qRT-PCR获得的差异倍数分别为4.96、1.74、0.89、2.10、1.60、1.06、-0.94、2.18,转录组测序获得的差异倍数分别为7.30、0.55、0.80、3.83、2.75、3.53、-0.86、1.93。qRT-PCR与转录组测序获得的转录趋势一致。结论 筛选获得刚地弓形虫慢性感染小鼠脑转录组DEG,弓形虫慢性感染小鼠中枢神经系统免疫持续激活,与抑郁症相关的KYN通路的7个DEG的转录水平上调。

关键词: 刚地弓形虫, 小鼠脑组织, 转录组, 抑郁症, 犬尿氨酸通路

Abstract:

Objective To screen the differentially expressed genes (DEGs) by comparing the transcriptome of the brain tissues between the mice chronically infected with Toxoplasma gondii and normal mice, to analyze the relative transcription level of DEGs in the depression-related kynurenine (KYN) pathway and to provide a theoretical basis for exploring the mechanism of depression-like symptoms caused by Toxoplasma gondii chronic infecttion in mice. Methods SV129 male mice (n = 18) were randomly and equally divided into the infection group and the control group. Mice in the infection group were intraperitoneally injected with 120 tachyzoites of T. gondii ME49 strain (200 μl), and mice in the control group were injected with the same volume of PBS. After 3 months post-infection, mice brain tissues of the two groups were collected for extraction of total RNA to undertake transcriptome sequencing for screening DEG. With the DEGs obtained, cluster analysis, gene ontology (GO) functional annotation analysis, Kyoto Encyclopedia of Genes and Genomes (KEGG) functional annotation and enrichment analysis were performed. Eight DEGs [interferon-γ (IFN-γ), indoleamine 2,3-dioxygenase 1 (IDO1), IDO2, kynurenine-3-monooxygenase (KYNU), kynurenine-3-monooxygenase (KMO), 3-hydroxyanthranilate 3,4-dioxygenase (3-HAO), vimentin (Vim) and brain-derived neurotrophic factor (BDNF)] related to KYN pathway associated with depression were selected to examine each gene’s relative transcription level by quantitative real-time PCR (qRT-PCR), using glyceraldehyde-3-phosphate dehydrogenase gene as an internal reference. Results Transcriptome sequencing found 2 295 DEGs in the brain of the mice from the infection and control groups, of which 2 016 were up-regulated and 279 were down-regulated. GO analysis showed that localisation was the most significantly enriched biological process, with a total of 257 DEGs. The most significantly enriched in cellular components was the protein-containing complex, with a total of 425 DEGs. The most significantly enriched molecular function was molecular transducer activity, with 177 DEGs. The largest number of DEGs enriched in biological process, cell component and molecular function were cell process, cell part and binding, with 1 039, 1 240 and 1 088 DEGs, respectively. KEGG analysis showed that the top three up-regulated metabolic pathways were the immune system, signaling transduction, and viral infectious disease, and the top three down-regulated pathways were signal transduction, signaling molecules and interaction and immune system. Functional enrichment analysis showed that 77 pathways were significantly enriched. The signaling pathways related to depression include tumor necrosis factor signaling pathway, neuroactive ligand-receptor interaction, NF-kappa B signaling pathway, JAK-STAT signaling pathway, necroptosis, apoptosis, chemokine signaling pathway, KYN pathway. The qRT-PCR results showed that the relative transcription levels of IFN-γ, IDO1, IDO2, KYNU, KMO, 3-HAO and Vim genes in the infection group were 3 023.08%, 355.52%, 190.17%, 496.55%, 339.92%, 212.74% and 507.34%, if the relative transcript level of control mice was taken as 100%. Compared with the control group, the transcription was significantly up-regulated (t = 3.782, 3.749, 3.226, 2.908, 2.533, 5.656, 2.948; all P < 0.05 or 0.01). The relative transcription level of BDNF was 63.32%, which was significantly down-regulated (t = 2.398, P < 0.05). The fold change of IFN-γ, IDO1, IDO2, KYNU, KMO, 3-HAO, BDNF, Vim obtained by qRT-PCR was 4.96, 1.74, 0.89, 2.10, 1.60, 1.06, -0.94, 2.18, respectively. The fold change obtained by transcriptome sequencing was 7.30, 0.55, 0.80, 3.83, 2.75, 3.53, -0.86 and 1.93, respectively. The transcriptional trend obtained by qRT-PCR was consistent with that obtained by transcriptome sequencing. Conclusion DEGs from brain tissues of mice chronically infected with T. gondii were screened. Transcriptome analysis revealed that the immune response of central nervous system of the mice with chronic infection of T. gondii was continuously activated. Seven DEGs in KYN pathway related to depression showed up-regulated transcription level.

Key words: Toxoplasma gondii, Mice brain tissue, Transcriptome, Depression, Kynurenine pathway

中图分类号: