王宇清，刀金威，朱麗慧，杜曉利，陸 珂，李 浩，劉金明，程國鋒

（中國農(nóng)業(yè)科學(xué)院上海獸醫(yī)研究所 農(nóng)業(yè)部動物寄生蟲學(xué)重點開放實驗室，上海200241）

兔感染日本血吸蟲后血漿中宿主源循環(huán)microRNAs變化研究

王宇清，刀金威，朱麗慧，杜曉利，陸 珂，李 浩，劉金明，程國鋒

（中國農(nóng)業(yè)科學(xué)院上海獸醫(yī)研究所 農(nóng)業(yè)部動物寄生蟲學(xué)重點開放實驗室，上海200241）

循環(huán)性microRNAs不僅參與生命體的重要活動調(diào)控，還有望成為某些疾病診斷的新型標(biāo)志分子。為深入探討宿主與血吸蟲之間的相互作用機(jī)制并發(fā)現(xiàn)新型診斷標(biāo)志分子，我們對前期獲得的感染和未感染新西蘭大白兔血漿小RNA的Solexia測序結(jié)果進(jìn)行生物信息學(xué)分析，對差異表達(dá)的microRNAs用實時定量PCR進(jìn)行驗證，同時對這些差異表達(dá)的microRNAs的靶基因進(jìn)行了分析。在剔除蟲源性序列后，我們發(fā)現(xiàn)了感染和未感染日本血吸蟲的新西蘭大白兔體內(nèi)的循環(huán)性宿主源microRNAs 173條，選取了10個差異表達(dá)的microRNAs利用實時定量PCR進(jìn)行驗證分析，結(jié)果表明所選取的10個microRNAs表達(dá)與Solex測序結(jié)果一致。靶基因的預(yù)測與分析結(jié)果表明感染血吸蟲兔血漿中差異表達(dá)microRNAs的靶基因主要參與細(xì)胞進(jìn)程、代謝調(diào)控、生命調(diào)控、發(fā)育進(jìn)程等生命進(jìn)程調(diào)控。本研究結(jié)果為深入開展日本血吸蟲與宿主的互作及其寄生的分子機(jī)理提供了基礎(chǔ)。

日本血吸蟲；循環(huán)microRNA；靶基因；功能分析

microRNAs（miRNAs）早在1990年就在線蟲體內(nèi)被發(fā)現(xiàn)［1-4］，隨后大量研究表明不同物種和組織中均有miRNAs參與細(xì)胞功能調(diào)控［5，6］。最近研究表明miRNAs可在人體的血液、淋巴液和尿液中存在，且這些miRNAs可通過循環(huán)系統(tǒng)，遞送到機(jī)體的其他組織和器官，進(jìn)而發(fā)揮生物學(xué)調(diào)控功能［7-9］。因此，循環(huán)性miRNAs分子不僅可成為某些疾病診斷的新型標(biāo)志分子，也可能成為疾病病理損害的關(guān)鍵調(diào)控分子，故而循環(huán)miRNAs本身以及其所調(diào)控的靶基因都可能作為潛在的治療疾病的靶位點。

我們課題組前期研究表明日本血吸蟲感染宿主后可在宿主血漿中檢測到血吸蟲特異miRNAs（Bantam、miR-3479、miR-10、miR-3096和miR-0001）［10］。隨后，Hoy等［11］研究也表明感染曼氏血吸蟲的小鼠血清中存在血吸蟲特異的miRNAs（miR-277、miR-3479-3p和bantam）。He等［12］在感染了日本血吸蟲的血清宿主中發(fā)現(xiàn)了miR-223表達(dá)量升高，經(jīng)過藥物治療后其表達(dá)水平降低，提示了miR-223可能是日本血吸蟲診斷標(biāo)志物。Zhu等［13］進(jìn)一步研究表明在感染了日本血吸蟲的小鼠血漿中有明顯的循環(huán)性miRNAs的變化，特別是miR-706以及miR-134-5p，并預(yù)測其參與Caspase-3、Creb1基因的調(diào)控，提示循環(huán)性microRNAs可能在血吸蟲的病理損害中發(fā)揮重要作用。故而，針對感染日本血吸蟲后宿主循環(huán)性miRNAs的研究不僅可發(fā)現(xiàn)血吸蟲病診斷新的標(biāo)志分子，而且可從miRNAs視角闡述病原與宿主相互作用機(jī)制［14］。為此，我們利用生物信息學(xué)對前期獲得的新西蘭大白兔感染血吸蟲后血漿miRNAs的Solexia測序結(jié)果進(jìn)行分析，研究在新西蘭大白兔體內(nèi)與血吸蟲感染相關(guān)的宿主源循環(huán)性miRNAs，并預(yù)測這些miRNAs所調(diào)控的靶基因及其功能。

