郭娟 曾燏 鄭羽晨 趙雨雙 鐘清

摘? ? 要：為探究實驗室養(yǎng)殖條件下麥穗魚（Pseudorasbora parva）eDNA濃度與其生物量間關(guān)系，根據(jù)麥穗魚線粒體COⅠ基因序列設(shè)計特異性引物，通過實時熒光定量PCR（qPCR），得到qPCR標準曲線，采用線性、指數(shù)和對數(shù)3種模型對麥穗魚生物量與eDNA濃度的相關(guān)性曲線進行擬合和比較分析。結(jié)果表明：基于CO I基因設(shè)計的引物對麥穗魚DNA具有特異性擴增，擴增效率為90.94%；qPCR的閾值循環(huán)數(shù)（Cycle threshold，Ct值）與已知濃度或拷貝數(shù)的樣品（標準品）間呈現(xiàn)良好的線性關(guān)系（R2=0.995 2），說明標準品的拷貝數(shù)與Ct值之間相關(guān)的可信度較好。通過3種曲線進一步擬合比較分析顯示：線性模型為麥穗魚生物量與eDNA濃度間的最佳擬合模型，其R2高達0.933 8，y=72.087x+81.978。這表明在本試驗條件下，eDNA濃度與麥穗魚生物量呈正相關(guān)性。

關(guān)鍵詞：eDNA濃度；生物量；麥穗魚

中圖分類號：S932? ? ? ? ? 文獻標識碼：A? ? ? ? ? DOI 編碼：10.3969/j.issn.1006－6500.2023.06.008

Abstract：In order to investigate the relationship between the environmental DNA （eDNA） concentration and biomass of Pseudorasbora parva reared in the laboratory， specific primers were designed based on the mitochondrial COI gene sequence of Pseudorasbora parva. The qPCR standard curve was obtained through quantitative real-time PCR （qPCR）. The correlation curves between Pseudorasbora parva biomass and eDNA concentration were fitted and compared using linear， exponential， and logarithmic models. The results showed that the primers designed based on the COI gene specifically amplified Pseudorasbora parva DNA， with an amplification efficiency of 90.94%. The cycle threshold （Ct） value of qPCR showed a good linear relationship with known concentration or copy number of samples （standards） （R2=0.995 2）， indicating a good credibility of the correlation between the copy number of standards and the Ct value. Further fitting and comparison of the three curves revealed that the linear model was the best fit for the relationship between Pseudorasbora parva biomass and eDNA concentration， with a high R2 of 0.933 8 and y=72.087x+81.978. This indicates a positive correlation between eDNA concentration and Pseudorasbora parva biomass under the experimental conditions.

Key words： eDNA concentration; biomass; Pseudorasbora parva

監(jiān)測魚類生物量是入侵物種種群管理工作的重要組成部分，然而，直接使用傳統(tǒng)捕撈調(diào)查方法會破壞水生生態(tài)環(huán)境，也難以準確評估個體小和群眾數(shù)量少的物種生物量[1]。近年來，隨著分子生物學的發(fā)展，eDNA技術(shù)作為一種新的監(jiān)測方法，憑借便捷、環(huán)境友好、經(jīng)濟高效等優(yōu)點在珍稀瀕危物種和入侵物種監(jiān)測領(lǐng)域發(fā)揮了重要作用[2-3]。環(huán)境DNA（Environmental DNA，eDNA）是指從環(huán)境樣本中（水、土壤、空氣等）提取到的總DNA[4]。目前，eDNA主要運用于4個方面。一是檢測目標物種的有無，通過檢測水體中是否有目標物種釋放的DNA來判斷水體中是否有該物種的存在[5]；二是生物多樣性評估，高通量測序的出現(xiàn)使得從一個eDNA樣品中同時檢測多種生物的DNA成為可能，現(xiàn)階段已被用于水生生態(tài)系統(tǒng)生物多樣性的檢測[6-7]；三是魚類分布，利用eDNA方法可以識別出一個特定物種在一個大區(qū)域（如海灣和大湖泊）上的分布，明確目標物種分布將有助于管理和監(jiān)測海洋和大型湖泊中的魚類[8]；四是生物量估算，對目標生物的生物量進行調(diào)查是水生生物資源監(jiān)測的重要內(nèi)容[9-10]。

