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

    環(huán)境DNA metabarcoding及其在生態(tài)學(xué)研究中的應(yīng)用

    2016-10-24 09:16:34徐海根
    生態(tài)學(xué)報(bào) 2016年15期
    關(guān)鍵詞:條形碼測序物種

    陳 煉,吳 琳,劉 燕,徐海根,*

    1 環(huán)境保護(hù)部南京環(huán)境科學(xué)研究所, 南京 210042 2 江蘇第二師范學(xué)院生命科學(xué)與化學(xué)化工學(xué)院, 南京 210013

    ?

    環(huán)境DNA metabarcoding及其在生態(tài)學(xué)研究中的應(yīng)用

    陳煉1, 2,吳琳2,劉燕1,徐海根1,*

    1 環(huán)境保護(hù)部南京環(huán)境科學(xué)研究所, 南京210042 2 江蘇第二師范學(xué)院生命科學(xué)與化學(xué)化工學(xué)院, 南京210013

    環(huán)境DNA metabarcoding(eDNA metabarcoding)是指利用環(huán)境樣本(如土壤、水、糞便等)中分離的DNA進(jìn)行高通量的多個物種(或高級分類單元)鑒定的方法。近年來,該方法引起了學(xué)者的廣泛關(guān)注,逐漸應(yīng)用于生物多樣性研究、水生生物監(jiān)測、珍稀瀕危物種和外來入侵物種檢測等生態(tài)學(xué)領(lǐng)域。介紹環(huán)境DNA metabarcoding的含義和研究方法;重點(diǎn)介紹環(huán)境DNA metabarcoding在物種監(jiān)測、生物多樣性研究和食性分析等生態(tài)學(xué)領(lǐng)域中的應(yīng)用;總結(jié)環(huán)境DNA metabarcoding應(yīng)用于生態(tài)學(xué)研究領(lǐng)域面臨的挑戰(zhàn)并對該方法的發(fā)展進(jìn)行展望。

    環(huán)境DNA;metabarcoding;生物多樣性;物種監(jiān)測

    在生態(tài)學(xué)研究中,要實(shí)現(xiàn)對物種有效的管理與保護(hù)首先要對物種分布有清晰而明確的認(rèn)識,這需要建立準(zhǔn)確、快速的物種監(jiān)測機(jī)制。近年來,環(huán)境DNA metebarcoding被廣泛地應(yīng)用于生態(tài)學(xué)的研究中。環(huán)境DNA metabarcoding是在DNA條形碼技術(shù)的基礎(chǔ)上研究出來的新興方法。該方法從環(huán)境樣品中直接提取DNA,當(dāng)結(jié)合第二代測序技術(shù)時可以實(shí)現(xiàn)對環(huán)境樣品中所有存在物種的有效檢測,在生物多樣性研究、物種監(jiān)測、動物食性分析和外來入侵物種檢測等方面都有較好的應(yīng)用前景。

    1 環(huán)境DNA metabarcoding

    1.1環(huán)境DNA

    環(huán)境DNA(environmental DNA, eDNA)的概念早在20世紀(jì)80年代末就已經(jīng)被提出,最初僅應(yīng)用于微生物研究領(lǐng)域,它是指從環(huán)境樣品(如土壤、沉積物、空氣、水體等)中提取的DNA,不對任何目標(biāo)生物進(jìn)行分離[1]。受該技術(shù)在微生物研究方面的啟發(fā),科學(xué)家們研究發(fā)現(xiàn)可以利用環(huán)境樣品中包含的DNA進(jìn)行物種鑒定及植食性動物的食性分析等方面的研究。

    環(huán)境樣品中通常含有動植物脫落的細(xì)胞或者是游離的DNA,生物排泄過程中脫落的皮膚,尿液和糞便等也可以是“環(huán)境DNA”的來源,并能提供一個環(huán)境系統(tǒng)中物種的存在記錄[2-3]。

    使用傳統(tǒng)采樣方法開展生態(tài)學(xué)研究,一些物種的樣本采集通常比較困難,并且耗時費(fèi)力。環(huán)境DNA允許從土壤、水體或是空氣中包含的動物毛發(fā)、脫落的細(xì)胞、糞便等進(jìn)行非損傷性取樣,從而成功進(jìn)行遺傳檢測。相比傳統(tǒng)的直接采樣,該方法成本更低,樣本采集受氣候條件影響較小[4],可以進(jìn)行更大數(shù)量的樣本收集。環(huán)境DNA還可以用來分析那些很難通過傳統(tǒng)方法監(jiān)測的物種,比如那些需要經(jīng)過相關(guān)部門許可才可以獲得的瀕?;蛘湎∥锓N[5]。目前,有關(guān)環(huán)境DNA在生態(tài)學(xué)中的應(yīng)用研究已有較多報(bào)道[6-12]。

    1.2環(huán)境DNA metabarcoding

    DNA 條形碼作為物種快速鑒定和發(fā)現(xiàn)新物種的方法已得到廣泛應(yīng)用。對于不同類別生物的DNA條形碼已確定了不同的擴(kuò)增片段,例如:線粒體細(xì)胞色素C氧化酶I 基因(COI) 已經(jīng)作為標(biāo)準(zhǔn)條形碼廣泛用于鑒定動物物種[13]; 核糖體16S rRNA 常用于細(xì)菌的鑒定[14]; 而ITS 基因則作為真菌的標(biāo)準(zhǔn)條形碼用于物種鑒定[15]。中國植物條形碼工作組建議ITS 基因(或ITS2)應(yīng)作為種子植物的核心條形碼之一[16]。目前通常采用rbcL+matK+ITS 基因組合作為植物的標(biāo)準(zhǔn)DNA 條形碼用于物種鑒定[17]。

    DNA條形碼技術(shù)雖然發(fā)展前景廣闊,但由于DNA 易發(fā)生降解,所以對于保存年代久遠(yuǎn)的博物館館藏標(biāo)本,要獲得標(biāo)準(zhǔn)的DNA 條形碼序列比較困難。對于DNA降解的標(biāo)本,相對而言,200bp左右的序列易于擴(kuò)增。由此,DNA 微型條形碼(DNA minibarcode)的概念被提出,它是指長度為100—200 bp的DNA 序列[18]。但在實(shí)際研究中,一方面,DNA易降解;另一方面,傳統(tǒng)測序技術(shù)僅適用于小樣品量的測序,當(dāng)需要研究的混合樣品包含上千個個體時,傳統(tǒng)測序技術(shù)不能滿足復(fù)雜的、大范圍的樣品測序。在第二代測序技術(shù)的發(fā)展和DNA條形碼存在局限性的雙重作用下,metabarcoding應(yīng)運(yùn)而生。Hajibabaei[19]認(rèn)為這兩種技術(shù)的結(jié)合極大地推動了“DNA 宏分類學(xué)(DNA metasystematics)”的發(fā)展。環(huán)境DNA metabarcoding是指利用宏條形碼對環(huán)境樣本(如土壤、水、糞便等)中分離的DNA進(jìn)行擴(kuò)增和高通量測序,從而達(dá)到同時對環(huán)境樣本中多個物種(或高級分類單元)鑒定的目的[1,20]。環(huán)境DNA metabarcoding克服傳統(tǒng)物種鑒定方法耗時費(fèi)力的缺點(diǎn),以更高的效率和更低的成本極大地拓展DNA 條形碼在生態(tài)學(xué)中的應(yīng)用。

    環(huán)境DNA metabarcoding整合DNA 條形碼和高通量測序技術(shù),通過提取環(huán)境樣品中的DNA,并使用特異性引物進(jìn)行PCR擴(kuò)增,對擴(kuò)增產(chǎn)物進(jìn)行測序后得到的可操縱分類單元(operational taxonomic units,OTUs)進(jìn)行物種鑒定[21]。該技術(shù)最大的優(yōu)勢在于高通量、低成本并能快速地鑒定物種[4- 5],最大的局限性在于PCR 的偏向性,導(dǎo)致條形碼的解析度和普適性及數(shù)據(jù)庫的完善度水平較低。要實(shí)現(xiàn)metabarcoding 技術(shù)高通量、準(zhǔn)確、快速物種鑒定,最為關(guān)鍵的是構(gòu)建完善的物種DNA 條形碼的標(biāo)準(zhǔn)數(shù)據(jù)庫、物種信息庫、信息共享和應(yīng)用平臺[22]。

    環(huán)境DNA metabarcoding擴(kuò)大了生物多樣性的研究范圍,其快速、可重復(fù)及高效性的特點(diǎn)讓更多的研究人員進(jìn)行物種鑒定和生物多樣性評估。在對混合樣品(水流樣品或土壤樣品等)的生物多樣性評估時,傳統(tǒng)方法中所遇的問題都會迎刃而解。

    2 環(huán)境DNA metabarcoding的研究方法

    2.1樣品采集

    目前,在生態(tài)學(xué)研究中環(huán)境DNA樣本的主要來源是水體、土壤和糞便等。水樣可以來源于實(shí)驗(yàn)室的水箱[23]、天然的池塘[24-25]、河流[26- 27]、溪流[8]、海水[6,28]等。池塘、溪流、河流等的采樣量從15mL到10L不等,對于溪流、河水和海水而言,一般標(biāo)準(zhǔn)的樣本量大小是1—2L[5]。在物種監(jiān)測中一些研究小組采取在水體的不同地點(diǎn)3次取樣的方法來提高水樣中的物種覆蓋度,增加物種檢測概率[7- 8,10,24]。對于流動水,應(yīng)在水流上游連續(xù)采樣以避免水樣的重復(fù)采集。采樣過程中始終佩戴一次性無菌手套,每采樣一次及時更換手套[29]。水樣的采集通常采用沉淀法(precipitation)或?yàn)V膜法(filtration)[10,24]。沉淀法是使用50mL的滅菌離心管,取15mL的水樣,加入1.5mL的3mol/L醋酸鈉和33.5mL的無水乙醇,放置于-20℃冰箱保存[24,30]。同沉淀法相比,濾膜法需要更多體積的水, 其優(yōu)點(diǎn)是可以通過大量的水收集DNA,但也會丟失一些能夠增加檢測幾率的可溶性DNA[31]。水通過裝有47mm直徑,0.45m孔徑的硝酸纖維素?zé)o菌過濾膜的一次性過濾器過濾。過濾后,將濾膜放置在無菌的2mL小瓶中,加入無水乙醇放置于-20℃保存[32]。對于濾膜規(guī)格和類型的選擇,并沒有統(tǒng)一的標(biāo)準(zhǔn),有些研究者選用25mm直徑,0.7m的玻璃纖維濾膜或1.5m的玻璃纖維素濾膜[27,33]。研究發(fā)現(xiàn)對于靜水和流水樣本的大型無脊椎動物,使用濾膜法檢測率最高。環(huán)境DNA采樣同傳統(tǒng)的方法比更容易以統(tǒng)一的標(biāo)準(zhǔn)實(shí)施,而且不易受到采樣強(qiáng)度和人員培訓(xùn)差異的影響而變化[29]。

    在對土壤的生物多樣性分析中,土壤樣品應(yīng)該盡可能地代表研究區(qū)域的生物多樣性,當(dāng)研究區(qū)域的生物分布具有異質(zhì)性時需要多份樣品混合取樣。在收集土壤樣品時需要注意兩點(diǎn):一是應(yīng)該隨機(jī)或有規(guī)律地采集來自不同深度的土樣,二是每個取樣地點(diǎn)應(yīng)該至少收集兩份樣品,同時每份樣品應(yīng)該是該取樣地點(diǎn)不同深度的土樣的混合品[31]。開展食性分析研究時,對采集到的1個或多個糞便樣本,應(yīng)從樣本混合物中進(jìn)行多次取樣,以增加檢測出樣品中所含被捕食動物的可能性[34]。

    2.2環(huán)境DNA的提取

    環(huán)境DNA提取方法也有很多種,包括氯仿萃取法[35]、 細(xì)胞溶解的物理破碎法[26]以及二氧化硅提取法[8]。Dveiner等[33]采集湖泊和河流水樣,比較了6種不同的樣本采集和DNA提取方法組合對水生態(tài)系統(tǒng)生物多樣性監(jiān)測的影響,結(jié)果表明不同的環(huán)境DNA采集和提取方法對DNA產(chǎn)量和高通量測序技術(shù)獲得的序列數(shù)量有顯著影響。采用濾膜法結(jié)合DNeasy 血液和組織試劑盒(Qiagen, 德國)提取DNA的方法對生活在激流環(huán)境中物種的檢測率最高。這也在檢測稀有物種或入侵物種的研究中得到證實(shí)[8,27,36]。目前研究中環(huán)境DNA提取使用較多的是DNeasy 血液及組織試劑盒(Qiagen, 德國), QIAamp 組織試劑盒(Qiagen德國), PowerWaterTMDNA 提取試劑盒 (MoBio Laboratories, 美國), PowerMaxTM土壤 DNA 提取試劑盒 (MoBio Laboratories, 美國)等[11,25,33]。

    2.3PCR擴(kuò)增與引物設(shè)計(jì)

    由于環(huán)境DNA樣本往往會混有不同物種,無法使用傳統(tǒng)的DNA條形碼進(jìn)行研究。根據(jù)研究目的的需要,通常會擴(kuò)增1個或多個基因。動物研究一般選擇擴(kuò)增線粒體COI基因,Cytb基因,12S rRNA和18S rDNA等,植物研究一般選擇擴(kuò)增葉綠體DNAtrnL 的P6環(huán)。Thomsen等[10]通過擴(kuò)增環(huán)境水樣中混合DNA的線粒體COI基因和Cytb基因,成功鑒別了該環(huán)境中包含的兩棲類、魚類、甲蟲、哺乳動物和昆蟲等。

    Metabarcoding技術(shù)依賴于設(shè)計(jì)理想的引物,通過PCR從環(huán)境樣品中擴(kuò)增出高質(zhì)量的混合DNA片段。由于環(huán)境中DNA易于受到紫外線和微生物分解作用的影響,物種的長DNA片段十分容易降解,因此使用短的片段(約90—120bp)十分重要[5]。在不同的研究中,選擇的宏條形碼不可能是完全相同的,在很多情況下需要針對特定的生物設(shè)計(jì)特異性的引物[38]。所以使用metabarcoding技術(shù)分析環(huán)境DNA的研究中,為了避免在不同分類種群中放大偏差,設(shè)計(jì)短小的、高度保守的引物是極其重要的。如果實(shí)驗(yàn)設(shè)計(jì)的引物保守度低和特異性不足的話,容易引起堿基錯配甚至引起實(shí)驗(yàn)結(jié)果出現(xiàn)假陽性和假陰性的結(jié)果[27]。因此理想的宏條形碼應(yīng)該具有序列短且高度保守、高辨別度等優(yōu)點(diǎn)。Leray等[37]設(shè)計(jì)了5組針對線粒體COI基因片段擴(kuò)增的特異性引物,研究天竺鯛科(Apogonidae)和金鱗魚科(Holocentridae)珊瑚魚的胃內(nèi)容物。結(jié)果表明其中一對引物(mlCOIintF/jgHCO2198)擴(kuò)增出的313bp的線粒體COI基因片段比傳統(tǒng)的COI基因通用引物(LCO1490/HCO2198)在后生動物多樣性研究中能夠獲得更多的OTUs,擴(kuò)增成功率更高。

    目前雖然有一些工具可以幫助生物學(xué)家設(shè)計(jì)引物,比如Primer5和QPrimer,但是還不能完全滿足metabarcoding設(shè)計(jì)引物的需求。Riaz等[38]通過研究證實(shí)ecoPrimers軟件可以滿足metabarcoding研究中設(shè)計(jì)引物的要求,盡管仍然存在一對引物覆蓋率較低的問題,但是可以通過設(shè)計(jì)極短的條形碼來提高引物的特異性和對目標(biāo)物種的覆蓋率。相信隨著公共數(shù)據(jù)庫中數(shù)據(jù)的不斷完善,可以根據(jù)生物的特異性和研究需求,精確地設(shè)計(jì)出更多更有效的條形碼。

    2.4第二代高通量測序技術(shù)

    1977年,Sanger等發(fā)明了雙脫氧核苷酸末端終止法,而Maxam和Gilbert發(fā)明的化學(xué)降解法標(biāo)志了第一代測序技術(shù)的產(chǎn)生。第一代測序技術(shù)完成了從噬菌體基因組到人類基因組草圖的測序工作,但是由于其測序成本高、速度慢等缺點(diǎn)并不是最理想的方法[39]。隨后,出現(xiàn)了可以邊合成邊測序,通過捕捉新合成的末端的標(biāo)記來確定DNA 的序列的第二代高通量測序技術(shù)(Next-generation sequencing,NGS),包括羅氏公司的454技術(shù)的GS FLX測序平臺、Illumina 公司的Solexa Genome Analyzer測序平臺和ABI公司的SOLiD測序平臺[40]。第二代測序技術(shù)的大致步驟為:將片斷化的基因組DNA 兩端連上特殊的接頭,通過不同的方式將每個片段結(jié)合在微珠或芯片上,形成幾百萬個可以同時進(jìn)行反應(yīng)的微型反應(yīng)池,每個單鏈DNA 片段在其微型反應(yīng)池中通過PCR 擴(kuò)增產(chǎn)生大量拷貝形成單克隆“DNA 簇群” (Polony),以作為測序的模板;再通過酶延伸反應(yīng)或寡核苷酸連接反應(yīng)同時對這幾百萬個微型反應(yīng)池中的模板進(jìn)行測序[41]。

    第二代測序技術(shù)最顯著的特征就是高通量,一次能對幾十萬到幾百萬條DNA分子進(jìn)行序列測序,目前SOLiD 3系統(tǒng)單次運(yùn)行可產(chǎn)生50GB的序列數(shù)據(jù),相當(dāng)于17倍人類基因組覆蓋度[42]。與第一代測序技術(shù)相比,雖然讀取的序列片段長度和準(zhǔn)確度不及Sanger測序,但是第二代測序技術(shù)不僅降低了測序成本并縮短了測序時間,同時又可以對復(fù)雜的混合樣品進(jìn)行測序。因此第二代測序技術(shù)的廣泛應(yīng)用可以在短時間內(nèi)對環(huán)境樣品中包含的上百萬條序列進(jìn)行測序,可以實(shí)現(xiàn)短時間對復(fù)雜環(huán)境樣品的研究。隨著第三代測序技術(shù)(如SMRT技術(shù),納米孔單分子測序技術(shù)等)的迅速發(fā)展,測序的長度和準(zhǔn)確性都會大幅度提高。同時第三代測序技術(shù)無需進(jìn)行PCR擴(kuò)增[43- 44],解決了第二代測序錯誤偏向性的問題。

    2.5生物信息學(xué)分析

    Metabarcoding使用通用PCR引物擴(kuò)增環(huán)境中的DNA,經(jīng)高通量測序后得到其中包含的大規(guī)模分類學(xué)信息基因。然而如何對測序得到的大量數(shù)據(jù)進(jìn)行處理是第二代測序的主要問題。由于得到的序列片段較短,而且環(huán)境樣品通常具有高度的復(fù)雜性,很難確定某條序列來自于哪個物種,因此面對海量的短序列,拼接和進(jìn)一步的分析具有相當(dāng)大的挑戰(zhàn)。一般來說,數(shù)據(jù)處理過程主要包括:(1)去除噪音,即對測序數(shù)據(jù)進(jìn)行預(yù)處理以去除低質(zhì)量的污染的序列。測序產(chǎn)物可能含有非特異性擴(kuò)增片段,利用特異性引物信息可以將其去除;序列中可能含有模糊堿基(ambiguous)、單堿基高重復(fù)區(qū)(homologous),長度過短的序列,及PCR過程中產(chǎn)生的一些嵌合體(chimera),將這些序列納入分析范圍會降低分析質(zhì)量,因此修剪、去除(trim)這部分序列,可得到供精準(zhǔn)分析的優(yōu)化序列。如:SeqClean、VecScree、UCHIME等可以用來去除噪音[39];(2)使用相關(guān)軟件如(USEARCH/CROP、Denoiser/UCLUST、OCTUPUS)根據(jù)相似度或序列組成對序列進(jìn)行歸類獲得OTUs。因?yàn)槊總€物種單獨(dú)貢獻(xiàn)許多DNA,而且通過PCR被復(fù)制了很多次,所以輸出的數(shù)據(jù)集需要使用軟件分配OTUs,其每一個理想的OTUs應(yīng)該只包含來自一個物種的序列[21]。最后,代表序列取自每個OTU和分配使用一個上面列出的數(shù)據(jù)庫的一個或多個分類法。通過多個對條形碼的分析可以獲得在該環(huán)境中各種物種的存在情況;(3)在獲得了OTUs后可以針對研究的具體內(nèi)容進(jìn)行后續(xù)分析,比如α-多樣性分析、β-多樣性分析、構(gòu)建系統(tǒng)發(fā)育樹、物種分類及物種豐度分析等[39]。

    3 環(huán)境DNA metabarcoding在生態(tài)學(xué)研究中的應(yīng)用

    3.1生物多樣性研究

    物種的準(zhǔn)確鑒定是探索、認(rèn)識物種和保護(hù)生物多樣性的基礎(chǔ)。但是,使用傳統(tǒng)的方法鑒定物種有很多局限性:(1)耗時較長,效率較低;(2)過于依賴專家的經(jīng)驗(yàn),鑒定與分類結(jié)果的正確性有待商榷;(3)在某些物種的特定階段或發(fā)育階段時,會影響鑒定結(jié)果。DNA 條形碼可以構(gòu)建植物與動物之間的網(wǎng)狀拓?fù)浣Y(jié)構(gòu),可以從群落水平探討物種多樣性、系統(tǒng)發(fā)育多樣性與群落穩(wěn)定性之間的相互關(guān)系[17]。

    近年來,環(huán)境DNA metabarcoding已成功應(yīng)用于動植物的多樣性研究。Yoccoz等[12]利用metabarcoding技術(shù)分析了不同植被類型的土壤DNA,獲得其中葉綠體trnL(UAA)內(nèi)含子P6環(huán)序列,結(jié)果表明利用土壤DNA檢測的植物多樣性與傳統(tǒng)地面調(diào)查估計(jì)的植物種類和生長型多樣性結(jié)果高度一致。Bienert等[45]利用metabarcoding技術(shù)結(jié)合線粒體16S基因,設(shè)計(jì)兩個DNA短序列的特異性引物(大約30bp和70bp)進(jìn)行PCR擴(kuò)增及測序,分析了法國阿爾卑斯山地生態(tài)系統(tǒng)土壤中蚯蚓組成。研究表明基于DNA方法分析的蚯蚓群落多樣性與傳統(tǒng)分類方法獲得的結(jié)果一致,闡明了環(huán)境DNA作為土壤動物多樣性評價(jià)工具的潛力。Thomsen等[6]從丹麥的一個海洋生態(tài)系統(tǒng)中采集海水,利用環(huán)境DNA metabarcoding技術(shù)研究海洋魚類生物多樣性。在獲得的環(huán)境DNA中發(fā)現(xiàn)了15種不同的魚類以及4種鳥類,其中包括一些用傳統(tǒng)監(jiān)測方法很難發(fā)現(xiàn)的稀有物種。同時與9種傳統(tǒng)的海洋魚類調(diào)查方法比較,環(huán)境DNA檢測出的魚類多樣性與其他傳統(tǒng)方法相當(dāng),甚至優(yōu)于傳統(tǒng)方法。Thomsen等[10]研究表明從湖泊、池塘和河流的少量水樣中直接獲得DNA,可以檢測并定量分析研究區(qū)域內(nèi)珍稀及瀕危淡水動物(兩棲類、魚類、哺乳類、昆蟲類和甲殼類)多樣性。同時表明通過池塘水中分離的DNA,利用第二代測序技術(shù)可以檢測到研究區(qū)域內(nèi)兩棲類和魚類群體,環(huán)境DNA可以作為監(jiān)測珍稀瀕危物種的工具,應(yīng)用于更廣泛的物種監(jiān)測中。

    應(yīng)用環(huán)境DNA metabarcoding技術(shù)不僅可以針對當(dāng)前的生物多樣性進(jìn)行有效的評估,也可以對過去的動植物群落的多樣性進(jìn)行評估[11,45]。Epp等[11]通過對過去和現(xiàn)在的環(huán)境DNA進(jìn)行metabarcoding研究,所有目標(biāo)類群都能通過metabarcoding的方法檢測出來,但是類群之間和不同年代的環(huán)境樣品所檢測的結(jié)果存在較大變異范圍,如更新世時期的樣品中真菌檢測成功率最高,同時也可以檢測到苔蘚植物、甲蟲和鳥類的序列。

    目前國內(nèi)在應(yīng)用環(huán)境DNA metabarcoding技術(shù)研究生物多樣性的報(bào)道還比較少。楊晨雪等[46]使用metabarcoding技術(shù)研究土壤樣品中的生物多樣性,通過三種不同的生物信息學(xué)分析方法進(jìn)行比較分析,探討了用18S rDNA基因擴(kuò)增土壤中后生動物的有效性,研究證明metabarcoding技術(shù)可以快速、全面、有力地調(diào)查土壤動物多樣性。中科院昆明動物所Yu等[47]利用馬氏網(wǎng)誘捕節(jié)肢動物,對整個混合樣本的DNA進(jìn)行線粒體COI基因擴(kuò)增,再進(jìn)行高通量測序,并結(jié)合生物信息學(xué)分析,探討metabarcoding方法對群落的差異度(β-多樣性)和群落內(nèi)系統(tǒng)發(fā)育(α-多樣性)評估的準(zhǔn)確性,結(jié)果表明metabarcoding方法可以在更廣的尺度上有效地測量生物多樣性。

    3.2食性分析

    食性分析及食物網(wǎng)絡(luò)的研究可直接體現(xiàn)出生態(tài)群落的功能結(jié)構(gòu),是生態(tài)系統(tǒng)研究的重要組成部分[34]。傳統(tǒng)的食性分析通常通過收集糞便并鑒別其中的食物殘留,研究人員通過肉眼直接觀察或利用早期的分子實(shí)驗(yàn)獲取食譜信息。由于食物殘留情況各異,這種方法耗時費(fèi)力,對分析人員的經(jīng)驗(yàn)依賴性大而且解析度不高,時間長,準(zhǔn)確度不高[48]。

    在動植物間的植食網(wǎng)絡(luò)關(guān)系中,當(dāng)觀察動物的行為存在困難時,可利用動物消化道中的食物殘?jiān)?,利用metabarcoding技術(shù),研究人員可快速準(zhǔn)確地進(jìn)行大量樣品的食性分析[34]。Shehzad等[49]基于第二代測序技術(shù)進(jìn)行了食肉動物豹貓(Prionailurusbengalensis)的食性分析。從38個豹貓的糞便提取出了混合DNA,通過擴(kuò)增約100bp片段的12S rRNA研究豹貓取食的食物種類,發(fā)現(xiàn)了18種獵物,包括八種哺乳動物、八種鳥類、一種兩棲類和一種魚類,并且確定豹貓主要食用嚙齒動物鼠類。

    3.3外來入侵物種的檢測

    目前外來入侵物種的情況日趨嚴(yán)重。當(dāng)外來物種入侵時,入侵物種很有可能破壞原本的生態(tài)平衡,使生態(tài)系統(tǒng)不完善,從而對本地物種的種類和數(shù)量產(chǎn)生影響[24]。但是在早期入侵階段外來入侵物種通常數(shù)量不多,使用傳統(tǒng)的檢測方法耗時且效率不高,難度高,不能及時被發(fā)現(xiàn)。這會導(dǎo)致外來入侵物種對本地生態(tài)平衡的危害達(dá)到最高值時才受到重視,為時已晚。環(huán)境DNA metabarcoding增強(qiáng)了對入侵物種早期監(jiān)測的可能性,可以有效地遏制物種入侵。

    研究發(fā)現(xiàn)利用環(huán)境DNA metabarcoding可以從淡水中提取環(huán)境DNA,檢測出外來入侵物種的存在。Ficetola等[24]最早將環(huán)境DNA metabarcoding應(yīng)用于物種監(jiān)測,他們使用特異性引物擴(kuò)增短線粒體DNA序列追蹤在受控環(huán)境和自然環(huán)境中一種入侵牛蛙(Ranacatesbeiana)的存在,發(fā)現(xiàn)對環(huán)境中低密度的物種,使用metabarcoding技術(shù)仍然可以進(jìn)行快速有效的物種檢測。這是首次結(jié)合大規(guī)模測序和DNA條形碼技術(shù)的發(fā)展來進(jìn)行物種識別,為利用環(huán)境DNA監(jiān)測外來物種提供了新方法和新思路。Jerde等[26]利用環(huán)境DNA研究美國芝加哥地區(qū)河流中的兩種入侵亞洲鯉科魚(Hypophthalmichthysmolitrix和H.nobilis),發(fā)現(xiàn)使用傳統(tǒng)電氣捕魚法檢測到一條魚需要93人d的努力,而環(huán)境DNA metabarcoding方法檢測到目標(biāo)物種只需要0.174人d。在亞洲鯉科魚種群密度很低的區(qū)域,只有環(huán)境DNA方法能夠檢測到這兩種魚[5]??梢?,環(huán)境DNA metabarcoding方法不僅在時間和資金成本上優(yōu)于傳統(tǒng)監(jiān)測方法,而且可以在種群密度很低的環(huán)境樣本中靈敏地檢測到入侵物種的存在。Piaggio等[23]捕獲入侵佛羅里達(dá)州的半水棲物種緬甸巨蟒(Pythonbivittatus),放置于90L的水箱中評估水中環(huán)境DNA的降解速率,結(jié)果表明緬甸巨蟒的DNA在96h以內(nèi)都可以檢測到。在佛羅里達(dá)州6個野外地點(diǎn)采集的水樣中,有5個地點(diǎn)的緬甸巨蟒環(huán)境DNA為陽性,這與緬甸巨蟒的野外分布結(jié)果一致。Takahara等[7]使用環(huán)境DNA方法采集日本海岸以及周圍島嶼的70個水樣對一種入侵藍(lán)腮太陽魚(Lepomismacrochirus)進(jìn)行檢測。利用藍(lán)腮太陽魚的特異性引物(100bp的Cytb片段)擴(kuò)增水環(huán)境DNA,結(jié)果表明70個池塘中有19個池塘遭受藍(lán)腮太陽魚的入侵。盡管魚的大小、習(xí)性等因素對實(shí)驗(yàn)結(jié)果有一定影響,但使用環(huán)境DNA metabacoding提高了監(jiān)測入侵物種的準(zhǔn)確度。

    3.4物種監(jiān)測

    通過水體中的環(huán)境DNA可以有效檢測水生脊椎動物,并且已在濕地和大型河流研究中得到證實(shí)[26]。利用環(huán)境DNA檢測河流中的稀有物種和隱秘物種可以增加調(diào)查結(jié)果準(zhǔn)確性,減少調(diào)查成本。Goldberg等[8]利用metabarcoding技術(shù),通過收集快速流動的水體環(huán)境DNA樣本,對尾突角蟾(Ascaphusmontanus)和愛達(dá)荷陸巨螈(Dicamptodonaterrimus)兩種珍稀物種進(jìn)行監(jiān)測。設(shè)計(jì)特異性引物尾突角蟾和愛達(dá)荷陸巨螈的線粒體Cytb區(qū)域的85bp和78bp片段,在5個不同物種密度的水樣中均成功進(jìn)行了PCR擴(kuò)增并靈敏地檢測到了目標(biāo)物種的存在。該研究還發(fā)現(xiàn)早春季節(jié)比早秋季節(jié)更難以檢測到兩棲類物種存在,推測可能是由于早春比早秋溫度更低,降低了生物代謝速率而導(dǎo)致。Foote等[28]從海水中分離DNA,使用特異性引物擴(kuò)增線粒體DNA區(qū)域60—80bp的片段檢測水體中是否存在港灣鼠海豚(Phocoenaphocoena),對控制條件下和自然條件下該物種進(jìn)行遺傳監(jiān)測。在一個樣地中檢測到長肢領(lǐng)航鯨(Globicephalamelas),而該物種極少在研究海域被發(fā)現(xiàn),結(jié)果表明環(huán)境DNA的方法可用于檢測稀有物種。Wilcox等[27]利用環(huán)境DNA結(jié)合定量PCR檢測了河流中的兩個稀有物種,美洲紅點(diǎn)鮭(Salvelinusfontinalis)和強(qiáng)壯紅點(diǎn)鮭(S.confluentus),同時發(fā)現(xiàn)定量PCR(qPCR)比傳統(tǒng)PCR更加靈敏,檢測概率更高。Rees等[25]利用環(huán)境DNA和傳統(tǒng)的實(shí)地調(diào)查方法開展對英國珍稀保護(hù)物種冠北螈(Trituruscristatus)的監(jiān)測,野外調(diào)查發(fā)現(xiàn)有冠北螈分布的地點(diǎn)中,采集的水樣DNA中成功檢測出冠北螈的概率為84%。使用環(huán)境DNA方法監(jiān)測物種,對物種及其生存環(huán)境不會造成損害。隨著測序成本的下降,環(huán)境DNA與metabarcoding技術(shù)結(jié)合應(yīng)用在許多情況下甚至優(yōu)于一些傳統(tǒng)監(jiān)測方法,成為物種監(jiān)測方法的又一選擇[4]。

    4 面臨的挑戰(zhàn)

    4.1樣品來源的影響

    使用環(huán)境DNA監(jiān)測水生生物時,必須保證能獲得該環(huán)境的完整樣品并能進(jìn)行很好的保存[27]。水體的溫度和光照情況[53-54],水體的流速[55],物種密度和在環(huán)境中的停留時間[54],生物體的發(fā)育階段[55]等均會對環(huán)境中所取的DNA的量造成影響。Maruyama等[56]對水體中DNA的半衰期進(jìn)行了研究和計(jì)算,發(fā)現(xiàn)魚類離開此水體后環(huán)境DNA的半衰期是6.3h,但對環(huán)境DNA降解速率及其影響因素,環(huán)境DNA在環(huán)境中存在時間的長短等問題仍缺乏充分的了解。

    在采集溪流水樣時,由于許多溪流水來自于地下水,所以接觸到目標(biāo)生物脫落DNA的可能性較低。此外,由于溪流水水位淺而且是充分混合的,因此脫落的DNA被光降解的可能性更大,導(dǎo)致在溪流中檢測到的目標(biāo)物種密度會比較低[27]。使用環(huán)境DNA metabarcoding技術(shù)對淡水中的物種進(jìn)行檢測已被證實(shí)快速有效。淡水中水流較平和且流量小,所以準(zhǔn)確度高。與淡水相比,海洋生態(tài)系統(tǒng)中有單位體積中的生物量比例更大,海流和波浪的影響作用,以及鹽度對保存和提取環(huán)境DNA的影響都需要加以考慮。這些因素可能意味著海水中的環(huán)境DNA不太集中,且分散更迅速,提取DNA時比較困難,導(dǎo)致在海水中使用時只能針對一小片區(qū)域,并且對外來入侵物種的檢測效率不高。因此在海水中的使用還需要進(jìn)一步的研究與完善[28]。

    4.2PCR擴(kuò)增的偏向性

    環(huán)境DNA宏條形碼技術(shù)的一個潛在的局限性是PCR引物設(shè)計(jì)。通用引物的使用是普遍的,對所有的序列都有擴(kuò)增的作用。非特異性的引物和探針會造成假陽性的結(jié)果,非靶標(biāo)模板的競爭會產(chǎn)生不準(zhǔn)確的DNA量估計(jì),從而導(dǎo)致物種豐富度的不準(zhǔn)確估計(jì)[53]??梢酝ㄟ^以下兩種方式來解決由于非特異性引物而造成的假陽性:(1)設(shè)計(jì)分析與非目標(biāo)物種最大化堿基對錯配的引物區(qū)域;(2)通過實(shí)驗(yàn)來測定靶DNA和非靶DNA的純樣品和混合樣品[27]。

    從長遠(yuǎn)角度看,在PCR擴(kuò)增時低密度的物種可能會因?yàn)橥ㄓ靡锏脑虮缓鲆暎@也就是PCR的偏向性。但是如果任由PCR的偏向性發(fā)生,會導(dǎo)致一些稀有的或者本身難以擴(kuò)增的物種有數(shù)據(jù)丟失的情況,或者引入嵌合體等錯誤產(chǎn)物[21]。目前的第三代測序技術(shù)就是采用單分子讀取技術(shù),有著更快的數(shù)據(jù)讀取速度。而且不依賴于PCR擴(kuò)增,不僅進(jìn)一步降低了測序成本,而且解決了在環(huán)境DNA metabarcoding的研究中PCR擴(kuò)增偏向的制約問題[43]。

    5 展望

    在第三代測序技術(shù)的推動下,不需要經(jīng)過PCR步驟的metabarcoding(PCR-free metabarcoding)技術(shù),相信在不久的將來就會實(shí)現(xiàn)。該技術(shù)的應(yīng)用可以對物種中存在線粒體DNA的總量進(jìn)行推斷,還可以在生態(tài)學(xué)的研究中作為一種新興方法對物種的生物量和相對豐度進(jìn)行評估,將極大地改變生態(tài)學(xué)家對生態(tài)多樣性檢測方法的認(rèn)知。所以在今后環(huán)境DNAmetabarcoding研究中, 一方面要大力發(fā)展PCR-free metabarcoding,不斷完善數(shù)據(jù)庫并且發(fā)展生物信息學(xué)工具增加數(shù)據(jù)處理的可信度[21]。另一方面應(yīng)盡快形成環(huán)境DNA metabarcoding對不同環(huán)境采樣的統(tǒng)一標(biāo)準(zhǔn)。雖然環(huán)境DNA metabarcoding并不能完全取代傳統(tǒng)的野外工作者和專業(yè)分類學(xué)家的工作,但是對于提高受限制和保護(hù)資源的有效利用而言是一個十分有效的工具[5]。

    [1]Taberlet P, Coissac E, Hajibabaei M, Rieseberg L H. Environmental DNA. Molecular Ecology, 2012, 21(8): 1789- 1793.

    [2]Haile J, Froese D G, Macphee R D, Roberts R G, Arnold L J, Reyes A V, Rasmussen M, Nielsen R, Brook B W, Robinson S, Demuro M, Gilbert M T, Munch K, Austin J J, Cooper A, Barnes I, M ller P, Willerslev E. Ancient DNA reveals late survival of mammoth and horse in interior Alaska. Proceedings of the National Academy of Sciences of the United States of America, 2009, 106(52): 22352- 22357.

    [3]Lydolph M C, Jacobsen J, Arctander P, Gilbert M T, Gilichinsky D A, Hansen A J, Willerslev E, Lange L. Beringian paleoecology inferred from permafrost-preserved fungal DNA. Applied and Environmental Microbiology, 2005, 71(2): 1012- 1017.

    [4]Thomsen P F, Willerslev E. Environmental DNA - An emerging tool in conservation for monitoring past and present biodiversity. Biological Conservation, 2015, 183: 4- 18.

    [5]Rees H C, Maddison B C, Middleditch D J, Patmore J R M, Gough K C. REVIEW: The detection of aquatic animal species using environmental DNA-a review of eDNA as a survey tool in ecology. Journal of Applied Ecology, 2014, 51(5): 1450- 1459.

    [6]Thomsen P F, Kielgast J, Iversen L L, Moller P R, Rasmussen M, Willerslev E. Detection of a diverse marine fish fauna using environmental DNA from seawater samples. PLoS One, 2012, 7(8): e41732.

    [7]Takahara T, Minamoto T, Doi H. Using environmental DNA to estimate the distribution of an invasive fish species in ponds. PLoS One, 2013, 8(2): e56584.

    [8]Goldberg C S, Pilliod D S, Arkle R S, Waits L P. Molecular detection of vertebrates in stream water: a demonstration using Rocky Mountain tailed frogs and Idaho giant salamanders. PLoS One, 2011, 6(7): e22746.

    [9]Mahon A R, Jerde C L, Galaska M, Bergner J L, Chadderton W L, Lodge D M, Hunter M E, Nico L G. Validation of eDNA surveillance sensitivity for detection of Asian carps in controlled and field experiments. PLoS One, 2013, 8(3): e58316.

    [10]Thomsen P F, Kielgast J, Iversen L L, Wiuf C, Rasmussen M, Gilbert M T, Orlando L, Willerslev E. Monitoring endangered freshwater biodiversity using environmental DNA. Molecular Ecology, 2012, 21(11): 2565- 2573.

    [11]Epp L S, Boessenkool S, Bellemain E P, Haile J, Esposito A, Riaz T, Erséus C, Gusarov V I, Edwards M E, Johnsen A, Sten?ien H K, Hassel K, Kauserud H, Yoccoz N G, Brathen K A, Willerslev E, Taberlet P, Coissac E, Brochmann C. New environmental metabarcodes for analysing soil DNA: potential for studying past and present ecosystems. Molecular Ecology, 2012, 21(8): 1821- 1833.

    [12]Yoccoz N G, Brathen K A, Gielly L, Haile J, Edwards M E, Goslar T, Von Stedingk H, Brysting A K, Coissac E, Pompanon F, S?nsteb? J H, Miquel C, Valentini A, De Bello F, Chave J, Thuiller W, Wincker P, Cruaud C, Gavory F, Rasmussen M, Gilbert M T, Orlando L, Brochmann C, Willerslev E, Taberlet P. DNA from soil mirrors plant taxonomic and growth form diversity. Molecular Ecology, 2012, 21(15): 3647- 3655.

    [13]Hebert P D N, Cywinska A, Ball S L, deWaard J R. Biological identifications through DNA barcodes. Proceeding of the Royal Society of London Series B: Biological Science, 2003, 270(1512): 313- 321.

    [14]Sogin M L, Morrison H G, Huber J A, Mark Welch D, Huse S M, Neal P R, Arrieta J M, Herndl G J. Microbial diversity in the deep sea and the underexplored "rare biosphere". Proceedings of the National Academy of Sciences of the United States of America, 2006, 103(32): 12115- 12120.

    [15]Lahaye R, van der Bank M, Bogarin D, Warner J, Pupulin F, Gigot G, Maurin O, Duthoit S, Barraclough T G, Savolainen V. DNA barcoding the floras of biodiversity hotspots. Proceeding of the National Academy of Sciences of the United States of America, 2008, 105(8): 2923- 2928.

    [16]Li D Z, Gao L M, Li H T, Wang H, Ge X J, Liu J Q, Chen Z D, Zhou S L, Chen S L, Yang J B, Fu C X, Zeng C X, Yan H F, Zhu Y J, Sun Y S, Chen S Y, Zhao L, Wang K, Yang T, Duan G W. Comparative analysis of a large dataset indicates that internal transcribed spacer (ITS) should be incorporated into the core barcode for seed plants. Proceeding of the National Academy of Sciences of the United States of America, 2011, 108(49): 19641- 19646.

    [17]羅亞皇, 劉杰, 高連明, 李德銖. DNA條形碼在生態(tài)學(xué)研究中的應(yīng)用與展望. 植物分類與資源學(xué)報(bào), 2013, 35(6): 761- 768.

    [18]Meusnier I, Singer G A C, Landry J F, Hickey D A, Hebert P D N, Hajibabaei M. A universal DNA mini-barcode for biodiversity analysis. BMC Genomics, 2008, 9: 214- 214.

    [19]Hajibabaei M. The golden age of DNA metasystematics. Trends in Genetics, 2012, 28(11): 535- 537.

    [20]Taberlet P, Coissac E, Pompanon F, Brochmann C, Willerslev E. Towards next-generation biodiversity assessment using DNA metabarcoding. Molecular Ecology, 2012, 21(8): 2045- 2050.

    [21]Ji Y Q, Ashton L, Pedley S M, Edwards D P, Tang Y, Nakamura A, Kitching R, Dolman P M, Woodcock P, Edwards F A, Larsen T H, Hsu W W, Benedick S, Hamer K C, Wilcove D S, Bruce C, Wang X Y, Levi T, Lott M, Emerson B C, Yu D W. Reliable, verifiable and efficient monitoring of biodiversity via metabarcoding. Ecology Letters, 2013, 16(10): 1245- 1257.

    [22]唐敏, 伊廷雙, 王欣, 譚美華, 周欣. Metabarcoding技術(shù)在植物鑒定和多樣性研究中的應(yīng)用. 植物分類與資源學(xué)報(bào), 2013, 35(6): 769- 773.

    [23]Piaggio A J, Engeman R M, Hopken M W, Humphrey J S, Keacher K L, Bruce W E, Avery M L. Detecting an elusive invasive species: a diagnostic PCR to detect Burmese python in Florida waters and an assessment of persistence of environmental DNA. Molecular Ecology Resources, 2014, 14(2): 374- 380.

    [24]Ficetola G F, Miaud C, Pompanon F, Taberlet P. Species detection using enviromental DNA from water samples. Biology Letters, 2008, 4(4): 423- 425.

    [25]Rees H C, Bishop K, Middleditch D J, Patmore J R, Maddison B C, Gough K C. The application of eDNA for monitoring of the Great Crested Newt in the UK. Ecology and Evolution, 2014, 4(21): 4023- 4032.

    [26]Jerde C L, Mahon A R, Chadderton W L, Lodge D M. “Sight-unseen” detection of rare aquatic species using environmental DNA. Conservation Letters, 2011, 4(2): 150- 157.

    [27]Wilcox T M, McKelvey K S, Young M K, Jane S F, Lowe W H, Whiteley A R, Schwartz M K. Robust detection of rare species using environmental DNA: the importance of primer specificity. PLoS One, 2013, 8(3): e59520.

    [28]Foote A D, Thomsen P F, Sveegaard S, Wahlberg M, Kielgast J, Kyhn L A, Salling A B, Galatius A, Orlando L, Gilbert M T P. Investigating the potential use of environmental DNA (eDNA) for genetic monitoring of marine mammals. PLoS One, 2012, 7(8): e41781.

    [29]Sigsgaard E E, Carl H, M?ller P R, Thomsen P F. Monitoring the near-extinct European weather loach in Denmark based on environmental DNA from water samples. Biological Conservation, 2015, 183: 46- 52.

    [30]Turner C R, Uy K L, Everhartb R C. Fish environmental DNA is more concentrated in aquatic sediments than surface water. Biological Conservation, 2015, 183: 93- 102.

    [31]Goldberg C S, Strickler K M, Pilliod D S. Moving environmental DNA methods from concept to practice for monitoring aquatic macroorganisms. Biological Conservation, 2015, 183: 1- 3.

    [32]Laramie M B, Pilliod D S, Goldberg C S. Characterizing the distribution of an endangered salmonid using environmental DNA analysis. Biological Conservation, 2015, 183: 29- 37.

    [33]Deiner K, Walser J C, M?chler E, Altermatt F. Choice of capture and extraction methods affect detection of freshwater biodiversity from environmental DNA. Biological Conservation, 2015, 183: 53- 63.

    [34]Pompanon F, Deagle B E, Symondson W O, Brown D S, Jarman S N, Taberlet P. Who is eating what: diet assessment using next generation sequencing. Molecular Ecology, 2012, 21(8): 1931- 1950.

    [35]Renshaw M A, Olds B P, Jerde C L, McVeigh M M, Lodge D M. The room temperature preservation of filtered environmental DNA samples and assimilation into a phenol-chloroform-isoamyl alcohol DNA extraction. Molecular Ecology Resources, 2015, 15(1): 168- 176.

    [36]Pilliod D S, Goldberg C S, Arkle R S, Waits L P. Estimating occupancy and abundance of stream amphibians using environmental DNA from filtered water samples. Canadian Journal of Fisheries and Aquatic Sciences, 2013, 70(8): 1123- 1130.

    [37]Leray M, Yang J Y, Meyer C P, Mills S C, Agudelo N, Ranwez V, Boehm J T, Machida R J. A new versatile primer set targeting a short fragment of the mitochondrial COI region for metabarcoding metazoan diversity: application for characterizing coral reef fish gut contents. Frontiers in Zoology, 2013, 10: 34- 34.

    [38]Riaz T, Shehzad W, Viari A, Pompanon F, Taberlet P, Coissac é. ecoPrimers: inference of new DNA barcode markers from whole genome sequence analysis. Nucleic Acids Research, 2011, 39(21): e145.

    [39]葉丹丹, 樊萌萌, 關(guān)瓊, 陳紅菊, 馬占山. 宏基因組研究的生物信息學(xué)平臺現(xiàn)狀. 動物學(xué)研究, 2012, 33(6): 574- 585.

    [40]Liu L, Li Y H, Li S L, Hu N, He Y M, Pong R, Lin D N, Lu L H, Law M. Comparison of next-generation sequencing systems. Journal of Biomedicine and Biotechnology, 2012, 2012: Article ID 251364.

    [41]賀紀(jì)正, 袁超磊, 沈菊培, 張麗梅. 土壤宏基因組學(xué)研究方法與進(jìn)展. 土壤學(xué)報(bào), 2012, 49(1): 155- 164.

    [42]解增言, 林俊華, 譚軍, 舒坤賢. DNA測序技術(shù)的發(fā)展歷史與最新進(jìn)展. 生物技術(shù)通報(bào), 2010, (8): 64- 70.

    [43]張得芳, 馬秋月, 尹佟明, 夏濤. 第三代測序技術(shù)及其應(yīng)用. 中國生物工程雜志, 2013, 33(5): 125- 131.

    [44]Roberts R J, Carneiro M O, Schatz M C. The advantages of SMRT sequencing. Genome Biology, 2013, 14(7): 405- 405.

    [45]Bienert F, De Danieli S, Miquel C, Coissac E, Poillot C, Brun J J, Taberlet P. Tracking earthworm communities from soil DNA. Molecular Ecology, 2012, 21(8): 2017- 2030.

    [46]楊晨雪, 季吟秋, 王曉陽, 楊春燕, Yu D W. 基于18S rDNA的metabarcoading技術(shù)分析土壤小型動物多樣性3種方法的比較. 中國科學(xué): 生命科學(xué), 2012, 42(12): 993- 1001.

    [47]Yu D W, Ji Y Q, Emerson B C, Wang X Y, Ye C X, Yang C Y, Ding Z L. Biodiversity soup: metabarcoding of arthropods for rapid biodiversity assessment and biomonitoring. Methods in Ecology and Evolution, 2012, 3(4): 613- 623.

    [48]Bohmann K, Monadjem A, Lehmkuhl N C, Rasmussen M, Zeale M R, Clare E, Jones G, Willerslev E, Gilbert M T. Molecular diet analysis of two african free-tailed bats (molossidae) using high throughput sequencing. PLoS One, 2011, 6(6): e21441.

    [49]Shehzad W, Riaz T, Nawaz M A, Miquel C, Poillot C, Shah S A, Pompanon F, Coissac E, Taberlet P. Carnivore diet analysis based on next-generation sequencing: application to the leopard cat (Prionailurusbengalensis) in Pakistan. Molecular Ecology, 2012, 21(8): 1951- 1965.

    [50]Soininen E M, Ravolainen V T, Br?then K A, Yoccoz N G, Gielly L, Ims R A. Arctic Small Rodents Have Diverse Diets and Flexible Food Selection. PLoS One, 2013, 8(6): e68128.

    [51]Quéméré E, Hibert F, Miquel C, Lhuillier E, Rasolondraibe E, Champeau J, Rabarivola C, Nusbaumer L, Chatelain C, Gautier L, Ranirison P, Crouau-Roy B, Taberlet P, Chikhi L. A DNA metabarcoding study of a primate dietary diversity and plasticity across its entire fragmented range. PLoS One, 2013, 8(3): e58971.

    [52]Hope P R, Bohmann K, Gilbert M T P, Zepeda-Mendoza M L, Razgour O, Jones G. Second generation sequencing and morphological faecal analysis reveal unexpected foraging behaviour byMyotisnattereri(Chiroptera, Vespertilionidae) in winter. Frontiers in Zoology, 2014, 11: 39- 39.

    [53]Takahara T, Minamoto T, Yamanaka H, Doi H, Kawabata Z. Estimation of fish biomass using environmental DNA. PLoS One, 2012, 7(4): e35868.

    [54]Pilliod D S, Goldberg C S, Arkle R S, Waits L P. Factors influencing detection of eDNA from a stream-dwelling amphibian. Molecular Ecology Resources, 2014, 14(1): 109- 116.

    [55]Jane S F, Wilcox T M, McKelvey K S, Young M K, Schwartz M K, Lowe W H, Letcher B H, Whiteley A R. Distance, flow and PCR inhibition: eDNA dynamics in two headwater streams. Molecular Ecology Resources, 2015, 15(1): 216- 227.

    [56]Maruyama A, Nakamura K, Yamanaka H, Kondoh M, Minamoto T. The release rate of environmental DNA from juvenile and adult fish. PLoS One, 2014, 9(12): e114639.

    Application of environmental DNA metabarcoding in ecology

    CHEN Lian1,2, WU Lin2, LIU Yan1, XU Haigen1,*

    1NanjingInstituteofEnvironmentalSciences,MinistryofEnvironmentalProtection,Nanjing210042,China2CollegeofLifeSciences,ChemistryandChemicalEngineering,JiangsuSecondNormalUniversity,Nanjing210013,China

    Environmental DNA metabarcoding is defined as high throughput DNA sequencing used for simultaneous molecular identification of multiple taxa in a complex environmental sample (i.e., soil, water, feces, etc.). In recent years, this method has aroused widespread concern among scientists and has been increasingly used in the field of ecology, such as biodiversity research, monitoring of aquatic, rare, or invasive species. This article describes the significance and the methods of environmental DNA metabarcoding, and highlights its application in the fields of species monitoring,detection of rare and invasive species, biodiversity, and diet analysis. Finally, this article summarizes the challenges in application of environmental DNA metabarcoding in ecological research and discusses the future prospects of its development.

    environmental DNA; metabarcoding; biodiversity; species detection

    中央級公益性科研院所基本科研業(yè)務(wù)費(fèi)專項(xiàng)(2013);江蘇省環(huán)境科技創(chuàng)新團(tuán)隊(duì)項(xiàng)目資助

    2015- 01- 15; 網(wǎng)絡(luò)出版日期:2015- 11- 17

    Corresponding author.E-mail: xhg@nies.org

    10.5846/stxb201501150125

    陳煉,吳琳,劉燕,徐海根.環(huán)境DNA metabarcoding及其在生態(tài)學(xué)研究中的應(yīng)用.生態(tài)學(xué)報(bào),2016,36(15):4573- 4582.

    Chen L, Wu L, Liu Y, Xu H G.Application of environmental DNA metabarcoding in ecology.Acta Ecologica Sinica,2016,36(15):4573- 4582.

    猜你喜歡
    條形碼測序物種
    吃光入侵物種真的是解決之道嗎?
    英語世界(2023年10期)2023-11-17 09:18:18
    杰 Sir 帶你認(rèn)識宏基因二代測序(mNGS)
    新民周刊(2022年27期)2022-08-01 07:04:49
    創(chuàng)意條形碼
    二代測序協(xié)助診斷AIDS合并馬爾尼菲籃狀菌腦膜炎1例
    傳染病信息(2021年6期)2021-02-12 01:52:58
    從條形碼到二維碼
    回首2018,這些新物種值得關(guān)注
    從條形碼到二維碼
    電咖再造新物種
    汽車觀察(2018年10期)2018-11-06 07:05:26
    條形碼大變身
    瘋狂的外來入侵物種
    又黄又爽又刺激的免费视频.| 97人妻精品一区二区三区麻豆| 亚洲一级一片aⅴ在线观看| 99九九线精品视频在线观看视频| 99热全是精品| 中国美白少妇内射xxxbb| 精品人妻偷拍中文字幕| 九九热线精品视视频播放| 午夜免费激情av| 蜜桃亚洲精品一区二区三区| 老熟妇乱子伦视频在线观看| 91av网一区二区| 久久久精品欧美日韩精品| 久久久色成人| 一级毛片aaaaaa免费看小| 成人特级黄色片久久久久久久| 草草在线视频免费看| av在线天堂中文字幕| 日韩 亚洲 欧美在线| 精品久久久噜噜| 在线a可以看的网站| 看非洲黑人一级黄片| 噜噜噜噜噜久久久久久91| 可以在线观看的亚洲视频| 亚洲在线观看片| 九九在线视频观看精品| 国产男人的电影天堂91| 人妻少妇偷人精品九色| 少妇的逼水好多| 12—13女人毛片做爰片一| 亚洲自偷自拍三级| 欧美xxxx性猛交bbbb| 亚洲丝袜综合中文字幕| 日产精品乱码卡一卡2卡三| 久久精品综合一区二区三区| 麻豆精品久久久久久蜜桃| a级一级毛片免费在线观看| 99久久久亚洲精品蜜臀av| 波多野结衣高清作品| 99久久精品国产国产毛片| 亚洲天堂国产精品一区在线| 麻豆久久精品国产亚洲av| 国产一区亚洲一区在线观看| 可以在线观看的亚洲视频| 成人欧美大片| 人妻少妇偷人精品九色| av专区在线播放| 99久久久亚洲精品蜜臀av| 少妇的逼好多水| 99在线视频只有这里精品首页| 99热这里只有是精品50| 国产女主播在线喷水免费视频网站 | 性色avwww在线观看| 蜜桃久久精品国产亚洲av| 国产精品乱码一区二三区的特点| 久久亚洲精品不卡| 久久韩国三级中文字幕| 国产亚洲av片在线观看秒播厂 | 国产av一区在线观看免费| 夫妻性生交免费视频一级片| 国产精品人妻久久久影院| 国产亚洲av嫩草精品影院| 国产色爽女视频免费观看| 国产亚洲精品久久久久久毛片| 菩萨蛮人人尽说江南好唐韦庄 | 国产淫片久久久久久久久| 午夜精品一区二区三区免费看| 日本-黄色视频高清免费观看| 久久中文看片网| 插阴视频在线观看视频| 国产色爽女视频免费观看| 亚洲一区二区三区色噜噜| 爱豆传媒免费全集在线观看| 能在线免费观看的黄片| 国产激情偷乱视频一区二区| 久久人人精品亚洲av| 免费大片18禁| av在线蜜桃| 99在线人妻在线中文字幕| 久久久国产成人免费| 免费观看精品视频网站| 久久人人精品亚洲av| 精品免费久久久久久久清纯| 亚洲不卡免费看| 在线观看66精品国产| 丰满人妻一区二区三区视频av| 18禁在线无遮挡免费观看视频| 色尼玛亚洲综合影院| 国产精品一区二区在线观看99 | 夫妻性生交免费视频一级片| 床上黄色一级片| 日韩 亚洲 欧美在线| 亚洲精华国产精华液的使用体验 | 成人av在线播放网站| 十八禁国产超污无遮挡网站| 日韩精品青青久久久久久| 久久久久久久午夜电影| 桃色一区二区三区在线观看| 老熟妇乱子伦视频在线观看| 长腿黑丝高跟| 国产亚洲欧美98| 日韩国内少妇激情av| 黄片wwwwww| 久久午夜福利片| 亚洲天堂国产精品一区在线| 此物有八面人人有两片| 亚洲欧美日韩卡通动漫| 日韩视频在线欧美| 国产精品福利在线免费观看| 卡戴珊不雅视频在线播放| 99久久人妻综合| 欧美区成人在线视频| 国产高清视频在线观看网站| 99热这里只有是精品在线观看| 伦精品一区二区三区| 黄片wwwwww| 亚洲精品久久国产高清桃花| 久久这里有精品视频免费| 欧美又色又爽又黄视频| 免费人成视频x8x8入口观看| 欧美一级a爱片免费观看看| 日本免费一区二区三区高清不卡| 国产精品三级大全| 干丝袜人妻中文字幕| 九草在线视频观看| 亚洲国产欧洲综合997久久,| 久久久国产成人精品二区| 一个人看视频在线观看www免费| 精品久久久久久久久亚洲| 欧美成人一区二区免费高清观看| 午夜视频国产福利| 国产不卡一卡二| 丝袜美腿在线中文| 一级毛片aaaaaa免费看小| 最好的美女福利视频网| 一本久久精品| 日韩欧美在线乱码| 国产高清激情床上av| 国内精品宾馆在线| 精品少妇黑人巨大在线播放 | 插阴视频在线观看视频| 久久久久免费精品人妻一区二区| 99热只有精品国产| 免费在线观看成人毛片| 美女cb高潮喷水在线观看| 久久中文看片网| 一进一出抽搐动态| 看十八女毛片水多多多| 一级毛片电影观看 | 久久久久久久久久久丰满| 国内精品美女久久久久久| 色综合站精品国产| 九九爱精品视频在线观看| av.在线天堂| 久久久久久国产a免费观看| 亚洲精品456在线播放app| 久久久久久久久久久免费av| 亚洲无线在线观看| 亚洲第一电影网av| 国产亚洲5aaaaa淫片| 成人午夜高清在线视频| 国产高清三级在线| 亚洲精品日韩在线中文字幕 | 精品一区二区三区人妻视频| 一个人免费在线观看电影| 亚洲第一电影网av| 久久热精品热| 嫩草影院入口| 欧美高清成人免费视频www| 亚洲精品成人久久久久久| 一本久久精品| 日本熟妇午夜| 99在线人妻在线中文字幕| 婷婷色av中文字幕| 两个人的视频大全免费| 一夜夜www| 国产精品麻豆人妻色哟哟久久 | 国产精品一区二区性色av| 偷拍熟女少妇极品色| 尤物成人国产欧美一区二区三区| 成人毛片a级毛片在线播放| 在线观看美女被高潮喷水网站| 别揉我奶头 嗯啊视频| 我的老师免费观看完整版| 超碰av人人做人人爽久久| 国内久久婷婷六月综合欲色啪| 悠悠久久av| 国产黄色小视频在线观看| 国产又黄又爽又无遮挡在线| 国产一区二区三区av在线 | 蜜臀久久99精品久久宅男| 亚洲精品影视一区二区三区av| 啦啦啦观看免费观看视频高清| 国产真实伦视频高清在线观看| 能在线免费看毛片的网站| 国产伦精品一区二区三区四那| 26uuu在线亚洲综合色| 久久久色成人| av天堂中文字幕网| 综合色丁香网| 国产在线男女| 久久精品夜夜夜夜夜久久蜜豆| 99在线视频只有这里精品首页| 两性午夜刺激爽爽歪歪视频在线观看| 国产国拍精品亚洲av在线观看| 日韩一区二区三区影片| 国产精品1区2区在线观看.| 成年女人永久免费观看视频| 中文在线观看免费www的网站| 蜜臀久久99精品久久宅男| 能在线免费看毛片的网站| 亚洲精华国产精华液的使用体验 | 少妇的逼水好多| 热99在线观看视频| 一进一出抽搐动态| 亚洲中文字幕一区二区三区有码在线看| 欧美性猛交╳xxx乱大交人| 亚洲国产精品久久男人天堂| 人人妻人人澡欧美一区二区| av女优亚洲男人天堂| 久久鲁丝午夜福利片| 亚洲国产欧美在线一区| 亚洲第一电影网av| 国产亚洲5aaaaa淫片| 日本一二三区视频观看| 欧美激情久久久久久爽电影| 又粗又硬又长又爽又黄的视频 | 亚洲精品亚洲一区二区| videossex国产| 97超碰精品成人国产| 亚洲综合色惰| 精品99又大又爽又粗少妇毛片| 麻豆久久精品国产亚洲av| 99国产精品一区二区蜜桃av| 亚洲国产日韩欧美精品在线观看| 色哟哟·www| 久久精品国产亚洲av天美| 亚洲精品亚洲一区二区| 麻豆成人午夜福利视频| 国产成人午夜福利电影在线观看| av又黄又爽大尺度在线免费看 | 男女那种视频在线观看| 国内揄拍国产精品人妻在线| 欧美3d第一页| 午夜福利在线在线| 精品人妻熟女av久视频| 身体一侧抽搐| 午夜老司机福利剧场| 中国美白少妇内射xxxbb| 亚洲18禁久久av| 特大巨黑吊av在线直播| av天堂在线播放| 亚洲四区av| 色综合色国产| 日本av手机在线免费观看| 国产av麻豆久久久久久久| 国产成人91sexporn| 桃色一区二区三区在线观看| 一进一出抽搐gif免费好疼| 国产精品久久久久久av不卡| 国产免费男女视频| 中文字幕熟女人妻在线| 日韩大尺度精品在线看网址| 日韩欧美三级三区| 中文字幕人妻熟人妻熟丝袜美| 美女大奶头视频| 成人性生交大片免费视频hd| 色哟哟哟哟哟哟| 亚洲aⅴ乱码一区二区在线播放| 成人综合一区亚洲| 99热网站在线观看| 亚洲人成网站在线播| 此物有八面人人有两片| 国产精品美女特级片免费视频播放器| 日本-黄色视频高清免费观看| 欧美激情在线99| 91精品一卡2卡3卡4卡| 国产精品一及| 欧美成人精品欧美一级黄| 亚洲av中文av极速乱| 午夜精品在线福利| 全区人妻精品视频| 日韩高清综合在线| 在线a可以看的网站| 亚洲七黄色美女视频| 久久久欧美国产精品| 亚洲18禁久久av| 别揉我奶头 嗯啊视频| 精品人妻视频免费看| 久久精品夜色国产| 国产乱人视频| 精品人妻熟女av久视频| 国产高潮美女av| 97热精品久久久久久| 国产精品乱码一区二三区的特点| 午夜亚洲福利在线播放| 久久久久久久久久成人| av在线天堂中文字幕| 美女被艹到高潮喷水动态| av免费在线看不卡| 日本三级黄在线观看| 97人妻精品一区二区三区麻豆| 51国产日韩欧美| 国产一级毛片七仙女欲春2| 亚洲精品影视一区二区三区av| 一个人看视频在线观看www免费| 成人特级av手机在线观看| 日韩大尺度精品在线看网址| 国产男人的电影天堂91| 黄色欧美视频在线观看| 国产精品一区二区三区四区久久| 国产日韩欧美在线精品| 波多野结衣高清作品| 精品人妻一区二区三区麻豆| 亚洲av男天堂| 国产av一区在线观看免费| 26uuu在线亚洲综合色| 成人二区视频| av天堂在线播放| 亚洲成人久久爱视频| 免费观看在线日韩| 国产亚洲av嫩草精品影院| 在线观看66精品国产| 免费人成视频x8x8入口观看| 免费无遮挡裸体视频| 最近中文字幕高清免费大全6| 日日摸夜夜添夜夜添av毛片| av在线天堂中文字幕| 少妇的逼水好多| 有码 亚洲区| 亚洲av成人精品一区久久| 久久精品国产清高在天天线| 99久久成人亚洲精品观看| eeuss影院久久| 国产精品伦人一区二区| www日本黄色视频网| 夜夜夜夜夜久久久久| 久久精品人妻少妇| 午夜激情欧美在线| 日日啪夜夜撸| 国产视频内射| 久久人人精品亚洲av| 麻豆精品久久久久久蜜桃| 亚洲精品自拍成人| 麻豆乱淫一区二区| av在线观看视频网站免费| 久久精品夜色国产| 国产av在哪里看| 精品无人区乱码1区二区| 欧美精品一区二区大全| 久久精品夜色国产| 免费看a级黄色片| 成人av在线播放网站| 国产免费一级a男人的天堂| 国产亚洲av嫩草精品影院| 全区人妻精品视频| 老司机影院成人| 你懂的网址亚洲精品在线观看 | 男女那种视频在线观看| 国产亚洲91精品色在线| 99精品在免费线老司机午夜| 亚洲内射少妇av| 99精品在免费线老司机午夜| 成人特级av手机在线观看| 午夜视频国产福利| 青春草亚洲视频在线观看| 人人妻人人看人人澡| 久久精品综合一区二区三区| 国产精品不卡视频一区二区| 免费看美女性在线毛片视频| 欧美激情国产日韩精品一区| 亚洲aⅴ乱码一区二区在线播放| 岛国毛片在线播放| 中国美女看黄片| 变态另类成人亚洲欧美熟女| 久久国内精品自在自线图片| 一进一出抽搐gif免费好疼| 成人av在线播放网站| 国产精品综合久久久久久久免费| 国产精品美女特级片免费视频播放器| 99国产极品粉嫩在线观看| 少妇人妻精品综合一区二区 | 高清毛片免费观看视频网站| 欧美三级亚洲精品| 简卡轻食公司| 久久99精品国语久久久| 特大巨黑吊av在线直播| 久久久成人免费电影| 亚洲人成网站在线播| 久久这里只有精品中国| 日韩欧美国产在线观看| 又粗又硬又长又爽又黄的视频 | 久久中文看片网| 久久久a久久爽久久v久久| 日韩三级伦理在线观看| 国产激情偷乱视频一区二区| 18禁在线播放成人免费| 国产精品99久久久久久久久| 99久久无色码亚洲精品果冻| 亚洲中文字幕日韩| 欧美色视频一区免费| 悠悠久久av| 亚洲人与动物交配视频| 国产精品无大码| 一级av片app| 亚洲精品日韩av片在线观看| 久久精品国产99精品国产亚洲性色| 18禁黄网站禁片免费观看直播| 中文字幕免费在线视频6| 久久精品国产亚洲av香蕉五月| 国内久久婷婷六月综合欲色啪| av免费在线看不卡| 国产精品电影一区二区三区| av在线天堂中文字幕| 两性午夜刺激爽爽歪歪视频在线观看| 亚洲精品影视一区二区三区av| 九九在线视频观看精品| 亚洲欧美中文字幕日韩二区| 免费看a级黄色片| 99久国产av精品| 亚洲欧洲国产日韩| 国产成人a区在线观看| 亚洲av成人av| 中文资源天堂在线| 卡戴珊不雅视频在线播放| 欧美日韩综合久久久久久| 男插女下体视频免费在线播放| 成人毛片60女人毛片免费| 51国产日韩欧美| 亚洲乱码一区二区免费版| 亚洲av不卡在线观看| 国产成人精品婷婷| 在线国产一区二区在线| 少妇裸体淫交视频免费看高清| 亚洲欧美日韩东京热| 美女内射精品一级片tv| 天堂中文最新版在线下载 | 寂寞人妻少妇视频99o| 国产精品美女特级片免费视频播放器| 日韩欧美精品v在线| 亚洲久久久久久中文字幕| 国产精品一及| 免费一级毛片在线播放高清视频| 免费看美女性在线毛片视频| 国产精品一区二区在线观看99 | 一边亲一边摸免费视频| 在线播放无遮挡| 国产高清有码在线观看视频| av天堂中文字幕网| 亚洲av免费高清在线观看| 男女啪啪激烈高潮av片| 色视频www国产| 黄片wwwwww| 乱码一卡2卡4卡精品| 精品久久久久久久久av| 搞女人的毛片| 99riav亚洲国产免费| 日韩欧美 国产精品| 日本撒尿小便嘘嘘汇集6| 亚洲精品日韩av片在线观看| 乱码一卡2卡4卡精品| 国产黄a三级三级三级人| 亚洲精华国产精华液的使用体验 | 亚洲中文字幕一区二区三区有码在线看| 国产成人aa在线观看| 2021天堂中文幕一二区在线观| 一卡2卡三卡四卡精品乱码亚洲| 免费av毛片视频| 可以在线观看毛片的网站| 99九九线精品视频在线观看视频| .国产精品久久| 亚洲av免费高清在线观看| 2022亚洲国产成人精品| 高清毛片免费看| 欧美精品国产亚洲| 日日摸夜夜添夜夜添av毛片| 亚洲高清免费不卡视频| 日本免费一区二区三区高清不卡| 中出人妻视频一区二区| 国产精品国产三级国产av玫瑰| 麻豆一二三区av精品| 精品免费久久久久久久清纯| 夜夜看夜夜爽夜夜摸| 热99re8久久精品国产| 3wmmmm亚洲av在线观看| 久久久精品大字幕| 青青草视频在线视频观看| 久久久久久久久久久免费av| 成人av在线播放网站| 亚洲国产精品sss在线观看| 国产在视频线在精品| 色哟哟·www| 国产精品乱码一区二三区的特点| 久久国内精品自在自线图片| 久久99蜜桃精品久久| 国内久久婷婷六月综合欲色啪| 国产黄a三级三级三级人| 一个人看视频在线观看www免费| 18禁裸乳无遮挡免费网站照片| 91久久精品电影网| 成熟少妇高潮喷水视频| 国产亚洲91精品色在线| 不卡视频在线观看欧美| 一本一本综合久久| 欧美日韩在线观看h| 国产亚洲av片在线观看秒播厂 | 日本与韩国留学比较| 直男gayav资源| 国产一区二区亚洲精品在线观看| 久久鲁丝午夜福利片| 国内精品宾馆在线| 国产真实伦视频高清在线观看| 中文在线观看免费www的网站| a级毛片免费高清观看在线播放| 天天一区二区日本电影三级| 久久午夜亚洲精品久久| 精品人妻视频免费看| 国产高清三级在线| 综合色丁香网| 亚洲中文字幕一区二区三区有码在线看| 成年免费大片在线观看| 久久久久久伊人网av| 日韩欧美在线乱码| 能在线免费看毛片的网站| 久久欧美精品欧美久久欧美| 国内久久婷婷六月综合欲色啪| 久久久久久久午夜电影| 国产精品一区二区在线观看99 | 久久精品久久久久久噜噜老黄 | 国产成人精品婷婷| 欧美在线一区亚洲| 美女脱内裤让男人舔精品视频 | 亚洲欧美精品综合久久99| www.色视频.com| 亚洲欧美精品综合久久99| 免费av毛片视频| 亚洲丝袜综合中文字幕| 最近视频中文字幕2019在线8| 国产免费男女视频| 性插视频无遮挡在线免费观看| 国产亚洲精品久久久久久毛片| 日韩亚洲欧美综合| 久久精品91蜜桃| 最近2019中文字幕mv第一页| 亚洲激情五月婷婷啪啪| 久久99热6这里只有精品| 免费观看a级毛片全部| 99久久成人亚洲精品观看| 长腿黑丝高跟| 日韩成人伦理影院| 免费人成视频x8x8入口观看| 国产精品久久久久久久电影| 国产亚洲av嫩草精品影院| 91久久精品电影网| 成年免费大片在线观看| 草草在线视频免费看| 大又大粗又爽又黄少妇毛片口| 淫秽高清视频在线观看| 亚洲在久久综合| 久久午夜福利片| av在线蜜桃| 可以在线观看的亚洲视频| 久久久午夜欧美精品| 久久精品久久久久久久性| av免费观看日本| 国产av一区在线观看免费| 亚洲欧洲日产国产| 成人永久免费在线观看视频| av在线播放精品| 在线观看66精品国产| 国产亚洲精品久久久com| 悠悠久久av| 久久人人爽人人片av| 国产真实伦视频高清在线观看| 精品久久久久久成人av| 91麻豆精品激情在线观看国产| 国产综合懂色| 国产黄a三级三级三级人| 久久久久久久久中文| 男女边吃奶边做爰视频| 久久久久久久久中文| 免费大片18禁| 白带黄色成豆腐渣| 春色校园在线视频观看| 欧美一区二区精品小视频在线| 少妇熟女欧美另类| 哪个播放器可以免费观看大片| 国内久久婷婷六月综合欲色啪| 国产精品一区二区三区四区久久| 少妇的逼水好多| 两个人视频免费观看高清| 精品不卡国产一区二区三区| 国产片特级美女逼逼视频| 久久久午夜欧美精品| 欧美成人免费av一区二区三区| 久久久成人免费电影| 99精品在免费线老司机午夜| 国产精品女同一区二区软件| 国语自产精品视频在线第100页| 日本欧美国产在线视频| 91午夜精品亚洲一区二区三区| 久久精品久久久久久久性| 九九久久精品国产亚洲av麻豆| 亚洲三级黄色毛片| 日韩强制内射视频| 国产精品.久久久| 久久久国产成人免费| 国产午夜精品久久久久久一区二区三区| 男人狂女人下面高潮的视频| 高清日韩中文字幕在线| a级一级毛片免费在线观看| 最近中文字幕高清免费大全6| 日本-黄色视频高清免费观看| 国产视频内射| www.色视频.com|