李龍召, 黃平平, 徐奎棟, 4, 趙 峰, 4

DNA/RNA提取及處理方法對(duì)評(píng)價(jià)沉積物中纖毛蟲(chóng)分子多樣性的影響

李龍召1, 2, 黃平平3, 徐奎棟1, 2, 4, 趙 峰1, 2, 4

(1. 中國(guó)科學(xué)院海洋研究所海洋生物分類(lèi)與系統(tǒng)演化實(shí)驗(yàn)室, 山東 青島 266071; 2. 中國(guó)科學(xué)院大學(xué), 北京 100049; 3. 山東省生物物理重點(diǎn)實(shí)驗(yàn)室, 德州學(xué)院生物物理研究院, 山東 德州 253023; 4. 中國(guó)科學(xué)院海洋大科學(xué)研究中心, 山東 青島 266071)

環(huán)境DNA(eDNA)技術(shù)結(jié)合高通量測(cè)序已廣泛應(yīng)用于評(píng)價(jià)微型生物多樣性與群落構(gòu)成。相較于eDNA, RNA在環(huán)境中易降解, 環(huán)境RNA(eRNA)可更準(zhǔn)確地反映群落近期生命活性狀態(tài)。理論上, 基于eRNA檢獲的生物多樣性要低于eDNA檢獲的總體多樣性。但前期已發(fā)表研究顯示同一站點(diǎn)eDNA獲得的多樣性低于eRNA檢獲量, 推測(cè)可能原因包括: 1)樣本量不同; 2)RNA中存在DNA污染。為驗(yàn)證假設(shè), 本研究以陸架區(qū)2個(gè)站位沉積物中的纖毛蟲(chóng)為實(shí)驗(yàn)對(duì)象, 系統(tǒng)性比較4種DNA/RNA提取方法: DNA直接提取(0.9 g沉積物), 大樣本量DNA直接提取(2 g), DNA洗脫法(2 g), RNA直接提取法(2 g); 以及3種RNA提取后的處理方法: 無(wú)DNA酶反轉(zhuǎn)錄, 含DNA酶反轉(zhuǎn)錄, 先純化RNA再反轉(zhuǎn)錄。研究結(jié)果表明, 大樣本量(2 g)DNA直接提取法所檢獲的可操作分類(lèi)單元(OTUs)約為小樣本量(0.9 g)的2倍。相同樣本量下, DNA洗脫法檢獲的OTUs最多, 而DNA直接提取檢獲的OTUs低于有/無(wú)DNA酶反轉(zhuǎn)檢獲的數(shù)量, 先純化再反轉(zhuǎn)錄所獲得OTUs最少。就群落構(gòu)成而言, DNA洗脫法和RNA純化可有效降低浮游類(lèi)群比例, 可更真實(shí)展現(xiàn)底棲群落的構(gòu)成。綜上, 建議使用DNA洗脫法評(píng)價(jià)沉積物中微型生物的總體多樣性; 采用eRNA研究具有生命活性的生物群落時(shí), 反轉(zhuǎn)錄之前可純化RNA, 以獲得更加準(zhǔn)確的具有生命活性的群落信息。

eDNA; eRNA; DNA/RNA提取; 纖毛蟲(chóng); 分子多樣性

真核微生物包含原生生物, 小型后生生物和單細(xì)胞真菌等所有自養(yǎng)和異養(yǎng)類(lèi)型[1-3]。海洋底棲真核微生物是食物網(wǎng)中的重要一環(huán), 在物質(zhì)循環(huán)和能量流動(dòng)中發(fā)揮重要作用, 并維持著海洋生態(tài)系統(tǒng)的穩(wěn)定[4-6]。纖毛蟲(chóng)作為海洋底棲真核微生物中重要一員, 作用廣泛, 例如作為指示物種和遺傳發(fā)育研究的模式生物[7], 以及穩(wěn)定微食物網(wǎng)[8]等功能。研究海洋底棲纖毛蟲(chóng)多樣性有利于深入了解生物地球化學(xué)循環(huán)過(guò)程[4-6, 8], 對(duì)認(rèn)知海洋生態(tài)系統(tǒng)具有重要意義。