近年來我國開展了大量與eDNA技術(shù)相關(guān)的研究，但主要運用的是檢測目標物種、生物多樣性評估和魚類分布3個方面，關(guān)于eDNA技術(shù)評估生物量的研究較少[1]。魚類資源的評估不僅要定性地評估魚類多樣性，也要定量地評估生物量，實現(xiàn)這些將有助于eDNA技術(shù)成為管理魚類資源的一般工具[11]。生物量是指某一特定時刻單位空間的個體數(shù)、質(zhì)量，如何科學地開展生物量評估是魚類種群生態(tài)的一個重難點問題，準確的生物量評估有助于魚類資源的合理開發(fā)利用和外來入侵物種的有效管理，建立eDNA濃度與生物量間關(guān)系是基于eDNA技術(shù)探究野生魚類生物量的重要依據(jù)[12]。

麥穗魚（Pseudorasbora parva）隸屬于鯉形目（Cypriniformes）、鯉科（Cyprinidae）、鮈亞科（Gobioninae）、麥穗魚屬（Pseudorasbora）。麥穗魚主要以浮游動物為食，分布廣，世代周期短，繁殖能力和適應(yīng)能力強，已成功入侵全球超過30個國家[13-14]。當前通過傳統(tǒng)種群生態(tài)學方法難以準確的評估麥穗魚種群的生物量，eDNA方法的大量應(yīng)用使得這一難題的順利解決成為了可能。目前關(guān)于麥穗魚的研究主要包括麥穗魚養(yǎng)殖、生活史和疾病等方面，而基于eDNA評估麥穗魚生物量未見報道[15-17]。

因此，本試驗以實驗室人工養(yǎng)殖的麥穗魚為目標物種，通過設(shè)計麥穗魚的特異性引物，實時熒光定量PCR建立標準曲線和熔解曲線，采用線性、指數(shù)和對數(shù)3種模型對麥穗魚生物量與eDNA濃度的相關(guān)性曲線進行了擬合和分析。研究結(jié)果有助于建立實驗室條件下eDNA濃度與其生物量的最佳相關(guān)性曲線，為eDNA濃度評估野外麥穗魚生物量提供理論支撐。

1 材料與方法

1.1 特異性引物設(shè)計

針對麥穗魚的線粒體CO I基因設(shè)計特異性引物。在GenBank上檢索麥穗魚的基因序列，根據(jù)已發(fā)表的麥穗魚序列（基因序列號MT805712.1），利用相關(guān)軟件Primer-Blast設(shè)計引物，使用引物分別對麥穗魚DNA和12種鮈亞科魚類DNA進行擴增，PCR試驗程序為：95 ℃預(yù)變性2 min；95 ℃變性20 s，55 ℃退火20 s，72 ℃延伸10 s，40個循環(huán)；72 ℃延伸7 min。擴增后對PCR產(chǎn)物凝膠電泳檢測，DNA的膠回收、連接和測序，最后比對測序結(jié)果是否是目的基因。

1.2 eDNA采集和提取

在5個魚缸（50 cm×45 cm×40 cm）中加入50 L曝氣水，溫度為5 ℃。分別在5個缸中放置0、2、4、8、16尾麥穗魚，平均每尾質(zhì)量4.5 g。試驗期為7 d，7 d內(nèi)不喂食，不換水。在試驗第7天采集樣本，使用消毒殺菌后的EP管取養(yǎng)殖缸表層水50 mL，每個缸做3次重復(fù)，設(shè)置1個蒸餾水陰性對照，共計18個樣本。采樣后立即用直徑為50 mm，孔徑為0.45 μm的混合纖維素濾膜進行無菌真空抽濾裝置抽濾，每次抽濾之前用無菌水沖洗過濾漏斗，防止交叉污染。抽濾完成后，用無菌鑷子將每張濾膜單獨卷起來，放入5 mL 的離心管中并做好標記，按照DNeasy PowerWater kit（Qiagen，Germany）試劑盒的操作方法提取eDNA，將收集到的DNA放置-20 ℃保存。

1.3 實時熒光定量PCR

