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

    藥用植物全基因組測序研究及應(yīng)用進(jìn)展

    2025-05-09 00:00:00陳虞超施文韜李瑩郭生虎袁媛楊建國
    寧夏農(nóng)林科技 2025年1期

    摘要:我國藥用植物資源豐富,應(yīng)用歷史悠久,是中醫(yī)藥事業(yè)發(fā)展的物質(zhì)基礎(chǔ)和人民生命健康的重要保障。全基因組測序(WGS)是揭示生物體基因組特征的主要研究手段,在藥用植物研究中受到高度關(guān)注并被廣泛應(yīng)用,已成為促進(jìn)該領(lǐng)域研究創(chuàng)新發(fā)展的關(guān)鍵驅(qū)動力之一。介紹了WGS基本原理及其測序技術(shù)更新迭代過程,按照時間順序整理了2010—2023 年完成的WGS藥用植物種類,歸納出藥用植物WGS 發(fā)展歷程;概述了WGS 在揭示藥用植物進(jìn)化與馴化、有效成分合成的分子機制,以及其在藥用植物種質(zhì)資源鑒定、分子輔助育種等方面的應(yīng)用現(xiàn)狀,探討了藥用植物WGS 研究及應(yīng)用之中存在的主要問題,并對其未來研究方向進(jìn)行了展望,以期為深入開展藥用植物全基因組相關(guān)理論研究與實踐應(yīng)用提供一定借鑒。

    關(guān)鍵詞:藥用植物;全基因組測序;泛基因組;種質(zhì)資源鑒定;分子輔助育種

    中圖分類號:S154.5 文獻(xiàn)標(biāo)識碼:A 文章編號:1002-204X(2024)01-0031-14

    doi:10.3969/j.issn.1002-204x.2025.01.005

    全基因組測序(Whole genome sequencing, WGS)為前沿生物技術(shù),其從整體性與系統(tǒng)性角度出發(fā),描述生物體所有DNA堿基序列的詳細(xì)特征,具有高通量、高靈敏度、系統(tǒng)性強、穩(wěn)定性好、準(zhǔn)確度高等特點,在現(xiàn)代生命科學(xué)領(lǐng)域研究過程中發(fā)揮著至關(guān)重要的作用[1-2]。我國藥用植物資源豐富,約有11 000種,其中常用的有500 多種,這些藥用植物是保障人民群眾生命健康及中醫(yī)藥現(xiàn)代化與大健康產(chǎn)業(yè)發(fā)展的物質(zhì)基礎(chǔ)[3]。在過去相當(dāng)一段時間內(nèi),藥用植物研究多集中于種質(zhì)資源調(diào)查與鑒別、種苗繁育與人工栽培、藥理藥效分析等方面,有關(guān)其全基因組測序的研究報道甚少,造成藥用植物學(xué)和現(xiàn)代生命科學(xué)之間缺乏“互通橋梁”,嚴(yán)重阻礙了該領(lǐng)域研究的創(chuàng)新發(fā)展[4]。

    2010年,CHAN A P等[5]繪制了蓖麻(Ricinus com原munis)全基因組草圖,由此拉開了藥用植物全基因組學(xué)測序研究的序幕。此后,YAN L等[6]、KIM N H等[7]、YUAN Y等[8]、PU X D等[9]、TU L C等[10]、XIONG XY等[11]、MA B等[12]先后報道了鐵皮石斛(Dendrobiumoffcinale)、人參(Panax ginseng)、天麻(Gastrodiaela原ta)、金銀花(Lonicera japonica)、雷公藤(Tripterygiumwilfordii)、紅豆杉(Taxus chinensis)、紫丁香(Syringaoblata)等藥用植物的全基因組測序研究工作。BORNOWSKIN等[13]、ZHANG X H等[14]分別報道了迷迭香(Rosmarinus offcinalis)、蘿卜(Raphanus sativus)等藥用植物泛基因組測序研究進(jìn)展。目前,這些測序數(shù)據(jù)已應(yīng)用于藥用植物研究領(lǐng)域諸多方面。例如基于WGS數(shù)據(jù)分析,KIM N H等[7]揭示了人參的全基因組復(fù)制(Whole genome duplications,WGD)事件等進(jìn)化特征,TU L C等[10]解析了雷公藤中雷公藤甲素生物合成途徑,GUO C等[15]開展了薏苡(Coix aquatica)“紙殼”品種分子輔助育種。

    雖然藥用植物WGS 研究及應(yīng)用取得了長足發(fā)展,但尚未見對其現(xiàn)狀系統(tǒng)性的總結(jié)報道。鑒于此,本文詳細(xì)梳理了藥用植物WGS 研究及應(yīng)用進(jìn)展,探討了其中存在的主要問題,并對其未來研究方向進(jìn)行了展望,以期為深入開展藥用植物全基因組相關(guān)理論研究與實踐應(yīng)用提供一定借鑒。

    1 WGS原理及測序技術(shù)更新迭代過程

    1.1 WGS的基本原理及應(yīng)用發(fā)展

    WGS 是指對一個生物體攜帶的所有遺傳信息,即細(xì)胞核染色體、細(xì)胞質(zhì)線粒體及葉綠體(植物)上所有的DNA 堿基序列進(jìn)行測序,旨在揭示堿基序列中嵌入的遺傳組成、結(jié)構(gòu)、組織、功能、多樣性和相互作用,是當(dāng)前轉(zhuǎn)錄組、蛋白質(zhì)組和代謝組等組學(xué)的研究基礎(chǔ)[16-17]。WGS 基本流程:①將基因組降解為長度不一且相對較短的片段,將短片段進(jìn)行測序,測序后獲得原始數(shù)據(jù);②將這些原始數(shù)據(jù)經(jīng)質(zhì)控操作后獲得高質(zhì)量的reads,根據(jù)reads之間的重疊部分進(jìn)行片段拼接,生成較長的conting片段;③確定conting片段的順序及位置關(guān)系,進(jìn)一步將conting 片段拼接成長度更長的scaffold 片段,借助Hi-C(High-through chromosomeconformation capture)等基因組輔助組裝技術(shù),對上述獲得的scaffold 片段進(jìn)行染色體組別的劃分、組內(nèi)的排序及定向操作,獲得生物體染色體水平的全基因組序列信息[18-19]。

    2000 年,第一個高等植物擬南芥(Arabidopsisthaliana)WGS的完成,揭開了植物WGS研究的序幕[17]。隨后,毛果楊(Populus trichocarpa)、葡萄(Vitis vinifer原a)、番木瓜(Carica papaya)等植物WGS 陸續(xù)完成,為植物生長性狀遺傳解析、分子育種等研究奠定了必要的遺傳信息基礎(chǔ)。目前,WGS 已成為植物科學(xué)發(fā)展最快的研究領(lǐng)域之一[18-20]。

    1.2 WGS 測序技術(shù)更新迭代過程

    DNA 測序技術(shù)是推動WGS 研究發(fā)展的關(guān)鍵手段。從1977年至今,DNA測序技術(shù)經(jīng)歷了第一代到第三代的迭代發(fā)展過程。第一代DNA 測序技術(shù)以學(xué)者SANGER F等發(fā)明的雙脫氧鏈終止法(亦稱Sanger法)為平臺,其測序讀長可達(dá)1 000 bp,堿基的讀取準(zhǔn)確率高達(dá)99.9%,在20 世紀(jì)90 年代之前被廣泛應(yīng)用[21]。21 世紀(jì)初,隨著相關(guān)技術(shù)的發(fā)展,以Illumina 測序等為代表的第二代測序技術(shù)憑借高通量、低成本等特點,逐漸占據(jù)測序領(lǐng)域主導(dǎo)地位[22]。2010 年后,以單分子測序為基礎(chǔ)的第三代測序技術(shù)如tSMSTM(truesingle-molecule sequencing)、ONT(oxford nanopore technologies)、SMRT(single-molecule real-time sequencing)等日趨成熟,適度降低了準(zhǔn)確率(為85%~90%),但測序讀長增長了數(shù)十倍(可達(dá)105 bp),具有實時測序、讀長更長、測序速度更快、成本較低等優(yōu)點(圖1)。目前,研究者們?yōu)榱双@得更準(zhǔn)確的測序數(shù)據(jù),通常將第二、三代測序技術(shù)整合應(yīng)用,借助第二代測序技術(shù)彌補第三代測序技術(shù)讀取準(zhǔn)確率低等缺陷,同時充分發(fā)揮第三代測序技術(shù)讀長較長、測序速度快等優(yōu)勢[23]。隨著DNA測序技術(shù)的不斷發(fā)展和完善,全基因組測序愈發(fā)精準(zhǔn)快捷,且成本越來越低,推動更多種類的生物體全基因組測序得以實現(xiàn)。

    2 藥用植物WGS 研究現(xiàn)狀

    2.1 WGS藥用植物種類逐年增加

    WGS對推動藥用植物品種改良、資源保護(hù)等具有極其重要的意義,但與模式植物、經(jīng)濟林木等相比,藥用植物WGS 發(fā)展一度較為緩慢[4]。2010 年,CHAN AP等[5]通過第一代測序技術(shù)繪制了蓖麻基因組框架圖,成為首個報道的藥用植物WGS;2011—2015 年,受第一代測序技術(shù)效率低、成本高等限制,藥用植物WGS研究進(jìn)展較為緩慢,僅有BAKEL H V等[24]、WANG ZW等[25]、HE N J等[26]、WANG Y等[27]、UPADHYAY AK等[28]報道了大麻(Cannabis sativa)、亞麻(Linum usi原tatissimum)、川桑(Morus notabilis)、蓮(Nelumbo nu原cifera)、圣羅勒(Ocimum tenuiforum)等10 多種植物的基因組。此后,因第二代測序技術(shù)逐步擴大應(yīng)用,該領(lǐng)域研究報道數(shù)量逐年加大。2016 年,GUAN R等[29]報道了銀杏(Ginkgo biloba)等6種藥用植物的基因組,2017 年,ZHAO D Y等[30]、WANG X等[31]、LIN YL等[32]報道了喜樹(Camptotheca acuminata)、柑橘(Cit原rus medica)、龍眼(Dimocarpus longan)等10 種藥用植物的基因組,2018 年,YUAN Y等[8]報道了天麻等14種藥用植物的基因組,2019 年,ZHANG Y H等[33]、QIN S S等[34]報道了虎杖(Polygonum cuspidatum)、雞血藤(Spatholobus suberectus)等12 種藥用植物的基因組。隨著第三代測序技術(shù)日趨成熟,以及第二代、第三代測序技術(shù)聯(lián)合應(yīng)用的普及,2020 年開始,藥用植物WGS 研究報道數(shù)量快速增長(41種),2021年報道數(shù)量為88 種,2022 年報道數(shù)量增加到129 種,2023 年報道數(shù)量達(dá)到136 種(圖2、表1)。目前,已有400 多種藥用植物開展了WGS,其中124 種為藥典收載的種類(表1)。

    2.2 藥用植物WGS研究內(nèi)容不斷豐富

    前期,藥用植物WGS側(cè)重于基因注釋(genome annotation)、多倍化進(jìn)化分析(polyploidization analysis)和關(guān)鍵化學(xué)成分合成基因家族(gene family)預(yù)測,如蓖麻[5]、大麻[24]、川桑[26]等。而后,該方面研究增加了基因組進(jìn)化(genome evolution)、基因家族進(jìn)化(genefamily evolution)、全基因組復(fù)制(whole genome duplication)、基因調(diào)控元件分析、基因組與其生物特性關(guān)聯(lián)分析等內(nèi)容,如鐵皮石斛[6]、天麻[8]、銀杏[35]等。當(dāng)前,隨著基因組測序深度與質(zhì)量的提升,藥用植物WGS研究進(jìn)一步擴展到染色體水平組裝(chromosome levelassembly)、藥效成分合成途徑及其參與調(diào)控的轉(zhuǎn)錄因子等內(nèi)容,如紫丁香[12]、菘藍(lán)(Isatis indigotica)[36]、艾(Artemisia argyi)[37]等。

    2.3 藥用植物WGS 已邁入泛基因組時代

    隨著藥用植物WGS研究不斷豐富與深入,其泛基因組(pan-genome)研究亦開始起步。泛基因組是指一個物種或者一類物種中所有基因組信息的總和,包括存在于群體中所有個體的核心基因(core genes)及某些個體中不存在的可變基因(variable genes)。泛基因組涵蓋了單一參考基因組無法體現(xiàn)的物種遺傳變異信息,尤其是大的結(jié)構(gòu)變異,而這些變異可能與藥用植物性狀密切關(guān)聯(lián)。目前,藥用植物泛基因組研究報道集中于藥食兩用的品種,如萊菔子(蘿卜,Raphanussativus)、淡豆豉(大豆,Glycine max)等。ZHANG X H等[38]采用PacBio、Illumina、Bionano 和Hi-C 技術(shù)對覆蓋蘿卜屬的11 份典型種質(zhì)進(jìn)行了高質(zhì)量基因組denovo 組裝,挖掘出極大量基因組間的SNP、indel、SV、PAV、易位、倒位等遺傳變異;通過系統(tǒng)進(jìn)化分析表明,蘿卜屬野生蘿卜和栽培蘿卜在1.8 百萬年(Mya)前分化,野生蘿卜兩個亞種的分化發(fā)生在1.67 Mya前。LIU Y C等[39]利用PacBio(96x)+Bionano(277x)+Hi-C(136x)等方法對2 898 份大豆樣本基因組進(jìn)行了重測序,從中選取26 份代表性樣本進(jìn)行基因組組裝,構(gòu)建泛基因組圖譜,進(jìn)一步對26 個新組裝的基因組和ZH13 基因組進(jìn)行全基因組水平的比較基因組分析,將27 個大豆基因組的所有基因分類為57 492 個家族,其中20 623 個基因家族被定義為核心基因,8 163個基因家族被定義為次核心基因,28 679 個基因家族被定義為非必需基因,只有1 個基因家族被定義為特異性基因。上述研究極大促進(jìn)蘿卜、大豆的進(jìn)化和功能基因組學(xué)研究,對深入挖掘其藥材性狀、藥效成分形成機制具有重要意義。

    3 藥用植物WGS的應(yīng)用

    3.1 揭示藥用植物的進(jìn)化與馴化歷程

    WGS 是揭示藥用植物進(jìn)化歷程的理想手段,而其中有關(guān)全基因組復(fù)制(Whole genome duplications,WGD)、全基因組三倍體(whole-genome triplication,WGT)的分析是揭示進(jìn)化歷程的關(guān)鍵環(huán)節(jié)。人參、三七、刺五加(Eleutherococcus senticosus)同屬于五加科(Araliaceae),其中三七為二倍體,而人參、刺五加為四倍體。WGD分析結(jié)果顯示,人參、三七、刺五加大約在29.6 Mya經(jīng)歷了一次WGD事件,而人參、刺五加分別在2.2 Mya、13.0 Mya 又經(jīng)歷了1 次WGD 事件[9,40-41]。JIANG Z Q等[40]通過對WGS 數(shù)據(jù)進(jìn)一步分析后發(fā)現(xiàn),三七、人參在大約62 Mya 于五加科中形成進(jìn)化分枝,而后兩者在大約4.2 Mya 形成彼此獨立的進(jìn)化。XU Z C等[42]、LI J R等[43]通過對唇形科(Labiatae)的4 種藥用植物黃芩(Scutellaria baicalensis)、半枝蓮(Scutellaria barbata)、丹參(Salvia miltiorrhiza)、薰衣草(Lavandula angustifolia)的全基因組測序研究后發(fā)現(xiàn),黃芩、半枝蓮、丹參在大約60.7 Mya 發(fā)生了1 次WGD 事件,造成了染色體擴張與重排;而薰衣草經(jīng)歷了2次WGD事件,導(dǎo)致其萜類生物合成相關(guān)基因家族的擴張。TU L C等[10]通過WGS 分析證實了雷公藤在大約21.0 Mya經(jīng)歷1次WGT事件,導(dǎo)致其雷公藤內(nèi)酯生物合成基因產(chǎn)生增倍進(jìn)化,使其雷公藤內(nèi)酯合成能力得以加強。此外,GUO C等[15]、LIU H B等[44]通過薏苡(Coixlacryma-jobi)全基因組測序結(jié)果發(fā)現(xiàn),薏苡栽培種的種殼為較薄的“紙殼”,而野生種的種殼為較厚的“石殼”,種殼硬度變化是薏苡的主要馴化性狀,顯示出其人工選擇馴化的歷程。

    3.2 解析藥用植物中有效成分合成的分子機制

    WGS 為闡明藥用植物中有效成分合成的分子機制提供了一條高效的途徑。人參中達(dá)瑪烷型人參皂苷是其主要活性成分,但有關(guān)其合成的分子機制長期處于未知狀態(tài)。KIM N H等[7]等揭示了人參全基因組結(jié)構(gòu),共注釋了59 352 個基因,并從中鑒定出254 個可能參與達(dá)瑪烷型人參皂苷合成的基因,初步解析了人參皂苷生物合成的分子網(wǎng)絡(luò)。TU L C等[10]在揭示雷公藤全基因組信息的基礎(chǔ)上,結(jié)合轉(zhuǎn)錄組、代謝組數(shù)據(jù),構(gòu)建了其藥用成分雷公藤甲素的基因-代謝網(wǎng)絡(luò)圖譜,預(yù)測了57個可能參與雷公藤甲素生物合成的CYP(cytochromeP450)基因,并鑒定出催化產(chǎn)生雷公藤甲素中間體Dehydroabietic acid 的CYP728B70 基因。MIAO Y H等[37]揭示了艾的基因組,解析了其黃酮合成途徑中PAL(phenylalanine ammonia-lyase)、C4H(4-hydroxylase)、HCT(hydroxycinnamoyl transferase)、CHS(chalcone synthase)、F3H(favanone hydroxylase)等44 個基因,以及萜類生物合成中DXS(1-deoxy-Dxylose-5-phosphate synthase)、GPPS(Geranyl diphos原phate synthetase)、TPS(terpene synthase)等66 個基因,并初步明確了廣泛的基因擴增和串聯(lián)重復(fù)是艾中黃酮類化合物和揮發(fā)油含量豐富的成因。由此可見,全基因組測序為深入解析藥用植物有效成分生物合成途徑及合成調(diào)控等分子機制方面的研究提供了可靠保障。

    3.3 鑒定藥用植物種質(zhì)資源

    藥用植物種質(zhì)資源鑒定是中藥材質(zhì)量提升的重要保障。傳統(tǒng)上主要通過表觀和顯微技術(shù)相結(jié)合進(jìn)行鑒定,存在經(jīng)驗依賴性強、結(jié)果精準(zhǔn)度低等局限。隨著藥用植物WGS 的不斷發(fā)展,積累的DNA分子信息不斷豐富,促使基于DNA分子信息的DNA分子標(biāo)記技術(shù)在藥用植物種質(zhì)資源鑒定中得到更廣泛應(yīng)用[45-46]。孫嘉苓等[47]利用三七葉綠體基因組開發(fā)分子標(biāo)記,不僅鑒定區(qū)分了人參屬Panax 不同物種,而且對三七不同栽培居群也進(jìn)行了有效區(qū)分。陳璇等[48]運用WGS技術(shù)對野生和栽培大麻SNP位點進(jìn)行分析,發(fā)現(xiàn)其中雜合度、二等位多態(tài)性SNP存在較大差異。李慧等[49]基于紅天麻、烏天麻全基因組重測序結(jié)果,篩選獲得了特異性SNP位點,實現(xiàn)了紅、烏天麻及其雜交天麻的精準(zhǔn)鑒定。HUA Z Y等[50]收集了來自3 899 個植物物種的4 356條葉綠體基因組序列,并基于高通量測序等方法開發(fā)并驗證了白術(shù)(Atractylodes macrocephala)、人參等165種藥用植物的DSS 標(biāo)記,實現(xiàn)了藥用植物種類的規(guī)?;b定。WGS 相關(guān)研究結(jié)果的應(yīng)用,有力促進(jìn)了藥用植物種質(zhì)資源的精準(zhǔn)鑒定,支撐了藥材“有序、安全、有效”生產(chǎn)。

    3.4 推動藥用植物分子輔助育種

    分子輔助育種技術(shù)將現(xiàn)代分子生物學(xué)與傳統(tǒng)育種手段相結(jié)合,借助DNA 分子標(biāo)記對育種材料及其子代進(jìn)行選擇,可減少育種盲目性,縮短育種年限,提高育種效率。藥用植物存在生長周期長、雜合度高、育種目標(biāo)多樣化等育種障礙,而分子輔助育種技術(shù)是突破這些障礙的有效方法[51-52]。隨著WGS 研究的不斷深化,眾多藥用植物分子輔助育種的發(fā)展進(jìn)入了快車道。GUO C等[15]通過薏苡野生、栽培品種的WGS比較分析,構(gòu)建了其遺傳連鎖圖譜,定位到與其種殼耐壓性關(guān)聯(lián)的兩個QTL(Ccph1、Ccph2),研究表明Ccph1調(diào)控種殼厚度,Ccph2 調(diào)控種殼顏色,為下一步分子設(shè)計選育“紙殼”薏苡品種奠定了基礎(chǔ)。董林林等[53]采用簡化基因組測序技術(shù)等手段,篩選出高產(chǎn)、優(yōu)質(zhì)、抗逆等表型關(guān)聯(lián)的DNA 標(biāo)記,結(jié)合系統(tǒng)選育方法,選育出三七抗病新品種“苗鄉(xiāng)抗七1 號”。范宏虹等[54]利用全基因組關(guān)聯(lián)分析挖掘菊花耐寒性相關(guān)優(yōu)異等位變異和候選基因,篩選出24 個與耐寒性顯著關(guān)聯(lián)的SNP 位點,并基于主效可加互作模型的雙標(biāo)圖選育出抗寒性強且穩(wěn)定遺傳的品種。

    4 問題與展望

    2000 年以來,隨著基因組測序技術(shù)的不斷發(fā)展及生物信息學(xué)分析方法的優(yōu)化,目前有400 多種藥用植物的全基因組被研究揭示,基因組數(shù)據(jù)在解析藥用植物的進(jìn)化與馴化歷程、有效成分合成分子機制、種質(zhì)資源鑒定、分子輔助育種等方面被廣泛應(yīng)用,對推動藥用植物研究領(lǐng)域發(fā)展發(fā)揮著至關(guān)重要的作用。但是,藥用植物WGS研究還存在不少問題與局限,主要包括①完成WGS的藥用植物種類還比較少,與其總量相比占比很低,亟待進(jìn)一步加強;②目前的測序往往會結(jié)合二、三代測序技術(shù)進(jìn)行優(yōu)勢互補,但測序錯誤率仍相對較高,且基因組較大、雜合度較高種類的測序質(zhì)量和精準(zhǔn)組裝等方面還存在較大障礙;③雖然全基因組的堿基序列信息被揭示,但基因組中大量基因功能尚未被注釋,而且對多數(shù)已知功能基因的研究還處于克隆分離及生物信息學(xué)預(yù)測階段,缺乏直接有效的實驗驗證;④藥用植物基因組應(yīng)用研究尚處于起步階段,研究廣度和深度多有不足,亟待系統(tǒng)性提升與加強。

    目前,藥用植物WGS 研究及應(yīng)用處在快速發(fā)展階段,其最新技術(shù)和成果將會不斷整合并推動結(jié)構(gòu)基因組、功能基因組、轉(zhuǎn)錄組、蛋白質(zhì)組、表觀基因組、宏基因組、合成生物學(xué)、代謝組、生物信息學(xué)等研究領(lǐng)域的發(fā)展。未來主要研究方向包括①藥用植物全基因組學(xué)測序規(guī)模將進(jìn)一步擴大?!吨袊幍洹罚?020)版收載了604 種藥用植物,其中約1/5 的種類進(jìn)行了全基因組測序,未測序的種類將會受到重點關(guān)注。②藥用植物全基因組大數(shù)據(jù)將逐步形成并得到完善。目前研究顯示藥用植物全基因組測序已經(jīng)從探索未知的階段進(jìn)入了大數(shù)據(jù)關(guān)聯(lián)研究的階段。如何運用大數(shù)據(jù)等分析手段從這些海量測序信息中挖掘篩選有效的信息將是一個有待持續(xù)攻關(guān)的方向。③藥用植物后基因組學(xué)研究與應(yīng)用將蓬勃發(fā)展起來。其中:基礎(chǔ)研究方面包括轉(zhuǎn)錄組、代謝組、蛋白質(zhì)組、表型組及合成生物學(xué)等將不斷深入,而且多組學(xué)聯(lián)合分析將成為新的研究趨勢;在應(yīng)用研究方面,藥材“優(yōu)形、優(yōu)質(zhì)”特征形成機制及目標(biāo)性狀關(guān)聯(lián)基因?qū)⒅鸩降靡越馕?,藥用植物分子輔助育種過程中目標(biāo)性狀精準(zhǔn)設(shè)計、優(yōu)異等位基因聚合等關(guān)鍵技術(shù)環(huán)節(jié)將有所突破,育種效率將會得到大幅提升。

    綜上而言,隨著全基因組測序技術(shù)日益革新與發(fā)展,大量藥用植物WGS 業(yè)已完成,相關(guān)基礎(chǔ)理論與實踐應(yīng)用的研究將得到進(jìn)一步深化和加強,藥用植物領(lǐng)域研究也將被推入空前繁榮的新時代。

    參考文獻(xiàn):

    [1] 劉景芳,李維林,王莉,等. 多組學(xué)技術(shù)及其在生命科學(xué)研究

    中應(yīng)用概述[J]. 生物工程學(xué)報,2022,38(10):3581-3593.

    [2] ORLOV Y L, ANASHKINA A A. Life: computational

    genomics applications in life sciences[J]. Life(Basel), 2021,

    11(11): 1211.

    [3] 鄧群,郎濤,夏建新. 藥用植物資源開發(fā)利用現(xiàn)狀及其

    發(fā)展[J]. 中央民族大學(xué)學(xué)報(自然科學(xué)版),2016,25(1):

    55-59.

    [4] CHENG Q Q, OUYANG Y, TANG Z Y, et al.

    Review on the development and applications of medicinal

    plant genomes[J]. Frontiersin Plant Science, 2021, 12:

    791219.

    [5] CHAN A P, CRABTREE J, ZHAO Q, et al. Draft

    genome sequence of the oilseed species Ricinus communis

    [J]. Nature Biotechnology, 2010, 28(9): 951-956.

    [6] YAN L, WANG X, LIU H, et al. The genome of

    Dendrobium officinale illuminates the biology of the

    important traditional Chinese orchid herb[J]. Molecular

    Plant, 2015, 8(6): 922-934.

    [7] KIM N H, JAYAKODI M, LEE S C, et al. Genome

    and evolution of the shade-requiring medicinal herb

    Panax ginseng[J]. Plant Biotechnology Journal, 2018, 16(11):

    1904-1917.

    [8] YUAN Y, JIN X H, LIU J, et al. The Gastrodia

    elata genome provides insights into plant adaptation to

    heterotrophy[J]. Nature Communications, 2018, 9(1): 1615.

    [9] PU X D, LI Z, TIAN Y, et al. The honeysuckle

    genome provides insight into the molecular mechanism

    of carotenoid metabolism underlying dynamic flower

    coloration[J]. The New Phytologist, 2020, 227(3): 930-943.

    [10] TU L C, SU P, ZHANG Z R, et al. Genome of

    Tripterygium wilfordii and identification of cytochrome

    P450 involved in triptolide biosynthesis [ J ]. Nature

    Communications, 2020, 11(1): 971.

    [11] XIONG X Y, GOU J B, LIAO Q G, et al. The Taxus

    genome provides insights into paclitaxel biosynthesis [J].

    Nature Plants, 2021, 7(8): 1026-1036.

    [12] MA B, WU J, SHI T L, et al. Lilac (Syringa oblata)

    genome provides insights into its evolution and molecular

    mechanism of petal color change[J]. Communications

    Biology, 2022, 5(1): 686.

    [13] BORNOWSKI N, HAMILTON J P, LIAO P, et al.

    Genome sequencing of four culinary herbs reveals

    terpenoid genes underlying chemodiversity in the

    Nepetoideae[J]. DNA Research, 2020, 27(3): dsaa016.

    [14] ZHANG X H, LIU T J, WANG J L, et al. Pan-genome

    of Raphanus highlights genetic variation and introgression

    among domesticated, wild, and weedy radishes

    [J]. Molecular Plant, 2021, 14(12): 2032-2055.

    [15] GUO C, WANG Y N, YANG A G, et al. The

    Coix genome provides insights into Panicoideae evolution

    and papery hull domestication[J]. Molecular Plant,

    2020, 13(2): 309-320.

    [16] VENTER J C, SMITH H O, ADAMS M D. The

    sequence of the human genome[J]. Clinical Chemistry,

    2015, 61(9): 1207-1208.

    [17] INITIATIVEA G. Analysis of the genome sequence

    of the flowering plant Arabidopsis thaliana[J]. Nature,

    2000, 408(6814): 796-815.

    [18] TUSKAN G A, DIFAZIO S, JANSSON S, et al. The

    genome of black cottonwood, Populus trichocarpa (Torr.

    amp; Gray)[J]. Science, 2006, 313(5793): 1596-1604.

    [19] JAILLON O, AURY J M, NOEL B, et al. The

    grapevine genome sequence suggests ancestral hexaploidization

    in major angiosperm phyla[J]. Nature, 2007,

    449(7052): 463-467.

    [20] MING R, HOU S B, FENG Y, et al. The draft

    genome of the transgenic tropical fruit tree papaya

    (Carica papaya Linnaeus) [J]. Nature, 2008, 452 (7190):

    991-996.

    [21] SANGER F, NICKLEN S, COULSON A R, et al.

    DNA sequencing with chain-terminating inhibitors [J].

    Proceedings of the National Academy of Sciences of

    the United States of America, 1977, 74(12): 5463-5467.

    [22] HU T S, CHITNIS N, MONOS D, et al. Next-generation

    sequencing technologies: An overview[J]. Human

    Immunology, 2021, 82(11): 801-811.

    [23] LIN B, HUI J N, MAO H J, et al. Nanopore technology

    and its applications in gene sequencing[J]. Biosensors,

    2021, 11(7): 214.

    [24] BAKEL H V, STOUT J M, COTE A G, et al. The

    draft genome and transcriptome of Cannabis sativa[J].

    Genome Biology, 2011, 12(10): R102.

    [25] WANG Z W, HOBSON N, GALINDO L, et al.

    The genome of flax (Linum usitatissimum) assembled de

    novo from short shotgun sequence reads[J]. The Plant

    Journal, 2012, 72(3): 461-473.

    [26] HE N J, ZHANG C, QI X W, et al. Draft genome

    sequence of the mulberry tree Morus notabilis[J]. Nature

    Communications, 2013, 4: 2445.

    [27] WANG Y, FAN G Y, LIU Y M, et al. The sacred lotus

    genome provides insights into the evolution of flowering

    plants[J]. The Plant Journal,2013,76(4):557-567.

    [28] UPADHYAY A K, CHACKO A R, GANDHIMATHI

    A, et al. Genome sequencing of herb Tulsi

    (Ocimum tenuiflorum) unravels key genes behind its strong

    medicinal properties[J]. BMC Plant Biology,2015,15(1):212.

    [29] GUAN R, ZHAO Y P, ZHANG H, et al. Draft genome

    of the living fossil Ginkgo biloba[J]. Gigascience, 2016, 5(1):

    49.

    [30] ZHAO D Y, HAMILTON J P, PHAM G M, et al. De

    novo genome assembly of Camptotheca acuminata, a

    natural source of the anti-cancer compound camptothecin[

    J]. Gigascience,2017, 6(9): 1-7.

    [31] WANG X, XU Y T, ZHANG S Q, et al. Genomic

    analyses of primitive, wild and cultivated citrus provide

    insights into asexual reproduction [J]. Nature Genetics,

    2017, 49(5): 765-772.

    [32] LIN Y L, MIN J M, LAI R L, et al. Genome-wide

    sequencing of longan (Dimocarpus longan Lour.) provides

    insights into molecular basis of its polyphenol-

    rich characteristics[J]. Gigascience, 2017, 6(5): 1-14.

    [33] ZHANG Y H, ZHENG L L, ZHENG Y, et al. Assembly

    and annotation of a draft genome of the medicinal

    plant Polygonum cuspidatum [J]. Frontiers in Plant

    Science, 2019, 10: 1274.

    [34] QIN S S, WU L Q, WEI K H, et al. A draft

    genome for Spatholobus suberectus[J]. Scientific Data,

    2019, 6(1): 1-9.

    [35] LIU H L, WANG X B, WANG G B, et al. The

    nearly complete genome of Ginkgo biloba illuminates

    gymnosperm evolution[J]. Nature Plants,2021,7(6):748-756.

    [36] KANG M H, WU H L, YANG Q, et al. A chromosome-

    scale genome assembly of Isatis indigotica, an

    important medicinal plant used in traditional Chinese

    medicine: An Isatis genome[J]. Horticulture Research,

    2020, 7: 18.

    [37] MIAO Y H, LUO D D, ZHAO T T, et al.

    Genome sequencing reveals chromosome fusion and

    extensive expansion of genes related to secondary

    metabolism in Artemisia argyi[J]. Plant Biotechnology

    Journal, 2022, 20(10): 1902-1915.

    [38] ZHANG X H, LIU T J, WANG J L, et al.

    Pan-genome of Raphanus highlights genetic variation

    and introgression among domesticated, wild, and weedy

    radishes[J]. Molecular Plant, 2021, 14(12): 2032-2055.

    [39] LIU Y C, DU H L, LI P C, et al. Pan-genome of wild

    and cultivated soybeans[J]. Cell,2020,182(1):162-176.

    [40] JIANG Z Q, TU L C, YANG W F, et al. The chromosome-

    level reference genome assembly for Panax

    notoginseng and insights into ginsenoside biosynthesis[J].

    Plant Communications, 2020, 2: 100113.

    [41] YANG Z J, CHEN S S, WANG S F, et al. Chromosomal-

    scale genome assembly of Eleutherococcus

    senticosus provides insights into chromosome evolution

    in Araliaceae[J]. Molecular Ecology Resources, 2021, 21

    (7): 2204-2220.

    [42] XU Z C, GAO R R, PU X D, et al. Comparative

    genome analysis of Scutellaria baicalensis and Scutel原

    laria barbata reveals the evolution of active flavonoid

    biosynthesis [J]. Genomics, Proteomics amp; Bioinformatics,

    2020, 18(3): 230-240.

    [43] LI J R, WANG Y M, DONG Y M, et al. The chromosome-

    based lavender genome provides new insights

    into Lamiaceae evolution and terpenoid biosynthesis [J].

    Horticulture Research, 2021, 8(1): 53.

    [44] LIU H B, SHI J P, CAI Z X, et al. Evolution and

    domestication footprints uncovered from the genomes

    of Coix[J]. Molecular Plant, 2020, 13(2): 295-308.

    [45] 黃璐琦,錢丹,鄧超. 雙分子標(biāo)記法的構(gòu)建及在中藥研

    究中的應(yīng)用[J]. 中國中藥雜志,2015,40(2):165-168.

    [46] KUMAR R, YADAV S S, MISHRA S K, et al.

    Molecular markers of aromatic and medicinal plants: a

    review[J]. Pharma Science Monitor, 2013, 4(4): 344-362.

    [47] 孫嘉苓,韓巖,崔秀明,等. 三七葉綠體分子標(biāo)記的開發(fā)與

    應(yīng)用研究[J]. 中國中藥雜志,2020,45(6):1342-1349.

    [48] 陳璇,郭蓉,王璐,等. 基于全基因組重測序的野生型

    大麻和栽培型大麻的多態(tài)性SNP 分析[J]. 分子植物

    育種,2018,16(3):893-897.

    [49] 李慧,錢潤,田娜,等. 紅天麻、烏天麻及其雜交天麻的

    PCR鑒別[J]. 中國中藥雜志,2020,45(15):3666-3671.

    [50] HUA Z Y, JIANG C, SONG S H, et al. Accurate

    identification of taxon-specific molecular markers in

    plants based on DNA signature sequence[J]. Molecular

    Ecology Resources, 2022, 23(1): 106-117.

    [51] 馬小軍,莫長明. 藥用植物分子育種展望[J]. 中國中藥

    雜志,2017,42(11):2021-2031.

    [52] 錢潤,周駿輝,楊健,等. 中藥材分子標(biāo)記輔助育種技術(shù)

    研究進(jìn)展[J]. 中國中藥雜志,2020,45(20):4812-4818.

    [53] 董林林,陳中堅,王勇,等. 藥用植物DNA 標(biāo)記輔助育

    種(一):三七抗病品種選育研究[J]. 中國中藥雜志,

    2017,42(1):56-62.

    [54] 范宏虹,徐婷婷,蘇江碩,等. 切花菊耐寒性相關(guān)SNP

    位點挖掘與候選基因分析[J]. 園藝學(xué)報,2019,46(11):

    2201-2212.

    [55] KRISHNAN N M, PATTNAIK S, JAIN P, et al. A

    draft of the genome and four transcriptomes of a

    medicinal and pesticidal angiosperm Azadirachta indica

    [J]. BMC Genomics, 2012, 13(1):464.

    [56]~[438]參考文獻(xiàn)略。

    責(zé)任編輯:李曉瑞

    国产亚洲精品第一综合不卡| 波野结衣二区三区在线| 在线免费观看不下载黄p国产| 大片电影免费在线观看免费| 亚洲色图 男人天堂 中文字幕| 久久99一区二区三区| 五月开心婷婷网| 精品国产一区二区三区四区第35| 亚洲国产欧美日韩在线播放| 99久国产av精品国产电影| 久久国产亚洲av麻豆专区| 国产精品免费大片| 中文字幕人妻熟女乱码| 午夜av观看不卡| 日韩av在线免费看完整版不卡| 丝袜美腿诱惑在线| 亚洲成av片中文字幕在线观看| 日韩一区二区视频免费看| 一个人免费看片子| 亚洲少妇的诱惑av| 精品国产超薄肉色丝袜足j| 777久久人妻少妇嫩草av网站| 日韩欧美精品免费久久| 亚洲成人手机| av国产精品久久久久影院| 在线观看免费午夜福利视频| 免费高清在线观看视频在线观看| 亚洲欧美一区二区三区黑人| 成人亚洲欧美一区二区av| 国产成人精品福利久久| 老汉色av国产亚洲站长工具| 午夜福利,免费看| 亚洲成人国产一区在线观看 | 亚洲精品成人av观看孕妇| 男人爽女人下面视频在线观看| 超碰成人久久| 精品久久久精品久久久| 午夜福利视频精品| 美女福利国产在线| 9色porny在线观看| 国产亚洲午夜精品一区二区久久| 在线观看三级黄色| 日韩一卡2卡3卡4卡2021年| 99香蕉大伊视频| 最黄视频免费看| av免费观看日本| 欧美精品一区二区免费开放| 少妇猛男粗大的猛烈进出视频| 久久青草综合色| 在线观看国产h片| 国产片内射在线| 2018国产大陆天天弄谢| 啦啦啦在线观看免费高清www| 一本—道久久a久久精品蜜桃钙片| 人人妻人人澡人人看| 日日撸夜夜添| 色吧在线观看| 亚洲av欧美aⅴ国产| 欧美激情高清一区二区三区 | 美女福利国产在线| 免费黄色在线免费观看| 亚洲精品国产av蜜桃| 日韩精品免费视频一区二区三区| 免费在线观看完整版高清| 亚洲第一青青草原| 亚洲少妇的诱惑av| 老司机亚洲免费影院| 人人澡人人妻人| 最近中文字幕高清免费大全6| 最近中文字幕高清免费大全6| 丝袜美足系列| 老司机深夜福利视频在线观看 | 大话2 男鬼变身卡| 免费黄网站久久成人精品| 丰满少妇做爰视频| 免费观看av网站的网址| 最近最新中文字幕大全免费视频 | 国产欧美日韩综合在线一区二区| 欧美日韩亚洲综合一区二区三区_| 精品国产乱码久久久久久男人| 男女免费视频国产| 久久久亚洲精品成人影院| 只有这里有精品99| 国产熟女午夜一区二区三区| 亚洲熟女精品中文字幕| 欧美人与性动交α欧美精品济南到| 日本av免费视频播放| 哪个播放器可以免费观看大片| 国产成人午夜福利电影在线观看| 亚洲视频免费观看视频| 黄片小视频在线播放| 19禁男女啪啪无遮挡网站| 成年人免费黄色播放视频| 精品酒店卫生间| 午夜福利在线免费观看网站| 亚洲人成77777在线视频| 人体艺术视频欧美日本| 免费观看人在逋| 美女大奶头黄色视频| 1024视频免费在线观看| 91国产中文字幕| 国产精品秋霞免费鲁丝片| 性少妇av在线| av在线观看视频网站免费| 在现免费观看毛片| 久热爱精品视频在线9| 亚洲欧美精品综合一区二区三区| 亚洲欧美精品综合一区二区三区| 欧美成人午夜精品| 亚洲三区欧美一区| 黄色 视频免费看| 免费黄网站久久成人精品| 天天躁狠狠躁夜夜躁狠狠躁| 成年动漫av网址| 在线精品无人区一区二区三| 午夜福利,免费看| 欧美最新免费一区二区三区| 男女之事视频高清在线观看 | 18禁观看日本| 成年av动漫网址| 亚洲图色成人| 免费少妇av软件| 国产伦理片在线播放av一区| 亚洲av日韩精品久久久久久密 | 国产免费又黄又爽又色| 人妻 亚洲 视频| 亚洲视频免费观看视频| 中文字幕人妻丝袜一区二区 | av一本久久久久| 校园人妻丝袜中文字幕| 国产成人a∨麻豆精品| 夜夜骑夜夜射夜夜干| 成年人午夜在线观看视频| a 毛片基地| 精品国产超薄肉色丝袜足j| 午夜精品国产一区二区电影| 高清在线视频一区二区三区| 日韩视频在线欧美| 两个人免费观看高清视频| av视频免费观看在线观看| 亚洲欧美成人精品一区二区| 黄色视频不卡| 欧美日韩一区二区视频在线观看视频在线| 免费女性裸体啪啪无遮挡网站| 最近最新中文字幕大全免费视频 | 国产精品一国产av| 看非洲黑人一级黄片| 国产国语露脸激情在线看| 视频区图区小说| 免费不卡黄色视频| 久久久精品94久久精品| 黄片播放在线免费| 午夜日韩欧美国产| 欧美日韩精品网址| 久久久久久人妻| 亚洲精品中文字幕在线视频| 一区二区三区激情视频| 日韩免费高清中文字幕av| 国产免费视频播放在线视频| 精品一区二区免费观看| 欧美国产精品一级二级三级| 色视频在线一区二区三区| 日本爱情动作片www.在线观看| av国产久精品久网站免费入址| 中文乱码字字幕精品一区二区三区| 国产精品99久久99久久久不卡 | 丝袜美腿诱惑在线| 在线看a的网站| 国产黄频视频在线观看| 亚洲一卡2卡3卡4卡5卡精品中文| 99久久精品国产亚洲精品| 激情五月婷婷亚洲| 一二三四中文在线观看免费高清| 亚洲精品国产av成人精品| 亚洲精品第二区| 亚洲av综合色区一区| 久久久久久久国产电影| 久久婷婷青草| 国产不卡av网站在线观看| 亚洲中文av在线| 制服丝袜香蕉在线| 欧美日韩成人在线一区二区| 国产精品亚洲av一区麻豆 | 午夜精品国产一区二区电影| 国产免费又黄又爽又色| 91精品三级在线观看| 久久女婷五月综合色啪小说| 国产成人精品无人区| 日本色播在线视频| 国产一区二区在线观看av| 一边亲一边摸免费视频| 久久热在线av| netflix在线观看网站| 99精国产麻豆久久婷婷| 久久久亚洲精品成人影院| 日本wwww免费看| 国产精品免费大片| 国产免费一区二区三区四区乱码| 亚洲欧美日韩另类电影网站| 新久久久久国产一级毛片| 国产精品二区激情视频| 男女之事视频高清在线观看 | 久久久精品区二区三区| 无遮挡黄片免费观看| 国产在线视频一区二区| 国产成人精品久久久久久| 最新的欧美精品一区二区| 日韩大片免费观看网站| 亚洲色图综合在线观看| 亚洲四区av| 嫩草影视91久久| 老汉色av国产亚洲站长工具| 亚洲精品aⅴ在线观看| 又大又黄又爽视频免费| 亚洲精品国产色婷婷电影| 夜夜骑夜夜射夜夜干| 欧美亚洲日本最大视频资源| 日韩欧美一区视频在线观看| 久久久久久久久久久免费av| 在线观看免费日韩欧美大片| 一本色道久久久久久精品综合| 久久久久久久精品精品| 亚洲国产精品成人久久小说| av线在线观看网站| 午夜福利,免费看| 国产免费福利视频在线观看| 欧美日韩福利视频一区二区| 国产精品 欧美亚洲| 99久久99久久久精品蜜桃| 丰满饥渴人妻一区二区三| 巨乳人妻的诱惑在线观看| 国产精品.久久久| 又大又爽又粗| 午夜激情久久久久久久| 午夜福利影视在线免费观看| 晚上一个人看的免费电影| 男人添女人高潮全过程视频| 女人久久www免费人成看片| 一区二区三区精品91| 一级毛片电影观看| 高清欧美精品videossex| 狠狠精品人妻久久久久久综合| 人妻 亚洲 视频| 亚洲七黄色美女视频| 精品国产一区二区久久| 久久久久人妻精品一区果冻| 国产老妇伦熟女老妇高清| www.自偷自拍.com| 成年人免费黄色播放视频| 在线天堂最新版资源| 亚洲欧美成人综合另类久久久| 两个人看的免费小视频| 午夜激情av网站| 婷婷色综合www| 国产乱人偷精品视频| 欧美日本中文国产一区发布| av在线观看视频网站免费| 最近最新中文字幕大全免费视频 | 男女边摸边吃奶| 99热网站在线观看| 久久热在线av| 免费日韩欧美在线观看| av在线老鸭窝| 欧美日韩综合久久久久久| 欧美人与性动交α欧美精品济南到| 多毛熟女@视频| 免费日韩欧美在线观看| 久久精品国产a三级三级三级| 欧美日本中文国产一区发布| 天天操日日干夜夜撸| av免费观看日本| 亚洲国产av新网站| 我的亚洲天堂| 亚洲熟女毛片儿| 亚洲精华国产精华液的使用体验| 久久亚洲国产成人精品v| 男女免费视频国产| 亚洲一码二码三码区别大吗| 亚洲激情五月婷婷啪啪| 又粗又硬又长又爽又黄的视频| 欧美日韩国产mv在线观看视频| 熟妇人妻不卡中文字幕| 五月开心婷婷网| 国产熟女午夜一区二区三区| 校园人妻丝袜中文字幕| 欧美亚洲日本最大视频资源| 深夜精品福利| 国产精品偷伦视频观看了| 男女高潮啪啪啪动态图| 国产精品久久久久成人av| 91成人精品电影| 欧美日韩综合久久久久久| 性少妇av在线| 最近手机中文字幕大全| 亚洲美女搞黄在线观看| 91精品三级在线观看| 只有这里有精品99| 精品第一国产精品| 久久久久网色| 国产精品麻豆人妻色哟哟久久| 日韩熟女老妇一区二区性免费视频| 欧美日韩福利视频一区二区| 久久 成人 亚洲| 女的被弄到高潮叫床怎么办| 人人妻人人爽人人添夜夜欢视频| 国产黄色免费在线视频| 精品第一国产精品| 久久精品久久久久久久性| 亚洲男人天堂网一区| 国产精品成人在线| 亚洲精品久久午夜乱码| 婷婷成人精品国产| 国产一区二区三区综合在线观看| 99香蕉大伊视频| 一二三四中文在线观看免费高清| 麻豆乱淫一区二区| 热99久久久久精品小说推荐| 亚洲精品,欧美精品| 国产精品女同一区二区软件| 蜜桃国产av成人99| 人人妻人人爽人人添夜夜欢视频| 制服诱惑二区| 久久青草综合色| 欧美亚洲日本最大视频资源| 午夜免费男女啪啪视频观看| 亚洲一码二码三码区别大吗| 久久久精品94久久精品| 在线免费观看不下载黄p国产| 80岁老熟妇乱子伦牲交| 男女边摸边吃奶| 成人漫画全彩无遮挡| 亚洲精品国产区一区二| 亚洲在久久综合| 国产国语露脸激情在线看| 国产成人a∨麻豆精品| 女人高潮潮喷娇喘18禁视频| 黄频高清免费视频| 精品久久蜜臀av无| 亚洲一区中文字幕在线| 日韩视频在线欧美| 午夜福利网站1000一区二区三区| 男女边吃奶边做爰视频| 亚洲欧美日韩另类电影网站| 日本黄色日本黄色录像| 亚洲欧洲国产日韩| 亚洲自偷自拍图片 自拍| 在线观看一区二区三区激情| 中文天堂在线官网| 久久影院123| 亚洲精品中文字幕在线视频| 成人免费观看视频高清| 久久久久视频综合| 曰老女人黄片| 人人妻人人爽人人添夜夜欢视频| 69精品国产乱码久久久| 9191精品国产免费久久| 午夜福利视频在线观看免费| 男女下面插进去视频免费观看| 国产男女内射视频| 国产亚洲一区二区精品| 精品酒店卫生间| 最近的中文字幕免费完整| 在线观看免费视频网站a站| 久久天躁狠狠躁夜夜2o2o | 看非洲黑人一级黄片| 王馨瑶露胸无遮挡在线观看| 无遮挡黄片免费观看| 欧美激情高清一区二区三区 | 男人舔女人的私密视频| 中文字幕另类日韩欧美亚洲嫩草| 国产男人的电影天堂91| 久久久久精品性色| 国产成人午夜福利电影在线观看| 亚洲成国产人片在线观看| 亚洲第一av免费看| 亚洲av欧美aⅴ国产| 国产亚洲最大av| 欧美成人精品欧美一级黄| 99精国产麻豆久久婷婷| 搡老乐熟女国产| 美女视频免费永久观看网站| 黑丝袜美女国产一区| 国产成人啪精品午夜网站| 亚洲精品中文字幕在线视频| 国产精品.久久久| 色综合欧美亚洲国产小说| 少妇 在线观看| 国产一区二区三区综合在线观看| 亚洲国产av影院在线观看| 啦啦啦视频在线资源免费观看| 丁香六月天网| 1024视频免费在线观看| 国产成人av激情在线播放| 久久精品国产a三级三级三级| 美女大奶头黄色视频| 精品少妇久久久久久888优播| 国产野战对白在线观看| 视频在线观看一区二区三区| 丝袜美足系列| 两个人免费观看高清视频| 午夜老司机福利片| 午夜老司机福利片| 啦啦啦在线免费观看视频4| 亚洲国产毛片av蜜桃av| 老鸭窝网址在线观看| 日日摸夜夜添夜夜爱| 亚洲国产看品久久| 一区二区三区精品91| 王馨瑶露胸无遮挡在线观看| 成人午夜精彩视频在线观看| 日本一区二区免费在线视频| 国产精品熟女久久久久浪| 欧美老熟妇乱子伦牲交| 国产麻豆69| 久热爱精品视频在线9| 日韩熟女老妇一区二区性免费视频| 国产亚洲欧美精品永久| 中文字幕人妻熟女乱码| 黑人欧美特级aaaaaa片| videos熟女内射| 超色免费av| 一区二区三区乱码不卡18| 久久久欧美国产精品| 精品国产一区二区久久| 2021少妇久久久久久久久久久| 欧美黑人欧美精品刺激| 操美女的视频在线观看| 亚洲欧美日韩另类电影网站| 亚洲av日韩精品久久久久久密 | 国产成人一区二区在线| av网站在线播放免费| 色婷婷久久久亚洲欧美| 熟女av电影| 日日撸夜夜添| 在现免费观看毛片| 久久av网站| 爱豆传媒免费全集在线观看| 黄网站色视频无遮挡免费观看| 90打野战视频偷拍视频| √禁漫天堂资源中文www| 美女视频免费永久观看网站| 高清黄色对白视频在线免费看| 久久久久久久精品精品| 人人妻,人人澡人人爽秒播 | 国产免费又黄又爽又色| 欧美精品一区二区免费开放| 中文字幕精品免费在线观看视频| 国产精品国产av在线观看| 色婷婷av一区二区三区视频| 亚洲精品国产色婷婷电影| 久久精品亚洲熟妇少妇任你| 国产成人一区二区在线| 自拍欧美九色日韩亚洲蝌蚪91| 欧美成人精品欧美一级黄| 久久久亚洲精品成人影院| 热re99久久国产66热| 成年动漫av网址| 一区二区av电影网| 免费黄色在线免费观看| 最近中文字幕2019免费版| 成人漫画全彩无遮挡| 国产色婷婷99| 日韩不卡一区二区三区视频在线| 侵犯人妻中文字幕一二三四区| 国产熟女欧美一区二区| 黄色怎么调成土黄色| 日韩精品免费视频一区二区三区| 久久99一区二区三区| 亚洲人成电影观看| 丰满少妇做爰视频| 中文字幕人妻丝袜一区二区 | 欧美日韩亚洲高清精品| 欧美在线一区亚洲| 嫩草影院入口| 在线观看人妻少妇| 男人操女人黄网站| 国产黄色视频一区二区在线观看| 美女脱内裤让男人舔精品视频| 日韩成人av中文字幕在线观看| 午夜福利乱码中文字幕| 日韩一卡2卡3卡4卡2021年| 最近2019中文字幕mv第一页| 成人亚洲欧美一区二区av| 成人18禁高潮啪啪吃奶动态图| 波多野结衣av一区二区av| 国产欧美日韩综合在线一区二区| 少妇人妻精品综合一区二区| 美女脱内裤让男人舔精品视频| 久久综合国产亚洲精品| 精品国产一区二区三区四区第35| 久久这里只有精品19| 免费看av在线观看网站| 色精品久久人妻99蜜桃| 欧美国产精品va在线观看不卡| 黑丝袜美女国产一区| 中文字幕av电影在线播放| 人人妻人人爽人人添夜夜欢视频| 丰满乱子伦码专区| 成年动漫av网址| 成人国产av品久久久| 中文字幕人妻熟女乱码| 天堂俺去俺来也www色官网| 天天躁夜夜躁狠狠久久av| 热99国产精品久久久久久7| 大陆偷拍与自拍| 嫩草影院入口| 亚洲av日韩精品久久久久久密 | 亚洲精品美女久久久久99蜜臀 | 国产精品成人在线| 国产亚洲av高清不卡| 久久久久久久久久久久大奶| 亚洲综合精品二区| 亚洲av日韩精品久久久久久密 | 免费不卡黄色视频| www.自偷自拍.com| 日韩伦理黄色片| 欧美日韩亚洲高清精品| 欧美精品一区二区免费开放| 亚洲欧美激情在线| 亚洲av中文av极速乱| 国产精品一区二区精品视频观看| 男女下面插进去视频免费观看| 亚洲自偷自拍图片 自拍| 欧美日韩综合久久久久久| 丁香六月欧美| 国产一区二区三区综合在线观看| 色视频在线一区二区三区| 人妻一区二区av| 久久热在线av| 国产伦理片在线播放av一区| 国产精品一区二区在线不卡| 久久久欧美国产精品| 亚洲精品国产色婷婷电影| 亚洲欧美成人精品一区二区| 搡老乐熟女国产| 国产成人午夜福利电影在线观看| svipshipincom国产片| 成人午夜精彩视频在线观看| 黑人猛操日本美女一级片| 国产av一区二区精品久久| 激情五月婷婷亚洲| 免费看不卡的av| 亚洲国产精品一区三区| 欧美黑人欧美精品刺激| 综合色丁香网| 久久久久精品人妻al黑| 久久青草综合色| 岛国毛片在线播放| 一区二区三区四区激情视频| 久久鲁丝午夜福利片| 黄色怎么调成土黄色| 亚洲国产日韩一区二区| 啦啦啦在线免费观看视频4| 五月开心婷婷网| 国产精品熟女久久久久浪| 欧美激情 高清一区二区三区| 丰满乱子伦码专区| 国产又色又爽无遮挡免| 国产精品欧美亚洲77777| 国产午夜精品一二区理论片| 久久ye,这里只有精品| a级片在线免费高清观看视频| 久久ye,这里只有精品| 午夜福利,免费看| 精品少妇久久久久久888优播| 欧美日本中文国产一区发布| 国产99久久九九免费精品| 亚洲中文av在线| 色播在线永久视频| 叶爱在线成人免费视频播放| 老司机深夜福利视频在线观看 | 天堂俺去俺来也www色官网| 两性夫妻黄色片| 欧美日韩亚洲高清精品| av不卡在线播放| 黑丝袜美女国产一区| 一本大道久久a久久精品| 高清黄色对白视频在线免费看| 亚洲欧美中文字幕日韩二区| 色精品久久人妻99蜜桃| 国产精品国产三级国产专区5o| 国产精品偷伦视频观看了| 亚洲精品中文字幕在线视频| 精品国产乱码久久久久久小说| 丝瓜视频免费看黄片| 国产97色在线日韩免费| 精品亚洲成国产av| 中文精品一卡2卡3卡4更新| 欧美黑人精品巨大| 国产视频首页在线观看| 国产精品二区激情视频| 中文字幕av电影在线播放| 天天添夜夜摸| 国产精品熟女久久久久浪| 母亲3免费完整高清在线观看| 99热全是精品| av免费观看日本| 亚洲av成人不卡在线观看播放网 | 涩涩av久久男人的天堂| 成年女人毛片免费观看观看9 | 十分钟在线观看高清视频www| 在线观看一区二区三区激情| √禁漫天堂资源中文www| 欧美老熟妇乱子伦牲交| av国产精品久久久久影院| 久久国产精品大桥未久av| av不卡在线播放| 亚洲精品自拍成人| 亚洲精品在线美女| 人人妻人人澡人人看| 免费黄频网站在线观看国产| 国产精品免费大片| 尾随美女入室| 久久韩国三级中文字幕| 免费观看人在逋|