Transcriptome analysis of human monocytic THP-1 cells overexpressing Toxoplasma gondii rhoptry protein 16

doi:10.12140/j.issn.1000-7423.2025.03.002

Abstract

Abstract:

Objective To investigate the changes in the gene expression profile of human monocyte THP-1 cells overexpressing Toxoplasma gondii rhoptry protein 16 (ROP16) using transcriptome sequencing, and to screen immune response-related driver genes. Methods Human monocytic THP-1 cells were seeded onto 96-well plates. Cells in the overexpression group was transfected with ROP16 overexpression lentivirus for 72 h, and cells in the control group were treated with an equal volume of culture medium. The efficiency of cell transfection was checked under a fluorescence microscope. Total RNA was extracted from THP-1 cells and THP-1 cells stably expressing ROP16 using the TRIzol reagent, followed by transcriptome sequencing. Differentially expressed genes (DEG) were screened, and a volcano plot was generated. DEGs were subjected to cluster analysis, Kyoto encyclopedia of genes and genomes (KEGG) metabolic pathway enrichment analysis, and gene ontology (GO) functional enrichment classification to screen immune response-related driver genes in human monocytic THP-1 cells following T. gondii infections, and gene expression was quantified using qPCR assay. Differences of means between groups were tested for statistical significance with independent sample t-test. Results Fluorescence microscopy displayed green fluorescence in more than 90% of the cells in the overexpression group in each field of view, and no green fluorescence was observed in the control group. qPCR assay quantified a higher relative expression level of ROP16 mRNA in the overexpression group (1 083.484 ± 68.990) than in the control group (1.000 ± 0) (t = 22.9, P < 0.01), indicating the successful generation of THP-1 cells that stably expressed ROP16. A total of 312 DEGs were identified in THP-1 cells stably expressing ROP16, including 193 upregulated genes and 119 downregulated genes. KEGG annotations showed that the highest proportion of DEGs were annotated to the organism system (24.2%), with 85 items significantly enriched, among which 47 genes were significantly enriched in the immune system. KEGG enrichment analysis showed that DEGs were significantly enriched in 20 signaling pathways, and 6, 7, and 8 DEG were significantly enriched in three pathways related to the immune system, including Th1 and Th2 cell differentiation, natural killer cell-mediated cytotoxicity, and IL-17 signaling pathways, respectively. GO functional annotations showed that a total of 309 DEGs were annotated to 55 secondary node classifications under three primary node classifications of biological processes, cellular components, and molecular functions. GO enrichment analysis showed that DEGs were significantly enriched in inflammatory response, negative regulation of tumor necrosis factor products, cytokine production, positive regulation of interferon-γ production, and immune response-related pathways. RT-qPCR assay detected that higher relative expression of MAN2B1, FOS, C1QA mRNA (25.994 ± 0.382、60.584 ± 2.968、 36.759 ± 0.180) in cells in the ROP16 overexpression group than in the control group (1.000 ± 0.039、1.000 ± 0.015、1.000 ± 0) (t = 92.00, 28.39, 280.7, P < 0.05 or 0.01), and lower relative expression of LMNB1, IL-6, and IL-12 mRNA (0.728 ± 0.054, 0.517 ± 0.073, 0.587 ± 0.015) in the overexpression group than in the control group (1.052 ± 0.027、1.000 ± 0.039、1.000 ± 0.010) (t = 7.64, 8.24, 33.62, P < 0.05 or 0.01). Conclusion The gene expression profile of human monocyte THP-1 cells changes significantly and may play important roles in the immune response following T. gondii infection.

Key words: Toxoplasma gondii, Rhoptry protein 16, Human monocytic THP-1 cell, Transcriptome, Immune response

CLC Number:

R382.5

LI Jiaming, CHEN Mei, DANG Tiantian, YIN He, ZHAO Zhijun. Transcriptome analysis of human monocytic THP-1 cells overexpressing Toxoplasma gondii rhoptry protein 16[J]. CHINESE JOURNAL OF PARASITOLOGY AND PARASITIC DISEASES, 2025, 43(3): 317-323.

Figures/Tables 5

Table 1

Fig. 1

Fig. 2

Fig. 3

Table 2

References 22

