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

    浙江長(zhǎng)興揚(yáng)子鱷種群遺傳多樣性研究

    2022-08-13 03:13:14方黎明付春正鄒衛(wèi)強(qiáng)任大斌徐菊敏韓群花萬(wàn)秋紅方盛國(guó)
    野生動(dòng)物學(xué)報(bào) 2022年3期
    關(guān)鍵詞:揚(yáng)子鱷多態(tài)微衛(wèi)星

    方黎明 付春正 李 慧 鄒衛(wèi)強(qiáng) 任大斌 徐菊敏 韓群花 萬(wàn)秋紅 方盛國(guó)*

    (1.長(zhǎng)興尹家邊揚(yáng)子鱷自然保護(hù)區(qū),長(zhǎng)興,313100;2.浙江大學(xué)生命科學(xué)學(xué)院,生命系統(tǒng)穩(wěn)態(tài)與保護(hù)教育部重點(diǎn)實(shí)驗(yàn)室(浙江大學(xué)),國(guó)家瀕危野生動(dòng)植物種質(zhì)基因保護(hù)中心(浙江大學(xué)),杭州,310058)

    生物多樣性是指地球上生存的所有動(dòng)物、植物、微生物和這些生物所擁有的所有遺傳物質(zhì),以及由這些生物構(gòu)成的生態(tài)系統(tǒng),即由遺傳多樣性、物種多樣性和生態(tài)系統(tǒng)多樣性三部分組成,其中遺傳多樣性是生物多樣性的重要組成部分,是自然界中物種或種群內(nèi)各種不同等位基因的總和,表現(xiàn)為在生物體遺傳上的千差萬(wàn)別[1-5]。遺傳多樣性能夠反應(yīng)物種或種群的遺傳變異水平,生物的遺傳多樣性越高,其對(duì)環(huán)境的適應(yīng)能力就越強(qiáng),生存能力也就越強(qiáng),因此遺傳多樣性是衡量物種或種群生存、進(jìn)化和環(huán)境條件適應(yīng)的關(guān)鍵指標(biāo)[6-7]。此外,遺傳多樣性也是保護(hù)生物學(xué)研究的重要內(nèi)容之一,它可以在一定范圍內(nèi)決定物種的數(shù)量及分布范圍,如果缺乏對(duì)生物的遺傳多樣性研究,就無(wú)法準(zhǔn)確地認(rèn)識(shí)物種種群的地理分布格局、進(jìn)化歷史、適應(yīng)機(jī)理、物種數(shù)量增長(zhǎng)和滅絕等問(wèn)題[8-9]。

    微衛(wèi)星標(biāo)記因多態(tài)性豐富、共顯性遺傳、易于檢測(cè)且有效性高和重復(fù)性好等特點(diǎn),在分析物種或種群的遺傳多樣性及系統(tǒng)進(jìn)化關(guān)系中表現(xiàn)出巨大的優(yōu)勢(shì),被認(rèn)為是研究物種或種群遺傳變異的最優(yōu)分子標(biāo)記,廣泛應(yīng)用于羊、馬、牛、豬和雞等家畜、家禽的遺傳多樣性研究[10-18]。此外,微衛(wèi)星標(biāo)記檢測(cè)因?qū)δ0錎NA要求低,且操作過(guò)程簡(jiǎn)捷、高效和自動(dòng)化,在對(duì)瀕危野生動(dòng)物如揚(yáng)子鱷(Alligatorsinensis)、東北馬鹿(Cervuscanadensis)、朱鹮(Nipponianippon)、大熊貓(Ailuropodamelanoleuca)、虎(Pantheratigris)、豹貓(Prionailurusbengalensis)、丹頂鶴(Grusjaponensis)、東北梅花鹿(Cervusnipponhortulorum)和西伯利亞狍(Capreoluspygargus)[19-27]等的保護(hù)遺傳學(xué)研究中,該檢測(cè)技術(shù)的應(yīng)用越來(lái)越廣泛。

    揚(yáng)子鱷是中國(guó)特有的珍稀瀕危物種,屬爬行綱(Reptilia),鱷形目(Crocodylia),鼉科(Alligatoridae),短吻鱷屬(Alligator)[28-29]。該物種曾廣泛分布于我國(guó)大部分區(qū)域,南起海南儋州,北到新疆呼圖壁;西起新疆呼圖壁,東至浙江余姚,均有揚(yáng)子鱷分布[30]。但隨著氣候變遷、人類對(duì)揚(yáng)子鱷的過(guò)度捕殺及人類對(duì)揚(yáng)子鱷棲息環(huán)境的破壞,導(dǎo)致?lián)P子鱷分布范圍大幅縮小,現(xiàn)今揚(yáng)子鱷主要分布于浙江長(zhǎng)興和安徽宣城[30-31]。浙江長(zhǎng)興尹家邊揚(yáng)子鱷保護(hù)區(qū)始建于1979年,最初的種群由11條奠基揚(yáng)子鱷個(gè)體組成,截至2019年9月,該保護(hù)區(qū)揚(yáng)子鱷數(shù)量已達(dá)6 691條[32-33]。浙江長(zhǎng)興尹家邊揚(yáng)子鱷保護(hù)區(qū)由核心區(qū)和野放區(qū)兩部分組成,核心區(qū)為揚(yáng)子鱷研究繁育中心,即最初建設(shè)的揚(yáng)子鱷保護(hù)區(qū)所在區(qū)域;野放區(qū)在“國(guó)家Ⅱ期工程”——浙江省揚(yáng)子鱷放歸自然建設(shè)項(xiàng)目的資助下于2009年開(kāi)始建設(shè),2011年完成[32,34]。2012年4月15日,保護(hù)區(qū)首次從核心區(qū)甄選120條揚(yáng)子鱷放歸至野放區(qū),作為野放區(qū)的奠基種群[34-35]。自此,保護(hù)區(qū)每年都從核心區(qū)內(nèi)挑選一定數(shù)量的7~8齡的揚(yáng)子鱷放歸至野放區(qū),截至2019年7月,已有960條揚(yáng)子鱷被放歸至野放區(qū)[34,36-37]。隨著野放區(qū)揚(yáng)子鱷的繁殖,野放區(qū)揚(yáng)子鱷已經(jīng)繁衍發(fā)育成一個(gè)數(shù)量較大的種群。

    已有研究[32,38-40]大多利用微衛(wèi)星標(biāo)記對(duì)浙江長(zhǎng)興尹家邊揚(yáng)子鱷保護(hù)區(qū)揚(yáng)子鱷種群個(gè)體進(jìn)行親子鑒定,或?qū)ΡWo(hù)區(qū)揚(yáng)子鱷種群遺傳多樣性進(jìn)行整體性分析,并未對(duì)核心區(qū)和野放區(qū)揚(yáng)子鱷種群遺傳多樣性進(jìn)行單獨(dú)分析和比較,導(dǎo)致保護(hù)區(qū)核心區(qū)和野放區(qū)揚(yáng)子鱷種群遺傳多樣性差異研究尚有不足。因此,本研究通過(guò)隨機(jī)采集浙江長(zhǎng)興尹家邊揚(yáng)子鱷保護(hù)區(qū)核心區(qū)和野放區(qū)的成年揚(yáng)子鱷血液樣本,利用微衛(wèi)星標(biāo)記技術(shù),對(duì)核心區(qū)和野放區(qū)種群進(jìn)行遺傳多樣性分析,從而了解核心區(qū)和野放區(qū)種群的遺傳多樣性及種群遺傳分化情況,為核心區(qū)揚(yáng)子鱷放歸至野放區(qū)提供分子水平上的理論依據(jù)。

    1 材料與方法

    1.1 樣品采集與DNA提取

    1.2 微衛(wèi)星標(biāo)記及引物信息

    13個(gè)多態(tài)性較高的微衛(wèi)星標(biāo)記信息見(jiàn)表1[38-39,43]。引物均由上海捷瑞生物工程有限公司合成,上游引物的5′端進(jìn)行熒光修飾(FAM、HEX或TAMRA)。

    表1 微衛(wèi)星標(biāo)記信息

    續(xù)表1

    1.3 PCR擴(kuò)增及STR擴(kuò)增產(chǎn)物掃描分析

    使用20.0 μL PCR擴(kuò)增反應(yīng)體系,包括2×TaqMaster Mix(上海捷瑞生物工程有限公司)10.0 μL,ddH2O 8.0 μL,模板DNA 1.0 μL,上、下游引物各0.5 μL。PCR反應(yīng)程序:95 ℃預(yù)變性5 min;95 ℃變性30 s,在各引物對(duì)應(yīng)的退火溫度下退火30 s,72 ℃延伸30 s,擴(kuò)增34個(gè)循環(huán);72 ℃終延伸5 min,4 ℃保存。PCR擴(kuò)增產(chǎn)物經(jīng)1%瓊脂糖凝膠電泳(120 V,30 min)檢測(cè),將擴(kuò)增成功的PCR產(chǎn)物送至Invitrogen(上海英俊)生物技術(shù)有限公司進(jìn)行STR擴(kuò)增產(chǎn)物掃描分析,確定擴(kuò)增產(chǎn)物片段大小。

    1.4 數(shù)據(jù)統(tǒng)計(jì)分析

    使用GeneMarker v1.5讀取STR擴(kuò)增產(chǎn)物掃描分析結(jié)果,確定微衛(wèi)星標(biāo)記片段大??;使用Excel Microsatellite Toolkit version 3.1計(jì)算等位基因頻率、等位基因大小范圍、多態(tài)信息含量和數(shù)據(jù)格式轉(zhuǎn)化;等位基因數(shù)量、有效等位基因數(shù)量、觀測(cè)雜合度、期望雜合度、香農(nóng)信息指數(shù)和F-統(tǒng)計(jì)值用POPGENE 1.32計(jì)算;分子方差分析(AMOVA)、FST矩陣計(jì)算和基因流計(jì)算由Arlequin 3.5.2.2完成。

    2 結(jié)果與分析

    2.1 等位基因分布及頻率

    核心區(qū)和野放區(qū)130條揚(yáng)子鱷在13個(gè)微衛(wèi)星標(biāo)記中共檢測(cè)到31個(gè)等位基因,其中含有2個(gè)稀有等位基因(基因頻率<5.00%[44])。13個(gè)微衛(wèi)星標(biāo)記在核心區(qū)和野放區(qū)揚(yáng)子鱷種群中均表現(xiàn)出多態(tài)性,微衛(wèi)星標(biāo)記CXA-9、CXA-34、CXA-35、CXA-43和CXA-136均具有3個(gè)等位基因,其余8個(gè)微衛(wèi)星標(biāo)記均具有2個(gè)等位基因。核心區(qū)和野放區(qū)揚(yáng)子鱷種群在13個(gè)微衛(wèi)星標(biāo)記上均具有相同的31個(gè)等位基因,平均等位基因數(shù)量均為2.38個(gè)。

    13個(gè)微衛(wèi)星標(biāo)記的等位基因頻率在保護(hù)區(qū)核心區(qū)和野放區(qū)揚(yáng)子鱷種群間存在差異。CXA-4微衛(wèi)星標(biāo)記的等位基因共有2個(gè),大小分別為123、128 bp,核心區(qū)揚(yáng)子鱷種群基因頻率較高的等位基因是128 bp(64.62%),而野放區(qū)為123 bp(52.31%);CXA-41微衛(wèi)星標(biāo)記的等位基因共有2個(gè),大小分別為156、160 bp,核心區(qū)揚(yáng)子鱷種群基因頻率較高的等位基因是160 bp(57.69%),而野放區(qū)為156 bp(53.08%);CXA-43微衛(wèi)星標(biāo)記的等位基因共有3個(gè),大小分別為113、117、129 bp,核心區(qū)揚(yáng)子鱷種群基因頻率較高的等位基因是為113 bp(56.92%),而野放區(qū)為117 bp(49.23%)。其余10個(gè)微衛(wèi)星標(biāo)記CXA-6、CXA-9、CXA-34、CXA-35、CXA-37、CXA-39、CXA-135、CXA-136、CXA-138和CXA-142,雖然每個(gè)微衛(wèi)星標(biāo)記的等位基因在核心區(qū)和野放區(qū)揚(yáng)子鱷種群間的基因頻率存在差異,但每個(gè)微衛(wèi)星標(biāo)記的較高基因頻率的等位基因是一致的,如CXA-35微衛(wèi)星標(biāo)記的等位基因共有3個(gè),大小分別為110、114、118 bp,核心區(qū)和野放區(qū)揚(yáng)子鱷種群基因頻率較高的等位基因均為114 bp,頻率分別為87.69%和90.00%。此外,基因頻率<5.00%的等位基因包括微衛(wèi)星標(biāo)記CXA-9的117 bp,核心區(qū)和野放區(qū)揚(yáng)子鱷種群中該基因的頻率分別為3.08%和0.77%;微衛(wèi)星標(biāo)記CXA-35的118 bp,核心區(qū)和野放區(qū)揚(yáng)子鱷種群該基因的頻率分別為3.08%和4.62%(表2)。

    表2 核心區(qū)和野放區(qū)揚(yáng)子鱷種群13個(gè)微衛(wèi)星標(biāo)記的等位基因長(zhǎng)度和頻率

    2.2 微衛(wèi)星標(biāo)記的遺傳多樣性

    保護(hù)區(qū)揚(yáng)子鱷種群中的各遺傳多樣性參數(shù)詳見(jiàn)表3。13個(gè)微衛(wèi)星標(biāo)記的等位基因數(shù)量為2或3個(gè),平均等位基因數(shù)量為2.38個(gè);有效等位基因數(shù)量為1.26~2.59個(gè),平均有效等位基因數(shù)量為1.88個(gè),平均等位基因數(shù)量與平均有效等位基因數(shù)量存在差異。有效等位基因數(shù)量最高的微衛(wèi)星標(biāo)記是CXA-34,為2.59個(gè);有效等位基因數(shù)量最低的微衛(wèi)星標(biāo)記是CXA-35,為1.26個(gè)。

    表3 保護(hù)區(qū)揚(yáng)子鱷種群13個(gè)微衛(wèi)星標(biāo)記的遺傳多樣性參數(shù)

    13個(gè)微衛(wèi)星標(biāo)記的觀測(cè)雜合度為0.21~0.75,平均觀測(cè)雜合度為0.50,觀測(cè)雜合度最高的微衛(wèi)星標(biāo)記為CXA-34;觀測(cè)雜合度最低的微衛(wèi)星標(biāo)記為CXA-9和CXA-35。13個(gè)微衛(wèi)星標(biāo)記的期望雜合度為0.20~0.62,平均期望雜合度為0.45,期望雜合度最高的微衛(wèi)星標(biāo)記為CXA-34;期望雜合度最低的微衛(wèi)星標(biāo)記為CXA-35。微衛(wèi)星標(biāo)記CXA-9和CXA-43的觀測(cè)雜合度低于其期望雜合度。

    13個(gè)微衛(wèi)星標(biāo)記的多態(tài)信息含量為0.19~0.54,平均多態(tài)信息含量為0.36,多態(tài)信息含量最高的微衛(wèi)星標(biāo)記為CXA-34;多態(tài)信息含量最低的微衛(wèi)星標(biāo)記為CXA-35。微衛(wèi)星標(biāo)記CXA-34和CXA-43為高度多態(tài)標(biāo)記(多態(tài)信息含量>0.50[16]);微衛(wèi)星標(biāo)記CXA-4、CXA-6、CXA-9、CXA-37、CXA-39、CXA-41、CXA-136、CXA-138和CXA-142為中度多態(tài)標(biāo)記(0.25<多態(tài)信息含量<0.50[16]);微衛(wèi)星標(biāo)記CXA-35和CXA-135為低度多態(tài)標(biāo)記(多態(tài)信息含量<0.25[16])。13個(gè)微衛(wèi)星標(biāo)記的香農(nóng)信息指數(shù)為0.42~1.01,平均香農(nóng)信息指數(shù)為0.68,香農(nóng)信息指數(shù)最高的微衛(wèi)星標(biāo)記為CXA-34,最低的為CXA-35。

    微衛(wèi)星標(biāo)記CXA-34的各遺傳多樣性參數(shù),包括有效等位基因數(shù)量、觀測(cè)雜合度、期望雜合度、多態(tài)信息含量和香農(nóng)信息指數(shù)的數(shù)值在13個(gè)微衛(wèi)星標(biāo)記中均最高,表明該微衛(wèi)星標(biāo)記的遺傳變異最大,遺傳多樣性最豐富;相反,微衛(wèi)星標(biāo)記CXA-35的各遺傳多樣性參數(shù)在13個(gè)微衛(wèi)星標(biāo)記中均最低,表明該微衛(wèi)星標(biāo)記的遺傳變異最小,遺傳多樣性最匱乏。

    2.3 核心區(qū)和野放區(qū)揚(yáng)子鱷種群遺傳多樣性

    核心區(qū)和野放區(qū)揚(yáng)子鱷種群具有相同的平均等位基因數(shù)量,但核心區(qū)揚(yáng)子鱷種群的平均有效等位基因數(shù)量高于野放區(qū);核心區(qū)揚(yáng)子鱷種群的平均觀測(cè)雜合度、平均期望雜合度、平均多態(tài)信息含量和平均香農(nóng)信息指數(shù)均高于野放區(qū),表明野放區(qū)揚(yáng)子鱷種群相較于核心區(qū)揚(yáng)子鱷種群的遺傳變異少,遺傳多樣性相對(duì)缺乏(表4)。

    表4 核心區(qū)和野放區(qū)揚(yáng)子鱷種群的遺傳多樣性參數(shù)

    2.4 保護(hù)區(qū)揚(yáng)子鱷種群的遺傳變異與分化

    在保護(hù)區(qū)揚(yáng)子鱷種群具有的總遺傳變異中,有96.83%的遺傳變異來(lái)源于各種群內(nèi)個(gè)體間的遺傳變異,僅有3.17%的遺傳變異來(lái)源于核心區(qū)和野放區(qū)種群間的遺傳變異,表明保護(hù)區(qū)種群內(nèi)的遺傳變異主要來(lái)源于個(gè)體之間的遺傳變異,核心區(qū)和野放區(qū)揚(yáng)子鱷種群間的遺傳變異較少,遺傳分化程度較低(表5)。

    表5 保護(hù)區(qū)揚(yáng)子鱷種群遺傳變異的分子方差分析結(jié)果

    13個(gè)微衛(wèi)星標(biāo)記在浙江長(zhǎng)興尹家邊揚(yáng)子鱷保護(hù)區(qū)揚(yáng)子鱷種群中的FIS值為-0.409 5(CXA-142)~0.289 9(CXA-9),平均值為-0.153 1,其中只有微衛(wèi)星標(biāo)記CXA-9的FIS值為正值,其余12個(gè)微衛(wèi)星標(biāo)記的FIS值均為負(fù)值;13個(gè)微衛(wèi)星標(biāo)記在保護(hù)區(qū)揚(yáng)子鱷種群中的FIT值為-0.387 6(CXA-41)~0.301 4(CXA-9),平均值為-0.130 0,其中微衛(wèi)星標(biāo)記CXA-9和CXA-43的FIT值為正值,其余11個(gè)微衛(wèi)星標(biāo)記的FIT值均為負(fù)值,且13個(gè)微衛(wèi)星標(biāo)記的FIS值均小于FIT值;13個(gè)微衛(wèi)星標(biāo)記在保護(hù)區(qū)揚(yáng)子鱷種群中的FST值為0.000 2(CXA-37、CXA-39和CXA-136)~0.056 7(CXA-135),平均值為0.020 0(表6)。

    13個(gè)微衛(wèi)星標(biāo)記在保護(hù)區(qū)核心區(qū)和野放區(qū)揚(yáng)子鱷種群中的FIS值分別為-0.607 3(CXA-142)~0.413 7(CXA-9)和-0.588 4(CXA-34)~0.0851(CXA-35);13個(gè)微衛(wèi)星標(biāo)記在保護(hù)區(qū)核心區(qū)和野放區(qū)揚(yáng)子鱷種群中的FIT值與其自身的FIS值相同,且13個(gè)微衛(wèi)星標(biāo)記在保護(hù)區(qū)核心區(qū)和野放區(qū)揚(yáng)子鱷種群中的FST值均為0(表6)。此外,保護(hù)區(qū)核心區(qū)與野放區(qū)的FST值為0.031 7,基因流為7.631 5(表7),表明保護(hù)區(qū)核心區(qū)和野放區(qū)揚(yáng)子鱷種群的遺傳分化程度較低。

    表6 13個(gè)微衛(wèi)星標(biāo)記的F-統(tǒng)計(jì)值

    表7 核心區(qū)和野放區(qū)揚(yáng)子鱷種群間的FST值和基因流

    3 討論

    在微衛(wèi)星標(biāo)記的等位基因中,若某一等位基因的頻率很高,那么該等位基因則很容易遺傳給子代;而若某一等位基因的頻率<5.00%,則稱為稀有等位基因,其在生物種群繁衍過(guò)程中,遺傳給子代的概率很小,很容易在種群繁衍過(guò)程中消失[26,44]。因此,在保護(hù)瀕危動(dòng)物種群繁衍過(guò)程中所發(fā)現(xiàn)的稀有等位基因,也是瀕危動(dòng)物遺傳多樣性研究的重要內(nèi)容[26]。本研究中的13個(gè)微衛(wèi)星標(biāo)記在浙江長(zhǎng)興尹家邊揚(yáng)子鱷保護(hù)區(qū)揚(yáng)子鱷種群中共檢測(cè)到稀有等位基因2個(gè),占所有等位基因的6.45%,分別為微衛(wèi)星標(biāo)記CXA-9,長(zhǎng)度為117 bp,該基因在核心區(qū)和野放區(qū)種群中的頻率分別為3.08%和0.77%;微衛(wèi)星標(biāo)記CXA-35,長(zhǎng)度為118 bp,該基因在核心區(qū)和野放區(qū)種群中的頻率分別為3.08%和4.62%。稀有等位基因是生物進(jìn)化的重要成果,同樣也是遺傳多樣性保護(hù)的重要內(nèi)容,特別是對(duì)于像浙江長(zhǎng)興尹家邊揚(yáng)子鱷保護(hù)區(qū)中人工繁殖的瀕危小種群,檢測(cè)到的2個(gè)稀有等位基因很可能是最初揚(yáng)子鱷種群中相對(duì)普遍存在的等位基因,但由于20世紀(jì)60年代的氣候變遷,人類對(duì)揚(yáng)子鱷的過(guò)度捕殺和人類對(duì)揚(yáng)子鱷棲息環(huán)境的破壞,造成了揚(yáng)子鱷種群的遺傳瓶頸以及在保護(hù)過(guò)程中人工選擇造成的基因頻率下降。因此,為了保護(hù)浙江長(zhǎng)興尹家邊揚(yáng)子鱷保護(hù)區(qū)揚(yáng)子鱷種群的遺傳多樣性,這些稀有等位基因更應(yīng)該得到重視[24,26,45]。在將來(lái)的揚(yáng)子鱷繁育保護(hù)過(guò)程中,研究人員應(yīng)該重點(diǎn)關(guān)注具有稀有等位基因的揚(yáng)子鱷個(gè)體,如有必要,可對(duì)其進(jìn)行集中飼養(yǎng)繁育,防止稀有等位基因在揚(yáng)子鱷隨機(jī)交配繁衍過(guò)程中丟失,造成保護(hù)區(qū)揚(yáng)子鱷種群遺傳多樣性的進(jìn)一步下降。此外,需對(duì)保護(hù)區(qū)揚(yáng)子鱷種群的遺傳多樣性進(jìn)行定時(shí)檢測(cè),監(jiān)測(cè)核心區(qū)和野放區(qū)揚(yáng)子鱷種群內(nèi)稀有等位基因動(dòng)態(tài)變化情況,如發(fā)現(xiàn)其中一個(gè)區(qū)域的揚(yáng)子鱷種群稀有等位基因頻率快速下降且呈消失趨勢(shì),應(yīng)及時(shí)在另一種群內(nèi)選取合適揚(yáng)子鱷個(gè)體補(bǔ)充至稀有等位基因即將消失的種群內(nèi)。

    衡量微衛(wèi)星標(biāo)記多態(tài)性的一個(gè)重要指標(biāo)是多態(tài)信息含量,當(dāng)一個(gè)標(biāo)記的多態(tài)信息含量>0.50時(shí)為高度多態(tài)標(biāo)記;0.25<多態(tài)信息含量<0.50時(shí)為中度多態(tài)標(biāo)記;多態(tài)信息含量<0.25時(shí)為低度多態(tài)標(biāo)記[6,46]。此外,微衛(wèi)星標(biāo)記的多態(tài)信息含量還與該標(biāo)記的可用性和使用效率相關(guān),微衛(wèi)星標(biāo)記的多態(tài)信息含量越大,在一個(gè)生物種群中該標(biāo)記的雜合子比例就越大,提供的遺傳信息就越多[46]。本研究使用的13個(gè)微衛(wèi)星標(biāo)記的多態(tài)信息含量為0.19~0.54,平均多態(tài)信息含量為0.36,其中2個(gè)微衛(wèi)星標(biāo)記為高度多態(tài)標(biāo)記,9個(gè)微衛(wèi)星標(biāo)記為中度多態(tài)標(biāo)記,2個(gè)微衛(wèi)星標(biāo)記為低度多態(tài)標(biāo)記,因此本研究使用的13個(gè)微衛(wèi)星標(biāo)記可作為研究浙江長(zhǎng)興尹家邊揚(yáng)子鱷保護(hù)區(qū)揚(yáng)子鱷種群遺傳多樣性的有效遺傳標(biāo)記[47]。浙江長(zhǎng)興尹家邊揚(yáng)子鱷保護(hù)區(qū)野放區(qū)揚(yáng)子鱷種群的平均有效等位基因數(shù)量、平均觀測(cè)雜合度、平均期望雜合度、平均多態(tài)信息含量和平均香農(nóng)信息指數(shù)均低于保護(hù)區(qū)核心區(qū)揚(yáng)子鱷種群,表明保護(hù)區(qū)野放區(qū)揚(yáng)子鱷種群的遺傳多樣性相對(duì)缺乏。這一結(jié)果可能與在野放過(guò)程中未對(duì)核心區(qū)揚(yáng)子鱷種群、野放區(qū)揚(yáng)子鱷種群以及挑選的野放個(gè)體的遺傳多樣性進(jìn)行系統(tǒng)性分析有關(guān)。

    本研究利用13個(gè)微衛(wèi)星標(biāo)記對(duì)浙江長(zhǎng)興尹家邊揚(yáng)子鱷保護(hù)區(qū)揚(yáng)子鱷種群進(jìn)行遺傳多樣性分析,并比較保護(hù)區(qū)核心區(qū)和野放區(qū)揚(yáng)子鱷種群的遺傳多樣性差異,了解了浙江長(zhǎng)興尹家邊揚(yáng)子鱷保護(hù)區(qū)揚(yáng)子鱷種群在生存和繁衍過(guò)程中可能存在低頻等位基因丟失的情況。本研究?jī)H檢測(cè)到2個(gè)稀有等位基因,因此在后續(xù)的揚(yáng)子鱷種群保護(hù)繁殖過(guò)程中,保護(hù)區(qū)應(yīng)注重對(duì)低頻等位基因及稀有等位基因的監(jiān)測(cè)和保護(hù)。此外,本研究還發(fā)現(xiàn),雖然2012—2019年保護(hù)區(qū)核心區(qū)揚(yáng)子鱷種群中已有960條揚(yáng)子鱷被放歸至野放區(qū),但保護(hù)區(qū)野放區(qū)揚(yáng)子鱷種群的遺傳多樣性與核心區(qū)揚(yáng)子鱷種群仍存在一定差異[36-37],野放區(qū)揚(yáng)子鱷種群的平均有效等位基因數(shù)量、平均觀測(cè)雜合度、平均期望雜合度、平均多態(tài)信息含量和平均香農(nóng)信息指數(shù)均低于核心區(qū)揚(yáng)子鱷種群,因此在后續(xù)的揚(yáng)子鱷野外放歸過(guò)程中,保護(hù)區(qū)應(yīng)注重分析核心區(qū)和野放區(qū)揚(yáng)子鱷種群的遺傳多樣性,保證保護(hù)區(qū)揚(yáng)子鱷種群的遺傳多樣性不會(huì)在野放過(guò)程中降低,從核心區(qū)揚(yáng)子鱷種群中挑選出合適的野放個(gè)體進(jìn)行野外放歸,從而豐富野放區(qū)揚(yáng)子鱷種群的遺傳多樣性。

    猜你喜歡
    揚(yáng)子鱷多態(tài)微衛(wèi)星
    揚(yáng)子鱷分四步過(guò)冬
    梅花鹿基因組微衛(wèi)星分布特征研究
    分層多態(tài)加權(quán)k/n系統(tǒng)的可用性建模與設(shè)計(jì)優(yōu)化
    成長(zhǎng)吧!揚(yáng)子鱷
    參差多態(tài)而功不唐捐
    浙江進(jìn)行揚(yáng)子鱷第7次野外放歸
    林麝全基因組微衛(wèi)星分布規(guī)律研究
    人多巴胺D2基因啟動(dòng)子區(qū)—350A/G多態(tài)位點(diǎn)熒光素酶表達(dá)載體的構(gòu)建與鑒定及活性檢測(cè)
    基于轉(zhuǎn)錄組測(cè)序的波紋巴非蛤微衛(wèi)星標(biāo)記研究
    揚(yáng)子鱷回家
    国产淫语在线视频| 小蜜桃在线观看免费完整版高清| 夜夜看夜夜爽夜夜摸| 国产精品无大码| 午夜激情久久久久久久| 99久国产av精品国产电影| 老熟女久久久| 国产一区二区三区综合在线观看 | 啦啦啦啦在线视频资源| 蜜桃在线观看..| 亚洲精品aⅴ在线观看| 搡女人真爽免费视频火全软件| 国产综合精华液| a 毛片基地| 国产亚洲av片在线观看秒播厂| 亚洲,欧美,日韩| 国产高潮美女av| 黑人高潮一二区| 97热精品久久久久久| 亚洲av不卡在线观看| 久久av网站| 国产乱人视频| 国产精品免费大片| 黄色视频在线播放观看不卡| 成人一区二区视频在线观看| 男人舔奶头视频| 汤姆久久久久久久影院中文字幕| 精品久久国产蜜桃| 熟女人妻精品中文字幕| 大又大粗又爽又黄少妇毛片口| 国产精品久久久久久久电影| 又粗又硬又长又爽又黄的视频| 婷婷色综合www| 久久精品国产鲁丝片午夜精品| 女人久久www免费人成看片| 少妇人妻一区二区三区视频| 亚洲aⅴ乱码一区二区在线播放| 免费看日本二区| 国产探花极品一区二区| 成人国产麻豆网| 免费观看无遮挡的男女| 高清在线视频一区二区三区| 人人妻人人澡人人爽人人夜夜| 婷婷色综合大香蕉| 国产免费一区二区三区四区乱码| 乱码一卡2卡4卡精品| 国产高清不卡午夜福利| 大码成人一级视频| 黄色视频在线播放观看不卡| 亚洲无线观看免费| 日韩伦理黄色片| 国产黄片视频在线免费观看| 老熟女久久久| 身体一侧抽搐| av一本久久久久| 狂野欧美激情性bbbbbb| 亚洲无线观看免费| av播播在线观看一区| 99九九线精品视频在线观看视频| 人妻系列 视频| 国产精品人妻久久久久久| 国产乱人视频| 国产高清国产精品国产三级 | 婷婷色综合大香蕉| 亚洲综合精品二区| 日韩欧美一区视频在线观看 | 一级a做视频免费观看| 国产老妇伦熟女老妇高清| 日韩成人av中文字幕在线观看| 18禁在线无遮挡免费观看视频| 国产在线一区二区三区精| 亚洲av二区三区四区| 亚洲怡红院男人天堂| 亚洲成人手机| 久久国产精品男人的天堂亚洲 | 极品教师在线视频| 色网站视频免费| 亚洲av男天堂| 亚洲va在线va天堂va国产| 在线观看人妻少妇| 插逼视频在线观看| 777米奇影视久久| 成年人午夜在线观看视频| 久久久精品免费免费高清| 精品人妻偷拍中文字幕| 国产黄片视频在线免费观看| 国产久久久一区二区三区| 大陆偷拍与自拍| 亚洲综合色惰| 高清欧美精品videossex| 又黄又爽又刺激的免费视频.| 欧美xxⅹ黑人| 多毛熟女@视频| 日韩一区二区三区影片| 女人久久www免费人成看片| 精品久久久久久电影网| 人妻少妇偷人精品九色| 精品少妇黑人巨大在线播放| 日日撸夜夜添| 最近的中文字幕免费完整| 80岁老熟妇乱子伦牲交| av国产免费在线观看| 毛片一级片免费看久久久久| 亚洲国产精品国产精品| 国产高潮美女av| 国产高清国产精品国产三级 | 国产伦精品一区二区三区视频9| 王馨瑶露胸无遮挡在线观看| 秋霞伦理黄片| 免费播放大片免费观看视频在线观看| 国产精品熟女久久久久浪| 一级毛片aaaaaa免费看小| 亚洲激情五月婷婷啪啪| 色视频在线一区二区三区| 99热这里只有精品一区| 国产免费福利视频在线观看| 精品少妇久久久久久888优播| 国产爽快片一区二区三区| 国产熟女欧美一区二区| 男女国产视频网站| 99热这里只有精品一区| 国产爱豆传媒在线观看| 观看美女的网站| 欧美极品一区二区三区四区| 夜夜看夜夜爽夜夜摸| 国产精品免费大片| 欧美极品一区二区三区四区| 男的添女的下面高潮视频| 91久久精品国产一区二区三区| 热re99久久精品国产66热6| 久久久久国产网址| xxx大片免费视频| 99热全是精品| 国产白丝娇喘喷水9色精品| 最近最新中文字幕免费大全7| 国产老妇伦熟女老妇高清| 亚洲欧美成人精品一区二区| 99热国产这里只有精品6| 成人无遮挡网站| 国产色婷婷99| 五月伊人婷婷丁香| 麻豆成人午夜福利视频| 色综合色国产| 国产亚洲欧美精品永久| 成年av动漫网址| 欧美三级亚洲精品| 久久99热这里只频精品6学生| 国产久久久一区二区三区| 寂寞人妻少妇视频99o| 丰满乱子伦码专区| 伊人久久国产一区二区| 一区二区三区精品91| 免费大片18禁| 在线观看人妻少妇| 国产精品成人在线| 99久国产av精品国产电影| 免费看不卡的av| 欧美日韩精品成人综合77777| 国语对白做爰xxxⅹ性视频网站| 久久久久久久久久久丰满| 亚洲欧美精品自产自拍| 久久6这里有精品| 最近最新中文字幕大全电影3| 久久久久精品久久久久真实原创| 18禁裸乳无遮挡免费网站照片| 国产淫语在线视频| 亚洲美女视频黄频| a 毛片基地| 久久鲁丝午夜福利片| 高清日韩中文字幕在线| 三级国产精品片| 久久女婷五月综合色啪小说| 搡老乐熟女国产| 99热这里只有精品一区| 午夜福利在线在线| 日本欧美国产在线视频| 日日摸夜夜添夜夜添av毛片| 亚洲综合色惰| 亚洲精品456在线播放app| 人体艺术视频欧美日本| 国国产精品蜜臀av免费| 九九爱精品视频在线观看| 交换朋友夫妻互换小说| 黄色怎么调成土黄色| 中文在线观看免费www的网站| 中文字幕av成人在线电影| 久久青草综合色| 国产 一区 欧美 日韩| 国内少妇人妻偷人精品xxx网站| 久久精品久久精品一区二区三区| 午夜免费观看性视频| 亚洲成人中文字幕在线播放| 国产精品嫩草影院av在线观看| 99re6热这里在线精品视频| 亚洲av男天堂| 中文字幕亚洲精品专区| 国产极品天堂在线| 日韩亚洲欧美综合| 黄色日韩在线| 成人影院久久| 精品人妻一区二区三区麻豆| 欧美日韩在线观看h| 国产在视频线精品| 亚洲最大成人中文| 男女国产视频网站| 黄色欧美视频在线观看| 精华霜和精华液先用哪个| av不卡在线播放| 欧美+日韩+精品| 丰满人妻一区二区三区视频av| 久久人人爽人人片av| 国产精品一区二区在线观看99| 午夜激情久久久久久久| 亚洲图色成人| 自拍偷自拍亚洲精品老妇| 最后的刺客免费高清国语| 国产又色又爽无遮挡免| 亚洲自偷自拍三级| 亚洲国产成人一精品久久久| 免费高清在线观看视频在线观看| 久久久久精品性色| 亚洲第一区二区三区不卡| 视频中文字幕在线观看| 亚洲三级黄色毛片| 99热网站在线观看| 国产v大片淫在线免费观看| 永久免费av网站大全| 久久 成人 亚洲| 国产av一区二区精品久久 | 久久人人爽人人片av| 99精国产麻豆久久婷婷| av国产精品久久久久影院| 日韩视频在线欧美| 亚洲国产精品专区欧美| 中文字幕免费在线视频6| 国产在线视频一区二区| 少妇被粗大猛烈的视频| 亚洲国产欧美在线一区| 精华霜和精华液先用哪个| 国国产精品蜜臀av免费| 韩国高清视频一区二区三区| 女人久久www免费人成看片| 国产成人精品婷婷| 国产精品一区www在线观看| a级一级毛片免费在线观看| 全区人妻精品视频| 色哟哟·www| 超碰97精品在线观看| 国产伦理片在线播放av一区| 深夜a级毛片| 嫩草影院新地址| 蜜桃亚洲精品一区二区三区| 伦精品一区二区三区| 国产午夜精品久久久久久一区二区三区| av在线播放精品| 成年av动漫网址| 肉色欧美久久久久久久蜜桃| 日韩av免费高清视频| 一级黄片播放器| 国产一区二区三区综合在线观看 | 国产成人91sexporn| 最近手机中文字幕大全| 成年女人在线观看亚洲视频| 免费高清在线观看视频在线观看| 女人十人毛片免费观看3o分钟| 97在线人人人人妻| 大陆偷拍与自拍| 中文精品一卡2卡3卡4更新| 午夜精品国产一区二区电影| 欧美亚洲 丝袜 人妻 在线| 婷婷色综合www| 国产精品.久久久| 国产亚洲最大av| 国产一区有黄有色的免费视频| 成人国产麻豆网| 热re99久久精品国产66热6| 亚洲av欧美aⅴ国产| 性色avwww在线观看| 亚洲久久久国产精品| 欧美少妇被猛烈插入视频| 中文在线观看免费www的网站| av在线观看视频网站免费| 超碰av人人做人人爽久久| 老师上课跳d突然被开到最大视频| 男女边摸边吃奶| 国产综合精华液| 精品视频人人做人人爽| 黑人高潮一二区| av线在线观看网站| 最新中文字幕久久久久| 如何舔出高潮| 亚洲熟女精品中文字幕| 国产乱来视频区| 男的添女的下面高潮视频| 中文字幕精品免费在线观看视频 | 国产人妻一区二区三区在| 久久久久久久久久人人人人人人| 国产综合精华液| 99久久精品一区二区三区| 一级毛片 在线播放| 特大巨黑吊av在线直播| 黄片无遮挡物在线观看| 少妇熟女欧美另类| 男的添女的下面高潮视频| 亚洲va在线va天堂va国产| 国产精品一区二区三区四区免费观看| 国产av国产精品国产| 亚洲丝袜综合中文字幕| 中文精品一卡2卡3卡4更新| 涩涩av久久男人的天堂| 十分钟在线观看高清视频www | 久久婷婷青草| 一级a做视频免费观看| 国产日韩欧美在线精品| 国产男人的电影天堂91| 美女中出高潮动态图| 爱豆传媒免费全集在线观看| 国产 一区 欧美 日韩| 亚洲怡红院男人天堂| 亚洲aⅴ乱码一区二区在线播放| 午夜视频国产福利| 国国产精品蜜臀av免费| 国产伦精品一区二区三区四那| 亚洲av福利一区| 精品亚洲成a人片在线观看 | 老熟女久久久| 欧美xxⅹ黑人| 国产精品人妻久久久久久| 亚洲一区二区三区欧美精品| 欧美另类一区| 国产黄色视频一区二区在线观看| 熟女av电影| 观看av在线不卡| 久久国产亚洲av麻豆专区| kizo精华| 男女国产视频网站| 在线播放无遮挡| 一级毛片电影观看| 久久毛片免费看一区二区三区| 亚洲欧洲日产国产| 久热这里只有精品99| 久久人人爽人人爽人人片va| 美女中出高潮动态图| 自拍欧美九色日韩亚洲蝌蚪91 | 精品少妇黑人巨大在线播放| 久久人人爽av亚洲精品天堂 | 国产综合精华液| 18+在线观看网站| 国产精品三级大全| 成人午夜精彩视频在线观看| 久久久精品94久久精品| 熟女av电影| 国产精品精品国产色婷婷| 亚洲欧美成人精品一区二区| 视频中文字幕在线观看| 国产成人精品福利久久| 91精品伊人久久大香线蕉| 久久久久国产网址| 美女内射精品一级片tv| 国产有黄有色有爽视频| 亚洲精品日本国产第一区| 高清黄色对白视频在线免费看 | 丝袜喷水一区| xxx大片免费视频| 久久国产精品男人的天堂亚洲 | 伦精品一区二区三区| 激情 狠狠 欧美| 国产伦理片在线播放av一区| 国产高清有码在线观看视频| 亚洲国产日韩一区二区| 久久韩国三级中文字幕| 高清午夜精品一区二区三区| 亚洲四区av| 我的老师免费观看完整版| 99久国产av精品国产电影| 大话2 男鬼变身卡| 久久久久国产精品人妻一区二区| 亚洲欧美精品专区久久| 亚洲精品456在线播放app| 婷婷色综合大香蕉| 尤物成人国产欧美一区二区三区| 亚洲一级一片aⅴ在线观看| 在线观看人妻少妇| 国产男女内射视频| 中国美白少妇内射xxxbb| 国产综合精华液| 国产在线视频一区二区| 18禁在线无遮挡免费观看视频| 最近中文字幕2019免费版| 亚洲国产最新在线播放| 三级国产精品欧美在线观看| 熟女人妻精品中文字幕| 久久99蜜桃精品久久| 干丝袜人妻中文字幕| 国产一区二区三区av在线| 久久久久久久大尺度免费视频| 婷婷色综合大香蕉| av在线app专区| 最近最新中文字幕免费大全7| 亚洲第一av免费看| 夫妻午夜视频| 国产成人一区二区在线| 亚洲欧美日韩卡通动漫| av天堂中文字幕网| av不卡在线播放| 国内少妇人妻偷人精品xxx网站| 一级黄片播放器| 国产午夜精品久久久久久一区二区三区| 蜜桃亚洲精品一区二区三区| 少妇人妻久久综合中文| 精品一区二区三卡| 亚洲精品国产色婷婷电影| 伦精品一区二区三区| 高清不卡的av网站| 最近中文字幕2019免费版| 日本欧美国产在线视频| 亚洲av中文av极速乱| 久久精品久久久久久久性| 观看av在线不卡| 久久精品国产亚洲av天美| 一本—道久久a久久精品蜜桃钙片| 在线观看一区二区三区激情| 多毛熟女@视频| 久久av网站| 亚洲aⅴ乱码一区二区在线播放| av福利片在线观看| 一本一本综合久久| 亚洲精品视频女| 欧美3d第一页| 亚洲欧美成人综合另类久久久| 亚洲精品一区蜜桃| 国产男女内射视频| 亚洲久久久国产精品| 男男h啪啪无遮挡| 国产熟女欧美一区二区| 狂野欧美白嫩少妇大欣赏| 只有这里有精品99| 丝瓜视频免费看黄片| 亚州av有码| 欧美变态另类bdsm刘玥| 国产av国产精品国产| 精品国产露脸久久av麻豆| 一本一本综合久久| 国产黄频视频在线观看| av线在线观看网站| 日本猛色少妇xxxxx猛交久久| 免费看不卡的av| 久久97久久精品| 九九在线视频观看精品| 乱系列少妇在线播放| 街头女战士在线观看网站| 中文字幕亚洲精品专区| 国产精品人妻久久久影院| 久久久精品免费免费高清| 校园人妻丝袜中文字幕| 国产精品久久久久久久电影| 少妇熟女欧美另类| 一区二区av电影网| 成人一区二区视频在线观看| 校园人妻丝袜中文字幕| 日本免费在线观看一区| 免费在线观看成人毛片| 少妇丰满av| 国产色爽女视频免费观看| 麻豆国产97在线/欧美| 国产伦理片在线播放av一区| 欧美成人一区二区免费高清观看| 亚洲第一区二区三区不卡| 亚洲最大成人中文| 一区二区三区乱码不卡18| 三级国产精品欧美在线观看| 干丝袜人妻中文字幕| www.色视频.com| 秋霞伦理黄片| 久久6这里有精品| 精品熟女少妇av免费看| 性色av一级| 草草在线视频免费看| 亚洲第一av免费看| 色视频在线一区二区三区| 99精国产麻豆久久婷婷| 久久人人爽人人片av| 成人国产麻豆网| 大又大粗又爽又黄少妇毛片口| 久久人人爽av亚洲精品天堂 | 精品人妻熟女av久视频| 黄片wwwwww| 久久国产精品男人的天堂亚洲 | 免费观看av网站的网址| 国产视频内射| 亚洲欧美日韩无卡精品| 色视频www国产| 深夜a级毛片| 91精品伊人久久大香线蕉| 亚洲图色成人| 少妇人妻 视频| 成年免费大片在线观看| 久久久久国产网址| 久久6这里有精品| 国产高清不卡午夜福利| 在线观看免费高清a一片| 黄色一级大片看看| 大片免费播放器 马上看| 精品酒店卫生间| 久久 成人 亚洲| 亚洲经典国产精华液单| 国产精品久久久久久av不卡| 波野结衣二区三区在线| 99久久精品国产国产毛片| 最近中文字幕2019免费版| 国产免费一区二区三区四区乱码| 国产男女超爽视频在线观看| 日日摸夜夜添夜夜爱| 人妻制服诱惑在线中文字幕| 成人18禁高潮啪啪吃奶动态图 | av国产免费在线观看| 极品教师在线视频| 少妇被粗大猛烈的视频| 亚洲不卡免费看| 少妇人妻久久综合中文| 乱码一卡2卡4卡精品| 亚洲精品色激情综合| 久久国产精品男人的天堂亚洲 | 人妻少妇偷人精品九色| 在线观看av片永久免费下载| 久久精品夜色国产| 国产精品偷伦视频观看了| 蜜桃亚洲精品一区二区三区| 精品久久久久久久久亚洲| 亚洲美女搞黄在线观看| 久久久久久久久大av| 老师上课跳d突然被开到最大视频| 久久久久视频综合| 欧美97在线视频| 色哟哟·www| 91狼人影院| 国产免费视频播放在线视频| 亚洲成人av在线免费| 最近2019中文字幕mv第一页| www.色视频.com| 日韩成人伦理影院| 一区二区三区乱码不卡18| 国产精品偷伦视频观看了| 国产精品麻豆人妻色哟哟久久| 观看美女的网站| 我的老师免费观看完整版| 18禁动态无遮挡网站| 国产成人精品婷婷| 亚洲色图av天堂| 亚洲精品第二区| 精品久久久久久电影网| 尾随美女入室| 色婷婷久久久亚洲欧美| 久久久久精品久久久久真实原创| 少妇的逼水好多| 国产免费一区二区三区四区乱码| 亚洲精品,欧美精品| 精品亚洲成国产av| 免费大片18禁| 午夜福利网站1000一区二区三区| 在线亚洲精品国产二区图片欧美 | 97超视频在线观看视频| 免费看光身美女| 黄色视频在线播放观看不卡| 少妇人妻久久综合中文| 亚洲,欧美,日韩| 亚洲国产色片| 亚洲伊人久久精品综合| 中国美白少妇内射xxxbb| 亚洲av免费高清在线观看| 3wmmmm亚洲av在线观看| 午夜免费观看性视频| av国产久精品久网站免费入址| 久久99蜜桃精品久久| 欧美97在线视频| 亚洲成人手机| 中文欧美无线码| 欧美一级a爱片免费观看看| 久久国产亚洲av麻豆专区| freevideosex欧美| 精品亚洲成国产av| 久久97久久精品| 亚洲精品日韩在线中文字幕| 亚洲国产毛片av蜜桃av| 91久久精品国产一区二区成人| 久久久午夜欧美精品| 国产成人精品一,二区| 亚洲色图综合在线观看| 婷婷色麻豆天堂久久| 亚洲人成网站高清观看| 少妇猛男粗大的猛烈进出视频| 汤姆久久久久久久影院中文字幕| 国产日韩欧美在线精品| 毛片一级片免费看久久久久| 久久99蜜桃精品久久| 亚洲丝袜综合中文字幕| 最后的刺客免费高清国语| av国产免费在线观看| 国产精品国产av在线观看| 国产亚洲精品久久久com| 免费播放大片免费观看视频在线观看| 91精品一卡2卡3卡4卡| 亚洲精品日本国产第一区| 久久99热这里只有精品18| 91狼人影院| 熟女人妻精品中文字幕| 丰满人妻一区二区三区视频av| 国产在线免费精品| 三级国产精品片| 九色成人免费人妻av| av天堂中文字幕网| 久热久热在线精品观看| 欧美另类一区| 制服丝袜香蕉在线| 精品亚洲乱码少妇综合久久|