生态毒理学报
Vol.1,No.3
Sep.2006
2006年9月
AsianJournalofEcotoxicology
环境中金属生物有效性的预测模型———生物配体模型研究进展
王学东
1,2
,马义兵
1,*
,华珞2,韦东普1,李波
1
1.中国农业科学院农业资源与农业区划研究所,北京100081
2.首都师范大学资源环境与旅游学院资源环境与地理信息系统北京市重点实验室,北京100037
摘要:生物配体模型(BLM)是一种用于预测环境中金属生物毒性的机理性模型.模型理论起源于自由离子活度模型(FIAM)和鱼鳃络合模型(GSIM),考虑了自由金属离子的活度以及自然环境存在的其他离子(如Ca2+、Na+、Mg2+、非生物配体(如可溶性有机质、氯化物、碳酸盐、硫酸盐)和生物配体的竞争.目前,在水生生态系统中,基于鱼H+)、
鳃络合模型的框架基础,通过生物化学实验手段并结合数学方法,建立了预测铜、锌、银、镍对Rainbowtrout(虹鳟鱼)、慢性毒性的BLM版本,并积极探索其在陆地生态系统Fatheadminnow(黑头呆鱼)和Daphniamagna(水蚤)的急、
中的应用.虽然生物配体模型在实验室模拟条件下取得了较为满意的结果,但其中包含着一些假设,在实际应用中还具有一定的局限性,尤其是陆地生态系统生物配体模型的发展还需要做许多研究工作.本文主要论述了生物配体模型的理论基础、实现手段和应用情况,讨论了生物配体模型的优势和局限性并对其未来研究方向进行了展望.关键词:生物配体模型;重金属;生物有效性/毒性
文章编号:1673-5897(2006)3-193-10
中图分类号:Q141,X171.5
文献标识码:A
AdvancesinBiotic-LigandModeltoPredicttheBioavailabilityofMetalsinEnvironments
WANGXue-dong1,2,MAYi-bing1,*,HUALuo2,WEIDong-pu1,LIBo1
1.InstituteofAgriculturalResourcesandRegionalPlanning,ChineseAcademyofAgriculturalSciences,Beijing100081
2.TheKeyLabofResourceEnvironmentandGIS,CollegeofResourceEnvironment&Tourism,CapitalNormalUniversity,Beijing100037
Received28June2006
receivedinrevisedform14August2006
accepted15August2006
Abstract:TheBiotic-ligandModel(BLM)isamechanisticmodelusedtopredicttheecotoxicityofmetalsinenvironments.ThetheoryoftheBLMevolvesfromFree-Ion-ActivityModel(FIAM)andGillSurfaceInteractionModel(GSIM),andincorporatesthecompetitionofthefreemetalionwithothernaturallyoccurringcations(e.g.Ca2+,Na+,Mg2+,H+),togetherwithcomplexationbyabioticligands(e.g.DOM,chloride,carbonates,sulfide)forbindingwiththebioticligand.Usingbiochemistryandmathematicsmethods,differentversionsofacuteandchronicBLMforthetoxicityofCu,Zn,AgandNitofish(Rainbowtrout,Fatheadminnow)andDaphniamagnainaquaticecosystemshadbeenpublished.AlthoughtheBLMwasfoundtobesatisfactoryinlaboratoryconditions,buttherearesomelimitationsinpractice,especiallyinterrestrialecosystems.Inthisarticle,theoreticalbasis,development,applicationandthelimitationsofBLMwerereviewed,meanwhiletheinvestigationinthefuturewereprospected.
Keywords:biotic-ligandmodel/BLM;heavymetals;bioavailability/ecotoxicity
收稿日期:2006-06-28修回日期:2006-08-14录用日期:2006-08-15
基金项目:国家自然科学基金项目(No.40571071);国际铜业协会项目(No.E-EU-05-03B;E-AS-05-04)作者简介:王学东(1978—),男,博士研究生;*通讯作者(Correspondingauthor),E-mail:ybma@caas.ac.cn
194
生态毒理学报第1卷
随着工业现代化、农村城镇化和养殖业集约化的发展,环境中重金属污染日趋严重,环境管理和监测所需的基于环境风险评价的环境质量标准制定越来越引起人们的重视.当前的重金属环境质量标准和风险评价方法主要是建立在金属总量的基础上(DeSchamphelaereetal.,2002a),忽视了环境介质(水、土壤)对金属生物毒性的影响.如何用比较简单的方法来准确地评估和预测环境中重金属的毒性呢?在过去几十年里,科学家开发出一些能够预测重金属生物毒性的数学模型,如自由离子活度模型(Free-Ion-ActivityModel,FIAM)(Morel,1983)、鱼鳃络合模型(GillSurface
困难,生物配体模型的研究还面临着很大的挑战.本文着重从生物配体模型的理论基础、模型的实现、模型的应用、局限性和未来的主要研究方向对生物配体模型进行了论述,目的是为制定适宜的环境质量标准和开展重金属污染风险评价提供理论依据和实现手段.
11.1
生物配体模型的理论基础重金属形态和生物毒性的关系
环境中重金属的形态与其生物有效性/毒性
有直接的关系,并不是所有形态的重金属都具有毒性,例如,水溶性铜的形态就有Cu2+,Cu(OH)+,
InteractionModel,GSIM)(Pagenkopf,1983)以及在此基础上发展而来的生物配体模型(Biotic-Ligand
Cu-DOC等,Cu-DOC占溶解性铜的大部分,但一般情况下不具有毒性(Paquinetal.,2000).(1995)根据大量相关文献指出金属毒性Campbell
直接和水中自由金属离子活度有关.在土壤溶液中,即使大部分铜以有机络合物的形式、大部分镉以无机络合物的形式存在,铜和镉对植物的毒性仍与溶液中自由金属离子活度显著相关(Bingham
Model,BLM)(Meyeretal.,1999).其中生物配体模型将生物受体位点作为生物配体,考虑了影响生物毒性的水化学性质,并把生物有效性的概念引入到水质标准中,在较宽的模拟水质范围内取得了较好的预测效果(Niyogietal.,2004).美国环保署对生物配体模型进行了评估,根据评估结果已经决定将该模型作为制定关于金属元素国家水质标准的基础.在欧洲,欧盟也在考虑把生物配体模型用于制定水质纲要等方面(黄圣彪,
etal.,1983;Sauvéetal.,1996).后来的研究指(OH)+,Me2(OH2)出,可溶性金属形态中的Me
2+
,
MeCl+等无机络合形态(Markichetal.,2003;McLaughlinetal.,1997;Smoldersetal.,1998;Saekietal.,2002)甚至一些有机络合形态(如铜)(Poldoski,1979)也可以被生物直接吸收而产生毒害.但是,关于Me(OH)+的毒害作用到底是由于(OH)+的直接作用目前还没有定pH的原因还是Me
论.总之,预测金属生物毒性的模型首先需要考虑金属的存在形态和生物毒性的关系.
2003).目前,生物配体模型已经被成功地用于解释铜对一系列具有不同敏感性生物体的毒性,模型在其他金属如银、镉、锌、镍、钴和铅上的应用也在建立和发展之中.同时研究人员正在积极探索生物配体模型在陆地生态系统和沉积物中的应用———即建立t-BLM(Santore,2006).
生物配体模型的提出具有重要的理论和现实意义.首先,它提供了一种预测水环境中金属铜毒性效应的简便、科学的方法,有利于实现对水体生态系统的监测和保护;其次,生物配体模型建立在铜对鱼的毒性机理基础上,不仅适用于铜对鱼的毒性,也适用于其他金属和其他生物,理论上也可以应用到底泥和陆地生态系统中,因此模型具有广阔的应用空间;最后,生物配体模型综合了金属吸收的化学、生理学、生物学等方面的成果,这些学科的发展也为模型的完善提供了基础.但是生物配体模型在实现过程中存在许多假设而且有些影响金属毒性的因素没有考虑,这使模型的应用具有一定的局限性,如陆地生态系统中土壤、生物、金属相互作用的复杂性使得模型的预测更加
1.2生物配体模型的理论框架
(1983)基于金属对生物的毒性取决于自Morel
由金属离子活度提出了预测金属毒性的自由离子活度模型,其机理为自由金属离子可以和细胞表面的物理活性位点结合,然后跨过细胞膜对生物产生毒性.该模型在理论上解释了金属的生物毒性取决于自由金属离子活度的原因,但研究表明:当水的硬度不同,即使水中自由金属离子活度相同,金属的生物毒性也不相同(Hunn,1985;Meyer(1983)提出的鱼鳃络合模etal.,1999).Pagenkopf
型成功地解释了这一现象.鱼鳃络合模型认为金属对鱼的毒性作用位点为鱼鳃,鱼鳃表面存在着带有负电荷的配体,金属离子与配体反应产生毒
第3期王学东等:环境中金属生物有效性的预测模型———生物配体模型研究进展
195
性,水中的阳离子可以和金属竞争鱼鳃的结合位点而减弱金属的毒性.许多试验也证实了阳离子竞争在减弱金属毒性中的作用(Ericksonetal.,
H+Ca2+Mg2+
1996;DeSchamphelaereetal.,2002b;2004a).同自由离子活度模型相比,鱼鳃络合模型不仅考虑了金属形态而且还考虑了水的硬度和H+浓度对金属生物毒性的影响,模型的应用范围更为广泛.但是鱼鳃络合模型仅把鱼鳃作为金属毒性作用位点,而没有计算出预测金属生物毒性的一个关键指标EC50/LC50值(半数抑制/致死浓度),不能替代水效应比(WaterEffectiveRatio,WER)来对水质标准进行修改.同自由离子活度模型和鱼鳃络合模型相比,生物配体模型更加全面地考虑了影响金属生物有效性的因素,能够准确地计算金属(如铜)在作用位点上的积累水平和EC50/LC50值,从而真实地预测环境中金属的毒性效应水平.
生物配体模型最早用来预测铜对鱼的毒性,模型认为铜对鱼的毒性是由于减少了鱼体中钠离子和氯离子的浓度引起的.这是因为铜离子和鱼腮表面的基团络合后累积在鱼鳃表面,占据了细胞表面的钠离子通道,从而抑制了细胞对氯离子和钠离子的吸收或细胞膜上钠/钾离子酶的活性(Janesetal.,1995).能够络合金属的生物部位(生物膜)称为生物配体,因此模型被称为生物配体模型.生物配体模型考虑了影响金属生物毒性的3个因素,即浓度、络合和竞争.金属的生物毒性取决于自由离子活度(后来发展的模型又考虑了其它的金属形态),而自由离子活度与总金属和溶解性金属的浓度有关,同时受有机无机配体络合的影响.硬度阳离子、H+等和自由金属离子竞争生物配体的作用位点.经过这些作用以后的自(M2+)和生物配体(BL)结合形成金属-由金属离子
生物配体络合物(MBL),MBL跨过生物膜后产生生物效应(见图1).
在生物配体模型的理论基础上,借助化学分析手段,
利用化学平衡模型(WHAM、CHESS、
DOC
金属M
自由金属
离子M2+
吸附解吸
MBL
同化
MeOH+MeCO3
图1生物配体模型中导致金属生物毒性的重要理化
(H+、反应过程Ca2+、Mg2+是参与竞争的阳离子;Me-DOC代表金属有机络合物;MeOH+、MeCO3代表金属
无机络合物;BL代表生物配体)
Fig.1
TheimportantphysicochemicalprocessesleadingtothetoxicityofheavymetaltobioticofBLM
自由金属离子的活度和环境水化学性质如
DOC、pH、Ca2+、Mg2+、CO3等有关,环境中这些有机(如DOC)无机(如CO3)配体和金属发生络合反应,而Ca2+、Mg2+等阳离子和金属竞争配体的结合位点.模型假设这些反应处于平衡状态,简单的(为简化方程,没有考虑离子的平衡方程表示如下电荷数):
M+L!ML
[ML]=KML[M][L]
(1)
2-
2-
无机配体;KML:M:溶液中自由金属离子;L:有机、
金属和配体反应的平衡常数.式中溶液中自由金属离子及配体的浓度单位通常为mol・L-1,平衡常数的单位为L・mol-1,后面公式中如无特别标注,离子的浓度单位均为mol・L-1,常数K值的单位均为L・mol-1.
当存在其他阳离子的竞争作用时:
Me+L!MeL
[MeL]=KMeL[Me][L]
(2)
Me:参与竞争的阳离子;KMeL:配体和阳离子反应的平衡常数.
目前,(Chemical
一些化学平衡模型如MINEQL
EquilibriumModelingSystem),CHESS
(ChemicalEquilibriumofSpeciesandSurfaces)和(Winder-mereHumicAcidModel)被包含于WHAM
生物配体模型中用于计算自由金属离子活度(表1).
金属生物毒性的预测(KS、EC50/LC50值的计算):
金属和生物配体的相互作用过程如下:
kfkint
M+BL!MBL
kd
Mint+BL
(3)MINEQL+)并结合数学方程(Michaelis-Menten、Langmuir),建立生物配体模型,计算出反映金属对生物毒性强弱的指标如金属和配体的络合常数(KS),用于预测重金属对生物的EC50/LC50值.
1.3生物配体模型的实现自由金属离子活度的计算:
{MBL}=Ks{BL}[Mn+]
BL:生物配体;kf、kd、kint:金属和生物配体的
196
生态毒理学报第1卷
表1金属形态计算模型、输入水质参数及应用举例
(Paquinetal.,2002)
程度(EC50/LC50)与结合金属的活性位点和金属总(fMBL)具有直接的关系,当考虑阳结合位点的比例
离子竞争作用时,它们之间的关系表示为:
[M]KS
fMBL={MBL}=
{BLTOT}1+KS[M]+\"iKC[iCi]
(9)
Table1
Inputparametersforspeciationcodesusedin
publishedBLMapplications
水质参数
应用举例
2-
模型
WHAM-V/CHESS
Ca2+,Mg2+,Na+,K+,SO4,Cl-,pH,碱度,DOC,%
腐殖酸
DiToroetal.,2001;Santore
etal.,2001;De
Schamphelaereetal.,2002a
McGeeretal.,2000;Schwartzetal.,2001
fMBL:金属和生物配体的结合部分;{MBL}:金
属-配体络合物浓度;{BLTOT}:络合容量(mol・g-1).
当对50%的生物产生影响时,方程可以写为:
(M)=EC50
50
MINEQL+
Ca,Mg,Na,pH,DOC,
2+
2+
+
S2O3,Cl-,SO4
2-2-
fMBL1-fMBLKS
501+\"K#i
Ci
[Ci]
$(10)
形成、离解和同化速率常数(mol・cm-2・s-1);Ks:金属和配体反应的平衡常数.当考虑阳离子竞争作用时:
Me+BL!MeBL
[MeBL]=KMeBL[Me][BL]
(4)
EC50:在一定时间范围内对50%的生物产生
影响时自由金属离子的浓度(mol・L-1);fMBL:对
50
50%的生物产生影响时金属-配体络合物浓度占络合容量的比例.
应用这个方法不需要知道金属和配体络合的浓度,所以当不能准确得到金属-配体络合物浓度时,它们之间的平衡常数通常使用这种方法来求.
金属累积试验结合Langmuir吸附方程也可以求得金属和配体的络合常数,在这个方法中和配
Me:参与竞争的阳离子;KMeBL:阳离子与生物配体反应的平衡常数.
KS、EC50/LC50和金属在生物配体的累积量通常通过金属膜通量、金属累积、金属毒性试验结合多Michaelis-Menten方程、Langmuir等温吸附方程、(Slaveykovaetal.,2005).元非线性回归方法来求
通过膜通量的方法可以求得金属和配体的平衡反应常数.其中细胞膜对金属的吸收速率通常通过对不同时间内金属累积曲线的斜率求得,穿过膜的重金属通过化学提取技术以区别于重金属的累积.它们之间的关系用Michaelis-Menten方程表示为:
Jint=Jmax
[M]
KM+[M]-2
-1
-1
体结合的总金属浓度{MBLTOT}通常通过生物累积或滴定的方法来求,考虑阳离子竞争的影响,方程表达式如下:
{MBLTOT}={MBLTOT}max
]K[SM
]+\"iKC[iCi]1+K[SM
(11)
{MBLTOT}:和配体结合的金属总浓度(mol・g-1);{MBLTOT}max:配体能够结合的金属最大浓度(mol・g-1).
通过这个方程可以求得金属和生物配体的总结合位点(络合容量)及它们之间反应的平衡常数.当考虑配体结合位点具有不同的特征,如具有多电子和多功能团的时候经常应用NICA-(5)
Jint:吸收速率(mol・cm・s);Jmax:最大膜通量(mol・cm・s);KM:半饱和常数(mol・L)
-2
-1
如果金属穿过膜受速率,那么金属和生物配体的络合常数(KS)可以通过KM的倒数来求:
KM=kint+kd
kf
Donnan方法.
利用上述方法在水生生态系统中获得了大量关于金属生物毒性的数据,并建立了不同金属的不同预测体系,而在陆地生态系统中主要累积了一些相关的毒性数据.下面对生物配体模型在水生和陆地生态系统中的应用进行简要的阐述.
(6)
如果kint<<kd,那么,当考虑其他阳离子竞争作用的影响,方程可以写为:
KM=kd=1
kfKS
(7)
则:
2
(8)
生物配体模型的应用水体生态系统
在水生生态系统中,通过生物配体模型的构
Jint=Jmax
KCi:
[M]KS
[M]+\"iKC[1+KSiCi]
2.1
阳离子和配体反应的平衡常数;[Ci]:参
建实现了不同金属毒性的预测体系,然后转为计算机软件加以推广和应用.由于预测金属毒性采取的方法和选取的生物可能不同,所以对于同一
与竞争的阳离子浓度.
在生物毒性终点和生物测定方法中假设毒性第3期王学东等:环境中金属生物有效性的预测模型———生物配体模型研究进展
197
金属出现了不同的预测版本.现已建立了多个铜、锌、银的BLM版本和一个镍的BLM版本用来预测金属对鱼(Rainbowtrout虹鳟鱼;Fathead(Daphniamagna)的毒性.minnow黑头呆鱼)或水蚤
目前公开出版的已有4个Cu-BLM版本,其中包括3个预测铜急性毒性版本:Ia(Santoreet
IIIa(Buryetal.,2002),没有慢性Ag-BLM版本.相对于铜和银,Zn-BLM发展较为滞后,当前有2个急性Zn-BLM版本Ia(Santoreetal.,2002)、IIa(Heij-ericketal.,2002)和1个慢性毒性Zn-BLM版本Ic(DeSchamphelaereetal.,2004b).Ni-BLM仅有1个急性毒性版本Ia(Meyeretal.,1999;Water
al.,2001;DiToroetal.,2001)、IIa(DeScham-phelaereetal.,2002b)、IIIa(DeSchamphelaereet(DeSchamphe-al.,2002a)和1个慢性毒性版本Ic
EnvironmentResearchFoundation,2002).还没有镉、铅、钴的BLM版本发行,但它们均有GSIM版本(Playleetal.,1993a;MacDonaldetal.,2002;Richards
laereetal.,2004a).急性Ag-BLM有3个版本:Ia(Paquinetal.,1999)、IIa(McGeeretal.,2000)和
表2
etal.,1998),这为建立相应的BLM提供了基础.所有这些版本选用的生物及相关K值见表2.
不同金属不同BLM和GSIM版本中金属、阳离子和配体的络合平衡常数(logK)
Table2
Affinityconstants(logK)forBL-CationandM-BLcomplexesusedinthedifferentversionsofBLMand
Gill-Mmodel
版本
生物
金属
KBL-Ca3.63.53.5NA3.32.32.32.34.83.33.64.05.04.04.7
2+
KBL-MgNA3.63.6NA----NA2.42.4NANA4.0NA
2+
KBL-Na3.03.23.22.94.72.32.92.9NA2.42.43.0NA3.53.2
+
KBL-H5.45.45.46.75.94.35.95.96.7NA6.37.56.74.06.2
+
KBL-Mn7.48.08.08.010.07.37.68.95.55.35.54.08.66.05.1
+
KBL-Cu(OH)
NA7.47.48.0-----------
+
KBL-CuCO
3
KBL-AgCl
----NA6.7NANA-------
-
CuⅠaⅡaⅢaⅠc
DMDMDMDMRBTDMRBTDMRBTDMFHMFHMFHMRBTRBT
NANA7.07.4-----------
Agg-AgⅠaⅡaⅢa
ZnⅠaⅡaⅠc
NiCdPbCo
Ⅰag-Cdg-Pbg-Co
注:表中所有的络合系数均为以10为底的对数值;a代表急性毒性;c代表慢性毒性;NA代表未检出;DM:Daphniamagna
(水蚤);RBT:Rainbowtrout(虹鳟鱼);FHM:Fatheadminnow(黑头呆鱼)
这些BLM版本的参数有的来源于之前建立的鱼鳃络合模型,有的则利用前边所述的方法求得.如铜的Ia中阳离子和配体的络合系数(logKBL-Me)来源于Playle等(1992;1993b)的鱼鳃络合模型,铜和有机物的络合反应平衡常数采用WHAM-V计算(Tipping,1994),
其他的常数取自CHESS
(Santoreetal.,1995).IIa、IIIa的logKBL-Me来源于金属累积和毒性的数据,铜和有机、无机的络合常数取自Martell等(1997)和WHAM-V数据库(Tipping,
和配体的络合常数来源于放射性同位素的方法(Alsopetal.,2000).对于同一种金属的不同版本它们的不同点主要集中在以下4个方面:对金属毒性机理的认识;模型包含的金属毒性形态;对
DOC的假设;与金属竞争配体络合位点的的阳离子种类,由此导致了模型预测结果的差异.在Cu-BLM的IIa和IIIa版本中考虑了Mg2+的竞争作用,
并求出了logKBL-Mg值,而Ia没有考虑.版本Ia在计
2+
算中假设所有的DOC都是活性的,其中包括90%的HA(胡敏酸)和10%的FA(腐殖酸),而版本IIa和IIIa假设DOC有50%是活性的,并以FA的形式和铜反应.另外版本Ia中,水的pH范围在7.0 ̄
1994).而Ag-BLM3个版本的参数基本上都来源于Janes等(1995)创立的鱼鳃络合模型.当然也可以采用其他方法获取参数,如锌的Ia版本中锌
198
生态毒理学报第1卷
(HCO3)+是水中铜的优势形态,其平8.5之间,Cu衡常数KCu(
HCO3)
+
年龄和体形大小影响镍对Fatheadminnow的毒性,这可能是由于镍影响了鱼的呼吸,而小鱼对氧的需求量更强烈(Paneetal.,2003;2004).对于不同金属的版本,共同的现象是当pH>8.0时,都可能会高估LC50值;DOC的取值在很大程度上了模型预测的准确性;大多数模型都没有在自然水体环境中进行验证.这都需要在以后的研究中进一步完善,但是模型在实验室条件下都取得了较好的效果.
大约是版本IIa和IIIa的100倍,由
+
于DOC和KCu
(HCO3)
取值不同导致了各版本LA50
(半致死累积量)值的不同.版本IIa和IIIa认为pH对铜毒性的影响不仅仅是由于与自由铜离子竞争结合位点的H+浓度的变化而引起的,铜的无机络合态(Cu(OH)+、CuCO3)也是有毒的.对于银的3个BLM版本,同另外2个相比,Ia考虑了AgCl的毒性,并求出了KBL-AgCl值.版本IIa认为银对鱼的毒性是由于抑制了Na+-K+-ATPase(腺苷三磷酸酶)的活性,而Na+-K+-ATPase是鱼吸收Na+、Cl-的能量来源(McGeeretal.,1998).因此IIa通过银对Na+-K+-ATPase的抑制率来求KBL-AgCl,当抑制率在85%时相当于96hLC50.Ia和IIIa尽管认识到银的生物毒性是由于阻碍了Na+、Cl-的吸收,但LC50是通过鱼鳃上银的累积量来求的(Woodetal.,
0
2.2陆地生态系统
生物配体模型在陆地生态系统中的应用要落
后于水生生态系统.目前已取得了一些金属对陆地植物(大麦、莴苣等)、动物(蚯蚓)、发光细菌的毒性数据,但是这些数据还非常有限,陆地生态系统生物配体模型(t-BLM)的研究还刚刚起步.下面概述一下t-BLM的研究进展及研究思路和方法(主要以植物为例).
过去生理学家发现并证实了植物根同鱼鳃一样,是最初的金属作用位点.在此基础上,又有人研究了竞争和络合机理在控制金属生物有效性中的作用.如Voigt等(2006)研究了金属和根的络合以及根吸附在控制金属生物有效性方面的作用,发现根具有和金属不连续的作用位点,而且Ca2+具有和金属竞争结合位点的作用.Zhao等(2006)通过在18种性质不同的土壤中添加不同浓度铜的试验研究了铜对大麦的毒性,发现大麦根的伸长和土壤溶液中自由铜离子活度、pH值(H+具有明显抑制铜毒性的作用)有关,和应用到水生生态系统中的生物配体模型理论一致.这些研究结果暗示了应用到水生生态系统中的生物配体模型经过修正并结合土壤理化性质后可以应用到陆地生态系统中.
虽然水生生态系统中的生物配体模型为t-1996;Morganetal.,1997).
由上可见,同一金属的不同BLM版本在一定程度上存在着差异,虽然后来的版本基本上是在前边版本的基础上演变而来,但并不一定版本越高预测效果就越好,只是它们适用的范围和条件不同,脱离这些条件,它们将达不到预测的效果.对于Cu-BLM,版本Ia在高K+浓度时会高估LC50值,
而且也会高估暴露在浅水湖泊中Fathead
minnow的LC50值(Ericksonetal.,1996);而版本IIa在pH>8时会高估EC50值(DeSchamphelaereetal.,2002a).Zn-BLM版本中Ia在预言Rainbowtrout和FatheadminnowEC50以及Daphniamagna48hLC50时取得了较好的效果,但是当水的pH>8时,模型预言的LC50值偏高,这可能和锌在水中的溶解性有关(Niyogietal.,2004).另外Ia在建立过程中没有对DOC水平进行,可能会导致不正确的DOC-Zn估计和自由锌活度的计算.目前的研究证明Ia对鱼的预测效果要比对Daphnia
magna的预测效果好,而IIa正好相反(Niyogietal.,2004).3个急性Ag-BLM版本在实验室里都进行了成功的预测,但是它们均不能用于银对(NiyogiAmphipod、Gammaruspulex的急性毒性预测
BLM提供了研究基础,但土壤或沉积物中的重金属生物毒害影响因子远比水体复杂,如土壤理化性质(土壤有机质、阳离子交换量、pH值等)和生物本身(如植物的根系分泌物)均影响土壤金属的生物毒害(McLaughlin,2000).应用生物配体模型预测土壤中金属的生物毒性首先需要从土壤中获得模型所需要输入的数据,目前的方法主要是提取用于分析的土壤溶液.用不同方法提取土壤孔隙水产生的效果不同.MacDonald等(2004a)提出
etal.,2004).对于Ni-BLM版本来说,在预测对Fatheadminnow的毒性时可能会高估LA50值,需要根据Fatheadminnow年龄和体形的大小调整LA50值(Groselletal.,2002).Hoang等(2004)也发现第3期王学东等:环境中金属生物有效性的预测模型———生物配体模型研究进展
199
淋洗柱淋洗然后用渗透来提取土壤溶液的方法,较好地模拟了田间土壤的化学性质.还有一些如离心、压力、抽气等方法提取的土壤溶液溶解性离子的浓度要大于渗透方法测得的离子浓度.
对金属累积和毒性终点的确定(以植物为例)首先需找到金属的受体位点.最初发展的生物配体模型把鱼鳃作为重金属的结合和毒性位点,而植物根质外体(包括细胞壁和Donnan自由空间)被认为是重金属在植物中的主要结合位点(Antunesetal.,2006).但对于植物学家来说,质外体结合不是真正意义上的植物对重金属的累积,因为质外体包括许多无生命物质如细胞壁和木质成分.在低的金属暴露浓度下,金属和这些无生命物质结合可能对植物没有任何毒性.在高的金属暴露浓度下,细胞壁发生改变有可能导致根际毒性(Marschner,1995).因此,发展土壤和植物的生物配体模型,需要考虑这些情况.Antunes等(2006)提出了通过比较不同植物部分、不同反应条件下的logKMe-plant
root
的能力减弱,金属毒性也会增强.2种说法都可以解释同一种结果,究竟那种说法正确还是2种情况都存在还需进一步研究.
生物配体模型在计算金属和DOC反应时通常假设DOC由一定比例的FA和HA或仅由FA组成(Niyogietal.,2004),而自然环境中DOC的组成和结构远比模型假设的要复杂,而且DOC是影响不同形态金属尤其是自由金属离子浓度的重要因素.研究发现,增加不同的络合物将会导致
Ag-DOC络合特征的不同(VanGenderenetal.,2003);在预测镍对生物的毒性时,镍和DOC络合以后其浓度可能会降至镍的毒性范围以下(0.1 ̄
1.0mg・L-1)(Niyogietal.,2004).由此可见,将模型应用到和模型中DOC假设不同的自然环境中就可能得出和实际情况相悖的结论.因此对DOC结构特征的研究也是未来需要解决的问题.
生物配体模型假设金属(通常为自由态金属离子)及其生物配体络合物的浓度代表了金属的毒性,但研究发现一些憎水性金属无机络合物如
(金属和植物质外体结合
shoot
的条件稳定常数)和logKMe-plant(金属和植物茎结HgCl2和CH3HgCl可以通过被动吸收的形式直接跨过生物膜而对生物产生毒害(Masonetal.,
合的条件稳定常数),以确定测试植物的种类及测定方法.但是其方法还有待标准化.总之,土壤中重金属生物毒害预测远比水体复杂,一些机理和方法还有待于进一步研究.
1996).在发展Ag-BLM时遇到的一个困难是不清楚Cl-在银的生物毒性过程中到底扮演何种角色,因为Cl-有减轻银毒性的作用,但并不阻碍鱼鳃上银的累积.观察还发现Cl-对Daphniamagna仅有轻微的保护作用,而对Fatheadminnow不存在保护作用(Niyogietal.,2004).在预测土壤中重金属对植物的有效性时这些情况可能会更明显,因为根际具有强烈的微生物活性和根系分泌物的累积,特别是这些根系分泌物可以和重金属形成络合物而被植物吸收,从而增加了金属的生物有效性(MacDonaldetal.,2004b).在重金属浓度较低时,金属络合物对金属毒性的贡献更需要考虑(Lehmannetal.,2000).
生物配体模型在实现的过程中还假设金属向生物配体移动时没有速率.实际上,当生物吸收较快时金属向生物配体移动就要受速率的,不过这种情况相对较少,但在底泥和土壤中这种情况就需要重视(Buffle,1988).
在金属对生物产生毒性的过程中生物配体模型假设配体的特征和数量没有改变.Heijerick等(2003)在预测锌对Daphniamagna毒性时提出了一个表面反应模型,考虑了这些变化.在较大
3生物配体模型的局限性和未来的研究方向生物配体模型的应用尤其是在水生生态系统
方面的应用取得了一定的进展,建立了不同金属的预测模型版本,但模型的应用范围还具有一定的局限性,这也是模型在未来的研究中需要解决的问题.
从模型在水生生态系统中的应用可知,当水体的pH>8时,大部分金属-BLM的预测效果都不好,即使pH不超出这个值,模型也存在着不能解释的情况.因为pH不但影响金属和有机、无机、生物配体的所有反应,而且会影响生物配体的某些特征,忽略了任何一方面,都可能得出不正确的预测结果.在预测铜的生物毒性时,除Ia外,其他版本都计算了KBL-Cu(加,
OH)
+
,有研究者认为Cu(OH)
+
对生物也有毒性,当pH增高,Cu(OH)+的浓度增
从而增加了金属对生物的毒性(De
Schamphelaereetal.,2002b;2005).从竞争原理来讲,高pH值下H+浓度减少,和金属竞争生物配体
200
生态毒理学报第1卷
硬度、(Ia、pH、DOC范围下比目前的Zn-BLMIIa)预(Ia、测效果要好,这可能是Zn-BLMIIa)没有考虑到
通讯作者简介:马义兵(1957—),男,研究员,博士生导师,中国农业科学院引进的杰出人才.现任中国农业科学院农业资源与农业区划研究所亚洲环境重金属研究中心主任.主要从事土壤重金属环境化学和生态毒理方面的研究工作.
pH的变化对生物配体特征和数量的影响.
痕量金属的吸收要受生物学的控制,如锌浓度增加或减少1000倍时,尽管Fatheadminnow鳃上结合位点增加,
但慢性半致死剂量减少
(Slaveykovaetal.,2005).这个结果说明生物配体模型对急性毒性的预测比对慢性毒性的预测要好,急性生物配体模型不能作为慢性模型的基础.在预测铜对Daphniamagna的毒性时也出现了这种情况,慢性版本Ic没有发现Ca2+、Mg2+竞争对铜毒性的影响.另外,自然界的重金属并不是单独存在的,模型没有考虑其他重金属的存在对这种金属毒性的影响.可见,要想达到更加满意的预测效果,生物配体模型还有很长的路要走.
在陆地生态系统中还有更多的情况需要考虑.目前采用提取土壤孔隙水的方法实际上简化了土壤对金属生物毒性的影响,这种方法忽视或简化了金属的老化过程、生物影响下金属的有效化过程、金属释放的动力学过程等对金属生物有效性的影响,尤其在研究金属慢性生物毒性时这些情况需更加重视.另外在预测金属毒性的指标(金属和配体的络合常数、配体的络合容量以及
References
AlsopDH,WoodCM.2000.Akineticanalysisofzinc
accumulationinthegillsofjuvenileRainbowtrout:Theeffectsofzincacclimationandimplicationsforbioticligandmodeling[J].EnvironToxicolChem,19:1911-1918
AntunesPMC,BerkelaarEB,BoyleD,HaleBA,HendershotW,VoigtA.2006.Thebioticligandmodelforplantsandmetals:Technicalchallengesforfieldapplications[J].EnvironToxicolChem,25:875-882
BinghamFT,StrongJE,SpositoG.1983.InfluenceofchloridesalinityoncadmiumuptakebySwisschard[J].SoilSci,35:160-165
BuffleJ.1988.Complexationreactioninaquaticsystems:Ananalyticalapproach[M].Chichester:EillsHorwood,692BuryNR,ShawJ,GloverC,HogstrandC.2002.Derivationofatoxicity-basedmodeltopredicthowwaterchemistryinfluences
[J].CompBiochemPhysiolC,133:silvertoxicitytoinvertebrates
259-270
CampbellPGC.1995.Interactionsbetweentracemetalsaquaticorganisms:Acritiqueofthefree-ionactivitymodel[In:TessierA,TurnerDR(Eds.).Metalspeciationbioavailabilityinaquaticsystems.NewYork:JohnWileySons,45-102
andM].andand
EC50等)时还没有较好的方法,目前较多采用的数学方法不能完全从机理方面解释金属的毒性机理.在金属对植物毒性的研究中,虽然对植物根茎重金属分配等研究为金属累积和毒性终点确定提供了思路,但由于根和茎生理过程的不同,金属在根茎中的分配还要受以下几个因素的影响:矿质营养、环境影响、金属之间的竞争作用、植物的健康状况等(Antunesetal.,2006),忽略了任何一个因素,预测结果都会产生偏差.
总之,生物配体模型是一个机理性的模型,它包含了生物有效性的概念,比较全面地考虑了影响金属生物有效性的因素,和建立在硬度基础上的模型相比,生物配体模型具有更好的预测能力,是一个能够替代生物毒性试验预测金属毒性的有用工具.尽管模型现在存在着局限性,但随着对痕量金属吸收过程中化学、生理学、生物学过程的深入理解,这些问题将会逐步得到解决,未来生物配体模型既面临挑战又有着广阔的应用和发展前景.
DeSchamphelaereKAC,HeijerickDG,JanssenCR.2002a.
Refinementandfield:ValidationofabioticligandmodelpredictingacutecoppertoxicitytoDaphniamagna[J].CompBiochemPhysiolC,133:243-258
DeSchamphelaereKAC,JanssenCR.2002b.AbioticligandmodelpredictingacutecoppertoxicityforDaphniamagna:Theeffectsofcalcium,magnesium,sodium,potassium,andpH[J].EnvironSciTechnol,36:48-54
DeSchamphelaereKAC,JanssenCR.2004a.DevelopmentandfieldvalidationofabioticligandmodelpredictingchroniccoppertoxicitytoDaphniamagna[J].EnvironToxicolChem,23:1365-1375
DeSchamphelaereKAC,JanssenCR.2004b.BioavailabilityandchronictoxicityofzinctojuvenileRainbowtrout(Oncorhynchusmykiss):Comparisonwithotherfishspeciesanddevelopmentofabioticligandmodel[J].EnvironSciTechnol,38:6201-6209DeforarepH
SchamphelaereKAC.2005.TowardabioticligandmodelfreshwaterGreenAlgae:Surface-boundandinternalcopperbetterpredictorsoftoxicitythanfreeCu2+ionactivitywhenisvaried[J].EnvironSciTechnol,39:2067-2072
DiToroDM,AllenHE,BergmanHL,MeyerJS,PaquinP第3期王学东等:环境中金属生物有效性的预测模型———生物配体模型研究进展
201
R,SantoreRC.2001.Bioticligandmodeloftheacutetoxicity
ofmetals.1Technicalbasis[J].EnvironToxicolChem,20:2383-2396
EricksonRJ,BenoitDA,MattsonVR,NelsonJr,HenryP,LeonardEN.1996.TheeffectsofwaterchemistryonthetoxicityofcoppertoFatheadminnows[J].EnvironToxicolChem,15:181-193
GrosellM,NielsenC,BianchiniA.2002.Sodiumturnoverratedeterminessensitivitytoacutecopperandsilverexposureinfreshwateranimals[J].CompBiochemPhysiolC,133:287-303HeijerickDG,DeSchamphelaereKAC,JanssenCR.2002.PredictingacutezinctoxicityforDaphniamagnaasafunctionofkeywaterchemistrycharacteristics:Developmentandvalidationof
(6):1309-abioticligandmodel[J].EnvironToxicolChem,21
1315
HeijerickDG,JanssenCR,DeCoenWM.2003.Thecombinedeffectsofhardness,pH,anddissolvedorganiccarbononthechronictoxicityofZntoD.magna:Developmentofasurfaceresponsemodel[J].ArchEnvironContamToxicol,44:210-217
HoangTC,TomassoJR,KlaineSJ.2004.InfluenceofwaterqualityandageonnickeltoxicitytoFatheadminnows(Pimephalespromelas)[J].EnvironToxicolChem,23:86-92HuangSB.2003.Measurementandpredictioninthetoxic/bioavailabileconcentrationofcopperinnaturalwater[D],Beijing:ResearchCenterforEco-EnvironmentalSciences,ChineseAcademyofSciences(inChinese)
HunnJB.1985.Roleofcalciumingillfunctioninfreshwaterfishes[J].CompBiochemPhysiol,82:543-547
JanesN,PlayleRC.1995.ModelingsilverbindingtogillsofRainbowtrout(Oncorhynchusmykiss)[J].EnvironToxicolChem,14:1847-1858
LehmannH,StelzerR,HolzamerS,KunzU.2000.Analyticalelectronmicroscopicalinvestigationsontheapoplasticpathwaysoflanthanumtransportinbarleyroots[J].Planta,211:816-822MacDonaldA,SilkL,SchwartzM,PlayleRC.2002.Alead-gillbindingmodeltopredictacuteleadtoxicitytoRainbowtrout(Oncorhynchusmykiss)[J].CompBiochemPhysiolC,133:227-242
MacDonaldJD,BelangerN,HendershotWH.2004a.Columnleachingusingdrysoiltoestimatesolid-solutionpartitioningobservedinzero-tensionlysimeters.1Methoddevelopment[J].SoilSedimentContam,13:361-374
MacDonaldJD,BelangerN,HendershotWH.2004b.Columnleachingusingdrysoiltoestimatesolid-solutionpartitioningobservedinzero-tensionlysimeters.2Tracemetals[J].SoilSedimentContam,13:375-390
MarkichSJ,BrownPL,JeffreeRA,LimRP.2003.TheeffectsofpHanddissolvedorganiccarbononthetoxicityofcadmiumandcoppertoafreshwaterBivalve:FurthersupportfortheextendedFreeIonActivityModel[J].ArchEnvironContam
(4):479-491Toxico,l45MarschnerH.1995.Mineralnutritionofhigherplants[M].
London:AcademicPress,889
MartellAE,SmithRM,MotekaitisRJ.1997.Criticalstabilityconstantsofmetalcomplexesdatabase,version4.0,NISTstandardreferencedatabase46[Z].GaithersburgMD:NationalInstituteofStandardsandTechnology
MasonRP,ReinfelderJR,MorelFMM.1996.Theuptake,toxicity,andtrophictransferofmercuryinacoastaldiatom[J].EnvironSciTechnol,30:1835-1845
McGeerJC,PlayleRC,WoodCM,GalvezF.2000.AphysiologicallybasedbioticligandmodelforpredictingtheacutetoxicityofwaterbornesilvertoRainbowtroutinfreshwaters[J].EnvironSciTechnol,34:4199-4207
McGeerJC,WoodCM.1998.ProtectiveeffectsofwaterCl-onphysiologicalresponsestowaterbornesilverinRainbowtrout[J].CanJFishAquatSci,55:2447-2454
McLaughlinMJ.2000.Bioavailabilityofmetalstoterrestrialplants[M].InAllenHE(Ed.).Bioavailabilityofmetalsinterrestrialecosystems.InfluenceofpartitioningforbioavailabilitytoInvertebrates.Pensacola,FL:SETACPress,39-67
McLaughlinMJ,TillerKG,SmartMK.1997.Speciationofcadmiuminsoilsolutionsofsaline-sodicsoilsandrelationshipwithcadmiumconcentrationsinpotatotubers(SolanumtuberosumL.)[J].AustJSoilRes,35:183-198
MeyerJS,SantoreRC,BobbittJP,DeBreyLD,BoeseCJ,PaquinPR,AllenHE,BergmanHL,DiToroDM.1999.Bindingofnickelandcoppertofishgillspredictstoxicitywhenwaterhardnessvaries,butfree-ionactivitydoesnot[J].EnvironSciTechnol,33:913-916
MorelFM.1983.Principlesofaquaticchemistry[M].NewYork:JohnWiley&Sons,300-309
MorganIJ,HenryRP,WoodCM.1997.ThemechanismofacutesilvertoxicityinfreshwaterRainbowtrout(Oncorhynchusmykiss)isinhibitionofgillNa+andCl-transport[J].AquatToxicol,38:145-163
NiyogiS,WoodCM.2004.Bioticligandmodel,aflexibletoolfordevelopingsite-specificwaterqualityguidelinesformetals[J].EnvironSciTechol,38:6177-6192
PagenkopfGK.1983.Gillsurfaceinteractionmodelfortrace-metaltoxicitytofishes:Roleofcomplexation,pHandwaterhardness[J].EnvironSciTechol,17:342-347
PaneEF,HaqueA,GossGG,WoodCM.2004.Thephysiologicalconsequencesofexposuretochronic,sublethalwaterbornenickelinRainbowtrout(Oncorhynchusmykiss):Exercisevsrestingphysiology[J].JExpBiol,207:1249-1261PaneEF,RichardsJG,WoodCM.2003.AcutewaterbornenickeltoxicityintheRainbowtrout(Oncorhynchusmykiss)occursbyarespiratoryratherthananionoregulatorymechanism[J].AquatToxicol,63:65-82
PaquinPR,DiToroDM,SantoreRC,TrivediD,WuKB.1999.AbioticligandmodeloftheacutetoxicityofmetalsIIIApplicationtofishandDaphniamagnaexposuretosilver.
202
生态毒理学报第1卷
Integratedapproachtoassessingthebioavailabilityandtoxicityofmetalsinsurfacewatersandsediments,Section3,asubmissiontotheEPAScienceAdvisoryBoard,OfficeofWater[R].USEPA-822-E-99-001.OfficeofResearchandDevelopment,WashingtonDC,USA,59-102
PaquinPR,GorsuchJW,ApteS,BatleyGE,BowlesKC,CampbellPG,DelosCG,DiToroDM,DwyerRL,GalvezF,GensemerRW,GossGG,HostrandC,JanssenCR,McGeerJC,NaddyRB,PlayleRC,SantoreRC,SchneiderU,StubblefieldWA,WoodCM,WuKB.2002.Thebioticligandmodel:Ahistoricaloverview[J].CompBiochemPhysiolC,133:3-35
PaquinPR,RobertC,SantoreRC,WuaKB,KavvadasCD,DiToroDM.2000.Thebioticligandmodel:Amodeloftheacutetoxicityofmetalstoaquaticlife[J].EnvironmentalScience&Policy3:S175-S182
PlayleRC,DixonDG,BumisonK.1993a.Copperandcadmiumbindingtofishgills:Modificationbydissolvedorganiccarbonandsyntheticligands[J].CanJFishAquatSci,50:2667-2677
PlayleRC,DixonDG,BurnisonK.1993b.Copperandcadmiumbindingtofishgills:Estimatesofmetal-gillstabilityconstantsandmodellingofmetalaccumulation[J].CanJFishAquatSci,50:2678-2687
PlayleRC,GensemerRW,DixonDG.1992.CopperaccumulationongillsofFatheadminnows:Influenceofwaterhardness,complexationandpHofthegillmicro-environment[J].EnvironToxicolChem,11:381-391
PoldoskiJE.1979.Cadmiumbio-accumulationassaystheirrelationshiptovariousionicequilibriainlakesuperiorwater[J].EnvironSciTechol,13:701
RichardsJG,PlayleRC.1998.Cobaltbindingtogillsof
(Oncorhynchusmykiss):Anequilibriummodel[J].Rainbowtrout
CompBiochemPhysiolC,119:185-197
SaekiK,KunitoT,OyaizuH,MatsumotoS.2002.Relationshipsbetweenbacterialtolerancelevelsandformsofcopperandzincinsoils[J].JEnvironQual,31:1570-1575
SantoreB.2006.Informationofsymposiumonbioticligandmodel&environmental(ecological)riskassessment[C],Beijing(inChinese)
SantoreRC,DiToroDM,PaquinPR,AllenHE,MeyerJS.2001.Bioticligandmodeloftheacutetoxicityofmetals.2.ApplicationtoacutecoppertoxicityinfreshwaterfishandDaphnia[J].EnvironToxicolChem,20:2397-2402
SantoreRC,DriscollCT.1995.Inchemicalequilibriumandreactionmodels[M].In:LoeppertR,SchwabAP,GoldbergS(Eds.).SoilsciencesocietyofAmericaspecialpublication42.
MadisonWI:AmericanSocietyofAgronomy,357-375
SantoreRC,MathewR,PaquinPR,DiToroD.2002.ApplicationofthebioticligandmodeltopredictingzinctoxicitytoRainbowtrout,FatheadminnowandDaphniamagna[J].CompBiochemPhysiolC,133:271-285
SauvéS,CookN,HendershotWH,McBrideMB.1996.Linkingplanttissueconcentrationsandsoilcopperpoolsinurbancontaminatedsoils[J].EnvironPollut,94:153-157
SchwartzML,PlayleRC.2001.Addingmagnesiumtothesilver-gillbindingmodelforRainbowtrout(Oncorhynchusmykiss)[J].EnvironToxicolChem,20:467-472
SlaveykovaVI,WilkinsonKJ.2005.Predictingthebioavailabilityofmetalcomplexes:Criticalreviewofthebioticligandmodel[J].EnvironChem,2:9-24
SmoldersE,LambregtsRM,McLaughlinMJ,TillerKG.1998.EffectofsoilsolutionchlorideoncadmiumavailabilitytoSwisschard[J].JEnvironQual,27:426-431
TippingE.1994.WHAM-Achemicalequilibriummodelandcomputercodeforwaters,sedimentsandsoilsincorporatingadiscretesite/electrostaticmodelofion-bindingbyhumicsubstances[J].ComputGeosci,20:973-1023
VanGenderenEJ,RyanAC,TomassoJR,KlaineSJ.2003.InfluenceofdissolvedorganicmattersourceonsilvertoxicitytoPimephalespromelas[J].EnvironToxicolChem,22:2746-2751VoigtA,HendershotWH,SunaharaGI.2006.Effect-plant,soil,invertebrate,microflorarhizotoxicityofcadmiumandcopperinsoilextracts[J].EnvironToxicolChem,25:692-701
WaterEnvironmentResearchFoundation.2002.DevelopmentofaBioticLigandModelfornickel[A].HydroqualProjectWERF0040.Hydroqual,MahwahNJ
WoodCM,HogstrandC,GalvezF,MungerRS.1996.ThephysiologyofwaterbornesilvertoxicityinfreshwaterRainbowtrout(Oncorhynchusmykiss)1.TheeffectsofionicAg+[J].AquatToxicol,35:93-109
ZhaoFJ,RooneyCP,ZhangH,McGrathSP.2006.Comparisonofsoilsolutionspeciationanddiffusivegradientsinthin-filmsmeasurementasanindicatorofcopperbioavailabilitytoplants[J].EnvironToxicolChem,25:733-742
中文参考文献
SantoreB.2006.生物配体模型和环境(生态)风险评价高级研讨会会议资料[C].北京
黄圣彪.2003.水环境中铜形态与其生物有效性/毒性关系及其预测模型研究[D].北京:中国科学院生态环境研究中心
◆
因篇幅问题不能全部显示,请点此查看更多更全内容
Copyright © 2019- sarr.cn 版权所有 赣ICP备2024042794号-1
违法及侵权请联系:TEL:199 1889 7713 E-MAIL:2724546146@qq.com
本站由北京市万商天勤律师事务所王兴未律师提供法律服务