• <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王莉莉阿力木江排爾哈提李莉偉居來提托合提
    中草藥 2021年22期

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

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

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

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

    揭示新疆塔里木盆地分布的野生脹果甘草的遺傳多樣性和種群遺傳結(jié)構(gòu)。利用ISSR分子標記方法,對來自新疆南疆不同地理方位的脹果甘草6個種群總計167個個體進行分析。對脹果甘草6個種群共篩選出46條引物,擴增獲得193條多態(tài)條帶,多態(tài)位點比率(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檢驗顯示種群間的遺傳分化與地理距離無相關(guān)性;種群間遺傳變異豐度和遺傳多樣性水平呈現(xiàn)BC>KEL>PS>QM>AQS>KC的遞減趨勢,聚類分析將6個種群在遺傳相似性系數(shù)為0.69處分為3類;6個種群內(nèi)個體間PPL的變化范圍為52.11%~81.87%,觀測等位基因數(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é)果對全面了解和掌握脹果甘草自然資源現(xiàn)狀、物種進化潛力等,制定資源利用與保護的正確策略提供了重要的研究基礎(chǔ)。

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

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

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

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

    1 材料和儀器

    1.1 材料

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

    1.2 儀器

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

    表1 脹果甘草的樣品來源信息

    2 方法

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

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

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

    2.2 引物的來源

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

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

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

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

    2.4 新疆脹果甘草遺傳多樣性ISSR檢測

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

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

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

    3 結(jié)果與分析

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

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

    3.1.2 循環(huán)次數(shù)的優(yōu)化 對脹果甘草的樣品進行ISSR擴增時,當PCR反應(yīng)的循環(huán)次數(shù)設(shè)置為30次時,電泳檢測結(jié)果顯示沒有條帶,進一步將循環(huán)次數(shù)設(shè)置為36次和38次,檢測結(jié)果無明顯差距,而循環(huán)次數(shù)為40次時,得到了條帶清晰度高、多態(tài)性較豐富的擴增結(jié)果,因此,40次擴增循環(huán)次數(shù)是進行脹果甘草葉片基因組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 脹果甘草的遺傳多樣性檢測

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

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

    3.2.2 種群內(nèi)個體間的遺傳多樣性 對脹果甘草6個種群進行遺傳結(jié)構(gòu)分析的結(jié)果顯示(圖3和表3),6個種群的PPL的變化范圍在52.11%~81.87%,觀測等位基因數(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ū)庫車市種群(KC)的變異相對貧乏,遺傳多樣性水平較低,群體間遺傳變異豐度和遺傳多樣性水平呈現(xiàn)BC>KEL>PS>QM>AQS>KC的遞減趨勢。

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

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

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

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

    根據(jù)6個脹果甘草種群間的Nei’s遺傳距離和遺傳一致度(表4)利用NTSYSpc 2.10e軟件UPGMA方法對脹果甘草6個種群的聚類分析(圖4)表明,6個種群在遺傳相似性系數(shù)為0.69處分成3個大類。首先,阿克蘇地區(qū)拜城縣種群(BC)單獨成為一類,表現(xiàn)出與來自其他5個地理分布區(qū)的脹果甘草種群之間較少的遺傳相似性,說明該種群與其他種群之間的遺傳分化程度較大;其次,阿克蘇地區(qū)庫車市種群(KC)與阿克蘇市種群(AQS)聚為第2類;第3,巴音郭楞蒙古州的庫爾勒市種群(KEL)與且末縣種群(QM)的遺傳距離最近,兩者先聚在一起,然后再進一步與和田地區(qū)皮山縣種群(PS)聚合為第3類。Mantel檢驗結(jié)果顯示群體間的遺傳距離與地理距離之間無相關(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]。遺傳多樣性豐富的物種具有更多的等位基因,對環(huán)境變化具有更好的適應(yīng)性和更強大的抵御能力[22]。影響物種遺傳多樣性和種群遺傳分化水平因素包括了地理隔離、生境片段化等外部因素,與繁育系統(tǒng)、遺傳漂變、基因突變和基因流等內(nèi)部因素[23-24]。

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

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

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

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

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

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

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

    [1] 中國藥典[S].一部.2020: 86.

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

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

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

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

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

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

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

    [9] 熊海錚.豇豆遺傳多樣性及若干農(nóng)藝性狀關(guān)聯(lián)分析 [D].杭州: 浙江大學, 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] 王剛, 曹佩, 韋學敏, 等.分子標記技術(shù)在藥用植物種質(zhì)資源研究中的應(yīng)用 [J].中國現(xiàn)代中藥, 2019, 21(11): 1435-1444.

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

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

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

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

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

    [17] 趙娜, 魯萍, 李景欣.東北地區(qū)野豌豆屬四個種間遺傳關(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] 燕雪飛.中國野生大豆遺傳多樣性及其分化研究 [D].沈陽: 沈陽農(nóng)業(yè)大學, 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分析中國原生董棕植物遺傳多樣性 [D].廈門: 廈門大學, 2009.

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

    [40] 張志南, 武高林, 王冬, 等.黃土高原半干旱區(qū)天然草地群落結(jié)構(gòu)與土壤水分關(guān)系 [J].草業(yè)學報, 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)新專項(科技支疆項目)(2017E0230)

    王莉莉,碩士研究生,研究方向為藥學。Tel: (021)34208145 E-mail: liliwangz@sjtu.edu.cn

    通信作者:劉 忠,副教授,研究方向為生藥學。Tel: (021)34208148 E-mail: liuzhong@sjtu.edu.cn

    #并列第一作者: 阿力木江·排爾哈提,講師,研究方向為生物科學。Tel: 13899478919 E-mail: 652343507@qq.com

    [責任編輯 時圣明]

    天堂动漫精品| 亚洲av不卡在线观看| 淫妇啪啪啪对白视频| 手机成人av网站| 中文字幕av成人在线电影| 国产精品,欧美在线| 久久久久久久精品吃奶| 免费看日本二区| 成人一区二区视频在线观看| 日本免费一区二区三区高清不卡| 久久精品影院6| 一区二区三区免费毛片| 国产精品久久电影中文字幕| 99国产精品一区二区三区| 国产av在哪里看| 久久香蕉国产精品| 操出白浆在线播放| 国产色爽女视频免费观看| 国产精品98久久久久久宅男小说| 国内揄拍国产精品人妻在线| 欧美zozozo另类| 亚洲精品456在线播放app | 亚洲乱码一区二区免费版| 国产欧美日韩一区二区三| 免费搜索国产男女视频| 丰满人妻熟妇乱又伦精品不卡| 国内精品一区二区在线观看| 搡女人真爽免费视频火全软件 | 18+在线观看网站| 久久久久久九九精品二区国产| 国产高清视频在线观看网站| 亚洲人成电影免费在线| 窝窝影院91人妻| 精品一区二区三区人妻视频| 色综合站精品国产| 在线观看免费午夜福利视频| 免费av观看视频| 国产乱人视频| 欧美一区二区亚洲| 久久久国产成人精品二区| netflix在线观看网站| 大型黄色视频在线免费观看| 女生性感内裤真人,穿戴方法视频| 黄片大片在线免费观看| 九色国产91popny在线| 女同久久另类99精品国产91| 长腿黑丝高跟| 免费在线观看亚洲国产| 日韩国内少妇激情av| 亚洲av成人不卡在线观看播放网| 亚洲av电影在线进入| 少妇人妻精品综合一区二区 | 久久久久国产精品人妻aⅴ院| 亚洲国产精品999在线| 欧美一区二区亚洲| 国产成人啪精品午夜网站| 国产黄片美女视频| 在线免费观看的www视频| 最后的刺客免费高清国语| 国产精品 欧美亚洲| 亚洲精品一卡2卡三卡4卡5卡| 搡女人真爽免费视频火全软件 | 国产成+人综合+亚洲专区| 男女午夜视频在线观看| 日本a在线网址| 美女高潮的动态| 亚洲天堂国产精品一区在线| 欧美最新免费一区二区三区 | 99久久精品热视频| 最近最新免费中文字幕在线| 高潮久久久久久久久久久不卡| 国产极品精品免费视频能看的| 精品乱码久久久久久99久播| 国产成人av教育| 看黄色毛片网站| 男人舔奶头视频| 国产一区二区三区在线臀色熟女| 国产v大片淫在线免费观看| 9191精品国产免费久久| 午夜福利在线在线| 69av精品久久久久久| 蜜桃久久精品国产亚洲av| 国产精品一区二区三区四区免费观看 | 18禁美女被吸乳视频| 免费看日本二区| 国产精品国产高清国产av| 午夜两性在线视频| 欧美黄色片欧美黄色片| 亚洲精品久久国产高清桃花| 蜜桃亚洲精品一区二区三区| 午夜两性在线视频| 国产一区二区三区在线臀色熟女| 一卡2卡三卡四卡精品乱码亚洲| 少妇高潮的动态图| 久久中文看片网| 亚洲av免费高清在线观看| 怎么达到女性高潮| 国产野战对白在线观看| 日日干狠狠操夜夜爽| 欧美日韩亚洲国产一区二区在线观看| 一级黄色大片毛片| 亚洲国产欧美人成| 99久久综合精品五月天人人| 国产极品精品免费视频能看的| 亚洲精品粉嫩美女一区| 欧美一级毛片孕妇| 成年女人看的毛片在线观看| 嫩草影院精品99| 国产成人av教育| 成人永久免费在线观看视频| 国产91精品成人一区二区三区| 青草久久国产| 美女被艹到高潮喷水动态| 久久久久久九九精品二区国产| 亚洲av成人av| www国产在线视频色| 亚洲va日本ⅴa欧美va伊人久久| 国产高清激情床上av| 一本久久中文字幕| 日日摸夜夜添夜夜添小说| 在线观看av片永久免费下载| 成人午夜高清在线视频| 国产色爽女视频免费观看| 我的老师免费观看完整版| 一二三四社区在线视频社区8| 国产精品美女特级片免费视频播放器| 久久精品国产亚洲av涩爱 | 免费人成视频x8x8入口观看| 亚洲av成人精品一区久久| 色哟哟哟哟哟哟| 久久香蕉国产精品| 国产v大片淫在线免费观看| 国产黄片美女视频| 女生性感内裤真人,穿戴方法视频| 国产视频内射| 黄色成人免费大全| 国产一区二区亚洲精品在线观看| 免费在线观看成人毛片| 老司机深夜福利视频在线观看| 亚洲精品成人久久久久久| 夜夜爽天天搞| 可以在线观看毛片的网站| av天堂在线播放| 在线观看美女被高潮喷水网站 | 欧美一级a爱片免费观看看| 精品人妻偷拍中文字幕| 国产主播在线观看一区二区| tocl精华| 久久精品亚洲精品国产色婷小说| 深夜精品福利| 可以在线观看的亚洲视频| 国产熟女xx| 男人舔女人下体高潮全视频| 精品一区二区三区视频在线 | 男女视频在线观看网站免费| 在线播放国产精品三级| 一级黄色大片毛片| 亚洲av电影在线进入| 国产v大片淫在线免费观看| 在线观看66精品国产| 精品乱码久久久久久99久播| 国产精品亚洲av一区麻豆| 在线十欧美十亚洲十日本专区| 在线观看一区二区三区| 亚洲人成电影免费在线| 老熟妇乱子伦视频在线观看| 99精品欧美一区二区三区四区| 老司机午夜十八禁免费视频| 变态另类成人亚洲欧美熟女| 国产精品爽爽va在线观看网站| 精品久久久久久成人av| 成人av一区二区三区在线看| 偷拍熟女少妇极品色| 校园春色视频在线观看| 此物有八面人人有两片| 久久天躁狠狠躁夜夜2o2o| 一级毛片女人18水好多| avwww免费| 最后的刺客免费高清国语| 可以在线观看毛片的网站| 亚洲精品日韩av片在线观看 | 内射极品少妇av片p| 久久99热这里只有精品18| 欧美成狂野欧美在线观看| 9191精品国产免费久久| 久久这里只有精品中国| 国产日本99.免费观看| 成年版毛片免费区| 久久久久久国产a免费观看| 亚洲精品456在线播放app | 欧美日韩精品网址| 国产一区二区激情短视频| www国产在线视频色| 欧美区成人在线视频| 成年女人永久免费观看视频| 国产久久久一区二区三区| 精品人妻偷拍中文字幕| 69av精品久久久久久| 国产色婷婷99| 热99在线观看视频| av片东京热男人的天堂| 国产一级毛片七仙女欲春2| 免费观看人在逋| 国产 一区 欧美 日韩| 好男人电影高清在线观看| 亚洲aⅴ乱码一区二区在线播放| 日本一本二区三区精品| 久久人人精品亚洲av| 最后的刺客免费高清国语| 少妇的逼水好多| 中文字幕人成人乱码亚洲影| 一进一出抽搐动态| 亚洲成人免费电影在线观看| 51午夜福利影视在线观看| 真实男女啪啪啪动态图| 欧美在线黄色| 嫁个100分男人电影在线观看| 少妇的逼水好多| 亚洲精品成人久久久久久| 毛片女人毛片| 岛国在线观看网站| 天天添夜夜摸| 99热精品在线国产| 日韩欧美精品v在线| 青草久久国产| 亚洲av五月六月丁香网| 极品教师在线免费播放| 免费观看精品视频网站| 搡女人真爽免费视频火全软件 | 三级国产精品欧美在线观看| 天天躁日日操中文字幕| 韩国av一区二区三区四区| 国产私拍福利视频在线观看| 欧美成人免费av一区二区三区| 亚洲电影在线观看av| 99国产综合亚洲精品| 亚洲人成网站高清观看| 一级a爱片免费观看的视频| 黄片大片在线免费观看| 亚洲精品一卡2卡三卡4卡5卡| 久久久久久大精品| 亚洲欧美日韩高清在线视频| 国产成人系列免费观看| 亚洲人成电影免费在线| 亚洲电影在线观看av| 观看美女的网站| 99热这里只有是精品50| 久久久久久久久大av| 97超级碰碰碰精品色视频在线观看| 久久国产精品影院| 亚洲av第一区精品v没综合| 99热精品在线国产| 久久中文看片网| 国产精品久久久久久久久免 | 搡老岳熟女国产| 欧美乱妇无乱码| 欧美激情久久久久久爽电影| 国产一区二区三区在线臀色熟女| 中文字幕久久专区| 91av网一区二区| 欧美色视频一区免费| 欧美不卡视频在线免费观看| 搡老妇女老女人老熟妇| 中文字幕高清在线视频| 老汉色av国产亚洲站长工具| 午夜久久久久精精品| 亚洲乱码一区二区免费版| 久久久国产成人免费| 亚洲av二区三区四区| 90打野战视频偷拍视频| 最新美女视频免费是黄的| 日本与韩国留学比较| 一本一本综合久久| 九九在线视频观看精品| 丝袜美腿在线中文| 国产伦在线观看视频一区| 亚洲精品色激情综合| 成人亚洲精品av一区二区| av天堂在线播放| 午夜精品在线福利| 欧美最新免费一区二区三区 | 久久性视频一级片| av天堂中文字幕网| 夜夜爽天天搞| 亚洲国产精品合色在线| 观看美女的网站| 亚洲中文日韩欧美视频| 99热只有精品国产| 级片在线观看| 国产伦精品一区二区三区视频9 | 日韩免费av在线播放| 一区二区三区免费毛片| 激情在线观看视频在线高清| 欧美日韩乱码在线| 免费av观看视频| 少妇的逼好多水| svipshipincom国产片| 亚洲真实伦在线观看| 精品一区二区三区视频在线 | 热99re8久久精品国产| 国产精品久久久久久久电影 | 国产老妇女一区| 国产精品爽爽va在线观看网站| 真实男女啪啪啪动态图| 国产精品免费一区二区三区在线| 欧美成人免费av一区二区三区| 日韩精品青青久久久久久| 十八禁人妻一区二区| 在线播放无遮挡| 哪里可以看免费的av片| 国产色爽女视频免费观看| av欧美777| 首页视频小说图片口味搜索| 99久久久亚洲精品蜜臀av| 久久久久九九精品影院| 欧美区成人在线视频| 母亲3免费完整高清在线观看| 成人永久免费在线观看视频| 制服丝袜大香蕉在线| 国产av麻豆久久久久久久| 亚洲国产欧洲综合997久久,| 在线播放国产精品三级| 久久久久久久久大av| 色av中文字幕| 男人舔女人下体高潮全视频| 欧美成人免费av一区二区三区| 亚洲美女黄片视频| 一级黄片播放器| 悠悠久久av| 久久性视频一级片| 精品国产亚洲在线| 久久九九热精品免费| 在线天堂最新版资源| 69av精品久久久久久| 日本熟妇午夜| 超碰av人人做人人爽久久 | 日本精品一区二区三区蜜桃| 观看免费一级毛片| 香蕉av资源在线| 欧美中文综合在线视频| eeuss影院久久| xxxwww97欧美| 村上凉子中文字幕在线| 在线观看免费午夜福利视频| 999久久久精品免费观看国产| 色播亚洲综合网| 成人永久免费在线观看视频| 制服丝袜大香蕉在线| 欧美乱码精品一区二区三区| 亚洲av电影在线进入| 日日摸夜夜添夜夜添小说| 哪里可以看免费的av片| 国产精品电影一区二区三区| 丰满乱子伦码专区| 日韩欧美在线乱码| 舔av片在线| 亚洲国产日韩欧美精品在线观看 | 欧美日韩精品网址| 嫩草影院入口| 免费av不卡在线播放| 一个人观看的视频www高清免费观看| 欧美绝顶高潮抽搐喷水| 日日干狠狠操夜夜爽| 内射极品少妇av片p| 听说在线观看完整版免费高清| 麻豆一二三区av精品| 长腿黑丝高跟| 美女被艹到高潮喷水动态| 麻豆国产97在线/欧美| 免费观看精品视频网站| 亚洲精品国产精品久久久不卡| 成人欧美大片| 亚洲无线在线观看| 国产精品亚洲美女久久久| 听说在线观看完整版免费高清| 色综合婷婷激情| 淫秽高清视频在线观看| 中文资源天堂在线| 日本免费a在线| 亚洲国产日韩欧美精品在线观看 | 亚洲av二区三区四区| 久久久久久久精品吃奶| 麻豆成人av在线观看| 99久久无色码亚洲精品果冻| 黄色日韩在线| 十八禁人妻一区二区| 99视频精品全部免费 在线| a在线观看视频网站| 三级男女做爰猛烈吃奶摸视频| 精品久久久久久久末码| 国产亚洲精品久久久久久毛片| 亚洲五月天丁香| 亚洲av二区三区四区| 一边摸一边抽搐一进一小说| 亚洲国产高清在线一区二区三| 亚洲成av人片免费观看| 两个人的视频大全免费| 国产美女午夜福利| 亚洲真实伦在线观看| 成年免费大片在线观看| 在线观看一区二区三区| 黄色丝袜av网址大全| 午夜免费观看网址| 国产成人av激情在线播放| 99热只有精品国产| 真人做人爱边吃奶动态| 亚洲国产欧美人成| 在线观看舔阴道视频| 国内毛片毛片毛片毛片毛片| 欧美精品啪啪一区二区三区| 中国美女看黄片| 午夜两性在线视频| 欧美日韩福利视频一区二区| 婷婷精品国产亚洲av在线| 非洲黑人性xxxx精品又粗又长| 国产精品影院久久| www国产在线视频色| 最新中文字幕久久久久| 中文在线观看免费www的网站| 91麻豆精品激情在线观看国产| 久久人妻av系列| 深夜精品福利| 在线国产一区二区在线| 午夜亚洲福利在线播放| av福利片在线观看| 久久性视频一级片| 日韩精品青青久久久久久| 日本撒尿小便嘘嘘汇集6| 男女午夜视频在线观看| 国产精品久久电影中文字幕| 国产91精品成人一区二区三区| 亚洲在线观看片| 久久香蕉国产精品| 国产熟女xx| 久久久久久久精品吃奶| 中文亚洲av片在线观看爽| 久9热在线精品视频| 一区二区三区高清视频在线| xxxwww97欧美| 少妇丰满av| av黄色大香蕉| 国产av不卡久久| 成人精品一区二区免费| 男女之事视频高清在线观看| 老熟妇乱子伦视频在线观看| 亚洲成av人片免费观看| 97碰自拍视频| 欧美日韩亚洲国产一区二区在线观看| 91麻豆精品激情在线观看国产| 日韩大尺度精品在线看网址| 久久久精品欧美日韩精品| www日本在线高清视频| 亚洲av第一区精品v没综合| 白带黄色成豆腐渣| 亚洲成av人片免费观看| 男人舔奶头视频| 婷婷丁香在线五月| 欧美午夜高清在线| 真人做人爱边吃奶动态| 搡老妇女老女人老熟妇| 黄色日韩在线| 久99久视频精品免费| 日韩欧美国产一区二区入口| 男插女下体视频免费在线播放| 国产黄片美女视频| 亚洲最大成人手机在线| 亚洲精品一区av在线观看| 亚洲欧美精品综合久久99| 国产私拍福利视频在线观看| 99久久精品一区二区三区| 日本免费a在线| 一个人看的www免费观看视频| 国产精品嫩草影院av在线观看 | 中文字幕av在线有码专区| 99国产综合亚洲精品| 激情在线观看视频在线高清| 真人一进一出gif抽搐免费| 男女那种视频在线观看| 国产精品爽爽va在线观看网站| 在线十欧美十亚洲十日本专区| 18禁裸乳无遮挡免费网站照片| 国内精品美女久久久久久| av女优亚洲男人天堂| 最好的美女福利视频网| www.熟女人妻精品国产| www.999成人在线观看| 亚洲成av人片在线播放无| 岛国在线免费视频观看| 三级男女做爰猛烈吃奶摸视频| 亚洲人成电影免费在线| 搡老岳熟女国产| 99热这里只有精品一区| 2021天堂中文幕一二区在线观| 成人亚洲精品av一区二区| 国产成年人精品一区二区| 午夜日韩欧美国产| 国产色爽女视频免费观看| 老汉色av国产亚洲站长工具| 亚洲国产日韩欧美精品在线观看 | 国产精品综合久久久久久久免费| 色视频www国产| 草草在线视频免费看| 怎么达到女性高潮| 最近视频中文字幕2019在线8| 国产精品98久久久久久宅男小说| 欧美成人一区二区免费高清观看| 精品电影一区二区在线| 国产v大片淫在线免费观看| 岛国视频午夜一区免费看| 免费人成在线观看视频色| 最近视频中文字幕2019在线8| av福利片在线观看| 国产成人aa在线观看| 在线观看一区二区三区| 日韩大尺度精品在线看网址| 亚洲七黄色美女视频| 九九在线视频观看精品| 99久久精品国产亚洲精品| 午夜a级毛片| av专区在线播放| 精品免费久久久久久久清纯| 亚洲 欧美 日韩 在线 免费| 在线免费观看不下载黄p国产 | 无遮挡黄片免费观看| svipshipincom国产片| 在线免费观看的www视频| 国产亚洲精品综合一区在线观看| 欧美激情在线99| 老汉色∧v一级毛片| 国产精品久久久久久久久免 | 成年女人永久免费观看视频| 欧美日韩瑟瑟在线播放| 欧美色视频一区免费| 亚洲乱码一区二区免费版| 国产色爽女视频免费观看| 乱人视频在线观看| 亚洲 欧美 日韩 在线 免费| 欧美中文日本在线观看视频| 97超视频在线观看视频| 国产午夜精品论理片| 性色avwww在线观看| 波多野结衣高清无吗| 亚洲国产精品久久男人天堂| 国产亚洲欧美在线一区二区| 国产精品98久久久久久宅男小说| 一进一出抽搐动态| 黄色女人牲交| 一级毛片高清免费大全| 欧美成人一区二区免费高清观看| 亚洲av成人不卡在线观看播放网| 中文字幕人妻熟人妻熟丝袜美 | 波野结衣二区三区在线 | 国内精品美女久久久久久| 久久久国产成人精品二区| 国产精品一区二区三区四区久久| 伊人久久精品亚洲午夜| 欧美日韩一级在线毛片| 国产真人三级小视频在线观看| 国产精品电影一区二区三区| 91久久精品国产一区二区成人 | 给我免费播放毛片高清在线观看| 国产av在哪里看| 黑人欧美特级aaaaaa片| 深爱激情五月婷婷| 成年女人看的毛片在线观看| 国产单亲对白刺激| 亚洲熟妇中文字幕五十中出| 欧美日韩瑟瑟在线播放| 51国产日韩欧美| 国产精品久久久人人做人人爽| 中文资源天堂在线| 免费观看的影片在线观看| 美女大奶头视频| 怎么达到女性高潮| 欧美成人a在线观看| 男人舔奶头视频| 国产 一区 欧美 日韩| 叶爱在线成人免费视频播放| 动漫黄色视频在线观看| 女人被狂操c到高潮| 国产国拍精品亚洲av在线观看 | 香蕉丝袜av| 国产99白浆流出| 熟女少妇亚洲综合色aaa.| 久久精品综合一区二区三区| 久久香蕉精品热| 性欧美人与动物交配| 亚洲av成人av| 免费无遮挡裸体视频| 中文字幕精品亚洲无线码一区| 99久久精品热视频| 国产精品99久久久久久久久| 久久午夜亚洲精品久久| 在线观看一区二区三区| 中文字幕高清在线视频| 男插女下体视频免费在线播放| 亚洲国产精品sss在线观看| 看免费av毛片| 老鸭窝网址在线观看| 国产野战对白在线观看| 怎么达到女性高潮| 亚洲精华国产精华精| 国产男靠女视频免费网站| 看免费av毛片| 一边摸一边抽搐一进一小说| 村上凉子中文字幕在线| av专区在线播放| 免费看十八禁软件| 亚洲欧美日韩无卡精品| 国产精品电影一区二区三区| 日韩欧美精品免费久久 | 51午夜福利影视在线观看| 久久久久久九九精品二区国产| svipshipincom国产片| 国产伦人伦偷精品视频|