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

    翹嘴鱖微衛(wèi)星標記及其與主要經濟性狀的相關分析

    2017-03-17 06:25:38孫際佳李桂峰劉麗王海芳胥鵬
    水產學雜志 2017年1期
    關鍵詞:微衛(wèi)星體長基因型

    孫際佳,李桂峰,劉麗,王海芳,胥鵬

    (1.華南農業(yè)大學海洋學院,廣東 廣州 510642;2.中山大學生命科學學院,水生經濟動物研究所,南海生物資源開發(fā)與利用協(xié)同創(chuàng)新中心,廣東 廣州 510275;3.佛山市南海百容水產良種有限公司,廣東省工程技術研究中心,廣東 佛山 528216)

    翹嘴鱖微衛(wèi)星標記及其與主要經濟性狀的相關分析

    孫際佳1,李桂峰2,3,劉麗1,王海芳2,胥鵬2,3

    (1.華南農業(yè)大學海洋學院,廣東 廣州 510642;2.中山大學生命科學學院,水生經濟動物研究所,南海生物資源開發(fā)與利用協(xié)同創(chuàng)新中心,廣東 廣州 510275;3.佛山市南海百容水產良種有限公司,廣東省工程技術研究中心,廣東 佛山 528216)

    利用篩選出的7對具有較高多態(tài)性的微衛(wèi)星標記,檢測106尾翹嘴鱖Siniperca chuatsi選育個體的基因組DNA,分析這些微衛(wèi)星標記與翹嘴鱖體長、體質量和體高的相關性。結果獲得67個等位基因,各位點的等位基因數(shù)為3~19,片段大小在147~530bp之間;期望雜合度0.5092~0.9207,均值為0.7574;各位點的多態(tài)信息含量在0.4639~0.9101之間,均值為0.7197,表明所選擇的SSR標記識別力較高,適用于翹嘴鱖選育群體遺傳分析和標記輔助育種研究。相關性分析結果表明,G4位點中含有的片段為219bp等位基因的基因型(229/219或219/219)個體的體質量、體長和體高的表型效應顯著高于其他的基因型,G10位點中246/246基因型的體質量和體高表型效應顯著高于其他基因型,可作為未來翹嘴鱖分子標記輔助育種的重要參考位點。

    翹嘴鱖;微衛(wèi)星標記;經濟性狀;相關性分析

    在分布于我國的2屬9種鱖類中,翹嘴鱖Siniperca chuatsi個體大、生長快、肉質鮮嫩而具有很高的經濟價值,是我國養(yǎng)殖的主要經濟鱖類。近年來,隨著翹嘴鱖養(yǎng)殖規(guī)模的擴大,其種質質量與生產要求之間的矛盾,如生長速度降低、種質退化和病害頻發(fā)等,嚴重制約了翹嘴鱖養(yǎng)殖業(yè)的持續(xù)健康發(fā)展[1-3]。因此,急需通過選擇育種來改善翹嘴鱖的種質質量,以滿足鱖養(yǎng)殖業(yè)的需求。

    近年來,隨著基于DNA多態(tài)性的分子標記技術的迅速發(fā)展和廣泛應用,分子標記輔助選擇(Marker Assisted Selection,MAS)給魚類育種提供了一種更有效的途徑[4]。分子標記輔助選擇育種的目標是尋找到與某個或某些數(shù)量性狀基因座位(Quantitative Trait Locus,QTL)連鎖的分子標記,并在育種過程中通過這些連鎖遺傳標記實現(xiàn)對目標數(shù)量性狀基因型的直接選擇,從而縮短育種周期,加快育種進程[5]。目前,在鯉Cyprinus carpio[6]、虹鱒Oncorhynchus mykiss[7]、大菱鲆Scophthalmus maximus[8]、美洲紅點鮭Salvelinus fontinalis[9]和紅鰭東方鲀Takifugu rubripes[10]等魚類中已有不少與性狀連鎖的分子標記的報道,而在鱖類研究中與性狀連鎖的分子標記鮮有報道。Wang等[11]利用SNP(Single Nucleotide Polymorphism)分子標記篩選到IGF-I基因的兩個SNP位點(g.A321G和g.A537G)與翹嘴鱖♀×斑鱖♂雜交F1的體質量和體寬生長性狀相關,但至今未見微衛(wèi)星標記與翹嘴鱖生長性狀間相關分析的報道。

    微衛(wèi)星DNA又稱簡單序列重復(Simple Sequence Repeat,SSR),是由Miesfeld等[12]首次發(fā)現(xiàn)的一種簡單串聯(lián)重復DNA序列,因具有數(shù)目多、分布廣、多態(tài)性高、信息含量豐富、突變快、雜合度高、共顯性遺傳和檢測方法簡便等優(yōu)點[13],成為第二代DNA遺傳標記的核心,已廣泛應用于QTL定位和分子標記輔助育種等方面[14,15]。近年來,微衛(wèi)星標記技術已廣泛應用于鱖選育群體遺傳多樣性分析[16,17]。

    本研究在已經獲得的翹嘴鱖微衛(wèi)星標記的基礎上,分析微衛(wèi)星標記與翹嘴鱖主要經濟性狀(體質量、體長和體高)之間的相關性,期望發(fā)現(xiàn)與主要經濟性狀相關聯(lián)的微衛(wèi)星標記,為翹嘴鱖的分子標記輔助選擇育種研究和實驗提供基礎參考數(shù)據(jù)。

    1 材料與方法

    1.1 材料

    本研究所用的野生翹嘴鱖親魚分別引自湖南省和江蘇省,按1∶1親本比例分別培育3個湖南翹嘴鱖家系、1個江蘇翹嘴鱖家系、以及2個湖南翹嘴鱖和江蘇翹嘴鱖雜交家系,構建由6個家系組成的翹嘴鱖選育群體。將同期孵出的6組家系各取1 000尾,共6 000尾全長約1cm的魚苗混合養(yǎng)殖在同一池塘中,以消除環(huán)境和飼料的影響。保證餌料魚充足的條件下,飼養(yǎng)至8月齡后,隨機選取106尾進行表型和基因型分析。

    1.2 微衛(wèi)星引物

    參考匡剛橋[18]分離的翹嘴鱖微衛(wèi)星序列,用軟件Primer Premier 5.0設計微衛(wèi)星引物。在設計的20對引物中,根據(jù)擴增穩(wěn)定、6%變性PAGE條帶清晰、多態(tài)性較高的選取原則,篩選出7對SSR引物(表1)分析群體遺傳多樣性和主要經濟性狀的相關性。所設計的引物由Invitrogen生物技術公司進行合成。

    表1 篩選的7對微衛(wèi)星引物的核心序列、退火溫度及其引物序列Tab.1 Core repeats,primer sequences and annealing temperatures of the 7 microsatellite DNA loci

    1.3 方法

    1.3.1 表型性狀的度量

    根據(jù)魚類學方法測定和計算翹嘴鱖的體質量、體長和體高等表型性狀。

    1.3.2 基因組DNA的提取及SSR分析

    提取106尾樣本魚的背部肌肉用于基因組DNA分析。魚類基因組DNA抽提方法參照傳統(tǒng)的苯酚-氯仿法基本原理,并根據(jù)本實驗實際情況對離心時間和抽提次數(shù)稍作更改。

    PCR反應體系總量25μL,包括2×PCR Reaction Mix12.5μL、上下游引物(10μmol/L)各1μL、模板DNA(30ng/μL)2μL、Taq DNA聚合酶(2.5U/μL)0.25μL和ddH2O8.25μL。擴增反應使用Bio-Rad公司的MyCycler型PCR儀,采用梯度PCR法篩選最佳退火溫度。PCR反應程序為:94℃預變性5min,94℃變性30s,退火45s(退火溫度依引物而異,見表1),72℃延伸1min,經35個循環(huán)后,72℃延伸10min,4℃保存。6%變性聚丙烯酰胺凝膠電泳,硝酸銀染色,掃描儀記錄電泳圖譜,用Gel-Proanalyzer32軟件手動分析SSR片段的大小。

    1.4 統(tǒng)計分析

    利用SPSS 17.0軟件統(tǒng)計體質量、體長、體高表型性狀的變異,計算各性狀間的相關系數(shù)。采用SAS 9.2軟件進行Shapiro-Wilk檢驗表型性狀的頻率分布是否顯著偏離正態(tài)分布。利用POPGENE version 1.32和PIC-CALC軟件統(tǒng)計各位點的等位基因數(shù)目、基因型、期望雜合度(He)和多態(tài)信息含量 PIC(Polymorphism Information Content)。利用SPSS17.0軟件對微衛(wèi)星標記與體重、體長和體高性狀進行單因素方差分析。利用SAS 9.2軟件進行基因型與性狀間的多重比較:對符合正態(tài)分布的性狀采用參數(shù)檢驗的Tukey-Cramer法進行比較,不符合的性狀則采用非參數(shù)檢驗的秩和檢驗法進行比較。由于部分基因型在微衛(wèi)星位點中的頻率太低,缺乏統(tǒng)計分析意義,因此在實際分析中僅考慮至少有2次觀測值的基因型。

    2 結果和分析

    2.1 表型性狀的統(tǒng)計分析

    翹嘴鱖選育群體的體質量范圍在88.4~648.1g,體長范圍在14.40~29.40cm,體高范圍在5.32~11.30 cm。體質量、體長和體高的變異系數(shù)分別為44.59%、13.86%和18.90%,體質量性狀的變異最大(表2)。

    表2 翹嘴鱖的體質量、體長和體高的表型統(tǒng)計量Tab.2 The phenotypic statistics of body weight,body length and body height in mandarin fish Siniperca chuatsi

    體質量、體長和體高的正態(tài)性檢驗Q-Q(Quantile-Quantile)圖如圖1所示。圖1中給出的體質量、體長和體高性狀的Shapiro-Wilk正態(tài)分布檢驗結果表明,體長符合正態(tài)分布(P=0.3930>0.05),而體質量和體高不符合正態(tài)分布(P<0.0001)。因此,基因型與體長的多重比較應采用Tukey-Kramer法檢驗,而基因型與體質量、體高的多重比較應采用秩和檢驗法檢驗。

    2.2 體質量、體長和體高的Pearson相關分析

    圖1 正態(tài)性檢驗的Q-Q圖(體質量(a),體長(b),體高(c))Fig.1 Normal quantile-quantile plots of body weight(a),body length(b)and body height(c)

    利用Pearson相關性檢驗方法分析3個性狀之間的相關程度(表3)。由表3可知,體質量、體長和體高的Pearson相關系數(shù)均達0.9以上,且相關極顯著(P<0.01)。

    表3 翹嘴鱖的體質量、體長和體高的Pearson相關系數(shù)Tab.3 Pearson correlation coefficient of body weight, body length and body height in mandarin fish Siniperca chuatsi

    2.3 電泳結果

    用7對微衛(wèi)星引物對翹嘴鱖基因組DNA進行PCR擴增和6%變性聚丙烯酰胺凝膠電泳檢測,各微衛(wèi)星引物均獲得了清晰、穩(wěn)定的DNA條帶,并在個體間表現(xiàn)出不同程度的多態(tài)性。

    2.4 微衛(wèi)星多態(tài)性檢測

    本研究中翹嘴鱖群體的遺傳多樣性信息及性狀F值見表4。在7對微衛(wèi)星標記中共檢測到67個等位基因,平均等位基因數(shù)為9.57個,等位基因片段大小在147~530bp之間。7個位點共檢測到基因型61個,平均每個位點8.71個;期望雜合度He變化在0.5092~0.9207之間,平均0.7574;多態(tài)信息含量PIC在0.4639~0.9101之間,平均0.7197。

    2.5 微衛(wèi)星位點與體質量、體長和體高的相關分析

    采用單因素方差分析方法檢驗了微衛(wèi)星位點與翹嘴鱖體質量、體長、體高性狀連鎖的顯著性。各微衛(wèi)星位點對體質量、體長和體高影響的F檢驗值見表4。7個微衛(wèi)星位點中有4個位點(G4、G10、G12和G13)與體質量、體長和體高均顯著相關(P<0.01)。對上述4個差異顯著的標記進行了不同基因型間不同性狀的多重比較(體長用Tukey-Kramer法,體質量和體高用秩和檢驗法),結果見表5。

    在標記G4中共檢測到4個基因型,其中基因型219/219的12個個體和基因型229/219的18個個體的平均體質量、體長和體高值均顯著高于其他基因型個體的平均值(P<0.05),且基因型219/219與基因型229/219個體的平均體質量、體長和體高值均差異不顯著(P>0.05),說明G4位點中219bp片段的等位基因可能對體質量、體長和體高性狀起正面影響。

    在標記G10中共檢測到4個基因型?;蛐?46/246的39個個體的體質量和體高均顯著高于其他基因型個體的平均值(P<0.05),說明該基因型與體質量和體高正相關(圖2)。基因型248/248個體的平均體高顯著小于基因型255/248和255/255的個體(P<0.05),因而推測248bp片段的等位基因可能對體高性狀起負面影響。

    在標記G12中共檢測到13個基因型?;蛐?29/279的2個個體和359/305的9個個體的體質量和體高顯著高于其他基因型個體(P<0.05),但基因型329/279與基因型359/305個體的平均體質量、體長和體高差異不顯著(P>0.05)。

    在標記G13中也檢測到13個基因型?;蛐?75/237的2個個體和304/268的9個個體的體質量和體高顯著高于其他基因型個體(P<0.05),但基因型275/237與基因型304/268個體的平均體質量、體長和體高差異均不顯著(P>0.05)。

    表4 所研究翹嘴鱖群體的遺傳信息及性狀F值Tab.4 Genetic information and F values in populations of mandarin fish Siniperca chuatsi

    表5 翹嘴鱖群體的4個微衛(wèi)星位點不同基因型體質量、體長和體高的平均值及多重比較Tab.5 Mean and multiple comparisons of bod weight,body length and body height in 4 microsatellite loci in mandarin fish Siniperca chuatsi

    3 討論

    3.1 分子標記輔助育種SSR標記的篩選

    動物遺傳育種的分子標記中,SSR標記已成為分析群體遺傳結構和與重要經濟性狀的遺傳連鎖關系最理想的標記之一,并已經獲得了廣泛的應用[19]。有關翹嘴鱖SSR標記的開發(fā)已有一些研究。Zhang等[20]、匡剛橋等[21]和Liu等[22]利用FIASCO(Fast Isolation by AFLP Sequences Containing Repeats)法分別篩選出18、18和40個具有多態(tài)性的SSR標記。Qu等2012年[23]和2013年[24]分析了翹嘴鱖♀×斑鱖♂雜交子代的轉錄組,分別篩選到適用于翹嘴鱖遺傳分析的43和29個多態(tài)SSR位點。但是,至今未見有應用SSR標記進行翹嘴鱖分子標記輔助育種的報道。本研究利用匡剛橋[18]分離的20個翹嘴鱖SSR序列,對選育群體進行擴增,篩選出7個具有多態(tài)性的4核苷酸重復的SSR標記。已有研究表明,核心序列過長的雙核苷酸重復的SSR標記易產生影子帶,對結果的判定有一定的影響,而3、4核苷酸重復的SSR標記具有多態(tài)性高,產生影子帶的幾率較低,結果判定更準確等優(yōu)點[25,26]。本研究篩選的SSR標記擴增條帶穩(wěn)定、清晰,PIC結果顯示,除G4位點表現(xiàn)為中度多態(tài)外(0.25<PIC<0.5),其余6個位點均為高度多態(tài)(PIC>0.5),表明這7個SSR標記具有較高的識別力,適用于翹嘴鱖群體遺傳分析和標記輔助育種的研究。

    3.2 選育群體的遺傳多樣性

    在分子標記輔助育種的選育體系中,維系選育群體的遺傳多樣性是選育體系能否可持續(xù)性發(fā)展的關鍵。如果選育群體的遺傳多樣性持續(xù)降低,極可能導致選育群體陷入近交衰退和較差選擇響應的惡性循環(huán)[27]。方展強[1]等利用RAPD(Randomly Amplified Polymorphic DNA)技術、楊宇暉等[28]利用線粒體控制區(qū)(D-loop)核苷酸序列、吳旭等[29]和梅秋蘭等[30]利用SSR標記技術對翹嘴鱖群體遺傳多樣性的研究結果均表明:翹嘴鱖養(yǎng)殖群體的遺傳多樣性水平顯著低于野生或原種群體。因此,開展翹嘴鱖的遺傳育種工作應維持連續(xù)的選育遺傳改良與高遺傳多樣性水平之間的平衡。本研究中,翹嘴鱖選育群體的平均期望雜合度(He=0.7574)和平均多態(tài)信息含量(PIC=0.7197)均高于吳旭等[29](He= 0.5210,PIC=0.4071)和梅秋蘭等[30](He=0.4109,PIC=0.5043)報道的翹嘴鱖養(yǎng)殖群體的遺傳多樣性。這是由于本研究所構建的翹嘴鱖選育群體中包含由湖南和江蘇野生翹嘴鱖親本遠緣雜交形成的子一代群體,因而能保持較高的遺傳多樣性水平,有利于后續(xù)進一步的遺傳選育。

    3.3 SSR標記與性狀的相關性分析

    分子標記輔助選擇育種是通過與目的基因緊密連鎖或共分離的分子標記,直接篩選DNA目標區(qū)域,因而不受環(huán)境條件的影響,可以提高選擇的準確性和可靠性[4]。分子標記輔助選擇育種的關鍵是尋找對目標性狀具有較大遺傳貢獻率的主效基因或標記。分子遺傳標記與數(shù)量性狀間存在一定的關聯(lián)關系,即標記與控制這些數(shù)量性狀的基因處于連鎖不平衡狀態(tài)是確定影響目標性狀的主效基因或分子標記的先決條件。如王宣朋等[31]利用SSR、EST-SSR(Expressed Sequence Tag-SSR)、SNP標記對鯉的3種性狀進行相關性分析,發(fā)現(xiàn)分別有27、35和27個標記與食物轉化率、體厚和體質量顯著相關(P<0.05),且標記SNP0041、SNP0044與食物轉化率相關極顯著(P<0.001)。Cnaani等[32]分別利用20個和6個SSR標記對莫桑比克羅非魚Oreochromis mossambica與奧利亞羅非魚Oreochromis aureus的2個雜交F2群體進行分子標記與耐寒能力和魚體規(guī)格的QTL分析,結果發(fā)現(xiàn)在同一個連鎖群中的兩個相距22cM的微衛(wèi)星位點與目標性狀相關顯著,其中UNH879與羅非魚的耐寒能力顯著相關(P<0.05),UNH130與羅非魚的體質量顯著相關(P<0.05)。本研究發(fā)現(xiàn),微衛(wèi)星位點G4、G10、G12和G13位點與翹嘴鱖選育群體的體質量、體長和體高性狀均極顯著相關(P<0.01),說明這些標記位點與相關的性狀間存在關聯(lián)。在這些關聯(lián)關系中,發(fā)現(xiàn)了一個性狀的表現(xiàn)受多個不同標記位點的基因影響,或者一個標記位點可以影響多個性狀的表現(xiàn),說明這些位點中存在“多因一效”和“一因多效”的現(xiàn)象。

    3.4 基因型與性狀的連鎖分析

    基因型-性狀連鎖分析是對個體的數(shù)量性狀值與標記位點的基因型進行顯著性檢驗,差異顯著則表明基因型與數(shù)量性狀存在關聯(lián)[33]。如孫效文等[34]對28個多態(tài)性SSR標記的基因型與208尾鏡鋰Cyprinus carpio var.specularis繁殖群體的體質量值進行連鎖分析,發(fā)現(xiàn)有12個基因型與鏡鯉體質量相關,且有3個基因型緊密相關。本研究對數(shù)量性狀相關的微衛(wèi)星位點進行不同基因型性狀均值的多重比較,結果表明在G4位點中,含有219bp片段的等位基因的基因型(229/219或219/219)個體的體質量、體長和體高的表型效應顯著高于其他基因型;在G10位點中246/246基因型的體質量和體高表型效應顯著高于其他的基因型。上述結果表明,G4位點的229/219和219/219基因型可作為體質量、體長和體高選擇的輔助標記,而G10位點的246/246基因型可作為體質量和體高選擇的輔助標記。

    [1]方展強,陳軍,鄭文彪,等.鱖野生群體與養(yǎng)殖群體的RAPD分析[J].大連水產學院學報,2005,20(1):16-19.

    [2]He J G,Zeng K,Weng S P,et al.Systemic disease caused by an iridovirus-like agent in cultured mandarinfish Siniperca chuatsi(Basilewsky)in China[J].Journal ofFishDiseases,2000,23(3):219-222.

    [3] Fu X Z,Li N Q,Liu L H,et al.Genotype and host range analysis of infectious spleen and kidney necrosis virus(ISKNV)[J].Virus Genes,2011,42(1):97-109.

    [4]Rothschild MF and Ruvinsky A.Chapter 12 Marker-Assisted Selection for Aquaculture Species,in Aquaculture Genome Technologies[M].Oxford,UK:Blackwell PublishingLtd,2007:199-214.

    [5]張義鳳,張研,魯翠云,等.鯉魚微衛(wèi)星標記與體重、體長和體高性狀的相關分析[J].遺傳,2008,30(5):613-619.

    [6]Sun X W and Liang L Q.A genetic linkage map ofcommon carp(Cyprinus carpio L.)and mappingofa locus associated with cold tolerance[J].Aquaculture,2004,238(1): 165-172.

    [7]DrewR E,Schwabl R,Wheeler P A,et al.Detection ofQTL influencing cortisol levels in rainbow trout(Oncorhynchus mykiss)[J].Aquaculture,2007,272(1):183-194.

    [8]Sánchez-Molano E,Cerna A,Toro MA,et al.Detection of growth-related QTL in turbot(Scophthalmus maximus)[J].BMCGenomics,2011,12(1):473.

    [9]SauvageC,VagnerM,Der?meN,etal.CodinggeneSNP mappingrevealsQTLlinkedtogrowthandstressresponseinbrook charr(Salvelinus fontinalis)[J].G3:Genes|Genomes|Genetic,2012,2(6):707-720.

    [10]萬玉美,王蕾,譚照君,等.紅鰭東方鲀兩個群體的遺傳結構及與經濟性狀的相關性分析[J].淡水漁業(yè),2011, 41(5):9-16.

    [11] Wang H F,Sun J J,Lu X,et al.Identification of insulin-like growth factor I gene polymorphisms using high-resolution melting and its effect on growth traits in sinipercid species[J].Fisheries Science,2013,79(3): 439-446.

    [12]Miesfeld R,Krystal Mand AmheimN.A member ofa new repeated sequence family which is conserved throughout eucaryotic evolution is found between the human δ and β globin genes[J].Nucleic Acids Research,1981,9(22):5931-5948.

    [13]O’Connell Mand Wright J M.Microsatellite DNA in fishes[J].Reviews in Fish Biologyand Fisheries,1997,7(3): 331-363.

    [14]Loukovitis D,Sarropoulou E,Batargias C,et al.Quantitative trait loci for body growth and sex determination in the hermaphrodite teleost fish Sparus aurata L.[J].Animal Genetics,2012,43(6):753-759.

    [15]Liu F L,Yao J T,Wang X L,et al.Genetic diversity and structure within and between wild and cultivated Saccharina japonica(Laminariales,Phaeophyta)revealed by SSR markers[J].Aquaculture,2012,358-359:139-145.

    [16]鄭荷子,易提林,梁旭方,等.翹嘴鱖連續(xù)4代選育群體遺傳多樣性及遺傳結構分析[J].淡水漁業(yè),2013,43(6): 8-12.

    [17]楊凱,高銀愛,丹成,等.翹嘴鱖養(yǎng)殖與野生群體及其家系的遺傳多樣性分析[J].淡水漁業(yè),2013,43(3):3-8.

    [18]匡剛橋.鱖魚微衛(wèi)星標記的開發(fā)及遺傳多樣性研究[D].長沙:湖南農業(yè)大學,2007.

    [19]孫效文,張曉鋒,趙瑩瑩,等.水產生物微衛(wèi)星標記技術研究進展及其應用[J].中國水產科學,2008,15(4):689-703.

    [20]Zhang B,Li Z J,Tong J G,et al.Isolation and characterization of18 polymorphic microsatellite markers in Chinese mandarin fish Siniperca chuatsi(Basilewsky)[J].Molecular EcologyNotes,2006,6(4):1216-1218.

    [21]Kuang G,Lu S,Zheng S,et al.Isolation and evaluation of 18microsatellitemarkersinSiniperca chuatsi(Basilewsky)[J].MolecularEcologyResources,2009,9(6):1473-1475.

    [22] Liu X L,Luo W,Zeng C,et al.Isolation of new 40 microsatellite markers in mandarin fish(Siniperca chuatsi)[J].International Journal ofMolecular Sciences,2011,12(7):4180-4189.

    [23]Qu C M,LiangX F,HuangW,et al.Isolation and characterization of 46 novel polymorphic EST-simple sequence repeats (SSR) markersin two sinipercine fishes(Siniperca)and cross-species amplification[J].International Journal of Molecular Sciences,2012,13(8): 9534-9544.

    [24]Qu C M,Liang X F,Huang W,et al.Development and characterization of twenty-nine novel polymorphic microsatellite loci in the mandarin fish Sininperca chuatsi[J].Journal ofGenetics,2013,92(1):19-23.

    [25] Lindqvist A K,Magnusson P K,Balciuniene J,et al. Chromosome-specific panels of tri-and tetranucleotide microsatellite markers for multiplex fluorescent detection and automated genotyping:evaluation of their utility in pathologyandforensics[J].GenomeResearch,1996,6(12): 1170-1176.

    [26]Pongsomboon S,Whanb V,Mooreb S S,et al.Characterization of tri-and tetranucleotide microsatellite in the black tiger prawn,Penaeus monodon[J].Population Studies,2000,26:1-8.

    [27]Su G S,Liljedahl L E and Gall G A E.Effects of inbreeding on growth and reproductive traits in rainbow trout(Oncorhynchus mykiss)[J].Aquaculture,1996,142(3):139-148.

    [28]楊宇暉,梁旭方,林群,等.廣東與江西翹嘴鱖養(yǎng)殖與天然種群的遺傳多態(tài)性分析[J].水產學報,2010,34(4): 515-520.

    [29]吳旭,嚴美姣,李鐘杰.長江中下游不同地理種群鱖遺傳結構研究[J].水生生物學報,2010,34(4):843-850.

    [30]梅秋蘭,劉臻,張學文,等.鱖原種群體和養(yǎng)殖群體遺傳多樣性的微衛(wèi)星分析[J].基因組學與應用生物學,2010, 29(2):266-272.

    [31]王宣朋,張曉峰,李文升,等.利用SSR、EST-SSR、SNP標記對鯉食物轉化率、體厚、體質量3種性狀的分析[J].中國水產科學,2011,18(3):565-573.

    [32]Cnaani A,Hallerman E M,Ron M,et al.Detection of a chromosomal region with two quantitative trait loci affecting cold tolerance and fish size,in an F2tilapia hybrid[J].Aquaculture,2003,223(1):117-128.

    [33]王高富,吳登俊.涼山半細毛羊微衛(wèi)星標記與羊毛性狀的相關分析[J].遺傳,2006,28(12):1505-1512.

    [34]孫效文,魯翠云,匡友誼,等.鏡鯉兩個繁殖群體的遺傳結構和幾種性狀的基因型分析[J].水產學報,2007,31(3):273-279.

    Analysis of Microsatellite DNA Markers Emphasis on Correlation with Some Economically Important Traits in Mandarin fish Siniperca chuatsi

    SUN Ji-jia1,LI Gui-feng2,3,LIU Li1,WANG Hai-fang2,XU Peng2,3
    (1.College of Marine Science,South China Agricultural University,Guangzhou 510640,China; 2.South China Sea Bio-Resource Exploitation and Utilization Collaborative Innovation Center,Institute of Aquatic Economic Animals,School of Life Sciences,Sun Yat-sen University,Guangzhou 510275,China;3.Guangdong Engineering Technology Research Center,Foshan Nanhai Bairong Aquatic Breeding Co.Ltd.,Foshan 528216,China)

    The genomic DNA properties and the correlation between the microsatellite markers and three economic traits including body weight,body length,and body height were analyzed in 106 individuals of mandarin fish Siniperca chuatsi based on the seven high polymorphism microsatellite markers.A total of 67 different alleles were observed,with alleles ranging from 3 to 19 per locus, and the DNA fragment length was varied from 147 bp to 530 bp.The expected heterozygosity was found to be from 0.5092 to 0.9207 with mean value of 0.7574,and the polymorphism information content(PIC)to be ranged from 0.4639 to 0.9101 with mean value of 0.7197,indicating that the selected microsatellite markers were practicable for the genetic analysis and molecular marker assisted breeding of the mandarin fish.The correlation analysis revealed that the genotype in locus G4 with DNA fragment length of 219(229/219 or 291/219)had significantly higher phenotype effect on body weight,body length and body height than others.However, the genotype(246/246)in locus G10 had significantly higher phenotype effect on body weight and body height.Both loci appear to offer some important reference information in the future research on molecular marker assisted breeding of the mandarin fish.

    Siniperca chuatsi;microsatellite marker;economic trait;correlation analysis

    S917

    A

    2016-07-29

    廣東省戰(zhàn)略性新興產業(yè)核心技術攻關項目(2012A020800001);廣東省教育廳項目(cxzd1104);佛山市引進科技創(chuàng)新團隊資助計劃(2014IT100122).

    孫際佳(1976-),女,博士,從事水產動物遺傳育種研究.E-mail:jjsun@scau.edu.cn

    李桂峰(1963-),男,教授,博士生導師.E-mail:liguif@mail.sysu.edu.cn

    1005-3832(2017)01-0011-08

    猜你喜歡
    微衛(wèi)星體長基因型
    梅花鹿基因組微衛(wèi)星分布特征研究
    特產研究(2022年6期)2023-01-17 05:05:06
    脖子占體長一半的巨龍——新疆巨龍
    天天愛科學(2022年9期)2022-09-15 01:12:28
    史上最長恐龍
    外泌體長鏈非編碼RNA在膀胱癌中的研究進展
    林麝全基因組微衛(wèi)星分布規(guī)律研究
    四川動物(2017年4期)2017-07-31 23:54:19
    基于轉錄組測序的波紋巴非蛤微衛(wèi)星標記研究
    西安地區(qū)育齡婦女MTHFRC677T基因型分布研究
    中國“一箭雙星”成功將“遙感衛(wèi)星二十一號”與“天拓二號視頻微衛(wèi)星”發(fā)射升空
    河北遙感(2014年3期)2014-07-10 13:16:48
    BAMBI基因敲除小鼠的繁育、基因型鑒定
    海州灣及鄰近海域大瀧六線魚的生長、死亡和最適開捕體長研究
    丝袜脚勾引网站| 亚洲人成网站在线观看播放| 久久精品aⅴ一区二区三区四区| 宅男免费午夜| 欧美日韩一区二区视频在线观看视频在线| 中文字幕人妻熟女乱码| 亚洲,一卡二卡三卡| 老司机在亚洲福利影院| 在线观看一区二区三区激情| 色网站视频免费| 久久精品亚洲av国产电影网| 丰满迷人的少妇在线观看| 国产精品熟女久久久久浪| 久久精品国产a三级三级三级| 在线精品无人区一区二区三| 如日韩欧美国产精品一区二区三区| 国产精品三级大全| 水蜜桃什么品种好| 99久久综合免费| 国产日韩欧美视频二区| a级片在线免费高清观看视频| 国产成人精品无人区| 操美女的视频在线观看| 久久精品国产综合久久久| 欧美日韩亚洲高清精品| 大香蕉久久网| 日韩伦理黄色片| 黄片播放在线免费| 一边摸一边做爽爽视频免费| 黄片小视频在线播放| 一级片免费观看大全| 国产精品人妻久久久影院| 精品卡一卡二卡四卡免费| 国产精品免费视频内射| 免费黄频网站在线观看国产| 国产免费福利视频在线观看| 国产精品偷伦视频观看了| 一区二区av电影网| 日韩欧美精品免费久久| 精品国产一区二区三区四区第35| 久久久久久久久久久久大奶| 成人影院久久| 亚洲欧美激情在线| 亚洲欧美成人精品一区二区| 老鸭窝网址在线观看| 自拍欧美九色日韩亚洲蝌蚪91| 一区二区三区乱码不卡18| 丝瓜视频免费看黄片| 岛国毛片在线播放| 日本午夜av视频| 日韩成人av中文字幕在线观看| 国产日韩欧美视频二区| 亚洲精品美女久久av网站| 国产精品成人在线| 国产99久久九九免费精品| 久久久久久免费高清国产稀缺| kizo精华| 久久性视频一级片| 久久影院123| 久热爱精品视频在线9| 日韩精品免费视频一区二区三区| 日韩av不卡免费在线播放| 美女扒开内裤让男人捅视频| 亚洲国产精品成人久久小说| 国产又色又爽无遮挡免| 黄片无遮挡物在线观看| 亚洲精品av麻豆狂野| 国产成人欧美| 国产欧美亚洲国产| 极品人妻少妇av视频| 欧美日韩国产mv在线观看视频| 亚洲av成人不卡在线观看播放网 | 国产精品麻豆人妻色哟哟久久| 国产男女超爽视频在线观看| 亚洲一级一片aⅴ在线观看| 我要看黄色一级片免费的| 一级爰片在线观看| 交换朋友夫妻互换小说| 亚洲av电影在线观看一区二区三区| 97精品久久久久久久久久精品| 久久天堂一区二区三区四区| 日本色播在线视频| 一区二区三区四区激情视频| 亚洲精品国产av蜜桃| 国产免费现黄频在线看| 十分钟在线观看高清视频www| 人妻人人澡人人爽人人| 精品视频人人做人人爽| 免费人妻精品一区二区三区视频| 中国国产av一级| 亚洲欧美一区二区三区国产| 久久天堂一区二区三区四区| 别揉我奶头~嗯~啊~动态视频 | 成人国产av品久久久| 夫妻性生交免费视频一级片| 成年女人毛片免费观看观看9 | 亚洲三区欧美一区| 亚洲成av片中文字幕在线观看| 制服诱惑二区| 精品国产国语对白av| 女人高潮潮喷娇喘18禁视频| 亚洲欧美一区二区三区黑人| 亚洲精品国产色婷婷电影| 麻豆乱淫一区二区| 国产一区二区激情短视频 | 欧美精品一区二区免费开放| 久久 成人 亚洲| 九九爱精品视频在线观看| 免费高清在线观看日韩| 啦啦啦在线观看免费高清www| 亚洲精品日本国产第一区| 青春草视频在线免费观看| 亚洲av日韩在线播放| av网站在线播放免费| 中文乱码字字幕精品一区二区三区| 在线观看免费视频网站a站| 久久久亚洲精品成人影院| www.熟女人妻精品国产| 亚洲国产成人一精品久久久| 欧美日韩一区二区视频在线观看视频在线| 中文字幕最新亚洲高清| 亚洲图色成人| 精品免费久久久久久久清纯 | 美女国产高潮福利片在线看| 亚洲精品美女久久久久99蜜臀 | 国产成人免费无遮挡视频| 亚洲美女搞黄在线观看| 天天操日日干夜夜撸| 国产精品一区二区在线不卡| 韩国精品一区二区三区| 永久免费av网站大全| 亚洲欧洲日产国产| 纯流量卡能插随身wifi吗| 欧美最新免费一区二区三区| 国产激情久久老熟女| 亚洲精品成人av观看孕妇| 人人妻人人爽人人添夜夜欢视频| 亚洲国产av新网站| 国产老妇伦熟女老妇高清| 80岁老熟妇乱子伦牲交| 别揉我奶头~嗯~啊~动态视频 | 在线免费观看不下载黄p国产| 午夜福利乱码中文字幕| 婷婷色综合大香蕉| 老司机影院成人| 国产精品久久久人人做人人爽| 街头女战士在线观看网站| 亚洲 欧美一区二区三区| 亚洲国产精品国产精品| 一级,二级,三级黄色视频| 人人澡人人妻人| 男人舔女人的私密视频| 一级,二级,三级黄色视频| 久久99热这里只频精品6学生| 亚洲av成人不卡在线观看播放网 | 午夜福利免费观看在线| 男女边吃奶边做爰视频| 男人舔女人的私密视频| 黄频高清免费视频| 色视频在线一区二区三区| 日韩制服丝袜自拍偷拍| 久久久久久久久久久久大奶| 午夜日韩欧美国产| 国产午夜精品一二区理论片| 性少妇av在线| 日韩中文字幕欧美一区二区 | 婷婷色综合大香蕉| 最近2019中文字幕mv第一页| 午夜福利免费观看在线| 久久影院123| 老熟女久久久| 国产精品国产三级专区第一集| 亚洲色图综合在线观看| 国产免费又黄又爽又色| 日本一区二区免费在线视频| 97人妻天天添夜夜摸| 国产不卡av网站在线观看| 日本一区二区免费在线视频| 国产亚洲av高清不卡| 国产在线免费精品| 精品一区在线观看国产| 国产乱人偷精品视频| 纵有疾风起免费观看全集完整版| 国产成人免费无遮挡视频| 欧美97在线视频| 国产黄频视频在线观看| 亚洲四区av| 深夜精品福利| 一本色道久久久久久精品综合| 久热爱精品视频在线9| 大香蕉久久网| 欧美精品av麻豆av| 免费黄频网站在线观看国产| 99精品久久久久人妻精品| 高清欧美精品videossex| 亚洲成人国产一区在线观看 | 女人爽到高潮嗷嗷叫在线视频| 午夜免费观看性视频| av福利片在线| 亚洲欧美精品综合一区二区三区| 夫妻性生交免费视频一级片| 日韩成人av中文字幕在线观看| av在线老鸭窝| 在线观看免费午夜福利视频| 少妇人妻精品综合一区二区| 亚洲成人av在线免费| 夫妻午夜视频| 亚洲国产毛片av蜜桃av| 中文字幕人妻熟女乱码| 欧美日韩av久久| 欧美日本中文国产一区发布| 国产精品欧美亚洲77777| 中文字幕亚洲精品专区| 国产免费福利视频在线观看| 久久97久久精品| 搡老岳熟女国产| 亚洲色图 男人天堂 中文字幕| 欧美人与性动交α欧美软件| 亚洲国产日韩一区二区| 国产视频首页在线观看| 日韩视频在线欧美| 人人妻人人澡人人看| av女优亚洲男人天堂| 啦啦啦 在线观看视频| 国产男人的电影天堂91| av在线app专区| 日韩大码丰满熟妇| 国产老妇伦熟女老妇高清| 19禁男女啪啪无遮挡网站| 国产精品三级大全| 一级爰片在线观看| 日韩伦理黄色片| 狂野欧美激情性bbbbbb| 肉色欧美久久久久久久蜜桃| 我的亚洲天堂| 国产一卡二卡三卡精品 | 亚洲国产成人一精品久久久| 在线观看免费日韩欧美大片| 亚洲欧美中文字幕日韩二区| 国产av码专区亚洲av| 久久热在线av| 精品一品国产午夜福利视频| 久久这里只有精品19| 久久久久网色| 国产高清不卡午夜福利| 久久 成人 亚洲| 赤兔流量卡办理| 久久久精品国产亚洲av高清涩受| 欧美 日韩 精品 国产| 国产精品 国内视频| 嫩草影院入口| 免费在线观看完整版高清| 久久久久国产一级毛片高清牌| 麻豆精品久久久久久蜜桃| 国产高清不卡午夜福利| 国产成人免费无遮挡视频| 亚洲人成网站在线观看播放| 亚洲精品国产一区二区精华液| 久久久久久久国产电影| 两性夫妻黄色片| 亚洲精品国产av成人精品| 午夜免费观看性视频| av片东京热男人的天堂| 国产深夜福利视频在线观看| 少妇被粗大猛烈的视频| 啦啦啦啦在线视频资源| 亚洲精品视频女| 一级毛片我不卡| 亚洲成av片中文字幕在线观看| 免费久久久久久久精品成人欧美视频| 国产精品久久久av美女十八| 精品久久久久久电影网| 亚洲av综合色区一区| 欧美在线一区亚洲| 免费在线观看视频国产中文字幕亚洲 | 久久久久人妻精品一区果冻| avwww免费| 日韩一卡2卡3卡4卡2021年| 亚洲精品中文字幕在线视频| 日韩一本色道免费dvd| 日韩欧美精品免费久久| 老司机影院毛片| 欧美日韩视频高清一区二区三区二| 国产精品香港三级国产av潘金莲 | 99九九在线精品视频| 亚洲成av片中文字幕在线观看| 亚洲欧洲日产国产| 久久精品熟女亚洲av麻豆精品| 一区二区av电影网| 免费黄网站久久成人精品| 免费高清在线观看视频在线观看| 国产片特级美女逼逼视频| 欧美日韩综合久久久久久| 国产成人av激情在线播放| 午夜免费观看性视频| 一本色道久久久久久精品综合| 18禁观看日本| 欧美日韩一级在线毛片| 晚上一个人看的免费电影| 在线观看一区二区三区激情| 精品久久久久久电影网| av线在线观看网站| 国产极品粉嫩免费观看在线| 国产精品久久久久久久久免| 日本午夜av视频| 亚洲欧美激情在线| 日韩大码丰满熟妇| 国产日韩欧美视频二区| 日韩大码丰满熟妇| 麻豆乱淫一区二区| 一区二区三区精品91| 一边摸一边抽搐一进一出视频| 国产高清不卡午夜福利| 中文字幕高清在线视频| 色94色欧美一区二区| 久久99热这里只频精品6学生| 亚洲欧美成人精品一区二区| 制服人妻中文乱码| 国产又色又爽无遮挡免| 搡老乐熟女国产| 久热这里只有精品99| 国产男人的电影天堂91| 精品国产一区二区三区四区第35| 午夜激情久久久久久久| 国产免费又黄又爽又色| 久久久国产一区二区| 嫩草影视91久久| 观看av在线不卡| 天堂俺去俺来也www色官网| 久久免费观看电影| 亚洲精品第二区| 老鸭窝网址在线观看| 老司机影院成人| 免费女性裸体啪啪无遮挡网站| 嫩草影院入口| 欧美日韩综合久久久久久| 国产成人a∨麻豆精品| 久久久久精品人妻al黑| 丁香六月欧美| 伦理电影大哥的女人| 天美传媒精品一区二区| 91精品三级在线观看| 久久久久人妻精品一区果冻| 欧美黄色片欧美黄色片| 国产精品女同一区二区软件| 中文字幕人妻丝袜一区二区 | 亚洲人成77777在线视频| 亚洲,一卡二卡三卡| 欧美日本中文国产一区发布| 一二三四中文在线观看免费高清| 精品人妻熟女毛片av久久网站| 国产精品一区二区在线不卡| 人妻一区二区av| 亚洲欧洲日产国产| 国产又色又爽无遮挡免| 9色porny在线观看| 男女床上黄色一级片免费看| 啦啦啦中文免费视频观看日本| 国产毛片在线视频| 最近中文字幕高清免费大全6| 99re6热这里在线精品视频| 青青草视频在线视频观看| 国产伦理片在线播放av一区| 久久人人爽av亚洲精品天堂| 男女床上黄色一级片免费看| 一区二区三区四区激情视频| 亚洲成色77777| 一二三四中文在线观看免费高清| 亚洲人成77777在线视频| 大片免费播放器 马上看| 久久精品人人爽人人爽视色| 丰满少妇做爰视频| 中文字幕人妻丝袜制服| 久久久久网色| 亚洲精品自拍成人| 亚洲精品在线美女| 看非洲黑人一级黄片| 天天躁夜夜躁狠狠躁躁| 操美女的视频在线观看| 搡老岳熟女国产| 中文精品一卡2卡3卡4更新| 国产色婷婷99| 亚洲欧洲日产国产| 黄片小视频在线播放| 欧美日韩视频高清一区二区三区二| 高清视频免费观看一区二区| 日本av手机在线免费观看| 亚洲精品美女久久久久99蜜臀 | 纯流量卡能插随身wifi吗| 操出白浆在线播放| 乱人伦中国视频| 美女主播在线视频| 亚洲av成人不卡在线观看播放网 | 亚洲天堂av无毛| 国产精品蜜桃在线观看| 菩萨蛮人人尽说江南好唐韦庄| 伊人久久大香线蕉亚洲五| 少妇人妻 视频| 伊人久久大香线蕉亚洲五| 极品人妻少妇av视频| 精品亚洲成国产av| 国产男人的电影天堂91| 午夜免费男女啪啪视频观看| 国产精品女同一区二区软件| 亚洲中文av在线| 男男h啪啪无遮挡| 黑丝袜美女国产一区| 咕卡用的链子| 在线看a的网站| 下体分泌物呈黄色| 国产精品偷伦视频观看了| 丝袜在线中文字幕| 精品一区二区免费观看| 美女午夜性视频免费| 熟女少妇亚洲综合色aaa.| 久久精品aⅴ一区二区三区四区| av一本久久久久| 亚洲精品一二三| 9色porny在线观看| 久久久久久久久久久免费av| 一区二区三区乱码不卡18| 欧美xxⅹ黑人| 免费黄色在线免费观看| 免费观看a级毛片全部| www.自偷自拍.com| 一本—道久久a久久精品蜜桃钙片| 亚洲欧美色中文字幕在线| www.精华液| a级片在线免费高清观看视频| 日韩成人av中文字幕在线观看| 伦理电影大哥的女人| 一级毛片我不卡| 国产福利在线免费观看视频| 亚洲熟女毛片儿| 黄片无遮挡物在线观看| 亚洲av成人不卡在线观看播放网 | 又大又黄又爽视频免费| 国产亚洲av片在线观看秒播厂| 最近最新中文字幕免费大全7| 欧美久久黑人一区二区| 亚洲第一青青草原| 91精品伊人久久大香线蕉| 女人久久www免费人成看片| 女人爽到高潮嗷嗷叫在线视频| 男女免费视频国产| 青春草国产在线视频| 欧美日韩视频精品一区| av在线app专区| 高清av免费在线| 亚洲精品一二三| 国产av国产精品国产| 国产成人精品在线电影| a级毛片黄视频| 免费高清在线观看视频在线观看| 熟妇人妻不卡中文字幕| 爱豆传媒免费全集在线观看| 亚洲视频免费观看视频| 啦啦啦视频在线资源免费观看| 精品国产露脸久久av麻豆| 亚洲欧美中文字幕日韩二区| 午夜av观看不卡| 日韩制服丝袜自拍偷拍| 熟女少妇亚洲综合色aaa.| 亚洲,欧美精品.| 亚洲,一卡二卡三卡| 国产精品av久久久久免费| 久久久久网色| 午夜av观看不卡| 伊人久久国产一区二区| 国产亚洲午夜精品一区二区久久| 精品人妻熟女毛片av久久网站| 精品国产国语对白av| 天天操日日干夜夜撸| 啦啦啦啦在线视频资源| 免费不卡黄色视频| av视频免费观看在线观看| 午夜免费鲁丝| 99精国产麻豆久久婷婷| 97精品久久久久久久久久精品| 99香蕉大伊视频| 午夜福利视频在线观看免费| 国产人伦9x9x在线观看| 成人亚洲精品一区在线观看| 亚洲国产精品成人久久小说| 亚洲熟女精品中文字幕| 国产一区有黄有色的免费视频| 多毛熟女@视频| 在线天堂中文资源库| 我要看黄色一级片免费的| 最新的欧美精品一区二区| 晚上一个人看的免费电影| 成人三级做爰电影| 人妻一区二区av| 亚洲欧美一区二区三区久久| 日韩av在线免费看完整版不卡| 国产精品久久久久久久久免| 可以免费在线观看a视频的电影网站 | 精品国产一区二区久久| 亚洲国产毛片av蜜桃av| 久久人人爽av亚洲精品天堂| 香蕉丝袜av| 欧美成人午夜精品| 少妇人妻久久综合中文| 午夜免费男女啪啪视频观看| 亚洲成人手机| 国产女主播在线喷水免费视频网站| 国产无遮挡羞羞视频在线观看| 日韩制服骚丝袜av| av视频免费观看在线观看| 国精品久久久久久国模美| 在线看a的网站| 超碰97精品在线观看| 成人国产av品久久久| 亚洲国产精品一区三区| 99精国产麻豆久久婷婷| 一边摸一边做爽爽视频免费| 日本vs欧美在线观看视频| 秋霞伦理黄片| 久久97久久精品| 国产成人精品久久二区二区91 | 婷婷色麻豆天堂久久| 国产欧美日韩一区二区三区在线| 一级片'在线观看视频| 高清视频免费观看一区二区| 中文字幕精品免费在线观看视频| 1024视频免费在线观看| 亚洲av在线观看美女高潮| 国产男女内射视频| 成年人免费黄色播放视频| www日本在线高清视频| 中文欧美无线码| 精品人妻一区二区三区麻豆| 黑人猛操日本美女一级片| videosex国产| 777久久人妻少妇嫩草av网站| av免费观看日本| 日本欧美视频一区| 精品少妇久久久久久888优播| 51午夜福利影视在线观看| 久久久精品区二区三区| 色吧在线观看| 免费看不卡的av| 99久久综合免费| 亚洲av国产av综合av卡| 伊人久久国产一区二区| 成人漫画全彩无遮挡| 欧美日韩综合久久久久久| 黄网站色视频无遮挡免费观看| 一级毛片 在线播放| 午夜福利,免费看| 一区二区三区乱码不卡18| 久久精品亚洲av国产电影网| 国产成人系列免费观看| 91老司机精品| 最近最新中文字幕大全免费视频 | 性高湖久久久久久久久免费观看| 日韩 亚洲 欧美在线| 日韩不卡一区二区三区视频在线| 国产精品久久久人人做人人爽| 国产免费一区二区三区四区乱码| 亚洲精品久久久久久婷婷小说| 精品亚洲成a人片在线观看| 久久久久久久久免费视频了| 久久精品国产亚洲av高清一级| 女的被弄到高潮叫床怎么办| 亚洲av福利一区| 精品酒店卫生间| 久久久久国产精品人妻一区二区| e午夜精品久久久久久久| 深夜精品福利| 国产老妇伦熟女老妇高清| 国产欧美日韩综合在线一区二区| 午夜免费观看性视频| 9191精品国产免费久久| 一区二区三区精品91| 亚洲第一av免费看| 国精品久久久久久国模美| 一级,二级,三级黄色视频| 涩涩av久久男人的天堂| 日本av手机在线免费观看| 交换朋友夫妻互换小说| 亚洲国产精品国产精品| 熟女少妇亚洲综合色aaa.| 欧美老熟妇乱子伦牲交| 中文字幕av电影在线播放| 丝袜美腿诱惑在线| 成年人午夜在线观看视频| 久久99精品国语久久久| 性色av一级| 日韩不卡一区二区三区视频在线| 亚洲精华国产精华液的使用体验| 亚洲精品国产av成人精品| 别揉我奶头~嗯~啊~动态视频 | 搡老岳熟女国产| 啦啦啦视频在线资源免费观看| 亚洲av国产av综合av卡| 在线观看一区二区三区激情| 久久国产亚洲av麻豆专区| 亚洲激情五月婷婷啪啪| tube8黄色片| 啦啦啦视频在线资源免费观看| 这个男人来自地球电影免费观看 | 亚洲情色 制服丝袜| 永久免费av网站大全| 日日啪夜夜爽| 日本爱情动作片www.在线观看| 久久久亚洲精品成人影院| 天天躁狠狠躁夜夜躁狠狠躁| 黄色毛片三级朝国网站| 超色免费av| 精品视频人人做人人爽| 日韩欧美精品免费久久| 观看av在线不卡| 蜜桃在线观看..| 色婷婷久久久亚洲欧美|