Blastocystis sp. infection in dairy cattle in Yunnan Province and its gene subtype analysis

doi:10.12140/j.issn.1000-7423.2021.06.012

Abstract

Abstract:

Objective To understand the prevalence of Blastocystis sp. in dairy cattle and its gene subtype in partial areas of Yunnan Province, thereby to provide information support for the prevention and control of the protozoan infection. Methods Fecal samples were collected from Holstein cows and dairy buffalos in Kunming, Dali, and Baoshan City from June 2019 to August 2020. DNA was extracted from the fecal samples, and the small subunit (SSU) rRNA of Blastocystis sp. was amplified by nested PCR and sequenced. The sequence alignment was performed with BLAST on NCBI database. Phylogenetic trees was constructed using neighbor-joining method of MEGA 6.0. SPSS 20.0 and SAS 9.1 software were used to analyze the infection rates of dairy cattle in different regions, seasons, farming patterns, breeds, genders and growth stages, and those factors with significant difference were further analysised by binary logistic regression. Results A total of 841 fecal samples were collected, of which 65 were SSU rRNA positive with the amplified fragment being 1 100 bp in length. The overall infection rate of Blastocystis sp. in dairy cattle was 7.7% (65/841, 95%CI: 5.92%-9.53%). The infection rates were 12.9% (32/248, 95%CI: 8.73%-17.08%), 8.1% (19/236, 95%CI: 4.58%-11.52%) and 3.9% (14/357, 95%CI: 1.91%~5.94%) in Kunming, Baoshan and Dali, with significant differences between places(χ² = 16.601, P < 0.01); the infection rate in summer, autumn and winter was 9.4% (56/599, 95%CI: 7.86%-37.30%), 8.3% (8/96, 95%CI: 2.80%-13.86%) and 0.7% (1/146, 95%CI: 0-2.02%), respectively, showing significant differences seasonally (χ² = 12.411, P < 0.05); the infection rate in captive farming and field farming was 8.3% (58/702, 95%CI: 6.23%-10.30%) and 5.0% (7/139, 95%CI: 1.40%-8.67%), respectively, with no significant difference (χ² = 1.693, P > 0.05); the infection rate among Holstein cows and dairy buffalos was 8.4% (41/490, 95%CI: 5.92%-10.82%) and 6.8% (24/351, 95%CI: 4.20%-9.48%), respectively, showing no significant difference (χ² = 0.671, P > 0.05); the rate in bull and cow was 11.1% (10/90, 95%CI: 4.62%-17.60%) and 7.3% (55/751, 95%CI: 5.46%-9.19%), repectively, presenting no significant difference (χ² = 1.617, P > 0.05); the rate in different growth stages of pre-weaned calves, post-weaned calves, juvenile cattle and adult cattle was 2/18 (95%CI: 0-25.63%), 26.2% (11/42, 95%CI: 12.89%-39.49%), 15.1% (11/73, 95%CI: 6.86%-23.28%) and 5.8% (41/708, 95%CI: 4.07%-7.51%), respectively, showing significant differences between growing stages(χ² = 29.604, P < 0.01). Binary logistic regression analysis showed that season (P < 0.05) and growth stage (P < 0.01) were independent influencing factor for Blastocystis sp. infection in dairy cattle. The sequence alignment found that the SSU rRNA sequence from one fecal sample was 96.3% homologous to reference sequence GenBank KX618192 (ST3 subtype), thirty-five samples were 97.3%-99.7% homologous to the reference sequences MH104059, KC148207 and MH358363 (ST10 subtype), twenty-three samples were 98.7%-100% homologous to the reference sequence MT453998 and MT453999 (ST12 subtype) and six samples were 98.7%-100% homologous to the reference sequence MH358361, MH358362, KC148205 and KC148206 (ST14 subtype). Phylogenetic tree analysis showed that the SSU rRNA sequences of Blastocystis sp. were clustered at to ST3, ST10, ST12 and ST14 subtype clade, respectively, having the constituent ratio of the subtypes being 1.5% (1/65), 53.9% (35/65), 35.4% (23/65) and 9.2% (6/65), respectively. Conclusion Blastocystis sp. infection are found in dairy cattle in Yunnan Province, and its ST3 and ST12 subtypes may have zoonotic transmission risk, thus, it is imperative to strengthen the control and prevention of the pathogen in the province.

Key words: Blastocystis sp., Zoonosis, Dairy cattle, Yunnan Province

CLC Number:

R531.19

