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

    方斑東風(fēng)螺3個選育世代遺傳多樣性和遺傳結(jié)構(gòu)的微衛(wèi)星分析

    2022-11-26 05:59:32付敬強沈銘輝游偉偉柯才煥
    海洋科學(xué) 2022年10期
    關(guān)鍵詞:分析

    梁 園, 付敬強, 沈銘輝, 駱 軒, 游偉偉, 柯才煥

    方斑東風(fēng)螺3個選育世代遺傳多樣性和遺傳結(jié)構(gòu)的微衛(wèi)星分析

    梁 園1, 3, 付敬強2, 3, 沈銘輝4, 駱 軒1, 3, 游偉偉1, 3, 柯才煥1, 3

    (1. 廈門大學(xué) 海洋與地球?qū)W院, 福建 廈門 361102; 2. 廈門大學(xué) 環(huán)境與生態(tài)學(xué)院, 福建 廈門 361102; 3. 廈門大學(xué) 福建省海洋經(jīng)濟生物遺傳育種重點實驗室, 福建 廈門 361102; 4. 海南省海洋與漁業(yè)科學(xué)院 海南省熱帶海水養(yǎng)殖技術(shù)重點實驗室, 海南 ???571126)

    采用12個多態(tài)性微衛(wèi)星標(biāo)記對方斑東風(fēng)螺()泰國選育系TF4~TF6連續(xù)3代選育群體的遺傳多樣性與遺傳結(jié)構(gòu)進(jìn)行了分析。3個群體的多態(tài)信息含量分別為0.730、0.717和0.708, 均高于0.5, 表現(xiàn)出較高的多態(tài)性。平均等位基因數(shù)()、有效等位基因數(shù)(e)、觀測雜合度(o)、期望雜合度(e)和Shannon’s信息指數(shù)()的變化范圍分別為11.667~12.583、6.334~6.915、0.658~0.672、0.737~ 0.787和1.837~2.003, 整體呈降低趨勢, 但差異并不顯著(>0.05)。分子方差分析結(jié)果顯示, 群體間的遺傳變異僅占總變異的0.95%; 群體間的遺傳分化系數(shù)(st)介于0.00569~0.01324, 處于低等分化水平(st0.05)。主坐標(biāo)分析和STRUCTURE分析結(jié)果顯示, 3個群體呈現(xiàn)出相似的遺傳背景和遺傳結(jié)構(gòu)。以上研究結(jié)果表明, 經(jīng)過連續(xù)多代的人工選育, 方斑東風(fēng)螺選育群體保持著較高的遺傳多樣性水平, 遺傳變異和遺傳分化水平較低, 遺傳結(jié)構(gòu)趨于穩(wěn)定, 仍具有較高的遺傳選育潛力。

    方斑東風(fēng)螺(); 選育群體; 微衛(wèi)星標(biāo)記; 遺傳多樣性; 遺傳結(jié)構(gòu)

    方斑東風(fēng)螺()俗稱花螺, 是生活于熱帶、亞熱帶海域的腐肉食性淺海底棲腹足類動物[1], 在中國主要分布于海南、廣東、廣西、福建和臺灣等東南沿海地區(qū)[2, 3]。方斑東風(fēng)螺不僅生長速度快, 肉質(zhì)鮮美, 氨基酸含量高, 還富含人體所需的EPA和DHA[1], 深受消費者的喜愛。目前, 方斑東風(fēng)螺已成為中國重要的經(jīng)濟貝類養(yǎng)殖品種[4, 5]。隨著其養(yǎng)殖業(yè)的快速發(fā)展, 養(yǎng)殖過程中個體大小差異大, 生長速度減緩, 抗逆性差以及病害頻發(fā)等問題日益突出[4, 6], 因此, 培育出生長速度快、抗逆性強的方斑東風(fēng)螺優(yōu)良養(yǎng)殖品種具有十分重要的意義。

    選擇育種作為一種傳統(tǒng)的育種方法, 能夠有效改善育種目標(biāo)性狀[7]。在人工選育過程中, 部分不利基因被淘汰, 與目標(biāo)性狀相關(guān)的基因逐漸趨向于純合, 選育群體的遺傳結(jié)構(gòu)逐漸趨向于穩(wěn)定, 遺傳基礎(chǔ)會逐漸得到純化[8-10]。但與此同時, 選育可能導(dǎo)致群體遺傳多樣性水平出現(xiàn)不同程度的下降, 而較低的遺傳多樣性則不利于物種的生存和適應(yīng)[11-13]。為了對方斑東風(fēng)螺進(jìn)行遺傳改良, 從2011年開始, 本實驗室對其開展了生長性狀相關(guān)的群體選擇育種研究。此前, FU等[14]對該選育系TF0-TF33代選育群體進(jìn)行了遺傳多樣性分析, 結(jié)果表明, 經(jīng)過連續(xù)3代的人工選育, 選育群體的遺傳多樣性水平略有下降, 但差異不顯著。但是, 再次經(jīng)過連續(xù)3代的選育, 該選育系各世代群體的遺傳多樣性水平變化如何, 目前并不清楚。

    微衛(wèi)星標(biāo)記具有種類多、分布廣、多態(tài)信息含量高、操作簡單等優(yōu)點[15], 已被廣泛應(yīng)用于大黃魚()[8]、刺參()[16]、太平洋牡蠣()[17]、文蛤()[18]、海灣扇貝()[19]、馬氏珠母貝()[20]等水產(chǎn)動物的遺傳多樣性分析。為了解人工選育對方斑東風(fēng)螺泰國選育群體TF4~TF6分子遺傳水平上產(chǎn)生的影響, 本研究采用12個微衛(wèi)星標(biāo)記對這3代選育群體的遺傳多樣性和遺傳結(jié)構(gòu)進(jìn)行分析, 以期為進(jìn)一步的良種選育和遺傳改良提供理論指導(dǎo)。

    1 材料與方法

    1.1 實驗材料

    從2011年開始, 本實驗室針對生長性狀開展了連續(xù)6代的群體選擇育種研究。本研究所用到的選育群體(TF4、TF5、TF6)是2017—2020年連續(xù)3代從生長快的個體中截留親本選育而成, 選擇強度為10%, 具體方法見FU[14]的報道。從每代選育群體中隨機選取30個6月齡的樣品, 取適量腹足肌于無水乙醇中脫水保存?zhèn)溆谩?/p>

    1.2 實驗內(nèi)容與方法

    1.2.1 基因組DNA的提取

    使用AxyPrep基因組DNA小量試劑盒提取方斑東風(fēng)螺的基因組DNA。利用1%的瓊脂糖凝膠電泳檢測DNA的完整性, 使用NanoDrop 2000紫外分光光度計檢測DNA的濃度與純度, 于–20℃保存?zhèn)溆谩?/p>

    1.2.2 微衛(wèi)星引物的篩選

    從已報道的文獻(xiàn)[14, 21-23]中篩選出12對多態(tài)性較好的微衛(wèi)星引物用以遺傳多樣性分析。利用FAM、HEX和ROX 3種不同顏色的熒光在上游引物的5’端進(jìn)行標(biāo)記, 熒光引物由上海生工生物工程有限公司合成。引物的詳細(xì)信息見表1。

    1.2.3 PCR擴增與檢測

    PCR反應(yīng)體系體積為25 μL, 包括: 10×PCR Buffer (含Mg2+)5.0 μL、dNTPs (10 mmol/L) 2.0 μL、上游引物(10 mmol/L) 1.0 μL、下游引物(10 mmol/L) 1.0 μL、rTaq DNA Polymerase (5 U/μL) 0.25 μL、模板DNA 2.0 μL和ddH2O 13.75 μL。

    表1 12對微衛(wèi)星引物信息表

    F. forward primer, 上游引物; R. reverse primer, 下游引物; N. pure, 連續(xù)的微衛(wèi)星序列; n. interrupted, 有間隔的微衛(wèi)星序列

    PCR反應(yīng)程序為: 94 ℃ 3 min; 94 ℃ 30 s, 退火溫度30 s(引物退火溫度見表1), 72 ℃ 1 min, 30個循環(huán); 72 ℃ 10 min; 4 ℃保溫。

    PCR產(chǎn)物經(jīng)1%的瓊脂糖凝膠電泳檢測后, 將其置于凝膠成像系統(tǒng)中觀察是否擴增出目的片段。將PCR擴增產(chǎn)物送往上海生工生物工程有限公司, 使用ABI 3500xl基因分析儀進(jìn)行毛細(xì)管電泳, 檢測熒光信號; 使用軟件GeneMapper ID v3.2分析微衛(wèi)星位點的片段長度[14]。

    1.2.4 數(shù)據(jù)分析

    使用GenAlEx v6.5[24]計算等位基因數(shù)(number of alleles,)、有效等位基因數(shù)(effective number of alleles,e)、觀測雜合度(observed heterozygosity,o)、期望雜合度(expected heterozygosity,e)、Shannon’s信息指數(shù)(Shannon’s information index,)和固定系數(shù)(fixation index,)。使用PIC-CALC 0.6軟件計算多態(tài)信息含量(polymorphism information content,)。使用NeEsti-mator v2.0.1[25], 采用連鎖不平衡方法(linkage disequi-li-brium method)[26]和卡方檢驗估算有效群體大小(effect-tive population size,e), 最低等位基因頻率取0.01。

    使用Arlequin v3.5[27]進(jìn)行分子方差分析(analysis of molecular variance, AMOVA), 計算群體間遺傳分化指數(shù)(st)。使用GenAlEx 6.5軟件[24]進(jìn)行主坐標(biāo)分析(principal coordinates analysis, PCoA)。使用STRUCTURE v2.3.4[28]分析群體遺傳結(jié)構(gòu), 采用混合模型(admixture model), 設(shè)定群體間等位基因頻率相關(guān); 設(shè)置值為1到3, 每個值重復(fù)估算10次; 設(shè)定運行參數(shù)MCMC的不作數(shù)迭代(lenth of Burnin period)為100 000次, 迭代后的MCMC (number of MCMC reps after Burnin)為500 000次。使用CLUMPP 1.1.2[29]進(jìn)行重復(fù)抽樣分析, 并使用Distruct 1.1[30]將其結(jié)果轉(zhuǎn)化為圖形。

    2 結(jié)果

    2.1 遺傳多樣性

    12對微衛(wèi)星引物在連續(xù)3代選育群體(TF4、TF5和TF6)中共擴增獲得188個等位基因, 每個群體檢測到的總等位基因數(shù)分別為146、151和140, 每個微衛(wèi)星位點的等位基因數(shù)介于4~34。3個群體的遺傳多樣性參數(shù)如表2所示, 多態(tài)信息含量()分別為0.730、0.717和0.708, 均表現(xiàn)出較高的多態(tài)性(>0.5)[31]。另外, 等位基因數(shù)()、有效等位基因數(shù)(e)、觀測雜合度(o)、期望雜合度(e)和Shannon’s信息指數(shù)()的范圍分別為11.667±2.061~12.583±1.869、6.334±1.366~6.915±1.274、0.658±0.072~0.672±0.059、0.737±0.058~0.787±0.039和1.837±0.224~2.003±0.187。3個群體的平均觀測雜合度(o)均低于對應(yīng)的平均期望雜合度(e), 且平均固定系數(shù)()均大于0(表2), 表明存在雜合子缺失的情況??偟膩碚f, 經(jīng)過連續(xù)3代的選育, 各遺傳多樣性參數(shù)的大小整體呈降低趨勢, 但世代間的差異并不顯著(>0.05)。

    表2 方斑東風(fēng)螺3個群體的遺傳多樣性參數(shù)

    注: a.同一字母表示群體間差異不顯著(>0.05)

    此外, 獲得3個群體的有效群體大小(e)如表3所示, TF4~TF6依次為無窮大、223.0和628.9。

    2.2 遺傳分化與變異

    利用分子方差分析(AMOVA)和遺傳分化系數(shù)(st)對3個群體間的遺傳變異與遺傳分化進(jìn)行分析。AMOVA結(jié)果見表4, 絕大部分變異來源于各群體內(nèi)部(99.05%), 群體間的變異僅占總變異的0.95%。群體間的遺傳分化系數(shù)(st)介于0.005 69~ 0.013 24(表5), 其中TF4與TF5的遺傳分化水平最高(st=0.013 24), TF5與TF6的遺傳分化水平最低(st=0.005 69)且差異不顯著。根據(jù)White對遺傳分化系數(shù)的界定[32], 3個群體的遺傳分化均處于低等水平(st0.05)。

    2.3 遺傳結(jié)構(gòu)變化

    利用主坐標(biāo)分析(PCoA)比較3個群體的遺傳聚類(圖1)。所有個體重疊交錯分布, 群體間未表現(xiàn)出明顯的分離現(xiàn)象。STRUCTURE分析結(jié)果(圖2)顯示,=2或=3時, 3個群體均表現(xiàn)出相似的遺傳背景。遺傳分化和遺傳結(jié)構(gòu)的分析結(jié)果保持一致, 即3個群體間的遺傳差異較小, 遺傳結(jié)構(gòu)趨于穩(wěn)定。

    表3 方斑東風(fēng)螺3個群體的有效群體大小(Ne)

    表4 方斑東風(fēng)螺3個群體的分子方差分析(AMOVA)

    ST= 0.00953;<0.01

    表5 方斑東風(fēng)螺3個群體的遺傳分化系數(shù)(Fst)

    **.<0.01; ns.>0.05

    圖1 方斑東風(fēng)螺3個群體的主坐標(biāo)分析(PCoA)

    圖2 方斑東風(fēng)螺3個群體的STRUCTURE分析

    3 討論

    遺傳多樣性是生物多樣性的核心[33], 一般來說, 種群的遺傳多樣性越高, 其適應(yīng)環(huán)境的能力越強, 相應(yīng)的選育潛力也越大[34]。在針對生長性狀的人工選育過程中, 目的性狀相關(guān)基因會不斷趨于純合[35-37], 同時部分與適應(yīng)性和抗逆性相關(guān)的等位基因很容易丟失, 選育群體的遺傳多樣性可能會有所下降, 甚至發(fā)生近交衰退而失去繼續(xù)選育的價值[38]。選育群體的遺傳多樣性水平隨著選育世代的推進(jìn)而不斷下降的現(xiàn)象在太平洋牡蠣[39]、馬氏珠母貝[40]、中國對蝦()[41]、羅氏沼蝦()[42]、尼羅羅非魚()[43-44]等水產(chǎn)動物中均有報道。此前FU等[14]對方斑東風(fēng)螺海南和泰國兩個選育系(HT0~HT3, TH0~TH3)的遺傳多樣性水平進(jìn)行了比較分析, 結(jié)果顯示選育群體存在微弱但不顯著的下降趨勢。因此, 有必要繼續(xù)對方斑東風(fēng)螺人工選育過程中各世代群體遺傳多樣性的變化進(jìn)行監(jiān)測。

    有效等位基因數(shù)和期望雜合度等參數(shù)是衡量選育群體遺傳多樣性水平的重要指標(biāo)[45]。本研究中, 由遺傳多樣性參數(shù)所揭示的3個群體的遺傳多樣性水平接近, 相比于TF4和TF5, TF6的遺傳多樣性水平略有下降, 但群體間的差異并不顯著(>0.05), 3個群體均保持著較高的遺傳多樣性水平。此前, FU等[14]利用10對微衛(wèi)星標(biāo)記獲得了TF0~TF34個群體的遺傳多樣性參數(shù), 其中等位基因數(shù)、有效等位基因數(shù)和Shannon’s信息指數(shù)的變化范圍分別為8.10±0.89~ 10.30±1.32、4.26±0.48~5.22±0.90和1.62±0.13~1.78±0.17;而本研究中TF4~TF6對應(yīng)參數(shù)的變化范圍分別為11.667±2.06~12.583±1.869、6.334±1.366~6.915±1.274和1.837±0.224~2.003±0.187, 高于TF0~TF3。本研究使用了12對微衛(wèi)星標(biāo)記, 與TF0~TF3相比有所增加, 微衛(wèi)星標(biāo)記數(shù)量和多態(tài)性的不同有可能是TF4~TF6部分遺傳多樣性參數(shù)升高的原因[46-48]。此外, 在養(yǎng)殖或選育群體中, 相比于等位基因數(shù)和有效等位基因數(shù), 短時間內(nèi)外界環(huán)境因素對雜合度的影響更小[39, 49]。NEI[50]認(rèn)為, 在小群體中, 期望雜合度是度量群體遺傳變異的最適參數(shù)。本研究中, TF4~TF6的期望雜合度(0.737±0.058~0.787±0.039)未表現(xiàn)出明顯的下降趨勢, 并且與TF0~TF3的水平較一致(0.726±0.040~ 0.759±0.028)。由此可見, 再次經(jīng)過連續(xù)3代的人工選育, 方斑東風(fēng)螺泰國選育群體仍保持著較高的遺傳多樣性水平, 具有一定的遺傳選育潛力。

    有效群體大小往往是影響群體近交率和遺傳漂變率的重要參數(shù)[12]。在連續(xù)多代的人工選育過程中, 不合理的育種操作常常會引起選育群體有效群體數(shù)量下降而近交幾率上升[51]。為了解有效群體大小對遺傳多樣性變化的影響, LIND等[52]在珠母貝(.)的3個野生群體與5個人工培育群體中進(jìn)行了比較分析, 發(fā)現(xiàn)兩種群體的有效群體大小變化范圍分別為109.6~432.4和3.5~9.2, 而等位基因豐富度的變化范圍分別為12.6±1.4~14.0±1.2和7.5±1.1~9.4±1.4, 人工培育群體的等位基因豐富度與野生群體相對比減少了29%~44%。另外, 在南非鮑()的選擇育種中, 當(dāng)有效群體大小從F1的471.6減少到F2的52.5時, F2的等位基因大約減少了30%[53]。在FU等[14]的前期研究中, TF0~TF3的有效群體大小e從最初的436.0逐漸降至151.8, 且選育群體遺傳多樣性的降低并不顯著。而本研究中TF4~TF6的有效群體大小的變化為無窮大、223.0和628.9, 相對于TF3有所增加, 這可能是3個群體保持較高遺傳多樣性的原因之一。維持一定的有效群體數(shù)量能夠在一定程度上降低近交和遺傳多樣性損失所帶來的負(fù)面影響[49], 因此, 在后續(xù)的選育過程中仍需重視有效群體數(shù)量的管理。

    眾多研究表明, 在選擇壓力的驅(qū)動下, 選育群體相鄰世代間的遺傳分化會逐漸減弱, 遺傳相似性逐漸升高, 遺傳結(jié)構(gòu)與選育性狀會逐漸穩(wěn)定[8, 37, 38]。本研究中3個群體整體的遺傳變異(0.95%)和遺傳分化水平(ST= 0.00953)很低, TF5與TF6的遺傳分化水平(st: 0.00569,0.05)與TF0~TF3中相鄰世代(st: 0.0103~0.0308)相比明顯下降[14]。此外, PCoA分析和STRUCTURE分析的結(jié)果均顯示出3個群體在遺傳結(jié)構(gòu)上的相似性。由此可知, 經(jīng)過連續(xù)多代的人工選育, 方斑東風(fēng)螺選育群體的遺傳分化和遺傳變異水平逐漸降低, 遺傳結(jié)構(gòu)逐漸趨于穩(wěn)定, 經(jīng)進(jìn)一步選育可望獲得較穩(wěn)定的品系。

    4 結(jié)論

    本研究利用12對微衛(wèi)星標(biāo)記對方斑東風(fēng)螺3代選育群體(TF4~TF6)的遺傳多樣性與遺傳結(jié)構(gòu)進(jìn)行了分析比較。結(jié)果顯示, 方斑東風(fēng)螺選育群體的遺傳多樣性水平較高且未表現(xiàn)出明顯的下降趨勢, 遺傳變異與遺傳分化水平較低, 遺傳結(jié)構(gòu)趨于穩(wěn)定。連續(xù)多代的人工選育并未對其遺傳多樣性和遺傳結(jié)構(gòu)造成不利影響, 該選育群體仍具有較高的遺傳選育潛力。在后續(xù)的選育進(jìn)程中, 可繼續(xù)保持一定的選擇強度, 同時應(yīng)注意增加有效親本的數(shù)量, 減少近交發(fā)生的幾率, 從而確保方斑東風(fēng)螺優(yōu)良性狀的穩(wěn)定遺傳。

    [1] 許貽斌, 沈銘輝, 魏永杰, 等. 兩種東風(fēng)螺的營養(yǎng)成分分析與評估[J]. 臺灣海峽, 2008, 27(1): 26-32.

    XU Yibin, SHEN Minghui, WEI Yongjie, et al. Analysis and evaluation of nutritional composition ofLink and[J]. Journal of Ocreanography in Taiwan Strait, 2008, 27(1): 26-32.

    [2] 陳利雄, 吳進(jìn)鋒. 東風(fēng)螺的增養(yǎng)殖技術(shù)及產(chǎn)業(yè)化前景[J]. 齊魯漁業(yè), 2004, 21(10): 9-11.

    CHEN Lixiong, WU Jinfeng. Culture technique and industrial prospect of ivory shell[J]. Shandong Fisheries, 2004, 21(10): 9-11.

    [3] 張漢華, 吳進(jìn)鋒, 陳利雄, 等. 東風(fēng)螺人工育苗、養(yǎng)殖及產(chǎn)業(yè)化發(fā)展前景[J]. 南方水產(chǎn)科學(xué), 2004, 11: 2-5.

    ZHANG Hanhua, WU Jinfeng, CHEN Lixiong, et al. Development prospect of artificial seedling, breeding and industrialization of[J]. South China Fisheries Science, 2004, 11: 2-5.

    [4] 沈銘輝, 符芳霞, 呂文剛, 等. 海南省東風(fēng)螺養(yǎng)殖產(chǎn)業(yè)現(xiàn)狀和展望[J]. 安徽農(nóng)業(yè)科學(xué), 2015, 43(26): 144- 145.

    SHEN Minghui, FU Fangxia, LV Wengang, et al. The current situation ofspp. farming industry in Hainan and prospect outlook[J]. Journal of Anhui Agriculture Sciences, 2015, 43(26): 144-145.

    [5] 梁飛龍, 毛勇, 余祥勇, 等. 方斑東風(fēng)螺人工育苗試驗[J]. 海洋湖沼通報, 2005, 1(14): 79-85.

    LIANG Feilong, MAO Yong, YU Xiangyong, et al. Experiment on artificial breeding of(Lamarck)[J]. Transactions of Oceanology and Limno-logy, 2005, 1(14): 79-85.

    [6] 王菁, 劉付柏, 許尤厚, 等. 基于轉(zhuǎn)錄組測序的方斑東風(fēng)螺單核苷酸多態(tài)性位點挖掘及功能注釋[J]. 廣東海洋大學(xué)學(xué)報, 2021, 21(1): 111-118.

    WANG Jing, LIU Fubo, XU Youhou, et al. SNP site biological analysis ofbased on RNA-seq technology[J]. Journal of Guangdong Ocean University, 2021, 21(1): 111-118.

    [7] 張守都. 海灣扇貝的選擇和雜交育種[D]. 青島: 中國科學(xué)院海洋研究所, 2013.

    ZHANG Shoudu. Selective and crossing breeding of Bay scallop,(Lamarck)[D]. Qingdao: Institute of Oceanology, Chinese Academy of Sciences, 2013.

    [8] 趙廣泰, 劉賢德, 王志勇, 等. 大黃魚連續(xù)4代選育群體遺傳多樣性與遺傳結(jié)構(gòu)的微衛(wèi)星分析[J]. 水產(chǎn)學(xué)報, 2010, 34(4): 500-507.

    ZHAO Guangtai, LIU Xiande, WANG Zhiyong, et al. Genetic structure and genetic diversity analysis of four consecutive breeding generations of large yellow croa-ker () using microsatellite marke-rs[J]. Journal of Fisheries of China, 2010, 34(4): 500- 507.

    [9] 郭紅軍, 羅潔, 洪一江, 等. 人工選育池蝶蚌的生長及不同世代遺傳分析[J]. 水生生物學(xué)報, 2008, 32(2): 220-224.

    GUO Hongjun, LUO Jie, HONG Yijiang, et al. Growth and genetic analysis on the selected breeding[J]. Acta Hydrobiologica Sinica, 2008, 32(2): 220-224.

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

    ZHENG Hezi, YI Tilin, LIANG Xufang, et al. Genetic structure and genetic diversity analysis of four cones-cutive breeding generations of[J]. Freshwater Fisheries, 2013, 43(6): 8-12.

    [11] 于愛清, 施永海, 鄧平平. 長江刀鱭選育和野生群體遺傳多樣性的微衛(wèi)星分析[J]. 水產(chǎn)科技情報, 2019, 46(3): 121-125.

    YU Aiqing, SHI Yonghai, DENG Pingping. Microsa-tellite analysis of genetic diversity of wild population and selection ofin the Yangtze River[J]. Aquatic Science and Technology Information, 2019, 46(3): 121-125.

    [12] CHEN N, LUO X, LU C K, et al. Effects of artificial selection practices on loss of genetic diversity in the Pacific abalone,[J]. Aquaculture Research, 2017, 48(9): 4923-4933.

    [13] XU L, LI Q, XU C X, et al. Genetic diversity and effe-ctive population size in successive mass selected generations of black shell strain Pacific oyster () based on microsatellites and mtDNA data[J]. Aquaculture, 2019, 500: 338-346.

    [14] FU J Q, LV W G, LI W D, et al. Comparative assess-ment of the genetic variation in selectively bred generations from two geographic populations of ivory shell ()[J]. Aquaculture Research, 2017, 48(8): 1-14.

    [15] 黎裕, 賈繼增, 王天宇. 分子標(biāo)記的種類及其發(fā)展[J]. 生物技術(shù)通報, 1999, 4: 19-22.

    LI Yu, JIA Jizeng, WANG Tianyu. Types of molecular markers and their development[J], Biotechnology Information, 1999, 4: 19-22.

    [16] 孫孝德, 孫國華, 袁延柱, 等. 刺參()野生及兩代選育群體間遺傳變異的微衛(wèi)星標(biāo)記研究[J]. 海洋與湖沼, 2016, 85(1): 92-99.

    SUN Xiaode, SUN Guohua, YUAN Yanzhu, et al. Genetic variations of wild and selectively bred popula-tions of a sea cucumber[J]. Ocea-no-lo-gia et Limnologia Sinica, 2016, 85(1): 92-99.

    [17] 王慶志, 孔令鋒, 李琪.長牡蠣3代人工選育群體的微衛(wèi)星分析[J]. 水產(chǎn)學(xué)報, 2012, 36(10): 1529-1536.

    WANG Qingzhi, KONG Lingfeng, LI Qi. Genetic variability assessed by microsatellites in mass selection lines of Pacific oyster ()[J]. Journal of Fisheries of China, 2012, 36(10): 1529-1536.

    [18] 鄭培. 文蛤三個選育世代的多樣性研究[D]. 上海: 上海海洋大學(xué), 2013.

    ZHENG Pei. Genetic diversity analisis of three gene-rations of selective breeding of[D]. Shanghai: Shanghai Ocean University, 2013.

    [19] WANG L L, ZHANG H, SONG L S, et al. Loss of allele diversity in introduced populations of the herma-phroditic bay scallop[J]. Aqua-culture, 2007, 271(1/4): 252-259.

    [20] 王愛民, 王嫣, 顧志峰, 等. 馬氏珠母貝()2個地理群體雜交子代的雜種優(yōu)勢和遺傳變異[J]. 海洋與湖沼, 2010, 41(1): 140-146.

    WANG Aimin, WANG Yan, GU Zhifeng, et al. Heterosis and genetic variation of hybrids from two geographical populations of pearl oyster,[J]. Oceanologia et Limnologia Sinica, 2010, 41(1): 140-146.

    [21] CHEN F, KE C H, WANG D X, et al. Isolation and characterization of microsatellite loci inand cross-species amplification in[J]. Molecular Ecology Resources, 2009, 9(2): 661-663.

    [22] WANG Y, LU H, ZHENG J, et al. Eight polymorphic microsatellite markers for the spotted babylon,(Buccinidae)[J]. Genetics & Molecular Research, 2011, 10(4): 3230-3235.

    [23] ZHANG N, YING Q, HUANG X Z, et al. Micro-satellite marker development and characterization in the spotted babylon,(Link, 1807): detection of duplicated loci at high frequency[J]. International Journal of Aquaculture, 2012, 2(2): 5-10.

    [24] PEAKALL R, SMOUSE P E. GenAlEx 6.5: genetic analysis in Excel. Population genetic software for tea-ching and research-an update[J]. Bioinformatics, 2012, 28(19): 2537-2539.

    [25] DO C, WAPLES R S, PEEL D, et al. NeEstimator v2: re-implementation of software for the estimation of contemporary effective population size (e) from genetic data[J]. Molecular Ecology Resources, 2014, 14(1): 209-214.

    [26] WAPLES R S. A bias correction for estimates of effe-ctive population size based on linkage disequili-brium at unlinked gene loci[J]. Conservation Genetics, 2006, 7(2): 167-184.

    [27] EXCOFFIER L, LISCHER H. Arlequin suite ver 3.5: a new series of programs to perform population genetics analyses under Linux and Windows[J]. Molecular Eco-logy Resources, 2010, 10(3): 564-567.

    [28] EARL D A, VONHOLDT B M. Structure harvester: a website and program for visualizing structure output and implementing the evanno method[J]. Conservation Genetics Resources, 2012, 4(2): 359-361.

    [29] JAKOBSSON M, ROSENBERG N A. Clumpp: a clu-ster matching and permutation program for dealing with label switching and multimodality in analysis of population structure[J]. Bioinformatics, 2007, 23(14): 1801-1806.

    [30] ROSENBERG N A. Distruct: a program for the gra-phical display of population structure[J]. Molecular Ecology Notes, 2003, 4(1): 137-138.

    [31] BOTSTEIN D, WHITE R L, SKOLNICK M, et al. Construction of a genetic linkage map in man using restriction fragment length polymorphisms[J]. Ame-rican Journal of Human Genetics, 1980, 32(3): 314.

    [32] WRIGHT S. Variability within and among natural popu-la-tions[M]. Chicago: University of Chicago Press, 1978.

    [33] 錢迎倩. 生物多樣性研究的原理與方法[M]. 北京: 中國科學(xué)技術(shù)出版社, 1994.

    QIAN Yingqian. Principles and methods of biodiversity research[M]. Beijing: Science and Technology of China Press, 1994.

    [34] 馬靜, 安永平, 王彩芬, 等. 遺傳多樣性研究進(jìn)展[J]. 陜西農(nóng)業(yè)科學(xué), 2010, 56(1): 126-130.

    MA Jing, AN Yongping, WANG Caifen, et al. Advances in genetic diversity research[J]. Shaanxi Agricultural Sciences, 2010, 56(1): 126-130.

    [35] 張進(jìn), 梁旭方, 楊敏, 等. 2個鱖魚選育群體遺傳多樣性分析[J]. 水產(chǎn)科學(xué), 2014, 33(7): 447-450.

    ZHANG Jin, LIANG Xufang, YANG Min, et al. Genetic structure and genetic diversity analysis of artifitial selection populations of mandarin fish[J]. Fisheries Science, 2014, 33(7): 447-450.

    [36] 馬春艷, 馬洪雨, 馬凌波, 等. 凡納濱對蝦引進(jìn)群體和2個養(yǎng)殖群體遺傳變異的微衛(wèi)星分析[J]. 海洋漁業(yè), 2011, 31(1): 1-8.

    MA Chunyan, MA Hongyu, MA Lingboet al. Micro-satellite analysis on genetic variation of importedpopulation and cultured stocks[J]. Marine Fisheries, 2011, 31(1): 1-8.

    [37] YI T L, GUO W J, LIANG X F, et al. Microsatellite analysis of genetic diversity and genetic structure in five consecutive breeding generations of mandarin fish(Basilewsky)[J]. Genetics and Molecular Research, 2015, 14(1): 2600-2607.

    [38] 馬冬梅, 蘇換換, 朱華平, 等. 華南鯉選育群體不同世代遺傳多樣性與遺傳結(jié)構(gòu)的微衛(wèi)星分析[J]. 水生生物學(xué)報, 2018, 42(5): 887-895.

    MA Dongmei, SU Huanhuan, ZHU Huaping, et al. Genetic diversity and genetic structure analysis of different selective breeding generations inusing microsatellite markers[J]. Acta Hydrobio-logica Sinica, 2018, 42(5): 887-895.

    [39] HAN Z Q, LI Q, LIU S K, et al. Genetic variability of an orange-shell line of the Pacific oysterduring artificial selection inferred from microsa-tellites and mitochondrial COI sequences[J]. Aquaculture, 2019, 508(May): 159-166.

    [40] 管云雁, 劉文廣, 何毛賢. 馬氏珠母貝選育群體4個世代的遺傳變異[J]. 中國水產(chǎn)科學(xué), 2013, 20(4): 764-770.

    GUAN Yunyan, LIU Wenguang, HE Maoxian. Genetic variation during four generations of selective breeding in the pearl oyster[J]. Journal of Fishery Sciences of China, 2013, 20(4): 764-770.

    [41] LI Z X, LI J, WANG Q Y, et al. The effects of selective breeding on the genetic structure of shrimppopulations[J]. Aquaculture, 2006, 258(1/4): 278-282.

    [42] 董丁健, 戴習(xí)林. 羅氏沼蝦不同群體世代遺傳多樣性的SSR分析[J]. 南方農(nóng)業(yè)學(xué)報, 2020, 51(2): 421-428.

    DONG Dingjian, DAI Xilin. SSR analysis on genetic diversity in different populations and generations of[J]. Journal of Southern Agriculture, 2020, 51(2): 421-428.

    [43] 唐首杰, 李思發(fā), 趙金良.“新吉富”羅非魚選育后期世代遺傳潛力的微衛(wèi)星分析[J]. 海洋湖沼通報, 2018, 3: 120-131.

    TANG Shoujie, LI Sifa, ZHAO Jinliang. Heritability analysis of the late selection generations of nile tilapia () NEW GIFT strain using micro-satellites[J]. Transactions of Oceanology and Limoolo-gy, 2018, 3: 120-131.

    [44] 劉志剛, 曹建萌, 高風(fēng)英, 等. 羅非魚“粵閩1號”母本選育群體世代間遺傳差異的微衛(wèi)星分析[J]. 大連海洋大學(xué)學(xué)報, 2021, 36(1): 16-22.

    LIU Zhigang, CAO Jianmeng, GAO Fengying, et al. Genetic differentiation analysis of maternal selective breeding generations of tilapia “Yuemin No. 1”using microsatellites[J]. Journal of Dalian Ocean University, 2021, 36(1): 16-22.

    [45] 楊弘, 李大宇, 曹祥, 等. 微衛(wèi)星標(biāo)記分析羅非魚群體的遺傳潛力[J]. 遺傳, 2011, 33(7): 768-775.

    YANG Hong, LI Dayu, CAO Xiang, et al. Genetic potential analysis of six Tilapia populations by micro-satellite DNA markers[J]. Hereditas (Beijing), 2011, 33(7): 768-775.

    [46] 閆路娜, 張德興. 種群微衛(wèi)星 DNA 分析中樣本量對各種遺傳多樣性度量指標(biāo)的影響[J]. 動物學(xué)報, 2004, 50(2): 279-290.

    YAN Luna, ZHANG Dexing. Effects of sample size on various genetic diversity measures in population genetic study with microsatellite DNA markers[J]. Acta Zoologica Sinica, 2004, 50(2): 279-290.

    [47] 戴習(xí)林, 劉潔, 李晶晶, 等. 羅氏沼蝦種群SSR分析中樣本量及標(biāo)記量對遺傳多樣性指標(biāo)的影響[J]. 水產(chǎn)學(xué)報, 2017, 41(7): 1083-1095.

    DAI Xilin, LIU Jie, LI Jingjing, et al. Effects of sample size and loci number on genetic diversity index inwith microsatellite markers[J]. Journal of Fisheries of China, 2017, 41(7): 1083-1095.

    [48] 李歐, 趙瑩瑩, 郭娜, 等. 草魚種群SSR分析中樣本量及標(biāo)記數(shù)量對遺傳多度的影響[J]. 動物學(xué)研究, 2009, 30(2): 121-130.

    LI Ou, ZHAO Yingying, GUO Na, et al. Effects of sample size and loci number on genetic diversity in wild population of grass carp revealed by SSR[J]. Zoological Research, 2009, 30(2): 121-130.

    [49] HILLEN J E J, COSCIA I, VANDEPUTTE M, et al. Estimates of genetic variability and inbreeding in experimentally selected populations of European sea bass[J/OL]. Aquaculture, 2017, 479(March): 742-749.

    [50] NEI M. Estimation of average heterozygosity and genetic distance from a small number of individuals[J]. Genetics, 1978, 89(3): 583-590.

    [51] BEAUMONT A, BOUDRY P, HOARE K. Biote-chnology and genetics in fisheries and aquaculture[M]. Singapore: Wiley-Blackwell, 2010.

    [52] LIND C E, EVANS B S, KNAUER J, et al. Decreased genetic diversity and a reduced effective population size in cultured silver-lipped pearl oysters()[J]. Aquaculture, 2009, 286(1/2): 12-19.

    [53] RHODE C, MADUNA S N, ROODT-WILDING R, et al. Comparison of population genetic estimates amon-gst wild, F1 and F2 cultured abalone ()[J]. Animal Genetics, 2014, 45(3): 456-459.

    Microsatellite analysis on genetic diversity and genetic structure of three successive selected generations of

    LIANG Yuan1, 3, FU Jing-qiang2, 3, SHEN Ming-hui4, LUO Xuan1, 3, YOU Wei-wei1, 3, KE Cai-huan1, 3

    (1. College of Ocean and Earth Sciences, Xiamen University, Xiamen 361102, China; 2. College of the Environment and Ecology, Xiamen University, Xiamen 361102, China; 3. Fujian Key Laboratory of Genetics and Breeding of Marine Organisms, Xiamen University, Xiamen 361102, China; 4. Hainan Academy of Ocean and Fisheries Sciences, Hainan Provincial Key Laboratory of Tropical Maricultural Technologies, Haikou 571126, China)

    We analyzed the genetic diversity and genetic structure among three successive selected generations (TF4–TF6) ofusing 12 polymorphic microsatellite markers. Polymorphism information content () was 0.730, 0.717, and 0.708, showing relatively high polymorphisms (>0.5). Additionally, number of alleles (), effective number of alleles (e), observed heterozygosity (o), expected heterozygosity (e), and Shannon’s information index () ranged from 11.667 to 12.583, 6.334 to 6.915, 0.658 to 0.672, 0.737 to 0.787, and 1.837 to 2.003, respectively, suggesting an overall decreasing but not significant trend (0.05). Furthermore, analysis of molecular variance (AMOVA) showed that genetic variation among generations only accounted for 0.95% of the total variation, whilestranged from 0.00569 to 0.01324, belonging to the low differentiation level (st0.05). Moreover, principal coordinates analysis (PCoA) and STRUCTURE analysis illustrated a similar genetic background and genetic structure among the three generations. These results suggested that the three generations maintained a high level of genetic diversity after consecutive generations of artificial breeding, with a low level of genetic variation and differentiation and a stable genetic structure. The results also indicated that a highly selective potential still existed in the breeding populations of..

    ; selective population; microsatellite marker; genetic diversity; genetic structure

    Oct. 7, 2021

    S917.4

    A

    1000-3096(2022)10-0085-09

    10.11759/hykx20211007003

    2021-10-07;

    2022-01-13

    海南省重點研發(fā)項目(ZDYF2022XDNY234); 福建省科技計劃項目高校產(chǎn)學(xué)合作項目(2020N5001); 國家自然科學(xué)基金資助項目(32202900); 海南省熱帶海水養(yǎng)殖技術(shù)重點實驗室項目(TMTOF-202003)

    [Hainan Province Science and Technology Special Fund, No. ZDYF2022XDNY234; Fujian Provincial S & T Project, No. 2020N5001; National Natural Science Foundation of China, No. 32202900; The Project of Hainan Provincial Key Laboratory of Tropical Maricultural Technologies, No. TMTOF-202003]

    梁園(1998—), 女,滿族, 遼寧鳳城人, 碩士, 主要從事東風(fēng)螺遺傳育種研究, E-mail: liangyuan@stu.xmu.edu.cn; 柯才煥(1962—),通信作者, E-mail: chke@xmu.edu.cn

    (本文編輯: 譚雪靜)

    猜你喜歡
    分析
    禽大腸桿菌病的分析、診斷和防治
    隱蔽失效適航要求符合性驗證分析
    電力系統(tǒng)不平衡分析
    電子制作(2018年18期)2018-11-14 01:48:24
    電力系統(tǒng)及其自動化發(fā)展趨勢分析
    經(jīng)濟危機下的均衡與非均衡分析
    對計劃生育必要性以及其貫徹實施的分析
    GB/T 7714-2015 與GB/T 7714-2005對比分析
    出版與印刷(2016年3期)2016-02-02 01:20:11
    網(wǎng)購中不良現(xiàn)象分析與應(yīng)對
    中西醫(yī)結(jié)合治療抑郁癥100例分析
    偽造有價證券罪立法比較分析
    国产美女午夜福利| 热99在线观看视频| 全区人妻精品视频| 黄色配什么色好看| 亚洲欧美日韩无卡精品| 欧美三级亚洲精品| 国产免费福利视频在线观看| 黄色日韩在线| 你懂的网址亚洲精品在线观看| 夫妻性生交免费视频一级片| 99re6热这里在线精品视频| 特级一级黄色大片| 亚洲最大成人中文| 午夜福利在线观看免费完整高清在| 最近中文字幕2019免费版| 全区人妻精品视频| 最新中文字幕久久久久| 大片免费播放器 马上看| 亚洲av电影不卡..在线观看| 国产成人一区二区在线| 神马国产精品三级电影在线观看| 欧美性猛交╳xxx乱大交人| 国产精品熟女久久久久浪| 禁无遮挡网站| 成人午夜高清在线视频| 91精品伊人久久大香线蕉| 亚洲综合精品二区| 免费看av在线观看网站| 亚洲国产精品国产精品| 亚洲av在线观看美女高潮| 免费人成在线观看视频色| 日本-黄色视频高清免费观看| 久久久久久伊人网av| 丝袜美腿在线中文| 久久久久久久久久久丰满| 国产毛片a区久久久久| 日本与韩国留学比较| 国产高清国产精品国产三级 | 国产黄片视频在线免费观看| 如何舔出高潮| av福利片在线观看| 黑人高潮一二区| 97热精品久久久久久| 狠狠精品人妻久久久久久综合| 国产成人freesex在线| 国产精品福利在线免费观看| 日韩av在线大香蕉| 尤物成人国产欧美一区二区三区| 夜夜看夜夜爽夜夜摸| 高清日韩中文字幕在线| 老司机影院成人| 99热这里只有是精品50| 日韩欧美 国产精品| 午夜激情久久久久久久| 青春草国产在线视频| 18+在线观看网站| 久久久亚洲精品成人影院| 高清av免费在线| 国产淫片久久久久久久久| 中文字幕久久专区| 男的添女的下面高潮视频| 精品久久国产蜜桃| 超碰97精品在线观看| 色5月婷婷丁香| 亚洲精品,欧美精品| 乱人视频在线观看| 麻豆av噜噜一区二区三区| 久久久久免费精品人妻一区二区| eeuss影院久久| 免费观看无遮挡的男女| 国产高清三级在线| 亚洲第一区二区三区不卡| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 一个人观看的视频www高清免费观看| 亚洲国产成人一精品久久久| 九九在线视频观看精品| 夫妻午夜视频| 欧美3d第一页| 精品国产一区二区三区久久久樱花 | 国产在线一区二区三区精| 国产成人精品福利久久| 欧美精品国产亚洲| 秋霞在线观看毛片| 大香蕉久久网| av又黄又爽大尺度在线免费看| 亚洲最大成人手机在线| 日韩在线高清观看一区二区三区| 亚洲自拍偷在线| 观看美女的网站| 亚洲av不卡在线观看| 成人亚洲精品一区在线观看 | 国产有黄有色有爽视频| 国产亚洲午夜精品一区二区久久 | 看黄色毛片网站| 亚洲精品久久久久久婷婷小说| 乱人视频在线观看| 久久草成人影院| 免费大片黄手机在线观看| 免费观看a级毛片全部| 欧美zozozo另类| av线在线观看网站| 22中文网久久字幕| 日日干狠狠操夜夜爽| 精品99又大又爽又粗少妇毛片| 一个人看的www免费观看视频| 精品一区二区免费观看| 国产精品久久久久久精品电影小说 | 人妻少妇偷人精品九色| 亚洲无线观看免费| 精品人妻视频免费看| 国产黄色小视频在线观看| 简卡轻食公司| 国产精品无大码| 自拍偷自拍亚洲精品老妇| 91aial.com中文字幕在线观看| 精品人妻视频免费看| 少妇丰满av| 日韩大片免费观看网站| or卡值多少钱| 亚洲国产最新在线播放| 男人爽女人下面视频在线观看| 久久人人爽人人片av| 成人综合一区亚洲| 天堂av国产一区二区熟女人妻| 欧美激情国产日韩精品一区| av福利片在线观看| 国产真实伦视频高清在线观看| 97人妻精品一区二区三区麻豆| 久久精品熟女亚洲av麻豆精品 | av在线亚洲专区| 国产探花极品一区二区| 亚洲欧美一区二区三区黑人 | 亚洲欧洲国产日韩| 亚洲av成人av| 国产黄片视频在线免费观看| 午夜精品在线福利| 中文字幕免费在线视频6| 99久久精品一区二区三区| 成人亚洲精品一区在线观看 | 九九爱精品视频在线观看| 国内精品美女久久久久久| 小蜜桃在线观看免费完整版高清| 在线a可以看的网站| 三级男女做爰猛烈吃奶摸视频| 久久久久国产网址| 国产熟女欧美一区二区| 国产伦精品一区二区三区视频9| 国产视频首页在线观看| 亚洲四区av| 久久久久网色| 日韩av免费高清视频| 久久99精品国语久久久| 亚洲精品第二区| 欧美xxⅹ黑人| 国产精品久久久久久久电影| 91精品一卡2卡3卡4卡| 美女脱内裤让男人舔精品视频| 亚洲性久久影院| 成人毛片60女人毛片免费| 国产一区亚洲一区在线观看| 国产淫片久久久久久久久| av线在线观看网站| 国产熟女欧美一区二区| 亚洲国产欧美在线一区| 大片免费播放器 马上看| 亚洲精品456在线播放app| 日韩,欧美,国产一区二区三区| 亚洲精品视频女| 久久久久久国产a免费观看| 最后的刺客免费高清国语| 亚洲精品456在线播放app| 看免费成人av毛片| 亚洲怡红院男人天堂| 女人十人毛片免费观看3o分钟| 免费看av在线观看网站| 久久久久久久久久久免费av| 免费不卡的大黄色大毛片视频在线观看 | 亚洲国产精品专区欧美| 国产老妇伦熟女老妇高清| 亚洲国产精品国产精品| 日本免费在线观看一区| 色综合站精品国产| 女人被狂操c到高潮| 亚洲国产日韩欧美精品在线观看| 精品久久久久久久久亚洲| 日本爱情动作片www.在线观看| 非洲黑人性xxxx精品又粗又长| 亚洲av中文字字幕乱码综合| 深夜a级毛片| 日日干狠狠操夜夜爽| 亚洲精品456在线播放app| 波野结衣二区三区在线| 国产高清国产精品国产三级 | 亚洲国产精品成人久久小说| 99re6热这里在线精品视频| 26uuu在线亚洲综合色| 尤物成人国产欧美一区二区三区| 人妻系列 视频| 亚洲精品自拍成人| 精品不卡国产一区二区三区| 免费观看精品视频网站| 熟妇人妻不卡中文字幕| 蜜臀久久99精品久久宅男| 天天一区二区日本电影三级| 欧美 日韩 精品 国产| 午夜激情福利司机影院| 最近2019中文字幕mv第一页| 亚洲国产色片| 春色校园在线视频观看| av在线播放精品| 嫩草影院入口| av国产久精品久网站免费入址| 天天躁日日操中文字幕| 婷婷色综合www| 国产成人精品久久久久久| av播播在线观看一区| 六月丁香七月| 免费看不卡的av| 男人舔女人下体高潮全视频| 国产成人aa在线观看| 日韩国内少妇激情av| 七月丁香在线播放| 婷婷色av中文字幕| 亚洲精品aⅴ在线观看| 在线观看美女被高潮喷水网站| 别揉我奶头 嗯啊视频| 中文天堂在线官网| 精品久久久久久久末码| 午夜免费男女啪啪视频观看| 一级毛片 在线播放| 97精品久久久久久久久久精品| 九草在线视频观看| 舔av片在线| 久久精品国产鲁丝片午夜精品| 22中文网久久字幕| 日韩欧美精品免费久久| 久久久久久久久久成人| 天堂√8在线中文| 男人狂女人下面高潮的视频| 国产女主播在线喷水免费视频网站 | 老司机影院成人| 国产午夜精品论理片| www.色视频.com| 三级毛片av免费| 网址你懂的国产日韩在线| 一边亲一边摸免费视频| 直男gayav资源| 日韩精品青青久久久久久| 91久久精品国产一区二区成人| 午夜精品一区二区三区免费看| 国产 一区精品| 国产淫片久久久久久久久| 伦理电影大哥的女人| 一边亲一边摸免费视频| 国产 一区 欧美 日韩| 日本与韩国留学比较| 99久久精品国产国产毛片| 伦精品一区二区三区| 亚洲国产欧美在线一区| 中文字幕人妻熟人妻熟丝袜美| 免费高清在线观看视频在线观看| 日韩亚洲欧美综合| 国产高清不卡午夜福利| 女人十人毛片免费观看3o分钟| 天堂av国产一区二区熟女人妻| 婷婷六月久久综合丁香| 亚洲国产成人一精品久久久| 国产精品国产三级专区第一集| 丰满少妇做爰视频| 久久鲁丝午夜福利片| 久久久久久久午夜电影| 国产一区有黄有色的免费视频 | 亚洲精品第二区| 国产精品久久久久久久电影| 18禁在线无遮挡免费观看视频| kizo精华| 精品久久久久久久久av| av在线老鸭窝| 99久久精品热视频| 极品教师在线视频| 亚洲真实伦在线观看| 99久久精品国产国产毛片| 99久久九九国产精品国产免费| 97在线视频观看| 天堂影院成人在线观看| 日韩精品有码人妻一区| 亚洲精品456在线播放app| 亚洲激情五月婷婷啪啪| 禁无遮挡网站| 精品国产露脸久久av麻豆 | 久久久久久久久中文| 日日撸夜夜添| 成年av动漫网址| 2022亚洲国产成人精品| 在线免费十八禁| 国产有黄有色有爽视频| 天天躁日日操中文字幕| 白带黄色成豆腐渣| 69人妻影院| 天堂av国产一区二区熟女人妻| 国产精品国产三级国产av玫瑰| 狠狠精品人妻久久久久久综合| 免费看a级黄色片| 国产精品熟女久久久久浪| 99久久九九国产精品国产免费| 18禁在线播放成人免费| 人人妻人人看人人澡| 午夜激情久久久久久久| 又黄又爽又刺激的免费视频.| eeuss影院久久| 日本黄色片子视频| 亚洲精品,欧美精品| 久99久视频精品免费| 美女脱内裤让男人舔精品视频| 国精品久久久久久国模美| 免费av观看视频| 天堂√8在线中文| 最近最新中文字幕大全电影3| eeuss影院久久| 久久精品久久久久久久性| 欧美精品一区二区大全| 五月天丁香电影| 成人亚洲欧美一区二区av| 亚洲精品久久午夜乱码| 国产成人一区二区在线| 我要看日韩黄色一级片| 午夜免费男女啪啪视频观看| 国产精品爽爽va在线观看网站| 欧美一级a爱片免费观看看| 国产伦精品一区二区三区视频9| 美女高潮的动态| 在线观看av片永久免费下载| 欧美xxxx性猛交bbbb| 永久网站在线| 亚洲四区av| av天堂中文字幕网| 亚洲精品影视一区二区三区av| 久久久久久久久久成人| 偷拍熟女少妇极品色| 亚洲精品一二三| 夫妻午夜视频| 精品久久久久久久人妻蜜臀av| 秋霞在线观看毛片| 国产一级毛片七仙女欲春2| 亚洲精品国产av成人精品| 国产成人精品福利久久| 69av精品久久久久久| 美女被艹到高潮喷水动态| 3wmmmm亚洲av在线观看| 日本一二三区视频观看| 99热这里只有是精品50| 黄色日韩在线| 亚洲怡红院男人天堂| 能在线免费看毛片的网站| 真实男女啪啪啪动态图| 精品国产一区二区三区久久久樱花 | 亚州av有码| 成人漫画全彩无遮挡| 国产成人aa在线观看| 国产伦精品一区二区三区四那| 91av网一区二区| 超碰av人人做人人爽久久| 搡女人真爽免费视频火全软件| 丰满少妇做爰视频| av专区在线播放| 2021天堂中文幕一二区在线观| 午夜福利视频精品| 国产精品久久久久久久久免| av福利片在线观看| 97人妻精品一区二区三区麻豆| 中文精品一卡2卡3卡4更新| 99久久精品国产国产毛片| 丝瓜视频免费看黄片| 精品少妇黑人巨大在线播放| 国产 一区精品| 亚洲va在线va天堂va国产| 大又大粗又爽又黄少妇毛片口| 国产精品综合久久久久久久免费| a级毛片免费高清观看在线播放| 成人亚洲精品av一区二区| freevideosex欧美| 女人被狂操c到高潮| 国产午夜精品久久久久久一区二区三区| 97精品久久久久久久久久精品| 国内少妇人妻偷人精品xxx网站| 成人av在线播放网站| 国产乱人视频| 亚洲精品第二区| 亚洲av国产av综合av卡| 久久99热6这里只有精品| 日韩精品青青久久久久久| 三级国产精品片| 免费高清在线观看视频在线观看| 99久国产av精品| av国产免费在线观看| 欧美成人精品欧美一级黄| 日韩中字成人| 久久精品国产亚洲av天美| 亚洲av成人精品一二三区| 精品久久久久久久久久久久久| 亚洲国产精品专区欧美| 成人综合一区亚洲| 国产国拍精品亚洲av在线观看| 精品一区二区三卡| 成人特级av手机在线观看| 波多野结衣巨乳人妻| 我要看日韩黄色一级片| 九草在线视频观看| 一区二区三区乱码不卡18| 在线天堂最新版资源| 2021少妇久久久久久久久久久| 三级男女做爰猛烈吃奶摸视频| 69人妻影院| 国产精品不卡视频一区二区| 九九久久精品国产亚洲av麻豆| 精品一区二区三区视频在线| 十八禁网站网址无遮挡 | 2022亚洲国产成人精品| 日本爱情动作片www.在线观看| 真实男女啪啪啪动态图| 久久精品人妻少妇| 卡戴珊不雅视频在线播放| 丝瓜视频免费看黄片| 亚洲欧美一区二区三区黑人 | 国产精品一区二区三区四区免费观看| 精品久久久久久久久亚洲| 久久久久精品性色| 免费高清在线观看视频在线观看| 女的被弄到高潮叫床怎么办| 久久久成人免费电影| 国产亚洲av片在线观看秒播厂 | 高清毛片免费看| 久久久成人免费电影| 亚洲在线自拍视频| 国产亚洲午夜精品一区二区久久 | 国产一区二区三区av在线| 国产91av在线免费观看| 国产欧美另类精品又又久久亚洲欧美| 中文字幕av在线有码专区| 欧美性感艳星| 人人妻人人澡人人爽人人夜夜 | 伦理电影大哥的女人| 一区二区三区乱码不卡18| 永久免费av网站大全| 51国产日韩欧美| 久久99精品国语久久久| 亚洲美女搞黄在线观看| 亚州av有码| 国产在线男女| 国精品久久久久久国模美| 国产av不卡久久| 国产免费又黄又爽又色| 激情五月婷婷亚洲| 久久精品国产亚洲av涩爱| 大片免费播放器 马上看| 亚洲欧美一区二区三区黑人 | 国产精品久久久久久精品电影| 欧美一级a爱片免费观看看| 综合色av麻豆| 免费看光身美女| 亚洲精品影视一区二区三区av| 欧美xxⅹ黑人| 韩国高清视频一区二区三区| 久久精品夜夜夜夜夜久久蜜豆| 乱码一卡2卡4卡精品| 色哟哟·www| 日日摸夜夜添夜夜添av毛片| 卡戴珊不雅视频在线播放| 日本-黄色视频高清免费观看| 国产精品一区www在线观看| 黄片无遮挡物在线观看| 午夜久久久久精精品| 夫妻午夜视频| 九草在线视频观看| 一本一本综合久久| 街头女战士在线观看网站| 国产91av在线免费观看| 久久人人爽人人爽人人片va| 内射极品少妇av片p| 亚洲婷婷狠狠爱综合网| 国产永久视频网站| 一级av片app| 色吧在线观看| 久久久久久久久中文| 成人午夜高清在线视频| 亚洲国产高清在线一区二区三| 精品人妻熟女av久视频| 看非洲黑人一级黄片| 成年免费大片在线观看| 又大又黄又爽视频免费| 久久精品熟女亚洲av麻豆精品 | 丝瓜视频免费看黄片| 男女下面进入的视频免费午夜| 狠狠精品人妻久久久久久综合| videos熟女内射| 日韩av不卡免费在线播放| 国产精品熟女久久久久浪| 国产有黄有色有爽视频| 观看美女的网站| 日日摸夜夜添夜夜添av毛片| 成人无遮挡网站| ponron亚洲| 亚洲国产成人一精品久久久| 天堂√8在线中文| 一个人看视频在线观看www免费| 一区二区三区免费毛片| 精品人妻偷拍中文字幕| 亚洲美女搞黄在线观看| 国产精品日韩av在线免费观看| 不卡视频在线观看欧美| 好男人在线观看高清免费视频| 自拍偷自拍亚洲精品老妇| av黄色大香蕉| 精品一区在线观看国产| 欧美不卡视频在线免费观看| 久久久久网色| 亚洲精品日韩av片在线观看| 又黄又爽又刺激的免费视频.| 欧美性感艳星| 亚洲电影在线观看av| 少妇熟女欧美另类| 国产成人91sexporn| 18+在线观看网站| av国产久精品久网站免费入址| 亚洲精品视频女| 国产在线一区二区三区精| 丰满乱子伦码专区| 国产精品.久久久| 午夜爱爱视频在线播放| 日韩精品有码人妻一区| 波多野结衣巨乳人妻| av卡一久久| 亚洲精品视频女| 久久久久久久久久久丰满| 男女边吃奶边做爰视频| 99久久九九国产精品国产免费| 久99久视频精品免费| 久久精品夜夜夜夜夜久久蜜豆| 天堂av国产一区二区熟女人妻| 91精品伊人久久大香线蕉| 91午夜精品亚洲一区二区三区| 中文字幕免费在线视频6| 国产亚洲午夜精品一区二区久久 | 亚洲激情五月婷婷啪啪| 精品亚洲乱码少妇综合久久| 精品一区在线观看国产| 欧美区成人在线视频| 青春草国产在线视频| 丰满人妻一区二区三区视频av| 97超碰精品成人国产| 欧美xxxx黑人xx丫x性爽| 国产av在哪里看| 真实男女啪啪啪动态图| 日本爱情动作片www.在线观看| 国产女主播在线喷水免费视频网站 | 永久网站在线| 欧美日韩亚洲高清精品| 国产亚洲午夜精品一区二区久久 | 国产精品一及| 欧美激情国产日韩精品一区| 亚洲最大成人手机在线| 久久精品国产自在天天线| a级一级毛片免费在线观看| 久久午夜福利片| 午夜福利在线在线| 国产精品麻豆人妻色哟哟久久 | 夜夜看夜夜爽夜夜摸| 男人狂女人下面高潮的视频| 少妇的逼好多水| 黄色日韩在线| 韩国高清视频一区二区三区| 777米奇影视久久| 欧美不卡视频在线免费观看| 一级黄片播放器| 日日撸夜夜添| 十八禁网站网址无遮挡 | 国产不卡一卡二| 国产午夜精品一二区理论片| 亚洲精品国产av成人精品| av免费观看日本| 免费人成在线观看视频色| 精品久久久久久久末码| 中文欧美无线码| 成人美女网站在线观看视频| 久久97久久精品| 22中文网久久字幕| 午夜久久久久精精品| 久久99热这里只频精品6学生| 国产永久视频网站| 中文字幕人妻熟人妻熟丝袜美| 搡老妇女老女人老熟妇| 亚洲精品aⅴ在线观看| 别揉我奶头 嗯啊视频| 青青草视频在线视频观看| 中文天堂在线官网| 又爽又黄a免费视频| 一夜夜www| 男女边摸边吃奶| 午夜久久久久精精品| 91在线精品国自产拍蜜月| 午夜福利高清视频| 亚洲精品一区蜜桃| 搡老妇女老女人老熟妇| 日韩av在线免费看完整版不卡| 国产高清三级在线| av在线播放精品| 久久久午夜欧美精品| 国产精品麻豆人妻色哟哟久久 | 国产爱豆传媒在线观看| 久久久久久久久久久丰满| 一个人免费在线观看电影| 22中文网久久字幕| 最近最新中文字幕大全电影3| 91久久精品国产一区二区成人| 日本与韩国留学比较|