近年來(lái), 環(huán)境DNA(eDNA)[9-11]和高通量測(cè)序技術(shù)的快速發(fā)展[12], 為研究纖毛蟲(chóng)等微型生物的分子多樣性提供新的思路與方法。eDNA是指從環(huán)境樣本中提取的總DNA, 無(wú)需分離目標(biāo)生物[13]。eDNA在環(huán)境樣本中存在時(shí)間長(zhǎng), 構(gòu)成復(fù)雜, 來(lái)源廣泛, 包含豐富的物種遺傳信息, 既可檢獲具有活性的生物, 又可檢獲休眠狀態(tài)生物以及來(lái)源于生物分泌物或者死亡個(gè)體的胞外DNA[14-15]。目前, 許多研究針對(duì)樣本的空間異質(zhì)性和反映生物群落的時(shí)效性, 以及實(shí)驗(yàn)流程的可重復(fù)性和穩(wěn)健性展開(kāi)探討, 以期得到統(tǒng)一規(guī)范的實(shí)驗(yàn)流程[16-20]?；緦?shí)驗(yàn)流程包含: 采樣及樣品保存→DNA/RNA提取→引物設(shè)計(jì), PCR擴(kuò)增, 測(cè)序獲得樣本序列信息→數(shù)據(jù)分析, 獲得生物群落信息[21]。目前, 提取方法的選用不盡相同, 已有多份研究探討不同的DNA提取方法對(duì)分子多樣性的影響[16, 22]。整體上, 試劑盒技術(shù)成熟, 在進(jìn)行DNA提取時(shí)操作規(guī)范, 流程統(tǒng)一, 獲得的DNA純度高, 系統(tǒng)誤差較小。目前大多數(shù)的分子研究都采用試劑盒提取[16-18, 22]。在使用試劑盒提取沉積物的研究中, 樣本量和提取方法的選用爭(zhēng)議較大[16, 18]。小的樣本量(0.2～0.5 g)有利于機(jī)械化流程的開(kāi)展[16], 而大樣本量(>2 g)能更真實(shí)反映生物群落信息[23]。

相較于eDNA, RNA在環(huán)境中存在時(shí)間短[24], 作為細(xì)胞代謝的過(guò)程產(chǎn)物, 環(huán)境RNA(eRNA)能更為準(zhǔn)確地反映群落活躍狀態(tài)[25-26]。eRNA作為一種新技術(shù), 尚缺乏一套標(biāo)準(zhǔn)化的方法來(lái)進(jìn)行樣本收集、處理和分析[27]。理論上, 同一樣品, 基于eRNA檢獲的生物多樣性要低于eDNA的檢獲量。然而, 課題組前期研究以及其他報(bào)道[16, 18]發(fā)現(xiàn), 采用目前國(guó)際主流的DNA直接提取法(0.9 g沉積物)獲得的多樣性低于RNA方法(2 g沉積物)獲得的多樣性[28]。推測(cè)可能原因包括1)樣本量不同; 2)在RNA反轉(zhuǎn)錄過(guò)程中, 污染基因組DNA的存在會(huì)對(duì)下游分析產(chǎn)生較大影響, 使得分子多樣性結(jié)果偏高[29-30]。

為驗(yàn)證本研究提出的假設(shè), 以陸架區(qū)2個(gè)站位沉積物中的纖毛蟲(chóng)為實(shí)驗(yàn)對(duì)象, 系統(tǒng)性比較4種DNA/RNA提取方法: DNA直接提取(0.9 g沉積物), 大樣本量DNA直接提取(2 g), DNA洗脫法(2 g), RNA直接提取法(2 g); 以及3種RNA提取后的處理方法: 無(wú)DNA酶反轉(zhuǎn)錄, 含DNA酶反轉(zhuǎn)錄, 先純化RNA再反轉(zhuǎn)錄。旨在評(píng)價(jià): 1) 不同沉積物樣本量對(duì)基于eDNA技術(shù)檢獲的纖毛蟲(chóng)分子多樣性的影響。2) 相同樣本量, 不同DNA/RNA提取方法對(duì)分子多樣性評(píng)價(jià)的影響。3) 不同RNA處理方法對(duì)評(píng)估纖毛蟲(chóng)分子多樣性的影響。以期提出利用eDNA/eRNA研究海洋沉積物中真核微生物的分子多樣性的最佳研究策略。

1 實(shí)驗(yàn)方法

1.1 樣品采集與樣品保存

2018年8月搭乘“科學(xué)三號(hào)”科學(xué)考察船于黃海3600-3(35°39′86″N, 112°29′71″E)和3600-6站位(36°00′07″N, 123°00′02″E), 利用0.1 m2改進(jìn)型的Gray-Ohara箱式采泥器, 進(jìn)行沉積物樣品采集。在每個(gè)站位采集含有上覆水的未受擾動(dòng)的沉積物, 刮取0～2 cm表層沉積物約20 g放入封口袋中, 3600-3站位采集4個(gè)重復(fù), 3600-6站位采集3個(gè)重復(fù)。置于–20 ℃超低溫冰箱中保存, 以待后續(xù)DNA提取。

1.2 分子生物學(xué)方法

1.2.1 DNA/RNA提取

1) 取0.9 g沉積物, 采用Power Soil DNA Isolation Kit(Qiagen, Germany)進(jìn)行DNA提取, 分3次進(jìn)行, 每次0.3 g (DNA直接提取0.9 g); 2) 取2 g沉積物, 采用Power Max Soil DNA Isolation kit(Qiagen, Germany)進(jìn)行DNA提取, 每個(gè)重復(fù)提取1份DNA (DNA直接提取2 g); 3) 取2 g沉積物, 采用Power Soil RNA Isolation Kit(Qiagen, Germany)提取RNA, 同時(shí)采用RNA Power Soil?DNA Elution Accessory Kit(Qiagen, Germany)于同一樣品洗脫DNA, 每個(gè)重復(fù)提取1份RNA和1份DNA(DNA洗脫法2 g)。