使用TOYOBO的SYBR Green Real-time PCR Master Mix 試劑盒和 Applied Biosystems StepOneTM qPCR儀。將10 μL的2× SYBR real-time PCR premixture和0.4 μL的正向引物與0.4 μL反向引物（10 μmol·L-1）混合為mixture A體系，上機體系為：8 μL的mixture A和8 μL模板稀釋樣。試驗程序為：95 ℃預(yù)變性2 min；95 ℃變性20 s，55 ℃退火20 s，72 ℃延伸10 s，40個循環(huán)；72 ℃延伸7 min，在每次運行退火之后采集熒光信號。首先，將目標片段用Taq酶擴增，連接T載體，裝入感受態(tài)細菌，挑單克隆，搖菌，經(jīng)測序確認片段正確后，在潔凈工作臺進行質(zhì)粒小抽和定量。將標準品依次稀釋，形成 101～107 copies·μL-1的7個梯度，制作熒光定量PCR標準曲線。每次試驗均設(shè)置陰性對照品。為了確定PCR反應(yīng)效率和定量未知樣本，計算并繪制熔解曲線和標準曲線。

1.4 擬合曲線模型的選擇

選用線性、指數(shù)和對數(shù)3種曲線模型進行擬合，利用R軟件繪制曲線圖，通過比較不同模型曲線擬合度R2的值大小來確定最佳擬合曲線。

2 結(jié)果與分析

2.1 提取DNA分析

從水體樣本提取的 eDNA和魚類肌肉中提取的 DNA 濃度分別為20～89? ng·μL-1和376～719 ng·μL-1。提取 eDNA 的OD值（A260/A280）為1.60～2.04，肌肉 DNA 的OD值為 1.79～2.13。盡管有部分eDNA濃度較低，但總體DNA的OD值能滿足后續(xù)試驗需求。陰性對照組和空白組均無擴增。

2.2 麥穗魚引物和標準曲線

得到麥穗魚引物F：CCGATTCTCATTGGGGGCTT，

R：CCTGCGAGAGGGGGATAAAC，片段長度188 bp。將麥穗魚eDNA擴增得到的PCR原液進行凝膠電泳并測序，目的條帶單一且明亮清晰，其測序結(jié)果在NCBI 網(wǎng)站上經(jīng)過 Blast 比對分析，其結(jié)果顯示與原始序列片段同源性達 100%。如圖1所示，目的條帶清晰明亮且無非特異性擴增。

通過實時熒光定量PCR擴增，閾值循環(huán)數(shù)（Cycle threshold，Ct值）平均值分別為4.83、12.75、14.53、20.09、24.56（圖2）。Ct值是指每個反應(yīng)管內(nèi)的熒光信號到達設(shè)定閾值時所經(jīng)歷的循環(huán)數(shù)。Ct值與模板的起始拷貝數(shù)的對數(shù)存在線性關(guān)系，起始模板濃度越高，Ct值越??；起始模板量濃度越低，Ct值越大。系統(tǒng)根據(jù)標準品濃度與擴增Ct值生成麥穗魚引物擴增的標準曲線。標準曲線用于確定實驗中未知樣本的目標DNA的相對或絕對量。標準曲線為：Ct=-3.56 x+38.89，（x為樣本中麥穗魚 eDNA 稀釋濃度的對數(shù)值），曲線的相關(guān)系數(shù) R2=0.995 2，擴增效率為90.94%，說明標準品的拷貝數(shù)與Ct值之間相關(guān)的可信度較好。反應(yīng)完成后，得到的熔解曲線如圖2所示。熔解曲線是指隨溫度升高 DNA 的雙螺旋結(jié)構(gòu)降解程度的曲線。本研究得到的溶解曲線出現(xiàn)單峰，說明反應(yīng)具有特異性。

2.3 麥穗魚eDNA濃度與生物量相關(guān)性曲線

