馬海璐，劉永宏，趙　麗，曹勝波，劉學(xué)鋒，王海國，李岳洋，馬世強(qiáng)，蔣秀梅，張保軍，薩依甫加瑪·卡依薩爾，努爾古麗·圖爾貢，廖秋萍，焦海宏

(1.塔里木大學(xué) 動(dòng)物科學(xué)學(xué)院，新疆生產(chǎn)建設(shè)兵團(tuán)塔里木畜牧科技重點(diǎn)實(shí)驗(yàn)室，新疆阿拉爾　843300；2.華中農(nóng)業(yè)大學(xué) 動(dòng)物醫(yī)學(xué)院，武漢　430070；3.新疆生產(chǎn)兵團(tuán) 第十四師一牧場，新疆和田　848306；4.新疆生產(chǎn)建設(shè)兵團(tuán) 第三師紅旗農(nóng)場生產(chǎn)科，新疆阿圖什　845350)

?

中國新疆小反芻獸疫病毒F基因序列分析

馬海璐1，劉永宏1，趙麗1，曹勝波2，劉學(xué)鋒3，王海國4，李岳洋1，馬世強(qiáng)1，蔣秀梅1，張保軍1，薩依甫加瑪·卡依薩爾1，努爾古麗·圖爾貢1，廖秋萍1，焦海宏1

(1.塔里木大學(xué) 動(dòng)物科學(xué)學(xué)院，新疆生產(chǎn)建設(shè)兵團(tuán)塔里木畜牧科技重點(diǎn)實(shí)驗(yàn)室，新疆阿拉爾843300；2.華中農(nóng)業(yè)大學(xué) 動(dòng)物醫(yī)學(xué)院，武漢430070；3.新疆生產(chǎn)兵團(tuán) 第十四師一牧場，新疆和田848306；4.新疆生產(chǎn)建設(shè)兵團(tuán) 第三師紅旗農(nóng)場生產(chǎn)科，新疆阿圖什845350)

摘要為分析中國新疆小反芻獸疫病毒(PPRV)毒株分子演變特點(diǎn)和疫情傳播路線，從GenBank數(shù)據(jù)庫下載所有國家全部的PPRV F基因全序列，應(yīng)用DNAStar軟件，結(jié)合疫情毒株時(shí)空分布進(jìn)行序列分析。結(jié)果顯示，中國新疆PPRV株F基因全序列，與2014年小反芻獸疫情中河南和北京毒株核苷酸的同源性最高，為99.8%～99.9%，與中國西藏3個(gè)毒株核苷酸的同源性為97.7%，與中國采用的疫苗株核苷酸同源性為92.9%，核苷酸同源性最低的是肯尼亞毒株，與國外毒株核苷酸同源性最高的是印度毒株。F基因系統(tǒng)進(jìn)化關(guān)系分析顯示，中國新疆PPRV株屬于基因Ⅳ系。F蛋白氨基酸同源性結(jié)果顯示，除與科特迪瓦1989年毒株ICV89氨基酸同源性最低外，與其他氨基酸同源性比對結(jié)果均與核苷酸的比對結(jié)果一致。中國新疆PPRV株存在一定程度的變異，可能是由周邊國家傳入。

關(guān)鍵詞中國；新疆；小反芻獸疫；F基因

小反芻獸疫(Peste des petits ruminants，PPR)，是由小反芻獸疫病毒(Peste des petits ruminants virus，PPRV)引起的一種小反芻動(dòng)物急性接觸性傳染疾病。有報(bào)道[1-2]顯示，牛、水牛、駱駝和獅均可感染該病毒，表現(xiàn)為發(fā)熱、眼鼻分泌物增多、口炎、腸炎和肺炎等。

1942年，西非科特迪瓦首次報(bào)道PPR[3]，之后PPR在多個(gè)國家和地區(qū)均有發(fā)生，呈全球蔓延趨勢。據(jù)報(bào)道，目前 至少有49個(gè)國家發(fā)生PPR疫情[2,4]，其中亞洲有23個(gè)國家，中國周邊有12個(gè)國家[4-5]。

