弓形虫毒力相关效应分子ROP16Ⅲ诱导人肺腺癌A549细胞的基因表达谱

doi:10.12140/j.issn.1000-7423.2021.04.009

中国寄生虫学与寄生虫病杂志 ›› 2021, Vol. 39 ›› Issue (4): 473-480.doi: 10.12140/j.issn.1000-7423.2021.04.009

弓形虫毒力相关效应分子ROP16Ⅲ诱导人肺腺癌A549细胞的基因表达谱

苏雅静¹(), 乔霞¹, 汪澎涛¹, 康宇婷¹, 杨宁爱¹, 贾伟¹^,², 赵志军¹^,²^,^*()

1 宁夏医科大学总医院宁夏临床病原微生物重点实验室,银川 750004
2 宁夏医科大学总医院医学实验中心,银川 750004

收稿日期:2020-10-05 修回日期:2020-11-07 出版日期:2021-08-30 发布日期:2021-06-11
通讯作者: 赵志军
作者简介:苏雅静（1989-）,女,硕士,助理研究员,从事弓形虫免疫研究。E-mail： suyajing8022664@126.com
基金资助:
宁夏自然科学基金(2019AAC03227);第四批宁夏青年科技人才托举工程项目(TJGC2019089)

Gene expression profiling of human lung adenocarcinoma of A549 cells induced by Toxoplasma virulence-related effector ROP16Ⅲ

SU Ya-jing¹(), QIAO Xia¹, WANG Peng-tao¹, KANG Yu-ting¹, YANG Ning-ai¹, JIA Wei¹^,², ZHAO Zhi-jun¹^,²^,^*()

1 Ningxia Key Laboratory of Clinical and Pathogenic Microbiology, General Hospital of Ningxia Medical University, Yinchuan 750004, China
2 Laboratory Center of Medicine, General Hospital of Ningxia Medical University, Yinchuan 750004, China

Received:2020-10-05 Revised:2020-11-07 Online:2021-08-30 Published:2021-06-11
Contact: ZHAO Zhi-jun
Supported by:
Natural Science Foundation of Ningxia(2019AAC03227);Forth Batch of Ningxia Young Scientific Talents Supporting Program(TJGC2019089)

摘要/Abstract

摘要：