根據(jù)公式Ct=-KlogX0+b計算可得（X0表示拷貝數(shù)，k 表示標準曲線斜率，b 表示標準曲線截距），0、2、4、8、16尾麥穗魚養(yǎng)殖缸中得到的平均拷貝數(shù)分別為0、275.56、354.03、596.67、1 264.28個·L-1。采用線性、對數(shù)和對數(shù)3種模型對麥穗魚生物量與eDNA濃度的相關(guān)性曲線進行了擬合。結(jié)果如圖3所示，3種曲線的擬合分析結(jié)果顯示：線性模型為麥穗魚生物量與eDNA濃度的最佳擬合模型，y = 72.087x+81.978，其R2高達0.933 8。eDNA濃度與麥穗魚生物量呈正相關(guān)性，即麥穗魚生物量越大，eDNA濃度也會越高。

3 討論與結(jié)論

3.1 引物和標準曲線

麥穗魚eDNA擴增后進行凝膠電泳并測序，目的條帶單一且明亮清晰，測序結(jié)果為麥穗魚。這說明基于COI基因設(shè)計的引物對麥穗魚DNA具有特異性擴增，且其擴增效率為90.94%，引物特異性好。

本研究得到的標準曲線為Ct=- 3.56 x + 38.89，曲線的相關(guān)系數(shù) R2=0.995 2，說明Ct值與標準品間呈現(xiàn)良好的線性關(guān)系，標準品的拷貝數(shù)與Ct值之間相關(guān)的可信度較好。結(jié)果表明，本研究在此麥穗魚eDNA濃度范圍內(nèi)具有良好的線性關(guān)系，所建立的標準曲線能夠正確反映出麥穗魚目的基因的擴增。

實時熒光定量PCR的SYBR Green Ⅰ熒光染料法其優(yōu)勢在于檢測方法簡便，檢測成本低。不足之處是染料與DNA雙鏈分子的結(jié)合是非特異性的，它可以和反應(yīng)體系中的所有DNA分子結(jié)合，易受到非特異性擴增和引物二聚體的影響，使定量結(jié)果不可靠[18]。解決這一問題可以借助熔解曲線分析法來區(qū)分非特異性產(chǎn)物和引物二聚體而排除假陽性[19]。熔解曲線如出現(xiàn)多峰，證明反應(yīng)非特異或存在二聚體。本研究得到的熔解曲線中出現(xiàn)單峰（圖2），說明定量結(jié)果是特異可靠的。

3.2 eDNA濃度與生物量相關(guān)性曲線

Takahara[9]等發(fā)現(xiàn)eDNA 濃度在第6天至第9天之間變化不大，得出eDNA濃度在第6天趨于穩(wěn)定的結(jié)論。因此，本研究在試驗第7天進行采樣時，eDNA濃度差異不明顯，得到的eDNA濃度與生物量的關(guān)系更為可靠。

通過3種曲線進一步擬合比較分析顯示：線性模型為麥穗魚生物量與eDNA濃度間的最佳擬合模型，其R2高達0.933 8，y=72.087x+81.978，表明在本試驗條件下，eDNA濃度與麥穗魚生物量呈正相關(guān)性。雖然在小型水族箱和人工池塘中eDNA濃度與生物量存在線性關(guān)系[20-22]，但在湖泊河流中eDNA與生物量的關(guān)系會更加的復(fù)雜，可能是線性或非線性。閆卉果[9]等在養(yǎng)殖桶中得到巖原鯉eDNA濃度與其生物量為線性關(guān)系；Coulter[23]等在河流中發(fā)現(xiàn)鰱魚eDNA濃度與其生物量存在非線性關(guān)系。本研究的擬合分析結(jié)果顯示：線性模型為麥穗魚生物量與eDNA濃度的最佳擬合模型，y=72.087x+81.978，其R2高達0.933 8，研究結(jié)果與養(yǎng)殖桶條件下目標魚eDNA濃度與其生物量間關(guān)系一致。

3.3 eDNA濃度的影響因素

eDNA技術(shù)作為一種新的監(jiān)測方法，具有便捷、環(huán)境友好、經(jīng)濟高效等優(yōu)點。eDNA產(chǎn)生后能在水體中保留一段時間，直至被降解不能檢測[24]。因此，它被認為提供了近似的實時數(shù)據(jù)。eDNA的降解與產(chǎn)生會導(dǎo)致eDNA濃度發(fā)生變化，但eDNA濃度的變化受多種因素的影響，進而影響eDNA濃度與其生物量的關(guān)系，主要包括生物因素和非生物因素2個方面。