1.2.2 RNA反轉(zhuǎn)錄

對(duì)提取的RNA采用3種不同的方式反轉(zhuǎn)錄處理獲取cDNA: 1) 采用不含DNA酶的反轉(zhuǎn)錄試劑盒Prime Script II 1st strand cDNA Synthesis Kit(TaKaRa, Japan)直接反轉(zhuǎn)錄(無(wú)DNA酶反轉(zhuǎn)); 2) 采用含DNA酶的反轉(zhuǎn)錄試劑盒5X All-In-One Master Mix (Gentaur, Europe)進(jìn)行反轉(zhuǎn)錄(含DNA酶反轉(zhuǎn)); 3) 先采用TURBO DNA-free? Kit(Thermo Fisher, USA)對(duì)提取的RNA進(jìn)行純化, 純化后的RNA采用Prime Script II 1st strand cDNA Synthesis Kit反轉(zhuǎn)錄(先純化再反轉(zhuǎn))(圖1)。

1.2.3 PCR與高通量測(cè)序

采用巢式PCR技術(shù)擴(kuò)增纖毛蟲(chóng)的18S rRNA基因的V4高變區(qū), 選用纖毛蟲(chóng)特異性引物CilF、Cil RⅠ、Cil RⅡ、Cil RⅢ[31], 以及真核微生物18S V4 高變區(qū)通用引物EukF、EukR[32]進(jìn)行擴(kuò)增, 每個(gè)樣品設(shè)置3份PCR重復(fù)及1份陰性空白對(duì)照, 具體巢式PCR流程參考文獻(xiàn)[33]。擴(kuò)增后使用瓊脂糖凝膠電泳檢測(cè)PCR產(chǎn)物質(zhì)量及片段長(zhǎng)度(480 bp), 將同一個(gè)樣品的3份PCR產(chǎn)物合并, 送至測(cè)序公司開(kāi)展測(cè)序。

1.3 生物信息分析

測(cè)序數(shù)據(jù)進(jìn)行拼接得到原始數(shù)據(jù)(raw data), 經(jīng)過(guò)質(zhì)量控制和去嵌合體得到可用于后續(xù)分析的有效數(shù)據(jù)(clean data)。質(zhì)量控制包括: 1) 將Raw Tags從連續(xù)低質(zhì)量值(≤19)堿基數(shù)達(dá)到設(shè)定長(zhǎng)度(3)的第一個(gè)低質(zhì)量堿基位點(diǎn)截?cái)? 2) 進(jìn)一步過(guò)濾掉其中連續(xù)高質(zhì)量堿基長(zhǎng)度小于Tags長(zhǎng)度75%的Tags。對(duì)有效數(shù)據(jù)進(jìn)行去冗余、去單一序列和 97%相似度的OUT (operational taxonomic units, 可操作分類(lèi)單元)聚類(lèi)分析, 根據(jù) OTU 聚類(lèi)結(jié)果, 對(duì)每個(gè) OTU的代表序列與 SILVA 數(shù)據(jù)庫(kù)比對(duì)作物種注釋, 得到對(duì)應(yīng)的物種信息和基于物種的豐度分布情況, 即OTU表格。利用EXCEL進(jìn)行統(tǒng)計(jì), 采用R進(jìn)行單因素方差分析, 并繪制箱線(xiàn)圖, 折線(xiàn)圖, 韋恩圖。使用PRIMER v6軟件中的CLUSTER應(yīng)用程序, 對(duì)同一樣本量(2 g)的所有樣本進(jìn)行聚類(lèi)分析。

2 實(shí)驗(yàn)結(jié)果

2.1 提取和處理方法對(duì)纖毛蟲(chóng)物種豐富度的影響

基于不同樣品量的DNA直接提取法比較結(jié)果顯示, 0.9 g樣本量所獲得的OTUs數(shù)明顯低于2 g樣本量所獲得的OTUs數(shù)。在3600-3站位, 直接提取DNA: 0.9 g樣本平均獲得96個(gè)OTUs, 2 g樣本量每個(gè)測(cè)序樣本平均獲得152個(gè)OTUs。在3600-6站位, 直接提取DNA: 0.9 g樣本量平均獲得83個(gè)OTUs, 2 g樣本量每個(gè)測(cè)序樣本平均獲得161個(gè)OTUs。