2005年，通過發(fā)病情況和血清學(xué)調(diào)查發(fā)現(xiàn)，中國西藏發(fā)生小反芻獸疫疫情；隨后，通過分子流行病學(xué)調(diào)查確證西藏于2007年7月發(fā)生PPR[6-8]。2008年和2010年，西藏至少發(fā)生3次PPR疫情[7-9]。2013年年底，新疆伊犁地區(qū)成為中國第2個(gè)發(fā)生PPR疫情的省區(qū)，隨后，新疆多地暴發(fā)該疫情[10]。2014年，PPR疫情呈大流行趨勢，在中國的至少22省(市、自治區(qū))256個(gè)縣發(fā)生，其中，云南、江蘇、安徽、湖北和貴州等省為重災(zāi)區(qū)，并于4月達(dá)到發(fā)病高峰。截止2014年9月，中國約有7萬多只綿羊和山羊因患PPR被撲殺，損失巨大[11]。

PPRV為副粘病毒科麻疹病毒屬的成員[12]，只有1個(gè)血清型[13]，該病毒核酸類型為單股負(fù)鏈RNA，基因組約為15 948 nt，從3′到5′依次為N、P、M、F、H和L6個(gè)基因，分別編碼核衣殼蛋白N、磷蛋白P、基質(zhì)蛋白M、融合蛋白F、血凝蛋白H、大蛋白L等6個(gè)結(jié)構(gòu)蛋白[14]。

2013年11月，新疆伊犁哈薩克自治州霍城縣三宮鄉(xiāng)某村發(fā)生PPR疫情，提取該村患病山羊淋巴結(jié)組織樣品，通過病毒RNA提取、RT-PCR、PCR產(chǎn)物純化測序及末端RACE法，獲得完整PPRV基因組序列，命名為PPRV China/XJYL/2013，為新奇變異株，GenBank登錄號為KM091959[15]。根據(jù)文獻(xiàn)[1,14]報(bào)道，以N或F基因序列可進(jìn)行遺傳演化分析。

本研究從GenBank數(shù)據(jù)庫下載各國PPRVF基因全序列，并結(jié)合每次疫情發(fā)生的時(shí)空分布，對中國新疆小反芻獸疫疫情毒株分子演變特點(diǎn)和疫情傳播路線進(jìn)行系統(tǒng)分析，以期為該病防治及疫病消滅提供依據(jù)。

1材料與方法

1.1材 料

中國新疆僅有的一個(gè)PPRV China/XJYL/2013株F基因全序列，中國其他地區(qū)和其他國家的全部PPRVF基因全序列，均來源于GenBank數(shù)據(jù)庫，各毒株的具體信息見表1。

表1　參考毒株

1.2方 法

搜索GenBank數(shù)據(jù)庫全部的PPRVF基因全序列，運(yùn)用DNAStar軟件的EditSeq程序?qū)⑵浜塑账岷桶被嵝蛄斜４鏋?.Seq和*.Pro格式文件，運(yùn)用該軟件MegAlign程序的Clustal W法進(jìn)行比對，然后在Veiw菜單中通過Sequence Distance和Phylogenetic Tree提取序列比對結(jié)果，并繪制進(jìn)化樹圖。同時(shí)，結(jié)合毒株時(shí)空分布，分析中國新疆伊犁地區(qū)PPRV來源及演變特征。

2結(jié)果與分析

2.1中國新疆PPRV F基因核苷酸序列分析

對不同時(shí)空分布的32個(gè)PPRVF基因全序列進(jìn)行核苷酸同源性分析，結(jié)果顯示，2013年新疆伊犁地區(qū)PPRV株China/XJYL/2013F基因與2014年P(guān)PP疫情大流行中河南和北京PPRV毒株F基因核苷酸同源性最高，達(dá)到99.8%～99.9%；與中國西藏阿里地區(qū)和那曲地區(qū)3個(gè)PPRV毒株F基因核苷酸同源性均達(dá)到97.7%；與中國疫苗株Nigeria/75/1F基因核苷酸同源性為92.9%；核苷酸同源性最低的是肯尼亞2011年KN5/2011株，與其他國家毒株核苷酸同源性最高的是印度1994年Izatnagar-94株；中國毒株與其他國家毒株核苷酸同源性為87.8%～98.3%(表2右上)。