生物因素方面體現(xiàn)在魚類的生物量會影響eDNA濃度。Michael等[25]發(fā)現(xiàn)，當大馬哈魚的密度越高，eDNA濃度也越高；Lacoursière等[26]研究表明，在同一溫度條件下，生物量越大，eDNA濃度越高。本研究結(jié)果發(fā)現(xiàn)，隨著養(yǎng)殖缸中麥穗魚的增加（0、2、4、8、16尾），eDNA濃度也呈正相關(guān)的線性關(guān)系。除生物因素外，其他非生物因素也會影響eDNA濃度，主要包括溫度、紫外線、pH等外界環(huán)境因素[27-30]。因此，在實際野外工作中需要綜合考慮各種因素，建立麥穗魚生物量和eDNA濃度關(guān)系模型才能準確評估麥穗魚的種群生物量。

本研究通過麥穗魚線粒體COI基因設(shè)計了麥穗魚特異性引物，F(xiàn)：CCGATTCTCATTGGGGGCTT，R： CCTGCGAGAGGGGGATAAAC，片段長度188 bp。確定了實驗室養(yǎng)殖麥穗魚eDNA濃度與其生物量的最佳擬合模型為線性模型，y=72.087x+81.978，R2=0.933 8。eDNA濃度與麥穗魚生物量呈正相關(guān)性。

參考文獻：

[1] 刁曹鋆， 王聞， 線薇薇， 等. 環(huán)境DNA技術(shù)在漁業(yè)資源生物量評估中的研究進展：現(xiàn)狀與展望[J]. 海洋科學， 2022， 46（2）： 135-144.

[2] 趙明， 趙夢迪， 馬春艷， 等. 環(huán)境DNA在水域生態(tài)中的研究進展[J]. 中國水產(chǎn)科學， 2018， 25（4）： 714-720.

[3] 趙彥偉， 陳家琪， 董麗， 等. 環(huán)境DNA技術(shù)在水生態(tài)領(lǐng)域應(yīng)用研究進展[J]. 農(nóng)業(yè)環(huán)境科學學報， 2021， 40（10）： 2057-2065.

[4] THOMSEN P F， KIELGAST J， IVERSEN L L， et al. Monitoring endangered freshwater biodiversity using environmental DNA[J]. Molecular Ecology， 2012， 21（11）： 2565-2573.

[5] FICETOLA G F， PANSU J， BONIN A， et al. Replication levels， false presences and the estimation of the presence/absence from eDNA metabarcoding data[J]. Molecular Ecology Resources， 2015， 15（3）： 543-556.

[6] JERDE C L， WILSON E A， DRESSLER T L. Measuring global fish species richness with eDNA metabarcoding[J]. Molecular Ecology Resources， 2019， 19（1）： 19-22.

[7] AHN H， KUME M， TERASHIMA Y， et al. Evaluation of fish biodiversity in estuaries using environmental DNA metabarcoding[J]. PLoS One， 2020， 15（10）： e0231127.

[8] DEJEAN T， VALENTINI A， DUPARC A， et al. Persistence of environmental DNA in freshwater ecosystems[J]. PLoS One， 2011， 6（8）： e23398.

[9] TAKAHARA T， MINAMOTO T， YAMANAKA H， et al. Estimation of fish biomass using environmental DNA[J]. PLoS One， 2012， 7（4）： e35868.

[10] KELLY R P， PORT J A， YAMAHARA K M， et al. Harnessing DNA to improve environmental management[J]. Science， 2014， 344（6191）： 1455-1456.

[11] 秦傳新， 左濤， 于剛， 等. 環(huán)境DNA在水生生態(tài)系統(tǒng)生物量評估中的研究進展[J]. 南方水產(chǎn)科學， 2020， 16（5）： 123-128.

[12] DOI H， INUI R， AKAMATSU Y， et al. Environmental DNA analysis for estimating the abundance and biomass of stream fish[J]. Freshwater Biology， 2017， 62（1）： 30-39.

