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

    猴頭菇DNA條形碼的構(gòu)建及驗證

    2017-10-16 05:44:07皇甫永冠閆寶松張躍新宗憲春
    中國林副特產(chǎn) 2017年5期

    皇甫永冠,閆寶松,張躍新,宗憲春

    (1.東北林業(yè)大學鹽堿地研究中心,哈爾濱 150040;2.黑龍江省林副特產(chǎn)研究所,黑龍江省非木質(zhì)林產(chǎn)品研發(fā)重點實驗室,黑龍江 牡丹江 157011;3 牡丹江師范學院,黑龍江 牡丹江 157000)

    猴頭菇DNA條形碼的構(gòu)建及驗證

    皇甫永冠1,閆寶松2,張躍新2,宗憲春3

    (1.東北林業(yè)大學鹽堿地研究中心,哈爾濱 150040;2.黑龍江省林副特產(chǎn)研究所,黑龍江省非木質(zhì)林產(chǎn)品研發(fā)重點實驗室,黑龍江 牡丹江 157011;3 牡丹江師范學院,黑龍江 牡丹江 157000)

    為構(gòu)建猴頭菇品種專屬的DNA條形碼。應用分子生物學方法、選定β-tubulin基因作為條形碼序列的目標基因,依據(jù)NCBI中猴頭菇β-tubulin基因序列的全長,設計了一對用于擴增目標基因的引物:Betu-F1:ATGCGTGAAATCGTCCACC;Betu-R2:GAAGACGGGGGAAAGGAAC。應用上述引物對現(xiàn)有的19個猴頭菇的品種進行PCR鑒定。測序后的序列結(jié)合BLAST、Clustal、MEGA、BioEdit等軟件分析、比對得到多態(tài)性位點,分析多態(tài)性位點設計出構(gòu)建條形碼序列的特異性引物。經(jīng)驗證表明,猴頭菇RT9品種的3個特有的條形碼成功獲得。試驗結(jié)果表明,實驗成功獲得的DNA條形碼序列可以將猴頭菇RT9品種與猴頭菇的其它品種區(qū)分鑒定。

    猴頭菇;DNA條形碼;內(nèi)含子多態(tài)性

    Abstract:The aim was to construct the DNA barcode of Hericium erinaceus. Using molecular biology method,Selected β-tubulin gene as the target gene of bar code sequence,A pair of primers for amplification of the target gene was designed according to the full-length of the β-tubulin gene sequence of NCBI:Betu-F1:ATGCGTGAAATCGTCCACC;Betu-R2:GAAGACGGGGGAAAGGAAC.The above primers were used to identify the 19 species of Hericium erinaceus of PCR.Sequencing of the sequence with BLAST, Clustal, MEGA, BioEdit and other software analysis, comparison of polymorphic loci,Analysis of polymorphic loci was designed to construct a specific primer sequence. The results showed that the three unique bar codes of Hericium erinaceus RT9 were successfully obtained.The experiment result shows:the DNA barcode sequence obtained from the experiment can distinguish the RT9 and other varieties of Hericium erinaceus

    Keywords:Hericiumerinaceus; DNA barcoding; Intron polymorphis

    0 緒論

    猴頭菇因營養(yǎng)豐富、味道鮮美及具抗癌等功效,近些年來成為相關行業(yè)的研究熱點,但對它的研究多集中在深加工等方面。中國現(xiàn)有的猴頭菇有幾十種到數(shù)百種,它們的營養(yǎng)成分及功效各不相同,而且名稱混亂,給引種、栽培和開發(fā)利用工作帶來較大的麻煩[1]。猴頭菇菌種的區(qū)分鑒定對于今后的育種工作及種質(zhì)資源的保護意義重大,但由于現(xiàn)在多數(shù)只是從形態(tài)學的角度來鑒定菌種,其準確性不是十分的可靠,研究表明,將分子標記應用在菌種鑒定上將會大幅度的增加鑒定效率。根據(jù)前人對物種的鑒定研究來看,DNA條形碼技術(shù)在區(qū)分物種時具簡單、快速、準確性高等優(yōu)點,因此本實驗結(jié)合內(nèi)含子多態(tài)性,力爭構(gòu)建出猴頭菇的專屬DNA條形碼,為今后食用菌行業(yè)的種質(zhì)資源鑒定提供行之有效的手段。自Hebert等[2-3]提議使用線粒體基因細胞色素C氧化酶(COI)基因片段作為動物的DNA條形碼后,DNA條形碼技術(shù)得以迅速發(fā)展,在動物、植物、微生物的鑒定及分類上開始廣泛的應用[4]。DNA barcoding技術(shù)在動物中的研究進展:COI基因通常作為動物的條形碼[5-6]目前主要運用于鳥類[7-16]、魚類[17-21]、昆蟲[22-24]的分類鑒定。DNA barcoding在植物中研究進展:2009年國際生命條形碼聯(lián)盟植物工作組(Consortiumfor the Barcode of Life Plant Working Group)初步確定,并推薦使用葉綠體基因片段的rbcL和matK作為植物的DNA條形碼基因[25]。DNA barcoding技術(shù)在微生物中的研究進展:目前主要在釀酒酵母(Saccharomyces cerevisiae)突變菌株的識別[26],毛霉屬(Mucor)[27]、青霉屬(Penicillium)[28]、口蘑屬(Tricholoma)[29]等物種鑒別中有應用。相比于動物、植物,菌物DNA barcoding技術(shù)的研究起步相對較晚,研究的也較少,到目前為止菌物的DNA條形碼的標準序列還沒有選定。依據(jù)Zhao等[30]、Samson et al.[31],分別對叢赤殼類真菌、青霉亞屬的代表性菌株進行的序列分析表明,β-tubulin基因在真菌鑒定中是一種理想的DNA barcoding標記基因,無論是在PCR鑒定還是測序上都表現(xiàn)出很高的成功率,而且不存在種間和種內(nèi)的交疊,鑒定效果優(yōu)于ITS。因此,本實驗選定β-tubulin 基因作為構(gòu)建猴頭菇品種的條形碼的目標基因,構(gòu)建猴頭菇品種的專屬條線碼。本研究針對19個猴頭菇的品種構(gòu)建猴頭菇品種的DNA條形碼,經(jīng)驗證成功構(gòu)建的條形碼序列在所有猴頭菇品種中均適用,因此,DNA條形碼技術(shù)在菌物鑒定中應用切實可行,為今后菌物菌種鑒定提供了技術(shù)手段。

    1 材料與方法

    1.1 實驗材料

    供試菌株見表1。

    表1 供試菌種

    續(xù)表1

    序號名稱采集地9野猴3黑龍江省牡丹江市三道關林區(qū)野生10RT9黑龍江省農(nóng)業(yè)經(jīng)濟學院11RT10黑龍江省農(nóng)業(yè)經(jīng)濟學院12RT11黑龍江省農(nóng)業(yè)經(jīng)濟學院13RT12黑龍江省農(nóng)業(yè)經(jīng)濟學院14RT13黑龍江省農(nóng)業(yè)經(jīng)濟學院15RT14黑龍江省農(nóng)業(yè)經(jīng)濟學院16RT16黑龍江省農(nóng)業(yè)經(jīng)濟學院17RT17黑龍江省農(nóng)業(yè)經(jīng)濟學院18RT18黑龍江省農(nóng)業(yè)經(jīng)濟學院19RT19黑龍江省農(nóng)業(yè)經(jīng)濟學院

    1.2 實驗方法

    1.2.1 DNA條形碼的構(gòu)建

    (1)選取材料:19個猴頭菇品種(見表1)。

    (2)獲得及培養(yǎng)19個猴頭菇品種的菌絲體。

    (3)DNA的提?。罕緦嶒灢捎肅TAB法對供試菌株的DNA進行提取。

    (4)設計引物及擴增基因

    設計引物

    依據(jù)NCBI中猴頭菇β-tubulin基因的全長設計了一對引物如下:Betu-F:ATGCGTGAAATCGTCCACC;Betu-R:GAAGACGGGGGAAAGGAAC。

    擴增目的基因

    PCR的體系:10×buffer 緩沖液 2 μL,d NTP(2.5mM each) 1.6 μL,上游引物(10μM) 1 μL,下游引物(10μM) 1 μL,r-Taq DNA 聚合酶(5U/ μL) 0.3 μL,基因組 DNA 1 μL,用滅過菌的去離子水補足至20μL。

    PCR反應程序:94℃預變性5min;94℃變性30s,55℃退火30s , 72℃延伸1min20s,35個循環(huán),72℃再延伸 20 min。

    (5)T-A克隆

    目的基因的回收及檢測

    利用瓊脂糖凝膠DNA回收試劑盒對目標條帶進行切膠回收,1%的凝膠濃度對回收的目的基因進行檢測。

    連接反應

    體系10μL:2×Rapid ligation Buffer 5μL,pGEM-T Easy vector 1μL,PCR product,1μL,T4 DNA Ligase,1μL,ddH2O,2μL,16℃過夜連接不超過16h。

    轉(zhuǎn)化

    熱激轉(zhuǎn)化:16℃連接過夜后,將10μL連接產(chǎn)物加入到感受態(tài)細胞中(DH5a)中,冰浴30min,42℃熱激90s,在冰浴2min,向其中加入500μL B液體培養(yǎng)基,37℃180轉(zhuǎn)震蕩培養(yǎng)40~60min,取震蕩培養(yǎng)后的菌液200μL涂布在加有AMP+IPTG+X-Gal的固體LB培養(yǎng)基的平板上(藍白斑篩選),為防止因立即倒置培養(yǎng),造成菌液流動現(xiàn)象,將涂布后的菌版在37℃正向培養(yǎng)1h,之后倒置過夜培養(yǎng),篩選陽性單菌落。

    篩選重組子

    挑取藍白斑篩選平板上白色的單菌落,至5mL含AMP 的LB液體培養(yǎng)基中,將震蕩培養(yǎng)箱的溫度調(diào)制37℃、轉(zhuǎn)速調(diào)至180r/min,培養(yǎng)約6h左右至菌液渾濁,進行菌液PCR檢測。

    體系20μL:10×buffer 緩沖液 2 μL,d NTP(2.5mM each) 1.6 μL,上游引物(10μM) 1 μL,下游引物(10μM) 1 μL,r-Taq DNA 聚合酶(5U/ μL) 0.3 μL,菌液 1 μL,用滅過菌的去離子水定容至20μL。

    PCR反應程序:94℃預變性5min;94℃變性 30s,60℃退火 30s,72℃延伸1min20s 35個循環(huán),72℃再延伸 5 min。

    PCR產(chǎn)物采用1%的瓊脂糖凝膠電泳對其進行檢測,陽性菌液保存?zhèn)溆谩?/p>

    重組子測序

    取出保存的陽性菌液送至生物公司(promega)進行測序。

    測序結(jié)果分析

    測序成功的序列結(jié)合primmer5.0對其進行整理,用BLAST軟件對單一序列進行比對,得到有效的序列,利用Clustal軟件對有效的序列進行多重比對,找出每個品種的代表性序列,再將代表性序列應用MEGA 5.1的鄰接法(Neighbor-Joining)構(gòu)建系統(tǒng)發(fā)育樹進、應用BioEdit軟件做出更直觀的多態(tài)性位點圖,分析整體多態(tài)性位點的分布情況,再將系統(tǒng)發(fā)育進化樹的每個大的分支的序列,應用BioEdit及Clustal軟件進行多重比對,分析每個大分支的多態(tài)性位點,進而找出成功設計DNA條形碼的多態(tài)性位點。

    (6)確定特異引物并合成

    分析多態(tài)性位點,結(jié)合primer軟件設計每個特異品種的特有引物(如果后續(xù)的條形碼驗證成功,此引物則為條形碼序列的特異引物)。將構(gòu)建的特異引物送至生物公司(上海生工)合成。

    1.2.2 DNA條形碼的驗證

    (1)現(xiàn)有菌種中的驗證:應用上述合成的特異性引物,分別對猴頭菇19個品種的基因組DNA進行PCR驗證,具體步驟如下:

    ①從樣品組織中提取基因組DNA;

    ②用步驟1所提取的基因組DNA為模板,用上述合成的特異性引物,進行聚合酶鏈式反應;

    ③對步驟2擴增得到的PCR產(chǎn)物進行電泳檢測,其特征是擴增片段長度,應與合成的特異性引物所擴增出目標條帶大小一致。

    (2)所有猴頭菇品種中的適用性驗證:把上述已構(gòu)建成功的DNA條形碼序列與GenBank數(shù)據(jù)庫中已發(fā)表的含有猴頭菇β-tubulin基因的序列(序列登錄號:>gi|386872545|>gi|63333592|>gi|84626253|>gi|84626251| >gi|84626249| >gi|84626247|>gi|84626245|>gi|84626243| >gi|84626241|>gi|84626239|>gi|84626237| >gi|84626235| >gi|84626233| >gi|84626231|>gi|84626229| >gi|84626227|>gi|84626225| >gi|84626223| >gi|84626221|>gi|84626219|>gi|84626217|>gi|84626215|>gi|84626213|>gi|84626212| >gi|84626210| >gi|84626208|>gi|84626206|>gi|84626204| >gi|84626202|>gi|84626200| >gi|84626198| >gi|84626196| >gi|84626194| >gi|84626192|>gi|84626190|>gi|84626188| >gi|84626186| >gi|84626184|>gi|84626182| )進行比對如果構(gòu)建成功的條形碼序列與已發(fā)表的序列間存在差異,說明該條形碼序列可以很好地將此品種與其他的猴頭菇品種區(qū)分開。

    2 結(jié)果與分析

    2.1 DNA提取

    提取的DNA采用1%的瓊脂糖凝膠電泳進行檢測,結(jié)果顯示所有樣本DNA均成功提取,見圖1。

    圖1 DNA提取檢測

    2.2 目的基因切膠回收

    依據(jù)NCBI中猴頭菇β-tubulin基因的全長設計的一對引物Betu-F1:ATGCGTGAAATCGTCCACC;Betu-R2:GAAGACGGGGGAAAGGAAC,用其對目的基因進行擴增,PCR產(chǎn)物采用1%的瓊脂糖凝膠電泳檢測。所有樣本均成功擴增,條帶大小約1.3kb,切膠回收電泳結(jié)果見圖2。

    圖2 目的基因切膠回收

    2.3 陽性鑒定

    在藍白斑篩選的平板上挑取白色菌落,放入LB液體培養(yǎng)基內(nèi)搖菌,做菌液PCR,PCR產(chǎn)物用1%瓊脂糖凝膠電泳檢測篩選出陽性菌。因為,所需要的陽性菌并非一次性檢出,所以實驗過程中將每次檢測出的陽性菌液均進行保存,待實驗結(jié)束后將所有檢測出的陽性菌的菌液PCR產(chǎn)物集中電泳(見圖3)。

    圖3 菌液陽性鑒定

    2.4 系統(tǒng)發(fā)育進化樹

    應用MEGA軟件對代表性序列進行分析得出系統(tǒng)發(fā)育進化樹圖(見圖4),從圖片上可直觀看出猴頭菇RT9品種是一條獨立的分支,說明RT9品種與其它猴頭菇品種間親緣關系較遠,構(gòu)建條形碼的成功率會增加。

    圖4 系統(tǒng)發(fā)育進化樹

    注:圖1中2.3、2.4、2.5為猴99的代表性序列;4.1、4.2為野猴3的代表性序列;RT16-2、RT16-3為RT16的代表性序列;HAI1、HAI3為海猴的代表性序列;RT17-1、RT17-2為RT17的代表性序列;RT18-2、RT18-3為RT18的代表性序列;RT19-1、RT19-2為RT19的代表性序列;RT10-3、RT10-4為RT10的代表性序列;RT13-1、RT13-2為RT13的代表性序列;RT12-1、RT12-2為RT12的代表性序列;RT9-1、RT9-3為RT9的代表性序列;RT11-4、RT11-5為RT11的代表性序列;RT14-3、RT14-4為RT14的代表性序列;HOU4-1、HOU4-2為猴4的代表性序列;HA2、HA5為HA的代表性序列;HOU2、HOU3為猴的代表性序列;HOU5-1、HOU5-2為猴5的代表性序列;HOU9-2、HOU9-3為猴9的代表性序列。

    2.5 多態(tài)性位點分析

    將測序后的所有序列應用BLAST、Clustal軟件篩查分析,得到每個品種具有代表性的序列,結(jié)合MEGA軟件構(gòu)建系統(tǒng)發(fā)育進化樹(見圖4)。但是由于序列過多,導致分析多態(tài)性位點、設計特異性引物難度增加,為此,將系統(tǒng)發(fā)育進化樹的每個大的分支的序列,應用Clustal軟件分別進行多重比對,每個分支的多重序列比對結(jié)果見附錄1,由于比對結(jié)果占篇幅龐大,附錄1中僅列出構(gòu)建成功的條形碼序列所在的分支的多重比對結(jié)果,根據(jù)多重比對的結(jié)果分析多態(tài)性位點,在此基礎上找出成功設計DNA條形碼的多態(tài)性位點(附錄1已標明),設計出每個品種特有的專屬條形碼引物。

    2.6 條形碼驗證

    分析多態(tài)性位點設計出了幾十對用于構(gòu)建條形碼的引物,經(jīng)逐一驗證發(fā)現(xiàn),猴頭菇RT9的3個條形碼成功構(gòu)建(見圖5-a;5-b),應用這三條條形碼序列可以將RT9與猴頭菇的其他18個品種區(qū)分開。對應的3個引物為CRT9-3F1:CGGTACGAGGATCGCGCAAAC;CRT9-3R1GCGACCGTCTT

    CCCACA;(擴增序列長度為353bp)、CRT9-3F1:CGGTACGAGGATCGCGCAAAC;CRT9-3R2:GTGTCAAGTCGCAAAAAT(擴增序列長度為483bp)、RT9-2F:TCCTCACACACACTGGTAC;RT9-2R:GTACGGATGTTGCAGTAG(擴增序列長度為723bp),應用這3對引物擴增出的條形碼序列。

    圖5-a DNA條形碼驗證

    圖5-b DNA條形碼驗證

    2.7 3個條形碼在整個猴頭菇品種中的適用性檢測

    把上述已構(gòu)建成功的3個猴頭菇RT9的DNA條形碼序列與GenBank數(shù)據(jù)庫中已發(fā)表的含有猴頭菇β-tubulin基因的序列(序列登陸號見1.2.2)進行比對,發(fā)現(xiàn)這3個條形碼序列與已發(fā)表的序列間存在差異,因此應用這3個序列可以很好地將RT9與其它的猴頭菇品種區(qū)分開,所以這3個序列可以作為鑒定猴頭菇RT9的專屬條形碼。

    3 討論

    3.1 構(gòu)建條形碼目標基因的選擇

    動物中通常作為條形碼片段的基因為COI基因[5-6],也是目前研究最為完善的。最初作為植物條形碼片段的基因為葉綠體基因片段的rbcL和matK[25],后來研究人員對高等植物、真菌、藻類的大量ITS2序列進行研究分析,結(jié)果顯示,ITS2在種級分辨上表現(xiàn)出很高的鑒定成功率。所以,建議以ITS2可以作為新的DNA條形碼,在鑒定綠色植物及部分真菌時使用[32-33]。

    相比于動物、植物,菌物DNA條形碼的研究及應用起步較晚,目前為止也沒有明確作為菌物條形碼片段的基因,前人研究發(fā)現(xiàn),β-tubulin基因在叢赤殼類真菌的研究中表現(xiàn)出高的 PCR 和測序成功率,且沒有種間及種內(nèi)的交疊,在一些菌種鑒定過程中,β-tubulin基因表現(xiàn)出了ITS2還要高的鑒定成功率[30-31],因此,本實驗選擇了β-tubulin基因作為條形碼片段的目標基因,并依據(jù)NCBI中猴頭菇β-tubulin基因的全長設計了一對引物用于擴增目的基因。

    3.2 假陽性分析

    依據(jù)多態(tài)性位點合成的幾十對用于擴增特異性條形碼的引物,在條形碼驗證時僅有RT9的3個條形碼成功獲得,可能導致此結(jié)果的原因如下:①在測序時相關的堿基錯讀或缺失,②菌液PCR過程中可能發(fā)生堿基錯配等原因。

    3.3 實驗體系優(yōu)化

    重組子鑒定過程中,開始陽性檢出率非常低,存在非特異性擴增現(xiàn)象,經(jīng)分析可能是如下原因:①可能實驗體系存在問題,如,Taq酶、T載體等的使用是否合適,②可能由于PCR的退火溫度過低,這樣雖然可以增加擴增的產(chǎn)量,但引物和模板之間的錯配現(xiàn)象也會相應增加,就會導致非特異性擴增。改進措施如下:①將實驗原來用的生工的T載體更換為大連寶生物的pGEM-T載體。②將最初使用的Taq酶(上海生工)更換為r-Taq酶(大連寶生物)。③通過提高PCR的退火溫度,由原來的55℃增至60℃。實驗體系經(jīng)上述優(yōu)化后,陽性檢出率低、非特異性擴增的現(xiàn)象均有所改進。

    4 結(jié)論

    本實驗應用分子生物學的研究手段,成功構(gòu)建出猴頭菇RT9品種的3個條形碼。經(jīng)驗證表明,這3個特異性條形碼在所有猴頭菇品種中均適用。因此,DNA條形碼技術(shù)在菌物的鑒定中應用切實可行。

    [1]Bunyard B A, Nicholson M S, Royse D J, et al. A systematic assessment of Morchella using RFLP analysis of the 28S ribosomal RNA gene[J]. Mycologia, 1994, 86(6):762-772.

    [2]Hebert P D, Cywinska A, Ball S L, et al. Biological identifications through DNA barcodes.[J]. Proceedings of The Royal Society B: Biological Sciences, 2003, 270(1512): 313-321.

    [3]Hebert P D, Ratnasingham S, De Waard J R, et al. Barcoding animal life: cytochrome c oxidase subunit 1 divergences among closely related species[J]. Proceedings of The Royal Society B: Biological Sciences, 2003,270:96-99.

    [4]裴男才,陳步峰.生物DNA條形碼:十年發(fā)展歷程、研究尺度和功能[J]. 生物多樣性,2013,2l(5):616-627.

    [5]Waugh J. DNA barcoding in animal species: progress, potential and pitfalls.[J]. BioEssays, 2007, 29(2):188-197.

    [6]Kwong S, Srivathsan A, Meier R, et al. An update on DNA barcoding: low species coverage and numerous unidentified sequences[J]. Cladistics, 2012, 28(6):639-644.

    [7]Hebert P D, Stoeckle M Y, Zemlak T S, et al. Identification of Birds through DNA Barcodes[J]. PLOS Biology, 2004, 2(10):1657-1663.

    [8]蔡延森, 張修月, 岳碧松,等. 我國8種猛禽的DNA條形碼技術(shù)研究[J]. 四川動物, 2009, 28(3):334-340.

    [9]Lijtmaer D A, Kerr K C, Barreira A S, et al. DNA Barcode Libraries Provide Insight into Continental Patterns of Avian Diversification[J]. PLOS ONE, 2011, 6(7).doi:10.1371/journal.pone.0020744.

    [10]馬明義, 閆穎, 王譯偉,等. 我國32種鳥類DNA條形碼分析[J]. 四川動物, 2012, 31(5):729-733.

    [11]Cai Y, Yue B, Jiang W, et al. DNA barcoding on subsets of three families in Aves[J]. Mitochondrial DNA, 2010,21(3-4):132-137.

    [12]戴傳銀, 張瑞瑩, 尹祚華,等. 部分山雀科鳥類的DNA條形碼與物種識別[J]. 動物分類學報, 2010, 35(4):835-841.

    [13]蔡延森, 張修月, 岳碧松,等. 我國8種猛禽的DNA條形碼技術(shù)研究[J]. 四川動物, 2009, 28(3):334-340.

    [14]屠云潔,高玉時,蘇一軍,等.我國部分地方雞種COⅠ基因多態(tài)性及其分子系統(tǒng)進化 研究[J].安徽農(nóng)業(yè)大學學報,2011,38(1):39-42.

    [15]徐向明. 我國3個地方品種鴨線粒體DNA COⅠ基因的DNA條形碼初步分析[J]. 畜牧與獸醫(yī), 2008, 40(11):51-53.

    [16]虞德兵, 陸應林, 徐昊翔,等. 基于線粒體COI基因序列分析家鴨系統(tǒng)發(fā)育關系[J]. 南京農(nóng)業(yè)大學學報, 2011, 34(6):109-114.

    [17]Ward R D, Zemlak T S, Innes B H, et al. DNA barcoding Australia's fish species.[J]. Philosophical Transactions of the Royal Society B, 2005, 360(1462):1847-1857.

    [18]Wong L L,Peatman E,L J G.DNA Barcoding of Cat?sh:Species Authentication and Phylogenetic Assessment[J].Plos One, 2011, 3(6):2578-2585.

    [19]Bhattacharjee M J, Laskar B A, Dhar B, et al. Identification and re-evaluation of freshwater catfishes through DNA barcoding.[J]. PLOS ONE, 2012, 7(11).doi:10.1371/journal.pone.0049950.

    [20]寧平, 劉靜. 基于線粒體COI基因的DNA條形碼探討印度-西太平洋沿岸6種金線魚屬魚類分類及親緣關系[C]. 中國海洋湖沼學會魚類學分會、中國動物學會魚類學分會.2012年學術(shù)研討會論文摘要匯編 2012:42-45.

    [21]彭居俐, 王緒禎, 王丁,等. 基于線粒體CO1基因序列的DNA條形碼在鯉科鲌屬魚類物種鑒定中的應用[J]. 水生生物學報, 2009, 33(2):271-276.

    [22]Hajibabaei M, Janzen D H, Burns J M, et al. DNA barcodes distinguish species of tropical Lepidoptera[J]. Proceedings of the National Academy of Sciences of the United States of America, 2006, 103(4):968-971.

    [23]Hajibabaei M, Singer G A, Clare E L, et al. Design and applicability of DNA arrays and DNA barcodes in biodiversity monitoring[J]. BMC Biology, 2007, 5(1):24-24.

    [24]潘程瑩, 胡婧, 張霞,等. 斑腿蝗科Catantopidae七種蝗蟲線粒體COI基因的DNA條形碼研究[J]. 昆蟲分類學報, 2006, 28(2):103-110.

    [25]Hollingsworth P M, Forrest L L, Spouge J L, et al. A DNA barcode for land plants[J]. Proceedings of the National Academy of Sciences of the United States of America, 2009, 106(31):12794-12797.

    [26]Eason R G, Pourmand N, Tongprasit W, et al. Characterization of synthetic DNA bar codes in Saccharomyces cerevisiae gene-deletion strains.[J]. Proceedings of the National Academy of Sciences of the United States of America, 2004, 101(30):11046-11051.

    [27]Druzhinina I S, Kopchinskiy A G, Komon M, et al. An oligonucleotide barcode for species identification in Trichoderma and Hypocrea[J]. Fungal Genetics and Biology, 2005, 42(10):813-828.

    [28]Seifert K A, Samson R A, Dewaard J R, et al. Prospects for fungus identification using CO1 DNA barcodes, with Penicillium as a test case[J]. Proceedings of the National Academy of Sciences of the United States of America, 2007, 104(10):3901-3906.

    [29]Mouhamadou B, Carriconde F, Gryta H, et al. Molecular evolution of mitochondrial ribosomal DNA in the fungal genus Tricholoma : Barcoding implications[J]. Fungal Genetics and Biology, 2008, 45(9):1219-1226.

    [30]Zhao P, Luo J, Zhuang W, et al. Practice towards DNA barcoding of the nectriaceous fungi[J]. Fungal Diversity, 2011, 46(1):183-191.

    [31]Samson R A, Seifert K A, Kuijpers A F A, et al. Phylogenetic analysis of Penicillium subgenus Penicillium using partial 尾-tubulin sequences. Stud Mycol[J]. Studies in Mycology, 2004, 2004(49):175-200.

    [32]Gao T, Chen S. Authentication of the Medicinal Plants in Fabaceae by DNA Barcoding Technique[J]. Planta Medica, 2009, 75(4).doi:10.1055/s-2009-1216451.

    [33]Chen S, Yao H, Han J, et al. Validation of the ITS2 Region as a Novel DNA Barcode for Identifying Medicinal Plant Species[J]. PLOS ONE, 2010, 5(1).doi:10.1371/journal.pone.0008613.

    ConstructionandVerificationofDNABarcodeforHericiumerinaceus

    Huangfu Yongguan1, Yan Baosong2, Zhang Yuexin2, Zong Xianchun3

    (1.Northeast Forestry University/Saline and alkaline land resources and Environment Research Center,Harbin 150040;2.Heilongjiang Forest Products Research Institute,Heilongjiang Provincial Key Laboratory of Non-wood Forest Product Development,Mudanjiang, Heilongjiang 157011;3.Mudanjiang Normal University/College of life science and technology,Mudanjiang 157000)

    2017-05-06

    林業(yè)公益性行業(yè)科研專項經(jīng)費(201404703-3)

    皇甫永冠(1990-),女,博士,從事植物抗逆基因與逆境脅迫應答的分子機理研究,E-mail:595287087@qq.com.

    S759.81

    A

    DOI.:10.13268/j.cnki.fbsic.2017.05.003

    桃花免费在线播放| 91老司机精品| 一级毛片女人18水好多 | 97在线人人人人妻| 久久久国产欧美日韩av| 午夜影院在线不卡| kizo精华| 岛国毛片在线播放| 欧美激情高清一区二区三区| 99国产精品一区二区蜜桃av | 国产精品欧美亚洲77777| 国产精品秋霞免费鲁丝片| 黄色片一级片一级黄色片| 99香蕉大伊视频| 少妇人妻 视频| 狂野欧美激情性xxxx| 久久久精品94久久精品| 精品第一国产精品| 国产99久久九九免费精品| 亚洲av片天天在线观看| 国产精品久久久久成人av| 欧美日韩精品网址| 久久久久久免费高清国产稀缺| 狠狠精品人妻久久久久久综合| 成人国产一区最新在线观看 | 狂野欧美激情性xxxx| 久久久久精品国产欧美久久久 | 亚洲av男天堂| 少妇的丰满在线观看| 丰满人妻熟妇乱又伦精品不卡| 人人妻人人添人人爽欧美一区卜| 少妇猛男粗大的猛烈进出视频| 不卡av一区二区三区| 丝袜脚勾引网站| 亚洲五月色婷婷综合| 韩国高清视频一区二区三区| 老熟女久久久| videos熟女内射| 两个人免费观看高清视频| 啦啦啦 在线观看视频| 狠狠婷婷综合久久久久久88av| 一级,二级,三级黄色视频| 亚洲一区二区三区欧美精品| 嫁个100分男人电影在线观看 | 久久久久久久大尺度免费视频| 精品亚洲成国产av| 国产av一区二区精品久久| 成年人免费黄色播放视频| 欧美激情高清一区二区三区| 中国美女看黄片| 人人妻人人澡人人看| 国产成人啪精品午夜网站| 国产成人免费无遮挡视频| 9热在线视频观看99| 婷婷成人精品国产| tube8黄色片| 日韩精品免费视频一区二区三区| 宅男免费午夜| 国产男女超爽视频在线观看| 91国产中文字幕| 色精品久久人妻99蜜桃| 视频在线观看一区二区三区| 亚洲一卡2卡3卡4卡5卡精品中文| 少妇的丰满在线观看| 制服人妻中文乱码| 99久久99久久久精品蜜桃| 亚洲精品日本国产第一区| 黄色 视频免费看| 精品人妻一区二区三区麻豆| 久久国产精品男人的天堂亚洲| 宅男免费午夜| 你懂的网址亚洲精品在线观看| 日韩大码丰满熟妇| 99久久精品国产亚洲精品| www.熟女人妻精品国产| 亚洲中文日韩欧美视频| 黄色视频在线播放观看不卡| 午夜视频精品福利| 久久狼人影院| 性高湖久久久久久久久免费观看| 人成视频在线观看免费观看| 欧美精品一区二区免费开放| 可以免费在线观看a视频的电影网站| 中国国产av一级| 丰满少妇做爰视频| 激情视频va一区二区三区| 少妇精品久久久久久久| 人人妻人人爽人人添夜夜欢视频| 黄频高清免费视频| 亚洲美女黄色视频免费看| 一本色道久久久久久精品综合| 国产精品一区二区免费欧美 | 亚洲国产欧美网| 欧美 亚洲 国产 日韩一| 激情视频va一区二区三区| 亚洲第一av免费看| 成年人免费黄色播放视频| 一级片免费观看大全| 18禁黄网站禁片午夜丰满| 婷婷丁香在线五月| 国精品久久久久久国模美| 9热在线视频观看99| 亚洲伊人色综图| 国产伦理片在线播放av一区| 成年美女黄网站色视频大全免费| 婷婷色综合大香蕉| 国产人伦9x9x在线观看| 午夜91福利影院| 超碰97精品在线观看| 国产精品久久久av美女十八| 一二三四社区在线视频社区8| 久久精品人人爽人人爽视色| 欧美精品一区二区大全| 免费在线观看黄色视频的| av一本久久久久| 中文字幕精品免费在线观看视频| 飞空精品影院首页| 五月开心婷婷网| 亚洲人成77777在线视频| 只有这里有精品99| 亚洲综合色网址| 亚洲国产欧美网| 超碰97精品在线观看| 午夜福利乱码中文字幕| 波野结衣二区三区在线| 久久精品亚洲av国产电影网| 桃花免费在线播放| 亚洲 欧美一区二区三区| 日韩 亚洲 欧美在线| 久久久国产欧美日韩av| 国产亚洲欧美在线一区二区| 丝袜人妻中文字幕| 成人三级做爰电影| 99久久综合免费| 精品国产超薄肉色丝袜足j| 美女中出高潮动态图| 人体艺术视频欧美日本| 99精品久久久久人妻精品| 国产视频首页在线观看| 曰老女人黄片| 精品人妻一区二区三区麻豆| 日本黄色日本黄色录像| 男女边摸边吃奶| 可以免费在线观看a视频的电影网站| 少妇猛男粗大的猛烈进出视频| 另类精品久久| 久久久欧美国产精品| 在线看a的网站| 亚洲成国产人片在线观看| 一本综合久久免费| 国产黄色视频一区二区在线观看| 欧美 日韩 精品 国产| 91成人精品电影| 久久狼人影院| 亚洲七黄色美女视频| 最新在线观看一区二区三区 | 精品亚洲成a人片在线观看| 看免费成人av毛片| 国产精品成人在线| 男女边吃奶边做爰视频| 亚洲精品国产av成人精品| 男人舔女人的私密视频| 亚洲久久久国产精品| 免费高清在线观看视频在线观看| 国产精品三级大全| 最近手机中文字幕大全| 老汉色∧v一级毛片| 国产精品国产av在线观看| 国产成人精品久久二区二区91| e午夜精品久久久久久久| 啦啦啦在线免费观看视频4| 亚洲欧美激情在线| 青春草视频在线免费观看| 香蕉丝袜av| 两人在一起打扑克的视频| 9色porny在线观看| 国产片特级美女逼逼视频| 午夜免费成人在线视频| 女人久久www免费人成看片| 欧美日韩视频精品一区| 国产精品一区二区在线不卡| 国产精品免费视频内射| 欧美在线黄色| 一级毛片黄色毛片免费观看视频| 日本wwww免费看| 丝袜美足系列| 久久人人97超碰香蕉20202| www.av在线官网国产| 99香蕉大伊视频| 日日夜夜操网爽| 在线亚洲精品国产二区图片欧美| www.999成人在线观看| 18在线观看网站| 亚洲欧洲精品一区二区精品久久久| 亚洲午夜精品一区,二区,三区| 日本欧美视频一区| 狠狠婷婷综合久久久久久88av| 最新在线观看一区二区三区 | 国产日韩欧美在线精品| 人体艺术视频欧美日本| 国产成人91sexporn| 热re99久久国产66热| 国产亚洲av片在线观看秒播厂| 国产精品.久久久| 男人添女人高潮全过程视频| 黄色视频不卡| 女人精品久久久久毛片| 国产片特级美女逼逼视频| 久久久精品免费免费高清| 久久精品亚洲av国产电影网| 高潮久久久久久久久久久不卡| 黄色 视频免费看| 国产精品.久久久| 中文字幕色久视频| 一级黄色大片毛片| 久久鲁丝午夜福利片| 亚洲伊人久久精品综合| 99九九在线精品视频| 国产淫语在线视频| 日韩人妻精品一区2区三区| 在线观看免费午夜福利视频| 香蕉丝袜av| 国产精品亚洲av一区麻豆| 夜夜骑夜夜射夜夜干| 在线观看人妻少妇| 肉色欧美久久久久久久蜜桃| 大码成人一级视频| 欧美精品啪啪一区二区三区 | 黄色视频在线播放观看不卡| 国产淫语在线视频| 欧美精品人与动牲交sv欧美| 男女国产视频网站| 一级黄色大片毛片| 久久久欧美国产精品| 国产黄频视频在线观看| 亚洲精品久久午夜乱码| 男人舔女人的私密视频| 国产av精品麻豆| 精品欧美一区二区三区在线| 久久久久网色| 妹子高潮喷水视频| 国产一区有黄有色的免费视频| 国产免费福利视频在线观看| 大片电影免费在线观看免费| 下体分泌物呈黄色| 亚洲成人国产一区在线观看 | 国产一级毛片在线| 久久鲁丝午夜福利片| 操出白浆在线播放| 精品卡一卡二卡四卡免费| 波多野结衣av一区二区av| 交换朋友夫妻互换小说| 欧美日韩亚洲国产一区二区在线观看 | 久久精品久久久久久噜噜老黄| 丁香六月欧美| 国产野战对白在线观看| 在线看a的网站| 国产高清videossex| 热re99久久精品国产66热6| 在线天堂中文资源库| 国产成人av教育| 一本综合久久免费| 日韩精品免费视频一区二区三区| 丝袜人妻中文字幕| av在线播放精品| 黑人巨大精品欧美一区二区蜜桃| 人体艺术视频欧美日本| av网站在线播放免费| 侵犯人妻中文字幕一二三四区| 男女床上黄色一级片免费看| 女警被强在线播放| 一本大道久久a久久精品| 欧美97在线视频| a级毛片在线看网站| 国产精品一区二区免费欧美 | 男女无遮挡免费网站观看| 国产成人影院久久av| 脱女人内裤的视频| 97在线人人人人妻| cao死你这个sao货| 欧美大码av| 亚洲国产成人一精品久久久| 黑丝袜美女国产一区| 狂野欧美激情性xxxx| 水蜜桃什么品种好| 深夜精品福利| 成人黄色视频免费在线看| 91精品伊人久久大香线蕉| 无遮挡黄片免费观看| 亚洲国产欧美一区二区综合| 一级毛片电影观看| 国产成人欧美| 久久国产精品大桥未久av| 国产男女内射视频| 亚洲精品国产区一区二| 国产高清不卡午夜福利| 又粗又硬又长又爽又黄的视频| 亚洲欧美激情在线| 国产免费现黄频在线看| av网站在线播放免费| 久久久久久免费高清国产稀缺| 国产精品.久久久| 一边亲一边摸免费视频| 妹子高潮喷水视频| 国产成人精品在线电影| 成人三级做爰电影| 91精品国产国语对白视频| 国产成人系列免费观看| 欧美日韩成人在线一区二区| 国产1区2区3区精品| 国产免费又黄又爽又色| 久久精品成人免费网站| 久久 成人 亚洲| 无遮挡黄片免费观看| 国产欧美日韩一区二区三 | 午夜福利乱码中文字幕| 亚洲伊人久久精品综合| 少妇 在线观看| 久久人妻福利社区极品人妻图片 | 日韩,欧美,国产一区二区三区| 91麻豆av在线| 国产黄色视频一区二区在线观看| 五月开心婷婷网| av一本久久久久| 午夜精品国产一区二区电影| 欧美精品人与动牲交sv欧美| 久久精品国产综合久久久| 1024香蕉在线观看| 国产亚洲欧美精品永久| 欧美成人精品欧美一级黄| 99国产综合亚洲精品| 国产精品一区二区免费欧美 | 五月开心婷婷网| 国产精品一区二区在线不卡| videos熟女内射| 成人黄色视频免费在线看| 精品福利观看| 成在线人永久免费视频| 亚洲情色 制服丝袜| 国产日韩欧美在线精品| 一区二区av电影网| 黄色视频在线播放观看不卡| 亚洲成色77777| 亚洲国产欧美在线一区| 超碰97精品在线观看| 亚洲欧美一区二区三区黑人| 精品高清国产在线一区| 欧美老熟妇乱子伦牲交| xxx大片免费视频| 国产欧美日韩一区二区三区在线| 日韩人妻精品一区2区三区| 久久久久久久精品精品| 一个人免费看片子| 成年人黄色毛片网站| 国产欧美日韩综合在线一区二区| 国产激情久久老熟女| 老汉色∧v一级毛片| 在线看a的网站| 人人妻人人添人人爽欧美一区卜| 日韩 欧美 亚洲 中文字幕| 99精品久久久久人妻精品| 欧美日韩精品网址| 国产免费现黄频在线看| 国产高清videossex| 在线观看人妻少妇| 大香蕉久久成人网| 国产精品久久久久成人av| 亚洲久久久国产精品| 亚洲精品国产一区二区精华液| 黄色 视频免费看| 国产福利在线免费观看视频| xxxhd国产人妻xxx| 一个人免费看片子| 一区二区三区精品91| av国产久精品久网站免费入址| 日韩电影二区| 亚洲伊人久久精品综合| 亚洲成人免费电影在线观看 | 亚洲精品国产一区二区精华液| 一区二区三区激情视频| 国产熟女欧美一区二区| 啦啦啦在线观看免费高清www| 一区二区三区激情视频| 久久精品国产亚洲av高清一级| 真人做人爱边吃奶动态| 国产高清视频在线播放一区 | 老司机午夜十八禁免费视频| 在线 av 中文字幕| 性高湖久久久久久久久免费观看| 精品一区在线观看国产| 国产野战对白在线观看| 欧美人与性动交α欧美精品济南到| 一边摸一边抽搐一进一出视频| 亚洲精品美女久久av网站| 亚洲欧洲日产国产| 精品人妻一区二区三区麻豆| 国产一区二区三区av在线| 国产老妇伦熟女老妇高清| www.av在线官网国产| 亚洲精品国产色婷婷电影| 精品亚洲成a人片在线观看| 亚洲av国产av综合av卡| 色婷婷久久久亚洲欧美| 99re6热这里在线精品视频| 男女边吃奶边做爰视频| 国产亚洲av片在线观看秒播厂| 久久九九热精品免费| 久久久久网色| 精品久久久久久久毛片微露脸 | 精品国产一区二区久久| 一级黄片播放器| 成在线人永久免费视频| 黄色怎么调成土黄色| 欧美黑人精品巨大| 99久久精品国产亚洲精品| 婷婷色av中文字幕| 日韩av不卡免费在线播放| 丝袜美足系列| 久久亚洲国产成人精品v| 精品久久蜜臀av无| 乱人伦中国视频| 国产一区二区三区av在线| 黄色视频在线播放观看不卡| 美女脱内裤让男人舔精品视频| 丝袜美足系列| 老汉色av国产亚洲站长工具| 欧美另类一区| 国产日韩欧美在线精品| 精品一区二区三区四区五区乱码 | 色94色欧美一区二区| 国产精品一区二区免费欧美 | 亚洲,欧美,日韩| 一级毛片电影观看| 国产免费又黄又爽又色| 99久久精品国产亚洲精品| 在线亚洲精品国产二区图片欧美| 亚洲人成电影观看| 9热在线视频观看99| 免费女性裸体啪啪无遮挡网站| 丝袜美足系列| 成年人免费黄色播放视频| 欧美久久黑人一区二区| 欧美黑人精品巨大| 首页视频小说图片口味搜索 | 操美女的视频在线观看| 亚洲专区中文字幕在线| 成人手机av| 久久久久国产精品人妻一区二区| 亚洲国产精品一区二区三区在线| 精品国产国语对白av| av有码第一页| √禁漫天堂资源中文www| 老熟女久久久| 在线观看免费午夜福利视频| 超碰97精品在线观看| 亚洲人成电影观看| a 毛片基地| 天堂中文最新版在线下载| 久久99热这里只频精品6学生| 午夜福利乱码中文字幕| 国产97色在线日韩免费| 亚洲精品成人av观看孕妇| 99国产综合亚洲精品| 精品卡一卡二卡四卡免费| 成年人免费黄色播放视频| 欧美久久黑人一区二区| 久久精品aⅴ一区二区三区四区| 国产精品人妻久久久影院| 女人高潮潮喷娇喘18禁视频| 操出白浆在线播放| 国产成人欧美| 亚洲色图 男人天堂 中文字幕| 国产精品国产三级国产专区5o| 母亲3免费完整高清在线观看| 久久久久网色| 大话2 男鬼变身卡| 中文字幕精品免费在线观看视频| 大陆偷拍与自拍| 亚洲中文日韩欧美视频| 国产亚洲欧美在线一区二区| 中文字幕av电影在线播放| 欧美人与善性xxx| 制服人妻中文乱码| 亚洲成人国产一区在线观看 | 国产精品国产av在线观看| 国产精品三级大全| 午夜福利乱码中文字幕| 又粗又硬又长又爽又黄的视频| 亚洲成人免费av在线播放| 一级毛片 在线播放| videos熟女内射| 久久这里只有精品19| av线在线观看网站| 国产成人一区二区三区免费视频网站 | 久久久久久久国产电影| 在线观看免费午夜福利视频| 国产高清视频在线播放一区 | 国产精品人妻久久久影院| 亚洲国产毛片av蜜桃av| 精品久久蜜臀av无| 免费少妇av软件| 又粗又硬又长又爽又黄的视频| 天天影视国产精品| 国产男人的电影天堂91| 久久狼人影院| 国产精品三级大全| 久久毛片免费看一区二区三区| 精品一区在线观看国产| 国产精品国产av在线观看| av一本久久久久| 1024香蕉在线观看| 9色porny在线观看| 99国产精品免费福利视频| 啦啦啦在线免费观看视频4| 国产三级黄色录像| 亚洲国产欧美网| 99国产综合亚洲精品| 欧美人与性动交α欧美软件| 51午夜福利影视在线观看| 国产老妇伦熟女老妇高清| 亚洲欧美一区二区三区黑人| 欧美 亚洲 国产 日韩一| 中文字幕人妻丝袜制服| 一本久久精品| 免费观看人在逋| 国产老妇伦熟女老妇高清| 电影成人av| 亚洲欧美日韩另类电影网站| 日韩一本色道免费dvd| 国产精品免费视频内射| 婷婷色麻豆天堂久久| 亚洲熟女毛片儿| 香蕉国产在线看| 亚洲国产成人一精品久久久| 人人妻人人澡人人看| 精品人妻一区二区三区麻豆| 一本—道久久a久久精品蜜桃钙片| av在线老鸭窝| 深夜精品福利| 欧美国产精品一级二级三级| 蜜桃在线观看..| 亚洲欧美日韩高清在线视频 | 精品高清国产在线一区| 色94色欧美一区二区| 久久精品国产综合久久久| 欧美日韩综合久久久久久| 国产在视频线精品| 大型av网站在线播放| 97精品久久久久久久久久精品| 在线观看免费日韩欧美大片| 黄色 视频免费看| 少妇人妻 视频| 在线观看人妻少妇| 欧美日韩亚洲国产一区二区在线观看 | av电影中文网址| 亚洲av日韩在线播放| 亚洲精品一卡2卡三卡4卡5卡 | 成在线人永久免费视频| 不卡av一区二区三区| 亚洲国产欧美一区二区综合| 亚洲熟女毛片儿| 欧美日韩福利视频一区二区| 人妻 亚洲 视频| 午夜91福利影院| 青草久久国产| 久久狼人影院| 精品少妇内射三级| 亚洲五月婷婷丁香| 日韩视频在线欧美| 亚洲国产中文字幕在线视频| 亚洲视频免费观看视频| 久久99热这里只频精品6学生| 日本黄色日本黄色录像| 久久久久精品国产欧美久久久 | 男人爽女人下面视频在线观看| 免费高清在线观看日韩| 欧美乱码精品一区二区三区| 最近最新中文字幕大全免费视频 | 亚洲av国产av综合av卡| 亚洲精品中文字幕在线视频| 99热国产这里只有精品6| av在线app专区| 久久久国产一区二区| 尾随美女入室| 日日爽夜夜爽网站| 免费观看人在逋| 日韩 亚洲 欧美在线| 成年人黄色毛片网站| 韩国精品一区二区三区| xxx大片免费视频| 亚洲精品一二三| 免费一级毛片在线播放高清视频 | 久久 成人 亚洲| 免费观看人在逋| 久久久久国产精品人妻一区二区| 男人操女人黄网站| 在线观看人妻少妇| 欧美少妇被猛烈插入视频| 久久精品国产亚洲av涩爱| 啦啦啦在线免费观看视频4| 精品国产乱码久久久久久男人| 国产片特级美女逼逼视频| 国产精品偷伦视频观看了| 好男人电影高清在线观看| 色综合欧美亚洲国产小说| 黄片小视频在线播放| 热re99久久国产66热| 99国产精品一区二区蜜桃av | 免费少妇av软件| 中国国产av一级| 女性生殖器流出的白浆| 欧美精品啪啪一区二区三区 | 亚洲成色77777| 国产精品久久久久久精品电影小说| 久久九九热精品免费|