[1]	Parlog A, Schlüter D, Dunay IR. Toxoplasma gondii-induced neuronal alterations[J]. Parasite Immunol, 2015, 37(3): 159-170.
[2]	Marín-García PJ, Planas N, Llobat L. Toxoplasma gondii in foods: Prevalence, control, and safety[J]. Foods, 2022, 11(16): 2542.
[3]	Fox BA, Bzik DJ. Nonreplicating, cyst-defective type Ⅱ Toxoplasma gondii vaccine strains stimulate protective immunity against acute and chronic infection[J]. Infect Immun, 2015, 83(5): 2148-2155.
[4]	Ye HM, Zhou XT, Zhu BK, et al. Toxoplasma gondii suppresses proliferation and migration of breast cancer cells by regulating their transcriptome[J]. Cancer Cell Int, 2024, 24(1): 144.
[5]	Du KG, Lu F, Xie CZ, et al. Toxoplasma gondii infection induces cell apoptosis via multiple pathways revealed by transcriptome analysis[J]. J Zhejiang Univ Sci B, 2022, 23(4): 315-327.
[6]	李佳铭, 王艺璇, 杨宁爱, 等. 刚地弓形虫ROP16蛋白对MH-S细胞极化和凋亡的影响及其相关机制[J]. 中国寄生虫学与寄生虫病杂志, 2022, 40(5): 579-586.
	Li JM, Wang YX, Yang NA, et al. Effects of ROP16 protein of Toxoplasma gondii on polarization and apoptosis of MH-S cells and their related mechanisms[J]. Chin J Parasitol Parasit Dis, 2022, 40(5): 579-586. (in Chinese)
[7]	Zhu LJ, Qi WJ, Yang G, et al. Toxoplasma gondii rhoptry protein 7 (ROP7) interacts with NLRP3 and promotes inflammasome hyperactivation in THP-1-derived macrophages[J]. Cells, 2022, 11(10): 1630.
[8]	Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2[J]. Genome Biol, 2014, 15(12): 550.
[9]	Saeij JJ, Coller S, Boyle JP, et al. Toxoplasma co-opts host gene expression by injection of a polymorphic kinase homologue[J]. Nature, 2007, 445(7125): 324-327.
[10]	Hu RS, He JJ, Elsheikha HM, et al. Transcriptomic profiling of mouse brain during acute and chronic infections by Toxoplasma gondii oocysts[J]. Front Microbiol, 2020, 11: 570903.
[11]	Yuan H, Zhang XX, Yang ZP, et al. Unveiling of brain transcriptome of masked palm civet (Paguma larvata) with chronic infection of Toxoplasma gondii[J]. Parasit Vectors, 2022, 15(1): 263.
[12]	Zhu LQ, Lei ZW, Xia XC, et al. Yeast shells encapsulating adjuvant AS04 as an antigen delivery system for a novel vaccine against Toxoplasma gondii[J]. ACS Appl Mater Interfaces, 2021, 13(34): 40415-40428.
[13]	Sprenkeler EGG, Zandstra J, van Kleef ND, et al. S100A8/A9 is a marker for the release of neutrophil extracellular traps and induces neutrophil activation[J]. Cells, 2022, 11(2): 236.
[14]	Wang S, Song R, Wang Z, et al. S100A8/A9 in inflammation[J]. Front Immunol, 2018, 9: 1298.
[15]	Purves-Tyson TD, Robinson K, Brown AM, et al. Increased macrophages and C1qA, C3, C4 transcripts in the midbrain of people with schizophrenia[J]. Front Immunol, 2020, 11: 2002.
[16]	Chen LH, Liu JF, Lu Y, et al. Complement C1q (C1qA, C1qB, and C1qC) may be a potential prognostic factor and an index of tumor microenvironment remodeling in osteosarcoma[J]. Front Oncol, 2021, 11: 642144.
[17]	高路, 姚瑞, 李亚彭, 等. 溶酶体组织蛋白酶B增加自噬保护缺氧诱导的心脏微血管内皮细胞损伤[J]. 中国药理学通报, 2022, 38(1): 53-60.
	Gao L, Yao R, Li YP, et al. Role of lysosomal cathepsin B in endothelial cell injury induced by hypoxia[J]. Chin Pharmacol Bull, 2022, 38(1): 53-60. (in Chinese)
[18]	Weiss-Sadan T, Maimoun D, Oelschlagel D, et al. Cathepsins drive anti-inflammatory activity by regulating autophagy and mitochondrial dynamics in macrophage foam cells[J]. Cell Physiol Biochem, 2019, 53(3): 550-572.
[19]	Wenzel TJ, Klegeris A. Novel multi-target directed ligand-based strategies for reducing neuroinflammation in Alzheimer’s disease[J]. Life Sci, 2018, 207: 314-322.
[20]	金宇扬, 陈光亮, 陈晓翔. Fosl2与免疫系统和自身免疫性疾病关系的研究进展[J]. 现代免疫学, 2018, 38(6): 509-512.
	Jin YY, Chen GL, Chen XX. Research progress on the relationship between Fosl2 and immune system and autoimmune diseases[J]. Curr Immunol, 2018, 38(6): 509-512. (in Chinese)
[21]	Finnsson J, Lubberink M, Savitcheva I, et al. Glucose metabolism in the brain in LMNB1-related autosomal dominant leukodystrophy[J]. Acta Neurol Scand, 2019, 139(2): 135-142.
[22]	Lin XL, Liu HW, Zhao HY, et al. Immune infiltration associated MAN2B1 is a novel prognostic biomarker for glioma[J]. Front Oncol, 2022, 12: 842973.

