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

    倒卵葉五加藥材的基原鑒定研究

    2021-09-03 03:37:15趙新禮舒泉湧崔蕾朱愛華
    安徽農(nóng)業(yè)科學(xué) 2021年15期

    趙新禮 舒泉湧 崔蕾 朱愛華

    摘要 [目的]應(yīng)用DNA條形碼鑒定技術(shù)對短柄五加等五加屬植物進行研究,以弄清倒卵葉五加和短柄五加的物種關(guān)系。[方法]通過對館藏標本與現(xiàn)地植物形態(tài)進行觀測,采用通用引物 ITS2 和 psbA-trnH對采自陜甘寧地區(qū)的短柄五加等50個樣本進行擴增后雙向測序,并從GenBank數(shù)據(jù)庫下載了五加屬等相關(guān)植物的序列,分別通過ITS2 database或其序列注釋,獲得ITS2和psbA-trnH序列,隨后進行BLAST比對鑒定分析;用MEGA 6.06軟件對所得序列進行比對,分析變異位點,計算GC含量、種內(nèi)和種間遺傳距離,構(gòu)建NJ系統(tǒng)發(fā)育樹。[結(jié)果]ITS2 和 psbA-trnH 序列的 PCR 擴增和測序成功率均為100%;倒卵葉五加和短柄五加的植物形態(tài)、ITS2序列和psbA-trnH序列完全相同。ITS2 序列:短柄五加和糙葉五加僅有一個堿基的差異,蜀五加與藤五加序列相同;psbA-trnH 序列:短柄五加與糙葉五加的差異顯著,短柄五加與蜀五加序列相同。[結(jié)論]同時使用ITS2、psbA-trnH這2個標記的DNA 條形碼可以準確鑒別短柄五加等五加屬植物;短柄五加與倒卵葉五加為同一種植物。

    關(guān)鍵詞 倒卵葉五加;短柄五加;基原鑒定;DNA條形碼;ITS2;psbA-trnH

    中圖分類號 R 282.5 ?文獻標識碼 A ?文章編號 0517-6611(2021)15-0166-07

    Abstract [Objective]The molecular identification technique of DNA barcode was used to identify the species of several plants of ?Eleutherococcus ?Maxim. in order to understand the relationship between ?Eleutherococcus brachypus ?and ?Eleutherococcus obovatus . [Method]By observing the morphology of collection specimens and onsite plants, using universal primers ITS2 and psbAtrnH, 50 samples including ?Eleutherococcus brachypus ?collected from Shaanxi, Gansu and Ningxia were amplified and sequenced in both directions,and downloaded the sequence of Araliacea and other related plants from GenBank database, obtained ITS2 and psbAtrnH sequence through ITS2 database or its sequence annotation, and then performed BLAST comparison and identification analysis.MEGA 6.06 software was used to compare the obtained sequences, analyze the mutation sites, calculate the GC content, intra and interspecies genetic distance, and construct the NJ phylogenetic tree.[Result]The success rate of PCR amplification and sequencing of ITS2 and psbAtrnH was 100%.The plant morphology, ITS2 sequence and psbAtrnH sequence of ?Eleutherococcus brachypus ?and ?Eleutherococcus obovatus ?were identical. ITS2 sequence: ?Eleutherococcus brachypus ?and ?Eleutherococcus henryi ?had only one base difference, and the sequence of ?Eleutherococcus leucorrhizus ?var. ?setchuenensis ?was the same as that of ?Eleutherococcus leucorrhizus .The sequence of psbAtrnH: ?Eleutherococcus brachypus ?and ?Eleutherococcus henryi ?sequence were significantly different, and the sequence of ?Eleutherococcus brachypus ?and ?Eleutherococcus leucorrhizus ?var. ?setchuenensis ?was the same.[Conclusion]At the same time, using ITS2 and psbAtrnH genetic markers, DNA barcode technology can accurately identify the 4 species of ?Eleutherococcus ?Maxim.; ?E. brachypus ?and ?E. obovatus ?should be the same plant.

    Key words ?Eleutherococcus obovatus;Eleutherococcus brachypus ;Base source identification;DNA barcoding;ITS2(internal transcribed spacer 2 sequence);psbAtrnH(psbAtrnH intergenic spacer sequence)

    基金項目 陜西省科技統(tǒng)籌創(chuàng)新工程項目(2016KTTSSF 01-01-02);陜西省咸陽市重大科技專項(2014K 01-17)。

    作者簡介 趙新禮(1961—),男,陜西扶風(fēng)人,主任藥師,從事中藥資源與栽培及分子生物學(xué)研究。

    收稿日期 2021-01-27

    倒卵葉五加與著名的中藥材刺五加具有相似的化學(xué)成分和藥理作用[1-2],倒卵葉五加與短柄五加兩者分布區(qū)域重疊,植物形態(tài)高度相似,在其產(chǎn)區(qū)內(nèi)均作為倒卵葉五加藥材使用,究竟兩者為同一種植物還是2種植物長期以來存在爭議[3-5]。自20世紀80年代倒卵葉五加被開發(fā)以來,由于無節(jié)制的采集使資源量銳減,為保證用藥的安全性和有效性,以及其野生撫育與栽培物種的準確性,有必要對該中藥材的基原進行鑒定研究。

    傳統(tǒng)的中藥材基原鑒定即依據(jù)其植物形態(tài)學(xué)和分類學(xué)原理,根據(jù)形態(tài)鑒別特征,確定到物種的生物學(xué)名,這種鑒定方式受到鑒定人員知識和經(jīng)驗等因素的影響,具有一定的局限性。DNA條形碼技術(shù)是利用基因組中一段或幾段公認的、標準的、有足夠變異的、易擴增且相對較短的DNA片段,以自身在物種種內(nèi)的特異性和種間的多樣性而創(chuàng)建的分子鑒定技術(shù),可以對物種實現(xiàn)準確、快速、簡便的鑒定,不受發(fā)育階段、供試部位、環(huán)境條件等的限制,彌補了傳統(tǒng)鑒定方法在這方面的不足。

    陳士林等[6]通過對大量中藥材進行DNA條形碼分子鑒定研究,建立以ITS2為核心、psbA-trnH為輔的植物類藥材DNA條形碼鑒定體系,并將該方法納入2015年版《中國藥典》。ITS2序列是去除5.8S和28S rRNA基因序列的核糖體DNA序列間隔區(qū),psbA-trnH序列是位于葉綠體編碼光合系統(tǒng)D1蛋白的psbA基因和編碼tRNA組氨酸的trnH基因之間的一段非編碼區(qū),兩者均具有進化速率快、引物通用性好、片段足夠短、易擴增與測序的優(yōu)點,已廣泛用于人參屬[7]、忍冬屬[8]、枸杞屬[9]、益母草[10]、覆盆子[11]、商陸[12]、活血丹[13]、金叉石斛[14]等藥材基原物種鑒定,五加皮[15]、竹節(jié)參[16]、當(dāng)歸[17]等中藥材真?zhèn)舞b別,人參[18]、西洋參[18]、蒼術(shù)[19]等中藥材種子種苗的鑒定。

    該研究采用植物形態(tài)鑒定與DNA條形碼分子鑒定技術(shù)相結(jié)合,利用通用引物ITS2和psbA-trnH序列對五加屬短柄五加等4種植物的50份標本進行鑒定分析,結(jié)合GenBank數(shù)據(jù)庫中的近緣物種序列構(gòu)建系統(tǒng)進化樹,弄清倒卵葉五加和短柄五加間的物種關(guān)系,為倒卵葉五加藥材DNA條形碼鑒定方法的建立提供依據(jù)。

    1 材料與方法

    1.1 材料

    1.1.1 試材。

    1.1.1.1 館藏標本。西北農(nóng)林科技大學(xué)和陜西中藥研究所館藏的倒卵葉五加和短柄五加標本。

    1.1.1.2 現(xiàn)地調(diào)查。對陜西、甘肅、寧夏的16個區(qū)縣進行五加資源野外調(diào)查、采集標本與分析樣品(表1)。標本經(jīng)西北農(nóng)林科技大學(xué)吳振海高級實驗師鑒定,現(xiàn)存放在陜西中藥研究所標本室。

    1.1.1.3 樣品采集。每個區(qū)縣選取2~4個五加居群,在每個居群內(nèi)選取相距2 m以上的10個植株,每株采集嫩葉一片放入自封袋內(nèi),加入硅膠快速干燥。待葉片徹底干燥后存放于-40 ℃的冰箱中備用。

    1.1.1.4 五加屬等物種的ITS2和psbA-trnH序列。五加屬等物種的ITS2和psbA-trnH序列均來源于GenBank數(shù)據(jù)庫(表2)。

    1.1.2 試劑。

    DNA Marker、植物基因組DNA提取試劑盒(北京天根生化);Tris、EDTA、冰乙酸、瓊脂糖等(西安海默生物);DNA凝膠回收試劑盒(北京擎科新業(yè)生物);Prime STAR HS DNA Polymerase(Takara,美國)。

    1.1.3 儀器。

    高速冷凍離心機(eppendorf 5810R,德國);PCR擴增儀(Life,美國);制冰機(IMS-50,西安儀創(chuàng));電泳儀(JUNYI-300,北京君意);凝膠成像儀(SAGECREATION,北京賽智);基因測序儀器(ABI 3730,美國)。

    1.2 方法

    1.2.1 植物標本鑒定。依照植物的形態(tài)特征,參考文獻資料,按傳統(tǒng)方法進行鑒定[4]。

    1.2.2 DNA提取。用干凈的剪刀將試樣剪碎,混勻,稱取適量,用植物基因組DNA提取試劑盒依照說明書提取,提取液 -20 ℃ 保存。

    1.2.3 引物。參見《中藥DNA條形碼分子鑒定》[20],由北京擎科新業(yè)生物技術(shù)有限公司合成(表3)。

    1.2.4 PCR擴增及測序。

    1.2.4.1 擴增體系。DNA模板1 μL,5X Prime STAR Buffer 5 μL,dNTP Mixture 2 μL,引物1與引物2各1 μL,Prime STAR HS DNA Polymerase ?0.25 μL,ddH2O 14.75 μL補足至 25 μL。

    1.2.4.2 擴增程序。94 ℃預(yù)變性 5 min;94 ℃變性30 s,60 ℃ 退火30 s,72 ℃ 延伸 35 s,33 個循環(huán);72 ℃ 5 min,16 ℃停止。

    1.2.4.3

    PCR產(chǎn)物電泳及DNA純化回收。PCR產(chǎn)物經(jīng)15%的瓊脂糖凝膠電泳后將目的片段割膠,并使用DNA凝膠回收試劑盒按說明書進行DNA片段的純化回收。

    1.2.4.4

    DNA序列的測序。將純化的產(chǎn)物送北京擎科新業(yè)生物技術(shù)有限公司進行雙向測序。

    1.3 數(shù)據(jù)處理

    用Chromas軟件對測序所得到的峰圖進行人工校對;用CodonCode Aligner V6.0.2軟件進行拼接。ITS2序列基于隱馬爾可夫模型的HMMer注釋方法采用ITS2 database網(wǎng)站(http://its2.bioapps.biozentrum.uniwuerzburg.de/)進行序列注釋,去除兩端的5.8S rRNA和28S rRNA區(qū)段,獲得完整的ITS2序列。根據(jù)GenBank數(shù)據(jù)庫中同屬物種psbA-trnH的注釋,去除序列兩端的psbA和trnH基因,獲得完整的psbA-trnH基因間隔區(qū)序列。用MEGA 6.06軟件進行序列對比,分析變異位點,計算GC含量,用Kimura 2-Parameter(K2P)法將所有數(shù)據(jù)計算種內(nèi)和種間遺傳距離,用NJ鄰接法(Neighbour 2 Joining Method,NJ)將比對過的序列分別構(gòu)建成五加屬植物的ITS2、pbsA-trnH序列的系統(tǒng)發(fā)育樹,各個分支的支持率使用bootstrap重復(fù)1 000次進行檢驗。最終將獲得序列在NCBI(https://blast.ncbi.nlm.nih.gov/Blast.cgi)中進行樣品序列的BLAST(basic local alignment search tool)用最高相似搜索法進行確認比對分析。

    2 結(jié)果與分析

    2.1 植物形態(tài)鑒定

    對館藏的倒卵葉五加和短柄五加標本進行了仔細觀測比對,未發(fā)現(xiàn)兩者的顯著區(qū)別。

    2018年7—8月到陜甘寧地區(qū)的16個區(qū)縣,對50個五加居群進行野外調(diào)查,采集標本與分析樣品。采集的標本經(jīng)鑒定,其中倒卵葉五加和短柄五加共42份,均鑒定為短柄五加;藤五加2份;蜀五加4份;糙葉五加2份。調(diào)查發(fā)現(xiàn)倒卵葉五加藥材的混淆品有藤五加和蜀五加(甘肅)、糙葉五加(陜西隴縣)。

    倒卵葉五加藥材主產(chǎn)地的人們對其原植物短柄五加的叫法不一,如在陜西的旬邑、黃陵、黃龍、寶塔區(qū)、志丹、隴縣,甘肅的華池、合水、正寧、靈臺等地稱倒卵葉五加或倒刨牛;在宜君、甘泉、涇源、清水、崆峒區(qū)等地稱短柄五加或甘肅刺五加,其植物特征為:直立灌木,莖圓柱形,多分枝,株高0.5~3.0 m。在灌叢的頂部常生出向上的長枝與側(cè)生出短枝,下部老枝及根莖發(fā)出的枝條全為長枝;幼枝無毛,葉柄的基部常生1~2個刺,細長,下彎,基部不膨大。掌狀復(fù)葉,小葉通常5枚;在長枝上葉互生,在短枝上葉簇生;幼枝的上部總?cè)~柄無或極短,無毛,無刺;從枝頂15~45 cm往下,總?cè)~柄漸長,達2~12 cm,具縱棱。從灌叢枝條上端側(cè)生的短枝上,葉近于無柄,無毛,無刺,小葉亦較小,常數(shù)片簇生。小葉片薄紙質(zhì),倒卵形、倒卵狀長圓形,稀菱形,葉尖圓形或短尖,基部狹尖,楔形,兩面均無毛,邊緣近全緣或上部邊緣具數(shù)對鈍形齒牙;側(cè)脈3~5對。傘形花序單生或數(shù)個形成圓錐花序頂生于枝端,徑1~4 cm,有花多數(shù)??偦üiL1~8 cm,具縱棱,無毛,無刺;花梗細長,長0.6~2.0 cm,無毛。由于受生長環(huán)境的影響,植株間在形態(tài)上有少許差異,如短枝上或光照不好的情況下傘形花序較小、總花梗較短,長枝上或光照較好時傘形花序較大,花梗也較長;花淺黃綠色。苞片卵形,紫色,長約1 mm,先端叢生銹毛,邊緣疏生纖毛;花萼無毛或有短柔毛,邊緣具5個三角形小裂齒?;ò?枚,三角狀卵形,先端尖,無毛,長約2 mm,開花時反曲。雄蕊5個,花絲長約2 mm?;ㄖ?個,全部合生成柱狀。子房5室。果實橢圓狀球形,未成熟時綠色、綠黑色、黑色,成熟后黑紫色,直徑5~8 mm,風(fēng)干時有5棱,宿存花柱長1.5~2.0 mm。花期7—8月。

    2.2 PCR擴增結(jié)果

    所有五加樣本ITS2序列與psbA-trnH序列PCR擴增均在序列長度500~750 bp出現(xiàn)一條清晰的DNA條帶(圖 1)。

    2.3 序列測定結(jié)果

    2.3.1

    ITS2序列。ITS2引物的擴增效率和測序成功率均為100%。短柄五加等樣本經(jīng)ITS2 ?F/R引物擴增后得到一段長455~498 bp的序列,去掉兩端的5.8S和28S區(qū)段后,獲得的ITS2序列長度為229~230 bp,GC含量為62.6%~63.0%。經(jīng)兩頭等齊對比,短柄五加在序列18位有一個C的缺失;糙葉五加與短柄五加相比在18位多一個C;藤五加與蜀五加的ITS2序列完全相同,與短柄五加相比在18位多一個C,在24位發(fā)生T/A的轉(zhuǎn)換,在25位發(fā)生G/T的顛換。

    短柄五加ITS2序列:

    1 CGCATCGCGT CGCCCCC-AA CCCTGCACTC CCTCATGGGA GTCATGACTG

    51 AGGGGCGGAT ACTGGCCTCC CGTGTCTCAC CGCGCGGTTG GCCCAAATGT

    101 GAGTCCTTGG CGACGGACGT CACGACAAGT GGTGGTTGTA AAAAGCCCTC

    151 TTCTCCTGTC GTGCGGTGGC CCGTCGCCAG CAAAAGCTCA TGTGACCCTG

    201 TTGTGCCGTC CTCGACGAGC ACTCCGACCG

    2.3.2

    psbA-trnH序列。psbA-trnH引物的擴增效率與測序成功率均為100%。短柄五加等樣本經(jīng)psbA-F/trnH-R引物擴增得到一段長487~522 bp的序列,去掉兩端的psbA和trnH基因后,獲得的psbA-trnH序列長度為392~399 bp,GC含量為28.3%~28.6%。psbA-trnH序列分別在短柄五加、蜀五加和糙葉五加種內(nèi)變異位點數(shù)差異為0,在藤五加種內(nèi)有一個變異位點的差異。短柄五加、蜀五加及甘肅清水產(chǎn)的藤五加三者psbA-trnH序列完全相同,甘肅清水產(chǎn)藤五加以3葉為主,少5葉;陜西隴縣產(chǎn)的藤五加全為3葉,與短柄五加相比在290位有一個G/C轉(zhuǎn)換。糙葉五加與短柄五加相比分別在149、294、397位有T/G、A/G、G/A的顛換,在132~133、393位發(fā)生AA/TT、C/G的轉(zhuǎn)換,在154~155位發(fā)生AT的缺失,在330~338位有一個CATAAGAAA序列片段的插入。

    短柄五加psbA-trnH 序列:

    1 GACCCGGTCT TAGTGTATAC GAGTTTTTTT GAAATAAAAA GGAGCAATAA

    51 CCTCTTTCTT GTTCTATCAA GAGGGCGGTA TTGCTCCTTT TTTTATTTAG

    101 TAGTATTTGT ACTTACCTAG TTTCTTTAAA TTTATTTAAA GAAACGTCGT

    151 TTTATTTTTA TTGGTTGGTT CATGAGTGAG TATCATGCTT TCGTTCTGTA

    201 TTAATCATTA ATCTAGAATT TAGCTACTTC TTCCCAATCT TTTGGGAAGT

    251 ATTTTTTAAA AAGAAAAGAT AAAAATGTTG GACTTTTTAC TTAGTTAATA

    301 CTTAATATTT TTATCTCGAA AATAAGAAAG AAATAATGAT GAATGGTAAA

    351 AGTTTAATCT TTTGAAACGT AAGGAAAAGA ATAGTAGAGG GG

    2.4 五加屬植物ITS2序列的差異與遺傳距離分析

    五加屬植物的ITS2序列長度為229~231 bp,GC含量61.1%~641%。當(dāng)空位(gap)作為缺失處理時,序列兩端對齊比對,屬內(nèi)保守位點210個,變異位點21個,簡約信息位點8個,自裔位點13個(表4)。短柄五加、糙葉五加、藤五加和蜀五加的ITS2序列在種內(nèi)完全相同,種內(nèi)遺傳距離均為0。屬內(nèi)種間遺傳距離為0.000~0.055,平均為0.020。短柄五加和糙葉五加與紅毛五加間的遺傳距離為0;藤五加與蜀五加間的遺傳距離為0,刺五加與無梗五加間的遺傳距離為0;烏蘞莓五加與狹葉五加間的遺傳距離為0??梢奍TS2序列在同物種之間高度保守,同屬不同種間、種內(nèi)變種之間遺傳距離差異有時不明顯,因此,ITS2序列只能作為鑒別五加屬物種的依據(jù)之一。

    2.5 五加屬植物psbA-trnH序列的差異與遺傳距離分析

    五加屬植物的psbA-trnH序列長度為392~399 bp,GC含量27.4%~28.8%。當(dāng)空位(gap)作為缺失處理時,序列兩端對齊比對(401 bp),屬內(nèi)保守位點387個,變異位點14個,簡約信息位點11個,自裔位點3個(表5)。糙葉五加、二歧五加、無梗五加和 E.japonicus 在154~155位有一個AT的缺失,在330~338位分別有一個CATAAGAAA的插入片段。短柄五加、蜀五加、糙葉五加3個種種內(nèi)遺傳距離均為0;甘肅產(chǎn)藤五加與陜西產(chǎn)藤五加的種內(nèi)遺傳距離為0.003。屬內(nèi)種間遺傳距離為0.000~0.026,平均為0.013。短柄五加與蜀五加、紅毛五加和甘肅產(chǎn)藤五加間的遺傳距離為0;糙葉五加與 E.japonicus 間的遺傳距離為0;二歧五加與無梗五加、刺五加之間的遺傳距離為0??梢?,psbA-trnH序列在同物種之間高度保守,同屬不同種間遺傳距離差異有時不明顯,因此psbA-trnH序列也只能作為鑒別五加屬物種的依據(jù)之一。

    將ITS2和psbA-trnH序列的特征參數(shù)匯總(表6)。K2P遺傳距離計算顯示,psbA-trnH序列種間遺傳距離為0.013,小于ITS2序列的種間遺傳距離(0.020);psbA-trnH序列的變異率為3.50%,小于ITS2序列的9.09%,說明psbA-trnH較ITS2序列保守,2種候選序列種內(nèi)遺傳距離為0或0~0003,小于種間遺傳距離,存在明顯的barcoding gap。

    2.6 序列比對鑒定

    對測定所獲得的ITS2、psbA-trnH序列在NCBI(https://blast.ncbi.nlm.nih.gov/Blast.cgi)中利用BLAST進行最高相似性搜索,進行確認比對分析。

    2.6.1 倒卵葉五加和短柄五加的鑒定。共42個樣本,從植物形態(tài)上看完全相同,統(tǒng)一鑒定為短柄五加;所有樣品其測序結(jié)果ITS2序列、psbA-trnH序列完全相同。其ITS2序列與紅毛五加(KR080440.1、KR080438.1、MH711360.1)和藤五加(MH711366.1)的相似性為99.57%,而psbA-trnH序列與紅毛五加(GQ435397.1)的相似性為100%,因此認為短柄五加和倒卵葉五加是同一種植物,且與紅毛五加的親緣關(guān)系較近。

    2.6.2 藤五加和蜀五加的鑒定。2個種共6個樣為一組,其ITS2序列完全相同,與藤五加(MH703148.1)、蜀五加(MH710849.1、MH703119.1)的相似性為100%。

    在GenBank中無藤五加和蜀五加psbA-trnH序列的相關(guān)數(shù)據(jù)。該研究所測定的蜀五加4個樣和甘肅產(chǎn)藤五加1個樣共5個樣的 psbA-trnH 與紅毛五加(GQ435397.1)的相似性為100%;陜西隴縣產(chǎn)藤五加的 psbA-trnH序列與紅毛五加(GQ435397.1)的相似性為99.74%。

    2.6.3 糙葉五加的鑒定。2個樣本為一組,其ITS2序列與紅毛五加(KR080440.1、KR080438.1、MH711360.1)和藤五加(MH711366.1)的相似性為100%,與刺五加(KU724205.1)和無梗五加(MF095901.1)的相似性為99.57%;其psbA-trnH序列與刺五加(MF348735.1)、無梗五加(KT153019.1)的相似性為99.75%,與紅毛五加(GQ435397.1)的相似性僅為96.01%。

    在GenBank中未見糙葉五加 ITS2和psbA-trnH序列的相關(guān)數(shù)據(jù)。

    2.7 序列聚類分析

    基于 K2P雙參數(shù)模型構(gòu)建的系統(tǒng)進化樹能顯示物種間的親緣關(guān)系,幫助判斷物種的進化關(guān)系。五加屬可能起源于鵝掌柴屬( Schefflera )[3],故選其作為外類群構(gòu)建了系統(tǒng)進化樹。

    2.7.1 ITS2序列的聚類分析。五加屬植物的ITS2序列經(jīng)過全局比對構(gòu)建的NJ系統(tǒng)進化樹見圖2。五加屬植物聚為一大支,其中所有測定的短柄五加、糙葉五加和紅毛五加聚為一小支,藤五加和蜀五加聚為另一小分枝。

    2.7.2 psbA-trnH序列的聚類分析。psbA-trnH序列構(gòu)建的五加屬植物NJ系統(tǒng)進化樹見圖3。五加屬植物獨聚為一支,其中所有測定的短柄五加、藤五加和蜀五加與紅毛五加等聚為一小枝,糙葉五加與無梗五加、刺五加等聚為另一小枝。

    3 結(jié)論

    目前 DNA 條形碼研究主要集中在如何提高近緣類群物種分辨率和構(gòu)建區(qū)域條形碼數(shù)據(jù)庫2個方向。受物種特異性影響,同屬物種較低的種間差異性以及分化產(chǎn)生的新物種的出現(xiàn),都會增加物種分類鑒定的難度。

    建立良好的DNA條形碼鑒定體系需要幾個序列進行相互補充,或與其他分子標記結(jié)合,才能提高DNA條形碼技術(shù)鑒定的效率。

    盡管 ITS2 序列對大多數(shù)植物可以進行鑒定,但由于其序列較短,遺傳信息有限,在鑒定近緣種或變種上仍存在一些不能鑒定的情況,增加psbA-trnH序列作為物種鑒別的輔助序列,可以將短柄五加、糙葉五加、藤五加和蜀五加準確鑒定,為鑒定短柄五加及其同屬植物提供新的分子鑒定方法,為控制藥材質(zhì)量提供科學(xué)依據(jù)。

    通過對現(xiàn)地42個短柄五加和倒卵葉五加居群的觀測,認為總?cè)~柄的長短、傘形花序的大小、總花梗的長短是與標本采集枝條的著生部位有關(guān),是植物適應(yīng)環(huán)境條件的結(jié)果,不能作為劃分種的依據(jù)。通過對42 份植物樣本的ITS2和 psbA-trnH 序列測定分析,其種內(nèi)序列差異為0,與五加屬其他物種間序列差異明顯,可以認定兩者為同一個物種,命名為短柄五加,即倒卵葉五加藥材的基原植物為短柄五加。

    通過對糙葉五加、藤五加和蜀五加共8份植物樣品的ITS2和psbA-trnH序列分析,其結(jié)果與植物形態(tài)分類相吻合,也得到較好的分類結(jié)果。

    在Genbank中未見糙葉五加ITS2和psbA-trnH序列,藤五加和蜀五加psbA-trnH序列的相關(guān)數(shù)據(jù)。

    應(yīng)用 ITS2 序列和 psbA-trnH 序列的聚類分析所建的NJ系統(tǒng)發(fā)育樹存在一定差異,在親緣關(guān)系的分析上不夠精準,還需有較多的數(shù)據(jù)積累才行。

    參考文獻

    [1] 衛(wèi)平,宋慶生.陜西產(chǎn)四種五加化學(xué)成分的含量測定[J].西北藥學(xué)雜志,1987,2(3):13-15.

    [2] 簡毓峰,胡浩斌.短柄五加的化學(xué)成分與藥理活性研究進展[J].中藥材,2011,34(8):1302-1306.

    [3] 中國科學(xué)院中國植物志編輯委員會.中國植物志:第54卷[M].北京:科學(xué)出版社,1978:86-116.

    [4] Flora of China Editorial Committee.Flora of China:V13[M].Beijing:Science Press & St Lious:Missouri Botanical Garden Press,2011:466-472.

    [5] 向其柏.五加科植物的新分類群及某些修訂[J].南京林業(yè)大學(xué)學(xué)報(自然科學(xué)版),1985,9(2):15-28.

    [6] 陳士林,姚輝,韓建萍,等.中藥材 DNA 條形碼分子鑒定指導(dǎo)原則[J].中國中藥雜志,2013,38(2):141-148.

    [7] 孫濤,孔德英,滕少娜,等.基于ITS2 條形碼的人參屬物種鑒定研究[J].湖北農(nóng)業(yè)科學(xué),2016,55(19):5072-5074.

    [8] 潘綠昌,李園園,劉紅昌,等.基于ITS2和psbA-trnH序列的忍冬屬5種中藥材DNA條形碼鑒定研究[J/OL].基因組學(xué)與應(yīng)用生物學(xué),2018-11-15[2020-10-25].http://kns.cnki.net/kcms/detail/45.1369.Q.20181115.0924.002.html.

    [9] 萬如,王亞軍,安巍,等.基于 psbA-trnH 序列條形碼鑒定21 份枸杞屬植物[J].江蘇農(nóng)業(yè)科學(xué),2019,47(1):56-59.

    [10] 李笑,譚朝陽,徐德宏,等.益母草植物DNA條形碼的鑒定研究[J].中國藥師,2020,23(3):442-447.

    [11] 方潔,呂群丹,陳正道,等.市售覆盆子藥材DNA條形碼鑒定研究[J].中國現(xiàn)代應(yīng)用藥學(xué),2020,37(4):437-442.

    [12] 呂瑞華,馮昭,馬添翼,等.陜西關(guān)中野生商陸資源的ITS2和psbA-trnH 條形碼序列研究[J].藥學(xué)學(xué)報,2020,55(8):1951-1956.

    [13] 鄭夢迪,張彥,張寒,等.活血丹及其同屬藥用植物的分子鑒定[J].時珍國醫(yī)國藥,2019,30(6):1397-1398.

    [14] 李園園,蔡莉,楊繼勇,等.基于ITS2和psbA-trnH序列的金釵石斛DNA條形碼鑒定研究[J].基因組學(xué)與應(yīng)用生物學(xué),2018,37(8):3516-3524.

    [15] 曾旭,李秋實,王鈐,等.基于 ITS2 條形碼的五加皮基原植物及其易混品的鑒定[J].中國現(xiàn)代中藥,2012,14(2):15-18.

    [16] 陳鏡安,楊璐,李榮釗,等.基于ITS2的竹節(jié)參及其近緣物種和混偽品鑒定評估[J].中草藥,2018,49(15):3672-3680.

    [17] 劉新星,歐巧明,石有太,等.基于ITS2 序列鑒別道地藥材岷縣當(dāng)歸及其混偽品[J].中草藥,2018,49(20):4877-4883.

    [18] 林鳳越,曹輝,任歡歡,等.基于ITS2序列的人參及同屬易混品西洋參種子的分子鑒定[J].中草藥,2019,50(9):2188-2193.

    [19] 劉金欣,李耿,陳彩霞,等.基于 ITS2 序列的中藥材蒼術(shù)種苗 DNA條形碼鑒定[J].中國實驗方劑學(xué)雜志,2018,24(2):34-38.

    [20] 陳士林.中藥DNA條形碼分子鑒定[M].北京:人民衛(wèi)生出版社,2012:494-495.

    久久精品国产99精品国产亚洲性色| 日韩欧美国产一区二区入口| 亚洲,欧美精品.| 国产亚洲精品久久久久5区| 欧美色视频一区免费| 美女扒开内裤让男人捅视频| 老司机在亚洲福利影院| 亚洲片人在线观看| 最近最新中文字幕大全免费视频| 欧美成人免费av一区二区三区| 亚洲av中文字字幕乱码综合| 欧美一区二区精品小视频在线| 国产高清激情床上av| 久久人妻福利社区极品人妻图片| 国产单亲对白刺激| 久久精品亚洲精品国产色婷小说| 成人三级黄色视频| 极品教师在线免费播放| 99热6这里只有精品| 脱女人内裤的视频| 999精品在线视频| 亚洲激情在线av| 69av精品久久久久久| 性欧美人与动物交配| 人人妻人人看人人澡| 每晚都被弄得嗷嗷叫到高潮| 久久久久久国产a免费观看| 日本一二三区视频观看| 两个人看的免费小视频| 成人特级黄色片久久久久久久| 亚洲美女视频黄频| 最近最新免费中文字幕在线| 91大片在线观看| av免费在线观看网站| 99re在线观看精品视频| 母亲3免费完整高清在线观看| 精品欧美一区二区三区在线| 99国产精品一区二区三区| 国产亚洲精品第一综合不卡| 一夜夜www| 成年免费大片在线观看| 亚洲va日本ⅴa欧美va伊人久久| 免费看a级黄色片| 亚洲无线在线观看| 日本一区二区免费在线视频| 免费一级毛片在线播放高清视频| av中文乱码字幕在线| 久久久国产精品麻豆| 亚洲一卡2卡3卡4卡5卡精品中文| 人人妻人人看人人澡| 此物有八面人人有两片| 欧美日韩中文字幕国产精品一区二区三区| 三级毛片av免费| 日日摸夜夜添夜夜添小说| 我的老师免费观看完整版| 十八禁网站免费在线| 亚洲色图 男人天堂 中文字幕| 两个人看的免费小视频| 成人国产综合亚洲| 99re在线观看精品视频| 一本久久中文字幕| 黄色毛片三级朝国网站| 午夜福利18| a级毛片a级免费在线| 国产成人一区二区三区免费视频网站| 国产免费男女视频| 五月伊人婷婷丁香| 亚洲黑人精品在线| 色尼玛亚洲综合影院| 国产精品一区二区精品视频观看| 黄色女人牲交| 亚洲欧美日韩高清专用| 99久久99久久久精品蜜桃| 波多野结衣高清无吗| 亚洲精品国产精品久久久不卡| 亚洲欧美日韩高清专用| 黑人操中国人逼视频| 国产男靠女视频免费网站| 久久中文字幕人妻熟女| 午夜成年电影在线免费观看| 日本黄大片高清| www.999成人在线观看| tocl精华| 亚洲国产日韩欧美精品在线观看 | 啦啦啦韩国在线观看视频| 欧美日韩黄片免| 天堂动漫精品| 可以在线观看的亚洲视频| 国产激情偷乱视频一区二区| 国产成人系列免费观看| 免费搜索国产男女视频| 淫妇啪啪啪对白视频| 久久久久久免费高清国产稀缺| 男女床上黄色一级片免费看| 国产视频内射| 亚洲欧美精品综合一区二区三区| 国产蜜桃级精品一区二区三区| 在线观看午夜福利视频| 亚洲欧美日韩高清在线视频| 人人妻,人人澡人人爽秒播| 天天添夜夜摸| 日本在线视频免费播放| 久久久久久人人人人人| 亚洲性夜色夜夜综合| 国产97色在线日韩免费| 岛国在线观看网站| 中文字幕人妻丝袜一区二区| 国产97色在线日韩免费| 小说图片视频综合网站| av有码第一页| 91麻豆av在线| 欧美日韩黄片免| 国产精品日韩av在线免费观看| 99久久国产精品久久久| 国产免费男女视频| 深夜精品福利| 成年免费大片在线观看| 夜夜看夜夜爽夜夜摸| www.999成人在线观看| 99久久国产精品久久久| 亚洲自偷自拍图片 自拍| 一a级毛片在线观看| 中亚洲国语对白在线视频| xxx96com| 久久国产乱子伦精品免费另类| 日日干狠狠操夜夜爽| 在线播放国产精品三级| 亚洲全国av大片| 最近视频中文字幕2019在线8| 91麻豆av在线| 悠悠久久av| 欧美性猛交╳xxx乱大交人| 亚洲第一欧美日韩一区二区三区| 欧美日韩精品网址| 19禁男女啪啪无遮挡网站| 亚洲成a人片在线一区二区| 日韩精品中文字幕看吧| 亚洲av成人av| 免费在线观看成人毛片| 国产精品1区2区在线观看.| 日韩欧美三级三区| 97超级碰碰碰精品色视频在线观看| 亚洲国产精品合色在线| 亚洲精品美女久久久久99蜜臀| 18禁黄网站禁片午夜丰满| 老司机深夜福利视频在线观看| avwww免费| 亚洲精品久久国产高清桃花| 色综合婷婷激情| 亚洲成人中文字幕在线播放| 给我免费播放毛片高清在线观看| 国产亚洲精品一区二区www| 中文字幕人成人乱码亚洲影| 一本久久中文字幕| 国产免费男女视频| 亚洲av成人av| 99热这里只有精品一区 | 国产精品影院久久| 午夜福利在线观看吧| 欧美一级a爱片免费观看看 | 亚洲国产精品久久男人天堂| 波多野结衣高清无吗| 大型av网站在线播放| 欧美大码av| 午夜福利高清视频| 夜夜爽天天搞| 国产精品电影一区二区三区| 男女做爰动态图高潮gif福利片| 给我免费播放毛片高清在线观看| 老司机深夜福利视频在线观看| 亚洲激情在线av| 国产成人一区二区三区免费视频网站| 国产精品国产高清国产av| 美女黄网站色视频| 又爽又黄无遮挡网站| 三级男女做爰猛烈吃奶摸视频| 九色成人免费人妻av| 欧美黄色淫秽网站| 亚洲中文字幕一区二区三区有码在线看 | 免费看美女性在线毛片视频| 成人18禁高潮啪啪吃奶动态图| 精品电影一区二区在线| 一个人免费在线观看电影 | 亚洲精品在线美女| 日本a在线网址| 亚洲天堂国产精品一区在线| 欧美av亚洲av综合av国产av| 久9热在线精品视频| 国产熟女xx| 婷婷六月久久综合丁香| 国产日本99.免费观看| 欧美精品亚洲一区二区| 亚洲人成77777在线视频| 欧美另类亚洲清纯唯美| 午夜免费成人在线视频| 久久九九热精品免费| 亚洲专区字幕在线| 国产黄a三级三级三级人| 免费看日本二区| 全区人妻精品视频| 日日干狠狠操夜夜爽| 国产成人av激情在线播放| 久久欧美精品欧美久久欧美| 青草久久国产| 色老头精品视频在线观看| 亚洲精品美女久久av网站| 欧美3d第一页| 亚洲全国av大片| 国产精品 国内视频| 欧美久久黑人一区二区| 99re在线观看精品视频| 国产日本99.免费观看| 一边摸一边做爽爽视频免费| 国产成人精品久久二区二区91| 国内精品久久久久精免费| 亚洲国产精品999在线| 一本久久中文字幕| 亚洲av五月六月丁香网| 国产精品精品国产色婷婷| 精品久久久久久久久久久久久| www.自偷自拍.com| 嫁个100分男人电影在线观看| 欧美日本视频| 免费电影在线观看免费观看| 特大巨黑吊av在线直播| 欧美精品亚洲一区二区| 国产亚洲精品久久久久5区| 50天的宝宝边吃奶边哭怎么回事| 一级a爱片免费观看的视频| 久久香蕉激情| 全区人妻精品视频| 性色av乱码一区二区三区2| 50天的宝宝边吃奶边哭怎么回事| 最好的美女福利视频网| 精品熟女少妇八av免费久了| 女同久久另类99精品国产91| 99国产精品一区二区蜜桃av| 日日爽夜夜爽网站| 亚洲av美国av| 欧洲精品卡2卡3卡4卡5卡区| 亚洲成人精品中文字幕电影| 操出白浆在线播放| 久久精品91蜜桃| 国产精品爽爽va在线观看网站| 国产精品永久免费网站| 天天躁狠狠躁夜夜躁狠狠躁| 日本一区二区免费在线视频| 999精品在线视频| 亚洲国产精品成人综合色| 中文字幕av在线有码专区| 久久99热这里只有精品18| 免费在线观看亚洲国产| 少妇粗大呻吟视频| 成年人黄色毛片网站| 欧美日韩瑟瑟在线播放| 淫秽高清视频在线观看| 国产激情欧美一区二区| 亚洲,欧美精品.| 欧美成人免费av一区二区三区| 999久久久国产精品视频| 亚洲五月天丁香| 中文资源天堂在线| 久久香蕉国产精品| 国产精品一区二区免费欧美| 99久久精品热视频| 首页视频小说图片口味搜索| 亚洲 国产 在线| 桃红色精品国产亚洲av| 日韩精品中文字幕看吧| а√天堂www在线а√下载| 亚洲精品av麻豆狂野| 亚洲激情在线av| 色综合欧美亚洲国产小说| 91字幕亚洲| a在线观看视频网站| 国内精品一区二区在线观看| 国产精品爽爽va在线观看网站| xxx96com| 久久婷婷人人爽人人干人人爱| 国产97色在线日韩免费| 国产av在哪里看| 午夜精品在线福利| 日韩中文字幕欧美一区二区| 成人午夜高清在线视频| 亚洲最大成人中文| 精品国产乱码久久久久久男人| 国产探花在线观看一区二区| 国产高清视频在线观看网站| 正在播放国产对白刺激| 国产不卡一卡二| 老司机午夜十八禁免费视频| 老司机靠b影院| 麻豆av在线久日| 又黄又粗又硬又大视频| videosex国产| 亚洲在线自拍视频| 在线免费观看的www视频| 欧美日韩亚洲国产一区二区在线观看| 97人妻精品一区二区三区麻豆| 操出白浆在线播放| 精品电影一区二区在线| 操出白浆在线播放| 久久久久久免费高清国产稀缺| 精品久久久久久久末码| 久热爱精品视频在线9| 精品久久久久久久毛片微露脸| 成人三级做爰电影| 国产精品亚洲一级av第二区| 亚洲精品色激情综合| 久久这里只有精品19| 国产精品免费视频内射| 久久久久久久精品吃奶| 亚洲午夜精品一区,二区,三区| 精品国产乱码久久久久久男人| 精品久久久久久久人妻蜜臀av| 亚洲精品美女久久久久99蜜臀| 一本精品99久久精品77| 欧美午夜高清在线| 国产又黄又爽又无遮挡在线| 亚洲成人国产一区在线观看| 99久久综合精品五月天人人| 久久人人精品亚洲av| 中文字幕精品亚洲无线码一区| 99久久久亚洲精品蜜臀av| 俺也久久电影网| 成年版毛片免费区| 人成视频在线观看免费观看| 欧美激情久久久久久爽电影| 两个人视频免费观看高清| 欧美乱码精品一区二区三区| 国产精品香港三级国产av潘金莲| 日韩大码丰满熟妇| 午夜久久久久精精品| 国产午夜福利久久久久久| 叶爱在线成人免费视频播放| www.熟女人妻精品国产| 免费看a级黄色片| 免费看十八禁软件| 中文亚洲av片在线观看爽| 色播亚洲综合网| 成人高潮视频无遮挡免费网站| 免费观看人在逋| 香蕉久久夜色| 欧美日韩亚洲综合一区二区三区_| 真人做人爱边吃奶动态| 亚洲人成77777在线视频| 麻豆一二三区av精品| 亚洲精品美女久久av网站| 午夜精品一区二区三区免费看| 日韩欧美三级三区| 啦啦啦观看免费观看视频高清| 亚洲av美国av| 啦啦啦观看免费观看视频高清| 欧美黑人欧美精品刺激| 国产成人av激情在线播放| 美女高潮喷水抽搐中文字幕| 亚洲av成人精品一区久久| 99热只有精品国产| 天堂影院成人在线观看| 国产91精品成人一区二区三区| 成人18禁在线播放| 国产av一区二区精品久久| 午夜成年电影在线免费观看| 亚洲激情在线av| 正在播放国产对白刺激| 99精品欧美一区二区三区四区| 制服人妻中文乱码| 岛国视频午夜一区免费看| 亚洲天堂国产精品一区在线| 亚洲18禁久久av| 国产精品一区二区免费欧美| 亚洲熟妇熟女久久| 欧美高清成人免费视频www| 精品久久久久久久毛片微露脸| 色综合婷婷激情| 嫩草影院精品99| 黑人操中国人逼视频| 叶爱在线成人免费视频播放| 成人高潮视频无遮挡免费网站| 国产高清有码在线观看视频 | 国产一区二区激情短视频| 日日夜夜操网爽| 国产欧美日韩精品亚洲av| 91字幕亚洲| 午夜激情福利司机影院| 村上凉子中文字幕在线| 亚洲av电影不卡..在线观看| 婷婷六月久久综合丁香| xxx96com| 国产三级黄色录像| 香蕉久久夜色| 听说在线观看完整版免费高清| 亚洲国产日韩欧美精品在线观看 | 久久久久久国产a免费观看| 波多野结衣巨乳人妻| a级毛片a级免费在线| 91在线观看av| 亚洲欧美一区二区三区黑人| 曰老女人黄片| 人成视频在线观看免费观看| 国产精品亚洲一级av第二区| 久久久国产成人免费| 精品电影一区二区在线| 欧美丝袜亚洲另类 | 精品免费久久久久久久清纯| 午夜激情av网站| 又大又爽又粗| 人妻夜夜爽99麻豆av| 欧美人与性动交α欧美精品济南到| 老司机靠b影院| 国产成人av教育| 女人高潮潮喷娇喘18禁视频| 一卡2卡三卡四卡精品乱码亚洲| 色噜噜av男人的天堂激情| 成人18禁在线播放| 亚洲av成人av| 91在线观看av| 亚洲aⅴ乱码一区二区在线播放 | 在线看三级毛片| 国产三级在线视频| 国产精品久久电影中文字幕| 亚洲av美国av| 久久久久久久久中文| 丝袜美腿诱惑在线| 国产精品98久久久久久宅男小说| 成年人黄色毛片网站| 日本免费一区二区三区高清不卡| 在线看三级毛片| 18禁黄网站禁片免费观看直播| 欧洲精品卡2卡3卡4卡5卡区| 天堂动漫精品| 国产亚洲精品av在线| 亚洲精品中文字幕一二三四区| 午夜激情av网站| 国产精华一区二区三区| www.自偷自拍.com| 色综合欧美亚洲国产小说| 日韩欧美一区二区三区在线观看| 精品福利观看| 久久久久国内视频| 两个人的视频大全免费| 久久久精品国产亚洲av高清涩受| 少妇人妻一区二区三区视频| 国产精品 欧美亚洲| 在线免费观看的www视频| 国产av麻豆久久久久久久| 99国产综合亚洲精品| 在线十欧美十亚洲十日本专区| 国产精品久久久久久人妻精品电影| 亚洲成人免费电影在线观看| 一级作爱视频免费观看| 天天躁狠狠躁夜夜躁狠狠躁| 正在播放国产对白刺激| 久久久久久国产a免费观看| 成人永久免费在线观看视频| 欧美成人一区二区免费高清观看 | 国产精品九九99| 久久精品国产综合久久久| 18禁美女被吸乳视频| 国产黄片美女视频| 一本一本综合久久| 欧美日本亚洲视频在线播放| 丰满人妻熟妇乱又伦精品不卡| 精品乱码久久久久久99久播| 真人做人爱边吃奶动态| 亚洲一码二码三码区别大吗| 久久久久久免费高清国产稀缺| 精品一区二区三区四区五区乱码| av欧美777| 淫妇啪啪啪对白视频| 亚洲av电影在线进入| 免费在线观看视频国产中文字幕亚洲| 色精品久久人妻99蜜桃| 我的老师免费观看完整版| 国产视频内射| 免费av毛片视频| 国产探花在线观看一区二区| 欧美日韩乱码在线| 熟女电影av网| 欧美色欧美亚洲另类二区| 久久精品夜夜夜夜夜久久蜜豆 | 国内揄拍国产精品人妻在线| 久久久精品大字幕| 亚洲九九香蕉| 国产精品免费视频内射| 欧美日韩一级在线毛片| 亚洲国产精品sss在线观看| 正在播放国产对白刺激| 淫秽高清视频在线观看| 欧美精品亚洲一区二区| 亚洲成人久久爱视频| 亚洲成人中文字幕在线播放| 精品一区二区三区av网在线观看| 国产一区二区三区在线臀色熟女| 黄色成人免费大全| 黄色片一级片一级黄色片| 午夜福利欧美成人| 国产成人av激情在线播放| 久久中文字幕人妻熟女| 在线观看免费视频日本深夜| 美女 人体艺术 gogo| 国产真人三级小视频在线观看| 亚洲天堂国产精品一区在线| 五月玫瑰六月丁香| 精品久久久久久久久久免费视频| 18禁国产床啪视频网站| 制服人妻中文乱码| 久久久久久久久久黄片| 久久人妻av系列| 国产精品,欧美在线| 中文在线观看免费www的网站 | 中文在线观看免费www的网站 | 国产精品一区二区免费欧美| 亚洲精品色激情综合| 日韩大尺度精品在线看网址| 亚洲精品国产精品久久久不卡| 国产1区2区3区精品| 国产亚洲欧美在线一区二区| 麻豆一二三区av精品| 精品久久久久久久毛片微露脸| 欧美高清成人免费视频www| 国产一区二区三区在线臀色熟女| 国产午夜精品论理片| 成人午夜高清在线视频| 99国产精品一区二区三区| 18禁美女被吸乳视频| www日本在线高清视频| 99国产精品一区二区蜜桃av| 手机成人av网站| 国产精品 国内视频| 三级男女做爰猛烈吃奶摸视频| 又大又爽又粗| 国产一级毛片七仙女欲春2| 国产欧美日韩一区二区精品| 老司机靠b影院| 狂野欧美白嫩少妇大欣赏| 波多野结衣高清无吗| 亚洲一区高清亚洲精品| 欧美在线黄色| 亚洲avbb在线观看| 国产精品久久久人人做人人爽| 中文在线观看免费www的网站 | 中文字幕最新亚洲高清| 天堂动漫精品| 亚洲av成人一区二区三| 我要搜黄色片| 欧美绝顶高潮抽搐喷水| 一个人观看的视频www高清免费观看 | www日本在线高清视频| 久久久国产成人免费| 在线观看一区二区三区| 成人精品一区二区免费| 这个男人来自地球电影免费观看| 久久欧美精品欧美久久欧美| 中文在线观看免费www的网站 | 熟女电影av网| 久久精品国产99精品国产亚洲性色| 国产高清视频在线观看网站| 日韩精品中文字幕看吧| 国产精品美女特级片免费视频播放器 | 日韩欧美三级三区| 欧美性长视频在线观看| 精品熟女少妇八av免费久了| 99久久精品热视频| 欧美日韩国产亚洲二区| 亚洲色图 男人天堂 中文字幕| 99re在线观看精品视频| 精品久久久久久久人妻蜜臀av| 1024视频免费在线观看| 国产男靠女视频免费网站| 精品久久久久久久久久免费视频| 岛国在线免费视频观看| 午夜视频精品福利| 国产av一区在线观看免费| 舔av片在线| 久久香蕉精品热| 久久草成人影院| 亚洲精品国产精品久久久不卡| av福利片在线观看| 亚洲一区高清亚洲精品| 亚洲熟妇熟女久久| 亚洲专区国产一区二区| 亚洲一区高清亚洲精品| 床上黄色一级片| 在线观看66精品国产| 欧美3d第一页| 麻豆国产97在线/欧美 | svipshipincom国产片| 国产精品精品国产色婷婷| 妹子高潮喷水视频| 亚洲欧美精品综合久久99| 久久精品综合一区二区三区| 一二三四在线观看免费中文在| 亚洲免费av在线视频| 哪里可以看免费的av片| 三级男女做爰猛烈吃奶摸视频| 欧美3d第一页| 亚洲精品国产精品久久久不卡| 日本成人三级电影网站| 可以免费在线观看a视频的电影网站| 19禁男女啪啪无遮挡网站| 精品少妇一区二区三区视频日本电影| 国产视频内射| 床上黄色一级片| 欧美成人午夜精品| 精品日产1卡2卡| 亚洲狠狠婷婷综合久久图片| 身体一侧抽搐| www.熟女人妻精品国产| 男女之事视频高清在线观看| 成人国产一区最新在线观看| 正在播放国产对白刺激| 日日干狠狠操夜夜爽| 久久久久亚洲av毛片大全| 黑人欧美特级aaaaaa片|