1 材料和方法

1.1 材料

1.1.1 實驗動物以及生物耗材 本研究的動物相關(guān)實驗由中國農(nóng)業(yè)科學(xué)院上海獸醫(yī)研究所動物保健及倫理委員會許可。本研究所用的釘螺由上海獸醫(yī)研究所釘螺室提供，新西蘭大白兔購自上海羅涇飛達(dá)實驗動物公司。實時定量PCR所用的引物均由Sunny有限責(zé)任公司負(fù)責(zé)合成。

1.1.2 生化試劑及儀器 PCR Mix（2X）、PrimeScript RT reagent Kit（DRRO37A）和SYBR Premix Ex Taq（RR420A）購自寶生物工程（大連）有限公司（大連TaKaRa公司）；血漿小RNA提取試劑盒mirVana Protein and RNA Isolation System購自Invitrogen公司；PCR擴(kuò)增儀購自伯樂生命醫(yī)學(xué)產(chǎn)品（上海）有限公司；離心機(jī)和Eppendorf 快速實時熒光定量PCR儀購自Eppendorf中國有限公司。

1.1.3 Solexia測序數(shù)據(jù)庫 本文分析的感染血吸蟲后新西蘭大白兔血漿miRNA數(shù)據(jù)庫P1和P2以及未感染血吸蟲的新西蘭大白兔血漿miRNA數(shù)據(jù)庫為本實驗室通過建庫后測序獲得［11］，其中P1為感染后14、28 d的新西蘭大白兔混合血漿中的小RNA庫；P2為感染后10、14、28、32 d的新西蘭大白兔混合血漿的小RNA庫。

1.2 方法

1.2.1 血液收集、血漿RNA的提取和小RNA庫的構(gòu)建血液材料的收集、血漿的分離以及血漿RNA的提取和小RNA庫的構(gòu)建參照文獻(xiàn)［11］。

1.2.2 感染和未感染新西蘭大白兔宿主源循環(huán)小RNA分析 將Solexia測序數(shù)據(jù)庫P1和P2以及未感染血吸蟲的新西蘭大白兔血漿miRNA數(shù)據(jù)庫數(shù)據(jù)利用SOAP（http∶//soap.genomics.org.cn/soap1/）與家兔基因組數(shù)據(jù)庫（ftp∶//ftp.ensembl.org/pub/release-65/fasta/ oryctolagus_cuniculus）進(jìn)行匹配，列出與基因組完全匹配的序列。在Rfam（Version 10.1，http∶//rfam. xfam.org/）中收集的非編碼RNA數(shù)據(jù)庫與NCBI中GeneBank數(shù)據(jù)（http∶//www.ncbi.nlm.nih.gov/）進(jìn)行BLAST分析，剔除所有與血吸蟲序列同源的序列后對這些小RNA進(jìn)行分類，比如rRNA，tRNA，snRNA以及snoRNA，對那些沒有匹配的小RNA進(jìn)一步與miRbase（http∶//www.mirbase.org/，Vision 20）進(jìn)行比對，以確定其中哪些是已知的miRNA或是同源的miRNA，并將沒有注釋的小RNAs進(jìn)行新miRNA的預(yù)測。然后將感染日本血吸蟲的新西蘭大白兔血漿的測序結(jié)果與未感染的大白兔血漿的小RNA測序結(jié)果進(jìn)行比較分析，找到P1和P2中特有的小RNA。

1.2.3 MicroRNAs靶基因的預(yù)測和功能分析 利用miRWalk 在線工具（http∶//www.umm.uni-heidelberg. de/apps/zmf/mirwalk/）預(yù)測miRNAs的靶基因；利用DAVID（http∶//david.abcc.ncifcrf.gov/）數(shù)據(jù)庫分析其靶基因的功能。