在相同樣本量(2 g)的情況下, 以3600-3站位為研究對(duì)象, 結(jié)果表明: 采用DNA洗脫法平均檢獲的OTUs最高; 其他方法排序: 依次為無(wú)DNA酶反轉(zhuǎn), 含DNA酶反轉(zhuǎn), DNA直接提取以及RNA先純化再反轉(zhuǎn)的方法。以3600-6站位為研究對(duì)象, DNA洗脫法平均所檢獲的OTUs亦最高; 其他依次為含DNA酶反轉(zhuǎn)、DNA直接提取、無(wú)DNA酶的反轉(zhuǎn)以及RNA先純化再反轉(zhuǎn)的方法(圖2)。在3600-3站位, 不同提取方法下的單因素方差分析(方差齊性檢驗(yàn)= 0.3, 符合正態(tài)分布)顯示, 不同提取方法之間存在顯著差異(= 0.023); 在3600-6站位, 不同提取方法下的單因素方差分析(方差齊性檢驗(yàn)= 0.1, 符合正態(tài)分布)顯示, 不同提取方法之間存在極顯著差異(= 0.007 5)。

圖2 3600-3與3600-6站位同一樣本量不同提取處理方法所獲得的OTU數(shù)量

2.2 提取處理方法與纖毛蟲(chóng)生物群落結(jié)構(gòu)的關(guān)系

總體來(lái)看, 直接DNA提取的方法獲得大量的浮游類(lèi)纖毛蟲(chóng)(Choreotrichia 以及Oligotrichia)序列。采用樣本量0.9 g, 在3600-3和3600-6站位平均獲得的浮游類(lèi)纖毛蟲(chóng)相對(duì)豐度最高(52%±0.155%; 51%±0.079%); 采用樣本量2 g, 平均獲得的浮游類(lèi)纖毛蟲(chóng)相對(duì)豐度次之(41%±0.145%; 38%±0.157%); 其余的方法(DNA洗脫法與RNA提取方法)所獲得的浮游類(lèi)群相對(duì)豐度較低。在3600-3站位, 無(wú)DNA酶反轉(zhuǎn)錄方法獲得的浮游類(lèi)纖毛蟲(chóng)較高, 平均相對(duì)豐度為20.9%±0.116%, 其余依次為洗脫DNA方法(20.2%± 0.156%)、含DNA酶的反轉(zhuǎn)錄方法(13.9%±0.028%)以及RNA先純化再反轉(zhuǎn)方法(11.5%±0.076%)。在3600-6站位, 無(wú)DNA酶反轉(zhuǎn)錄方法獲得的浮游類(lèi)纖毛蟲(chóng)亦較高, 平均相對(duì)豐度為13.4%±0.073%, 其余依次為含DNA酶的反轉(zhuǎn)錄方法(7.9%±0.054%)、DNA洗脫方法(7.0%±0.015%)以及RNA先純化再反轉(zhuǎn)方法(1.5%±0.002%) (圖3)。

圖3 每個(gè)測(cè)序樣本浮游類(lèi)群(舞毛亞綱和寡毛亞綱)相對(duì)豐度折線(xiàn)圖

在樣本量相同的情況下, 綜合所有重復(fù), 分析不同方法獲得的纖毛蟲(chóng)總體上以及刪除浮游類(lèi)群后所共有和特有的OTUs。結(jié)果表明, DNA洗脫法檢獲的特有OTUs數(shù)最多, 3600-3站位總體上51個(gè), 刪除浮游類(lèi)群后47個(gè); 3600-6站位總體上71個(gè), 刪除浮游類(lèi)群后57個(gè)。而DNA提取法在兩個(gè)站位檢獲的特異OTUs數(shù)波動(dòng)較大(總體上, 3600-3: 6個(gè); 3600-6: 33個(gè))。有或無(wú)DNA酶反轉(zhuǎn)所檢獲的特異OTUs數(shù)亦較高, 先純化再反轉(zhuǎn)所檢獲的特異OTUs數(shù)最低。以不同的RNA反轉(zhuǎn)錄方法為研究對(duì)象, 總體上看, 無(wú)論是在3600-3, 還是在3600-6站位; 含DNA酶反轉(zhuǎn)錄所獲得的分子多樣性和無(wú)DNA酶反轉(zhuǎn)錄所共檢的OTUs比例較高, 均在50%以上。在3600-3站位, 總體上含DNA酶和無(wú)DNA酶共檢獲OTUs占比67.8%; 去除浮游類(lèi)群后為65.6%。在3600-6站位, 總體上含DNA酶和無(wú)DNA酶共檢獲OTUs占比52%。去除浮游類(lèi)群后發(fā)現(xiàn)為52.7%。綜合兩個(gè)站位, 對(duì)比含DNA酶反轉(zhuǎn)和先純化再反轉(zhuǎn)的方法發(fā)現(xiàn), 無(wú)論是否去除浮游類(lèi)群, 先純化再反轉(zhuǎn)所檢獲的特異OTUs占比較少, 均在10%以下, 分別為6.9%, 7.8%, 2.9%, 3.3%(圖4)。