采用PPRVF基因核苷酸全序列繪制基因進(jìn)化樹(圖1)，結(jié)果顯示，1989年西非科特迪瓦的1個(gè)毒株和1969年塞內(nèi)加爾的1個(gè)毒株構(gòu)成1個(gè)小分支(Ⅱ)；2個(gè)西非尼日利亞的毒株、1個(gè)塞內(nèi)加爾2013年毒株和1個(gè)加納2010年毒株構(gòu)成1個(gè)小分支(Ⅰ)；亞洲鄰國阿聯(lián)酋和阿曼各1個(gè)毒株、非洲東北部埃塞俄比亞1994年1個(gè)毒株、非洲東部的肯尼亞和烏干達(dá)各1個(gè)毒株構(gòu)成為1個(gè)分支(Ⅲ)；亞洲18個(gè)毒株、非洲東北部埃塞俄比亞2010年1個(gè)毒株和非洲西北部摩洛哥2個(gè)毒株構(gòu)成1個(gè)分支(Ⅳ)。F基因系統(tǒng)進(jìn)化樹顯示，埃塞俄比亞和塞內(nèi)加爾不同年代的毒株被劃分到不同的分支；2013年中國新疆伊犁地區(qū)毒株與中國西藏毒株和2014年P(guān)PP疫情中國其他毒株被劃分在同一個(gè)大分支(Ⅳ)，但2013-2014年毒株與中國西藏毒株被劃分為2個(gè)不同的小分支，與孟加拉國2010年毒株系統(tǒng)進(jìn)化關(guān)系最近。

圖1　不同地區(qū)和不同年代的小反芻獸疫病毒毒株F基因系統(tǒng)進(jìn)化樹

2.2中國新疆PPRV F蛋白氨基酸序列分析

對不同時(shí)空分布的32個(gè)PPRV F蛋白全序列進(jìn)行氨基酸同源性分析，結(jié)果顯示，2013年中國新疆伊犁地區(qū)PPRV株China/XJYL/2013 F蛋白與2014年P(guān)PP疫情大流行中國河南和北京PPRV毒株F蛋白氨基酸同源性最高，達(dá)到99.8%～100%；與中國西藏阿里地區(qū)和那曲地區(qū)3個(gè)PPRV毒株F蛋白氨基酸同源性為98.7%～98.9%；與中國疫苗株Nigeria/75/1 F蛋白氨基酸同源性為96.7%；與科特迪瓦1989年毒株ICV89氨基酸同源性最低，與國外毒株氨基酸同源性最高為印度1994年Izatnagar-94株、1996年Sungri/96株和2003年Jhansi/2003株，同源性均為99.1%；中國毒株與其他國家毒株氨基酸同源性為94.0%～99.3%(表2左下)。

3討 論

3.12014年中國PPR疫情來源分析

通過對F基因核苷酸序列進(jìn)行遺傳演化分析[1]，可反映出PPRV更多的以時(shí)間為軸地演化過程，從而為該病毒在不同地理區(qū)域間的傳播途徑提供線索[12]。研究者們用PPRVF基因部分核苷酸序列進(jìn)行分析，將不同來源的PPRV毒株劃分為4系，其中，Ⅰ系由塞內(nèi)加爾和尼日利亞的毒株組成，Ⅱ系由幾內(nèi)亞和科特迪瓦的毒株組成，Ⅲ系包括埃塞俄比亞、也門、蘇丹以及1株1992年分離自印度南部的毒株，Ⅳ系均為亞洲毒株[16]。

本研究采用PPRVF基因全序列進(jìn)行系統(tǒng)進(jìn)化分析，將GenBank數(shù)據(jù)庫的32個(gè)PPRV流行株F基因全序列進(jìn)行分系，結(jié)果顯示，新疆China/XJYL/2013株與西藏3個(gè)毒株均劃分為Ⅳ系，與前人報(bào)道相符[12,17]。但是，中國2014年P(guān)PP疫情中的PPRV毒株與早前中國西藏毒株未劃分在同一個(gè)小分支，且基因序列也有差異。由于2010-2013年11月中國未發(fā)生過小反芻獸疫疫情；同時(shí)，考慮新疆少數(shù)民族生活習(xí)慣及消費(fèi)羊產(chǎn)品的群體、新疆豐富的野生小反芻獸資源、新疆的地理位置和國家對新疆大力發(fā)展養(yǎng)羊業(yè)的政策等導(dǎo)致羊只流動(dòng)頻繁、接觸性大等因素，分析得出，2014年P(guān)PP疫情為西藏疫情的延續(xù)或者是西藏毒株為適應(yīng)中國環(huán)境和易感動(dòng)物而發(fā)生一些基因變異的可能性極小。

