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

    草魚(yú)載脂蛋白A-I-1基因3'非編碼區(qū)SNPs篩選及其與生長(zhǎng)性狀的關(guān)聯(lián)分析

    2012-03-20 14:52:20劉小獻(xiàn)白俊杰于凌云李勝杰
    關(guān)鍵詞:生長(zhǎng)

    劉小獻(xiàn),白俊杰,于凌云,李勝杰

    (1.大連海洋大學(xué)水產(chǎn)與生命學(xué)院,遼寧大連116023;2.中國(guó)水產(chǎn)科學(xué)研究院珠江水產(chǎn)研究所,農(nóng)業(yè)部熱帶亞熱帶水產(chǎn)資源利用與養(yǎng)殖重點(diǎn)實(shí)驗(yàn)室,廣州510380)

    草魚(yú)載脂蛋白A-I-1基因3'非編碼區(qū)SNPs篩選及其與生長(zhǎng)性狀的關(guān)聯(lián)分析

    劉小獻(xiàn)1、2,白俊杰2,于凌云2,李勝杰2

    (1.大連海洋大學(xué)水產(chǎn)與生命學(xué)院,遼寧大連116023;2.中國(guó)水產(chǎn)科學(xué)研究院珠江水產(chǎn)研究所,農(nóng)業(yè)部熱帶亞熱帶水產(chǎn)資源利用與養(yǎng)殖重點(diǎn)實(shí)驗(yàn)室,廣州510380)

    采用PCR產(chǎn)物直接測(cè)序法和PCR-RFLP法對(duì)草魚(yú)EST文庫(kù)中載脂蛋白A-I-1基因(Apoprotein A-I -1,apoA-I-1)3'非編碼區(qū)進(jìn)行SNPs篩選,在珠江水系草魚(yú)Ctenopharyngodon idella養(yǎng)殖群體144個(gè)樣本中發(fā)現(xiàn)2個(gè)SNPs位點(diǎn)。C792T位點(diǎn)的CC型占46.53%,CT型占50.69%,TT型占2.78%;G851A位點(diǎn)的GG型占77.08%,GA型占20.83%,AA型占2.08%。采用一般線性模型對(duì)2個(gè)SNPs位點(diǎn)與草魚(yú)4個(gè)生長(zhǎng)性狀——體質(zhì)量、體長(zhǎng)、體高和頭長(zhǎng)進(jìn)行關(guān)聯(lián)分析,結(jié)果表明:C792T位點(diǎn)TT型個(gè)體的體質(zhì)量、體長(zhǎng)、體高和頭長(zhǎng)均值高于CT型和CC型個(gè)體的生長(zhǎng)性狀指標(biāo)值,但未達(dá)到顯著水平(P>0.05);G851A位點(diǎn)GG型個(gè)體的體質(zhì)量、體長(zhǎng)、體高和頭長(zhǎng)均值高于GA型和AA型個(gè)體的生長(zhǎng)性狀指標(biāo)值,且GG型個(gè)體的體長(zhǎng)和頭長(zhǎng)均值與GA型和AA型個(gè)體均存在顯著差異(P<0.05)。將2個(gè)SNPs位點(diǎn)不同基因型組合成6種雙倍型,關(guān)聯(lián)分析結(jié)果表明,雙倍型D3(TTC792T-GGG851A)個(gè)體的體質(zhì)量均值最高,D6型(CCC792TAAG851A)個(gè)體的體質(zhì)量均值最低,且二者在體長(zhǎng)、體高、頭長(zhǎng)均值間存在顯著差異(P<0.05)。推測(cè)C792T位點(diǎn)TT型為有利基因型;G851A位點(diǎn)GG型為有利基因型,AA型為不利基因型。本研究結(jié)果為草魚(yú)分子輔助育種提供了候選標(biāo)記,為加快草魚(yú)育種進(jìn)程奠定了基礎(chǔ)。

    草魚(yú);SNP;載脂蛋白A-I;apoA-I-1;生長(zhǎng)性狀;雙倍型;關(guān)聯(lián)分析

    草魚(yú)Ctenopharyngodon idella隸屬于鯉形目、鯉科、雅羅魚(yú)亞科、草魚(yú)屬,是中國(guó)淡水養(yǎng)殖四大家魚(yú)之一,有著重要的經(jīng)濟(jì)價(jià)值。由于草魚(yú)親魚(yú)的個(gè)體大、性成熟周期長(zhǎng),給其良種選育帶來(lái)了困難,目前尚未獲得遺傳改良的品種或品系。利用DNA遺傳標(biāo)記進(jìn)行輔助選育已成為動(dòng)物育種過(guò)程中的重要方法[1]。篩選并利用與草魚(yú)經(jīng)濟(jì)性狀相關(guān)的分子標(biāo)記,可有效提高選擇效率,縮短育種年限,對(duì)加快培育具有優(yōu)良經(jīng)濟(jì)性狀的草魚(yú)養(yǎng)殖品系具有重要的意義。

    魚(yú)類(lèi)能量主要來(lái)源于脂類(lèi),脂類(lèi)代謝在魚(yú)類(lèi)的能量代謝及生長(zhǎng)過(guò)程中有著重要的作用。載脂蛋白A-I(Apoprotein A-I,簡(jiǎn)稱(chēng)apoA-I)是硬骨魚(yú)類(lèi)血漿蛋白中高密度脂蛋白(High-density lipoprotein,HDL)的重要組成成分,在逆向轉(zhuǎn)運(yùn)膽固醇過(guò)程中具有重要作用[2],對(duì)魚(yú)類(lèi)生長(zhǎng)存在潛在的作用。一些學(xué)者在對(duì)鯉、虹鱒和斑點(diǎn)叉尾鮰的研究中發(fā)現(xiàn),apoA-I還存在潛在的抗菌活性,具有部分免疫功能,對(duì)魚(yú)類(lèi)的健康生長(zhǎng)具有一定的保護(hù)作用[2-4]。本試驗(yàn)中,作者采用PCR產(chǎn)物直接測(cè)序法和PCR-RFLP法對(duì)草魚(yú)EST文庫(kù)中載脂蛋白A-I-1基因(Apoprotein A-I-1,apoA-I-1)3'非編碼區(qū)SNPs位點(diǎn)進(jìn)行篩選,以期為草魚(yú)分子標(biāo)記輔助育種提供候選標(biāo)記。

    1 材料與方法

    1.1 材料

    試驗(yàn)用草魚(yú)取自廣東省佛山市三水白金水產(chǎn)種苗有限公司,為同批繁殖且同塘養(yǎng)殖的20月齡草

    魚(yú)群體,樣本數(shù)量為144尾,體質(zhì)量均值為1 292.5 g。親本來(lái)自珠江水系。

    Taq DNA聚合酶、buffer和dNTP購(gòu)自華美生物工程公司,限制性?xún)?nèi)切酶BseGI和MspI購(gòu)自Fermentas公司,瓊脂糖、去離子甲酰胺丙烯酰胺和N,N-亞甲雙丙烯酰胺等購(gòu)自廣州威佳科技有限公司。

    1.2 方法

    1.2.1 DNA的提取 采用尾靜脈活體取血,抗凝劑(ACD)與血液體積比為6∶1,使用天根生化科技(北京)有限公司(TIANGEN)血液基因組提取試劑盒提取DNA,樣本于-20℃下保存。

    1.2.2 引物設(shè)計(jì)及合成 根據(jù)草魚(yú)EST-SNPs數(shù)據(jù)庫(kù)中apoA-I-1基因cDNA序列設(shè)計(jì)引物,擴(kuò)增含SNPs位點(diǎn)的序列,引物序列如下:

    PF:GGAGAACTGAGGAAAAAGATGACCG;

    PR:TTTATTCATATCTGTGTTTGCCCATTAG,由上海英駿生物技術(shù)有限公司合成。

    1.2.3 SNPs篩選及基因型的檢測(cè) 從144個(gè)草魚(yú)DNA樣本中隨機(jī)挑選20個(gè)樣本,用引物PF和PR對(duì)其進(jìn)行PCR擴(kuò)增。PCR反應(yīng)體積為20 μL,包括10×PCR buffer 2.0 μL,dNTP(10 mmol/L)0.3 μL,Taq酶(5 U/μL)0.2 μL,上下游引物(20 μmol/L)各0.4 μL,基因組DNA(80 ng/μL)1 μL,ddH2O 15.8 μL。反應(yīng)程序?yàn)?94℃下預(yù)變性4 min;94℃下變性30 s,56℃下退火30 s,72℃下延伸30 s,共進(jìn)行32個(gè)循環(huán);最后在72℃下延伸7 min,12℃下保存。將擴(kuò)增產(chǎn)物在10 g/L的瓊脂糖凝膠中進(jìn)行電泳檢測(cè)。PCR產(chǎn)物經(jīng)純化后送廣州博尚生物技術(shù)有限公司進(jìn)行測(cè)序,用Vector NIT Suite 8.0軟件分析測(cè)序結(jié)果,查找該片段上的SNPs位點(diǎn)。

    1.2.4 草魚(yú)apoA-I-1基因3'UTR序列分析 采用Vector NIT Suite 8.0軟件和Clustalx 1.83軟件將草魚(yú)apoA-I-1基因3'UTR序列與GenBank中已登錄的白鰱(登錄號(hào)為HM124749)、稀有鮈鯽(登錄號(hào)為EU327775)、花魚(yú)骨(登錄號(hào)為FJ170104, FJ170105)3種鯉科魚(yú)類(lèi)的apoA-I-1基因3'UTR序列進(jìn)行分析比較。

    1.2.5 采用PCR-RFLP法對(duì)SNPs位點(diǎn)進(jìn)行分型PCR產(chǎn)物經(jīng)10 g/L的瓊脂糖凝膠電泳檢測(cè)后,用限制性?xún)?nèi)切酶進(jìn)行酶切。酶切體系為:PCR產(chǎn)物0.20 μg,10×fast digest buffer 0.8 μL,限制性?xún)?nèi)切酶0.3 μL,ddH2O 6 μL。酶切產(chǎn)物用10%的聚丙烯酰胺凝膠電泳分離(Acr與Bis的質(zhì)量比為29∶1),以120 V恒壓電泳1.5 h后,用銀染法顯色。

    1.3 數(shù)據(jù)處理

    采用Popgene(Version 3.2)軟件分析基因型頻率和等位基因頻率。由于樣本是同一批繁育且同池養(yǎng)殖,采樣時(shí)間也一致,不存在時(shí)間、環(huán)境及人工飼養(yǎng)水平的差別,所以在建立模型時(shí)不考慮季節(jié)、環(huán)境及人工飼養(yǎng)技術(shù)差別,利用SPSS 15.0軟件一般線性模型(General Linear Model,GLM)對(duì)SNPs位點(diǎn)基因型與草魚(yú)主要生長(zhǎng)性狀之間的相關(guān)性進(jìn)行最小二乘分析。統(tǒng)計(jì)分析模型為

    Yij=μ+Bi+eij,

    其中:Yij為某個(gè)性狀第i個(gè)標(biāo)記第j個(gè)個(gè)體的觀測(cè)值;μ為試驗(yàn)觀測(cè)所有個(gè)體的平均值(即總體平均值);Bi為第i個(gè)標(biāo)記的效應(yīng)值;eij為對(duì)應(yīng)于觀察值的隨機(jī)誤差效應(yīng)。試驗(yàn)數(shù)據(jù)用平均值±標(biāo)準(zhǔn)誤差表示。

    2 結(jié)果

    2.1 草魚(yú)apoA-I-1基因cDNA序列及SNPs位點(diǎn)的初步分析

    由草魚(yú)EST文庫(kù)獲得草魚(yú)apoA-I-1基因cDNA序列1 291 bp,與GenBank已登錄序列分析比較,該序列含CDS序列774 bp,5'UTR序列177 bp,3'UTR序列346 bp,共編碼257個(gè)氨基酸。經(jīng)序列比對(duì)分析,發(fā)現(xiàn)在第792個(gè)堿基處(相對(duì)于起始密碼子ATG)存在C-T突變,命名為C792T位點(diǎn);在第851個(gè)堿基處(相對(duì)于起始密碼子ATG)存在G-A突變,命名為G851A位點(diǎn)(圖1)。

    2.2 草魚(yú)apoA-I-1基因3'UTR SNPs位點(diǎn)基因型的檢測(cè)

    采用限制性?xún)?nèi)切酶BseGI對(duì)C792T位點(diǎn)進(jìn)行基因型分型,該314 bp序列存在2個(gè)BseGI酶切位點(diǎn)。經(jīng)酶切,CC基因型個(gè)體具有173、105、36 bp 3條帶,CT基因型個(gè)體具有173、121、105、52、36 bp 5條帶,TT基因型個(gè)體具有121、105、52、36 bp 4條帶,較小的片段位于圖譜下方,因電泳時(shí)間長(zhǎng),顏色較淡,在此僅標(biāo)出可鑒別基因型的173 bp和121 bp兩個(gè)片段(圖2)。用限制性?xún)?nèi)切酶MspI對(duì)G851A位點(diǎn)進(jìn)行基因型分型,經(jīng)酶切, GG基因型個(gè)體具有177、137 bp 2條帶,GA基因型個(gè)體具有314、177、137 bp 3條帶,AA基因型

    個(gè)體只有314 bp 1條帶(圖3)。

    圖1 草魚(yú)apoA-I-1基因cDNA部分序列及SNPs位點(diǎn)Fig.1 Partical cDNA of apoA-I-1 gene and SNPs sites in grass carp

    圖2 草魚(yú)apoA-I-1基因C792T位點(diǎn)BseGI酶切電泳圖譜Fig.2 Electrophoresis of C792T site in apoA-I-1 gene demonstrated by BseGI PCR-RFLP in grass carp

    圖3 草魚(yú)apoa-I-1基因G851A位點(diǎn)MspI酶切電泳圖譜Fig.3 Electrophoresis of G851A site in apoA-I-1 gene demonstrated by MspI PCR-RFLP in grass carp

    2.3 草魚(yú)apoA-I-1基因3'UTR SNPs位點(diǎn)在群體中的分布

    草魚(yú)apoA-I-1基因C792T位點(diǎn)和G851A位點(diǎn)3種基因型及其等位基因頻率的統(tǒng)計(jì)結(jié)果見(jiàn)表1和表2。從表1、表2可見(jiàn):在該試驗(yàn)群體中,C792T位點(diǎn)的CT型為優(yōu)勢(shì)基因型,等位基因C頻率明顯高于等位基因T,C為優(yōu)勢(shì)等位基因;G851A位點(diǎn)的GG型為優(yōu)勢(shì)基因型,GG型所占比例是GA型的3.7倍,G為優(yōu)勢(shì)等位基因。

    表1 草魚(yú)apoA-I-1基因3'UTR C792T位點(diǎn)在群體中的基因頻率分布Tab.1 Frequencies of alleles and genotypes of C792T in 3'UTR of apoA-I-1 gene in grass carp population%

    表2 草魚(yú)apoA-I-1基因3'UTR G851A位點(diǎn)在群體中的基因頻率分布Tab.2 Frequencies of alleles and genotypes of G851A in 3'UTR of apoA-I-1 gene in grass carp population%

    2.4 草魚(yú)apoA-I-1基因3'UTR SNPs位點(diǎn)與生長(zhǎng)性狀的關(guān)聯(lián)分析

    對(duì)兩個(gè)SNPs位點(diǎn)不同基因型與草魚(yú)4個(gè)主要生長(zhǎng)性狀的關(guān)聯(lián)分析結(jié)果見(jiàn)表3。從表3可見(jiàn): C792T位點(diǎn)TT型個(gè)體的4個(gè)主要生長(zhǎng)性狀均值最高,CT型的體質(zhì)量和體高值高于CC型個(gè)體,T為有利等位基因;G851A位點(diǎn)GG型個(gè)體的4個(gè)生長(zhǎng)性狀均值最高,AA型個(gè)體的各項(xiàng)值最小,且GG型個(gè)體的體長(zhǎng)和頭長(zhǎng)與GA型、AA型均存在顯著差異(P<0.05)。

    表3 草魚(yú)apoA-I-1基因C792T位點(diǎn)和G851A位點(diǎn)不同基因型與草魚(yú)生長(zhǎng)性狀的關(guān)聯(lián)分析Tab.3 Association of C792T site and G851A site of apoA-I-1 gene polymorphisms with growth traits in grass carp

    將兩個(gè)多態(tài)位點(diǎn)的不同基因型進(jìn)行組合,得到6種雙倍型(表4)。從表4可見(jiàn):D3型個(gè)體的4個(gè)生長(zhǎng)性狀均值最高,D6型個(gè)體的4個(gè)生長(zhǎng)性狀均值最低。其中D3型個(gè)體的體長(zhǎng)、體高均值顯著高于D4、D6型個(gè)體(P<0.05),D6型個(gè)體的體長(zhǎng)均值顯著低于D1~D5型個(gè)體;D3型個(gè)體的頭長(zhǎng)均值極顯著高于D4、D6型個(gè)體(P<0.01),D6型個(gè)體的頭長(zhǎng)均值顯著低于D1型個(gè)體(P<0.05)。

    3 討論

    SNPs主要是指在基因組水平上由于單個(gè)核苷酸的變異所引起的DNA序列多態(tài)性,即同一物種不同個(gè)體間染色體上遺傳密碼子單個(gè)堿基的變化。常表現(xiàn)為單堿基的轉(zhuǎn)換、顛換、插入和缺失[5],通常發(fā)生轉(zhuǎn)換與顛換的SNPs之比為2∶1[6]。本試驗(yàn)中兩個(gè)多態(tài)位點(diǎn)都屬于轉(zhuǎn)換,且只有兩個(gè)等位基因。目前,國(guó)內(nèi)外采用SNPs方法對(duì)魚(yú)類(lèi)生長(zhǎng)方面的研究已有很多報(bào)道,Tao等[7]對(duì)北極嘉魚(yú)Salvelinus alpinus的研究發(fā)現(xiàn),GHRH/PACAP2基因上的SNPs位點(diǎn)與北極嘉魚(yú)早期生長(zhǎng)速率極顯著相關(guān); Case等[8]對(duì)大西洋鱈Gadus morhua和其它海域鱈panⅠ基因的研究結(jié)果表明,大西洋鱈panⅠ基因的變異對(duì)其生長(zhǎng)性狀有著重要的影響;Li等[9]對(duì)大口黑鱸Micropterus salmoides IGF-I基因5'非編碼區(qū)進(jìn)行了研究,結(jié)果表明,啟動(dòng)子區(qū)存在多態(tài)位

    點(diǎn),且與大口黑鱸的生長(zhǎng)性狀有相關(guān)性。

    表4 草魚(yú)apoA-I-1基因不同雙倍型與草魚(yú)生長(zhǎng)性狀的關(guān)聯(lián)分析Tab.4 Association of apoA-I-1 gene diplotype with growth traits in grass carp

    本研究中,在草魚(yú)apoA-I-1基因的3'非編碼區(qū)發(fā)現(xiàn)兩個(gè)SNPs位點(diǎn),其中C792T位點(diǎn)TT型個(gè)體的生長(zhǎng)性狀優(yōu)于CC型和CT型個(gè)體;G851A位點(diǎn)GG型個(gè)體的生長(zhǎng)性狀優(yōu)于GA型和AA型個(gè)體。將草魚(yú)、白鰱、稀有鮈鯽和花魚(yú)骨4種鯉科魚(yú)類(lèi)apoA-I-1基因3'UTR序列進(jìn)行比較,發(fā)現(xiàn)C792T位點(diǎn)處存在兩種堿基類(lèi)型,白鰱和稀有鮈鯽的堿基為C,花魚(yú)骨的堿基為T(mén);G851A位點(diǎn)處白鰱、稀有鮈鯽和花魚(yú)骨3種魚(yú)的堿基類(lèi)型都為G,推測(cè)草魚(yú)C966T位點(diǎn)的TT型為突變型,且該突變對(duì)草魚(yú)的生長(zhǎng)有利,而G851A位點(diǎn)的AA型也為突變型,對(duì)草魚(yú)的生長(zhǎng)不利。

    魚(yú)類(lèi)的生長(zhǎng)性狀屬于數(shù)量性狀,由多個(gè)基因控制,不同的SNPs位點(diǎn)間可能存在相互作用,由單倍型構(gòu)成的雙倍型基因比單個(gè)SNPs能提供更多的基因型頻率信息[10]。本試驗(yàn)中對(duì)草魚(yú)apoA-I-1基因不同雙倍型個(gè)體與草魚(yú)的生長(zhǎng)性狀進(jìn)行關(guān)聯(lián)分析,發(fā)現(xiàn)D3型(TTC792T-GGG851A)個(gè)體的體質(zhì)量、體長(zhǎng)和體高均值最高,D6型(CCC792TAAG851A)個(gè)體的體質(zhì)量、體長(zhǎng)和體高均值最低,可見(jiàn)D3型是有利基因型。由此可見(jiàn),當(dāng)兩個(gè)位點(diǎn)都為有利基因型時(shí),其生長(zhǎng)性狀均值最高;當(dāng)兩個(gè)多態(tài)位點(diǎn)中存在一個(gè)有利基因型或者兩個(gè)位點(diǎn)都為不利基因型時(shí),其生長(zhǎng)性狀均值最低。C792T位點(diǎn)的T為有利等位基因,且在群體中所占的比例不高,在草魚(yú)的分子輔助育種中具有較大的選擇潛力。在選育時(shí)可以考慮以apoA-I-1基因作為候選基因,應(yīng)用于草魚(yú)的分子輔助育種,選擇C792T位點(diǎn)為T(mén)T型以及G851A為GG型的個(gè)體,淘汰G851A位點(diǎn)為AA型的個(gè)體。

    目前對(duì)基因表達(dá)調(diào)控的研究主要集中在基因的5'端,其實(shí)3'端序列亦含有豐富的信息量,基因的3'非翻譯區(qū)涉及了基因轉(zhuǎn)錄后水平的調(diào)控過(guò)程,主要參與mRNA的穩(wěn)定性及降解速率、翻譯時(shí)間和位點(diǎn)、控制翻譯起始及效率等[11-14]。有關(guān)魚(yú)類(lèi)apoA-I-1基因的3'UTR調(diào)控元件的研究還很少,其區(qū)域堿基突變而引起生長(zhǎng)性狀改變的作用機(jī)理尚不清楚,還有待深入研究。

    本研究中對(duì)草魚(yú)apoA-I-1基因cDNA序列進(jìn)行了分析,并在3'UTR發(fā)現(xiàn)了兩個(gè)SNPs位點(diǎn),通過(guò)一般線性模型分析,了解到這些突變位點(diǎn)不同基因型與草魚(yú)的生長(zhǎng)性狀存在著一定的關(guān)聯(lián),這將有利于解決草魚(yú)性成熟周期長(zhǎng)等問(wèn)題,并為培育出草魚(yú)的優(yōu)良品系提供參考資料。

    [1] 尹紹武,黃海,雷從改,等.DNA標(biāo)記技術(shù)在魚(yú)類(lèi)遺傳育種中的應(yīng)用[J].水產(chǎn)科學(xué),2005,24(11):42-45.

    [2] Amthauer R,Villanueva J,Vera M,et al.Characterisation of the major plasma apolipoproteins of the high density lipoprotein in the carp(Cyprinus carpio)[J].Comparative Biochemistry&Physiology,1989,92B:787-793.

    [3] Concha M I,Smith V J,Castro K,et al.Apolipoproteins A-I and A -II are potentially important effectors of innate immunity in the teleost fish Cyprinus carpio[J].Eur J Biochem,2004,271:2984-2990.

    [4] Villarroel F,Bastías A,Casado A,et al.Apolipoprotein A-I,an antimicrobial protein in Oncorhynchus mykiss:evaluation of its expression in primary defense barriers and plasma levels in sick and healthy fish[J].Fish&Shellfish Immunol,2007,23:197-209.

    [5] Alain V,Denis M,Magali S C,et al.A review on SNP and other types of molecular markers and their use in animal genetics[J]. Genet Sel Evol,2002,34:275-305.

    [6] Brookes A J.The essence of SNPs[J].Gene,1999,234:177-186.

    [7] Tao W J,Boulding E G.Association between single nucleotide polymorphisms in candidate genes and growth rate in Arctic charr (Salvelinus alpinus L.)[J].Heredity,2003,91:60-69.

    [8] Case R A J,Hutchinson W F.Association between growth and Pan I genotype within Atlantic cod full-sibling families[J].American Fisheries Society,2006,135:241-250.

    [9] Li Xiaohui,Bai Junjie,Ye Xing,et al.Polymorphisms in the 5'flanking region of the insulin-like growth factor I gene are associated with growth traits in large mouth bass Micropterus salmoides [J].Fisheries Science,2009,75(2):351-358.

    [10] Stephens M,Smith N,Donnelly P.A new statistical method for haplotype reconstruction from population data[J].Am J Hum Genet,2001,68:978-989.

    [11] 江元清,凌毅,趙武玲.真核mRNA的3'非翻譯區(qū)轉(zhuǎn)錄調(diào)控作用研究進(jìn)展[J].植物學(xué)通報(bào),2001,18(1):3-10.

    [12] Xie X H,Lu J,Kulbokas E J,et al.Systematic discovery of regulatory motifs in human promoters and 3'-UTRs by comparison of several mammals[J].Nature,2005,434:338-345.

    [13] Wilkie G S,Dickson K S,Gray N K.Regulation of mRNA translation by 5'-and 3'-UTR binding factors[J].Trends in Biochemical Sciences,2003,28(4):182-188.

    [14] 陳淑華.真核生物mRNA 3'-非翻譯區(qū)的調(diào)控功能[J].國(guó)外醫(yī)學(xué):生理、病理科學(xué)與臨床分冊(cè),2003,23(12):611-614.

    SNPs screening of 3'UTR in Apoprotein A-I-1 gene and its association with growth traits in grass carp

    LIU Xiao-xian1,2,BAI Jun-jie2,YU Ling-yun2,LI Sheng-jie2
    (1.College of Fisheries and Life Science,Dalian Ocean University,Dalian 116023,China;2.Key Laboratory of Tropical&Subtropical Fishery Resource Application&Cultivation,Ministry of Agriculture,Pearl River Fisheries Research Institute,Chinese Academy of Fishery Sciences,Guangzhou 510380,China)

    Screening and genotype SNPs in 3'UTR of Apoprotein A-I-1 gene were conducted in grass carp(Ctenopharyngodon idellus)through PCR sequencing and PCR-RFLP method from EST database.Two SNPs(C792T and G851A)were obtained in the cultured population of grass carp from Pearl River.The CC genotype in C792T site was found 46.53%,CT genotype 50.69%,and TT genotype 2.78%.The GG genotype in G851A was found 77.08%,GA genotype 20.83%,and AA genotype 2.08%.A general linear model was used to analyze associations of apoA-I-1 gene with 4 major growth traits(body weight,body length,body height and head length).The results showed that the mean body weight of BB genotype was higher than that of CC and CT genotype in C792T site (P>0.05).The mean weight,body length,body height and head length of GG genotype in G851A was higher than GA and AA genotype(P<0.05).Six diplotypes were constructed based on two SNPs in the experimental population.The mean body weight of diplotype D3(TTC792T-GGG851A)was higher than other diplotypes.The diplotype D6 (CCC792T-AAG851A)had the minimal mean body weight,and had significant differences in body length,body height and head length(P<0.05).TT genotype in C792T site was advantage genotype.GG genotype in G851A site was advantage genotype,and AA genotype disadvantage genotype.It is suggested that apoA-I-1 gene be a candidate modifier gene for grass carp growth traits,indicating that the apoA-I-1 gene could be useful for grass carp molecular breeding.

    grass carp;SNP;Apoprotein A-I-1;apoA-I-1;growth trait;diplotype;association analysis

    Q789

    A

    2095-1388(2012)01-0012-06

    2011-04-06

    現(xiàn)代農(nóng)業(yè)產(chǎn)業(yè)技術(shù)體系建設(shè)專(zhuān)項(xiàng)資金資助項(xiàng)目(CARS-46-03);國(guó)家高技術(shù)研究發(fā)展計(jì)劃項(xiàng)目(2011AA100403);國(guó)家科技支撐計(jì)劃項(xiàng)目(2012BAD26B02)

    劉小獻(xiàn)(1986-),女,碩士研究生。E-mail:littlebearliu@126.com

    白俊杰(1957-),男,研究員。E-mail:jjbai@163.net

    猜你喜歡
    生長(zhǎng)
    野蠻生長(zhǎng)
    碗蓮生長(zhǎng)記
    小讀者(2021年2期)2021-03-29 05:03:48
    生長(zhǎng)的樹(shù)
    自由生長(zhǎng)的家
    美是不斷生長(zhǎng)的
    快速生長(zhǎng)劑
    共享出行不再“野蠻生長(zhǎng)”
    生長(zhǎng)在哪里的啟示
    野蠻生長(zhǎng)
    NBA特刊(2018年21期)2018-11-24 02:48:04
    生長(zhǎng)
    文苑(2018年22期)2018-11-19 02:54:14
    欧美日韩视频高清一区二区三区二| 一本大道久久a久久精品| 黄色配什么色好看| 久久久a久久爽久久v久久| 丰满饥渴人妻一区二区三| 老司机影院成人| 色婷婷久久久亚洲欧美| 免费观看性生交大片5| 亚洲av.av天堂| 免费在线观看黄色视频的| 亚洲色图 男人天堂 中文字幕 | 999精品在线视频| 99久国产av精品国产电影| 美女内射精品一级片tv| 国产1区2区3区精品| 亚洲内射少妇av| 久久久久视频综合| 你懂的网址亚洲精品在线观看| 免费看光身美女| 人妻一区二区av| 人妻 亚洲 视频| 18+在线观看网站| 人妻一区二区av| 亚洲,欧美精品.| 精品人妻在线不人妻| 亚洲性久久影院| 午夜日本视频在线| 夫妻午夜视频| 人妻少妇偷人精品九色| 一区二区日韩欧美中文字幕 | 国产精品久久久久久久电影| 久久久久国产网址| 另类亚洲欧美激情| 亚洲av在线观看美女高潮| 国产亚洲av片在线观看秒播厂| videos熟女内射| 国产免费又黄又爽又色| 免费少妇av软件| 日韩欧美一区视频在线观看| 成年女人在线观看亚洲视频| 午夜视频国产福利| 亚洲国产精品999| 黄色一级大片看看| 成人免费观看视频高清| 免费女性裸体啪啪无遮挡网站| 在线观看免费日韩欧美大片| 精品一区在线观看国产| 一边摸一边做爽爽视频免费| 久久精品熟女亚洲av麻豆精品| 亚洲第一区二区三区不卡| 精品久久久久久电影网| 一本久久精品| 国产色婷婷99| 高清黄色对白视频在线免费看| 少妇人妻 视频| 亚洲综合色惰| 99久国产av精品国产电影| 十八禁高潮呻吟视频| 久久久久久久久久久免费av| 视频区图区小说| 高清不卡的av网站| 在线亚洲精品国产二区图片欧美| 国产探花极品一区二区| 久久精品久久久久久噜噜老黄| 免费观看性生交大片5| 热re99久久精品国产66热6| 少妇的逼水好多| 人成视频在线观看免费观看| 午夜激情久久久久久久| 久久人人爽人人片av| av视频免费观看在线观看| av网站免费在线观看视频| 亚洲综合色网址| 中文字幕精品免费在线观看视频 | 少妇的丰满在线观看| av电影中文网址| 国产亚洲av片在线观看秒播厂| 大陆偷拍与自拍| 国产精品秋霞免费鲁丝片| 一二三四中文在线观看免费高清| 久久久久久久国产电影| 黑丝袜美女国产一区| 搡女人真爽免费视频火全软件| 久久97久久精品| 欧美成人精品欧美一级黄| 久久亚洲国产成人精品v| 日韩av不卡免费在线播放| 有码 亚洲区| 女性被躁到高潮视频| 久热久热在线精品观看| 成年女人在线观看亚洲视频| 91精品三级在线观看| 老司机影院毛片| 久久99一区二区三区| 18在线观看网站| 免费观看a级毛片全部| 一区二区av电影网| 夫妻午夜视频| 一本久久精品| 欧美成人午夜精品| 少妇精品久久久久久久| a级片在线免费高清观看视频| 欧美日本中文国产一区发布| 免费女性裸体啪啪无遮挡网站| 国产黄色免费在线视频| 大码成人一级视频| 视频中文字幕在线观看| 欧美日韩av久久| 赤兔流量卡办理| 国产极品天堂在线| 亚洲成人手机| 亚洲国产看品久久| 汤姆久久久久久久影院中文字幕| 国产综合精华液| 欧美成人精品欧美一级黄| 日韩熟女老妇一区二区性免费视频| 亚洲国产av新网站| 在线 av 中文字幕| 国产又色又爽无遮挡免| 国产极品粉嫩免费观看在线| 丰满少妇做爰视频| 人妻 亚洲 视频| 91国产中文字幕| 久久久久久久精品精品| 久久久国产欧美日韩av| 亚洲av国产av综合av卡| 国产1区2区3区精品| 日日啪夜夜爽| 边亲边吃奶的免费视频| 亚洲 欧美一区二区三区| 国产在线一区二区三区精| 狠狠婷婷综合久久久久久88av| 亚洲国产欧美在线一区| 国产精品99久久99久久久不卡 | 日韩视频在线欧美| 免费高清在线观看视频在线观看| 99热网站在线观看| 男人爽女人下面视频在线观看| 精品人妻在线不人妻| 免费人妻精品一区二区三区视频| 久久精品久久精品一区二区三区| av黄色大香蕉| 丝袜人妻中文字幕| 国产高清国产精品国产三级| 最近中文字幕高清免费大全6| 久久久久精品性色| 女的被弄到高潮叫床怎么办| 天天影视国产精品| 久久久久精品人妻al黑| 婷婷色麻豆天堂久久| 久久久久视频综合| 精品卡一卡二卡四卡免费| 黄色 视频免费看| 久久毛片免费看一区二区三区| 最近中文字幕2019免费版| 日本与韩国留学比较| 少妇的丰满在线观看| 99久久综合免费| 国产成人欧美| 18禁国产床啪视频网站| 国产成人精品无人区| 桃花免费在线播放| 国精品久久久久久国模美| 黄网站色视频无遮挡免费观看| 精品视频人人做人人爽| 中文字幕人妻丝袜制服| 久久久久精品久久久久真实原创| 人人澡人人妻人| 少妇 在线观看| 国产xxxxx性猛交| 美女脱内裤让男人舔精品视频| 国产熟女午夜一区二区三区| 在线观看国产h片| 一级黄片播放器| 乱人伦中国视频| 日韩,欧美,国产一区二区三区| 亚洲一区二区三区欧美精品| 啦啦啦在线观看免费高清www| 久久午夜福利片| 黄色 视频免费看| 人人妻人人爽人人添夜夜欢视频| 亚洲激情五月婷婷啪啪| 国产男女内射视频| 我要看黄色一级片免费的| 十八禁高潮呻吟视频| 成人亚洲精品一区在线观看| 婷婷色综合大香蕉| 综合色丁香网| 久久久精品区二区三区| kizo精华| 亚洲精品久久久久久婷婷小说| 午夜激情久久久久久久| 99久国产av精品国产电影| 在线观看人妻少妇| 国产精品一区二区在线观看99| av福利片在线| 日韩成人av中文字幕在线观看| 韩国高清视频一区二区三区| 精品人妻在线不人妻| 国产高清不卡午夜福利| 又大又黄又爽视频免费| 日本爱情动作片www.在线观看| 国产无遮挡羞羞视频在线观看| 18禁国产床啪视频网站| 午夜日本视频在线| 97在线视频观看| www日本在线高清视频| 新久久久久国产一级毛片| 久久久久久久久久久免费av| 99香蕉大伊视频| 久久青草综合色| 超碰97精品在线观看| 宅男免费午夜| 亚洲精品乱码久久久久久按摩| 国产激情久久老熟女| 国产老妇伦熟女老妇高清| 又黄又爽又刺激的免费视频.| 日本免费在线观看一区| 大码成人一级视频| 亚洲欧美一区二区三区国产| 国产欧美日韩综合在线一区二区| xxx大片免费视频| 色婷婷久久久亚洲欧美| 一区二区日韩欧美中文字幕 | 亚洲精品色激情综合| www日本在线高清视频| 成年女人在线观看亚洲视频| 国产成人精品无人区| 男女边吃奶边做爰视频| 久久国产亚洲av麻豆专区| 高清毛片免费看| 国产不卡av网站在线观看| 亚洲天堂av无毛| 最后的刺客免费高清国语| 青春草视频在线免费观看| 久久免费观看电影| 最近中文字幕高清免费大全6| 国产日韩欧美亚洲二区| 成人亚洲欧美一区二区av| 又粗又硬又长又爽又黄的视频| 狂野欧美激情性bbbbbb| 免费大片黄手机在线观看| av免费在线看不卡| 狂野欧美激情性bbbbbb| 亚洲精品美女久久久久99蜜臀 | 国产精品久久久久成人av| 在线观看免费视频网站a站| 2022亚洲国产成人精品| 欧美性感艳星| 天天躁夜夜躁狠狠躁躁| 一级毛片黄色毛片免费观看视频| 一级片'在线观看视频| 国产精品一二三区在线看| 久久鲁丝午夜福利片| freevideosex欧美| 卡戴珊不雅视频在线播放| 美女视频免费永久观看网站| 亚洲第一av免费看| 一级,二级,三级黄色视频| 黄色怎么调成土黄色| 韩国精品一区二区三区 | 高清不卡的av网站| 在线亚洲精品国产二区图片欧美| 免费高清在线观看日韩| 免费久久久久久久精品成人欧美视频 | 一级爰片在线观看| 精品熟女少妇av免费看| 一边亲一边摸免费视频| 国产精品久久久久久精品古装| 国产精品.久久久| av免费在线看不卡| 久久午夜福利片| 香蕉精品网在线| 久久精品国产亚洲av涩爱| 汤姆久久久久久久影院中文字幕| 两个人看的免费小视频| 国产亚洲一区二区精品| 大香蕉久久网| 国产欧美另类精品又又久久亚洲欧美| 99热6这里只有精品| 九草在线视频观看| 七月丁香在线播放| 高清欧美精品videossex| 亚洲精品久久成人aⅴ小说| 久久99精品国语久久久| 久久精品久久精品一区二区三区| 丰满饥渴人妻一区二区三| 伦理电影大哥的女人| 一区在线观看完整版| 黄色怎么调成土黄色| 午夜福利视频在线观看免费| 人体艺术视频欧美日本| 亚洲av.av天堂| 最后的刺客免费高清国语| 久久精品aⅴ一区二区三区四区 | 久久久久网色| 十八禁网站网址无遮挡| 一本色道久久久久久精品综合| 下体分泌物呈黄色| 欧美老熟妇乱子伦牲交| 国产极品粉嫩免费观看在线| 在线 av 中文字幕| 少妇被粗大的猛进出69影院 | √禁漫天堂资源中文www| 亚洲五月色婷婷综合| 在线观看免费视频网站a站| 免费人妻精品一区二区三区视频| 丰满迷人的少妇在线观看| 久久久精品94久久精品| 欧美亚洲日本最大视频资源| 在线天堂中文资源库| 精品一区二区三区四区五区乱码 | 春色校园在线视频观看| 国产又爽黄色视频| 2022亚洲国产成人精品| 高清不卡的av网站| 乱码一卡2卡4卡精品| 欧美精品一区二区大全| 亚洲国产精品一区三区| 香蕉精品网在线| 在线天堂中文资源库| 午夜av观看不卡| 久久97久久精品| 高清在线视频一区二区三区| 一级,二级,三级黄色视频| av卡一久久| 又黄又爽又刺激的免费视频.| 亚洲精品日韩在线中文字幕| 卡戴珊不雅视频在线播放| h视频一区二区三区| 只有这里有精品99| 最近中文字幕高清免费大全6| 久久人人97超碰香蕉20202| 肉色欧美久久久久久久蜜桃| 国产精品免费大片| 纯流量卡能插随身wifi吗| 国产成人精品无人区| 午夜免费鲁丝| 精品亚洲成a人片在线观看| av.在线天堂| 两个人看的免费小视频| 精品亚洲乱码少妇综合久久| 最新的欧美精品一区二区| 久久午夜福利片| 97在线人人人人妻| 91精品三级在线观看| 母亲3免费完整高清在线观看 | 啦啦啦啦在线视频资源| 男人添女人高潮全过程视频| 精品人妻熟女毛片av久久网站| 街头女战士在线观看网站| 欧美国产精品一级二级三级| 国产精品一区www在线观看| 高清在线视频一区二区三区| 国产精品不卡视频一区二区| 国产福利在线免费观看视频| 少妇被粗大猛烈的视频| 观看美女的网站| 午夜精品国产一区二区电影| 亚洲成av片中文字幕在线观看 | 在线看a的网站| 看免费av毛片| 精品视频人人做人人爽| 人人妻人人澡人人看| 国产精品久久久久久精品电影小说| 精品一区二区三卡| 永久免费av网站大全| 亚洲国产av影院在线观看| 午夜91福利影院| 国产色婷婷99| 伊人亚洲综合成人网| 一级毛片电影观看| 国产精品一二三区在线看| 国产在线视频一区二区| 亚洲av成人精品一二三区| 麻豆乱淫一区二区| 人人妻人人澡人人爽人人夜夜| 欧美少妇被猛烈插入视频| 男人舔女人的私密视频| 视频在线观看一区二区三区| 亚洲 欧美一区二区三区| 精品亚洲乱码少妇综合久久| 午夜福利网站1000一区二区三区| 观看美女的网站| 97在线人人人人妻| 国产精品无大码| 国产老妇伦熟女老妇高清| 日韩,欧美,国产一区二区三区| 欧美精品国产亚洲| 欧美成人精品欧美一级黄| 交换朋友夫妻互换小说| 美女xxoo啪啪120秒动态图| 9191精品国产免费久久| 妹子高潮喷水视频| 中国美白少妇内射xxxbb| 91午夜精品亚洲一区二区三区| 亚洲综合色网址| 成年人午夜在线观看视频| 国产成人精品一,二区| av在线观看视频网站免费| 久久ye,这里只有精品| 麻豆乱淫一区二区| 深夜精品福利| 69精品国产乱码久久久| 91在线精品国自产拍蜜月| 女的被弄到高潮叫床怎么办| 大香蕉久久网| 菩萨蛮人人尽说江南好唐韦庄| 热re99久久国产66热| 伦精品一区二区三区| 欧美少妇被猛烈插入视频| 夜夜爽夜夜爽视频| 国产精品国产三级国产av玫瑰| 天天躁夜夜躁狠狠躁躁| 欧美日韩精品成人综合77777| 亚洲精品成人av观看孕妇| 天天操日日干夜夜撸| 在线观看免费高清a一片| 免费在线观看完整版高清| 欧美另类一区| 日韩视频在线欧美| 一级毛片 在线播放| 热99国产精品久久久久久7| 国产免费现黄频在线看| 人人澡人人妻人| 久久精品aⅴ一区二区三区四区 | 全区人妻精品视频| 久久99热这里只频精品6学生| 午夜免费观看性视频| 18禁裸乳无遮挡动漫免费视频| 18禁观看日本| 国产精品偷伦视频观看了| 少妇人妻精品综合一区二区| 免费观看在线日韩| 精品少妇内射三级| 午夜视频国产福利| 亚洲成国产人片在线观看| 香蕉精品网在线| 欧美日韩视频精品一区| 韩国av在线不卡| 亚洲av国产av综合av卡| 水蜜桃什么品种好| 美女大奶头黄色视频| 我要看黄色一级片免费的| 老司机影院毛片| 亚洲成人手机| 亚洲性久久影院| 伦理电影免费视频| 亚洲三级黄色毛片| av又黄又爽大尺度在线免费看| 熟妇人妻不卡中文字幕| 国产一区二区在线观看日韩| 少妇人妻精品综合一区二区| 日韩av不卡免费在线播放| 熟女人妻精品中文字幕| 寂寞人妻少妇视频99o| av天堂久久9| 18禁在线无遮挡免费观看视频| 国产精品不卡视频一区二区| 午夜久久久在线观看| 又黄又爽又刺激的免费视频.| 国产一区二区在线观看av| 成人二区视频| 国产精品成人在线| 亚洲欧美中文字幕日韩二区| 亚洲第一区二区三区不卡| 中文字幕精品免费在线观看视频 | 国产成人免费无遮挡视频| 国产欧美亚洲国产| 一二三四中文在线观看免费高清| 久久精品国产亚洲av涩爱| 免费大片18禁| 蜜桃在线观看..| 欧美性感艳星| 少妇被粗大的猛进出69影院 | 国产精品99久久99久久久不卡 | 夜夜骑夜夜射夜夜干| kizo精华| 欧美日韩精品成人综合77777| 国产免费一区二区三区四区乱码| videos熟女内射| 边亲边吃奶的免费视频| 伦理电影免费视频| 少妇人妻精品综合一区二区| 黄色一级大片看看| 国产片内射在线| 久久精品国产a三级三级三级| 纵有疾风起免费观看全集完整版| 精品卡一卡二卡四卡免费| 欧美日本中文国产一区发布| 午夜激情久久久久久久| 国产无遮挡羞羞视频在线观看| 性高湖久久久久久久久免费观看| 久久热在线av| 精品一区二区三卡| xxx大片免费视频| 国产精品久久久av美女十八| 各种免费的搞黄视频| 极品少妇高潮喷水抽搐| 国产精品无大码| 九色亚洲精品在线播放| 国产福利在线免费观看视频| 亚洲av免费高清在线观看| 在线看a的网站| 黄色视频在线播放观看不卡| 肉色欧美久久久久久久蜜桃| 一二三四中文在线观看免费高清| 伦理电影大哥的女人| 大片免费播放器 马上看| 亚洲成色77777| 一级毛片我不卡| 国产av码专区亚洲av| 国语对白做爰xxxⅹ性视频网站| 香蕉丝袜av| 国产老妇伦熟女老妇高清| 国产免费一级a男人的天堂| 国产成人aa在线观看| 久久99精品国语久久久| 久久久久精品性色| 亚洲精品自拍成人| 亚洲国产色片| 日韩人妻精品一区2区三区| 日韩 亚洲 欧美在线| 伊人亚洲综合成人网| 久久人人爽人人片av| 成人手机av| 亚洲五月色婷婷综合| av在线观看视频网站免费| 性色av一级| 国产黄频视频在线观看| 欧美3d第一页| 国产精品一国产av| 国产亚洲一区二区精品| av免费观看日本| 麻豆乱淫一区二区| 汤姆久久久久久久影院中文字幕| 热re99久久精品国产66热6| 少妇人妻 视频| 少妇精品久久久久久久| 国产激情久久老熟女| 国产又色又爽无遮挡免| 黄片播放在线免费| 欧美激情极品国产一区二区三区 | 久久免费观看电影| 侵犯人妻中文字幕一二三四区| 久久99一区二区三区| 国产精品99久久99久久久不卡 | 亚洲成人一二三区av| 午夜av观看不卡| 晚上一个人看的免费电影| 精品久久久久久电影网| 欧美精品高潮呻吟av久久| 亚洲av男天堂| 亚洲国产毛片av蜜桃av| 极品人妻少妇av视频| 国产乱人偷精品视频| 国产精品国产三级国产av玫瑰| av黄色大香蕉| 国产一区亚洲一区在线观看| 另类精品久久| av免费观看日本| 丝瓜视频免费看黄片| 我的女老师完整版在线观看| 久久国产精品男人的天堂亚洲 | 欧美日韩亚洲高清精品| 如日韩欧美国产精品一区二区三区| 国产精品一二三区在线看| 18在线观看网站| 激情视频va一区二区三区| av福利片在线| 少妇精品久久久久久久| 国产日韩欧美亚洲二区| 欧美日韩成人在线一区二区| 熟妇人妻不卡中文字幕| 亚洲一级一片aⅴ在线观看| 成年女人在线观看亚洲视频| 丰满乱子伦码专区| 亚洲成av片中文字幕在线观看 | 久久精品久久久久久噜噜老黄| 国产高清三级在线| 丰满乱子伦码专区| 国产午夜精品一二区理论片| 国产综合精华液| 日本黄色日本黄色录像| 久久精品久久精品一区二区三区| 狠狠婷婷综合久久久久久88av| 国产日韩欧美视频二区| 日本vs欧美在线观看视频| 久久婷婷青草| 国产高清国产精品国产三级| 国产一区二区三区综合在线观看 | 国产伦理片在线播放av一区| 中文字幕人妻丝袜制服| 男人爽女人下面视频在线观看| 日韩av不卡免费在线播放| 亚洲av综合色区一区| 精品国产一区二区久久| 国产精品麻豆人妻色哟哟久久| 亚洲内射少妇av| 欧美激情 高清一区二区三区| 人人澡人人妻人| 亚洲久久久国产精品| 亚洲成人手机| 精品熟女少妇av免费看| 熟女人妻精品中文字幕| 国产一区二区三区av在线| av线在线观看网站| 久久久国产一区二区| 国产又色又爽无遮挡免| 日日摸夜夜添夜夜爱| 欧美激情极品国产一区二区三区 | 国产精品久久久av美女十八| xxx大片免费视频| 热99久久久久精品小说推荐| xxxhd国产人妻xxx|