分別對(duì)3600-3和3600-6站位相同樣本量(2 g)的測(cè)序樣本聚類(lèi)分析, 結(jié)果顯示: 在3600-3站位, 無(wú)DNA酶反轉(zhuǎn)以及含DNA酶反轉(zhuǎn)獲得的纖毛蟲(chóng)群落結(jié)構(gòu)較為相似, 然后與DNA洗脫法獲得的樣本群落聚合; 在3600-6站位, 含DNA酶反轉(zhuǎn)獲得的樣本群落首先與DNA洗脫法獲得的樣本群落聚合, 其次與無(wú)DNA酶反轉(zhuǎn)獲得的樣本群落聚合, 最后與DNA直接提取獲得的樣本群落聚合。結(jié)果顯示, 無(wú)論是在哪個(gè)站位, 先純化再反轉(zhuǎn)所獲得的樣本群落信息與其他方法均在相似度50%以下聚合(圖5)。

3 討論

3.1 樣本量對(duì)纖毛蟲(chóng)分子多樣性的影響

研究表明, 不同樣本量對(duì)重建土壤微生物群落(細(xì)菌, 真菌等)信息存在影響[34-36]。Nascimento等[18]利用不同的巖芯采樣量評(píng)估了樣本量對(duì)于真核生物多樣性及群落結(jié)構(gòu)的影響, 分析得出隨著樣本量的增加, 群落結(jié)構(gòu)信息隨之穩(wěn)定; 而在Pearman等[16]的研究中發(fā)現(xiàn)在低重復(fù)樣品數(shù)的情況下, 不同樣本量檢獲的群落差異不顯著。為降低沉積物樣品異質(zhì)性的影響, 我們的研究將每種處理方式設(shè)置3～4個(gè)重復(fù), 結(jié)果顯示不同樣本量對(duì)纖毛蟲(chóng)分子多樣性的評(píng)價(jià)具有重要影響。采用同種提取方法, 0.9 g樣本量平均所獲得的OTUs數(shù)約為2 g樣本量所獲得的OTUs數(shù)的1/2。從側(cè)面證明樣品量小是導(dǎo)致從0.9 g沉積物直接提取DNA獲得的多樣性低于從2 g沉積物提取RNA獲得的多樣性的重要因素。

圖4 3600-3與3600-6站位總體上和刪除浮游類(lèi)群(舞毛亞綱以及寡毛亞綱)后, 不同RNA處理方法之間共有和特有OTUs

生物群落構(gòu)成分析結(jié)果顯示0.9 g樣本量所檢測(cè)到的浮游類(lèi)纖毛蟲(chóng)占比最多, 說(shuō)明小樣本量的DNA提取方法檢測(cè)能力有限, 富集于海底的浮游物種的DNA會(huì)干擾底棲群落的研究。提取DNA時(shí), 由于小樣本量抽樣不全[18], 會(huì)導(dǎo)致多樣性的漏檢, 進(jìn)而形成生物群落構(gòu)成的偏好性等。

3.2 不同提取方法對(duì)纖毛蟲(chóng)分子多樣性的影響

利用eDNA作為物種檢測(cè)手段已經(jīng)十分成熟, 但目前仍沒(méi)有一套統(tǒng)一規(guī)范的實(shí)驗(yàn)流程。現(xiàn)有的研究結(jié)果表明, 從環(huán)境樣品中提取DNA, 不同的提取處理方法會(huì)導(dǎo)致分析結(jié)果的差異[16, 22]。本研究發(fā)現(xiàn)相較DNA提取法, DNA洗脫法可檢獲更多的OTUs, 這可能歸因于直接提取DNA時(shí), 腐殖質(zhì)含量高, 更易與DNA的結(jié)合, 進(jìn)而使得獲得的生物遺傳信息降低[37]。DNA洗脫法可更高效除去腐殖質(zhì)的影響, 獲得更為全面的分子多樣性信息, 可為 eDNA研究提供可靠的模板[38]。

相較于eDNA, RNA在環(huán)境中存在時(shí)間短[24], 作為細(xì)胞代謝的過(guò)程產(chǎn)物, 環(huán)境RNA(eRNA)能更準(zhǔn)確地反映群落近期的生命活性狀態(tài)[25-26]。相較于有/無(wú)DNA酶反轉(zhuǎn)錄, 先純化再反轉(zhuǎn)錄所獲得的OTUs數(shù)量明顯偏低; 有DNA酶反轉(zhuǎn)錄所獲得的OTUs比無(wú)DNA酶反轉(zhuǎn)錄所獲得的少, 推測(cè)可能是樣本中殘留的基因組DNA所致。這表明在RNA反轉(zhuǎn)錄過(guò)程中去除基因組污染是必要的[39]。