PU Li-hua, ZHANG Xing-ze, GUAN Shao-jun, CHENG Wen-jie, ZOU Feng-cai, MAO Hua-ming, YANG Jian-fa. Blastocystis sp. infection in dairy cattle in Yunnan Province and its gene subtype analysis[J]. CHINESE JOURNAL OF PARASITOLOGY AND PARASITIC DISEASES, 2021, 39(6): 809-815.

Figures/Tables 4

Fig. 1

Table 1

Table 2

Fig. 2

References 29

[1]	Liu LK, Wang PL, Han H, et al. Current status of Blastocystis infection in China[J]. Chin J Zoonoses, 2021, 37(6): 548-555, 562. (in Chinese)
	(刘林科, 王朋林, 韩涵, 等. 我国芽囊原虫感染现状[J]. 中国人兽共患病学报, 2021, 37(6): 548-555, 562.)
[2]	Zhang X, Qiao JY, Dong XH, et al. Study on morphology of Blastocystis hominis in culture and from diarrhea patients[J]. Chin J Parasitol Parasit Dis, 2003, 21(2): 116-118. (in Chinese)
	(张旭, 乔继英, 董小惠, 等. 腹泻患者及培养基中人芽囊原虫的各期形态研究[J]. 中国寄生虫学与寄生虫病杂志, 2003, 21(2): 116-118.)
[3]	Tian CL, Liu DY, Lu ZC, et al. Stool examination and morphological observation of Blastacystis hominis[J]. Chin Trop Med, 2005, 5(2): 221-223. (in Chinese)
	(田春林, 刘登宇, 卢作超, 等. 人芽囊原虫粪便检查与形态观察[J]. 中国热带医学, 2005, 5(2): 221-223.)
[4]	Li TC, Pu LH, Li HX. Advances in studies on non-human primates infected with Blastocystis[J]. Chin J Wild, 2019, 40(4): 1050-1054. (in Chinese)
	(李廷翠, 普丽花, 李红霞. 非人灵长类动物感染芽囊原虫的研究进展[J]. 野生动物学报, 2019, 40(4): 1050-1054.)
[5]	Qian K, Hu ZH, Wu XP, et al. Prevalence and genotype analysis of Blastocystis hominis among HIV/AIDS patients in Nanchang City[J]. Chin J Parasitol Parasit Dis, 2021, 39(1): 55-60. (in Chinese)
	(钱科, 胡主花, 吴秀萍, 等. 南昌市HIV/AIDS患者人芽囊原虫感染现况及基因型分析[J]. 中国寄生虫学与寄生虫病杂志, 2021, 39(1): 55-60.)
[6]	Maloney JG, Lombard JE, Urie NJ, et al. Zoonotic and genetically diverse subtypes of Blastocystis in US pre-weaned dairy heifer calves[J]. Parasitol Res, 2019, 118(2): 575-582. doi: 10.1007/s00436-018-6149-3 pmid: 30483890
[7]	Ramírez JD, Sánchez A, Hernández C, et al. Geographic distribution of human Blastocystis subtypes in South America[J]. Infect Genet Evol, 2016, 41: 32-35. doi: S1567-1348(16)30086-7 pmid: 27034056
[8]	Li LH, Zhang XP, Lv S, et al. Cross-sectional surveys and subtype classification of human Blastocystis isolates from four epidemiological settings in China[J]. Parasitol Res, 2007, 102(1): 83-90. doi: 10.1007/s00436-007-0727-0
[9]	Tao W. Molecular characterization and genetic features of three enteric protozoa with zoonotic importance in dairy cattle in partial areas of Heilongjiang Province[D]. Harbin: Northeast Agricultural University, 2016: 1-66. (in Chinese)
	(陶威. 黑龙江省部分地区奶牛三种重要人兽共患原虫的分子鉴定及其遗传学特征[D]. 哈尔滨: 东北农业大学, 2016: 1-66.)
[10]	Zhu W, Tao W, Gong B, et al. First report of Blastocystis infections in cattle in China[J]. Vet Parasitol, 2017, 246: 38-42. doi: 10.1016/j.vetpar.2017.09.001
[11]	Wang J, Gong B, Yang F, et al. Subtype distribution and genetic characterizations of Blastocystis in pigs, cattle, sheep and goats in northeastern China’s Heilongjiang Province[J]. Infect Genet Evol, 2018, 57: 171-176. doi: 10.1016/j.meegid.2017.11.026
[12]	Wang SS. Study on prevalence and subtypes of Blastocystis sp. from calves in partial areas of Shaanxi Province[D]. Yangling: Northwest A & F University, 2018: 1-40. (in Chinese)
	(王莎莎. 陕西部分地区犊牛芽囊原虫的感染情况及亚型研究[D]. 杨凌: 西北农林科技大学, 2018: 1-40.)
