Transcriptomic analysis of mouse lungs ofter acute infection with Toxoplasma gondii
CAO Wen-qiao1, LI Kai1, YU Yun1, TAN Qing-qing1, LU He1,2, ZHANG Jing1,2
1.Molecular Medicine and Cancer Research Center of Chongqing Medical University, Chongqing 400016, China; 2.Department of pathogenic biology of basic medicine college Chongqing medical university, Chongqing 400016,China
Abstract:The purpose of this study was to investigate the transcriptomic changes in mouse lungs infected with Toxoplasma gondii. Total RNA was extracted from infected and uninfected mouse lung samples, and a cDNA library was constructed from the mRNA, and then subjected to high throughput RNA sequencing (RNA-seq) with the BGISEQ-500 platform. Ten differentially expressed genes (DEGs) were randomly selected for verification with quantitative real-time PCR (q-PCR). Gene Ontology (GO) enrichment analysis, Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis and key driver gene analysis (KDA) were performed on the sequencing data. A total of 286 DEGs were identified, among which 234 were up-regulated and 52 were down-regulated. Q-PCR experiments showed consistent trends with the RNA-seq results, indicating that the RNA-seq data were reliable. Most of the GO terms were associated with immunity and cytokines; KEGG enrichment analysis showed that the terms cytokine-cytokine receptor interaction and chemokine signaling pathway were significantly enriched, i.e., Top1 and Top4, respectively. KDA analysis was performed on the 24 DEGs that were significantly enriched in cytokine receptor interaction and chemokine signaling pathway, and ten key driver genes were identified from the 24 DEGs: Tnf, Cxcl10, Il10, Ccl2, Ifng, Cxcl9, Cxcr2, Ccr5, Ccl4, and Ccl7. Collectively, our research revealed ten key driver genes that may play an important role in mouse lungs after acute Toxoplasma gondii infection.
[1] Liu Q, Wang ZD, Huang SY,et al.Diagnosis of toxoplasmosis and typing of Toxoplasma gondii[J]. ParasitVector, 2015, 8:292.DOI: 10.1186/s13071-015-0902-6 [2] Geita S, Majid G.A review on human toxoplasmosis[J].Scand J Infect Dis, 2012, 44(11): 805-814. DOI: 10.3109/00365548.2012.693197 [3] Schmidt M, Sonneville R, Schnell D, et al.Clinical features and outcomes in patients with disseminated toxoplasmosis admitted to intensive care: a multicenter study[J]. Clin Infect Dis, 2013, 57: 1535-1541. DOI: 10.1093/cid/cit557 [4] Ludlam GB, EdinMD,et al. Pulmonary Toxoplasmosis?[J]. Lancet, 1963, 2(7318):1136-1138.DOI: 10.1016/s0140-6736(63)90789-x [5] Montoya JG, Liesenfeld O.Toxoplasmosis[J]. Lancet, 2004, 363: 1965-1976. DOI: 10.1016/S0140- 6736(04)16412-X [6] Suzuki Y, Orellana M A, Schreiber R D, et al.Interferon-gamma: the major mediator of resistance against Toxoplasma gondii[J]. Science, 1988, 240: 516-518. DOI: 10.1126/science.3128869 [7] López-Yglesias AH, Camanzo E, Martin AT, et al.TLR11-independent inflammasome activation is critical for CD4+ T cell-derived IFN-γ production and host resistance to Toxoplasma gondii[J]. PLoS Pathog, 2019, 15: e1007872. DOI: 10.1371/journal.ppat.1007872 [8] Safronova A, Araujo A,Camanzo ET, et al.Alarmin S100A11 initiates a chemokine response to the human pathogen Toxoplasma gondii[J]. Nat Immunol, 2019, 20: 64-72. DOI: 10.1038/s41590-018-0250-8 [9] Kobayashi K,Umeda K,Ihara F, et al.Transcriptome analysis of the effect of C-C chemokine receptor 5 deficiency on cell response to Toxoplasma gondii in brain cells[J].BMC Genomics, 2019, 20: 705.DOI: 10.1186/s12864-019-6076-4 [10] Burg JL, Grover CM, Pouletty P, et al.Direct and sensitive detection of a pathogenic protozoan, Toxoplasma gondii, by polymerase chain reaction.[J]. J Clin Microbiol, 1989, 27: 1787-1792. [11] Li R,Yu C,Li Y, et al.SOAP2: an improved ultrafast tool for short read alignment[J]. Bioinformatics, 2009, 25: 1966-1973. DOI: 10.1093/bioinformatics/btp336 [12] Kim D,Langmead B,Salzberg SL.HISAT: a fast spliced aligner with low memory requirements[J]. Nat Methods, 2015, 12: 357-360. DOI: 10.1038/nmeth.3317 [13] Langmead B,Salzberg SL.Fast gapped-read alignment with Bowtie 2[J]. Nat Methods, 2012, 9: 357-359. DOI: 10.1038/nmeth.1923 [14] Li B, Dewey CN.RSEM: accurate transcript quantification from RNA-Seq data with orwithout a reference genome[J]. BMC Bioinformatics, 2011, 12: 323. DOI: 10.1186/1471-2105-12-323 [15] Love MI,Huber W,Anders S.Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2[J]. Genome Biol, 2014, 15: 550. DOI: 10.1186/s13059-014-0550-8 [16] Ryan MC, Stucky M, Wakefield C, et al. Interactive Clustered Heat Map Builder: An easy web-based tool for creating sophisticated clustered heat maps[J].F1000 Res,2019,8:ISCB Comm J-1750. DOI: 10.12688/f1000research.20590.2. [17] Bando H,Lee Y,Sakaguchi N,et al.Toxoplasma Effector GRA15-dependent suppression of IFN-γ-induced antiparasitic response in human neurons[J]. Front Cell Infect Microbiol, 2019, 9: 140. DOI: 10.3389/fcimb.2019.00140 [18] Yamamoto M,Okuyama M,Ma JS,et al.A cluster of interferon-γ-inducible p65 GTPases plays a critical role in host defense against Toxoplasma gondii[J]. Immunity, 2012, 37: 302-313. DOI: 10.1016/j.immuni.2012. 06.009 [19] Hunn JP,Koenen-Waisman S,Papic N, et al.Regulatory interactions between IRG resistance GTPases in the cellular response to Toxoplasma gondii[J]. EMBO J, 2008, 27: 2495-2509. DOI: 10.1038/emboj.2008.176 [20] Martens S,Parvanova I,Zerrahn J,et al.Disruption of Toxoplasma gondii parasitophorousvacuoles by the mouse p47-resistance GTPases[J]. PLoS Pathog, 2005, 1: e24. DOI: 10.1371/journal.ppat.0010024 [21] Wen X,Kudo T,Payne L,et al.Predominant interferon-γ-mediated expression of CXCL9, CXCL10, and CCL5 proteins in the brain during chronic infection with Toxoplasma gondii in BALB/c mice resistant to development of toxoplasmicencephalitis[J]. J Interferon Cytokine Res, 2010, 30: 653-660. DOI: 10.1089/ jir.2009.0119 [22] Ochiai E,Sa Q, Brogli M,et al.CXCL9 is important for recruiting immune T cells into the brain and inducing an accumulation of the T cells to the areas of tachyzoite proliferation to prevent reactivation of chronic cerebral infection with Toxoplasma gondii[J]. Am J Pathol, 2015, 185: 314-324. DOI: 10.1016/j.ajpath.2014. 10.003 [23] Norose K,Kikumura A,Luster AD,et al.CXCL10 is required to maintain T-cell populations and to control parasite replication during chronic ocular toxoplasmosis[J]. Invest Ophthalmol Vis Sci, 2011, 52: 389-398. DOI: 10.1167/iovs.10-5819 [24] 李文东, 唐娟, 卓卫. 血清趋化因子CCL5和CXCL10在刚地弓形虫急性感染早期诊断中的应用[J].医学理论与实践, 2019, (18):2879-2882. DOI: 10.19381/j.issn.1001-7585.2019.18.008 [25] He JJ,Ma J,Wang JL,et al.Toxoplasma gondii global transcriptome profiling of multiple porcine organs reveals-induced transcriptional landscapes[J]. Front Immunol, 2019, 10: 1531. DOI: 10.3389/fimmu.2019. 01531 [26] Khan IA,Murphy PM,Casciotti L, et al.Mice lacking the chemokine receptor CCR1 show increased susceptibility to Toxoplasma gondii infection[J]. J Immunol, 2001, 166: 1930-1937. DOI: 10.4049/jimmunol. 166.3.1930 [27] Del Rio L,Bennouna S,Salinas J,et al.CXCR2 deficiency confers impaired neutrophil recruitment and increased susceptibility during Toxoplasma gondii infection[J]. J Immunol, 2001, 167: 6503-6509. DOI: 10.4049/ jimmunol.167.11.6503 [28] BonfáG,Benevides L,Souza Mdo C,et al. CCR5 controls immune and metabolic functions during Toxoplasma gondii infection[J].PLoS One, 2014, 9: e104736. DOI: 10.1371/journal.pone.0104736 [29] Gazzinelli RT,Wysocka M,Hieny S,et al.In the absence of endogenous IL-10, mice acutely infected with Toxoplasma gondii succumb to a lethal immune response dependent on C<inline-formula><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="Mml1-1002-2694-37-3-203"><mml:msup><mml:mrow><mml:msub><mml:mrow><mml:mi mathvariant="normal">D</mml:mi></mml:mrow><mml:mrow><mml:mn>4</mml:mn></mml:mrow></mml:msub></mml:mrow><mml:mrow><mml:mo>+</mml:mo></mml:mrow></mml:msup></mml:math></inline-formula> T cells and accompanied by overproduction of IL-12, IFN-gamma and TNF-alpha[J]. J Immunol, 1996, 157: 798-805.