另外，China/XJYL/2013株與國外毒株中的2010年孟加拉國毒株遺傳演化關(guān)系最近，這與王志亮等[17]報(bào)道不一致，可能與該文獻(xiàn)PPRV進(jìn)化分析中采用的是F基因的部分序列有關(guān)，采用具有代表性基因的全序列分析遺傳發(fā)生關(guān)系，可能更具代表性；若能結(jié)合宿主、發(fā)生時(shí)間和地點(diǎn)進(jìn)行分析則更有意義。

F蛋白氨基酸同源性分析顯示，2013年中國新疆PPRV毒株F蛋白與2014年P(guān)PR疫情大流行中國其他毒株氨基酸同源性最高，與之前中國西藏毒株氨基酸同源性也較高，與國外毒株氨基酸同源性最高為印度不同年代的3個(gè)毒株，且同源性高于其與中國西藏毒株的同源性。F蛋白氨基酸同源性比對結(jié)果與核苷酸比對結(jié)果相似。

綜合分析，此次PPR大流行疫情可能是周邊國家毒株的跨境傳入。至于更全面更準(zhǔn)確的結(jié)果，有待于不同時(shí)空分布的更多PPRV毒株F基因全序列的報(bào)道。

3.2國外PPR疫情特征

本研究顯示，Ⅰ系和Ⅱ系毒株分類與前人報(bào)道相符[1,16]；Ⅲ系為亞洲鄰國阿聯(lián)酋和阿曼各1個(gè)毒株、非洲東北部埃塞俄比亞1994年1個(gè)毒株、非洲東部的肯尼亞和烏干達(dá)各1個(gè)毒株構(gòu)成；Ⅳ系包括21個(gè)毒株，其中18個(gè)亞洲毒株，1個(gè)非洲東北部和2個(gè)非洲西北部毒株。由此可知，在系統(tǒng)分類中出現(xiàn)一個(gè)地區(qū)或一個(gè)國家存在多個(gè)PRRV基因系的情況。若一些國家出現(xiàn)多個(gè)基因系毒株或新奇毒株，可能與其在自身進(jìn)化中受到獨(dú)特的易感動(dòng)物、生存環(huán)境等的影響有關(guān)。

綜上所述，周邊國家疫情、養(yǎng)殖特點(diǎn)、獸醫(yī)衛(wèi)生制度、動(dòng)物及動(dòng)物產(chǎn)品生產(chǎn)條件等均對中國畜牧業(yè)發(fā)展和疫情發(fā)生有一定影響，尤其對新疆等周臨多國的中國邊境地區(qū)有明顯影響。此外，中國的引種制度、活畜頻繁交易而檢疫實(shí)際操作不嚴(yán)格或不作為等，均可能與中國多種疫情的發(fā)生有至關(guān)重要的關(guān)系。

4結(jié) 論

2013年中國新疆伊犁地區(qū)PPRV株China/XJYL/2013F基因與2014年P(guān)PR疫情大流行中國河南和北京PPRV毒株F基因核苷酸同源性最高，與肯尼亞2011年毒株KN5/2011核苷酸同源性最低，與國外毒株核苷酸同源性最高為印度1994年Izatnagar-94株。2013年中國新疆伊犁地區(qū)毒株屬于Ⅳ系毒株，該毒株可能是由周邊國家傳入。

參考文獻(xiàn)Reference：

[1]BALAMURUGAN V,SEN A,VENKATESAN G,etal.Peste des petits ruminants virus detected in tissues from an Asiatic lion (Panthera Leo Persica) belongs to Asian lineage iv[J].JournalofVeterinaryScience,2012,13(2):203-206.

[2]劉永宏,曹勝波,趙麗,等.中國小反芻獸疫疫情分析[J].西北農(nóng)業(yè)學(xué)報(bào),2014,23(9):19-26.

LIU Y H,CAO SH B,ZHAO L,etal.Analysis of epidemic Chinese Peste des Petits Ruminants[J].ActaAgriculturaeBoreali-OccidentalisSinica,2014,23(9):19-26(in Chinese with English abstract).

[3]GARGADENNEC L,LALANNE A.La peste des petits ruminants[J].BulletindesServicesZooTechniqueetdesEpizootiedel’AfriqueOccidentaleFranfiaise,1942,5:16-21.

[4]羅靜,何宏軒.小反芻獸疫病毒的分子生物學(xué)特性及其在全球的流行[J].河北師范大學(xué)學(xué)報(bào)(自然科學(xué)版),2009,3(4):533-540.