[13]	Quílez J, Sánchez-Acedo C, Clavel A, et al. Occurrence of Blastocystis sp. in cattle in Aragón, northeastern Spain[J]. Parasitol Res, 1995, 81(8): 703-705. pmid: 8570589
[14]	Masuda A, Sumiyoshi T, Ohtaki T, et al. Prevalence and molecular subtyping of Blastocystis from dairy cattle in Kanagawa, Japan[J]. Parasitol Int, 2018, 67(6): 702-705. doi: S1383-5769(18)30098-9 pmid: 30009956
[15]	Fayer R, Santin M, Macarisin D. Detection of concurrent infection of dairy cattle with Blastocystis, Cryptosporidium, Giardia, and Enterocytozoon by molecular and microscopic methods[J]. Parasitol Res, 2012, 111(3): 1349-1355. doi: 10.1007/s00436-012-2971-1
[16]	Ramírez JD, Sánchez LV, Bautista DC, et al. Blastocystis subtypes detected in humans and animals from Colombia[J]. Infect Genet Evol, 2014, 22: 223-228. doi: 10.1016/j.meegid.2013.07.020 pmid: 23886615
[17]	Badparva E, Sadraee J, Kheirandish F. Genetic diversity of Blastocystis isolated from cattle in Khorramabad, Iran[J]. Jundishapur J Microbiol, 2015, 8(3): e14810.
[18]	Sharifi Y, Abbasi F, Shahabi S, et al. Comparative genotyping of Blastocystis infecting cattle and human in the south of Iran[J]. Comp Immunol Microbiol Infect Dis, 2020, 72: 101529. doi: 10.1016/j.cimid.2020.101529
[19]	Udonsom R, Prasertbun R, Mahittikorn A, et al. Blastocystis infection and subtype distribution in humans, cattle, goats, and pigs in central and western Thailand[J]. Infect Genet Evol, 2018, 65: 107-111. doi: 10.1016/j.meegid.2018.07.007
[20]	Aynur ZE, Güçlü Ö, Yıldız İ, et al. Molecular characterization of Blastocystis in cattle in Turkey[J]. Parasitol Res, 2019, 118(3): 1055-1059. doi: 10.1007/s00436-019-06243-8
[21]	Greige S, El Safadi D, Khaled S, et al. First report on the prevalence and subtype distribution of Blastocystis sp. in dairy cattle in Lebanon and assessment of zoonotic transmission[J]. Acta Trop, 2019, 194: 23-29. doi: 10.1016/j.actatropica.2019.02.013
[22]	Suwanti LT, Susana Y, Hastutiek P, et al. Blastocystis spp. subtype 10 infected beef cattle in Kamal and Socah, Bangkalan, Madura, Indonesia[J]. Vet World, 2020, 13(2): 231-237. doi: 10.14202/vetworld.
[23]	Wang LY, Cao LT, Wang YG, et al. Infection and subtype analysis of Blastocystis in some dairy farms in Guangdong[J]. Chin J Vet Sci, 2020, 40(2): 345-348. (in Chinese)
	(王璐阳, 曹乐天, 王艳歌, 等. 广东省部分奶牛场芽囊原虫感染情况及亚型分析[J]. 中国兽医学报, 2020, 40(2): 345-348.)
[24]	Fu Y, Zhao HY, Chen YC, et al. Molecular epidemiology of Enterocytozoon bieneusi and Blastocystis sp. in dairy cattle in Henan Province[J]. Chin J Zoonoses, 2021, 37(10): 960-964. (in Chinese)
	(付寅, 赵洪育, 陈远才, 等. 河南省部分地区奶牛毕氏肠微孢子虫和芽囊原虫分子流行病学调查[J]. 中国人兽共患病学报, 2021, 37(10): 960-964.)
[25]	Ren M, Song JK, Yang F, et al. First genotyping of Blastocystis in yaks from Qinghai Province, northwestern China[J]. Parasit Vectors, 2019, 12(1): 171. doi: 10.1186/s13071-019-3436-5
[26]	Liang XX, Zou Y, Li TS, et al. Molecular epidemiological studies of Enterocytozoon bieneusi and Blastocystis sp. in Yunling cattle in Yunnan Province[J]. Chin J Animal Infect Dis, 2021, 29(5): 75-82. (in Chinese)
	(梁霞霞, 邹扬, 李涛善, 等. 云岭牛毕氏肠微孢子虫和芽囊原虫分子流行病学调查研究[J]. 中国动物传染病学报, 2021, 29(5): 75-82.)