基因名称 Gene name	引物序列（5'→3'） Primer sequence（5'→3'）
甘油醛-3-磷酸脱氢酶GAPDH	F：CAGGAGGCATTGCTGATGAT
甘油醛-3-磷酸脱氢酶GAPDH	R：GAAGGCTGGGGCTCATTT
棒状体蛋白 ROP-16	F：TGGGCTCCTGAACTTGCGAAATC
棒状体蛋白 ROP-16	R：AGACGAACTCGAAGATTGCCAACC
组织蛋白酶B CTSB	F：GGTGGCCTCTAYGAATCCCA
组织蛋白酶B CTSB	R：GCCAGGTCACAGATCTTGCTA
原癌基因 Fos FOS	F：GACTGATACACTCCAAGCGG
原癌基因 Fos FOS	R：CAGGGATCTTGCAGGCAGGT
核纤层蛋白B1 LMNB1	F：CTGAGGCGGAATGAGAGATGC
核纤层蛋白B1 LMNB1	R：CTGAGGCGGAATGAGAGATGC
甘露糖苷酶α 2B成员1 MAN2B1	E：GGAGGCGGAATGAGAGATGC
甘露糖苷酶α 2B成员1 MAN2B1	R：CTGTCTGTGTAGAAGCGTCCC
补体C1q B链 C1QB	F：ACAGTAGGGCTTGGTGAATG
补体C1q B链 C1QB	R：GCTGGGACACAGAGTTGAAT
组织蛋白酶A CTSA	F：CATCTTCACTCGCCTGCCAC
组织蛋白酶A CTSA	R：CGGGTTGTTGAGGTAGGTGG
S100钙结合蛋白A8 S100A8	F：CTCTATCATCGACGTCTACC
S100钙结合蛋白A8 S100A8	R：ATCATGGAGGACCTGGACAC
S100A9	F：ATCATGGAGGACCTGGACAC
S100A9	R：ACCCTCGTGCATCTTCTCGT
C1QA	F：CCGGCTACTACTACTTCACC
C1QA	R：CCTTGTTGGTGGTGTCACAG
C1QC	F：TGCACCCAACAGCCTGATCA
C1QC	R：GCGACGCGTGGTAGACAAAG
白细胞介素1β IL-1β	F：ATGGCTTATTACAGTGGCAATGAG
白细胞介素1β IL-1β	R：GTAGTGGTGG TCGGAGATTCG
IL-18	F：GCTGAAGATGATGAAAACCTGG
IL-18	R：CAAATAGAGGCCGATTTCCTTG
IL-10	F：GCCAAGCCTTGTCTGAGATGATCC
IL-10	R：AATCGATGACAGCGCCGTAGC
IL-12	F：GAGGTCCCTCCAAACCGTTGTC
IL-12	R：GATGTAATAGTCCCATCCTTCTTT
IL-6	F：GCCAGAGCTGTGCAGATGAGTAC
IL-6	R：AGTGGCATTGCATCCCTGAGTTG
IL-17	F：TCCCCATCCAGCAAGAGATCCTG
IL-17	R：AGCCCACGGACACCAGTATCTTC
肿瘤坏死因子α TNF-α	F：TCTCCTTCCTGATCGTGGCA
肿瘤坏死因子α TNF-α	R：GGCTACAGGCTTGTCACTCG
转化生长因子β TGF-β	F：AGCAACAATTCCTGGCGATACCTC
转化生长因子β TGF-β	R：TCAACCACTGCCGCACAACTC

基因 Gene	ROP16过表达组 ROP16 overexpression group	对照组 Control group
组织蛋白酶B CTSB	32.226 ± 0.632	1.000 ± 0.015
原癌基因Fos FOS	60.584 ± 2.968	1.000 ± 0.015
核纤层蛋白B1 LMNB1	0.728 ± 0.054	1.052 ± 0.027
甘露糖苷酶α 2B成员 MAN2B1	25.994 ± 0.382	1.000 ± 0.039
补体C1q B链 C1QB	52.166 ± 0.511	1.000 ± 0
C1QA	36.759 ± 0.180	1.000 ± 0
CTSA	3.580 ± 0.035	1.001 ± 0.044
S100钙结合蛋白A8 S100A8	13.693 ± 0.403	1.003 ± 0.108
S100A9	41.434 ± 3.651	1.001 ± 0.064
C1QC	23.837 ± 0.467	1.001 ± 0.059
白细胞介素1β IL-1β	2.462 ± 0.277	1.000 ± 0.010
IL-18	1.454 ± 0.036	1.002 ± 0.093
IL-10	1.718 ± 0.076	1.000 ± 0
IL-12	0.587 ± 0.014	1.000 ± 0.010
IL-6	0.517 ± 0.073	1.000 ± 0.039
IL-17	1.366 ± 0.007	1.000 ± 0.025
肿瘤坏死因子α TNF-α	1.235 ± 0.012	1.000 ± 0.014
转化生长因子β TGF-β	14.124 ± 0.138	1.016 ± 0.251