1.2.4 差異miRNAs的qRT-PCR驗證分析 將新西蘭大白兔感染了約1000條日本血吸蟲尾蚴，參考文獻(xiàn)［11］方法收集感染日本血吸蟲后7、14、28、32 d的新西蘭大白兔血液并提取血漿RNA，反轉(zhuǎn)錄所使用的特異性引物以及實時定量PCR引物見表1。實時定量PCR通用下游引物：5'-GCAGGGTCCGAGGTATTC-3'。逆轉(zhuǎn)錄所使用的反應(yīng)體系：5×PrimeScriptTMBuffer （for Real Time）2μL、PrimeScriptTMRT Enzyme MixⅠ0.5 μL、Specific primer 0.5 μL、Total RNA 25 ng，然后用RNase free H2O 補(bǔ)充至10 μL；反應(yīng)條件：42℃ 30 min，85℃ 5 s，4℃保存。實時定量PCR的反應(yīng)體系：cDNA 1μL，SBYR?Premix Ex TaqTMⅡ（2X）10 μL，上游特異性引物0.5 μL，通用下游引物0.5 μL，然后用RNase free H2O 定容至20 μL；反應(yīng)條件：95℃ 30 s，然后（95℃ 5 s，60.5℃ 30 s）45個循環(huán)。

表1 逆轉(zhuǎn)錄和實時定量PCR引物Table 1 Primers for reverse PCR and Real-Time PCR

2 結(jié)果

2.1 感染和未感染兔差異表達(dá)循環(huán)性宿主源已知小

RNA 將兩個生物學(xué)重復(fù)P1和P2的Solexa數(shù)據(jù)分別與未感染的陰性對照血漿小RNA測序結(jié)果比較分析。結(jié)果表明：P1中有721條已知miRNAs，P2中有1101條。其中，在P1和P2兩次實驗中共有的miRNAs 為173條。進(jìn)一步分析表明感染后有64個已知miRNAs在P1和P2呈現(xiàn)差異表達(dá)（差異倍數(shù)在兩倍以上），其中12個下調(diào)，52個上調(diào)（圖1）。

2.2 感染后宿主源特異性循環(huán)miRNAs分布 將P1與P2分別與未感染的家兔血漿小RNAs庫進(jìn)行比對，獲得感染后特有的小RNAs并進(jìn)行分類分析（圖2A和2B）。對未注釋的特異性的小RNA進(jìn)行新miRNAs預(yù)測，結(jié)果表明P1有191條新miRNAs，P2中有190條新miRNAs。

圖1 宿主血漿中已知miRNAs的差異表達(dá)Fig. 1 Differential levels of circulating miRNAs associated with Schistosoma japonicum

2.3 感染兔特異性循環(huán)miRNAs靶基因預(yù)測 對感染兔特異性表達(dá)的宿主源循環(huán)miRNAs以人的數(shù)據(jù)庫為參考進(jìn)行靶基因預(yù)測，按五個軟件都能預(yù)測到的靶基因為標(biāo)準(zhǔn)，結(jié)果在124個miRNAs中，共獲得了18個miRNAs以及它們所調(diào)控的共416條靶基因（表2）。

圖2 P1（A）和P2（B）中宿主源特異性小RNA分布Fig. 2 The distribution of specifi c miRNAs from P1（A） and P2（B）

表2 差異表達(dá)的miRNAs靶基因預(yù)測Table 2 Target genes prediction for differentially expressed miRNAs

2.4 靶基因功能的生物信息學(xué)分析 利用David數(shù)據(jù)庫提供的在線工具Pather進(jìn)行了靶基因GO聚類分析，結(jié)果他們表明主要參與包括細(xì)胞進(jìn)程、代謝調(diào)控、生命調(diào)控、發(fā)育進(jìn)程等生物學(xué)進(jìn)程。分子功能分析表明他們主要參與結(jié)合、催化活性等分子功能（圖3A、B）。KEGG代謝調(diào)控分析表明，這些miRNAs主要與信號調(diào)節(jié)通路相關(guān)，13個靶基因參與了整聯(lián)素信號通路，9個參與了p53 通路，8個參與了促性腺釋放激素受體通路等（圖4）。