分析不同提取和處理方法所得群落信息, 發(fā)現(xiàn)采用RNA反轉(zhuǎn)錄方法所獲得的浮游生物占比較少, 且純化后RNA獲得浮游生物占比最少。先純化再反轉(zhuǎn)錄所獲得的群落單獨(dú)聚合, 并與其他方法之間存在差異(相似度小于50%)?？紤]到RNA作為一種易降解的活性遺傳物質(zhì)[40], 我們推測(cè)先純化再反轉(zhuǎn)的方法可去除環(huán)境樣本中較長(zhǎng)時(shí)間內(nèi)的基因組DNA等污染, 因此導(dǎo)致RNA先純化再反轉(zhuǎn)獲得的OTUs也較少, 且群落構(gòu)成與其他方法存在差異, 這也從側(cè)面證實(shí)RNA先純化再反轉(zhuǎn)可更加準(zhǔn)確的展現(xiàn)具有生命活性的生物群落。

圖5 3600-3和3600-6站位相同樣本量(2 g)測(cè)序樣本聚類(lèi)分析圖

4 結(jié)論

增加樣本量(0.9 g至2 g范圍內(nèi))可檢獲更多的OTUs。相較于DNA提取法, DNA洗脫法可更加全面的反映微型生物的群落結(jié)構(gòu)。RNA在環(huán)境樣本中存在時(shí)間短, 可反映具有生命活性的微生物群落信息情況, 但直接提取的eRNA中存在DNA污染, 提取后直接反轉(zhuǎn)錄, 會(huì)檢獲較高比例的浮游類(lèi)群。RNA提取純化后再反轉(zhuǎn)錄可有效降低浮游類(lèi)群比例, 可更加真實(shí)展現(xiàn)具有生命活性的底棲生物群落。綜上, 我們推薦采用DNA洗脫法評(píng)價(jià)沉積物中微型生物的整體多樣性; 采用eRNA研究具有生命活性的生物群落時(shí), 反轉(zhuǎn)錄之前可先純化RNA, 以獲得更加準(zhǔn)確的具有生命活性的底棲生物群落信息。

致謝: 中國(guó)科學(xué)院海洋研究所徐雨協(xié)助樣品采集, 謹(jǐn)致謝忱。

[1] MCGRATH C L, KATZ L A. Genome diversity in microbial eukaryotes[J]. Trends in Ecology & Evolution, 2004, 19(1): 32-38.

[2] BIK H M, SUNG W, DELEY P, et al. Metagenetic community analysis of microbial eukaryotes illuminates biogeographic patterns in deep-sea and shallow water sediments[J]. Molecular Ecology Resources, 2012, 21(5): 1048-1059.

[3] 徐奎棟. 海洋微型底棲生物的多樣性與地理分布[J]. 生物多樣性, 2011, 19(6): 661-675.

XU Kuidong. Biodiversity and biogeography of marine microbenthos: progress and prospect[J]. Biodiversity Science, 2011, 19(6): 661-675.

[4] SHERR E B, SHERR B F. Significance of predation by protists in aquatic microbial food webs[J]. Antonie Van Leeuwenhoek International Journal of General and Molecular Microbiology, 2002, 81(1/4): 293-308.

[5] ARISTEGUI J, GAASOL J M, DUARTE C M, et al. Microbial oceanography of the dark ocean’s pelagic realm[J]. Limnology and Oceanography, 2009, 54(5): 1501-1529.

[6] SNELGROVE P, BLACKBURN T H, HUTCHINGS P A, et al. The importance of marine sediment biodiversity in ecosystem processes[J]. Ambio, 1997, 26(8): 578-583.

[7] PUCCIARELLI S, BALLARININ P, BARCHETTA S, et al. Ciliates as models to study diversities in microtubule functions by functional genomics[J]. New Biotechnology, 2010, 27: S32.

[8] 陳伊凡. 海洋纖毛蟲(chóng)攝食微型藻類(lèi)過(guò)程中磷的動(dòng)力學(xué)研究[D]. 廈門(mén): 廈門(mén)大學(xué), 2014.

Chen Yifan. Kinetics of phosphorus uptake by marine ciliates on microalgae[D]. Xiamen: Xiamen University, 2014.

[9] JERDE C L, MAHON A R, CHADDERTON W L, et al. “Sight-unseen” detection of rare aquatic species using environmental DNA[J]. Conservation Letters, 2011, 4(2): 150-157.

[10] DEJEAN T, VALENTINI A, DUPARC A, et al. Persistence of Environmental DNA in Freshwater Ecosystems[J]. PloS One, 2011, 6(8): e23398.

[11] CORINALDESI C, BARUCCA M, LUNA G M, et al. Preservation, origin and genetic imprint of extracellular DNA in permanently anoxic deep-sea sediments[J]. Molecular Ecology, 2011, 20(3): 642-654.

[12] LOMAN N J, MISRA R V, DALLMAN T J, et al. Performance comparison of benchtop high-throughput sequencing platforms[J]. Nature Biotechnology, 2012, 30(5): 434-439.