[13] 楊麗萍， 胡俊儀， 秦超彬， 等. 基于線粒體Cytb分析入侵云南四大水系的麥穗魚群體遺傳結(jié)構(gòu)[J]. 水產(chǎn)學報， 2020， 44（3）： 339-350.

[14] 丁瑞華. 四川魚類志[M]. 成都： 四川科學技術(shù)出版社， 1994.

[15] 張鳳玉， 劉柳， 張靜， 等. 齊齊哈爾地區(qū)麥穗魚吸蟲囊蚴感染情況及其種類鑒定[J]. 中國寄生蟲學與寄生蟲病雜志， 2023， 41（1）： 112-116.

[16] 孫歡歡， 朱仁， 馮秀， 等. 西藏濕地外來麥穗魚生活史特征的適應(yīng)性研究[J]. 水生生物學報， 2022， 46（12）： 1780-1787.

[17] 劉東姣， 杜新科， 王世興， 等. 水蕹菜和羅莎生菜對麥穗魚養(yǎng)殖系統(tǒng)水質(zhì)的影響[J]. 水產(chǎn)學雜志， 2023， 36（1）： 89-94.

[18] 袁亞男， 劉文忠. 實時熒光定量PCR技術(shù)的類型、特點與應(yīng)用[J]. 中國畜牧獸醫(yī)， 2008， 35（3）： 27-30.

[19] 陳旭， 齊鳳坤， 康立功， 等. 實時熒光定量PCR技術(shù)研究進展及其應(yīng)用[J]. 東北農(nóng)業(yè)大學學報， 2010， 41（8）： 148-155.

[20] LACOURSI？魬RE-ROUSSEL A， C？魺T？魪G， LECLERC V， et al. Quantifying relative fish abundance with eDNA： a promising tool for fisheries management[J]. Journal of Applied Ecology， 2016， 53（4）： 1148-1157.

[21] SALTER I， JOENSEN M， KRISTIANSEN R， et al. Environmental DNA concentrations are correlated with regional biomass of Atlantic cod in oceanic waters[J]. Communications Biology， 2019， 2（1）： 461.

[22] 閆卉果， 董智玲， 馬婷婷， 等. 基于環(huán)境DNA的巖原鯉檢測及生物量評估[J]. 水產(chǎn)學報， 2022， 46（6）： 1018-1026.

[23] COULTER D P， WANG P， COULTER A A， et al. Nonlinear relationship between silver carp density and their eDNA concentration in a large river[J]. PLoS One， 2019， 14（6）： e0218823.

[24] PIAGGIO A J， ENGEMAN R M， HOPKEN M W， et al. Detecting an elusive invasive species： a diagnostic PCR to detect Burmese python in Florida waters and an assessment of persistence of environmental DNA[J]. Molecular Ecology Resources， 2014， 14（2）： 374-380.

[25] TILLOTSON M D， KELLY R P， DUDA J J， et al. Concentrations of environmental DNA （eDNA） reflect spawning salmon abundance at fine spatial and temporal scales[J]. Biological Conservation， 2018， 220： 1-11.

[26] LACOURSI？魬RE-ROUSSEL A， ROSABAL M， BERNATCHEZ L. Estimating fish abundance and biomass from eDNA concentrations： variability among capture methods and environmental conditions[J]. Molecular Ecology Resources， 2016， 16（6）： 1401-1414.

[27] JO T， MURAKAMI H， YAMAMOTO S， et al. Effect of water temperature and fish biomass on environmental DNA shedding， degradation， and size distribution[J]. Ecology and Evolution， 2019， 9（3）： 1135-1146.

[28] EICHMILLER J J， BEST S E， SORENSEN P W. Effects of temperature and trophic state on degradation of environmental DNA in lake water[J]. Environmental Science & Technology， 2016， 50（4）： 1859-1867.

[29] STRICKLER K M， FREMIER A K， GOLDBERG G S. Quantifying effects of UV-B， temperature， and pH on eDNA degradation in aquatic microcosms[J]. Biological Conservation， 2015， 183： 85-92.

[30] LANCE R F， KLYMUS K E， RICHTER C A， et al. Experimental observations on the decay of environmental DNA from bighead and silver carps[J]. Management of Biological Invasions， 2017， 8（3）： 343-359.