圖 3 生物學(xué)進(jìn)程GO聚類分析（A）和分子功能GO聚類分析（B）Fig. 3 GO analysis of biological process （A） and molecular function （B）

圖4 KEGG代謝通路分析Fig.4 The KEGG pathway analysis

2.5 qRT-PCR的驗證 將選取的10個（9個上調(diào)和1個下調(diào)miRNAs）與血吸蟲感染相關(guān)的宿主源miRNAs利用qRT-PCR進(jìn)行驗證分析。結(jié)果表明這些miRNAs的表達(dá)變化均與Solexia測序結(jié)果相一致（圖5）。

3 討論

圖5 差異表達(dá)的miRNAs的qRT-PCR驗證分析Fig. 5 The qRT- PCR result of selected differentially expressed miRNAs

宿主在感染了日本血吸蟲之后其體內(nèi)會發(fā)生一系列的變化，其中宿主體內(nèi)循環(huán)miRNA的變化就是其中之一［12，15］。通過生物信息學(xué)分析我們篩選出感染兔和未感染兔差異表達(dá)的miRNAs以及特異表達(dá)的miRNAs，總共64條miRNAs，其中上調(diào)52條，下調(diào)12條。通過實時定量PCR，進(jìn)一步驗證了一些miRNAs的確在感染后宿主體內(nèi)存在差異表達(dá)。對這些miRNAs進(jìn)行了靶基因的預(yù)測，發(fā)現(xiàn)其中18個miRNA的預(yù)測靶基因共416個。Zhu等［13］研究了感染日本血吸蟲后小鼠體內(nèi)循環(huán)性miRNAs的變化，發(fā)現(xiàn)總共有11個miRNAs下調(diào)，25個miRNAs上調(diào)，其中miR-125b-1-3p在感染后上調(diào)，但本研究發(fā)現(xiàn)它在感染血吸蟲的新西蘭大白兔中為下調(diào)，二者剛好相反，提示宿主感染血吸蟲后循環(huán)性miRNAs可能存在宿主間的差異性，這有待后續(xù)研究進(jìn)一步查明。

為了更好的推測這些miRNAs的功能，我們預(yù)測了miRNAs的靶基因并進(jìn)行了KEGG代謝通路分析，發(fā)現(xiàn)這些差異表達(dá)的循環(huán)性miRNAs主要參與許多重要的信號傳導(dǎo)調(diào)控，比如p38 MAPK、P53 pathway等，這些信號通路與多種疾病的發(fā)展進(jìn)程密切相關(guān)。測序結(jié)果和qRT-PCR驗證均表明在感染日本血吸蟲的新西蘭大白兔的循環(huán)性miRNAs中miR-10是上調(diào)的，有研究指出在多種癌癥病理過程中miR-10是上調(diào)的，比如miR-10a在結(jié)腸癌中被發(fā)現(xiàn)，miR-10b在胰腺癌和乳腺癌中被發(fā)現(xiàn)［16］。因此對感染日本血吸蟲的宿主體內(nèi)循環(huán)microRNAs的研究，不僅為發(fā)現(xiàn)血吸蟲病的新標(biāo)志物提供線索，進(jìn)而被應(yīng)用于診斷，同時也加深了我們對宿主與寄生蟲之間相互作用機(jī)制的理解。

［1］Lee R C， Feinbaum R， Ambros V. The C. elegans heterochronic gene lin-4 encodes small RNAs with antisense complementarity to lin-14［J］. Cell， 1993， 75（5）∶843-854.

［2］Lagos-Quintana M， Rauhut R， Lendeckel W， et al. Identification of novel genes coding for small expressed RNAs［J］. Science， 2001， 294（5543）∶ 853-858.

［3］Lau N C， Lim L P， Weinstein E G， et al. An abundant class of tiny RNAs with probable regulatory roles in Caenorhabditis elegans［J］. Science， 2001， 294（5543）∶858-862.

［4］Ambros V， Lee R C， Lavanway A， et al. MicroRNAs and other tiny endogenous RNAs in C. elegans［J］. Curr Biol，2003， 13（10）∶ 807-818.