[13] TABERLET P, COISSAC E, HAJIBABAEI M, et al. Environmental DNA[J]. Molecular Ecology, 2012, 21(8): 1789-1793.

[14] HARRISON J B, SUNDAY J M, ROGERS S M. Predicting the fate of eDNA in the environment and implications for studying biodiversity[J]. Proceedings of the Royal Society B-Biological Sciences, 2019, 286(1915): 20191409.

[15] PIETRAMELLARA G, ASCHER J, BORGOGNI F, et al. Extracellular DNA in soil and sediment: fate and ecological relevance[J]. Biology and Fertility of Soils, 2009, 45(3): 219-235.

[16] PEARMAN J K, KEELEY N B, WOOD S A, et al. Comparing sediment DNA extraction methods for assessing organic enrichment associated with marine aquaculture[J]. PeerJ, 2020, 8: e10231.

[17] GUERRIERI A, BONIN A, MUNKEMULLER T, et al. Effects of soil preservation for biodiversity monitoring using environmental DNA[J]. Molecular Ecology Resources, 2020, 30(13): 3313-3325.

[18] NASCIMENTO F J A, LALLIAS D, BIK H M, et al. Sample size effects on the assessment of eukaryotic diversity and community structure in aquatic sediments using high-throughput sequencing[J]. Scientific Reports, 2018, 8(1): 11737.

[19] CALDERONSANOU I, MUNKEMULLER T, BOYER F, et al. From environmental DNA sequences to ecological conclusions: How strong is the influence of methodological choices?[J]. Journal of Biogeography, 2019, 47(1): 193-206.

[20] PEREZBROCAL V, MAGNE F, RUIZ S, et al. Optimized DNA extraction and purification method for characterization of bacterial and fungal communities in lung tissue samples[J]. Scientific Reports, 2020, 10(1): 17377.

[21] ZINGER L, BONIN A, ALSOS I G, et al. DNA metabarcoding-Need for robust experimental designs to draw sound ecological conclusions[J]. Molecular Eco-logy, 2019, 28(8): 1857-1862.

[22] CHANGEY F, BLAUD A, PANDO A, et al. Monitoring soil microbial communities using molecular tools: DNA extraction methods may offset long-term management effects[J]. European Journal of Soil Science, 2020, 31(1): 134-148.

[23] LEAR G, DICKIE I, BANKS J, et al. Methods for the extraction, storage, amplification and sequencing of DNA from environmental samples[J]. New Zealand Journal of Ecology, 2018, 42(1): 10, 1A-50A.

[24] WOOD S A, BIESSY L, LATCHFORD J L, et al. Release and degradation of environmental DNA and RNA in a marine system[J]. Science of the Total Environment, 2020, 704: 135314.

[25] GUARDIOLA M, WANGENSTEEN O S, TABERLET P, et al. Spatiotemporal monitoring of deep-sea communities using metabarcoding of sediment DNA and RNA[J]. PeerJ, 2016, 4: e2807.

[26] PAWLOWSKI J, ESLING P, LEJZEROWICZ F, et al. Environmental monitoring through protist next-generation sequencing metabarcoding: assessing the impact of fish farming on benthic foraminifera communities[J]. Molecular Ecology Resources, 2014, 14(6): 1129-1140.

[27] WOOD S A, POCHON X, MING W, et al. Considerations for incorporating real-time PCR assays into routine marine biosecurity surveillance programmes: a case study targeting the Mediterranean fanworm (Sabella spallanzanii) and club tunicate (Styela clava)[J]. Genome, 2019, 62(3): 137-146.

[28] HUANG P, ZHAO F, XU K. Complementary DNA sequencing (cDNA): an effective approach for assessing the diversity and distribution of marine benthic ciliates along hydrographic gradients[J]. Journal of Oceanology and Limnology, 2020, 39(1): 208-222.

[29] LARMAN M G, KATZJAFFE M G, MCALLIE B, et al. Analysis of global gene expression following mouse blastocyst cryopreservation[J]. Human Reproduction, 2011, 26(10): 2672-2680.

[30] ALKHALIL A, HAMMAMIEH R, HARDICK J, et al. Gene expression profiling of monkeypox virus-infected cells reveals novel interfaces for host-virus interactions[J]. Virology Journal, 2010, 7: 173. DOI: 10.1186/ 1743-422X-7-173.

[31] LARA E, BERNEY C, HARMS H, et al. Cultivation-independent analysis reveals a shift in ciliate 18S rRNA gene diversity in a polycyclic aromatic hydrocarbon-polluted soil[J]. FEMS Microbiology Ecology, 2007, 62(3): 365-73.

[32] STOECK T, BASS D, NEBEL M, et al. Multiple marker parallel tag environmental DNA sequencing reveals a highly complex eukaryotic community in marine anoxic water[J]. Molecular Ecology, 2010, 19: 21-31.

