 * g1 g& N: |0 ?; o8 W( T
中文题目：人类活动对中国沿海鱼类生物多样性的整体影响评价 Q8 j! X( O; c) C' E. {. U3 v
英文题目：Holistic Impact Evaluation of Human Activities on the Coastal Fish Biodiversity in Chinese Coastal Environment
* p* x) y0 q- K2 y& G O' L 期刊：Environ. Sci. Technol. 2022 - h4 V0 o1 ^/ J( V9 g3 K) n
第一作者：Wenjun Zhong 5 H. ]5 @- j2 t' R
通讯作者：Hongjun Li, Xiaowei Zhang
- r* A9 C# c1 a4 y+ m 作者单位：南京大学环境学院；国家海洋环境监测中心
: e. c" T* Q p/ d 方法解读：张政
& D1 ?$ h& `2 Z1 A1 l' x 正文部分转自公众号【国家海洋环境监测中心】 * S* \) d9 Y; s1 K& o" P( V
在此基础上，我们对文章的研究方法进行了详细解读  ' R1 J5 y* W+ l; q
正文解读 1 f9 f- }. y+ ^
环境DNA（Environmental DNA, eDNA）是指从环境样品（如土壤、沉积物、空气、水体等）中提取的DNA，不对任何目标生物进行分离。在海洋生态系统中，通过采集水样和沉积物，收集微型生物的生物个体（个体DNA）和大型生物的碎片或游离DNA（非个体DNA），得到各类海洋生物的环境DNA，从而进行物种或群落的分析。海洋中心李宏俊研究员团队联合南京大学和汕头大学等单位在鸭绿江口和南澳岛海域开展了鱼类环境DNA研究，为评估人类活动对我国近岸鱼类多样性影响提供案例，相关研究成果发表于环境领域顶级期刊Environmental Science & Technology。
6 t2 L+ b% P. ?2 l! b; f D0 P( E- Z0 j1 q, q- R. O O
随着人类高强度的活动和气候变化，各种人类活动严重威胁近岸海域生态系统质量和稳定性，但是由于缺乏有效的生物多样性监测和评价方法，极大限制了针对不同海洋生态系统环境压力影响的评估。本研究通过整合eDNA技术和生态功能特性，建立了一套评估环境污染和水产养殖对我国近岸海域鱼类生物多样性影响的综合性方法。通过建立我国鱼类12S rRNA条形码参考数据库，提高了eDNA技术对我国海洋鱼类群落分类学注释的准确性，同时分别在序列、分类学和功能尺度计算鱼类群落的多样性，在多个维度提供鱼类群落多样性的分析结果。选择黄海鸭绿江口和南海南澳岛作为代表性研究区域，对本文开发的eDNA生物监测评价工作流程的适用性进行了验证（图1）。 $ Z, o/ [4 g) k- Y5 Q n4 W5 Y
; D: L- Q! o( X9 Q8 h6 V8 A 图1 研究框架 * P. a1 m* t& [7 _" f, t
本研究利用便携式水环境DNA过滤仪搭配富集器采集表层海水，通过DNA提取、PCR扩增、DNA文库构建和高通量测序比对等步骤，在鸭绿江口检测到鱼类物种14目24科43属46种，环境DNA所检出鱼类有76.1%可被鸭绿江口历史渔获物种验证，监测结果还表明，处于海洋保护区中的鱼类丰富度显著高于非保护区。南澳岛区域鱼类物种的丰富度是鸭绿江口的1.8倍，南澳岛中低营养级的鱼类的检出率是鸭绿江口地区的1.64倍（图2）。 . _9 L7 z+ v2 h ?
* c y+ Y* Q9 U c2 _
图2 研究区域不同营养级鱼类物种多样性 ; g! N5 [: v2 Z x1 s7 F
; l8 p9 L, o/ j: Z 随着采样站位与岸边的距离逐渐增加，鸭绿江口地区鱼类序列多样性首先快速升高随后保持稳定，而南澳岛地区鱼类序列多样性呈现持续升高趋势。鱼类分类学多样性和功能多样性在鸭绿江口随着离岸距离的增加逐渐降低，而在南澳岛地区，其呈现先增加后降低之后再增加的变化趋势（图3）。此外，无论是在南澳岛还是鸭绿江口地区，鱼类分类学多样性和功能多样性具有显著相关的一致性变化趋势。在鸭绿江口地区，河口污染因素（COD和Zn）导致鱼类群落序列多样性明显下降。在南澳岛地区，养殖区域的鱼类群落序列多样性要低于非养殖区域（图4）。
8 i* d3 a% T- Z' k! P 
% U! c' c& p& O! E; W7 } 图3 研究区域鱼类多样性空间分布规律
$ a! i7 s% X2 B4 Q, i$ o 
" t. k- H: ~8 c' y, x x 图4 人类活动对研究区域鱼类多样性影响机制 4 V& l* |6 V* y
本研究认为，目前eDNA研究海洋生物多样性，不能局限在注释到的具体物种，eDNA也不可能检测出来所有的物种。首先参考数据库目前还并不完整，有很多物种目前还没有对应的barcode序列，参考数据库中没有序列就可能检测不到；其次是无论选择哪一段目标序列，都无法做到准确区分所有物种，肯定会存在一些在测序区域序列相似度特别高甚至序列一致而无法得到有效区分的物种。构建针对目标研究区域的个性化参考数据库可能会部分解决注释错误的问题，但是肯定也没办法做到100%的准确，就像本研究专门建立了一个我国鱼类12S rRNA的参考数据库，确实能够提高分类学注释的准确性，但是从结果中也可以看出，eDNA结果和历史观测数据的一致性在鸭绿江口海域为76.09%，但是在南澳岛海域也仅有35.8%。
) k! _* w8 l5 d 因此，在进行大型生物eDNA的研究中，应该采取弱化具体物种的研究思路，使用一些综合性的评价指数，比如alpha多样性指数、beta多样性指数、功能多样性，或者是一些生态功能分类，比如营养级、生活习性、食性等，来代替具体的物种以达到研究目的。 : u- A+ o& b. \$ I7 {' y1 N
下一步，课题组将开展我国近海环境DNA“生命之树”扫描式监测，为摸清我国海洋生物多样性本底提供基础数据，开展不同生物类群生态响应和阈值研究，为构建我国近海生态质量评价标准和业务化监测提供科学依据。
9 O: g1 U0 ^# y9 u0 X, W7 o: D3 C 方法详解 3 H" o" I8 }: q; T
1. 对引物Teleo进行in silico评价 ! j+ K/ x* J+ w. ~
（1）从EMBL下载序列信息
9 c' ^: Z; B+ l% I5 ]3 ?+ u （2）使用引物Teleo针对下载并整合分类信息的序列进行in silico分析，只保留可以被该引物比对出的鱼类信息
7 ^0 Z1 H6 f6 t% x) ]& D) E( i) X （3）分析这些鱼类的世界和中国海域的分布情况，从而说明该引物在中国沿海地区的鱼类检测效力
& D* M6 w! u0 j+ h" P 2. 构建中国鱼类12S-DNA条形码参考数据库
3 H4 p7 w! G/ T, _: ^, e$ W （1）从fishbase.org获得鱼类分类信息
# R( |9 f/ n% N& x+ Z （2）根据分类信息，从NCBI Taxonomy Browser获得每种鱼类的“taxid” " Q& h4 Y# f7 E* t0 S
（3）根据“taxid”搜索获得每种鱼的12S-rDNA序列
. l# f7 C+ h; n2 o. U （4） 如果某鱼种没有公开的12S-rDNA序列但有全基因组信息，从该鱼种的全基因组中选择相应的12S-rDNA；如果一个物种有2个及以上的12S序列时，比较后去重复 - N3 S: ~, w, D# |2 U
3. 获得中国鱼类性状数据
/ d) l* \0 z& k1 m$ O; a2 r7 j0 w0 \ （1）从fishbase.org 获得迁移类型、摄食类型和营养级指数等性状信息；在 globalbioticinteractions.org上筛选鱼类捕食关系，并将其映射为鱼类食物网
- }8 j, E, A: C& O2 `0 I' m3 J （2）从rfishbase 获得洄游等信息
4 F/ P' z; `$ X4 t* z （3）从fishbase.org 下载具有最大体长数据的整条成鱼的高质量图片和科学侧视图，然后是针对体型变量使用对侧视图图片进行评估。
P* R0 u4 N" a4 ~# z! n: z! s 4. eDNA取样、扩增和测序 & k/ O2 T! W8 l0 Z- W8 c
（1）在鸭绿江口和南澳岛海域收集海水，过膜
+ ~: c: P7 k# R0 } （2）使用DNeasy Blood & Tissue Kit提取DNA（66个海水样品、3个空白滤膜、3个提取空白和3个阴性对照）
5 M- o9 T9 |/ {3 M9 i （3）使用12S-rRNA引物Teleo扩增 & n# u; c5 a5 N6 I
（4）使用Ion Torrent S5测序仪进行测序 4 ~4 q% o' L7 C! W2 B. |' O
5. 生物信息学分析和分类注释
, H) }+ t% U: }6 t: ^8 P1 x" j) X7 E （1）针对测序数据使用USEARCH的--unoise3命令生成ASV
; X& n7 w% O9 l0 @ （2）只保留reads大于10的ASV，并只保留注释到种水平的ASV 5 q9 d2 L+ w; K1 D: ~
（3）只有在每个站点至少2个重复中出现的ASV才被认作存在，并且对该站点的重复进行平均后四舍五入来合并该站点的ASV reads
! Y! i$ I4 m! N 6. 序列、分类和功能多样性的计算 . I$ A4 }. s# S) e, E
（1）对特定区域内同一物种的所有ASV序列进行多序列比对，仅保留重叠率>50%的序列 , \* i/ H0 x5 q7 w
（2）通过跨地区重叠物种的ASV数量和物种内的平均序列变异来进行序列多样性的比较 . y( I/ b7 s% @' @7 n- e- B
（3）根据 计算核苷酸多样性 ( }3 C1 a) @8 {! N# o( u( Z
（4）mFD package计算每个采样点的functional richness (FRic)
, q2 A) B# Q- D* e2 T$ Y6 B 7. 比较区域之间和区域内的多样性变化 1 l- n' z' m" L- O
（1）序列多样性差异进行Wilcoxon秩和检验 ) c; X& R+ N: }% o, N
（2）β多样性比较
6 M2 ^# M; |9 N7 C0 {, ~ 8. 针对不同养殖区域评价河口污染对海岛养殖的影响 0 Y) J- c8 c; j, m
（1）Mantel检验用于确定与鸭绿江口不同鱼类多样性方面显著相关的环境因素 ( U8 ]( U5 P8 F) m, ^
（2）采用Wilcoxon秩和检验对南澳海域3类养殖区的序列、分类和功能多样性进行检验 (Zhong et al., 2022)
1 l" z8 h( q7 z1 f 参考文献
- E" S$ |7 \8 ~8 D7 D z/ C 参考文献： 4 N6 Y2 T$ o* f# u
Zhong, W., Zhang, J., Wang, Z., Lin, J., Huang, X., Liu, W., ... Zhang, X., 2022. Holistic Impact Evaluation of Human Activities on the Coastal Fish Biodiversity in the Chinese Coastal Environment. Environ Sci Technol. 2 y; X# x F( ~3 v" N1 f' y
原文链接： % U1 o3 [3 y; a, v3 f; w* t' W C
https://pubs.acs.org/doi/10.1021/acs.est.2c01339#
; e7 Q2 Q+ P! B+ x! i( O! h  
) z& K T {% m6 d- f' F' U9 \ 中国科学院生态环境研究中心 , v- X9 |' n+ J
8 M! R8 R, J* C L& @$ g$ t
环境生物技术重点实验室
! O2 }3 R" D x5 m0 g- }8 ~ 邓晔 研究员课题组发布 
r' {% }2 l/ b' n0 }8 K 公众号: 5 d' T$ C" P! J( `4 w. U' A
微生态笔记
Q& L9 U# O, F' ^ 微信号:
! i7 G4 k7 e3 u7 r5 d mem_denglab
! L; ^6 y$ h" U; @ ● 扫码关注我们
6 L* N+ E! P/ ?. j1 ^, B9 t
' ?; g q' n1 E5 \# r
& T) Q+ d. I5 V( d( t% o: j! S, s( E4 I! N
! y, h. `$ p y, B% G8 p: Z2 j | 