[27]	Teng XJ, Chu YH, Zhai CC, et al. The epidemiological characteristics and influencing factors for Blastocystis hominis infection among human immunodeficiency virus seropositive individuals in Tengchong of Yunnan Province[J]. Chin J Parasitol Parasit Dis, 2018, 36(2): 129-134. (in Chinese)
	(滕雪娇, 储言红, 翟铖铖, 等. 云南腾冲人类免疫缺陷病毒携带者人芽囊原虫感染的流行病学特征及影响因素分析[J]. 中国寄生虫学与寄生虫病杂志, 2018, 36(2): 129-134.)
[28]	Liu Q. Studies on prevalence and genotypes of G. lamblia, B. hominis and E. bieneusi in humans in Dali[D]. Dali: Dali University, 2018: 1-67. (in Chinese)
	(刘晴. 大理地区人群贾第虫、芽囊原虫和微孢子虫分子流行病学调查研究[D]. 大理: 大理大学, 2018: 1-67.)
[29]	Lee H, Lee SH, Seo MG, et al. Occurrence and genetic diversity of Blastocystis in Korean cattle[J]. Vet Parasitol, 2018, 258: 70-73. doi: 10.1016/j.vetpar.2018.06.010

类别 Category	样本量 No. tested	阳性数 No. positive	感染率/% Infection rate/%	95%置信区间/% 95%CI/%	χ²	P值 P-vlaule	基因亚型（样品数） Subtype （No. sample）
地区Region
昆明Kunming	248	32	12.9	8.73~17.08	16.601	< 0.01	ST10（23）、ST12（6）、ST14（3）
大理Dali	357	14	3.9	1.91~5.94			ST10（6）、ST12（6）、ST14（2）
保山Baoshan	236	19	8.1	4.58~11.52			ST3（1）、ST10（6）、ST12（11）、ST14（1）
季节Season
夏季Summer	599	56	9.4	7.86~37.30	12.411	< 0.05	ST3（1）、ST10（31）、ST12（20）、ST14（4）
秋季Autumn	96	8	8.3	2.80~13.86			ST10（3）、ST12（3）、ST14（2）
冬季Winter	146	1	0.7	0~2.02			ST10（1）
品种Breed
荷斯坦奶牛Holstein cow	490	41	8.4	5.92~10.82	0.671	> 0.05	ST10（27）、ST12（9）、ST14（5）
奶水牛Dairy buffalo	351	24	6.8	4.20~9.48			ST3（1）、ST10（8）、ST12（14）、ST14（1）
养殖模式Farming pattern
圈养Captive	702	58	8.3	6.23~10.30	1.693	> 0.05	ST3（1）、ST10（31）、ST12（20）、ST14（6）
放牧Graze	139	7	5.0	1.40~8.67			ST10（4）、ST12（3）
性别Gender
公牛 Bull	90	10	11.1	4.62~17.60	1.617	> 0.05	ST3（1）、ST10（31）、ST12（20）、ST14（6）
母牛 Cow	751	55	7.3	5.46~9.19			ST10（29）、ST12（21）、ST14（5）
发育阶段 Growth stage
断奶前Pre-weaned	18	2	11.1	0~25.63	29.604	< 0.01	ST3（1）、ST10（1）
断奶后Post-weaned	42	11	26.2	12.89~39.49			ST10（10）、ST14（1）
育成牛Juvenile	73	11	15.1	6.86~23.28			ST10（7）、ST12（3）、ST14（1）
成年牛Adult	708	41	5.8	4.07~7.51			ST10（17）、ST12（20）、ST14（4）
总计Total	841	65	7.7	5.92~9.53			ST3（1）、ST10（35）、ST12（23）、ST14（6）

因素 Factor	回归系数 Regression coefficien	标准误差 S.E.	瓦尔德 Wald	自由度 df	P值 P-vlaule	比值比 OR	95%CI
因素 Factor	回归系数 Regression coefficien	标准误差 S.E.	瓦尔德 Wald	自由度 df	P值 P-vlaule	比值比 OR	上限Lower	下限Upper
地区Region	-0.036	0.190	0.036	1	0.849	0.964	0.664	1.400
季节Season	-0.722	0.304	5.650	1	0.017	0.486	0.268	0.881
发育阶段 Growth stage	-0.496	0.147	11.314	1	0.001	0.609	0.456	0.813
常量Constant	0.319	0.826	0.149	1	0.700	1.375