一種太白貝母的分子生物學鑒定方法
【專利摘要】本發(fā)明公開一種太白貝母的分子生物學鑒定方法，該方法包括提取待鑒定藥材的葉綠體基因組DNA并純化，進行測序與序列拼接，得到完整的葉綠體基因組序列，將該序列與太白貝母的葉綠體基因組序列進行比對，從而實現(xiàn)對太白貝母進行快速和準確的鑒定。
【專利說明】—種太白貝母的分子生物學鑒定方法
【技術領域】
[0001]本發(fā)明屬于中藥鑒定【技術領域】，具體公開了一種太白貝母的分子生物學鑒定方法。
【背景技術】
[0002]太白貝母,又稱太貝、尖貝、秦貝，為百合科Liliaceae貝母屬Fritillaria L.植物太白貝母F.taipaiensis P.Y.Li的干燥鱗莖，主治肺熱燥咳、干咳少痰、陰虛勞嗽、咯痰帶血，為2010版《中華人民共和國藥典》收錄品種。太白貝母和瓦布貝母在藥典中均作為川貝母入藥，加之川貝母、梭砂貝母、甘肅貝母和暗紫貝母，川貝母基源植物達到6個，但其外形極為相似，僅有微小變異，傳統(tǒng)經(jīng)典分類學主要依據(jù)花特征進行鑒定，但由于花器官特征性差別小及過渡類型較多，定性非定量的方法難以客觀的進行區(qū)分。以川貝母入藥的源植物與貝母屬其它常見藥典載品種同樣難以區(qū)分，如平貝母、浙貝母和伊犁貝母等。生物條形碼對物種的鑒定科學、快速、有效，易于形成標準，對貝母屬藥用植物的鑒別具有重要意義。
[0003]葉綠體基因作為DNA Barcoding標準序列的地位已經(jīng)得到公認，但是在DNA條形碼鑒定研究工作中發(fā)現(xiàn)，盡管該方法具有操作簡單、可重復性強、適用于多個分類階元等諸多優(yōu)點，但DNA條形碼長度多介于300-700bp之間，有時包含序列的變異信息不充足，在區(qū)分和鑒定部分藥用植物(如貝母屬)種內(nèi)和種間樣品時存在困難；同時植物基因組的高變異性導致在植物中至今還未找到公認的通用DNA條形碼序列，以上從理論上證明一些常用的DNA條形碼在通用性上存在局限性。葉綠體基因組在植物中廣泛存在，長度達到120-160kb，與DNA條形碼相比，具有更好的通用性和更強的分辨力。但是由于測序模板提取還存在一定的困難，在沒有近緣物種的序列時對測序數(shù)據(jù)進行拼接組裝仍然是復雜的。因此，現(xiàn)有技術中認為采用葉綠體基因組進行物種的鑒定目前還不是一個好的選擇。本發(fā)明通過對太白貝母葉綠體基因組的獲得及同屬近緣物種的分析，首次提出可以利用葉綠體全基因組對太白貝母進行快速和準確的鑒定。

【發(fā)明內(nèi)容】