LUO J,HE H X.The molecular biological characterization of peste des petits ruminants virus and its prevalence in the world [J].JournalofHebeiNormalUniversity(NaturalScienceEdition),2009,3(4):533-540(in Chinese with English abstract).

[5]蔣梅,楊仕標(biāo),張念祖.小反芻獸疫的流行趨勢與防控[J].動(dòng)物醫(yī)學(xué)進(jìn)展,2007,28(S1):88-91.

JIANG M,YANG SH B,ZHANG N Z.Epidemic trend and control of peste des petits ruminants[J].ProgressinVeterinaryMedicine,2007,28(S1):88-91(in Chinese).

[6]BAO J,LI L,WANG Z,etal.Development of one-step real-time RT-PCR assay for detection and quantitation of peste des petits ruminants virus[J].JournalofVirologicalMethods,2008,148(1-2):232-236.

[7]次仁羅布．西藏阿里地區(qū)小反芻獸疫的流行情況與綜合防控[J].山東畜牧獸醫(yī),2013,5(34):44-45.

CIRLB.The epidemic situation and comprehensive prevention and control of Tibet Ali region of peste des petits ruminants[J].ShandongJournalofAnimalScienceandVeterinaryMedicine,2013,5(34):44-45(in Chinese).

[8]WANG Z,BAO J,WU X,etal.Peste des petits ruminants virus in Tibet,China[J].EmergingInfectiousDiseases,2009,15(2):299-301.

[9]王樂元,次真,吳國珍,等.中國西藏小反芻獸疫的發(fā)生狀況與防控[J].畜牧獸醫(yī)學(xué)報(bào),2011,42(5):717-720.

WANG L Y,TSED,WU G ZH,etal.Occurrence and control strategy of peste des petits ruminants in Tibet Autonomous region of China[J].ChineseJournalofAnimalandVeterinarySciences,2011,42(5):717-720(in Chinese with English abstract).

[10]劉永宏,趙麗,曹勝波,等.小反芻獸疫在中國的流行趨勢及應(yīng)對措施[J].動(dòng)物醫(yī)學(xué)進(jìn)展,2015,36(1):110-113.

LIU Y H,ZHAO L,CAO SH B,etal.Popular trends and countermeasures of peste des petits ruminants in China[J].ProgressinVeterinaryMedicine,2015,36(1):110-113(in Chinese with English abstract).

[11]WU X,LI L,LI J,etal.Peste des petits ruminants viruses re-emerging in China,2013-2014[J].Transboundary&EmergingDiseases,2015.

[12]包靜月,王志亮,李林,等.我國西藏小反芻獸疫病毒China/Tib/Gej/07-30核衣殼蛋白基因和基因組啟動(dòng)子區(qū)的分子特征分析[J].病毒學(xué)報(bào),2008,24(6):464-471.

BAO J Y,WANG ZH L,LI L,etal.Sequence analysis of the nucleocapsid gene and genome promoter region of peste des petits ruminants virus of Chinese origin[J].ChineseJournalofVirology,2008,24(6):464-471(in Chinese with English abstract).

[13]包靜月,李林,王志亮,等．一步法實(shí)時(shí)定量RT-PCR檢測小反芻獸疫病毒方法的建立[J].中國動(dòng)物檢疫,2007,24(8):27-29.

BAO J Y,LI L,WANG ZH L,etal.One-step real-time quantitative reverse transcriptase PCR (qRT- PCR) for the detection of peste des petits ruminants virus[J].ChineseJournalofAnimalHealthInspection,2007,24(8):27-29(in Chinese with English abstract).

[14]BAILEY D,BANYARD A,DASH P,etal.Full genome sequence of peste des petits ruminants virus,a member of the Morbillivirus genus[J].VirusResearch,2005,110(1/2):119-124.

[15]BAO J Y,WANG Q H,ZHANG Y Q,etal.Complete genome sequence of a novel variant strain of peste des petits ruminants virus,China/XJYL/2013[J].GenomeAnnouncements,2014,2(5):1-2.

[16]SHAILA M S,SHAMAKI D,FORSYTH M A,etal.Geographic distribution and epidemiology of peste des petits ruminants virus[J].VirusResearch,1996,43 (2):149-153.

[17]王志亮,包靜月,吳曉東,等.我國首例小反芻獸疫診斷報(bào)告[J].中國動(dòng)物檢疫,2007,24(8):24-26.