目的利用生物信息学方法分析刚地弓形虫（Toxoplasma gondii）棒状体蛋白16Ⅲ （ROP16Ⅲ）重组表达载体转染人肺腺癌A549细胞后的基因表达谱,筛选致病候选基因。方法构建真核表达载体pEGFP-N1-ROP16Ⅲ。将A549细胞分为未转染组、pEGFP-N1组、pEGFP-N1-ROP16Ⅲ组,后两组分别用pEGFP-N1、pEGFP-N1-ROP16Ⅲ质粒进行瞬时转染,48 h后收集各组细胞,提取总RNA,纯化后进行芯片杂交分析,筛选差异表达基因（DEG）。将DEG进行基因本体（GO）功能富集分类和京都基因与基因组百科全书(KEGG)通路分析。利用ToppGene从DEG中筛选弓形虫ROP16Ⅲ致病候选基因。利用STRING在线分析软件对候选基因所编码的蛋白质进行蛋白间相互作用(PPI）网络分析。结果 pEGFP-N1-ROP16Ⅲ转染A549细胞后共有566个基因表达差异显著,其中382个基因上调,184个基因下调。对DEG进行生物信息学分析,共获得997条GO注释和50条KEGG信号通路。这些DEG参与了多种生物学过程,包括防御反应、I型干扰素信号通路、先天免疫应答、对外界生物刺激的反应等;参与的信号通路主要集中在抗原处理及呈递,细胞质DNA传感通路,Toll样受体信号通路,NOD样受体信号通路等。ToppGene筛选共获得14个弓形虫ROP16Ⅲ致病候选基因,其中C-X-C基序趋化因子配体11（CXCL11）、Toll样受体3（TLR3）、C-C基序趋化因子26（CCL26）、人类白细胞抗原-E（HLA-E）、早幼粒白血病基因（PML）及信号转导子与转录激活因子2（STAT2）等6个基因在弓形虫方面的研究尚未见报道。PPI网络构建分析发现,DEG编码的蛋白质PPI网络由376个节点（蛋白质）和1 152条边（相互作用关系）组成,STAT2、HLA-E和PML位于PPI网络的中心节点处,删除这些关键节点后PPI网络结构涣散。通路富集分析结果显示,CXCL11等6个基因共涉及12条信号通路,主要与趋化因子信号通路、Toll样受体信号通路、程序性坏死、内吞作用及细胞因子-细胞因子受体相互作用等相关,这些生物信号通路可能参与弓形虫ROP16Ⅲ的致病过程。结论本研究筛选发现CXCL11、TLR3、CCL26、HLA-E、PML和STAT2致病候选基因,为进一步探索刚地弓形虫病的致病机制提供了理论依据。

关键词: 弓形虫ROP16Ⅲ, 表达谱, 生物信息学, 蛋白互作网络

Abstract:

Objective To analyze the gene expression profile in the A549 human lung adenocarcinoma cell line transfected with the expression vector of Toxoplasma gondii rhoptry protein 16Ⅲ (ROP16Ⅲ), using bioinformatics methods, in order to identify candidate pathogenic genes. Methods The eukaryotic expression vector pEGFP-N1-ROP16Ⅲ was constructed. The A549 cells were assigned to three groups: untransfected, pEGFP-N1 and pEGFP-N1-ROP16Ⅲ groups. Cells in the latter two groups were instantantaneously transfected with pEGFP-N1 and pEGFP-N1-ROP16 plasmid, respectively. After 48 h, the cells were harvested and RNA was extracted. After RNA purification, microarray hybridization analysis was performed to identified differentially expressed genes (DEGs) which were further analyzed by Gene Ontology(GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) for functional enrichment classification and pathway analysis. ToppGene was used to screen for Toxoplasma ROP16Ⅲ pathogenic candidate genes from DEGs. The STRING online analysis software was used to analyze the protein-protein interaction (PPI) network for protein products of the candidate genes. Results A total of 566 genes were differentially expressed in A549 transfected with Toxoplasma ROP16Ⅲ, of which 382 genes were up-regulated and 184 genes down-regulated. Bioinformatics analysis of the DEGs revealed 997 GO annotations and 50 KEGG signaling pathways. The DEGs were significantly enriched in defense response, type I interferon signaling pathway, innate immune response, and response to external biotic stimulus. KEGG pathway enrichment analysis suggested the DEGs were most significantly enriched in antigen processing and presentation, cytosolic DNA-sensing pathway, Toll-like receptor signaling pathway, and NOD-like receptor signaling pathway. In addition, 14 candidate genes were screened out from the DEGs by the ToppGene Suite. Among them, C-X-C motif chemokine 11 (CXCL11), toll-like receptor 3 (TLR3), C-C motif chemokine ligand 26 (CCL26), human leukocyte antigen E(HLA-E), promyelocytic leukemia protein (PML) and signal transducer and activator of transcription 2 (STAT2) had not been reported in Toxoplasma studies. The PPI network of the DEGs comprised 376 nodes and 1 152 edges, in which STAT2, HLA-E and PML were just located at the core position of the PPI network. Deletion of these key nodes led to the PPI network structure disorganized. Pathway enrichment analysis of 6 genes including CXCL11 resulted in 12 signaling pathways, which were mainly related to chemokine signaling pathway, Toll-like receptor signaling pathway, necroptosis, endocytosis and cytokine-cytokine receptor interaction. These bio-signal pathway may be related to the pathogenic mechanism of Toxoplasma ROP16Ⅲ. Conclusion CXCL11, TLR3, CCL26, HLA-E, PML and STAT2 were identified as candidate pathogenic genes, providing theoretical basis for further investigation of pathogenic mechanism of toxoplasmosis.

Key words: Toxoplasma ROP16Ⅲ, Gene expression profiling, Bioinformatics, Protein-protein interaction network

中图分类号:

R382.5

苏雅静, 乔霞, 汪澎涛, 康宇婷, 杨宁爱, 贾伟, 赵志军. 弓形虫毒力相关效应分子ROP16Ⅲ诱导人肺腺癌A549细胞的基因表达谱[J]. 中国寄生虫学与寄生虫病杂志, 2021, 39(4): 473-480.

SU Ya-jing, QIAO Xia, WANG Peng-tao, KANG Yu-ting, YANG Ning-ai, JIA Wei, ZHAO Zhi-jun. Gene expression profiling of human lung adenocarcinoma of A549 cells induced by Toxoplasma virulence-related effector ROP16Ⅲ[J]. Chinese Journal of Parasitology and Parasitic Diseases, 2021, 39(4): 473-480.

图/表 4

表1

pEGFP-N1-ROP16Ⅲ转染的人肺腺癌A549细胞FC值（上调/下调）表达量最大的部分基因

基因标志Gene symbol	基因名称Gene name	变化倍数Fold change	表达情况Expression
CYP2C18	细胞色素P450家族2亚家族C成员18 Cytochrome P450 family 2 subf amily C member 18	18.70	上调 Up
CMPK2	胞苷/尿苷酸激酶2 Cytidine/uridine monophosphate kinase 2	18.40	上调 Up
IFITM1	干扰素诱导的跨膜蛋白1 Interferon induced transmembrane protein 1	17.41	上调 Up
IFIT2	三角形四肽重复干扰素诱导蛋白2 Interferon induced protein with tetratricopeptide repeats 2	17.17	上调 Up
TRIM22	三结构域分子22 Tripartite motif-containing 22	14.88	上调 Up
OAS2	2′,5′-寡腺苷酸合成酶 2 2′,5′-oligoadenylate synthetase 2	13.65	上调 Up
IFIT3	扰素诱导的三角形四肽重复蛋白3 Interferon induced protein with tetratricopeptide repeats 3	12.98	上调 Up
OASL	2′,5′-寡腺苷酸合成酶样蛋白 2′,5′-oligoadenylate synthetase like	12.67	上调 Up
RSAD2	含基S-腺苷蛋氨酸结构域的蛋白2 Radical S-adenosyl methionine domain containing 2	12.63	上调 Up
XAF1	X连锁凋亡抑制蛋白相关因子1 X-linked inhibitor of apoptosis-associated factor 1	12.57	上调 Up
OR51G2	嗅觉受体家族51亚家族G成员2 Olfactory receptor family 51 subfamily G member 2	-127.96	下调Down
RGS3	G蛋白信号转导调节因子3 Regulator of G protein signaling 3	-5.84	下调Down
LIPG	内皮脂肪酶 Lipase G, endothelial type	-5.71	下调Down
OR51F2	嗅觉感受器家族51亚家族F成员2 Olfactory receptor family 51 subfamily F member 2	-4.65	下调Down
DIO3	碘甲腺原氨酸脱碘酶3 Iodothyronine deiodinase 3	-4.63	下调Down
NCRNA00204B	非蛋白编码RNA 204B Non-protein coding RNA 204B	-4.61	下调Down
TGIF2LY	重组人TGF-β诱导的转录因子2样蛋白 TGF-β induced factor homeobox 2 like Y-linked	-4.51	下调Down
C18orf26	18号染色体开放阅读框26 Chromosome 18 open reading frame 26	-4.21	下调Down
ZFYVE9	锌指FYVE域含蛋白9 Zinc finger FYVE-type containing 9	-4.13	下调Down
APCDD1	腺瘤性息肉病大肠埃希菌下调蛋白1 Adenomatosis polyposis coli down-regulated 1	-3.92	下调Down

表1

表2

表3

表4

参考文献 21

[1]	Furtado JM, Smith JR, Belfort R, et al. Toxoplasmosis: a global threat[J]. J Glob Infect Dis, 2011, 3(3):281-284.
[2]	Xu YW, Xing RX, Zhang WH, et al. Toxoplasma ROP16Ⅰ/Ⅲ ameliorated inflammatory bowel diseases via inducing M2 phenotype of macrophages[J]. World J Gastroenterol, 2019, 25(45):6634-6652.
[3]	Cui W, Wang C, Luo QL, et al. Toxoplasma gondii ROP16I deletion: the exacerbated impact on adverse pregnant outcomes in mice[J]. Front Microbiol, 2019, 10:3151.
[4]	Tuladhar S, Kochanowsky JA, Bhaskara A, et al. The ROP16Ⅲ-dependent early immune response determines the sub-acute CNS immune response and type Ⅲ Toxoplasma gondii survival[J]. PLoS Pathog, 2019, 15(10):e1007856.
[5]	Wang C, Cheng WS, Yu Q, et al. Toxoplasma Chinese 1 strain of WH3Δrop16Ⅰ/Ⅲ/gra15Ⅱ genetic background contributes to abnormal pregnant outcomes in murine model[J]. Front Immunol, 2018, 9:1222.
[6]	Wang L, Shen JL. Research progress on genotype and genotype-associated pathogenesis of Toxoplasma gondii[J]. Chin J Parasitol Parasit Dis, 2013, 31(4):319-324. (in Chinese)
[6]	(王林, 沈继龙. 刚地弓形虫基因型和与基因型相关的致病机制研究进展[J]. 中国寄生虫学与寄生虫病杂志, 2013, 31(4):319-324.)
[7]	Liu JF, LV FL. Advances in research on toxoplasmosis of the lung[J]. J Trop Med, 2014, 14(4):541-544. (in Chinese)
[7]	(刘晋锋, 吕芳丽. 肺弓形虫病研究进展[J]. 热带医学杂志, 2014, 14(4):541-544.)
[8]	Su YJ, Dong H, Wang YJ, et al. Differential expression profiling of A549 cells induced by Toxoplasma polymorphism effector ROP16[J]. Chin J Zoonoses, 2018, 34(3):200-206. (in Chinese)
[8]	(苏雅静, 董辉, 王云杰, 等. 弓形虫多态性ROP16效应分子诱导宿主A549细胞基因表达谱的差异研究[J]. 中国人兽共患病学报, 2018, 34(3):200-206.)
[9]	Chen J, Bardes EE, Aronow BJ, et al. ToppGene Suite for gene list enrichment analysis and candidate gene prioritization[J]. Nucleic Acids Res, 2009, 37(web server issue):W305-W311.
[10]	Shuldiner AR, O’Connell JR, Bliden KP, et al. Association of cytochrome P450 2C19 genotype with the antiplatelet effect and clinical efficacy of clopidogrel therapy[J]. JAMA, 2009, 302(8):849-857.
[11]	Su YJ, Dong H, Qiao X, et al. Difference of expression profiling of A549 cells induced by Toxoplasma effector ROP16Ⅱ[J]. Chin J Zoonoses, 2018, 34(4):323-329. (in Chinese)
[11]	(苏雅静, 董辉, 乔霞, 等. 弓形虫ROP16Ⅱ效应分子对宿主A549细胞基因表达谱的影响[J]. 中国人兽共患病学报, 2018, 34(4):323-329.)
[12]	Su YJ, Qiao X, Liu XL, et al. Effects of gene expression profiles of A549 cells induced by Toxoplasma effector ROP16Ⅰ by using bioinformatics analysis[J]. Chin J Zoonoses, 2019, 35(10):957-963, 972. (in Chinese)
[12]	(苏雅静, 乔霞, 刘学雷, 等. 基于生物信息学探讨弓形虫ROP16Ⅰ对A549细胞表达谱的影响[J]. 中国人兽共患病学报, 2019, 35(10):957-963, 972.)
[13]	Eusemann TN, Willmroth F, Fiebich B, et al. Adenosine receptors differentially regulate the expression of regulators of G-protein signalling (RGS) 2, 3 and 4 in astrocyte-like cells[J]. PLoS One, 2015, 10(8):e0134934.
[14]	Melo MB, Nguyen QP, Cordeiro C, et al. Transcriptional analysis of murine macrophages infected with different Toxoplasma strains identifies novel regulation of host signaling pathways[J]. PLoS Pathog, 2013, 9(12):e1003779.
[15]	Lin QW, Zhang S, Lu WQ. Research progress of NOD-like signaling pathways and the relationship between NOD and tumor[J]. China Oncol, 2019, 29(3):223-228. (in Chinese)
[15]	(林巧卫, 张思, 陆维祺. NOD样受体介导的信号转导通路及其与肿瘤关系的研究进展[J]. 中国癌症杂志, 2019, 29(3):223-228.)
[16]	Gao YJ, Liu DL, Li S, et al. Down-regulation of CXCL11 inhibits colorectal cancer cell growth and epithelial-mesenchymal transition[J]. Oncotargets Ther, 2018, 11:7333-7343.
[17]	Beiting DP, Peixoto L, Akopyants NS, et al. Differential induction of TLR3-dependent innate immune signaling by closely related parasite species[J]. PLoS One, 2014, 9(2):e88398.
[18]	Arimoto KI, Löchte S, Stoner SA, et al. STAT2 is an essential adaptor in USP18-mediated suppression of type I interferon signaling[J]. Nat Struct Mol Biol, 2017, 24(3):279-289.
[19]	Kraemer T, Blasczyk R, Bade-Doeding C. HLA-E: a novel player for histocompatibility[J]. J Immunol Res, 2014, 2014:352160.
[20]	Guan D, Factor D, Liu Y, et al. The epigenetic regulator UHRF1 promotes ubiquitination-mediated degradation of the tumor-suppressor protein promyelocytic leukemia protein[J]. Oncogene, 2013, 32(33):3819-3828.
[21]	Ahn HJ, Kim JY, Nam HW. IL-4 independent nuclear translocalization of STAT6 in HeLa cells by entry of Toxoplasma gondii[J]. Korean J Parasitol, 2009, 47(2):117-124.

序号Rank	类别Category	GO号 GO ID	GO注释 GO annotation	基因数Count	P值P value
上调Up-regulated
1	BP	0051607	病毒防御反应 Defense response to virus	46	< 0.01
2	BP	0009615	对病毒的应答反应 Response to virus	52	< 0.01
3	BP	0006952	防御反应 Defense response	94	< 0.01
4	BP	0060337	I型干扰素信号通路 Type I interferon signaling pathway	32	< 0.01
5	BP	0071357	细胞对I型干扰素的反应 Cellular response to type I interferon	32	< 0.01
6	BP	0034340	对I型干扰素的反应 Response to type I interferon	32	< 0.01
7	BP	0045087	先天免疫反应 Innate immune response	73	< 0.01
8	BP	0098542	对其他生物的防御反应 Defense response to other organism	52	< 0.01
9	BP	0043207	对外界生物刺激的反应 Response to external biotic stimulus	65	< 0.01
10	BP	0051707	对其他有机体的反应 Response to other organism	65	< 0.01
下调Down-regulated
1	CC	0005886	细胞质膜 Plasma membrane	40	< 0.01
2	CC	0071944	细胞外围 Cell periphery	40	< 0.01
3	CC	0009986	细胞表面 Cell surface	12	< 0.01
4	CC	0005887	细胞膜的组成成分 Integral component of plasma membrane	16	< 0.01
5	BP	0006814	钠离子运输 Sodium ion transport	6	< 0.01
6	CC	0044459	细胞膜部分 Plasma membrane part	22	< 0.01
7	CC	0031226	细胞膜的固有成分 Intrinsic component of plasma membrane	16	< 0.01
8	CC	0016021	膜组成部分 Integral component of membrane	40	< 0.01
9	BP	0032501	多细胞生物过程 Multicellular organismal process	46	< 0.01
10	MF	0004872	受体活性Receptor activity	17	< 0.01

序号Rank	通路号 Pathway ID	信号通路 Signaling pathway	基因数Gene count	P值 P value
上调 Up-regulated
1	Hsa 05168	单纯疱疹病毒感染-智人（人类） Herpes simplex infection-Homo sapiens (human)	26	< 0.01
2	Hsa 05164	甲型流感-智人（人类） Influenza A-Homo sapiens (human)	21	< 0.01
3	Hsa 05162	麻疹-智人（人类） Measles-Homo sapiens (human)	15	< 0.01
4	Hsa 04612	抗原处理及呈递-智人（人类） Antigen processing and presentation-Homo sapiens (human)	11	< 0.01
5	Hsa 05160	丙型肝炎-智人（人类） Hepatitis C-Homo sapiens (human)	13	< 0.01
6	Hsa 04623	细胞质DNA传感通路-智人（人类） Cytosolic DNA-sensing pathway-Homo sapiens (human)	9	< 0.01
7	Hsa 04620	Toll样受体信号通路-智人（人类） Toll-like receptor signaling pathway-Homo sapiens (human)	11	< 0.01
8	Hsa 05332	移植物抗宿主病-智人（人类） Graft-versus-host disease-Homo sapiens (human)	7	< 0.01
9	Hsa 05330	同种异体移植排斥反应-智人（人类） Allograft rejection-Homo sapiens (human)	6	< 0.01
10	Hsa 04621	NOD样受体信号通路-智人（人类） NOD-like receptor signaling pathway-Homo sapiens (human)	7	< 0.01
下调 Down-regulated
1	Hsa 00020	柠檬酸循环(TCA循环)-智人（人类） Citrate cycle (TCA cycle)-Homo sapiens (human)	2	< 0.05
2	Hsa 04024	cAMP信号通路-智人（人） cAMP signaling pathway-Homo sapiens (human)	4	< 0.05
3	Hsa 04360	轴突导向-智人（人类）Axon guidance-Homo sapiens (human)	3	< 0.05

序号 Rank	ToppGene
序号 Rank	基因标志 Gene symbol	平均得分 Average score（× 10^-1）	P值/× 10^-5 P value/× 10^-5	基因标志 Gene symbol	相互作用基因数 Interactant count	得分/× 10^-4 Score/× 10^-4
1	IL6	9.84	0.16	CCL5	41	12.90
2	MYH11	7.30	2.43	CCL2	21	10.10
3	STAT1	8.79	2.75	STAT1	199	9.72
4	CXCR4	9.71	3.05	CCL4	11	8.90
5	B2M	8.25	3.60	CCL3	13	8.52
6	CCL3	9.13	3.91	CXCL11	5	6.76
7	CCL2	9.09	4.47	CCL26	3	6.68
8	HLA-A	8.76	4.97	BIRC3	102	5.04
9	CCL5	9.02	7.59	AMOT	104	4.05
10	CXCL11	8.98	7.62	HDAC5	355	3.94
11	NOD2	8.58	8.10	NMI	70	3.72
12	HLA-B	8.66	8.23	PML	257	3.61
13	TLR3	8.31	11.10	IFI16	358	3.48
14	HLA-C	8.32	12.20	IRF7	48	3.39
15	HLA-G	8.61	16.20	TNFAIP3	99	3.35
16	TNFSF10	8.90	16.50	STAT2	51	3.30
17	CCL4	8.96	16.90	CXCR4	53	3.21
18	FGA	7.94	17.70	TLR3	32	3.14
19	LTB	8.81	18.90	LMO2	208	3.02
20	IFNB1	8.90	20.10	PARP14	14	2.91
21	TNFAIP3	7.77	22.10	CISH	34	2.86
22	FGFR2	7.82	24.80	CD93	35	2.63
23	CCL26	8.88	28.90	OPTN	80	2.45
24	ERAP1	7.78	30.40	BST2	13	2.39
25	DDX58	7.98	36.90	TRIM21	99	2.20
26	HLA-E	8.44	37.70	HLA-E	54	2.14
27	IFNA16	8.68	48.60	ISG15	189	1.88
28	RET	7.13	56.90	PCDH7	40	1.82
29	PML	7.53	58.30	EIF2AK2	187	1.68
30	STAT2	6.56	62.70	IL6	13	1.65

弓形虫毒力相关效应分子ROP16Ⅲ诱导人肺腺癌A549细胞的基因表达谱

Gene expression profiling of human lung adenocarcinoma of A549 cells induced by Toxoplasma virulence-related effector ROP16Ⅲ

RichHTML

PDF (PC)

可视化

摘要/Abstract

引用本文

使用本文

图/表 4

参考文献 21

相关文章 15

编辑推荐

Metrics

本文评价

[1]	李启松, 杨李, 腾飞翔, 马桂芳. 屋尘螨过敏原Der p 4的克隆、表达、免疫原性鉴定及生物信息学分析[J]. 中国寄生虫学与寄生虫病杂志, 2024, 42(1): 42-47.
[2]	仲顺虎, 孙玥, 郭小腊, 郑亚东, 陈轶霞. 多房棘球蚴感染小鼠脾淋巴细胞中差异表达miRNA的鉴定及其生物信息学分析[J]. 中国寄生虫学与寄生虫病杂志, 2022, 40(3): 288-294.
[3]	何成, 潘帅, 许妹珍, 原飞, 何静妹, 刘转转. 刚地弓形虫蛋白磷酸酶的全蛋白质组鉴定及生物信息学分析[J]. 中国寄生虫学与寄生虫病杂志, 2021, 39(1): 120-124.
[4]	范俊杰, 韩秀敏, Nur Fazleen Binti Idris, 李锴, 谭青青, 曹纹侨, 李想, 廖鹏, 叶彬. 细粒棘球绦虫蛋白激酶A的生物信息学特征及免疫反应性[J]. 中国寄生虫学与寄生虫病杂志, 2020, 38(6): 682-687.
[5]	郭楚桐, 欧阳春艳, 何俊贤, 季树宇, 杨礼腾, 刘晓宇. 屋尘螨第十类变应原的克隆、表达、纯化及免疫原性分析[J]. 中国寄生虫学与寄生虫病杂志, 2020, 38(6): 723-729.
[6]	郑文琪, 冯学敏, 曹雅明, 韩艳秋, 王俊瑞. 伯氏疟原虫静息巯基氧化酶的生物信息学分析和体外酶活检测[J]. 中国寄生虫学与寄生虫病杂志, 2020, 38(2): 159-165.
[7]	古霞, 欧阳春艳, 袁如意, 袁谢芳, 李国平, 刘志刚. 屋尘螨3类变应原基因的克隆、表达及免疫原性分析[J]. 中国寄生虫学与寄生虫病杂志, 2018, 36(6): 627-631.
[8]	李想, 罗波, 周必英. 重要绦虫功能基因组学研究进展[J]. 中国寄生虫学与寄生虫病杂志, 2018, 36(5): 520-525.
[9]	姚嘉青，刘丛珊，薛剑，陶奕，张皓冰*. 细粒棘球蚴原头节β4微管蛋白的可溶性表达、纯化及生物信息学分析[J]. 中国寄生虫学与寄生虫病杂志, 2018, 36(3): 6-231-237.
[10]	蔡玉春，陈韶红，李浩，卢艳，艾琳，储言红，沈慧敏，陈家旭*. 田鼠巴贝虫可溶性虫体蛋白的质谱鉴定及生物信息学分析[J]. 中国寄生虫学与寄生虫病杂志, 2018, 36(3): 8-246-252.
[11]	程金芝, 刘鉴, 翟慧, 吕清巧, 颜凤, 商正玲, 吴家红. 埃及伊蚊丝氨酸蛋白酶抑制蛋白基因家族的筛选和表达谱分析[J]. 中国寄生虫学与寄生虫病杂志, 2017, 35(6): 527-535.
[12]	王楠1, 滕飞翔1, 俞黎黎1, 杨李1, 张承伯1, 周鹰1,2, 崔玉宝1,2*. 粉尘螨变应原Der f 8全长基因克隆及表达[J]. 中国寄生虫学与寄生虫病杂志, 2017, 35(4): 8-351-356.
[13]	魏超君1，贾彦娟1，徐辉1，吴玲1，卢思奇2，李永红1*. 蓝氏贾第鞭毛虫α-8贾第素特异性优势抗原表位肽的抗原性分析[J]. 中国寄生虫学与寄生虫病杂志, 2016, 34(5): 1-389-393.
[14]	颜凤，吕清巧，程金芝，牟小会，翟慧，杨迅，商正玲，吴家红. 埃及伊蚊yellow*基因家族的鉴别和表达谱分析[J]. 中国寄生虫学与寄生虫病杂志, 2016, 34(1): 1-1-10.
[15]	史世俊，崔勇，李瑾，王洪法，尹昆，魏庆宽，黄炳成，孙慧，刘功振*. 弓形虫TgROP21基因克隆及生物信息学分析[J]. 中国寄生虫学与寄生虫病杂志, 2015, 33(4): 16-294-296.