[0004]本發(fā)明要解決的技術問題是以葉綠體基因組全序列作為太白貝母的分子鑒定標記，對太白貝母進行快速和準確的鑒定。
[0005]為解決上述技術問題，本發(fā)明采用下述技術方案:
[0006]一種太白貝母的分子生物學鑒定方法，該方法包括提取待鑒定藥材的葉綠體基因組DNA并純化，進行測序與序列拼接，得到完整的葉綠體基因組序列，將該序列與太白貝母的葉綠體基因組序 列進行比對，得到分析鑒定結果，其中，太白貝母葉綠體基因組序列的NCBI 序列號為 KC543997。
[0007]上述方法中，所述提取待鑒定藥材的葉綠體基因組DNA并純化的具體方法為:
[0008](I)取新鮮幼嫩葉片，4°C暗處理24h，洗凈，去葉脈，撕成碎片，加4°C預冷的緩沖液A置于勻漿器中，充分勻漿，然后過濾，收集濾液，棄渣；[0009](2)將濾液，于4°C，1500g，離心15min，棄上清；
[0010](3)用緩沖液I清洗沉淀；
[0011](4)用洗滌緩沖液懸浮沉淀，在一個新的離心管中加入30%的蔗糖溶液，將沉淀懸浮液覆蓋在30%的蔗糖溶液的表面，于4°C，2500g，離心15min，沉淀葉綠體；
[0012](5)沉淀重懸于洗滌緩沖液中，于4°C，100g離心5min，上清轉(zhuǎn)移到新管并在380g離心10min ；
[0013](6)沉淀重懸于洗滌緩沖液中，將該溶液覆蓋于蔗糖梯度的頂部，于4°C，20000g離心70min ;其中，鹿糖梯度由7mL30%和16.5~18mL52%的鹿糖溶液組成；
[0014](7 )去掉上層，從52%和30%的蔗糖溶液的中間相，將葉綠體DNA收集到離心管中，然后加入洗滌緩沖液，于4°C，1500g離心15min沉淀葉綠體DNA ；
[0015](8)用緩沖液2清洗沉淀，于4°C，1500g，離心15min，取沉淀，加入緩沖液2和濃度為10mg/mL的蛋白酶K，室溫靜置2min，加入1/5體積的裂解液，顛倒混勻，然后在50°C保溫至少Ih，冰上冷卻；
[0016](9)轉(zhuǎn)移到滅菌離心管中，加入等體積的酚:氯仿:異戊醇，顛倒混勻，靜置5min，2600g室溫離心20min ;其中，酚:氯仿:異戊醇的體積比為25:24:1 ；
[0017](10)上清轉(zhuǎn)移到新的離心管中，加入2/5體積的5mol ^r1NaCl, 2/5體積的預熱至650C的CTAB/NaCl，并68°C加熱15分鐘，加入等體積的氯仿:異戊醇，顛倒混勻，靜置5min，2600g室溫離心20min,取上清；其中，氯仿:異戍醇的體積比為24:1 ;
[0018](11)重復(10)`
[0019](12)加等體積的氯仿:異戍醇，顛倒混勻，靜置5min, 2600g室溫離心20min ;其中，氯仿:異戊醇的體積比為24:1 ;
[0020](13)取上清，加入0.1倍體積的3mol.L4NaAc和0.8倍體積的預冷異丙醇或是加入0.1倍體積的3mol.T1NaAc和2倍體積的無水乙醇，_20°C放置至少20min沉淀葉綠體 DNA ；
[0021](14)4°C,13500g離心20min，將沉淀用70%的乙醇室溫清洗兩次，13500g離心IOmin,風干，用含有RNA酶的TE重懸，37°C保溫lh，得到葉綠體DNA ；
[0022]其中，所述緩沖液A 含有 50mM Tris,25mM EDTA, 1.25M NaCl，0.25M 抗壞血酸，1%PVP40,0.1%BSA,ImM DTT，ρΗ3.6 ；
[0023]所述緩沖液I 含有 1.25Μ NaCl，ρΗ8.0 的 50mM Tris-HCl, ρΗ8.0 的 7mMEDTA,5%PVP40,0.1%BSA,ImM DTT ；
[0024]所述洗滌緩沖液含有0.35M山梨醇，pH8.0的50mM Tris-HCl, pH8.0的25mMEDTA ；
[0025]所述緩沖液2 含有 ρΗ8.0 的 IOmM Tris-HCl, ρΗ8.0 的 5mM EDTA ；
[0026]所述裂解液含有20%十二烷基肌氨酸鈉，pH8.0的50mM Tris-HCl, pH8.0的25mMEDTA0
[0027]進一步地，葉綠體基因組DNA測序與序列拼接的具體方法為:采用第二代高通量測序技術，通過Roche454GS FLX高通量測序儀完成葉綠體基因組的測序和序列拼接，獲得葉綠體全基因組序列。有時，拼接后的葉綠體全基因組序列上可能留有少量的空隙，根據(jù)空隙兩邊的序列設計引物，進行PCR擴增，并用Sanger法測序，填充空隙后得到完整的葉綠體基因組序列。
[0028]本申請通過將待鑒定藥材的葉綠體基因組全序列與GenBank中的太白貝母中葉綠體全基因的序列進行比對，更加準確的鑒定藥材，若其基因序列與之相一致，即同源性達到98%以上，該樣品就是太白貝母，否則就不是該植物。
[0029]本申請還對太白貝母葉綠體全序列進行基因組注釋(采用在線注釋軟件DOGMA (http: //dogma, ccbb.utexas.edu/)),并繪制了葉綠體全基因組物理圖譜(采用在線繪圖工具(http://wolfe.gen.tcd.1e/GenomeVx/))。
[0030]本發(fā)明的有益效果如下:
[0031 ] (I)與總DNA提取方法相比較而言，葉綠體DNA提取方法研究較少，且目前沒有廉價、高效的商品試劑盒。針對一個物種進行優(yōu)化提取流程往往需要數(shù)天或數(shù)十天，阻礙了科研項目的推進。應用文獻資料的方法進行提取葉綠體DNA往往需要30-50個小時，遠遠超過常規(guī)DNA1-2個小時的提取時間，且產(chǎn)量不穩(wěn)定。本發(fā)明提供的葉綠體基因組的提取方法改良了 cpDNA提取流程，太白貝母cpDNA提取效率為206.0 μ g/100g，0D260/0D280為1.94。
[0032](2)本發(fā)明首次提出了葉綠體基因組更適合作為太白貝母和其它貝母屬植物物種的分子鑒定標記。通過對待測樣品葉綠體全基因組全序列進行測定，可以在分子水平對太白貝母進行鑒別，還可以有效地區(qū)分太白貝母及同屬近緣植物。對于傳統(tǒng)分類學中難以區(qū)分的類群，采用葉綠體基因組序列分析的方法可以拋開形態(tài)相似的假象，從基因水平上提供新的分類依據(jù)。
【專利附圖】

【附圖說明】
[0033]下面結合附圖對本發(fā)明的【具體實施方式】作進一步詳細的說明。
[0034]圖1為太白貝母葉綠體全基因組物理圖譜。
【具體實施方式】
[0035]為了更清楚地說明本發(fā)明，下面結合優(yōu)選實施例和附圖對本發(fā)明做進一步的說明。本領域技術人員應當理解，下面所具體描述的內(nèi)容是說明性的而非限制性的，不應以此限制本發(fā)明的保護范圍。
[0036]實驗試劑與實驗儀器
[0037](1)緩沖液 A (pH=3.6):50mM Tris,25mM EDTA, 1.25M NaCl，0.25M 抗壞血酸(VC)，1%PVP40 (聚乙稀吡咯烷酮)，0.1%BSA, ImM DTT ；
[0038](2)緩沖液 I:1.25M NaCl,50mM Tris-HCl (ρΗ8.0)，7mM EDTA (ρΗ8.0),5%PVP40,
0.1%BSA,ImM DTT ；
[0039](3)緩沖液 2:10mM Tris-HCl (pH8.0),5mM EDTA (pH8.0)；
[0040](4)洗滌緩沖液:0.35M 山梨醇，50mM Tris-HCl (ρΗ8.0)，25mMEDTA (ρΗ8.0)；
[0041](5) 52% 的蔗糖:52% 的蔗糖溶于 50mM 的 Tris-HCl (pH8.0),25mMEDTA (ρΗ8.0)；
[0042](6) 30% 的蔗糖:30% 的蔗糖溶于 50mM 的 Tris-HCl (ρΗ8.0)，25mMEDTA (ρΗ8.0)；
[0043](7)裂解液:20% 十二烷基肌氨酸鈉，50mM Tris-HCl (pH8.0),25mMEDTA (ρΗ8.0); [0044](8 ) CTAB/NaCl 溶液:10%CTAB，0.7M NaCl。
[0045](9)瓊脂糖:Biowest Agrose (Spain)。[0046](10)含有 RNA 酶的 TE:200uL TE 加入 luL50mg/mL 的 RNAase
[0047]其它試劑耗材:10ml注射器和平頭針、異丙醇、漂白劑、無水乙醇、70%乙醇，以上試劑組分均購自北京碧橙藍生物科技有限責任公司？。
[0048]天根植物組織DNA提取試劑盒:目錄號DP305，天根生化科技(北京)有限公司(中國)。
[0049]PCR材料:TaKaRa Ex Taq?,Code drrooia,寶生物工程(大連)有限公司(中國)。
[0050]脫脂棉紗布:非滅菌，山東省曹縣華魯衛(wèi)生材料有限公司(中國)；MiraCl0th:cat#475855, CalbiOChem?，Merck (德國);
[0051]各種規(guī)格離心管:海門市三和新華玻璃實驗儀器廠。
[0052]實驗儀器
[0053]低溫高速離心機:SIGMALaboratory Centrifuges3K3O (德國)；
[0054]制冰機:三洋(日本)；
[0055]臺式高速冷凍大容量離心機(水平轉(zhuǎn)子):Eppendorf Centrifuge5810R (德國)；
[0056]高速勻漿機:IKA?T18 basic MLTRA-TURRAX^ (德國)；
[0057]超高速低溫離心機:Beckm`anCOMLTER 0ptima?L-80XPMLtracentrifUge (美國)；
[0058]凝膠成像儀:Biorad (美國)；
[0059]電熱恒溫水浴鍋:DK_S24型，上海森信實驗儀器有限公司(中國)；
[0060]電子天平:METTLER TOLEDO PL203，梅特勒-托利多儀器(上海)有限公司(中國)；
[0061]電泳儀:DYY_8C型，北京市六一儀器廠(中國)；
[0062]微量分光光度計:ThermoScientific NanoDrop2000 (美國)；
[0063]PCR 儀:Applied Biosystems? GeneAmp? pcr System9700 (美國);
[0064]Applied Biosystems2720Thermal Cycler (美國)；
[0065]生物安全柜:蘇凈安泰(中國)；
[0066]渦旋儀:IKA' MS3 (德國);
[0067]高壓滅菌鍋:三洋(日本)。
[0068]實施例1太白貝母葉綠體基因組的獲得
[0069]( I)太白貝母葉綠體基因組DNA提取與純化
[0070]取太白貝母新鮮幼嫩葉片(太白貝母的新鮮葉片采自重慶市巫溪縣文峰店紅池壩(緯度:31.4°，經(jīng)度:109.6° )，經(jīng)中國醫(yī)學科學院藥用植物研究所鑒定為太白貝母Fritillaria.taipaiensis P.Y.Li,),采取改良的鹿糖密度梯度離心法分離提取并純化葉綠體基因組(cpDNA),步驟如下:
[0071]①取新鮮幼嫩葉片100g，4°C暗處理24h，洗凈，去葉脈，撕成lcm2左右大小的碎片，加500mL4°C預冷的緩沖液A置于勻漿器中，使葉片基本浸沒，充分勻漿。然后用4層紗布過濾I次，4層Miracloth濾布過濾I次,收集濾液,棄洛。
[0072]②將濾液分裝到50ml離心管中，4°C下1500g, 15min離心,棄上清。[0073]③在沉淀中加入20mL緩沖液1，用玻璃棒充分懸浮沉淀，4°C下1500g離心15min，棄上清；20mL緩沖液I重懸沉淀一次，離心，以徹底清洗沉淀。
[0074]④分別用20mL洗滌緩沖液懸浮各管沉淀。在一個新的50mL離心管中加入20mL30%的蔗糖溶液，將20mL的沉淀懸浮液小心覆蓋在蔗糖溶液的表面，于4°C，2500g水平離心15min,沉淀葉綠體。將離心后的沉淀合并,再次重懸于20mL洗滌緩沖液，并重復小心覆蓋在20mL30%的蔗糖溶液上，2500g水平離心10min。將沉淀合并，重懸于20mL洗滌緩沖液，再次小心覆蓋在20mL30%的蔗糖溶液上，2500g水平離心5min。
[0075]⑤最后的沉淀物重懸于40mL洗漆緩沖液，4°C, 100g離心5min。上清轉(zhuǎn)移到新管并在380g離心10min。如果使用了大量的起始材料,進行一次850g的第三輪離心。如果仍然沒有離心下來，則可適當加大轉(zhuǎn)數(shù)和延長離心時間。
[0076]⑥將沉淀重懸于5mL洗滌緩沖液，將此葉綠體溶液小心覆蓋于蔗糖梯度的頂部(蔗糖梯度由7mL30%和16.5~18mL52%的蔗糖溶液組成，提前配置，4°C放置至少一個小時)，立即上水平轉(zhuǎn)子離心機，4°C下20000g離心lhlOmin。
[0077]⑦去掉上層，從52%和30%的蔗糖溶液的中間相，將葉綠體DNA收集到50mL的離心管中，然后用洗滌緩沖液稀釋到45mL，取1.2mL4°C保存?zhèn)溆谩?5mL管子在4°C，1500g離心15min沉淀葉綠體DNA。
[0078]⑧沉淀在IOmL的緩沖液2中重懸兩次(1500g，15min，4°C離心)，取沉淀，加入2mL緩沖液2和200uL蛋白酶K (10mg/mL)，室溫靜置2min，加入1/5體積(440uL)的裂解液，輕微顛倒管子混合數(shù)次，然后在50°C保溫至少Ih (—般2~3h)，中間溫和震蕩混合，冰上冷卻。
[0079]⑨轉(zhuǎn)移到2mL滅菌離心管中，加入等體積的酚:氯仿:異戊醇(25:24:1 )，顛倒混勻，靜置5min,2600g室溫離心20min。
[0080]⑩上清轉(zhuǎn)移到新2mL離心管中，加入2/5體積的5mol *^01,2/5體積的預熱至65°C的CTAB/NaCl，并68°C加熱15分鐘。加入等體積的氯仿:異戊醇(24:1)，顛倒混勻，靜置5min，2600g室溫離心20min，取上清。
[0081]O重復⑩
[0082]⑩取上清，加等體積的氯仿:異戊醇(24:1)，顛倒混勻，靜置5min，2600g室溫離心 20min。
[0083]O取上清，加入0.1倍體積的3mol.L-1NaAc70.8倍體積的預冷異丙醇(或是0.1倍體積的3mol -T1NaAc, 2倍體積的無水乙醇)，沉淀葉綠體DNA (_20°C放置至少20min，一般過夜)。
[0084]4 °C, 13500g離心20min，將沉淀用70%的乙醇室溫清洗兩次，13500g離心lOmin，風干，用含有RNA酶的TE重懸至30~IOOuL (—般50uL)，37°C保溫lh，0.8%的瓊脂糖凝膠電泳檢測DNA。所得到的葉綠體DNA可以在4°C短期保存，-20°C長期保存。利用NanoDrop2000微量分光光度計檢測cpDNA的濃度并用0.8%的瓊脂糖電泳檢測cpDNA的質(zhì)量。米用 invitrogen 生化 DNA 產(chǎn)物純化試劑盒(Purelink Genomic DNA Kit, K182001)對所提取的cpDNA進一步純化以符合高通量測序要求。
[0085]經(jīng)改良的cpDNA提取流程，太白貝母cpDNA提取效率為206.0 μ g/100g, 260/280為 1.94。
[0086](2)測序和序列拼接
[0087]經(jīng)純化后的太白貝母cpDNA樣品采用第二代高通量測序平臺Roche454GS FLXTitanium進行測序。先將雙鏈DNA打斷為SOObp的小片段，在片段兩端添加特異性銜接子，使之變性為單鏈連接到磁珠上，經(jīng)Emulsion PCR對樣本cpDNA進行擴增，富集到PicoTiterPlate 板,上機測序。具體方法可以參見 Margulies M, Egholm M, Altman WE, etal.Genome sequencing in microfabricated high-density picolitre reactors[J].Nature, 2005，437:3767380 中的記載。
[0088](3)葉綠體序列中空隙的填充
[0089]測序產(chǎn)生的數(shù)據(jù)去掉低質(zhì)量序列后，采用454測序儀自帶Newbler軟件進行denovo拼裝。采用454FLX測序平臺對太白貝母葉綠體基因組測序共獲得41296個CCS序列(14.5Mbp)，覆蓋度約76倍。拼接后獲得45條contig，最長的為42，767bp。其中20條contig的覆蓋度較高(>=20倍)，總長126.398bp。對這些contig進一步拼接,并確定其位置關系，最終獲得完整太白貝母葉綠體基因組序列。
[0090](4)基因組序列 注釋和提交
[0091]采用在線注釋軟件DOGMA (http://dogma, ccbb.utexas.edu/)對太白貝母葉綠體全序列進行基因組注釋。根據(jù)起始密碼子和終止密碼子序列手工調(diào)整DOGMA初步注釋的編碼蛋白基因范圍，將注釋的結果用NCBI網(wǎng)站指定的Seqin軟件提交太白貝母cpDNA序列數(shù)據(jù)和注釋信息，生成太白貝母序列的sqn文件。將太白貝母序列的Genbank格式文件用GenomeVx在線繪圖工具(http://wolfe.gen.tcd.1e/GenomeVx/)繪制太白貝母葉綠體全基因組物理圖譜(圖1)。
[0092]太白貝母葉綠體DNA(chloroplast DNA, cpDNA)為環(huán)狀分子結構，其大小為151693bp(NCBI序列號KC543997)。一對反向重復區(qū)(IRs)長26335bp，大單拷貝區(qū)(LSC)和小單拷貝區(qū)(SSC)分別為81471bp、17552bp。太白貝母cpDNA編碼132個基因，包括86個編碼蛋白基因，8個rRNA基因和38個tRNA基因。蛋白編碼區(qū)占整個基因組的44.88%，rRNA占5.97%, tRNA占6.61%，基因間隔區(qū)和內(nèi)含子占42.54%。其中18個基因(6個tRNA和12個蛋白編碼基因)包含內(nèi)含子。太白貝母葉綠體基因組的GC含量為42.49%，其IR區(qū)GC含量(42.49%)高于LSC (34.82%)和SSC (30.36%)，IR區(qū)的高GC含量是由于IR區(qū)包含4個 rRNA 基因:rrnl6> rrn23> rrn4.5、rrn5 (GC 含量 55.07%)。
[0093]實施例2基于葉綠體基因組全序列進行太白貝母鑒定的可行性
[0094]( I)基于單基因序列的太白貝母和川貝母生物條形碼鑒定
[0095]根據(jù)現(xiàn)有技術中已知的方法提取藥材的基因組DNA，利用葉綠體trnh-psbA序列進行物種鑒定。trnh-psbA區(qū)擴增使用的引物為:
[0096]上游引物:5，-GTTATGCATGACGTAATGCTC-3’
[0097]下游引物:5’ -CGCGCATGGTGATTCACAATCC-3，。
[0098]PCR擴增條件:94 V預變性5min，94 V變性lmin，54 °C退火lmin，72 V延伸1.5min,30個循環(huán)，最后72°C延伸7min。
[0099]測序后，川貝母trnH-psbA序列:
[0100]ACAATTTCCCTCTAGACCTAGCTGCTATTGAAGTTCCATCTACAAATGGATAAGACTTACATTCTGTCTTAGTGTATACGAATCCTTGATGGAGCAATATCTTTGTATCTTGCTAAAGATATTGCTCCATCAAGGATTTTTTTTTA
TTGAACGTGTCACAGGGGATTACTCCTTTATTACATTACATTTTTAAAATTGGCATTCTATGCCCAATATATCGATC
TAAGTATGAAGGTAAGAATAAATACAATAATGATGAATGGAAAAATAAAAAAGCTAAAATCCTTTAGCTAGATAAGG
GGCGGATGTAGCCAAGTGGATCAAGGCAGT
[0101 ] 太白貝母trnH-psbA序列:
[0102]ACAATTTCCCTCTAGACCTAGCTGCTATTGAAGTTCCATCTACAAATGGATAAGACTTACATTCTGTCTTAGTGTATACGAATCCTTGATGGAGCAATATCTTTGTATCTTGCTAAAGATATTGCTCCATCAAGGATTTTTTTTTATTGAACGTGTCACAGGGGATTACTCCTTTATTACATTACATTTTTAAAATTGGCATTCTATGCCCAATATATCGATCTAAGTATGAAGGTAAGAATAAATACAATAATGATGAATGGAAAAATAAAAAAGCTAAAATCCTTTAGCTAGATAAGGGGCGGATGTAGCCAAGTGGATCAAGGCAGT
[0103]從目前對trnH-psbA序列研究的結果看,該序列在對川貝母和太白貝母分析時，序列相似度為100%，提示我們兩個物種在進化時親緣關系較近，對太白貝母采用單基因進行分子鑒定存在較大困難。
[0104](2)基于葉綠體基因組全序列對太白貝母進行鑒定
[0105]a高通量測序及數(shù)據(jù)提取
[0106]根據(jù)實施例1中太白貝母的高通量測序及注釋，得到86個蛋白編碼序列?；谠?6個蛋白基因序列，獲得其它近緣物種(表1) 14個樣本的86個共有基因的序列。
[0107]b數(shù)據(jù)處理
[0108]對14個樣本的共有基因分別進行序列比對，組成序列矩陣，進行wilcoxon twosample test種間種內(nèi)的遺傳距離分析,應用paup軟件并采用blast和distance的方法進行不同分類級別(genus, species, population)的比較分析。為驗證葉綠體基因組序列對貝母屬近緣物種的鑒定能力，采用國際上認可的最優(yōu)單基因組合與基因組序列進行對比。
[0109]c結果分析
[0110]在葉綠體全基因組序列的基礎上，提取所有樣本的共有蛋白編碼基因序列進行太白貝母和同屬近緣物種的核苷酸序列差異分析，以獲得的遺傳差異作為比較對象，開展物種鑒定(結果見表2和表3)。
[0111]從種間種內(nèi)分析看，以葉綠體基因組作為鑒定標記的遺傳距離是國際條形碼組織推薦序列matK+rbcL的2.15倍，葉綠體基因組平均種間差異是最小種間差異的1.46倍，matK+rbcL與平均距離幾乎一樣(1.05倍)。基因組水平，平均種間差異是種內(nèi)差異的16.2倍，matK+rbcL只有4倍，因此葉綠體基因組更適合作為分子鑒定標記。
[0112]從不同分類級別上看，葉綠體基因組無論是blast還是distance統(tǒng)計方法,鑒定效率均為100%，matK+rbcL為78.6，為進一步驗證葉綠體基因組的鑒定效力，在matK+rbcL的基礎上，我們增加了公認的具有更大變異位點的非編碼區(qū)psbA-trnH，分析結果顯示兩種統(tǒng)計方法平均鑒定效率為85.8，增加不顯著，這與國際條形碼組織的結論相一致，即多基因組合不能明顯提高DNA條形碼的鑒定能力。在種下等級，葉綠體基因組的鑒定優(yōu)勢更加明顯，兩種方法的鑒定效率分別為92.9%和100%,平均比Matk+rbcL和Matk+rbcL+psbA-trnH分別高出50.05%和28.6%。因此，葉綠體基因組適合作為貝母屬近緣物種的鑒定標記，可以作為太白貝母和其它貝母屬植物物種的分子條碼。
[0113]表1實驗材料[0114]
【權利要求】
1.一種太白貝母的分子生物學鑒定方法，其特征在于，該方法包括提取待鑒定藥材的葉綠體基因組DNA并純化，進行測序與序列拼接，得到完整的葉綠體基因組序列，將該序列與太白貝母的葉綠體基因組序列進行比對，得到分析鑒定結果。
2.根據(jù)權利要求1所述的太白貝母的分子生物學鑒定方法，其特征在于，太白貝母葉綠體基因組序列的NCBI序列號為KC543997。
3.根據(jù)權利要求1或2所述的太白貝母的分子生物學鑒定方法，其特征在于，提取待鑒定藥材的葉綠體基因組DNA并純化的具體方法為: (1)取新鮮幼嫩葉片，4°C暗處理24h，洗凈，去葉脈，撕成碎片，加4°C預冷的緩沖液A置于勻漿器中，充分勻漿，然后過濾，收集濾液，棄渣； (2)將濾液，于41:，1500g，離心15min，棄上清； (3)用緩沖液1清洗沉淀； (4)用洗滌緩沖液懸浮沉淀，在一個新的離心管中加入30%的蔗糖溶液，將沉淀懸浮液覆蓋在30%的蔗糖溶液的表面，于4°C，2500g,離心15min，沉淀葉綠體； (5)沉淀重懸于洗滌緩沖液中，于4°C，100g離心5min，上清轉(zhuǎn)移到新管并在380g離心10min ； (6)沉淀重懸于洗滌緩沖液中，將該溶液覆蓋于蔗糖梯度的頂部，于4°C，20000g離心70min ;其中，蔗糖梯度由7mL30%和16.5~18mL52%的蔗糖溶液組成； (7)去掉上層，從52%和30%的蔗糖溶液的中間相，將葉綠體DNA收集到離心管中，然后加入洗滌緩沖液，于4°C，1500g離心15min沉淀葉綠體DNA ； (8)用緩沖液2清洗沉淀，于4°C，1500g，離心15min，取沉淀，加入緩沖液2和濃度為10mg/mL的蛋白酶K，室溫靜置2min，加入1/5體積的裂解液，顛倒混勻，然后在50°C保溫至少Ih，冰上冷卻； (9)轉(zhuǎn)移到滅菌離心管中，加入等體積的酚:氯仿:異戊醇，顛倒混勻，靜置5min，2600g室溫離心20min ;其中，酹:氯仿:異戍醇的體積比為25:24:1 ； (10)上清轉(zhuǎn)移到新的離心管中，加入2/5體積的5mol.^01,2/5體積的預熱至650C的CTAB/NaCl，并68°C加熱.15分鐘，加入等體積的氯仿:異戊醇，顛倒混勻，靜置5min，2600g室溫離心20min,取上清；其中，氯仿:異戍醇的體積比為24:1 ; (11)重復(10) (12)加等體積的氯仿:異戍醇,顛倒混勻，靜置5min,2600g室溫離心20min ;其中，氯仿:異戊醇的體積比為24:1 ; (13)取上清，加入0.1倍體積的3mol.T1NaAc和0.8倍體積的預冷異丙醇或是加入0.1倍體積的3mol.L-1NaAc和2倍體積的無水乙醇，-20°C放置至少20min沉淀葉綠體DNA ； (14)4°C,13500g離心20min，將沉淀用70%的乙醇室溫清洗兩次，13500g離心10min，風干，用含有RNA酶的TE重懸，37°C保溫lh，得到葉綠體DNA ； 其中，所述緩沖液A含有50mM Tris, 25mM EDTA, 1.25M NaCl,0.25M抗壞血酸，1%PVP40，0.1%BSA,1mM DTT，ρΗ3.6 ；
所述緩沖液 1 含有 1.25Μ NaCl,ρΗ8.0 的 50mM Tris_HCl，pH8.0 的 7mMEDTA, 5%PVP40,0.1%BSA,ImM DTT ； 所述洗滌緩沖液含有0.35M山梨醇，pH8.0的50mM Tris-HCl，pH8.0的25mM EDTA ；所述緩沖液 2 含有 pH8.0 的 IOmM Tris-HCl，pH8.0 的 5mM EDTA ； 所述裂解液含有20%十二烷基肌氨酸鈉，pH8.0的50mM Tris-HCl，pH8.0的25mM EDTA。
4.根據(jù)權利要求1所述的太白貝母的分子生物學鑒定方法，其特征在于，葉綠體基因組DNA測序與序列拼接是采用第二代高通量測序技術。
5.根據(jù)權利要求3所述的太白貝母的分子生物學鑒定方法，其特征在于，所述第二代高通量測序技術是通過Roche454GS FL`X高通量測序儀完成。
【文檔編號】C12Q1/68GK103849687SQ201410064966
【公開日】2014年6月11日 申請日期:2014年2月25日 優(yōu)先權日:2014年2月25日
【發(fā)明者】盧其福, 陳士林, 黃喜梅, 李西文, 劉薇, 黃林芳, 高麗, 周立 申請人:廣州白云山潘高壽藥業(yè)股份有限公司