[33] STOCK A, EDGCOMB V, ORSI W, et al. Evidence for isolated evolution of deep-sea ciliate communities through geological separation and environmental selection[J]. BMC Microbiology, 2013, 13: 150. DOI: 10.1186/1471-2180-13-150.

[34] KANG S, MILLS A L. The effect of sample size in studies of soil microbial community structure[J]. Journal of Microbiological Methods, 2006, 66(2): 242-250.

[35] ELLINGSOE P, JOHNSEN K. Influence of soil sample sizes on the assessment of bacterial community structure[J]. Soil Biology & Biochemistry, 2002, 34(11): 1701-1707.

[36] PENTON C R, GUPTA V S R, YU J, et al. Size matters: assessing optimum soil sample size for fungal and bacterial community structure analyses using high throu-ghput sequencing of rRNA gene amplicons[J]. Frontiers in Microbiology, 2016, 7: 824. DOI: 10.3389/fmicb. 2016.00824.

[37] LAKAY F M, BOTHA A, PRIOR B A. Comparative analysis of environmental DNA extraction and purification methods from different humic acid-rich soils[J]. Journal of Applied Microbiology, 2007, 102(1): 265- 273.

[38] LEJZEROWICZ F, ESLING P, PILLET L, et al. High- throughput sequencing and morphology perform equally well for benthic monitoring of marine ecosystems[J]. Scientific Reports, 2015, 5: 13932. DOI: https://doi.org/ 10.1038/srep13932.

[39] NORHAZLIN J, NORASHIKIN M N, HOH B P, et al. Effect of DNase treatment on RNA extraction from preimplantation murine embryos[J]. Genetics and Molecular Research, 2015, 14(3): 10172-10184.

[40] LI Y F, BREAKER R R. Kinetics of RNA degradation by specific base catalysis of transesterification involving the 2'-hydroxyl group[J]. Journal of the American Chemical Society, 1999, 121(23): 5364-5372.

Effects of DNA and RNA extraction methods for the evaluation of ciliate diversity in marine sediments

LI Long-zhao1, 2, HUANG Ping-ping3, XU Kui-dong1, 2, 4, ZHAO Feng1, 2, 4

(1. Department of Marine Organism Taxonomy and Phylogeny Institute of Oceanology, Chinese Academy of Sciences, Qingdao 266071, China; 2. University of Chinese Academy of Sciences, Beijing 100049, China; 3. Shandong Key Laboratory of Biophysics, Institute of Biophysics, Dezhou University, Dezhou 253023, China; 4. Center for Ocean Mega-Science, Chinese Academy of Sciences, Qingdao 266071, China)

Environmental DNA (eDNA) has been used along with high-throughput sequencing to evaluate microbial diversity and community composition. Compared with eDNA, environmental RNA (eRNA) degrades easily and can be used to investigate the community of active taxon. The species richness obtained by eRNA is lower than that by eDNA in theory. However, previous studies have shown that the diversity obtained by eDNA at the same site is lower than that detected by eRNA. This unexpected finding may be attributed to the difference in the sample size and the DNA contamination in RNA. To test this hypothesis, we estimated the effects of DNA and RNA extraction methods on ciliate diversity in sediments obtained from the continental shelf area. Four extraction methods were performed: DNA direct extraction with 0.9 g sediments, DNA direct extraction with 2 g sediments, DNA elution, and RNA direct extraction from 2 g sediments. Meanwhile, three types of eRNA treatments before the reverse transcription, including the addition of DNA enzyme, no addition of DNA enzyme, and purifying RNA, were also compared. The results revealed that the number of operational taxonomic units (OTUs) detected from 2 g sediments was about twice that from the 0.9 g sediments by DNA direct extraction. Using the same sample size, the highest number of OTUs was detected via DNA elution and the lowest number by purified RNA. The number of OTUs detected by DNA direct extraction was lower than that by eRNA, irrespective of the addition of DNA enzyme. In terms of community composition, DNA elution and RNA purification can effectively reduce the proportion of planktonic groups, thereby improving the estimates of the composition of the benthic community. In summary, the DNA elution method is recommended for evaluating the total diversity of microorganisms in marine sediments. eRNA can be used after purification to investigate active communities.

eDNA: eRNA; DNA / RNA extraction; ciliate; molecular diversity

Sep. 22, 2021

[National Natural Science Foundation of China, No. 41876171]

Q958

A

1000-3096(2022)07-0052-09

10.11759/hykx20210922001

2021-09-22;

2021-10-28

國(guó)家自然科學(xué)基金資助項(xiàng)目(41876171)

李龍召(1997—), 男, 安徽定遠(yuǎn)人, 碩士研究生, 主要從事海洋真核微生物多樣性研究, E-mail: 15375544229@163.com; 趙峰(1987—),通信作者, 副研究員, E-mail: fzhao@qdio.ac.cn

(本文編輯: 趙衛(wèi)紅)