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

    新疆脹果甘草遺傳多樣性的ISSR分析

    2021-11-21 11:15:06王莉莉阿力木江排爾哈提李莉偉居來(lái)提托合提
    中草藥 2021年22期

    王莉莉,阿力木江·排爾哈提,李莉偉,居來(lái)提·托合提,劉 忠

    新疆脹果甘草遺傳多樣性的ISSR分析

    王莉莉1,阿力木江·排爾哈提2#,李莉偉1,居來(lái)提·托合提2,劉 忠*

    1.上海交通大學(xué)藥學(xué)院,上海 200240 2.新疆維吾爾醫(yī)學(xué)??茖W(xué)校,新疆 和田 848000

    揭示新疆塔里木盆地分布的野生脹果甘草的遺傳多樣性和種群遺傳結(jié)構(gòu)。利用ISSR分子標(biāo)記方法,對(duì)來(lái)自新疆南疆不同地理方位的脹果甘草6個(gè)種群總計(jì)167個(gè)個(gè)體進(jìn)行分析。對(duì)脹果甘草6個(gè)種群共篩選出46條引物,擴(kuò)增獲得193條多態(tài)條帶,多態(tài)位點(diǎn)比率(PPL)為69.68%,Nei’s遺傳多樣性指數(shù)()為0.272 2,Shannon多樣性指數(shù)()為0.401 0;Nei’s總種群遺傳多樣性指數(shù)(t)為0.513 0,種群間遺傳分化系數(shù)(st)為0.530 7,基因流(m)為0.438 8,Mantel檢驗(yàn)顯示種群間的遺傳分化與地理距離無(wú)相關(guān)性;種群間遺傳變異豐度和遺傳多樣性水平呈現(xiàn)BC>KEL>PS>QM>AQS>KC的遞減趨勢(shì),聚類分析將6個(gè)種群在遺傳相似性系數(shù)為0.69處分為3類;6個(gè)種群內(nèi)個(gè)體間PPL的變化范圍為52.11%~81.87%,觀測(cè)等位基因數(shù)(a)和有效等位基因數(shù)(e)的變化范圍分別為1.521 1~1.818 7和1.307 9~1.498 9,平均為0.240 8,平均為0.357 1。脹果甘草有較高的遺傳多樣性水平,但種群間的分化程度較種群內(nèi)大,遺傳多樣性更多的存在于種群水平。生境片段化可能是造成其遺傳多樣性空間分布現(xiàn)狀的主要原因,土壤水分含量可能是影響其種群遺傳多樣性的制約因素。研究結(jié)果對(duì)全面了解和掌握脹果甘草自然資源現(xiàn)狀、物種進(jìn)化潛力等,制定資源利用與保護(hù)的正確策略提供了重要的研究基礎(chǔ)。

    脹果甘草;ISSR標(biāo)記;遺傳多樣性;遺傳結(jié)構(gòu);資源保護(hù)

    甘草et是常用中藥材,被譽(yù)為百藥之首。甘草也是優(yōu)良的甜味劑、乳化劑和矯味劑,廣泛應(yīng)用于食品、飲料、煙草、日化等工業(yè),市場(chǎng)需求量巨大。正品甘草的基原植物有3種,即豆科甘草屬L.植物烏拉爾甘草Fisch.、脹果甘草Bat.和光果甘草L.[1]。脹果甘草主要分布在新疆南疆,少部分經(jīng)東疆?dāng)U展至甘肅疏勒河和額濟(jì)納河流域呈零星分布[2-3]。脹果甘草自然分布區(qū)域的生態(tài)環(huán)境均極為嚴(yán)酷,因而也極為脆弱,一旦遭到破壞,則難以恢復(fù),野生資源總量十分有限,在商品甘草中的占比遠(yuǎn)小于烏拉爾甘草。但是,脹果甘草具有不同于烏拉爾甘草的品質(zhì)和開發(fā)利用價(jià)值。脹果甘草不僅在總黃酮和多糖含量上均高于烏拉爾甘草,而且在黃酮化合物的組成上也不同于烏拉爾甘草,顯示出與烏拉爾甘草有所不同的藥理活性[4-5]。在生態(tài)適應(yīng)能力上,脹果甘草也遠(yuǎn)勝烏拉爾甘草,具有極強(qiáng)的抗旱、耐鹽堿、耐沙埋等抗極端不利環(huán)境的能力[6]。

    遺傳多樣性是生物多樣性的核心內(nèi)容[7-8]。對(duì)物種遺傳多樣性的了解是揭示其資源水平、生境適應(yīng)性、進(jìn)化潛力的基礎(chǔ),也是合理開發(fā)利用與保護(hù)的依據(jù)[9]。藥用植物遺傳多樣性分析可為新型、特異新成分的發(fā)掘提供有益的參考和指導(dǎo),在天然藥物的研發(fā)中具有積極意義。簡(jiǎn)單重復(fù)序列間片段多態(tài)性(inter-simple sequence repeat,ISSR)分子標(biāo)記方法是在簡(jiǎn)單重復(fù)序列(simple sequence repeat,SSR)分子標(biāo)記方法基礎(chǔ)上發(fā)展起來(lái)的檢測(cè)DNA片段多態(tài)性的技術(shù)。ISSR標(biāo)記具有擴(kuò)增穩(wěn)定性好、檢測(cè)能力強(qiáng)、分辨力高等優(yōu)點(diǎn),被廣泛應(yīng)用于遺傳多樣性評(píng)價(jià)、優(yōu)良種質(zhì)發(fā)掘、基因定位、遺傳作圖、輔助育種等方面的研究[10-12]。

    脹果甘草的基礎(chǔ)研究十分欠缺,對(duì)其資源蘊(yùn)藏量、遺傳多樣性水平、品質(zhì)變異格局等資源狀況都所知甚少[13]。本研究采用ISSR分子標(biāo)記技術(shù),對(duì)新疆南疆不同地理來(lái)源的脹果甘草野生種群進(jìn)行遺傳多樣性分析,評(píng)估脹果甘草遺傳多樣性水平,為脹果甘草生態(tài)適應(yīng)機(jī)制、品質(zhì)變異特征、優(yōu)良種質(zhì)篩選、資源保護(hù)的研究提供理論指導(dǎo),同時(shí)也為脹果甘草規(guī)范化種植優(yōu)良栽培品系的培育奠定基礎(chǔ),是脹果甘草開發(fā)利用中不可或缺的基礎(chǔ)性研究工作。

    1 材料和儀器

    1.1 材料

    材料來(lái)自于新疆南疆6個(gè)不同縣市自然分布的脹果甘草野生種群(表1)。經(jīng)上海交通大學(xué)藥學(xué)院劉忠副教授鑒定為豆科甘草屬植物脹果甘草Bat.,憑證標(biāo)本保存于新疆維吾爾醫(yī)學(xué)??茖W(xué)校藥用植物標(biāo)本室(新疆和田市)。每個(gè)種群采集20個(gè)以上個(gè)體,6個(gè)種群共167個(gè)個(gè)體。采集時(shí),選取新鮮、健康、幼嫩的葉片放入裝有變色硅膠顆粒的自封袋中,標(biāo)明采集信息。及時(shí)更換變色硅膠,使樣品快速脫水,直至完全干燥。

    1.2 儀器

    JXFSPRP-24型組織研磨儀(凈信實(shí)業(yè)有限公司,上海),Mastercycler nexus GX2 PCR儀(Eppendorf,德國(guó)),1600 凝膠成像分析系統(tǒng)(天能有限公司,上海),NanoDrop ONE(Thermo,美國(guó)),DYY-6C型電泳儀(六一儀器廠,北京)。

    表1 脹果甘草的樣品來(lái)源信息

    2 方法

    2.1 脹果甘草基因組DNA提取與質(zhì)量檢測(cè)

    在進(jìn)行種群遺傳多樣性分析時(shí),采取混合取樣法,即每個(gè)種群均稱取其所有個(gè)體的等量葉片,充分混合后用于總DNA的提取,將種群整體上作為一個(gè)分析樣品。在進(jìn)行種群遺傳結(jié)構(gòu)檢測(cè)時(shí),采取個(gè)體取樣法,即每一個(gè)種群均對(duì)其內(nèi)每一個(gè)個(gè)體單獨(dú)提取基因組DNA,用于ISSR-PCR擴(kuò)增,每個(gè)個(gè)體均為一個(gè)樣品。

    利用CTAB法進(jìn)行脹果甘草基因組DNA的提取。參照相關(guān)文獻(xiàn)并加以適當(dāng)改進(jìn)[14]。主要步驟如下:稱取0.1 g研磨成粉末狀的植物材料,加入600 μL 3% CTAB提取液,65 ℃水浴中溫育30 min;加入600 μL氯仿-異戊醇(24∶1,現(xiàn)配),靜置,離心,取上清;加入預(yù)冷無(wú)水乙醇,靜置,離心,棄上清;室溫風(fēng)干,除去殘留乙醇,加入滅菌蒸餾水溶解。以1%瓊脂糖凝膠電泳檢測(cè)所提DNA質(zhì)量。Nanodrop ONE(Thermo公司,美國(guó))測(cè)定樣品DNA的濃度。

    2.2 引物的來(lái)源

    基于哥倫比亞大學(xué)(UBC)公布的100條ISSR通用引物序列信息,以及豆科植物相關(guān)文獻(xiàn)報(bào)道的引物序列信息[15-18],委托生工生物(上海)公司進(jìn)行引物合成。

    2.3 ISSR分析檢測(cè)體系的優(yōu)化

    ISSR-PCR反應(yīng)體系為:2×Taq PCR Master Mix(12.5 μL)、引物10 pmol/μL、DNA模板2 ng/μL,共25 μL。擴(kuò)增程序:94 ℃預(yù)變性5 min;94 ℃變性45 s,53 ℃退火60 s,70 ℃延伸90 s,30個(gè)循環(huán);72 ℃保持7 min。本實(shí)驗(yàn)對(duì)該擴(kuò)增程序的引物退火溫度和擴(kuò)增循環(huán)次數(shù)進(jìn)行了優(yōu)化。

    預(yù)實(shí)驗(yàn)提示,退火溫度不僅影響ISSR-PCR擴(kuò)增條帶的清晰度,而且可造成條帶數(shù)目的變化。實(shí)驗(yàn)中,以每條引物m值的±5 ℃為上下限設(shè)置溫度梯度,進(jìn)行梯度PCR擴(kuò)增,篩選引物最適退火溫度;設(shè)置30、32、35、38、40 5種不同擴(kuò)增循環(huán)次數(shù),考察擴(kuò)增循環(huán)次數(shù)的優(yōu)化設(shè)置。以巴音郭楞蒙古自治州且末縣種群(QM)和庫(kù)爾勒市種群(KEL)為樣本,采用優(yōu)化后的擴(kuò)增程序進(jìn)行擴(kuò)增,驗(yàn)證上述優(yōu)化的ISSR-PCR擴(kuò)增程序的實(shí)效性,并進(jìn)行有效引物篩選。

    2.4 新疆脹果甘草遺傳多樣性ISSR檢測(cè)

    以篩選得到的條帶明晰、多態(tài)性豐富,并且穩(wěn)定性良好的引物,分別對(duì)6個(gè)種群進(jìn)行種群間遺傳多樣性和種群遺傳結(jié)構(gòu)的ISSR分析。

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

    對(duì)擴(kuò)增產(chǎn)物進(jìn)行人工判讀,有條帶記錄為“1”,無(wú)條帶則記錄為“0”,形成0、1排列的原始數(shù)據(jù)矩陣。利用Popgene1.32軟件計(jì)算脹果甘草遺傳多樣性參數(shù)。利用NTSYSpc 2.10e軟件構(gòu)建種群遺傳關(guān)系聚類圖。利用Mathematica軟件計(jì)算6個(gè)種群間的地理距離,Mantel test分析6個(gè)種群遺傳距離和地理距離的相關(guān)性。

    3 結(jié)果與分析

    3.1 ISSR反應(yīng)條件的優(yōu)化

    3.1.1 引物退火溫度的篩選 以引物UBC807、UBC808、UBC809、UBC810為例,對(duì)比不同退火溫度對(duì)ISSR-PCR擴(kuò)增結(jié)果的影響(圖1)。結(jié)果顯示,當(dāng)引物UBC807、UBC808、UBC809、UBC810的退火溫度分別為50.1、46.0、47.0、45.0 ℃時(shí),擴(kuò)增獲得的條帶不夠清晰、多態(tài)性不豐富,但當(dāng)每個(gè)引物對(duì)應(yīng)的退火溫度分別為46.1、52.0、51.0、49.0 ℃時(shí),擴(kuò)增產(chǎn)物不僅條帶清晰,而且條帶數(shù)目多、多態(tài)性豐富,擴(kuò)增效果良好。其他引物的最佳退火溫度見表2。

    3.1.2 循環(huán)次數(shù)的優(yōu)化 對(duì)脹果甘草的樣品進(jìn)行ISSR擴(kuò)增時(shí),當(dāng)PCR反應(yīng)的循環(huán)次數(shù)設(shè)置為30次時(shí),電泳檢測(cè)結(jié)果顯示沒(méi)有條帶,進(jìn)一步將循環(huán)次數(shù)設(shè)置為36次和38次,檢測(cè)結(jié)果無(wú)明顯差距,而循環(huán)次數(shù)為40次時(shí),得到了條帶清晰度高、多態(tài)性較豐富的擴(kuò)增結(jié)果,因此,40次擴(kuò)增循環(huán)次數(shù)是進(jìn)行脹果甘草葉片基因組DNA ISSR-PCR反應(yīng)的適宜條件。

    M-Marker 1~2-UBC807,50.1、46.1 ℃ 3~4-UBC808,46、52 ℃ 5~6-UBC809,47、51 ℃ 7~8-UBC810,49、45 ℃

    表2 各引物的退火溫度

    3.2 脹果甘草的遺傳多樣性檢測(cè)

    3.2.1 種群間的遺傳多樣性 在對(duì)脹果甘草6個(gè)種群間的遺傳多樣性分析中(圖2),篩選出46條引物,擴(kuò)增獲得的總條帶數(shù)為277條,其中,有193條多態(tài)條帶,多態(tài)條帶比率(PPL)為69.68%,平均每個(gè)引物可擴(kuò)增得到4.2個(gè)多態(tài)條帶。Nei’s遺傳多樣性指數(shù)()為0.272 2,Shannon多樣性指數(shù)()為0.401 0,顯示脹果甘草在基因水平上具有較豐富的遺傳多樣性。

    圖2 引物UBC811的ISSR-PCR的結(jié)果分析

    3.2.2 種群內(nèi)個(gè)體間的遺傳多樣性 對(duì)脹果甘草6個(gè)種群進(jìn)行遺傳結(jié)構(gòu)分析的結(jié)果顯示(圖3和表3),6個(gè)種群的PPL的變化范圍在52.11%~81.87%,觀測(cè)等位基因數(shù)(a)和有效等位基因數(shù)(e)的變化范圍分別為1.521 1~1.818 7和1.307 9~1.498 9,變化范圍為0.183 1~0.288 6(平均0.240 8),變化范圍為0.274 1~0.429 1(平均0.357 1)。其中,阿克蘇地區(qū)拜城縣種群(BC)遺傳變異豐富,遺傳多樣性水平高,而阿克蘇地區(qū)庫(kù)車市種群(KC)的變異相對(duì)貧乏,遺傳多樣性水平較低,群體間遺傳變異豐度和遺傳多樣性水平呈現(xiàn)BC>KEL>PS>QM>AQS>KC的遞減趨勢(shì)。

    圖3 引物UBC844對(duì)KC種群的ISSR擴(kuò)增結(jié)果

    表3 種群內(nèi)個(gè)體間的遺傳多樣性分析

    3.3 種群間遺傳分化與聚類分析

    Nei’s種群間遺傳多樣性指數(shù)(st)為0.272 2,Nei’s種群內(nèi)遺傳多樣性指數(shù)(s)為0.240 8,計(jì)算可知Nei’s總遺傳多樣性指數(shù)(t)為0.513 0,種群間遺傳分化系數(shù)(st)為0.530 7,即脹果甘草種群間的遺傳變異水平占總遺傳變異水平的53.07%,基因流(m)為0.438 8,小于1,顯示脹果甘草種群間的分化程度大于種群內(nèi)的分化程度。脹果甘草的遺傳多樣性更多的存在于種群水平。

    根據(jù)6個(gè)脹果甘草種群間的Nei’s遺傳距離和遺傳一致度(表4)利用NTSYSpc 2.10e軟件UPGMA方法對(duì)脹果甘草6個(gè)種群的聚類分析(圖4)表明,6個(gè)種群在遺傳相似性系數(shù)為0.69處分成3個(gè)大類。首先,阿克蘇地區(qū)拜城縣種群(BC)單獨(dú)成為一類,表現(xiàn)出與來(lái)自其他5個(gè)地理分布區(qū)的脹果甘草種群之間較少的遺傳相似性,說(shuō)明該種群與其他種群之間的遺傳分化程度較大;其次,阿克蘇地區(qū)庫(kù)車市種群(KC)與阿克蘇市種群(AQS)聚為第2類;第3,巴音郭楞蒙古州的庫(kù)爾勒市種群(KEL)與且末縣種群(QM)的遺傳距離最近,兩者先聚在一起,然后再進(jìn)一步與和田地區(qū)皮山縣種群(PS)聚合為第3類。Mantel檢驗(yàn)結(jié)果顯示群體間的遺傳距離與地理距離之間無(wú)相關(guān)性(=0.158 7,=0.254 7)。

    表4 不同地理方位的脹果甘草的遺傳距離和遺傳一致度

    左下角表示Nei’s遺傳距離,右上角表示遺傳一致度

    Figures in the lower left half part of the table indicates Nei ’s genetic distance, whereas those in the upper right half part indicates genetic consistency

    圖4 脹果甘草種群間的聚類分析

    4 討論

    4.1 伴隨生境片段化的種群隔離分化可能是造成脹果甘草遺傳多樣性空間分布現(xiàn)狀的主要原因

    遺傳多樣性反映了物種在起源演化和環(huán)境適應(yīng)上的特征[19-21]。遺傳多樣性豐富的物種具有更多的等位基因,對(duì)環(huán)境變化具有更好的適應(yīng)性和更強(qiáng)大的抵御能力[22]。影響物種遺傳多樣性和種群遺傳分化水平因素包括了地理隔離、生境片段化等外部因素,與繁育系統(tǒng)、遺傳漂變、基因突變和基因流等內(nèi)部因素[23-24]。

    脹果甘草在生境適應(yīng)性方面具有突出的特點(diǎn),脹果甘草是一種生長(zhǎng)在極端干旱區(qū)綠洲至沙漠過(guò)渡帶的荒漠植被類群,塔里木盆地是其主要分布區(qū)域。根據(jù)古地磁學(xué)研究結(jié)果,新疆南疆地區(qū)經(jīng)歷了多次海浸海退過(guò)程,最終整體成陸于三疊紀(jì)至侏羅紀(jì)時(shí)期,因此,當(dāng)被子植物在白堊紀(jì)開始繁盛和演化時(shí),這里應(yīng)該是氣候溫和濕潤(rùn),條件適宜植被發(fā)育的[25-26]。至第四紀(jì)青藏高原開始隆升,塔里木盆地逐漸形成,周緣的高原和高山阻擋了大氣環(huán)流及其攜帶的水汽,導(dǎo)致盆地日益變得干旱,土壤沙漠化程度和范圍不斷加深擴(kuò)大,期間雖然由于冰期與間冰期的相互交替使得沙漠與綠洲的范圍發(fā)生震蕩,但總的趨勢(shì)仍然是沙漠化不斷加劇[27-28]。目前,南疆綠洲在盆地周緣主要呈片狀間斷分布,部分區(qū)域片狀綠洲可沿河流相互聯(lián)結(jié)呈串珠狀[29]。遺傳多樣性分析結(jié)果表明新疆南疆脹果甘草具有較豐富的遺傳多樣性,且遺傳多樣性更多地體現(xiàn)在種群之間,同時(shí),種群間的遺傳分化與地理距離無(wú)相關(guān)性。結(jié)合地質(zhì)變遷、氣候變化和植物演化的研究資料,推測(cè)在塔克拉瑪干沙漠形成和不斷擴(kuò)大的過(guò)程中,生境片段化帶來(lái)的種群隔離分化是脹果甘草種群間遺傳分化和遺傳多樣性基本特征形成的主要原因。

    蒺藜科白刺屬植物大白刺Kom.在地理分布和生境上與脹果甘草十分近似。對(duì)大白刺遺傳多樣性的研究結(jié)果很好地支持了這一推論。大白刺主要分布于塔里木盆地、河西走廊和內(nèi)蒙古阿拉善。陸偉等[30]利用葉綠體基因組DNA片段的序列變異分析了大白刺的種群遺傳結(jié)構(gòu)和譜系地理特征,所得研究結(jié)果顯示,中國(guó)西北干旱區(qū)更新世以來(lái)持續(xù)沙漠化擴(kuò)張導(dǎo)致大白刺棲息地生境破碎化,種群被隔離,使種群間基因流受到限制,致使種群間表現(xiàn)出較高的遺傳多樣性。在對(duì)另一種塔里木盆地中的沙生植物沙生檉柳M.T.Liu所進(jìn)行的遺傳多樣性分析中[31],結(jié)果也表明,分布于盆地邊緣的一些沙生檉柳種群發(fā)生了明顯的遺傳分化,地理隔離可能是這些群體分化的主要因素。

    而在種群內(nèi),由于地理隔離帶來(lái)的種群間基因交流受阻,通過(guò)雜交增加遺傳多樣性的效應(yīng)變小,而在另一方面,脹果甘草較高的自交率[32-34]又使得脹果甘草個(gè)體間基因交流的頻率顯著降低。一些由于生境片段化導(dǎo)致的種群規(guī)模變小的種群還會(huì)顯著偏離哈迪-溫伯格平衡狀態(tài),進(jìn)一步加大遺傳漂變引起種群遺傳多樣性降低的效應(yīng)[35-36]。在諸多因素的綜合作用下,脹果甘草種群內(nèi)遺傳多樣性趨于丟失和貧乏化,使得脹果甘草的遺傳變異更多地存在于種群之間。

    4.2 土壤水分含量可能是影響脹果甘草維持或增加種群遺傳多樣性的制約因素

    在所檢測(cè)的脹果甘草6個(gè)野生種群中,阿克蘇地區(qū)拜城縣種群(BC)的遺傳多樣性水平最高,并且與其它5個(gè)種群的遺傳距離均很大,在聚類分析中單獨(dú)成為一類,而同為阿克蘇地區(qū)的庫(kù)車市種群(KC)卻顯示出最低的遺傳多樣性水平。目前尚沒(méi)有直接的依據(jù)來(lái)闡述這一結(jié)果,還有待進(jìn)一步的研究。

    土壤水分含量是土壤的一個(gè)關(guān)鍵性環(huán)境變量。土壤是植物生長(zhǎng)必須依賴的基本條件,土壤水分含量是影響干旱半干旱地區(qū)沙漠生態(tài)系統(tǒng)制約性環(huán)境因素,與生態(tài)系統(tǒng)的物質(zhì)循環(huán)、能量流動(dòng)、信息傳遞過(guò)程都有直接或間接的關(guān)系,對(duì)沙地土壤的發(fā)生、演化和土地生產(chǎn)力起著決定性的作用,在土地沙漠化區(qū)域作為維持生態(tài)系統(tǒng)穩(wěn)定、正常發(fā)揮生態(tài)系統(tǒng)結(jié)構(gòu)和功能的關(guān)鍵因子而存在,對(duì)整個(gè)生態(tài)系統(tǒng)的水熱平衡狀態(tài)的維護(hù)起決定性作用[37]。辛倩[38]在甘草屬植物中開展研究的結(jié)果表明土壤因素與遺傳多樣性趨勢(shì)存在一致性,脹果甘草分布區(qū)的土壤含水量越低,遺傳多樣性越低,干旱環(huán)境的選擇壓力使脹果甘草耐旱性強(qiáng)的基因型適應(yīng)極端環(huán)境形成,反過(guò)來(lái),土壤含水量越高,另一親本的出現(xiàn)有可能促進(jìn)雜交和漸滲種的產(chǎn)生,導(dǎo)致遺傳多樣性升高。宋雪梅[39]通過(guò)研究表明內(nèi)蒙古高原紅砂(Pall.) Maxim.種群間遺傳變異的形成是生境中水熱組合共同作用的結(jié)果,紅砂種群的遺傳差異在很大程度上是由其生長(zhǎng)環(huán)境中的土壤平均含水量決定的。張志南等[40]對(duì)黃土高原半干旱區(qū)天然草地的群落結(jié)構(gòu)研究發(fā)現(xiàn),物種多樣性和物種豐富度均隨著土壤含水量的增加而增加,土壤水分是限制半干旱地區(qū)群落生產(chǎn)力和物種多樣性的重要因子,而禾本科和豆科作為指示物種可以反映土壤的水分狀況,含水量越高,群體的多樣性越高。因此,在適宜的生境生長(zhǎng),脹果甘草的基因庫(kù)不至于遭受破壞,避免了遺傳漂變的發(fā)生,種群內(nèi)的遺傳多樣性也會(huì)相應(yīng)地升高。拜城縣的脹果甘草種群生長(zhǎng)在河灘地帶,而其他5個(gè)種群均生長(zhǎng)在村落旁,在土壤水分含量方面有著明顯不同。土壤水分含量可能是影響脹果甘草遺傳分化的核心因素。

    利益沖突 所有作者均聲明不存在利益沖突

    [1] 中國(guó)藥典[S].一部.2020: 86.

    [2] 韓路, 王家強(qiáng), 王海珍, 等.塔里木荒漠綠洲過(guò)渡帶主要種群生態(tài)位與空間格局分析 [J].植物科學(xué)學(xué)報(bào), 2016, 34(3): 352-360.

    [3] 李慧, 王歡, 宋經(jīng)元, 等.新疆地區(qū)藥用植物地理分布模式和氣候特征研究 [J].干旱區(qū)地理, 2015, 38(1): 36-42.

    [4] 張振巍, 張娜娜, 李月梅.不同品種甘草中總黃酮的含量考察 [J].中國(guó)藥師, 2013, 16(1): 49-51.

    [5] 李薇, 宋新波, 孫成榮, 等.三個(gè)不同品種甘草多糖的含量測(cè)定 [J].天津中醫(yī)藥, 2013, 30(1): 47-49.

    [6] 陸嘉惠.三種藥用甘草耐鹽性及耐鹽機(jī)制研究 [D].石河子: 石河子大學(xué), 2014.

    [7] Ellegren H, Galtier N.Determinants of genetic diversity [J]., 2016, 17(7): 422-433.

    [8] 彭羽, 卿鳳婷, 米凱, 等.生物多樣性不同層次尺度效應(yīng)及其耦合關(guān)系研究進(jìn)展 [J].生態(tài)學(xué)報(bào), 2015, 35(2): 577-583.

    [9] 熊海錚.豇豆遺傳多樣性及若干農(nóng)藝性狀關(guān)聯(lián)分析 [D].杭州: 浙江大學(xué), 2016.

    [10] Ge X J, Liu M H, Wang W K,.Population structure of wild bananas,, in China determined by SSR fingerprinting and cpDNA PCR-RFLP [J]., 2005, 14(4): 933-944.

    [11] 王剛, 曹佩, 韋學(xué)敏, 等.分子標(biāo)記技術(shù)在藥用植物種質(zhì)資源研究中的應(yīng)用 [J].中國(guó)現(xiàn)代中藥, 2019, 21(11): 1435-1444.

    [12] 胡一凡, 張雪梅, 石乃星, 等.云南草果種質(zhì)資源DNA條形碼序列分析 [J].中草藥, 2019, 50(24): 6091-6097.

    [13] 李文斌, 羅琳, 崔潔, 等.3種藥用甘草的遺傳多樣性研究進(jìn)展 [J].中藥材, 2019, 42(2): 463-469.

    [14] 王艷, 成金樂(lè).改進(jìn)的CTAB法提取32種中藥破壁飲片DNA及物種鑒定 [J].中國(guó)中醫(yī)藥現(xiàn)代遠(yuǎn)程教育, 2018, 16(22): 89-91,130.

    [15] 強(qiáng)正澤.基于ISSR標(biāo)記技術(shù)的甘肅地產(chǎn)藥材紅芪的遺傳多樣性研究 [D].蘭州: 甘肅中醫(yī)藥大學(xué), 2016.

    [16] 邱亞玲, 秦民堅(jiān), 戴軼群, 等.不同產(chǎn)地貓豆遺傳多樣性的ISSR分析 [J].中國(guó)實(shí)驗(yàn)方劑學(xué)雜志, 2015, 21(6): 86-90.

    [17] 趙娜, 魯萍, 李景欣.東北地區(qū)野豌豆屬四個(gè)種間遺傳關(guān)系的ISSR分析 [J].作物雜志, 2011(3): 93-96.

    [18] Shah A, Li D Z, Gao L M,.Genetic diversity within and among populations of the endangered species(Taxaceae) from Pakistan and implications for its conservation [J]., 2008, 36(3): 183-193.

    [19] 燕雪飛.中國(guó)野生大豆遺傳多樣性及其分化研究 [D].沈陽(yáng): 沈陽(yáng)農(nóng)業(yè)大學(xué), 2014.

    [20] Szczecińska M, Sramko G, Wo?osz K,.Genetic diversity and population structure of the rare and endangered plant species(L.) mill in east central Europe [J]., 2016, 11(3): e0151730.

    [21] Li E X, Yi S, Qiu Y X,.Phylogeography of two East Asian species in(Stemonaceae) inferred from chloroplast DNA and ISSR fingerprinting variation [J]., 2008, 49(3): 702-714.

    [22] Li B J, Wang J Y, Liu Z J,.Genetic diversity and ex situ conservation of, an endangered species endemic to China [J]., 2018, 19(1): 12.

    [23] 趙春磊.利用ISSR分析中國(guó)原生董棕植物遺傳多樣性 [D].廈門: 廈門大學(xué), 2009.

    [24] 吳濤, 陳少瑜, 寧德魯, 等.怒江州深紋核桃種質(zhì)資源SSR遺傳多樣性分析 [J].福建農(nóng)林大學(xué)學(xué)報(bào): 自然科學(xué)版, 2019, 48(2): 252-258.

    [25] 潘裕生, 方愛民.中國(guó)青藏高原特提斯的形成與演化 [J].地質(zhì)科學(xué), 2010, 45(1): 92-101.

    [26] 林暢松, 楊海軍, 劉景彥, 等.塔里木早古生代原盆地古隆起地貌和古地理格局與地層圈閉發(fā)育分布 [J].石油與天然氣地質(zhì), 2008, 29(2): 189-197.

    [27] 葛肖虹, 劉俊來(lái), 任收麥, 等.青藏高原隆升對(duì)中國(guó)構(gòu)造-地貌形成、氣候環(huán)境變遷與古人類遷徙的影響 [J].中國(guó)地質(zhì), 2014, 41(3): 698-714.

    [28] 馬志偉, 孫東懷, 易治宇, 等.晚新生代塔里木盆地西部沙漠形成與演化的粒度記錄 [J].蘭州大學(xué)學(xué)報(bào): 自然科學(xué)版, 2017, 53(2): 152-162.

    [29] 李錦軼.新疆東部新元古代晚期和古生代構(gòu)造格局及其演變 [J].地質(zhì)論評(píng), 2004, 50(3): 304-322.

    [30] 陸偉.基于cpDNA對(duì)大白刺(Nitraria roborowskii)的譜系地理研究 [D].石河子: 石河子大學(xué), 2015.

    [31] 蘇志豪, 卓立, 李文軍, 等.塔里木盆地沙生檉柳()遺傳多樣性 [J].中國(guó)沙漠, 2018, 38(5): 1041-1048.

    [32] 郭爭(zhēng)爭(zhēng), 程瑞婧, 楊麗, 等.不同繁殖模式下甘草藥用成分含量積累速度差異研究 [J].現(xiàn)代生物醫(yī)學(xué)進(jìn)展, 2015, 15(10): 1816-1819.

    [33] 謝良碧, 陸嘉惠, 李曉嵐, 等.三種甘草屬植物的種間雜交親和性及雜交種子活力 [J].植物分類與資源學(xué)報(bào), 2014, 36(3): 342-348.

    [34] 田潤(rùn)煒, 陸嘉惠, 謝良碧, 等.光果甘草與烏拉爾甘草開花與傳粉方式對(duì)生殖及種間關(guān)系的影響 [J].西北植物學(xué)報(bào), 2012, 32(10): 2004-2008.

    [35] 郭洪英.四川榿木優(yōu)樹群體遺傳多樣性分析及育種群體構(gòu)建 [D].北京: 北京林業(yè)大學(xué), 2019.

    [36] 俞文灝, 吳保鋒, 劉勇波.生境破碎化對(duì)動(dòng)植物遺傳多樣性的影響研究進(jìn)展 [J].應(yīng)用與環(huán)境生物學(xué)報(bào), 2019, 25(3): 743-749.

    [37] 馬芳.毛烏素沙地油蒿種群遺傳多樣性與環(huán)境相關(guān)性研究 [D].呼和浩特: 內(nèi)蒙古大學(xué), 2015.

    [38] 辛倩.基于轉(zhuǎn)錄組SSR分子標(biāo)記的甘草屬植物雜交區(qū)組成及種間基因漸滲研究 [D].石河子: 石河子大學(xué), 2019.

    [39] 宋雪梅.內(nèi)蒙古高原跨地帶分布紅砂種群遺傳多樣性研究 [D].呼和浩特: 內(nèi)蒙古大學(xué), 2012.

    [40] 張志南, 武高林, 王冬, 等.黃土高原半干旱區(qū)天然草地群落結(jié)構(gòu)與土壤水分關(guān)系 [J].草業(yè)學(xué)報(bào), 2014, 23(6): 313-319.

    Genetic diversity ISSR analysis offrom Xinjiang

    WANG Li-li1, ALIMJAN·PARHAT2, LI Li-wei1, JURAT·TOHTI2, LIU Zhong1

    1.School of Pharmacy, Shanghai Jiao Tong University, Shanghai 200240, China 2.College of Xinjiang Uyghur Medicine, Hotan 848000, China

    To evaluate the genetic diversity of, which is one of the source plants ofet, abundantly distributed in the Xinjiang.ISSR molecular marker method was used to analyze 167 individuals from six populations offrom different geographical locations of southern Xinjiang.In the six populations, 46 primers were used to amplify 193 polymorphic bands, accounting for 69.68% of the percentage of polymorphic loci (PPL).Moreover, Nei's genetic diversity index () was 0.272 2, and Shannon's diversity index () was 0.401 0.Nei's genetic diversity index at the level of total populations (t) was 0.513 0; The genetic differentiation coefficient (st) was 0.530 7, showing 0.4388 of the gene flow (m).In contrast, among these populations, PPLs ranged from 52.11% to 81.87%; the observed allele numbers (a) and the effective allele numbers (e) varied from 1.521 1 to 1.818 7 and from 1.307 9 to 1.498 9, respectively; Likewise, the averaged Nei's genetic diversity index () was 0.240 8, while the averaged Shannon's diversity index () was 0.3571.Finally, Mantel test showed that there was no correlation among genetic differentiations and geographical distances.Besides, the level of genetic diversity among populations showed a decreasing tendency of BC> KEL> PS> QM> AQS> KC.In cluster calculation, the six populations were gathered into three groups (clades) at the genetic correlation coefficient of 0.69.The PPL range of six populations was 52.11%—81.87%; The range ofaandewas 1.521 1—1.818 7 and 1.307 9—1.498 9, the averagevalue was 0.240 8, the averagevalue was0.357 1.It is well known thatpossesses comparatively rich genetic diversity.The genetic diversity of inter-populations was higher than that of intra-populations in this species.Habitat fragmentation was considered as the principle factor giving rise to the spatial distribution pattern of the genetic diversity, and the water content of habitat soils might be the crucial factor impacting the differentiation of the genetic diversity of the species.In sum, this study would provide a serviceable reference for a comprehensive and better understanding of the natural resource and their evolutionary potential, contributing to the conservation and sustainable utilization of such unadequately-studied medicinal plants.

    Bat.; ISSR; genetic diversity; genetic structure; natural resource conservation

    R282.12

    A

    0253 - 2670(2021)22 - 6975 - 08

    10.7501/j.issn.0253-2670.2021.22.024

    2021-05-06

    新疆維吾爾自治區(qū)區(qū)域協(xié)同創(chuàng)新專項(xiàng)(科技支疆項(xiàng)目)(2017E0230)

    王莉莉,碩士研究生,研究方向?yàn)樗帉W(xué)。Tel: (021)34208145 E-mail: liliwangz@sjtu.edu.cn

    通信作者:劉 忠,副教授,研究方向?yàn)樯帉W(xué)。Tel: (021)34208148 E-mail: liuzhong@sjtu.edu.cn

    #并列第一作者: 阿力木江·排爾哈提,講師,研究方向?yàn)樯锟茖W(xué)。Tel: 13899478919 E-mail: 652343507@qq.com

    [責(zé)任編輯 時(shí)圣明]

    99精品欧美一区二区三区四区| 国模一区二区三区四区视频 | 俺也久久电影网| 女人被狂操c到高潮| 精品国内亚洲2022精品成人| 在线观看www视频免费| 男女做爰动态图高潮gif福利片| 我要搜黄色片| 非洲黑人性xxxx精品又粗又长| 窝窝影院91人妻| 亚洲中文字幕日韩| 国产高清激情床上av| 丝袜美腿诱惑在线| 亚洲免费av在线视频| 中文字幕av在线有码专区| 久久婷婷成人综合色麻豆| 国产主播在线观看一区二区| 国产精品 国内视频| 久久九九热精品免费| av片东京热男人的天堂| 午夜精品一区二区三区免费看| 中文字幕高清在线视频| 亚洲真实伦在线观看| 国产一区二区激情短视频| 毛片女人毛片| 免费观看精品视频网站| 日韩欧美 国产精品| 亚洲18禁久久av| 一本久久中文字幕| 精品一区二区三区av网在线观看| 亚洲电影在线观看av| 国产免费av片在线观看野外av| 俄罗斯特黄特色一大片| 可以免费在线观看a视频的电影网站| 日韩欧美在线乱码| 人人妻,人人澡人人爽秒播| 亚洲国产欧美人成| 日本撒尿小便嘘嘘汇集6| 精品久久久久久久人妻蜜臀av| 不卡一级毛片| 长腿黑丝高跟| 久久精品国产亚洲av高清一级| videosex国产| 9191精品国产免费久久| 亚洲 欧美 日韩 在线 免费| 午夜福利18| 日韩精品免费视频一区二区三区| 色综合亚洲欧美另类图片| 黄片小视频在线播放| 国产精品一区二区三区四区久久| 十八禁网站免费在线| or卡值多少钱| 免费观看人在逋| 啪啪无遮挡十八禁网站| 日本 欧美在线| 丁香欧美五月| 亚洲男人天堂网一区| 免费在线观看完整版高清| 哪里可以看免费的av片| 99热6这里只有精品| 国产成人精品久久二区二区免费| 国产三级中文精品| 日韩欧美免费精品| 亚洲一区高清亚洲精品| 欧美绝顶高潮抽搐喷水| 丰满的人妻完整版| 99在线视频只有这里精品首页| 亚洲av成人精品一区久久| 日本黄色视频三级网站网址| 久久香蕉国产精品| 最新在线观看一区二区三区| 久久久久国产精品人妻aⅴ院| 搡老岳熟女国产| 69av精品久久久久久| 欧美日韩亚洲国产一区二区在线观看| 国产视频一区二区在线看| 国产精品一及| 日本一二三区视频观看| 黄色毛片三级朝国网站| 午夜免费成人在线视频| 日本a在线网址| www国产在线视频色| 99国产精品一区二区三区| 午夜免费成人在线视频| 日韩av在线大香蕉| 九色国产91popny在线| 国产精品亚洲一级av第二区| 亚洲精品美女久久久久99蜜臀| 制服诱惑二区| 国产av麻豆久久久久久久| 成人av一区二区三区在线看| 久久久精品欧美日韩精品| 久久久精品国产亚洲av高清涩受| 91在线观看av| 国产精品香港三级国产av潘金莲| 日韩大码丰满熟妇| 天堂av国产一区二区熟女人妻 | 久久精品综合一区二区三区| 美女 人体艺术 gogo| 波多野结衣巨乳人妻| 亚洲全国av大片| 亚洲国产欧美人成| 91麻豆av在线| 亚洲电影在线观看av| 长腿黑丝高跟| 女人被狂操c到高潮| 亚洲va日本ⅴa欧美va伊人久久| 少妇人妻一区二区三区视频| 精品久久久久久久末码| 啦啦啦观看免费观看视频高清| videosex国产| 欧美在线一区亚洲| 人妻夜夜爽99麻豆av| 2021天堂中文幕一二区在线观| 久久久久久久久久黄片| 搡老妇女老女人老熟妇| 免费在线观看亚洲国产| 久久久精品大字幕| 全区人妻精品视频| 成人三级黄色视频| 国产精品九九99| 俺也久久电影网| 亚洲免费av在线视频| 色综合婷婷激情| 男人舔奶头视频| 国产一区二区在线av高清观看| 国产av又大| 日韩大尺度精品在线看网址| 精品国产超薄肉色丝袜足j| 日韩精品免费视频一区二区三区| 丝袜美腿诱惑在线| 美女 人体艺术 gogo| 一级毛片精品| 国产精品免费视频内射| 日韩中文字幕欧美一区二区| 国产精品综合久久久久久久免费| 欧美日本亚洲视频在线播放| avwww免费| 俄罗斯特黄特色一大片| 香蕉丝袜av| 非洲黑人性xxxx精品又粗又长| 亚洲欧美日韩高清在线视频| 欧美性猛交╳xxx乱大交人| 欧美日韩乱码在线| 亚洲国产精品成人综合色| 无人区码免费观看不卡| 中文在线观看免费www的网站 | 麻豆成人av在线观看| 女同久久另类99精品国产91| 亚洲精品一区av在线观看| 亚洲av成人av| 18美女黄网站色大片免费观看| 少妇的丰满在线观看| 久久 成人 亚洲| 国产精品99久久99久久久不卡| 天天躁夜夜躁狠狠躁躁| 亚洲精品久久成人aⅴ小说| 国产亚洲欧美在线一区二区| 欧美成人一区二区免费高清观看 | 欧美极品一区二区三区四区| 欧美日韩国产亚洲二区| 好男人电影高清在线观看| 一区福利在线观看| 久久久国产成人免费| 国产亚洲精品一区二区www| 精华霜和精华液先用哪个| 这个男人来自地球电影免费观看| 国产精品亚洲美女久久久| 他把我摸到了高潮在线观看| 国产熟女午夜一区二区三区| 久久久久国产精品人妻aⅴ院| 亚洲,欧美精品.| 国产成人av教育| 亚洲专区字幕在线| 国产高清视频在线观看网站| 色综合欧美亚洲国产小说| 精品一区二区三区视频在线观看免费| 日韩欧美免费精品| 亚洲性夜色夜夜综合| 久久久国产成人精品二区| 一二三四在线观看免费中文在| 亚洲精品美女久久久久99蜜臀| 欧美丝袜亚洲另类 | 亚洲精品中文字幕在线视频| 日韩av在线大香蕉| 成人精品一区二区免费| 免费在线观看完整版高清| 欧美色视频一区免费| 色综合婷婷激情| 国产不卡一卡二| 欧美人与性动交α欧美精品济南到| 香蕉久久夜色| 99精品欧美一区二区三区四区| 免费在线观看亚洲国产| 熟女电影av网| 国产精品av久久久久免费| 精品久久久久久久久久久久久| 国产精品综合久久久久久久免费| 亚洲aⅴ乱码一区二区在线播放 | 搡老妇女老女人老熟妇| 一区二区三区国产精品乱码| 18禁国产床啪视频网站| 久久久久性生活片| 国产亚洲av高清不卡| 日本精品一区二区三区蜜桃| 又爽又黄无遮挡网站| 亚洲精品在线美女| 国产成人精品无人区| 日本一本二区三区精品| av国产免费在线观看| 真人一进一出gif抽搐免费| 日本熟妇午夜| 亚洲,欧美精品.| 国产69精品久久久久777片 | 国产av不卡久久| 亚洲av电影不卡..在线观看| 一区二区三区高清视频在线| a在线观看视频网站| 9191精品国产免费久久| 哪里可以看免费的av片| 午夜久久久久精精品| 亚洲精品在线观看二区| 国产欧美日韩一区二区精品| 又紧又爽又黄一区二区| 狠狠狠狠99中文字幕| 怎么达到女性高潮| 国产精品一及| www.精华液| 好男人电影高清在线观看| 午夜福利在线观看吧| 看黄色毛片网站| 国产欧美日韩一区二区三| 97超级碰碰碰精品色视频在线观看| 亚洲成人免费电影在线观看| 亚洲av第一区精品v没综合| 脱女人内裤的视频| 国产又色又爽无遮挡免费看| aaaaa片日本免费| 亚洲中文字幕日韩| 狠狠狠狠99中文字幕| 国产激情偷乱视频一区二区| 在线免费观看的www视频| 女人被狂操c到高潮| 两个人的视频大全免费| 国模一区二区三区四区视频 | 亚洲成人精品中文字幕电影| 热99re8久久精品国产| 午夜两性在线视频| 啪啪无遮挡十八禁网站| 亚洲最大成人中文| 国内毛片毛片毛片毛片毛片| 成人精品一区二区免费| 88av欧美| 久久久久精品国产欧美久久久| 亚洲美女黄片视频| 亚洲人成伊人成综合网2020| 免费在线观看影片大全网站| 日韩大码丰满熟妇| 亚洲精品一卡2卡三卡4卡5卡| 国产精品爽爽va在线观看网站| 国产人伦9x9x在线观看| 一a级毛片在线观看| 18禁黄网站禁片午夜丰满| 精品福利观看| 一夜夜www| 成人精品一区二区免费| 19禁男女啪啪无遮挡网站| 日韩欧美在线二视频| 国产三级在线视频| 国产成人一区二区三区免费视频网站| 亚洲电影在线观看av| 亚洲一码二码三码区别大吗| 国产成人精品久久二区二区免费| 国产成年人精品一区二区| 色综合欧美亚洲国产小说| 神马国产精品三级电影在线观看 | 精品人妻1区二区| 国产av又大| av福利片在线观看| cao死你这个sao货| 欧美激情久久久久久爽电影| 国内毛片毛片毛片毛片毛片| 91九色精品人成在线观看| 亚洲熟妇熟女久久| 中文字幕久久专区| 88av欧美| 最新美女视频免费是黄的| 国产av在哪里看| 一本一本综合久久| 夜夜看夜夜爽夜夜摸| 亚洲av熟女| 色综合婷婷激情| 波多野结衣高清作品| 亚洲国产中文字幕在线视频| 男人舔女人下体高潮全视频| 成年女人毛片免费观看观看9| 成人av一区二区三区在线看| а√天堂www在线а√下载| 91大片在线观看| 两个人的视频大全免费| 舔av片在线| 欧美另类亚洲清纯唯美| 亚洲自偷自拍图片 自拍| 在线观看一区二区三区| 午夜精品在线福利| 天堂av国产一区二区熟女人妻 | 大型av网站在线播放| 男女那种视频在线观看| 精品不卡国产一区二区三区| 亚洲专区字幕在线| 午夜a级毛片| 可以在线观看的亚洲视频| 啦啦啦韩国在线观看视频| 国产三级在线视频| 日本一二三区视频观看| netflix在线观看网站| 一a级毛片在线观看| 1024视频免费在线观看| 淫妇啪啪啪对白视频| 成人特级黄色片久久久久久久| 亚洲精品美女久久久久99蜜臀| 99精品久久久久人妻精品| 老司机深夜福利视频在线观看| 在线观看一区二区三区| 十八禁人妻一区二区| 亚洲av中文字字幕乱码综合| 日韩精品青青久久久久久| 岛国在线免费视频观看| 午夜成年电影在线免费观看| 动漫黄色视频在线观看| 亚洲国产欧美人成| 可以免费在线观看a视频的电影网站| 一二三四在线观看免费中文在| 美女扒开内裤让男人捅视频| 国产午夜福利久久久久久| 最近视频中文字幕2019在线8| www.精华液| 日日摸夜夜添夜夜添小说| 最新美女视频免费是黄的| 亚洲熟女毛片儿| 色综合站精品国产| 精品欧美一区二区三区在线| 精品一区二区三区视频在线观看免费| cao死你这个sao货| 中文字幕久久专区| 久久精品夜夜夜夜夜久久蜜豆 | 亚洲人与动物交配视频| 搡老熟女国产l中国老女人| 日韩 欧美 亚洲 中文字幕| 久久九九热精品免费| 夜夜躁狠狠躁天天躁| 亚洲自拍偷在线| 99国产精品一区二区蜜桃av| 国内精品久久久久精免费| 999久久久精品免费观看国产| 黄色视频不卡| 国内久久婷婷六月综合欲色啪| 久久欧美精品欧美久久欧美| av在线天堂中文字幕| 一本精品99久久精品77| 人妻夜夜爽99麻豆av| 国产黄色小视频在线观看| 久久国产精品人妻蜜桃| www.自偷自拍.com| 夜夜爽天天搞| 国产精品电影一区二区三区| 五月伊人婷婷丁香| 久久午夜亚洲精品久久| netflix在线观看网站| 色综合亚洲欧美另类图片| e午夜精品久久久久久久| 777久久人妻少妇嫩草av网站| av有码第一页| 欧美日本视频| 国产伦人伦偷精品视频| 9191精品国产免费久久| www.999成人在线观看| 国产成人欧美在线观看| 三级男女做爰猛烈吃奶摸视频| 一二三四社区在线视频社区8| 性欧美人与动物交配| 久久香蕉精品热| 亚洲欧美一区二区三区黑人| 99国产精品99久久久久| 黄色视频不卡| 久久热在线av| 91麻豆精品激情在线观看国产| 久久草成人影院| 久久精品91蜜桃| 国产高清有码在线观看视频 | 变态另类成人亚洲欧美熟女| 宅男免费午夜| 两人在一起打扑克的视频| 色哟哟哟哟哟哟| 亚洲欧美一区二区三区黑人| 日本一本二区三区精品| 免费看a级黄色片| 91国产中文字幕| 两性午夜刺激爽爽歪歪视频在线观看 | 国产成人系列免费观看| 欧美成狂野欧美在线观看| 国产伦一二天堂av在线观看| 亚洲av成人av| √禁漫天堂资源中文www| 成人国产一区最新在线观看| 人成视频在线观看免费观看| 精品久久久久久久人妻蜜臀av| 国产免费av片在线观看野外av| 成年免费大片在线观看| 日韩 欧美 亚洲 中文字幕| 国产成人啪精品午夜网站| 丰满的人妻完整版| 高潮久久久久久久久久久不卡| www.999成人在线观看| 久久 成人 亚洲| 舔av片在线| 搞女人的毛片| 久久天躁狠狠躁夜夜2o2o| 国产精品久久久久久久电影 | 中国美女看黄片| 夜夜看夜夜爽夜夜摸| 麻豆成人午夜福利视频| 亚洲精品一区av在线观看| 可以免费在线观看a视频的电影网站| 男人舔女人下体高潮全视频| 久久人妻av系列| 久久久水蜜桃国产精品网| 亚洲av中文字字幕乱码综合| 国产黄色小视频在线观看| 成年版毛片免费区| 成人手机av| 精品电影一区二区在线| 两人在一起打扑克的视频| 国产一区在线观看成人免费| 国产又黄又爽又无遮挡在线| 国产一区二区在线av高清观看| 亚洲午夜理论影院| 大型av网站在线播放| 欧美精品亚洲一区二区| 麻豆国产av国片精品| 日韩 欧美 亚洲 中文字幕| 丝袜人妻中文字幕| 老司机靠b影院| 操出白浆在线播放| 色哟哟哟哟哟哟| 嫁个100分男人电影在线观看| 色精品久久人妻99蜜桃| 人人妻人人看人人澡| 无遮挡黄片免费观看| 国产欧美日韩一区二区精品| 老司机在亚洲福利影院| 成在线人永久免费视频| 国产99久久九九免费精品| 久久婷婷人人爽人人干人人爱| 夜夜夜夜夜久久久久| 亚洲精品中文字幕在线视频| 欧美在线一区亚洲| 国产成+人综合+亚洲专区| АⅤ资源中文在线天堂| 国产人伦9x9x在线观看| 午夜福利欧美成人| 狂野欧美白嫩少妇大欣赏| 国产一区二区三区视频了| 中亚洲国语对白在线视频| 黄色毛片三级朝国网站| 免费在线观看完整版高清| 国产精品 欧美亚洲| 18禁黄网站禁片免费观看直播| aaaaa片日本免费| 日本免费一区二区三区高清不卡| 禁无遮挡网站| 男人舔女人的私密视频| 中出人妻视频一区二区| 88av欧美| 国产精品亚洲美女久久久| 1024视频免费在线观看| 91字幕亚洲| 可以免费在线观看a视频的电影网站| 国产精品久久久人人做人人爽| 一级毛片女人18水好多| 一夜夜www| 中文亚洲av片在线观看爽| 成人av在线播放网站| 精品久久久久久,| 中文字幕最新亚洲高清| 亚洲熟妇中文字幕五十中出| 9191精品国产免费久久| 色精品久久人妻99蜜桃| 午夜视频精品福利| 精品日产1卡2卡| 男女床上黄色一级片免费看| 成人永久免费在线观看视频| 99riav亚洲国产免费| 国产三级中文精品| 精品久久久久久成人av| 亚洲精品美女久久久久99蜜臀| 国产午夜精品论理片| 久久精品成人免费网站| 欧美又色又爽又黄视频| 国内久久婷婷六月综合欲色啪| 国产精品九九99| 757午夜福利合集在线观看| 欧美乱妇无乱码| 久久久国产欧美日韩av| 成人18禁高潮啪啪吃奶动态图| 国产一区二区在线av高清观看| 亚洲最大成人中文| tocl精华| 精华霜和精华液先用哪个| 精品久久久久久久久久免费视频| 我要搜黄色片| 亚洲中文字幕一区二区三区有码在线看 | 亚洲精华国产精华精| 久久久久久九九精品二区国产 | www.999成人在线观看| 成人永久免费在线观看视频| 男插女下体视频免费在线播放| 日韩精品青青久久久久久| √禁漫天堂资源中文www| 国产一区二区激情短视频| 又黄又爽又免费观看的视频| 好男人在线观看高清免费视频| 亚洲美女黄片视频| 久久久久久久久免费视频了| 久久人妻福利社区极品人妻图片| 久久精品91蜜桃| 首页视频小说图片口味搜索| 国产精品免费一区二区三区在线| 91九色精品人成在线观看| 国产成人影院久久av| 亚洲精品一区av在线观看| 欧美大码av| 丝袜美腿诱惑在线| 国产99白浆流出| 精品欧美一区二区三区在线| 国产av不卡久久| 色老头精品视频在线观看| 制服诱惑二区| 午夜福利欧美成人| 看片在线看免费视频| 亚洲五月天丁香| 伊人久久大香线蕉亚洲五| 国产精品永久免费网站| 久久久久久大精品| 欧美成人午夜精品| 亚洲人与动物交配视频| 色综合亚洲欧美另类图片| 国产亚洲精品av在线| 99riav亚洲国产免费| 国产野战对白在线观看| 国产成人av激情在线播放| 日本免费a在线| 人妻夜夜爽99麻豆av| 欧美日韩亚洲国产一区二区在线观看| 国产精品九九99| 一二三四在线观看免费中文在| 免费在线观看成人毛片| 少妇被粗大的猛进出69影院| 在线观看舔阴道视频| 午夜免费激情av| 成人三级黄色视频| 国产激情欧美一区二区| 久久天躁狠狠躁夜夜2o2o| www.自偷自拍.com| 国产aⅴ精品一区二区三区波| 在线观看美女被高潮喷水网站 | 禁无遮挡网站| 五月伊人婷婷丁香| 国产1区2区3区精品| 国模一区二区三区四区视频 | 日本a在线网址| 老熟妇仑乱视频hdxx| 丰满人妻熟妇乱又伦精品不卡| 亚洲av电影在线进入| 欧美三级亚洲精品| 亚洲av电影不卡..在线观看| 久久久久久大精品| 欧美三级亚洲精品| 欧美日韩亚洲国产一区二区在线观看| 午夜福利免费观看在线| 夜夜看夜夜爽夜夜摸| 精品日产1卡2卡| 欧美性猛交黑人性爽| 黑人巨大精品欧美一区二区mp4| 久久久久久亚洲精品国产蜜桃av| 亚洲熟妇中文字幕五十中出| 国产乱人伦免费视频| e午夜精品久久久久久久| www.999成人在线观看| 午夜激情av网站| 免费在线观看成人毛片| 国产高清视频在线观看网站| 国产乱人伦免费视频| 最新在线观看一区二区三区| 亚洲精品美女久久av网站| 一二三四在线观看免费中文在| 欧美3d第一页| 久久久水蜜桃国产精品网| 久久国产精品影院| 免费看日本二区| 国产午夜精品久久久久久| 精品不卡国产一区二区三区| 中出人妻视频一区二区| 欧美性猛交黑人性爽| 国产午夜福利久久久久久| 国产成人啪精品午夜网站| 精品国产乱码久久久久久男人| 欧美一区二区精品小视频在线| 丰满人妻熟妇乱又伦精品不卡| 一个人免费在线观看电影 | 最近视频中文字幕2019在线8| 亚洲精品中文字幕在线视频| 欧美国产日韩亚洲一区| 欧美性猛交╳xxx乱大交人| 亚洲成人久久性|