［5］Bartel D P. MicroRNAs∶ genomics， biogenesis，mechanism， and function［J］. Cell， 2004， 116（2）∶ 281-297.

［6］Friedman R C， Farh K K， Burge C B， et al. Most mammalian mRNAs are conserved targets of microRNAs［J］. Genome Res， 2009， 19（1）∶ 92-105.

［7］Taylor D D， Gercel-Taylor C. MicroRNA signatures of tumor-derived exosomes as diagnostic biomarkers of ovarian cancer［J］. Gynecol Oncol， 2008， 110（1）∶ 13-21.

［8］Valadi H， Ekstr?m K， Bossios A. Exosome-mediated transfer of mRNAs and microRNAs is a novel mechanism of genetic exchange between cells［J］. Nat Cell Biol， 2007，9（6）∶ 654-659.

［9］Denzer K， Kleijmeer M J， Heijnen H F， et al. Exosome∶from internal vesicle of the multivesicular body to intercellular signaling device［J］. J Cell Sci， 2000，113Pt19∶ 3365-3374.

［10］Cheng G， Luo R， Hu C， et al. Deep sequencing-based identification of pathogen-specific microRNAs in the plasma of rabbits infected with Schistosoma japonicum［J］. Parasitology， 2013， 140（14）∶ 1751-1761.

［11］Hoy A M， Lundie R J， Ivens A， et al. Parasite-Derived MicroRNAs in Host Serum As Novel Biomarkers of Helminth Infection ［J］. Plos Negl Trop Dis， 2014， 8（2）∶85-88.

［12］He X， Sai X， Chen C， et al. Host serum miR-223 is a potential new biomarker for Schistosoma japonicum infection and the response to chemotherapy ［J］. Parasit Vectors， 2013， 6∶ 272.

［13］Zhu L， Dao J， Du X， et al. Altered levels of circulating miRNAs are associated Schistosoma japonicum infection in mice［J］. Parasites Vectors， 2015， 8∶196.

［14］Xue X， Sun J， Zhang Q， et al. Identification and characterization of novel microRNAs from Schistosoma japonicum［J］. PLoS One， 2008， 3（12）∶ e4034.

［15］Liu F， Lu J， Hu W， et al. New perspectives on hostparasite interplay by comparative transcriptomic and proteomic analyses of Schistosoma japonicum［J］. PLoS Pathog， 2006， 2（4）∶ e29.

［16］Lund A H. miR-10 in development and cancer［J］. Cell Death Differ， 2010， 17（2）∶ 209-214.

STUDIES ON RABBIT CIRCULATING MICRORNAS ASSOCIATED WITH SCHISTOSOMA JAPONICUM INFECTION

WANG Yu-qing， DAO Jin-wei， ZHU Li-hui， DU Xiao-li， LU Ke， LI Hao，LIU Jin-ming， CHENG Guo-feng
（Key Laboratory of Animal Parasitology， Ministry of Agriculture， Shanghai Veterinary Research Institute， CAAS， Shanghai 200241， China）

Circulating microRNAs are not only involved in the regulation of many life processes， but also considered to be a novel class of biomarkers for diagnosis of diseases. To better understand the interaction mechanisms between hosts and Schistosoma japonicum，the microRNA populations in infected New Zealand rabbits were further analyzed based on our previous Solexia data from infected and uninfected rabbits. Total 173 rabbit microRNAs associated with Schistosoma japonicum infection were obtained， from which 10 microRNAs different between infected and uninfected rabbits were verifi ed in a qRT-PCR. Bioinformatical analysis also indicated that the target genes of these altered microRNAs were putatively involved in many biological processes such as cellular process， metabolic process， biological regulation， developmental process and catalytic activity. The results obtained here provided a basis for further research on interaction between hosts and schistosomes.

Schistosoma japonicum； circulating microRNAs； target genes； functional analysis

S852.735

A

1674-6422（2015）05-0058-07

2015-05-22

國家自然科學(xué)基金（30901068）

王宇清，男，碩士研究生，預(yù)防獸醫(yī)學(xué)專業(yè)

程國鋒，E-mail：chenggfeng@yahoo.com.cn