WANG ZH L,BAO J Y,WU X D,etal.Diagnosis of the first outbreak of peste des petits ruminants in Tibet [J].ChineseJournalofAnimalHealthInspection,2007,24(8):24-26(in Chinese).

Received 2015-07-30Returned2015-11-20

Foundation itemCollege Students Innovation and Entrepreneurship Training Program of China (No.201510757011); College Students Innovation Program of College of Animal Science, Tarim University (No.DKY2014004,No.DKY2014009); Joint Fund Project of Huangzhou Agricultural University and Tarim University (No.HNTDLH1404).

First authorMA Hailu, male,undergraduate student.Research area:infectious diseases and immune pathology.E-mail: 420723205@qq.com

ZHAO Li, female,Ph.D,associate professor.Research area:parasitic diseases and immune pathology.E-mail: zhaolidky@126.com

(責(zé)任編輯：郭柏壽Responsible editor:GUO Baishou)

Sequence Analysis ofFGene of Peste Des Petits Ruminants in Xinjiang

MA Hailu1, LIU Yonghong1, ZHAO Li1, CAO Shengbo2, LIU Xuefeng3,WANG Haiguo4, LI Yueyang1, MA Shiqiang1,JIANG Xiumei1, ZHANG Baojun1,Sayifujiama Kayisaer1, Nuerguli Tuergong1, LIAO Qiuping1and JIAO Haihong1

(1.College of Animal Science, Tarim University, Key Laboratory of Tarim Animanl Husbandry Science and Technology of Xinjiang Production & Construction Corps, Alar Xinjiang843300, China; 2.College of Veterinary Medicine,Huazhong Agricultural University, Wuhan430070, China; 3.The 1st Pasturage Farm, Agricultural Division 14 of Xinjiang Production and Construction Corps, Hetian Xinjiang848306, China; 4.Hongqi Farm Production Department, Third Divisions, Xinjiang Production and Construction Corps, Artush Xinjiang845350, China)

AbstractIn order to analyze the molecular evolution characteristic of peste des petits ruminants virus (PPRV) and the epidemic disseminated route in Xinjiang.We downloaded the complete sequences of F gene of PPRV from GenBank of all the countries in the world.We used DNAStar software and combined it with spatial and temporal distribution of strains from NCBI to analyze its sequence.The results showed that the nucleotide homology of Xinjiang strain was the highest with the epidemic strainsof other places in China, the nucleotide homology with Tibet strain which was epidemic before, and with Chinese vaccine strain used were 97.7% and 92.9% respectively.The nucleotide homology with Kenya strain was the lowest, with India strain was the highest among foreign strains.The phylogenetic analysis of F gene showed that Xinjiang strain belongs to the gene series Ⅳ.The amino acid homology of F protein was consistent nearly with the nucleotide homology, but was the lowest with ICV89 strain from Ivory Coast.There are certain variation in the PPRV strain, it may be introduced into the neighboring countries.

Key wordsChina; Xinjiang; PPR; F gene

Corresponding authorLIU Yonghong,male,associate professor.Research area: infectious diseases and immune pathology.E-mail: lyhdky@126.com

中圖分類號S855.3

文獻(xiàn)標(biāo)志碼A

文章編號1004-1389(2016)03-0328-07

通信作者：劉永宏，男，副教授，研究方向?yàn)閭魅静∨c免疫病理學(xué)。E-mail: lyhdky@126.com

基金項(xiàng)目：國家級大學(xué)生創(chuàng)新創(chuàng)業(yè)訓(xùn)練計(jì)劃(201510757011)；塔里木大學(xué)動(dòng)物科學(xué)學(xué)院大學(xué)生創(chuàng)新項(xiàng)目(DKY2014004，DKY2014009)；華中農(nóng)業(yè)大學(xué)塔里木大學(xué)科研聯(lián)合基金(HNTDLH1404)。

收稿日期：2015-07-30修回日期：2015-11-20

網(wǎng)絡(luò)出版日期：2016-03-06

網(wǎng)絡(luò)出版地址：http://www.cnki.net/kcms/detail/61.1220.S.20160306.1610.004.html

第一作者：馬海璐，男，本科生，研究方向?yàn)閭魅静∨c免疫病理學(xué)。E-mail：420723205@qq.com

趙麗，女，博士，副教授，研究方向?yàn)榧纳x病與免疫病理學(xué)。E-mail: zhaolidky@126.com