• <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基因敲除小鼠的繁育、基因型鑒定
    海州灣及鄰近海域大瀧六線魚的生長、死亡和最適開捕體長研究
    成年av动漫网址| 国产精品香港三级国产av潘金莲 | 丰满迷人的少妇在线观看| xxx大片免费视频| 亚洲免费av在线视频| 成年人黄色毛片网站| 久久精品亚洲av国产电影网| 久久性视频一级片| 婷婷色综合大香蕉| 女人精品久久久久毛片| 色婷婷av一区二区三区视频| 国产精品亚洲av一区麻豆| 老司机午夜十八禁免费视频| 一本色道久久久久久精品综合| 女警被强在线播放| av片东京热男人的天堂| 欧美中文综合在线视频| 欧美大码av| 国产成人av激情在线播放| 少妇猛男粗大的猛烈进出视频| 一级毛片电影观看| 9色porny在线观看| 色播在线永久视频| 老司机影院成人| 女性被躁到高潮视频| 日韩中文字幕欧美一区二区 | 欧美精品一区二区大全| 亚洲黑人精品在线| 在线观看免费午夜福利视频| 精品久久蜜臀av无| 亚洲精品国产一区二区精华液| 亚洲精品国产色婷婷电影| 99久久99久久久精品蜜桃| 免费日韩欧美在线观看| 国产精品国产三级国产专区5o| 亚洲欧美清纯卡通| 亚洲成人免费电影在线观看 | 免费高清在线观看视频在线观看| 免费人妻精品一区二区三区视频| 亚洲熟女精品中文字幕| 性少妇av在线| av网站在线播放免费| 91老司机精品| 亚洲成人手机| 一级a爱视频在线免费观看| 啦啦啦啦在线视频资源| 日韩制服丝袜自拍偷拍| 午夜福利一区二区在线看| 视频区欧美日本亚洲| 久久久久久久久久久久大奶| 欧美日韩一级在线毛片| 多毛熟女@视频| 在线观看免费高清a一片| 曰老女人黄片| 国产爽快片一区二区三区| 国产在线观看jvid| 亚洲三区欧美一区| 午夜视频精品福利| 国产精品香港三级国产av潘金莲 | 久久久久久久久免费视频了| 亚洲欧美清纯卡通| 亚洲一卡2卡3卡4卡5卡精品中文| 成人亚洲精品一区在线观看| 欧美精品人与动牲交sv欧美| 波多野结衣一区麻豆| 国产精品免费大片| 精品高清国产在线一区| 一二三四在线观看免费中文在| 亚洲五月色婷婷综合| 好男人电影高清在线观看| 亚洲第一av免费看| 黄片小视频在线播放| 精品国产一区二区三区久久久樱花| 天天操日日干夜夜撸| 少妇人妻 视频| 中文字幕最新亚洲高清| 午夜福利视频精品| 水蜜桃什么品种好| 新久久久久国产一级毛片| 国产91精品成人一区二区三区 | 在线观看www视频免费| 午夜福利视频在线观看免费| 免费在线观看完整版高清| 日本a在线网址| 99香蕉大伊视频| 色综合欧美亚洲国产小说| 人人妻人人添人人爽欧美一区卜| 最近手机中文字幕大全| 丰满饥渴人妻一区二区三| 精品人妻1区二区| a级片在线免费高清观看视频| 看免费成人av毛片| 亚洲av国产av综合av卡| 欧美黄色淫秽网站| 国产免费一区二区三区四区乱码| 老司机在亚洲福利影院| 老熟女久久久| 欧美黄色片欧美黄色片| 免费观看人在逋| 日本色播在线视频| 多毛熟女@视频| 又黄又粗又硬又大视频| 免费在线观看影片大全网站 | 波野结衣二区三区在线| 久久精品亚洲av国产电影网| 各种免费的搞黄视频| 黄色怎么调成土黄色| 成年人午夜在线观看视频| 日韩制服丝袜自拍偷拍| 赤兔流量卡办理| 啦啦啦 在线观看视频| 天天操日日干夜夜撸| 午夜老司机福利片| 亚洲av电影在线观看一区二区三区| 国产精品一区二区免费欧美 | 老司机深夜福利视频在线观看 | 激情视频va一区二区三区| 久久精品国产a三级三级三级| 搡老岳熟女国产| 女性生殖器流出的白浆| 中国美女看黄片| 国产男人的电影天堂91| 日日夜夜操网爽| 亚洲欧美中文字幕日韩二区| 999精品在线视频| 男女午夜视频在线观看| 国产精品久久久人人做人人爽| 黑丝袜美女国产一区| 又大又爽又粗| 久久久久国产一级毛片高清牌| 亚洲欧美日韩高清在线视频 | 国产一区二区三区av在线| 一本—道久久a久久精品蜜桃钙片| 多毛熟女@视频| 丰满少妇做爰视频| 亚洲精品第二区| 丝袜喷水一区| 成人亚洲精品一区在线观看| 婷婷色麻豆天堂久久| 在线av久久热| 久久 成人 亚洲| 欧美成人午夜精品| 精品国产一区二区三区四区第35| 日韩,欧美,国产一区二区三区| 久久人妻福利社区极品人妻图片 | 久久久精品区二区三区| 两个人免费观看高清视频| 男女之事视频高清在线观看 | svipshipincom国产片| 一级毛片我不卡| 国产高清视频在线播放一区 | 亚洲国产毛片av蜜桃av| 国产免费现黄频在线看| 国产精品成人在线| 2018国产大陆天天弄谢| 91九色精品人成在线观看| 狠狠精品人妻久久久久久综合| 男女无遮挡免费网站观看| 国产伦人伦偷精品视频| 精品国产一区二区久久| 日本vs欧美在线观看视频| 黄片播放在线免费| 亚洲精品第二区| avwww免费| 99精国产麻豆久久婷婷| 亚洲人成电影免费在线| 两个人看的免费小视频| 亚洲国产精品成人久久小说| 精品国产超薄肉色丝袜足j| 两个人免费观看高清视频| 大话2 男鬼变身卡| 一本—道久久a久久精品蜜桃钙片| 久久精品久久精品一区二区三区| 国产片特级美女逼逼视频| 97在线人人人人妻| 老司机影院毛片| 亚洲精品自拍成人| 国产男女内射视频| 成年人午夜在线观看视频| 免费少妇av软件| 日韩熟女老妇一区二区性免费视频| 电影成人av| 国产亚洲欧美精品永久| 欧美 亚洲 国产 日韩一| 无遮挡黄片免费观看| 国产三级黄色录像| 亚洲人成电影免费在线| 亚洲精品av麻豆狂野| 天堂俺去俺来也www色官网| 精品一区二区三卡| 免费观看av网站的网址| 在线av久久热| 99久久99久久久精品蜜桃| 欧美亚洲日本最大视频资源| 国产在线观看jvid| 成年美女黄网站色视频大全免费| 男人操女人黄网站| 久久ye,这里只有精品| 丝袜脚勾引网站| 欧美在线黄色| 国产亚洲精品久久久久5区| 亚洲精品乱久久久久久| 女人爽到高潮嗷嗷叫在线视频| 国产成人欧美| 国产伦理片在线播放av一区| 久久精品亚洲av国产电影网| 精品一区二区三区av网在线观看 | 在线天堂中文资源库| 丝袜人妻中文字幕| www.熟女人妻精品国产| 国产精品久久久久久精品电影小说| 日本91视频免费播放| 日本黄色日本黄色录像| 国产精品一国产av| 伦理电影免费视频| 国产精品一区二区精品视频观看| 日韩制服丝袜自拍偷拍| 大香蕉久久网| 久9热在线精品视频| 激情五月婷婷亚洲| 手机成人av网站| avwww免费| 日韩av在线免费看完整版不卡| 日韩中文字幕视频在线看片| 一本大道久久a久久精品| 亚洲成人免费av在线播放| 日本a在线网址| 日本猛色少妇xxxxx猛交久久| 女人爽到高潮嗷嗷叫在线视频| 少妇粗大呻吟视频| 日韩熟女老妇一区二区性免费视频| 国产精品偷伦视频观看了| 国产亚洲午夜精品一区二区久久| 欧美国产精品va在线观看不卡| 中文字幕亚洲精品专区| 看免费成人av毛片| 午夜老司机福利片| 男男h啪啪无遮挡| 亚洲图色成人| 九草在线视频观看| 69精品国产乱码久久久| 一区二区三区四区激情视频| 黄色 视频免费看| 亚洲精品美女久久久久99蜜臀 | 一二三四在线观看免费中文在| 亚洲第一av免费看| 亚洲av日韩在线播放| 国产男人的电影天堂91| 欧美国产精品va在线观看不卡| 极品少妇高潮喷水抽搐| 精品少妇黑人巨大在线播放| 婷婷色综合大香蕉| 咕卡用的链子| 欧美在线黄色| 视频在线观看一区二区三区| 在线观看国产h片| 97精品久久久久久久久久精品| 精品一区二区三卡| 久久久久久久久久久久大奶| www.熟女人妻精品国产| 国产亚洲av高清不卡| 熟女av电影| 一区在线观看完整版| 免费一级毛片在线播放高清视频 | 精品福利永久在线观看| 欧美少妇被猛烈插入视频| 水蜜桃什么品种好| 久久亚洲国产成人精品v| 亚洲五月婷婷丁香| 久久热在线av| 亚洲精品中文字幕在线视频| 老司机影院毛片| 久久久久国产一级毛片高清牌| 中文字幕精品免费在线观看视频| 精品亚洲成a人片在线观看| 欧美日韩综合久久久久久| 人人妻人人澡人人看| 我要看黄色一级片免费的| 巨乳人妻的诱惑在线观看| 1024视频免费在线观看| 国产精品 欧美亚洲| 一级黄片播放器| 免费在线观看完整版高清| 中文字幕人妻熟女乱码| 一区二区三区乱码不卡18| 免费少妇av软件| 日韩av在线免费看完整版不卡| 欧美国产精品一级二级三级| 欧美日韩黄片免| 国产免费现黄频在线看| 精品久久蜜臀av无| 一本久久精品| 五月天丁香电影| 国产欧美日韩一区二区三 | 午夜免费男女啪啪视频观看| 久久久久久久久久久久大奶| 精品第一国产精品| 欧美成人精品欧美一级黄| 精品久久久久久久毛片微露脸 | 午夜福利在线免费观看网站| 一本—道久久a久久精品蜜桃钙片| 欧美乱码精品一区二区三区| 91成人精品电影| 纵有疾风起免费观看全集完整版| 一本色道久久久久久精品综合| 免费日韩欧美在线观看| 久热这里只有精品99| 免费人妻精品一区二区三区视频| 色视频在线一区二区三区| 精品人妻一区二区三区麻豆| 大话2 男鬼变身卡| 中文字幕最新亚洲高清| 久久久久久人人人人人| 丝袜美腿诱惑在线| 美女扒开内裤让男人捅视频| 久久久久久免费高清国产稀缺| 久久av网站| 又大又黄又爽视频免费| 久久国产亚洲av麻豆专区| 91成人精品电影| 十八禁人妻一区二区| 91精品伊人久久大香线蕉| 天天躁狠狠躁夜夜躁狠狠躁| 一级片免费观看大全| 国产99久久九九免费精品| 丝瓜视频免费看黄片| 国产女主播在线喷水免费视频网站| 亚洲国产中文字幕在线视频| 超碰97精品在线观看| 一二三四社区在线视频社区8| 女警被强在线播放| 在现免费观看毛片| 久久毛片免费看一区二区三区| 超色免费av| 99久久综合免费| h视频一区二区三区| 国产黄频视频在线观看| av国产精品久久久久影院| 国产亚洲av片在线观看秒播厂| 交换朋友夫妻互换小说| 亚洲,一卡二卡三卡| 在线观看人妻少妇| 亚洲精品久久久久久婷婷小说| 国产精品一区二区在线不卡| 久久久久久久大尺度免费视频| 久久国产精品大桥未久av| 国产激情久久老熟女| 欧美日韩黄片免| 久久人妻熟女aⅴ| 午夜免费鲁丝| tube8黄色片| 国产精品一国产av| 婷婷色综合www| 老汉色av国产亚洲站长工具| 精品欧美一区二区三区在线| 亚洲av综合色区一区| 亚洲七黄色美女视频| 一区二区三区四区激情视频| 久久狼人影院| 黄色片一级片一级黄色片| 国产黄色视频一区二区在线观看| 国产在线视频一区二区| 可以免费在线观看a视频的电影网站| 日韩 亚洲 欧美在线| 女性生殖器流出的白浆| 免费av中文字幕在线| 亚洲精品国产av蜜桃| 90打野战视频偷拍视频| 激情五月婷婷亚洲| 日本五十路高清| 久久精品成人免费网站| 色播在线永久视频| 少妇被粗大的猛进出69影院| 精品国产一区二区久久| 精品少妇一区二区三区视频日本电影| 精品人妻熟女毛片av久久网站| 国产色视频综合| 国产精品久久久久成人av| 18禁国产床啪视频网站| 亚洲av男天堂| 日韩一区二区三区影片| 国产高清不卡午夜福利| 丰满人妻熟妇乱又伦精品不卡| 悠悠久久av| 欧美精品一区二区大全| 免费在线观看日本一区| 在线观看免费视频网站a站| 国产精品国产三级专区第一集| 少妇猛男粗大的猛烈进出视频| 丝袜脚勾引网站| 岛国毛片在线播放| 精品久久蜜臀av无| 999久久久国产精品视频| 国产亚洲精品久久久久5区| 久久久久久久精品精品| 日韩 欧美 亚洲 中文字幕| 99精品久久久久人妻精品| 美女国产高潮福利片在线看| 捣出白浆h1v1| 欧美日韩亚洲国产一区二区在线观看 | 久久亚洲国产成人精品v| 七月丁香在线播放| 激情视频va一区二区三区| 国产免费福利视频在线观看| 国产1区2区3区精品| 久久久欧美国产精品| 色婷婷av一区二区三区视频| 国产精品一二三区在线看| 亚洲欧美一区二区三区黑人| 亚洲中文日韩欧美视频| 欧美变态另类bdsm刘玥| 亚洲伊人久久精品综合| 午夜两性在线视频| 国产女主播在线喷水免费视频网站| 精品一区二区三区av网在线观看 | 交换朋友夫妻互换小说| 好男人视频免费观看在线| 久久天堂一区二区三区四区| 在线av久久热| 91字幕亚洲| 久久99精品国语久久久| 欧美人与性动交α欧美软件| 欧美黄色片欧美黄色片| 亚洲综合色网址| 韩国精品一区二区三区| 成年美女黄网站色视频大全免费| 欧美大码av| 女人高潮潮喷娇喘18禁视频| 中文字幕人妻熟女乱码| 久久人人爽人人片av| 国产熟女欧美一区二区| av国产精品久久久久影院| 巨乳人妻的诱惑在线观看| 国产极品粉嫩免费观看在线| 精品国产超薄肉色丝袜足j| 国产又爽黄色视频| 婷婷色综合大香蕉| 亚洲五月色婷婷综合| 天堂中文最新版在线下载| 国产1区2区3区精品| 中文字幕人妻丝袜一区二区| 中文精品一卡2卡3卡4更新| 亚洲精品国产av蜜桃| 一级,二级,三级黄色视频| 天天躁日日躁夜夜躁夜夜| e午夜精品久久久久久久| 色精品久久人妻99蜜桃| av电影中文网址| 视频区欧美日本亚洲| 国产日韩一区二区三区精品不卡| www.自偷自拍.com| 国产色视频综合| 99热全是精品| 男的添女的下面高潮视频| 我要看黄色一级片免费的| 久久精品国产亚洲av涩爱| 久久影院123| 考比视频在线观看| 精品久久蜜臀av无| 日本五十路高清| 日韩伦理黄色片| 日韩中文字幕欧美一区二区 | 久久综合国产亚洲精品| 男女无遮挡免费网站观看| 人人妻人人澡人人看| 婷婷色av中文字幕| 国产精品国产三级专区第一集| 日本欧美国产在线视频| 色网站视频免费| 中文字幕av电影在线播放| 欧美人与善性xxx| av一本久久久久| 久热这里只有精品99| xxx大片免费视频| 国产野战对白在线观看| 亚洲国产av影院在线观看| 悠悠久久av| 两性夫妻黄色片| 精品少妇一区二区三区视频日本电影| 久久女婷五月综合色啪小说| 久久精品久久精品一区二区三区| 两个人免费观看高清视频| 欧美日韩亚洲国产一区二区在线观看 | 一本—道久久a久久精品蜜桃钙片| 超碰97精品在线观看| 十八禁人妻一区二区| 亚洲五月色婷婷综合| 纵有疾风起免费观看全集完整版| 久久精品国产综合久久久| 亚洲中文av在线| 国产99久久九九免费精品| 国产成人一区二区在线| 伦理电影免费视频| 赤兔流量卡办理| 精品一区在线观看国产| 精品一品国产午夜福利视频| 中文精品一卡2卡3卡4更新| 搡老岳熟女国产| 精品久久蜜臀av无| 91成人精品电影| 一级毛片 在线播放| 新久久久久国产一级毛片| 一区二区av电影网| 每晚都被弄得嗷嗷叫到高潮| 日本av免费视频播放| 一边摸一边抽搐一进一出视频| 波野结衣二区三区在线| 大片免费播放器 马上看| 亚洲中文av在线| 国产成人免费无遮挡视频| av有码第一页| 真人做人爱边吃奶动态| 亚洲成色77777| 精品少妇内射三级| 女性生殖器流出的白浆| 亚洲色图综合在线观看| 中文字幕人妻丝袜一区二区| 高潮久久久久久久久久久不卡| 涩涩av久久男人的天堂| 女人精品久久久久毛片| 欧美av亚洲av综合av国产av| 一本—道久久a久久精品蜜桃钙片| 久久精品国产a三级三级三级| 亚洲欧洲国产日韩| av线在线观看网站| 桃花免费在线播放| 国产成人系列免费观看| 午夜激情久久久久久久| 午夜免费观看性视频| 久久ye,这里只有精品| 亚洲,欧美精品.| 看免费av毛片| 亚洲中文av在线| 精品卡一卡二卡四卡免费| 亚洲人成77777在线视频| 亚洲国产最新在线播放| 香蕉丝袜av| www日本在线高清视频| 男的添女的下面高潮视频| 一级黄色大片毛片| 一区二区三区精品91| 亚洲精品中文字幕在线视频| 最黄视频免费看| 亚洲av在线观看美女高潮| 欧美黄色淫秽网站| 久久国产精品男人的天堂亚洲| 欧美黄色淫秽网站| 亚洲国产看品久久| 国产人伦9x9x在线观看| 欧美黄色片欧美黄色片| 午夜免费鲁丝| 日韩熟女老妇一区二区性免费视频| 啦啦啦中文免费视频观看日本| 免费看不卡的av| 中国国产av一级| 91精品三级在线观看| 免费一级毛片在线播放高清视频 | 熟女少妇亚洲综合色aaa.| 天天影视国产精品| 高清av免费在线| 免费看不卡的av| 高潮久久久久久久久久久不卡| 老汉色av国产亚洲站长工具| 国产成人免费无遮挡视频| 亚洲成色77777| 高清欧美精品videossex| 国产精品av久久久久免费| 亚洲一区中文字幕在线| 久久久亚洲精品成人影院| 一级黄色大片毛片| 国产色视频综合| 亚洲图色成人| 大码成人一级视频| h视频一区二区三区| 国产精品一区二区在线不卡| 亚洲成av片中文字幕在线观看| 一级毛片女人18水好多 | 啦啦啦在线免费观看视频4| 国产男女超爽视频在线观看| 美国免费a级毛片| 伊人亚洲综合成人网| 国产亚洲一区二区精品| 操出白浆在线播放| 91成人精品电影| 国产男女内射视频| 精品少妇黑人巨大在线播放| 大片免费播放器 马上看| 黄频高清免费视频| 侵犯人妻中文字幕一二三四区| 欧美精品一区二区大全| 精品国产超薄肉色丝袜足j| 日韩av不卡免费在线播放| 黄色视频在线播放观看不卡| 男的添女的下面高潮视频| 亚洲av综合色区一区| 极品人妻少妇av视频| 久久精品aⅴ一区二区三区四区| 国产高清国产精品国产三级| 成人免费观看视频高清| av不卡在线播放| 18在线观看网站| 欧美日韩亚洲高清精品| 最新在线观看一区二区三区 | 91精品国产国语对白视频| 国产在线观看jvid| 婷婷丁香在线五月| 国产一卡二卡三卡精品| 久久国产亚洲av麻豆专区| 欧美日韩综合久久久久久| 国产精品熟女久久久久浪| 国产黄频视频在线观看| 老熟女久久久| 多毛熟女@视频| 久久鲁丝午夜福利片| 无限看片的www在线观看|