• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    上海地方豬種抗腹瀉基因與免疫相關因子分析

    2024-06-28 21:30:07周金勇孫玲偉張克勤戴建軍張德福王振國主性吳彩鳳
    中國農(nóng)業(yè)科技導報 2024年4期

    周金勇 孫玲偉 張克勤 戴建軍 張德?!⊥跽駠≈餍浴遣束P

    摘要:為鑒定上海地方豬種上海白豬、楓涇豬和沙烏頭豬中抗腹瀉基因抗原加工相關轉(zhuǎn)運體1(transporterassociatedwith antigen processing 1, TAP1)、α(1, 2)巖藻糖轉(zhuǎn)移酶1[α( 1, 2)-fucosyltransferase 1, FUT1]、天然抗性相關的巨噬細胞蛋白1(natural resistance-associated macrophage protein 1, NRAMP1)、黏蛋白4(mucin 4,MUC4)和黏蛋白13(mucin 13, MUC13)基因的有效遺傳標記,為上海地方豬種資源特性提供參考,應用聚合酶鏈式反應-限制性片段長度多態(tài)性(polymerase chain reaction restriction fragment length polymorphism, PCRRFLR)和序列測序分析上述基因的多態(tài)性,結合部分免疫相關因子,探討對3個上海地方豬種免疫力的影響。結果顯示,3個豬種中TAP1 和MUC4 基因均具有抗腹瀉基因型GG,而僅在楓涇豬和沙烏頭豬中檢測到MUC13基因抗腹瀉基因型GG,在3個豬種中均未檢測到FUT1 和NRAMP1 基因抗腹瀉基因型AA。上海白豬、楓涇豬和沙烏頭豬TAP1 基因處于中度多態(tài),上海白豬的MUC4 基因處于低度多態(tài),楓涇豬和沙烏頭豬的MUC4 基因處于中度多態(tài),上海白豬MUC13 基因處于中度多態(tài),其中上海白豬和沙烏頭豬的TAP1 基因不滿足Hardy-Weinberg平衡,楓涇豬的MUC4 基因不滿足Hardy-Weinberg平衡。上海白豬MUC13 基因AA型的白細胞介素12(interleukin-12, IL-12)水平顯著高于AG 型,楓涇豬TAP1 基因AA 型的腫瘤壞死因子α(tumor necrosisfactor, TNF-α)指標顯著高于GG和AG型,沙烏頭豬TAP1 基因AG型的IL-12指標顯著高于GG型。以上研究結果對上海白豬、楓涇豬和沙烏頭豬的抗腹瀉育種和分子選育具有一定指導意義,為今后上海各地方豬種抗腹瀉育種工作奠定了基礎。

    關鍵詞:上海白豬;楓涇豬;沙烏頭豬;抗腹瀉基因;免疫因子

    doi:10.13304/j.nykjdb.2023.0229

    中圖分類號:S828.2 文獻標志碼:A 文章編號:10080864(2024)04012816

    根據(jù)2022年10月國家統(tǒng)計局數(shù)據(jù),我國生豬存欄44 394萬頭,擁有相當高的生豬存欄量,而豬腹瀉病常年頻發(fā),致死率高,是嚴重制約我國乃至全世界養(yǎng)豬業(yè)發(fā)展的一大難題[1]。當前,對豬群的抗腹瀉相關基因通過分子育種技術進行其等位基因頻率的檢測與分析,從而用于指導豬群育種,從育種角度提高豬群整體的抗腹瀉能力,進而從根本上解決生豬腹瀉病帶來的經(jīng)濟損失。

    產(chǎn)腸毒素大腸桿菌(enterotoxi-genic escherichiacoli, ETEC)F18的候選基因主要是α (1, 2)巖藻糖轉(zhuǎn)移酶基因1[α (1, 2)-fucosyltransferase1, FUT1]和抗原處理相關轉(zhuǎn)運體基因(transporter-associatedwith antigen processing 1, TAP1)[23]。Meijerink等[4]研究發(fā)現(xiàn),瑞士長白豬FUT1 基因的A 型可抗ETEC F18的感染。Kim等[5]進一步證明,F(xiàn)UT1 基因的AA基因型個體的仔豬存活率幾乎是GG個體的2倍。因此,F(xiàn)UT1 基因多態(tài)性可用作選擇程序的有效標記,以提高斷奶后仔豬的存活率。趙喬輝等[6]研究發(fā)現(xiàn),在淮豬、榮昌豬和長白豬等豬種中TAP1 基因的GG型個體有更強的抗ETEC F18感染的能力,可將TAP1 基因作為豬抗腹瀉育種的有效遺傳標記之一。天然抗性相關的巨噬蛋白1(natural resistance-associated macrophage protein 1,NRAMP1)基因能夠影響動物機體的抗病性能,豬NRAMP1 基因AA基因型為抗腹瀉優(yōu)勢基因型[7]。劉艷冬等[8]通過觀測記錄哺乳仔豬的腹瀉情況發(fā)現(xiàn),NRAMP1 在各品種豬中表現(xiàn)出豐富的多態(tài)性,進一步證明NRAMP1 可作為豬抗病育種候選基因。粘蛋白4(mucin 4, MUC4)和粘蛋白13(mucin 13, MUC13)基因是ETEC F4ab/F4ac 受體的重要候選基因[9-11]。Peng 等[12]研究表明,豬對ETEC F4ab/ac 感染的易感性/抗性與MUC4 基因內(nèi)含子17中的G243位點A→G突變顯著相關,且GG 型是抗腹瀉基因型。任軍等[13] 也發(fā)現(xiàn),MUC13 基因的rs319 699 771位點上存在G→A突變,GG型為抗腹瀉基因型,且通過鑒定該位點突變能夠鑒別抗腹瀉個體。因此,MUC4 和MUC13基因被視為F4ab/F4ac 受體最有希望的候選基因。

    本研究對上海白豬、楓涇豬和沙烏頭豬部分抗腹瀉基因(TAP1、FUT1、NRAMP1、MUC4 和MUC13)的基因多態(tài)性進行篩查,并結合部分血液免疫因子,如免疫球蛋白(immunoglobulin, IgG)、白細胞介素2(interleukin-2, IL-2)、白細胞介素6(interleukin-6, IL-6)、白細胞介素10(interleukin-10, IL-10)、白細胞介素12(interleukin-12, IL-12)、腫瘤壞死因子ɑ(tumor necrosis factor, TNF-ɑ)、干擾素γ(interferon γ, IFN-γ)等,進一步分析上述豬種的抗腹瀉能力,以期為提高上述豬種抗腹瀉育種工作提供數(shù)據(jù)支撐。

    1 材料與方法

    1.1 采樣

    本研究中使用的動物為90 頭后備種公豬(200~300日齡),包括3種中國上海本土豬:上海白豬(30頭)、楓涇豬(30頭)和沙烏頭豬(30頭),分別來自上海農(nóng)系白豬種原種場、上海金山楓涇豬保種場和上海崇明沙烏頭農(nóng)業(yè)科技有限公司。試驗豬種圍欄中單獨飼養(yǎng),每天按照標準飼喂,并隨意飲用水,按照《中國豬飼喂標準指南》[14]進行病害控制。每個豬種分別采集血液樣本各30份抗凝血用于提取DNA,30份未抗凝血用于測定部分免疫因子。采集的未抗凝血經(jīng)4 000 r· min-1離心20 min后分離血清,置于-80 ℃冰箱保存。

    1.2 DNA 提取

    按照血液全基因組DNA提取試劑盒[天根生化科技(北京)有限公司]使用說明書,對豬血液樣品進行DNA 提取。利用紫外分光光度計(eppendorf公司)對提取的DNA含量進行測定,選擇含量在80~250 μg· μL-1的樣品進行后續(xù)試驗。用1.5% 瓊脂糖凝膠電泳檢測提取DNA 的質(zhì)量,統(tǒng)一稀釋至50 ng·μL-1,-20 ℃冰箱保存?zhèn)溆谩?/p>

    1.3 聚合酶鏈式反應-限制性酶切酶段長度多態(tài)性分析。

    在NCBI 上查找相關序列,使用PrimerPremier 6.0 軟件設計FUT1、NRAMP1、TAP1、MUC13 和MUC4 基因的引物(表1),由中國上海生工生物工程有限公司合成。

    PCR 反應體系為20 μL,包含template DNA2μL,2× M5 Taq hifi PCR mix 10 μL,primer 1(10 μmol·L-1)0.5 μL,primer 2(10 μmol·L-1)0.5 μL,nuclease-free ddH2O 補充至20 μL。PCR 反應條件:預變性94 ℃ 4 min;94 ℃變性30 s,退火30 s,72 ℃延伸(TAP1 和FUT1 為60 s,循環(huán)30 次;NRAMP1 為45 s,循環(huán)29 次;MUC4 為15 s,循環(huán)34次;MUC13 為45 s,循環(huán)35次);72 ℃延伸5 min,并置于4 ℃保存待用(PCR儀購自美國ABI公司)。

    TAP1 基因選用MboI [寶如億(北京)生物技術有限公司]對PCR產(chǎn)物進行消化,酶切反應總體系為10 μL,包括PCR 擴增產(chǎn)物0.4 μL,restrictionenzyme 0.4 μL(10 U·μL-1),10× NE buffer 1 μL,ddH2O 8.2 μL,置37 ℃恒溫反應2 h,65 ℃溫育20 min后經(jīng)1.5%的瓊脂糖凝膠電泳后分析。

    FUT1 與MUC4 基因選用HhaⅠ [寶如億(北京)生物技術有限公司],NRAMP1 基因選用NdeⅠ[寶如億(北京)生物技術有限公司]對PCR產(chǎn)物進行消化。酶切反應總體系為10 μL,包括PCR擴增產(chǎn)物1 μL(500 ng· μL-1),RE10× buffer 1 μL,acetylated BSA 0.1 μL(10 μg·μL-1),restrictionenzyme 0.25 μL(10 U·μL-1),ddH2O 7.65 μL,置37 ℃恒溫反應2 h,后添加4 μL 6×加載緩沖液經(jīng)1.5%的瓊脂糖凝膠電泳后分析。由于MUC13 無特異性內(nèi)切酶位點,故不進行酶切。

    將上述PCR產(chǎn)物切膠回收純化后,送至上海生工生物工程技術服務有限公司單向測序,測序結果用DNAstar軟件分析。

    1.4 血液免疫因子測定

    采用ELISA免疫因子檢測試劑盒(elabsciencebiotechnology Co., Ltd.)測定樣品吸光度,吸光度值用于繪制標準曲線,用于測量以下免疫因子:IgG(DG50 132P-96T)、IL-2(DG50 032P-96T)、IL-6(DG50 031P-96T)、IL-10(DG50 034P-96T)、IL-12(DG50 055P-96T)、TNF- ɑ(DG50 014P-96T)和INF-γ(DG50 141P-96T)[15-17]。所測指標作為上海白豬、楓涇豬和沙烏頭豬免疫應答和一般抗病力的評價。按照使用說明具體操作方法如下:首先,從室溫平衡20 min后的鋁箔袋中取出所需板條,剩余板條用自封袋密封放回4 ℃;設置標準品孔和樣本孔,標準品孔各加不同含量的標準品50 μL;樣本孔先加待測樣本10 μL,再加樣本稀釋液40 μL;空白孔不加;除空白孔外,標準品孔和樣本孔中每孔加入辣根過氧化物酶(horseradishperoxidase, HRP)標記的檢測抗體100 μL,用封板膜密封反應孔,37 ℃水浴鍋或恒溫箱溫育60 min;棄去液體,吸水紙上拍干,每孔加滿洗滌液,靜置1 min,甩去洗滌液,吸水紙上拍干,如此重復洗板5次(也可用洗板機洗板);每孔加入底物A、B各50 μL,37 ℃避光孵育15 min;每孔加入終止液50 μL,15 min內(nèi)在450 nm波長處測定各孔的OD值。繪制標準曲線:在Excel工作表中,以標準品含量作橫坐標,對應OD值作縱坐標,繪制標準品線性回歸曲線,按曲線方程計算各樣本含量(pg·mL-1)。所有免疫因子指標檢測方法相同。

    1.5 統(tǒng)計分析

    測序結果利用NCBI 系統(tǒng)中BLAST 進行比對,鑒定SNP 位點,根據(jù)Hardy-Weinberg 平衡定律,計算基因的基因型頻率。采用適合性檢驗法,進行Hardy-Weinberg 平衡狀態(tài)、卡方(χ2)獨立性檢驗分析,計算各基因的等位基因和基因型頻率。群體內(nèi)遺傳變異和分化的相關參數(shù)通過popgene3.1 軟件計算獲得。采用cervus 3.0 軟件計算基因的多態(tài)信息含量(polymorphism informationcontent, PIC)、觀測雜合度(observed heterozygosity,Ho)和期望雜合度(expected heterozygosity, He);其中PIC≥0.50時為高度多態(tài),0.25≤PlC<0.50時為中度多態(tài),PIC<0.25時為低度多態(tài)。各性狀值與基因型的關系采用SPSS 14.0軟件進行最小二乘分析,其中P<0.05為顯著相關。

    2 結果與分析

    2.1 基因PCR 擴增結果

    根據(jù)TAP1、FUT1、NRAMP1、MUC4 和MUC13基因的引物對上海白豬、楓涇豬和沙烏頭豬提取的DNA 進行擴增,PCR 產(chǎn)物用1.5% 的瓊脂糖電泳檢測。由圖1 可見,TAP1、FUT1、NRAMP1、MUC4 和MUC13 基因在不同豬種中均擴增出清晰的特異性條帶,片段大小分別為767、421、483、220和230 bp,與預測的擴增片段大小一致。

    2.2 PCR 產(chǎn)物的酶切結果

    由圖2可知,上海白豬、楓涇豬和沙烏頭豬的TAP1 基因PCR產(chǎn)物經(jīng)限制性內(nèi)切酶MboⅠ酶切,經(jīng)過1.5%的瓊脂糖凝膠電泳后得到上海白豬有2種不同基因型(AG和GG),楓涇豬和沙烏頭豬有3種不同的基因型條帶,分別為AA、GG和AG。根據(jù)酶切結果,對PCR產(chǎn)物進行測序發(fā)現(xiàn),在測序峰圖中部分豬種個體的TAP1 基因有1個G→A的單堿基突變,將該位點未突變呈現(xiàn)單峰的基因型定義為AA型,將該位點突變呈現(xiàn)單峰的基因型定義為GG型,呈現(xiàn)雙峰的基因型定義為AG型(圖2D)。GG 型個體的序列與NCBI 中AK396698.1的序列一致,定義為野生型,相比之下發(fā)現(xiàn)GG 型在第729位置有1個G→A 的單堿基突變,這與酶切結果一致(圖2E)。

    上述豬種FUT1 基因的PCR產(chǎn)物經(jīng)限制性內(nèi)切酶MboⅠ酶切,經(jīng)1.5%的瓊脂糖凝膠電泳后得到3個豬種的FUT1 基因都只有1種基因型條帶,為GG表型(241/93/87 bp)(圖3)。根據(jù)酶切結果對其進行驗證,對PCR產(chǎn)物進行測序,對測序峰圖分析發(fā)現(xiàn),在測序峰圖中所有豬種個體的FUT1 基因都不存在單堿基突變,將該位點未突變呈現(xiàn)單峰的基因型定義為GG型(圖3D)。所有豬種GG型個體的序列與NCBI中U70 883.2的序列一致,這與酶切結果一致,沒有發(fā)現(xiàn)突變的基因位點。

    3個豬種NRAMP1 基因第6內(nèi)含子的PCR產(chǎn)物經(jīng)限制性內(nèi)切酶NdeⅠ酶切,經(jīng)過1.5%的瓊脂糖凝膠電泳后得到3個豬種的NRAMP1 基因都只有1種基因型條帶,為BB表型(373/110 bp)(圖4)。根據(jù)酶切結果對其進行驗證,對PCR產(chǎn)物進行測序發(fā)現(xiàn),在測序峰圖中所有豬種個體的NRAMP1基因都不存在單堿基突變,將該位點未突變呈現(xiàn)單峰的基因型定義為BB型(圖4D)。所有豬種個體的序列與NCBI中上AY368 473.1的序列一致,這與酶切結果一致,沒有發(fā)現(xiàn)突變的基因位點(圖4E)。

    對上海白豬、楓涇豬的MUC4 基因擴增得到的所有PCR產(chǎn)物進行測序分析,對所得測序峰圖(圖5)分析發(fā)現(xiàn),部分豬種的MUC4 基因有1個G→A的單堿基突變,將該位點未突變呈現(xiàn)單峰的基因型定義為GG型,將該位點突變呈現(xiàn)單峰的基因型定義為AA型,呈現(xiàn)雙峰的基因型定義為AG型(圖5B)。對沙烏頭豬的MUC4 基因的PCR產(chǎn)物經(jīng)限制性內(nèi)切酶HhaⅠ酶切后,經(jīng)過1.5%的瓊脂糖凝膠電泳后都得到了3種不同的基因型條帶,分別命名為AA、GG和AG,抗性等位基因(A)不能被HhaI 消化,為AA型(367 bp),而敏感等位基因被消化成2個片段,為GG型(151/216 bp)。2個基因型同時存在時為AG型(151/216/367 bp),如圖5所示。根據(jù)酶切結果,對PCR產(chǎn)物進行測序發(fā)現(xiàn),在測序峰圖中沙烏頭豬部分豬種的MUC4 基因有1 個G→A 的單堿基突變。將上述豬種的MUC4 基因所測序列與NCBI中KT966 749的序列對比,AA型個體的序列與其一致,定義為野生型,相比之下發(fā)現(xiàn)GG型有1個G→A的單堿基突變,這與酶切結果一致。

    對上述豬種的MUC13 基因擴增得到的所有PCR產(chǎn)物進行測序分析,得到測序峰圖與基因序列,對所得的測序峰圖進行分析發(fā)現(xiàn),上海白豬中有部分豬種存在G→A的單堿基突變,楓涇豬和沙烏頭豬的所有個體都存在G→A 的單堿基突變,將該位點未突變呈現(xiàn)單峰的基因型定義為AA型,將該位點突變呈現(xiàn)單峰的基因型定義為GG型,呈現(xiàn)雙峰的基因型定義為AG型(圖6A)。得到的序列與NCBI上JN613 418.1對比,所有個體的對比結果與測序峰圖一致(圖6B)。

    2.3 TAP1、FUT1、NRAMP1、MUC4 和MUC13基因型和等位基因頻率分析

    利用PCR-RFLP方法對3個豬群個體進行檢驗,分型后計算各豬群的基因型和等位基因頻率,結果如表2所示。上海白豬中TAP1 基因有1個等位基因,2種基因型(AG和GG),其中AG、GG基因型數(shù)量和基因型頻率分別為17、13和0.57、0.43;抗病等位基因G的頻率為0.72,敏感等位基因A的頻率為0.28。楓涇豬TAP1 基因的AA、AG、GG基因型數(shù)量及基因型頻率分別為3、11、16 和0.10、0.37、0.53,抗病等位基因G和敏感等位基因A的頻率分別為0.72和0.28。在檢測的30個沙烏頭豬樣本中,TAP1 基因的AA、AG、GG 基因型數(shù)量分別為1、19、10,基因型頻率分別為0.03、0.63、0.33,抗病等位基因G的頻率為0.65,敏感等位基因A的頻率為0.35。在3個豬群中,AA型的頻率不高,總體上等位基因G為優(yōu)勢等位基因。

    在檢測的所有豬種的樣本中,F(xiàn)UT1 基因只檢測出敏感基因型GG型,等位基因G為優(yōu)勢等位基因。

    對NRAMP1 基因的檢測表明,3 個地方豬種中沒有出現(xiàn)抗病基因型AA型,只檢測出敏感基因型BB型,等位基因B為優(yōu)勢等位基因。

    如表2所示,共檢測出MUC4 基因的1對等位基因,3個豬種均檢測出3種基因型,但在各個豬種中所含有的基因型和基因頻率有差異;在上海白豬中,MUC4 基因的AA、GG、AG基因型數(shù)量及頻率分別為1、25、4 和0.03、0.83、0.13,其中G 為優(yōu)勢等位基因。在楓涇豬中MUC4 基因的AA、GG、AG 基因型數(shù)量及頻率分別為10、1、19 和0.33、0.03、0.63,其中A為優(yōu)勢等位基因。在沙烏頭豬中MUC4 基因的AA、GG、AG 基因型數(shù)量及頻率分別為6、10、14和0.20、0.33、0.47,其中G為優(yōu)勢等位基因。

    3個地方豬種中MUC13 基因共檢測出1對等位基因,檢測出AA、GG和AG共3種基因型,上海白豬中檢測出2種基因型(AA和AG),楓涇豬和沙烏頭豬只檢測出了1種基因型(GG),在上海白豬中AA和AG基因型頻率分別為0.87和0.13,其中A為優(yōu)勢等位基因。

    2.4 TAP1、FUT1、NRAMP1、MUC4 和MUC13基因的遺傳多態(tài)性分析

    由表3 可知,上海白豬、楓涇豬和沙烏頭豬TAP1 基因的多態(tài)信息含量(PIC)處于0.25~0.50,都屬于中度多態(tài)。沙烏頭豬的Ho、He和PIC值均比上海白豬和楓涇豬高,表明其遺傳多樣性高于其他2個品種。上海白豬和沙烏頭豬的χ20.05(1)<χ2<χ20.01(1)(χ20.05(1)=3.84,χ20.01(1)=6.63),0.010.05,差異不顯著,滿足Hardy-Weinberg平衡。

    對上海白豬、楓涇豬和沙烏頭豬FUT1 和NRAMP1 基因的酶切與測序結果顯示,只有GG和BB 基因型,無法進行遺傳多態(tài)性和Hardy-Weinberg平衡分析。

    上海白豬的MUC4 基因的PIC<0.25,屬于低度多態(tài),楓涇豬和沙烏頭豬的MUC4 基因PIC 處于0.25~0.50,屬于中度多態(tài)。沙烏頭豬的He 與PIC值均高于楓涇豬和沙烏頭豬,表明沙烏頭豬的遺傳多態(tài)性高于其他2個品種,經(jīng)過比對,3個品種中上海白豬遺傳多態(tài)性最低。分析發(fā)現(xiàn),上海白豬和沙烏頭豬χ2<χ20.05(1),P>0.05,差異不顯著,滿足Hardy-Weinberg 平衡。楓涇豬χ20.05(1)<χ2<χ20.01(1),0.01

    上海白豬MUC13 基因PIC<0.25,屬于低度多態(tài)。χ2<χ20.05(1),P>0.05,差異不顯著,滿足Hardy-Weinberg平衡。然而,楓涇豬和沙烏頭豬MUC13基因僅有GG基因型,無法繼續(xù)進行遺傳多態(tài)性和Hardy-Weinberg平衡分析。

    2.5 上海白豬、楓涇豬和沙烏頭豬免疫因子的測定分析

    不同豬種免疫因子測定結果見表4,上海白豬、楓涇豬和沙烏頭豬之間血液免疫因子平均水平無顯著差異(P>0.05)。分析發(fā)現(xiàn),沙烏頭豬的IL-12和TNF-α免疫因子均比上海白豬和楓涇豬高。其中,IL-6免疫因子水平與其正好相反,上海白豬血液中IL-6水平高于楓涇豬和沙烏頭豬。楓涇豬的IFN-γ、IL-10和IgG 的水平比上海白豬和楓涇豬高,其中IL-10和IgG的水平較低的是沙烏頭豬。上海白豬IL-2的水平較高,楓涇豬IL-2的水平低于其他2個豬種。

    2.6 TAP1、FUT1、NRAMP1、MUC4 和MUC13基因多態(tài)性與免疫因子相關性分析

    由于FUT1 和NRAMP1 基因只檢測出1種基因型,所以這2個基因不與免疫因子進行相關性分析。如表5所示,上海白豬的TAP1 基因各基因型個體表達的免疫因子的量差異不顯著(P>0.05),但GG基因型的IFN-γ、IL-12、IgG、TNF-α免疫因子水平高于AG基因型,AG基因型的IL-2、IL-6和IL-10的免疫因子水平高于GG基因型。在上海白豬中,MUC4 基因由于AA 型只檢測出1個,AA型免疫因子與GG和AG不作比較,通過最小二乘法分析,GG和AG型各免疫因子不存在顯著差異(P>0.05)。MUC13 基因由于沒有檢測到GG 基因型,GG 型免疫因子與AA 和AG 不作比較,通過最小二乘法分析,AA型的IL-12水平顯著高于AG型(P<0.05)。其他免疫因子AA和AG型之間不存在顯著差異(P>0.05)。

    對楓涇豬TAP1 基因的各基因型的免疫因子的分析可知,TNF-α在3個基因型豬中存在明顯差異,其中AA 型的TNF-α 水平顯著高于GG 和AG型(P<0.05)。其他免疫因子在各基因型間不存在差異,但從結果看,AA型中IL-2的水平高于GG和AG型,其中AG型IL-2水平較低。GG型中IFN-γ和IL-10水平高于AA和AG型,其中AA型的IFN- γ 和IL-10 水平較低。IL-6 水平AG 型較高,GG 型較低。AG 型的IL-12 水平高于AA 和GG,AA型的IL-12水平較低。還可以看出,GG的IgG 水平高于AA 和AG,其中AG 型水平相對較低。楓涇豬MUC4 基因由于GG型只檢測出1個,GG型免疫因子與AA和AG不作比較,通過最小二乘法分析,GG 和AG 型免疫因子不存在差異(P>0.05)。

    在沙烏頭豬中,TAP1 基因由于AA型只檢測出1個,AA型免疫因子與GG和AG不作比較,通過最小二乘法分析,GG 型的IL-12水平與AG 型之間存在顯著差異(P<0.05),其他免疫因子的GG和AG型之間不存在顯著差異。從表4可知,GG型的IFN-γ、IL-2和IgG高于AG型。AG型的IL-6和IL-12 高于GG 型。GG 型與AG 型的TNF-α 水平趨于一致。沙烏頭豬中MUC4 各基因型的免疫因子沒有顯著差異(P>0.05)。

    3 討論

    3.1 抗腹瀉基因的多態(tài)性

    本研究對上海地區(qū)上海白豬、楓涇豬和沙烏頭豬3個地方豬種的TAP1 基因外顯子3多態(tài)性分析發(fā)現(xiàn),上海白豬檢測到AG和GG型,楓涇豬和沙烏頭豬都檢測到AA、AG和GG 3種基因型,其中G為優(yōu)勢等位基因。趙喬輝等[6]研究發(fā)現(xiàn),中國本地豬種與亞洲野豬及國外引進豬品種的基因頻率存在明顯差異,外來引進豬的等位基因G頻率相對較高,而中國地方豬種(淮豬和榮昌豬等)的A頻率相對較高,這與本研究結果一致,本研究中上海白豬沒有檢測出AA型基因,可能是因為上海白豬是與外國豬種培育的結果。Zhao 等[18]報道,蘇太豬TAP1 基因G729A突變顯著影響mRNA表達,且GG基因型個體抗腹瀉能力顯著高于其他2種基因型的個體。因此,TAP1 基因的GG基因型可作為抗腹瀉感染的候選基因型。

    本研究中,F(xiàn)UT1 基因M307位點結果顯示,3個地方豬種均僅存在1 種不抗病基因型(GG)。張雄等[19]發(fā)現(xiàn),國內(nèi)地方品種只含有GG型,而國外及國內(nèi)雜交培育品種豬FUT1 基因M307位點中含有AA、AG 和GG 基因型,其中GG 為絕對優(yōu)勢基因型。推測FUT1 基因抗病性等位基因A在中國地方豬種中呈極端偏態(tài)分布,AA抗病基因型可能來源于外國豬種。這一結果也與Yan等[20]和吳華莉等[21]的研究結果一致,可能是上海的幾個地方豬種與國內(nèi)外豬種基因交流較少,具有上海本地獨特的遺傳資源性。也有學者在中國上海的浦東白豬中檢測出了抗病基因GG型[22],可能是由于浦東白豬作為本土種豬和國外豬種雜交,使AA抗病基因型得到遺傳。

    作為相對比較保守的基因,NRAMP1 基因會對動物自身免疫產(chǎn)生影響,且對沙門氏菌及多種胞內(nèi)寄生病原菌具有抗病[23]。本研究利用PCRRFLP技術,在3個豬種中均未檢測到AA型,只檢測到BB 型。楊軍等[24] 在多個豬種中檢測到NRAMP1 基因的AB和BB型,均未檢測到抗病基因型AA型。但劉艷冬[25]研究發(fā)現(xiàn),貴州地方豬種NRAMP1 基因內(nèi)含子6具有豐富多態(tài)性,其中香豬抗病基因型AA型比例為26%、糯米豬為20%、宗地花豬為7.5%、大約克豬為4.3%。其他學者在中國的其他地方豬種中也檢測到有AA型[22]??赡苁巧虾0棕i、楓涇豬和沙烏頭豬在長期的進化選擇中,該基因的個體逐漸被淘汰,群體該位點的遺傳平衡從而被打破。

    研究表明,MUC4 基因可以作為F4ab/F4ac抗性的候選基因,且G為控制F4ac大腸菌抗性等位基因[2627]。本研究對上海地區(qū)3 個地方豬種的MUC4 基因第17 內(nèi)含子243 位點多態(tài)性鑒定發(fā)現(xiàn),在該位點共檢測到了1 個等位基因及AA、AG、GG共3種基因型,其中上海白豬、楓涇豬和沙烏頭豬中都檢測到3種基因型。在上海白豬和沙烏頭豬中G為優(yōu)勢等位基因,楓涇豬中A為優(yōu)勢等位基因。李聰聰?shù)萚28]研究表明,淮南豬和太白豬的MUC4 基因檢測出3 種基因型AA、AG 和GG,這與本研究結果相似,在Liu 等[27]和李聰聰?shù)萚28]研究中,GG型個體對抗ETEC F4ab/F4ac感染能力高于其他基因型。本研究中,3個地方豬種都存在GG型抗病基因,可以通過選育來實現(xiàn)群體抗ETEC F4ab/F4ac感染。

    MUC13 基因被定位于豬13號染色體長臂41區(qū)的Sw207和S0075之間,研究表明,MUC13 基因是決定仔豬F4ac腹瀉易感/抗性的基因之一[29-31]。本研究對上海地區(qū)3個地方豬種的MUC13 基因第2內(nèi)含子多態(tài)性研究發(fā)現(xiàn),在該位點共檢測到了1個等位基因和AG、GG 2種基因型。其中上海白豬檢測到2種基因型,楓涇豬和沙烏頭豬只檢測出抗病基因型GG型,其中G為優(yōu)勢等位基因。與本研究的結果相似,楊明等[32]研究杜洛克豬和皮特豬中MUC13 基因的分子標記輔助選擇發(fā)現(xiàn),2個豬種中G為優(yōu)勢等位基因;任軍等[13]發(fā)現(xiàn),該基因在rs319699771 位點的GG 型為優(yōu)勢抗腹瀉基因型,而AG型和AA型為易感腹瀉基因型。本研究發(fā)現(xiàn),楓涇豬和沙烏頭豬只存在抗病基因型GG型,可能是這2個豬種在長期的進化選擇中,該基因的個體逐漸被保留,形成了抗病基因型的種群。

    3.2 遺傳多態(tài)性分析

    作為群體內(nèi)遺傳變異的檢測參數(shù),PIC、Ho和He 數(shù)值反映了群體內(nèi)個體的均質(zhì)度,其數(shù)值越高,說明遺傳變異性就越大,對環(huán)境的適應能力越強,具有較大的選擇潛力,應用于遺傳育種效果越好,反之則說明該群體雜合度小、基因純合度高,可以看作純系加以利用[33-35]。上海白豬、楓涇豬和沙烏頭豬TAP1 基因的PIC處于0.25~0.50,都屬于中度多態(tài)。沙烏頭豬的TAP1 基因Ho和He值高于上海白豬和楓涇豬,說明沙烏頭豬遺傳多樣性豐富,對環(huán)境的適應能力較強,遺傳潛力大。經(jīng)過分析,TAP1 基因在上海白豬和沙烏頭豬中不滿足Hardy-Weinberg平衡,這一結果與在浦東白豬、皖南黑豬、霍壽黑豬和梅山豬等地方豬種中的研究結果一致,TAP1 基因偏離了Hardy-Weinberg 平衡[6, 22]。楓涇豬中滿足Hardy-Weinberg平衡,說明TAP1 基因在楓涇豬中基本處于自由交配的狀態(tài),也就是群體的選育對TAP1 基因沒有造成影響,對現(xiàn)在的選育性狀沒有作用。這與前人在杜洛克、大白和長白中的研究結果一致[6]。

    上海白豬的MUC4基因?qū)儆诘投榷鄳B(tài)(<0.25),而楓涇豬和沙烏頭豬的MUC4 基因?qū)儆谥卸榷鄳B(tài)。沙烏頭豬He值比其他2個豬種高,與PIC值的變化趨勢一致,說明沙烏頭豬的遺傳多態(tài)性較高,可作為抗病育種的標記輔助性選擇。上海白豬MUC4 基因的PIC較低,說明遺傳變異性較低,但分析發(fā)現(xiàn)上海白豬和沙烏頭豬滿足Hardy-Weinberg平衡,楓涇豬不滿足Hardy-Weinberg平衡。說明人工選擇影響了MUC4 基因抗腹瀉基因型在楓涇豬中的分布。上海白豬MUC13 基因?qū)儆诘投榷鄳B(tài),滿足Hardy-Weinberg平衡。此外,本試驗中有部分豬種基因位點偏離了Hardy-Weinberg平衡??赡苁怯捎谌后w研究的樣本量不夠或者群體內(nèi)近交積累過高[36]。本研究的上海白豬、楓涇豬和沙烏頭豬品系樣本量都達到了統(tǒng)計學意義,但分析3個品種的基因多態(tài)性以及雜合度,結果表明,群體中出現(xiàn)近交積累過高和部分等位基因缺失的問題。進一步分析群體出現(xiàn)連鎖不平衡,主要是由于目前上海白豬、楓涇豬和沙烏頭豬群體在保種過程中出現(xiàn)了高度近交或遺傳漂變,而這與上述豬種后期的品種改良有關。

    3.3 抗腹瀉基因多態(tài)性與免疫因子相關性分析

    IgG是在針對外來抗原的免疫應答過程中產(chǎn)生的主要免疫球蛋白類別,通過其雙功能特性有效地提供保護,血清中IgG含量可以代表機體體液的免疫狀態(tài)[37]。IL-2、IL-6、IL-10 和IL-12 在機體免疫系統(tǒng)中具有重要作用,其中IL-2主要是激活免疫,IL-6 和IL-10 主要參與炎癥反應[38-40]。IL-12主要是促進T細胞和NK細胞的增殖,也被證明可以促進Thl細胞的分化和增殖,阻止瘧疾的傳播[4142]。作為一種強烈抗病毒、抗腫瘤和免疫調(diào)節(jié)作用的細胞因子,IFN-γ在機體免疫系統(tǒng)中具有重要作用[43]。TNF-α屬于腫瘤壞死因子超家族的一員,在機體內(nèi)可抑制病毒復制、病原體的擴散,并預防腫瘤發(fā)生,是重要的抗病因子[44]。上述血清學免疫因子,都是機體免疫不可缺少,它們之間互相影響、互相調(diào)節(jié),存在密切聯(lián)系。本研究對3個豬種的上述免疫因子分析顯示,各豬種間的免疫因子水平均差異不顯著,但在沙烏頭豬中IL-6指標低于其他2個豬種。研究表明,低水平IL-6更具有抗病優(yōu)勢[45],說明沙烏頭豬的抗病性可能更高。此外,對上海地區(qū)3個地方豬種的部分抗腹瀉基因的基因型與部分免疫因子進行關聯(lián)性分析,在上海白豬中僅有MUC13 基因AA型的IL-12免疫因子顯著高于AG型;楓涇豬中TAP1 基因的AA基因型個體TNF-a表達水平顯著高于GG和AG型,這與薛偉偉[15]對圩豬的研究結果一致。由于TNF-a不僅能夠調(diào)控促炎細胞因子,也能控制趨化因子和黏附因子的活化和表達[46]。因此,推測在正常值范圍內(nèi),AA型豬的抗病能力與AG利GG型存在差異。在沙烏頭豬中,TAP1 基因GG型的IL-12水平顯著低于AG型,研究表明,IL-12水平過高會抑制Th2細胞的產(chǎn)生從而抑制IL-10的產(chǎn)生,且IL-10的表達與IL-12的表達早現(xiàn)負相關[47]。雖然沙烏頭豬的TAP1 基因GG型與AG型的IL-10 沒有顯著差異,但AG 型的IL-10 水平低于GG型。本研究總結了上海白豬,楓涇豬和沙烏頭豬部分抗腹瀉基因與部分免疫因子間的的關系,通過分析,對各豬種的抗腹瀉有了初步的認識,可為后期試驗提供一定的理論依據(jù)。

    參 考 文 獻

    [1] CANIBE N, H?JBERG O, KONGSTED H, et al .. Review onpreventive measures to reduce post-weaning diarrhoea inpiglets [J/OL]. Animals, 2022, 12(19):2585 [2023-02-23]. https://doi.org/10.3390/ani12192585.

    [2] MEIJERINK E, NEUENSCHWANDER S, FRIES R, et al .. ADNA polymorphism influencing α(1,2) fucosyltransferase activityof the pig FUT1 enzyme determines susceptibility of smallintestinal epithelium to Escherichia coli F18 adhesion [J].Immunogenetics, 2000, 52:129-136.

    [3] 方晨, 胡瑞舉, 楊明華, 等. 迪慶藏豬FUT1 基因遺傳變異分析[J]. 中國畜牧雜志, 2020, 56(12):29-35.

    FANG C, HU R J, YANG M H, et al .. Genetic variationanalysis of FUT1 gene in Diqing Tibetan pig [J]. Chin. J. Anim.Sci., 2020, 56 (12):29-35.

    [4] MEIJERINK E, FRIES R, VOGELI P, et al .. Two α(1, 2)fucosyltransferase genes on porcine chromosome 6q11 areclosely linked to the blood group inhibitor (S) and Escherichiacoli F18 receptor (ECF18R) loci [J]. Mamm. Genome, 1997, 8(10):736-741.

    [5] KIM K, NGUYEN D T, CHOI M, et al .. Alpha (1, 2) -fucosyltransferase M307A polymorphism improves pigletsurvival [J]. Anim. Biotechnol., 2013, 24(3):243-250.

    [6] 趙喬輝, 劉穎, 董文華, 等. TAP1 基因外顯子3在9個不同豬群體中的遺傳多態(tài)性分析[J]. 中國畜牧雜志, 2014, 50(7):1-5.

    ZHAO Q H, LIU Y, DONG W H, et al .. Analysis of geneticdiversity in TAP1 gene exon 3 from 9 different pig populations [J].Chin. J. Anim. Sci., 2014, 50(7):1-5.

    [7] 劉艷冬, 許厚強, 稽辛勤, 等. 香豬NRAMP1 基因多態(tài)性與仔豬腹瀉的研究[J]. 中國畜牧獸醫(yī), 2010, 37(4):125-127.

    LIU Y D, XU H Q, JI X Q, et al .. NRAMP1 gene polymorphismand pig piglet diarrhea research [J]. China Anim. Husb. Vet.Med., 2010, 37(4):125-127.

    [8] 劉艷冬, 許厚強, 林家棟, 等. 貴州地方豬NRAMP1 基因多態(tài)性研究[J]. 貴州畜牧獸醫(yī), 2009, 33(2):1-2.

    LIU Y D, XU H Q, LIN J D, et al .. Genetic polymophism ofNRAMP1 genc in Guizhou local pig [J]. Guizhou J. Anim.Husb. Vet. Med., 2009, 33(2):1-2.

    [9] JOLLER D, JORGENSEN C B, BERTSCHINGER H U, et al ..Refined localization of the Escherichia coli F4ab/F4ac receptorlocus on pig chromosome 13 [J]. Anim. Genet., 2009, 40(5):749-752.

    [10] JACOBSEN M, KRACHT S S, ESTESO G, et al .. Refinedcandidate region specified by haplotype sharing for Escherichiacoli F4ab/F4ac susceptibility alleles in pig [J]. Anim. Genet.,2010, 41(1):21-25.

    [11] 洪淵, 韓云珍, 薛永欽, 等. 長白豬MUC13基因多態(tài)性及其與生產(chǎn)性狀的關聯(lián)性分析[J]. 福建畜牧獸醫(yī), 2022, 44(4):7-9.

    HONG Y, HAN Y Z, XUE Y Q, et al .. Association betweenpolymorphism of MUC13 gene and production traits inAmerican Landrace pig [J]. Fujian J. Anim. Husb. Vet. Med.,2022, 44(4):7-9.

    [12] PENG Q L, REN J, YAN X M, et al .. The g. 243A> G mutationin intron 17 of MUC4 is significantly associated withsusceptibility/resistance to ETEC F4ab/ac infection in pig [J].Anim. Genet., 2007, 38(4):397-400.

    [13] 任軍, 晏學明, 艾華水, 等. 仔豬斷奶前腹瀉抗病基因育種技術的創(chuàng)建及應用[J]. 豬業(yè)科學, 2012, 29(1):44-48.

    REN J, YAN X M, AI H S, et al .. Establishment andapplication of gene breeding technology for resistance todiarrhea in piglets before weaning [J]. Swin. Ind. Sci., 2012, 29 (1):44-48.

    [14] 劉榮榮,徐勇,周慧姝,等.規(guī)?;i場免疫程序的制定標準[J]. 畜牧獸醫(yī)科學, 2017 (10):87.

    [15] 薛瑋瑋. 豬TAP1 基因多態(tài)性及其與斷奶仔豬免疫指標的相關性分析[D]. 合肥:安徽農(nóng)業(yè)大學, 2014.

    XUE W W. Poylmorphisms of porcine TAP1 gene and theirassociationswith immune indicator in weaned piglets [D]. Hefei:Anhui Agricultural University, 2014.

    [16] 劉穎, 董文華, 蘇先敏,等. 梅山豬Mx1 基因第14外顯子多態(tài)性及其與部分免疫指標和早期生長性能的關聯(lián)分析[J].中國畜牧雜志, 2014, 50(5):18-20,64.

    LIU Y, DONG W H, SU X M, et al .. Polymorphisms in exon 14of Mx1 gene and its relationship with partial cytokines leveland early growth performance in Meishan pig [J]. Chin. J.Anim. Sci., 2014, 50(5):18-20,64.

    [17] 馬小軍. 不同豬種免疫功能及免疫相關基因表達差異研究[D]. 蘭州:甘肅農(nóng)業(yè)大學, 2005.

    MA X J. Studies on immune function and differentialexpressionof immune-related gene in different pig breeds [D].Lanzhou: Gansu Agricultural University, 2005.

    [18] ZHAO Q, LIU Y, DONG W, et al .. Genetic variations of TAP1gene exon 3 affects gene expression and Escherichia coli F18resistance in piglets [J]. Int. J. Mol. Sci., 2014, 15(6):11161-11171.

    [19] 張雄, 王婧, 張靜, 等. 3個群體豬FUT1 基因M229與M307位點遺傳多樣性研究[J]. 基因組學與應用生物學, 2020, 39(6):2536-2541.

    ZHANG X, WANG J, ZHANG J, et al .. The polymorphismanalysis between M229 and M307 locus of FUT1 gene in threegroup pig [J]. Genomics Appl. Biol., 2020, 39(6):2536-2541.

    [20] YAN X M, GUO Y M, DING N S, et al .. Study on the geneticvariation ofal-fucosytransferase gene in different pig breed [J].Chin. J. Anim. Sci., 2004, 40(3): 8-10.

    [21] 吳華莉, 涂尾龍, 曹建國, 等. FUT1 基因在上海地區(qū)5個豬群體內(nèi)的多態(tài)性分析[J]. 畜牧與飼料科學, 2017, 38(1):18-20.

    WU H L, TU W L, CAO J G, et al .. Polymorphism analysis ofFUT1 gene in five pig herds in Shanghai [J]. Anim. Husb. FeedSci., 2017, 38(1):18-20

    [22] 張艷. 浦東白豬重要抗病基因多態(tài)性位點及其可能相關免疫學指標檢測分析[D]. 上海:上海交通大學, 2016.

    ZHANG Y. Polymorphism sites of important disease resistantgenes in Pudong white pig and detection and analysis ofimmune indices possibly related to the genes [D]. Shanghai:Shanghai Jiaotong University, 2016.

    [23] BLACKWELL J M. Structure and function of the naturalresistance-associated macrophage protein (NRAMP1), a candidateprotein for infectious and autoimmune disease susceptibility [J].Mol. Med. Today, 1996, 2(5):205-211.

    [24] 楊軍, 陳紅頌, 丁山河,等. 3個外種豬NRAMP1 基因多態(tài)性研究[J]. 長江大學學報(自然科學版), 2009, 6(4):40-43,120.

    YANG J, CHEN H S, DING S H, et al .. Study on polymorphismof NRAMP1 gene in three exotic pig [J]. J. Yangtze Univ. (Nat.Sci.), 2009, 6(4):40-43,120.

    [25] 劉艷冬. 香豬NRAMP1 基因多態(tài)性及與仔豬腹瀉相關性研究[D]. 貴陽:貴州大學, 2009.

    LIU Y D. Study of Xiang pig NRAMP1 gene polymorphism andcorrelation with the piglet diarrhea [D]. Guiyang: GuizhouUniversity, 2009.

    [26] 欽偉云, 甘麗娜, 魏宗友, 等. MUC4 和MUC13 基因在ETECF18抗性型和易感型梅山斷奶仔豬間的差異表達分析[J]. 中國畜牧獸醫(yī), 2017, 44(4):959-964.

    QIN W Y, GAN L N, WEI Z Y, et al .. Differential expressionanalysis of MUC4 and MUC13 genes between resistant andsensitive weaned Meishan piglets to ETEC F18 [J]. ChinaAnim. Husb. Vet. Med., 2017, 44(4):959-964.

    [27] LIU Y, YIN X M, XIA R W, et al .. Association between theMUC4 g. 243A> G polymorphism and immune andproductiontraits in Large white pig [J]. Turk. J. Vet. Anim. Sci., 2015, 39(2):141-146.

    [28] 李聰聰, 劉思雨, 吳姣, 等. 不同品種豬SLA-DQA、FUT1 和MUC4 基因遺傳變異初析[J]. 黑龍江畜牧獸醫(yī), 2020 (2):55-61.

    LI C C, LIU S Y, WU J, et al .. Preliminary analysis of geneticvariation of SLA-DQA, FUT1 and MUC4 genes in differentbreeds of pig [J]. Heilongjiang J. Anim. sci., 2020 (2):55-61.

    [29] JACOBSEN M, CIRERA S, JOLLER D, et al .. Characterisationof five candidate genes within the ETEC F4ab/ac candidateregion in pig [J/OL]. BMC Res. Notes, 2011, 4: 225 [2023-02-23]. https://doi.org/10.1186/1756-0500-4-225.

    [30] JACOBSEN M, KRACHT S S, ESTESO G, et al .. Refinedcandidate region specified by haplotype sharing for Escherichiacoli F4ab/F4ac susceptibility alleles in pig [J]. Anim. Genet.,2010, 41(1):21-25.

    [31] REN J, YAN X, AI H, et al .. Susceptibility towardsenterotoxigenic Escherichia coli F4ac diarrhea is governed bythe MUC13 gene in pig [J/OL]. PLoS One, 2012, 7(9):e44573[2023-02-23]. https://doi.org/10.1371/journal.pone.0044573.

    [32] 楊明, 王青來, 劉敬順, 等. MUC13、FUT1 基因在2個種豬核心群中的分子標記輔助選擇研究[J]. 華南農(nóng)業(yè)大學學報,2015, 36(6):1-8.

    YANG M, WANG Q L, LIU J S, et al .. Molecular maker-assistedselections of MUC13 and FUT1 genes in the two swine nucleuspopulations [J]. J. South China Agric. Univ., 2015, 36(6):1-8.

    [33] 蔣欽楊, 韋英明, 陳寶劍, 等. 馬生長激素基因多態(tài)性與體尺指標之間的關聯(lián)性分析[J]. 中國畜牧雜志, 2013, 49(3):1-4.

    JIANG Q Y, WEI Y M, CHEN B J, et al .. Association ofpolymorphisms of growth hormone gene with body sizesindexes in horses [J]. Chin. J. Anim. Sci., 2013, 49(3):1-4.

    [34] 嚴燕, 張陳華, 王陽, 等. 圩豬OB 基因SNPs檢測及其與產(chǎn)仔性能的關系[J]. 中國農(nóng)業(yè)大學學報, 2010, 15(6):78-83.

    YAN Y, ZHANG C H, WANG Y, et al .. SNPs detection of OBgene and its association with reproduction in Wei pig [J]. J.China Agric. Univ., 2010, 15(6):78-83

    [35] 相德才, 張斌, 劉邵娜, 等. 迪慶藏豬NRAMP1 基因第二內(nèi)含子和第二外顯子多態(tài)性位點遺傳分析[J]. 養(yǎng)豬, 2020, 169(2):60-62.

    XIANG D C, ZHANG B, LIU S N, et al .. Genetic analysis ofthe polymorphism in the second intron and the second exon ofNRAMP1 gene in diqing tibetan pig [J]. Swin. Prod., 2020, 169(2):60-62.

    [36] ZERGER K R, RICHARDERSON B J. A rapid poulationexpansion retains genetic diversity within european rabbit inaustralia [J]. Mol. Ecol., 2003, 12(3):789-794.

    [37] BOURNAZOS S, RAVETCH J V. Diversification of IgGeffector functions [J]. Int. Immunol., 2017, 29(7):303-310.

    [38] KOGUT M, ROTHWELL L, KAISER P. Differential effects ofage on chicken heterophil functional activation by recombinantchicken interleukin-2 [J]. Dev. Comp. Immunol., 2002, 26(9):817-830.

    [39] ZIZZO G, TAMBURELLO A, CASTELNOYO L, et al ..Immunotherapy of COVID-19: inside and beyond IL-6signalling [J/OL]. Front. Immunol., 2022, 13:795315 [2023-02-23]. https://doi.org/10.3389/fimmu.2022.795315.

    [40] WU Y R, HSING C H, CHIU C J, et al .. Roles of IL-1 and IL-10 family cytokines in the progression of systemic lupuserythematosus: friends or foes [J]. IUBMB Life, 2022, 74(2):143-156.

    [41] XUE D, MOON B, LIAO J, et al .. A tumor-specific pro-IL-12activates preexisting cytotoxic T cells to control establishedtumors [J/OL]. Sci. Immunol., 2022, 7(67):6899 [2023-02-23].https://www.science.org/doi/10.1126/sciimmunol.abi6899.

    [42] 姜素華, 唐慧, 邢建新. IL-12在約氏瘧原蟲感染中作用的初步實驗研究[J]. 農(nóng)墾醫(yī)學, 2008, 30(5):360-362.

    JIANG S H, TANG H, XING J X. Primary research oninfluence of IL-12 on plasmodium yoelii infection [J]. J.NongkenMed., 2008, 30(5):360-362.

    [43] 郭紫晶, 陳飛, 張志雄, 等. 白細胞介素-10對口蹄疫病毒感染小鼠T 細胞增殖及其表達TNF- α、IFN- γ 和IL-2 的影響[J]. 畜牧獸醫(yī)學報, 2023, 54(2):694-705.

    GUO Z J, CHEN F, ZHANG Z X, et al .. Effects of interleukin-10 on T cell proliferation and its expression of TNF-a, IFN- γand IL-2 in mice infected with foot and-mouth disease virus [J].J. Anim. Vet. Sci., 2023, 54 (2):694-705.

    [44] KOPER-LENKIEWICZ O M, SUTKOWSKA K,WAWRUSIEWICZ-KURYLONEK N, et al .. Proinflammatorycytokines (IL-1,-6,-8,-15,-17,-18,-23, TNF-α) single nucleotidepolymorphisms in rheumatoid arthritis-a literature review [J/OL].Int. J. Mol. Sci., 2022, 23(4):2106 [2023-02-23]. https://doi.org/10.3390/ijms23042106.

    [45] 呂瓊霞, 張書霞, 趙茹茜. 運輸應激對豬淋巴結IL-2、IL-6和IL-10 及其受體mRNA 轉(zhuǎn)錄的影響[J]. 南京農(nóng)業(yè)大學學報, 2011, 34(1):95-100.

    LYU Q X, ZHANG S X, ZHAO R X. Effects of transportationstress on transcription of IL-2, IL-6 and IL-10 and the relativereceptors mRNA of lymph nodes in pig [J]. J. Nanjing Agric.Univ., 2011, 34 (1):95-100.

    [46] PENNICA D, KOHR W J, FENDLY B M, et al ..Characterization of a recombinant extracellular domain of thetype 1 tumor necrosis factor receptor: evidence for tumornecrosis factor. alpha. induced receptor aggregation [J].Biochemistry, 1992, 31(4):1134-1141.

    [47] OGARRA A, MURPHY K M. From IL-10 to IL-12: howpathogens and their products stimulate APCs to induce TH1development [J]. Nat. Immunol., 2009, 10(9):929-932.

    (責任編輯:胡立霞)

    基金項目:國家重點研發(fā)計劃項目(2021YFD120303);上海市農(nóng)業(yè)科學院助跑計劃人才項目(ZP22171);重慶市技術創(chuàng)新與應用發(fā)展專項(CSTC2021JSCX-DXWTBX0004);寧波科技創(chuàng)新2025重大專項(2019B10023)。

    黑丝袜美女国产一区| 午夜日韩欧美国产| 欧美日韩成人在线一区二区| 亚洲一码二码三码区别大吗| 欧美精品一区二区大全| 久久人妻熟女aⅴ| 国产高清不卡午夜福利| 欧美激情极品国产一区二区三区| 亚洲成av片中文字幕在线观看| 深夜精品福利| 九草在线视频观看| 久久女婷五月综合色啪小说| 制服丝袜香蕉在线| 99热网站在线观看| 日本欧美国产在线视频| 亚洲伊人色综图| 在线天堂中文资源库| 一区二区三区精品91| 激情五月婷婷亚洲| 一二三四在线观看免费中文在| 精品一区二区三卡| 51午夜福利影视在线观看| 亚洲国产av影院在线观看| 久久99热这里只频精品6学生| 日韩一区二区三区影片| 国产精品.久久久| 黄片无遮挡物在线观看| 丰满迷人的少妇在线观看| 久久久久精品性色| 久久天堂一区二区三区四区| 18在线观看网站| 一级片'在线观看视频| 亚洲国产欧美网| 另类精品久久| 女性生殖器流出的白浆| 欧美精品一区二区免费开放| 我的亚洲天堂| 久久久久久人妻| 精品少妇黑人巨大在线播放| 日韩精品有码人妻一区| 天天躁日日躁夜夜躁夜夜| 自线自在国产av| 啦啦啦 在线观看视频| 一区在线观看完整版| 久久人人爽av亚洲精品天堂| 久久久久久久国产电影| 99国产综合亚洲精品| 精品一区二区三区四区五区乱码 | 日韩制服骚丝袜av| 亚洲国产欧美网| 成年av动漫网址| 不卡视频在线观看欧美| 亚洲av成人精品一二三区| av在线播放精品| 热re99久久精品国产66热6| 18禁国产床啪视频网站| 91aial.com中文字幕在线观看| 欧美日韩亚洲国产一区二区在线观看 | 国产熟女午夜一区二区三区| 免费观看性生交大片5| 蜜桃国产av成人99| 人人妻人人澡人人爽人人夜夜| 看非洲黑人一级黄片| 人人澡人人妻人| 精品国产超薄肉色丝袜足j| 免费在线观看完整版高清| 婷婷色麻豆天堂久久| www.熟女人妻精品国产| 国产日韩欧美视频二区| 亚洲精品美女久久av网站| av在线观看视频网站免费| 国产日韩一区二区三区精品不卡| 久久久久久人妻| 亚洲av电影在线进入| 菩萨蛮人人尽说江南好唐韦庄| av卡一久久| av线在线观看网站| 亚洲欧美一区二区三区黑人| 男女高潮啪啪啪动态图| 婷婷色av中文字幕| 国产男女超爽视频在线观看| 高清av免费在线| 欧美日韩av久久| 男女无遮挡免费网站观看| 最近中文字幕2019免费版| 桃花免费在线播放| 丝袜脚勾引网站| 午夜福利乱码中文字幕| 欧美日韩视频精品一区| 色播在线永久视频| 日韩电影二区| 欧美精品一区二区免费开放| 午夜91福利影院| tube8黄色片| 最近最新中文字幕大全免费视频 | 免费日韩欧美在线观看| 欧美精品高潮呻吟av久久| 精品免费久久久久久久清纯 | 一二三四在线观看免费中文在| 1024香蕉在线观看| 国产精品无大码| 啦啦啦中文免费视频观看日本| 久久人人爽人人片av| 午夜免费男女啪啪视频观看| 亚洲一卡2卡3卡4卡5卡精品中文| 美国免费a级毛片| 国产精品久久久av美女十八| 亚洲精品久久久久久婷婷小说| 黑人欧美特级aaaaaa片| 精品亚洲乱码少妇综合久久| av网站免费在线观看视频| 激情五月婷婷亚洲| 亚洲精品国产av蜜桃| 欧美日韩综合久久久久久| 女人精品久久久久毛片| 超碰成人久久| 高清视频免费观看一区二区| 少妇人妻久久综合中文| 777久久人妻少妇嫩草av网站| av在线app专区| av在线老鸭窝| 亚洲色图综合在线观看| av电影中文网址| 不卡av一区二区三区| 飞空精品影院首页| 搡老岳熟女国产| 亚洲成av片中文字幕在线观看| 欧美日韩综合久久久久久| 天天躁夜夜躁狠狠久久av| 亚洲欧洲国产日韩| 激情五月婷婷亚洲| 男女国产视频网站| 久久 成人 亚洲| 国产 精品1| 精品人妻一区二区三区麻豆| 亚洲第一青青草原| 制服丝袜香蕉在线| 少妇人妻久久综合中文| 亚洲精品久久午夜乱码| 国产无遮挡羞羞视频在线观看| 1024香蕉在线观看| 美女脱内裤让男人舔精品视频| 十分钟在线观看高清视频www| 97人妻天天添夜夜摸| 久久天堂一区二区三区四区| 观看av在线不卡| 色播在线永久视频| 宅男免费午夜| 精品一区二区免费观看| 一级黄片播放器| av福利片在线| 青春草亚洲视频在线观看| 在线亚洲精品国产二区图片欧美| 成人毛片60女人毛片免费| 中国国产av一级| 亚洲av日韩精品久久久久久密 | 国产成人欧美在线观看 | 色婷婷av一区二区三区视频| av国产久精品久网站免费入址| 成年av动漫网址| a级片在线免费高清观看视频| 狠狠婷婷综合久久久久久88av| 一区二区av电影网| 午夜福利视频精品| 无遮挡黄片免费观看| 欧美xxⅹ黑人| 日本猛色少妇xxxxx猛交久久| 亚洲欧美一区二区三区国产| 男人添女人高潮全过程视频| 日韩av在线免费看完整版不卡| 久久久久精品国产欧美久久久 | 考比视频在线观看| 亚洲成色77777| 中文精品一卡2卡3卡4更新| 在线免费观看不下载黄p国产| 欧美黑人欧美精品刺激| 日韩不卡一区二区三区视频在线| 狂野欧美激情性bbbbbb| 香蕉丝袜av| 中文字幕最新亚洲高清| 午夜影院在线不卡| 两性夫妻黄色片| 丰满少妇做爰视频| 中文字幕人妻丝袜一区二区 | 午夜福利视频在线观看免费| 91aial.com中文字幕在线观看| 黄色怎么调成土黄色| 精品亚洲成国产av| 国产精品偷伦视频观看了| 99久久99久久久精品蜜桃| 秋霞伦理黄片| 性少妇av在线| 少妇 在线观看| 看十八女毛片水多多多| 久久影院123| 老熟女久久久| 久久99精品国语久久久| 热99国产精品久久久久久7| 黄片小视频在线播放| 一级毛片 在线播放| 秋霞在线观看毛片| 国产福利在线免费观看视频| 91精品伊人久久大香线蕉| 国产精品香港三级国产av潘金莲 | 久久午夜综合久久蜜桃| 精品亚洲成国产av| 91aial.com中文字幕在线观看| av卡一久久| 欧美97在线视频| 超色免费av| 啦啦啦视频在线资源免费观看| 久久久国产欧美日韩av| 日韩人妻精品一区2区三区| 国产成人av激情在线播放| www.自偷自拍.com| 亚洲av中文av极速乱| 啦啦啦视频在线资源免费观看| 亚洲欧美色中文字幕在线| 午夜日本视频在线| 欧美乱码精品一区二区三区| 国产不卡av网站在线观看| 欧美人与性动交α欧美精品济南到| 高清不卡的av网站| 赤兔流量卡办理| 精品国产一区二区久久| 韩国高清视频一区二区三区| 亚洲一码二码三码区别大吗| 9191精品国产免费久久| 欧美人与善性xxx| 在线观看三级黄色| 亚洲精品一区蜜桃| 在线观看免费午夜福利视频| 在线免费观看不下载黄p国产| 99精国产麻豆久久婷婷| 高清不卡的av网站| 高清视频免费观看一区二区| 男女高潮啪啪啪动态图| 日韩视频在线欧美| 精品国产一区二区久久| 欧美成人精品欧美一级黄| 天天操日日干夜夜撸| 在线观看免费高清a一片| 一区二区日韩欧美中文字幕| 成年女人毛片免费观看观看9 | 久久 成人 亚洲| 国产日韩欧美亚洲二区| 你懂的网址亚洲精品在线观看| 亚洲一级一片aⅴ在线观看| 2021少妇久久久久久久久久久| 青青草视频在线视频观看| 青春草亚洲视频在线观看| 国产精品99久久99久久久不卡 | 又大又爽又粗| 欧美少妇被猛烈插入视频| 一级黄片播放器| 老汉色av国产亚洲站长工具| 男女下面插进去视频免费观看| 一级片'在线观看视频| 欧美日韩国产mv在线观看视频| 老汉色∧v一级毛片| 亚洲视频免费观看视频| 久久人妻熟女aⅴ| 日韩av免费高清视频| 黄色视频在线播放观看不卡| 欧美在线黄色| 日韩大码丰满熟妇| 人人妻人人添人人爽欧美一区卜| av线在线观看网站| 在线 av 中文字幕| 极品少妇高潮喷水抽搐| 亚洲av中文av极速乱| 又粗又硬又长又爽又黄的视频| 国产片特级美女逼逼视频| 午夜免费观看性视频| 亚洲精品aⅴ在线观看| 精品第一国产精品| 成年av动漫网址| 国产一区二区激情短视频 | xxxhd国产人妻xxx| 久久久精品区二区三区| 五月开心婷婷网| 国产一区二区三区av在线| 最近最新中文字幕免费大全7| 777久久人妻少妇嫩草av网站| 天天躁狠狠躁夜夜躁狠狠躁| av在线观看视频网站免费| 亚洲情色 制服丝袜| 97精品久久久久久久久久精品| 80岁老熟妇乱子伦牲交| 男女床上黄色一级片免费看| 亚洲一码二码三码区别大吗| 伦理电影免费视频| 久久久久久免费高清国产稀缺| 亚洲国产毛片av蜜桃av| 国产 一区精品| 久久精品亚洲熟妇少妇任你| videos熟女内射| 九草在线视频观看| 观看av在线不卡| 不卡视频在线观看欧美| 色吧在线观看| 国产精品女同一区二区软件| 久久久精品94久久精品| 午夜老司机福利片| 欧美激情极品国产一区二区三区| 成人亚洲精品一区在线观看| 中文字幕另类日韩欧美亚洲嫩草| 国产精品 国内视频| 街头女战士在线观看网站| 亚洲国产精品一区二区三区在线| 国产欧美日韩综合在线一区二区| 人人妻人人添人人爽欧美一区卜| 国产精品免费视频内射| 九色亚洲精品在线播放| 美女国产高潮福利片在线看| 日本欧美国产在线视频| 老司机深夜福利视频在线观看 | videos熟女内射| 午夜精品国产一区二区电影| 欧美黄色片欧美黄色片| 人体艺术视频欧美日本| av福利片在线| xxx大片免费视频| 久久 成人 亚洲| 亚洲自偷自拍图片 自拍| 免费黄色在线免费观看| 国产成人系列免费观看| 免费高清在线观看日韩| 制服诱惑二区| 性少妇av在线| 人人澡人人妻人| 日本爱情动作片www.在线观看| 好男人视频免费观看在线| 国产乱来视频区| 99久久人妻综合| 一本一本久久a久久精品综合妖精| 日韩人妻精品一区2区三区| 欧美少妇被猛烈插入视频| 亚洲精品国产av成人精品| 亚洲国产精品成人久久小说| 老司机影院成人| av不卡在线播放| 日韩免费高清中文字幕av| 热re99久久国产66热| 一二三四中文在线观看免费高清| 精品一区二区三区av网在线观看 | 亚洲欧美成人精品一区二区| 久久精品熟女亚洲av麻豆精品| 精品亚洲成a人片在线观看| 午夜福利视频在线观看免费| 国产色婷婷99| 人妻 亚洲 视频| 欧美激情 高清一区二区三区| 久久精品国产亚洲av涩爱| 国产在线视频一区二区| 国产伦人伦偷精品视频| 纵有疾风起免费观看全集完整版| 亚洲色图 男人天堂 中文字幕| www.熟女人妻精品国产| 国产老妇伦熟女老妇高清| 国产精品久久久久久精品电影小说| 久久久久国产精品人妻一区二区| 欧美黑人欧美精品刺激| 亚洲国产av影院在线观看| 看免费av毛片| 女的被弄到高潮叫床怎么办| 亚洲男人天堂网一区| av网站免费在线观看视频| 亚洲av成人精品一二三区| 我的亚洲天堂| 久久久久国产精品人妻一区二区| 久久久国产欧美日韩av| 90打野战视频偷拍视频| 亚洲欧洲国产日韩| 精品第一国产精品| 成年人午夜在线观看视频| 欧美 日韩 精品 国产| 亚洲国产精品999| 最近最新中文字幕大全免费视频 | 欧美激情 高清一区二区三区| 亚洲精品自拍成人| 亚洲av中文av极速乱| 欧美国产精品一级二级三级| 老司机在亚洲福利影院| 日韩一区二区视频免费看| 久久久久久久久久久久大奶| 激情五月婷婷亚洲| 亚洲精品成人av观看孕妇| 国产精品久久久av美女十八| 亚洲欧美成人精品一区二区| 女人久久www免费人成看片| 久久女婷五月综合色啪小说| 午夜福利网站1000一区二区三区| 大香蕉久久成人网| 最黄视频免费看| 日韩av免费高清视频| 人人澡人人妻人| 久久韩国三级中文字幕| 色婷婷久久久亚洲欧美| 一区二区三区激情视频| 黑人巨大精品欧美一区二区蜜桃| 精品亚洲成国产av| 日韩,欧美,国产一区二区三区| 大码成人一级视频| 黄网站色视频无遮挡免费观看| 在线精品无人区一区二区三| 亚洲av日韩精品久久久久久密 | 精品卡一卡二卡四卡免费| 国产精品一二三区在线看| 日韩一区二区视频免费看| 51午夜福利影视在线观看| 精品国产乱码久久久久久小说| 一级毛片电影观看| 99久国产av精品国产电影| 新久久久久国产一级毛片| 18禁国产床啪视频网站| 亚洲av国产av综合av卡| 久久久国产一区二区| 丝袜美足系列| 一个人免费看片子| 国产精品 国内视频| 精品第一国产精品| 亚洲av在线观看美女高潮| 满18在线观看网站| 国产激情久久老熟女| 777久久人妻少妇嫩草av网站| 蜜桃在线观看..| 制服诱惑二区| 9191精品国产免费久久| 亚洲欧美色中文字幕在线| 1024视频免费在线观看| 少妇猛男粗大的猛烈进出视频| 亚洲精品久久成人aⅴ小说| 日本色播在线视频| 国产精品久久久久成人av| 精品人妻在线不人妻| 香蕉国产在线看| 咕卡用的链子| 青春草视频在线免费观看| 视频在线观看一区二区三区| 国产精品女同一区二区软件| 黑人巨大精品欧美一区二区蜜桃| 男女免费视频国产| 亚洲精品国产色婷婷电影| 亚洲精品国产av蜜桃| 免费女性裸体啪啪无遮挡网站| 男女边吃奶边做爰视频| 啦啦啦啦在线视频资源| 90打野战视频偷拍视频| 亚洲精品国产一区二区精华液| 激情视频va一区二区三区| 日韩熟女老妇一区二区性免费视频| 中文精品一卡2卡3卡4更新| 国产片特级美女逼逼视频| 我的亚洲天堂| 日韩制服丝袜自拍偷拍| 男的添女的下面高潮视频| 久久青草综合色| 免费黄频网站在线观看国产| 久久免费观看电影| 在线观看免费高清a一片| 秋霞伦理黄片| 精品国产乱码久久久久久小说| 国产片特级美女逼逼视频| 国产探花极品一区二区| 成人国语在线视频| 国产亚洲av高清不卡| 成人亚洲欧美一区二区av| 久久久久精品国产欧美久久久 | 久久久久久久国产电影| 成人毛片60女人毛片免费| 久久人人爽av亚洲精品天堂| 女人高潮潮喷娇喘18禁视频| 色精品久久人妻99蜜桃| 日日爽夜夜爽网站| 精品视频人人做人人爽| av免费观看日本| 中文乱码字字幕精品一区二区三区| 精品人妻熟女毛片av久久网站| 中文精品一卡2卡3卡4更新| 久久久久久久大尺度免费视频| 日韩av不卡免费在线播放| 啦啦啦中文免费视频观看日本| 在线观看人妻少妇| 亚洲天堂av无毛| 五月天丁香电影| 99国产综合亚洲精品| 色播在线永久视频| 欧美激情极品国产一区二区三区| 久久久久久久精品精品| 欧美老熟妇乱子伦牲交| 久久久国产一区二区| 成年av动漫网址| 国产成人精品久久二区二区91 | 国产日韩欧美视频二区| 岛国毛片在线播放| 久久久精品94久久精品| 男女午夜视频在线观看| 黄片无遮挡物在线观看| 亚洲精品日本国产第一区| 90打野战视频偷拍视频| 中国国产av一级| avwww免费| 精品一区二区三区四区五区乱码 | 国产精品三级大全| 嫩草影视91久久| 成人手机av| 国产精品久久久久久精品电影小说| 大陆偷拍与自拍| 国产成人a∨麻豆精品| 18禁国产床啪视频网站| 中文字幕色久视频| av在线app专区| 最黄视频免费看| 精品国产超薄肉色丝袜足j| 欧美日韩一区二区视频在线观看视频在线| 国产一区二区 视频在线| 国产色婷婷99| 极品人妻少妇av视频| 建设人人有责人人尽责人人享有的| 中文天堂在线官网| 婷婷色综合www| 亚洲精品美女久久久久99蜜臀 | 免费观看a级毛片全部| 亚洲欧美精品综合一区二区三区| 久久女婷五月综合色啪小说| 黑人欧美特级aaaaaa片| 国产淫语在线视频| 又大又黄又爽视频免费| 久久久精品区二区三区| 男的添女的下面高潮视频| 免费日韩欧美在线观看| 在线看a的网站| 久久国产亚洲av麻豆专区| 中文字幕人妻丝袜一区二区 | 久久久久久久国产电影| 亚洲国产欧美一区二区综合| 欧美日韩国产mv在线观看视频| 捣出白浆h1v1| 免费不卡黄色视频| a 毛片基地| 免费观看人在逋| 精品一区在线观看国产| 久久久精品区二区三区| 日韩av在线免费看完整版不卡| 99九九在线精品视频| 女人久久www免费人成看片| 久久久久久久精品精品| 久久亚洲国产成人精品v| 丰满乱子伦码专区| 成人午夜精彩视频在线观看| 多毛熟女@视频| 久久毛片免费看一区二区三区| 国产激情久久老熟女| 日韩电影二区| 亚洲免费av在线视频| 999久久久国产精品视频| 久久久久精品久久久久真实原创| 黄色视频在线播放观看不卡| 国产精品久久久久久人妻精品电影 | 亚洲国产日韩一区二区| 久久ye,这里只有精品| 精品一品国产午夜福利视频| 国产成人免费观看mmmm| 亚洲图色成人| 一区二区av电影网| 九九爱精品视频在线观看| 国产 一区精品| 啦啦啦啦在线视频资源| 亚洲精品久久久久久婷婷小说| a 毛片基地| 精品福利永久在线观看| 免费日韩欧美在线观看| 男女之事视频高清在线观看 | 一本色道久久久久久精品综合| 伊人久久大香线蕉亚洲五| 精品国产一区二区久久| 伊人久久大香线蕉亚洲五| 国产一区二区在线观看av| 一个人免费看片子| 狂野欧美激情性xxxx| 国产日韩欧美在线精品| 久久久久视频综合| 丝袜喷水一区| 精品久久蜜臀av无| 一级爰片在线观看| 成年女人毛片免费观看观看9 | 搡老乐熟女国产| 国产午夜精品一二区理论片| 午夜久久久在线观看| 久热这里只有精品99| 免费av中文字幕在线| 欧美精品一区二区免费开放| 欧美精品高潮呻吟av久久| 999久久久国产精品视频| 国产免费一区二区三区四区乱码| 久久久久久久久免费视频了| 亚洲熟女毛片儿| 夫妻性生交免费视频一级片| 亚洲欧洲精品一区二区精品久久久 | 在线观看免费午夜福利视频| 成年人午夜在线观看视频| 国产一卡二卡三卡精品 | 精品国产国语对白av| 亚洲第一青青草原| 亚洲国产中文字幕在线视频| 美女中出高潮动态图| 99热网站在线观看| 久久久久久久国产电影| 日韩伦理黄色片| 自线自在国产av| 久久ye,这里只有精品| 精品人妻熟女毛片av久久网站| 黄色怎么调成土黄色| 大香蕉久久网|