CN109583662B - 大气污染二元混合物健康风险评价方法 - Google Patents
大气污染二元混合物健康风险评价方法 Download PDFInfo
- Publication number
- CN109583662B CN109583662B CN201811500291.5A CN201811500291A CN109583662B CN 109583662 B CN109583662 B CN 109583662B CN 201811500291 A CN201811500291 A CN 201811500291A CN 109583662 B CN109583662 B CN 109583662B
- Authority
- CN
- China
- Prior art keywords
- reference body
- ref
- concentration
- relative
- binary mixture
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 239000000203 mixture Substances 0.000 title claims abstract description 110
- 230000036541 health Effects 0.000 title claims abstract description 55
- 238000011156 evaluation Methods 0.000 title claims abstract description 24
- 231100000419 toxicity Toxicity 0.000 claims abstract description 66
- 230000001988 toxicity Effects 0.000 claims abstract description 66
- 239000000178 monomer Substances 0.000 claims abstract description 58
- 239000003344 environmental pollutant Substances 0.000 claims abstract description 49
- 231100000719 pollutant Toxicity 0.000 claims abstract description 49
- 238000000338 in vitro Methods 0.000 claims abstract description 43
- 230000000694 effects Effects 0.000 claims abstract description 42
- 238000000034 method Methods 0.000 claims abstract description 33
- 238000012360 testing method Methods 0.000 claims abstract description 30
- 230000009471 action Effects 0.000 claims abstract description 11
- 238000012544 monitoring process Methods 0.000 claims abstract description 4
- 108090000623 proteins and genes Proteins 0.000 claims description 32
- 230000037361 pathway Effects 0.000 claims description 28
- 230000014509 gene expression Effects 0.000 claims description 14
- 230000004044 response Effects 0.000 claims description 13
- 231100000331 toxic Toxicity 0.000 claims description 12
- 230000001186 cumulative effect Effects 0.000 claims description 11
- 238000012502 risk assessment Methods 0.000 claims description 11
- 231100000673 dose–response relationship Toxicity 0.000 claims description 10
- 239000000356 contaminant Substances 0.000 claims description 8
- 230000002588 toxic effect Effects 0.000 claims description 8
- 230000008859 change Effects 0.000 claims description 7
- 238000004088 simulation Methods 0.000 claims description 7
- 230000004913 activation Effects 0.000 claims description 6
- 238000004364 calculation method Methods 0.000 claims description 6
- 230000036387 respiratory rate Effects 0.000 claims description 5
- 238000012614 Monte-Carlo sampling Methods 0.000 claims description 4
- 230000000996 additive effect Effects 0.000 claims description 4
- 238000004113 cell culture Methods 0.000 claims description 4
- 238000012512 characterization method Methods 0.000 claims description 4
- 238000002474 experimental method Methods 0.000 claims description 4
- 239000002245 particle Substances 0.000 claims description 4
- 230000008021 deposition Effects 0.000 claims description 3
- 230000003042 antagnostic effect Effects 0.000 claims description 2
- 238000012937 correction Methods 0.000 claims description 2
- 230000000977 initiatory effect Effects 0.000 claims description 2
- 230000002195 synergetic effect Effects 0.000 claims description 2
- 229910002056 binary alloy Inorganic materials 0.000 claims 1
- 230000007613 environmental effect Effects 0.000 abstract description 4
- 239000012298 atmosphere Substances 0.000 abstract description 2
- 230000008878 coupling Effects 0.000 abstract description 2
- 238000010168 coupling process Methods 0.000 abstract description 2
- 238000005859 coupling reaction Methods 0.000 abstract description 2
- 239000000463 material Substances 0.000 abstract description 2
- 238000000342 Monte Carlo simulation Methods 0.000 abstract 1
- 231100000048 toxicity data Toxicity 0.000 abstract 1
- DXBHBZVCASKNBY-UHFFFAOYSA-N 1,2-Benz(a)anthracene Chemical compound C1=CC=C2C3=CC4=CC=CC=C4C=C3C=CC2=C1 DXBHBZVCASKNBY-UHFFFAOYSA-N 0.000 description 68
- FMMWHPNWAFZXNH-UHFFFAOYSA-N Benz[a]pyrene Chemical compound C1=C2C3=CC=CC=C3C=C(C=C3)C2=C2C3=CC=CC2=C1 FMMWHPNWAFZXNH-UHFFFAOYSA-N 0.000 description 46
- 125000005575 polycyclic aromatic hydrocarbon group Chemical group 0.000 description 15
- 238000009826 distribution Methods 0.000 description 10
- 108091006146 Channels Proteins 0.000 description 9
- 239000000126 substance Substances 0.000 description 7
- 238000003915 air pollution Methods 0.000 description 6
- 238000011160 research Methods 0.000 description 5
- BBEAQIROQSPTKN-UHFFFAOYSA-N pyrene Chemical compound C1=CC=C2C=CC3=CC=CC4=CC=C1C2=C43 BBEAQIROQSPTKN-UHFFFAOYSA-N 0.000 description 4
- 239000013558 reference substance Substances 0.000 description 4
- 231100000820 toxicity test Toxicity 0.000 description 4
- 238000004166 bioassay Methods 0.000 description 3
- 230000000711 cancerogenic effect Effects 0.000 description 3
- 150000001875 compounds Chemical class 0.000 description 3
- 230000005183 environmental health Effects 0.000 description 3
- 238000003912 environmental pollution Methods 0.000 description 3
- 230000008821 health effect Effects 0.000 description 3
- 230000007246 mechanism Effects 0.000 description 3
- MWUXSHHQAYIFBG-UHFFFAOYSA-N nitrogen oxide Inorganic materials O=[N] MWUXSHHQAYIFBG-UHFFFAOYSA-N 0.000 description 3
- 239000013618 particulate matter Substances 0.000 description 3
- 230000029058 respiratory gaseous exchange Effects 0.000 description 3
- 108010057466 NF-kappa B Proteins 0.000 description 2
- 102000003945 NF-kappa B Human genes 0.000 description 2
- 239000000654 additive Substances 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 230000008827 biological function Effects 0.000 description 2
- 231100000315 carcinogenic Toxicity 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- GVEPBJHOBDJJJI-UHFFFAOYSA-N fluoranthrene Natural products C1=CC(C2=CC=CC=C22)=C3C2=CC=CC3=C1 GVEPBJHOBDJJJI-UHFFFAOYSA-N 0.000 description 2
- 231100000003 human carcinogen Toxicity 0.000 description 2
- 210000004072 lung Anatomy 0.000 description 2
- 150000003071 polychlorinated biphenyls Chemical class 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 231100000027 toxicology Toxicity 0.000 description 2
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 1
- 108020005124 DNA Adducts Proteins 0.000 description 1
- 241000282412 Homo Species 0.000 description 1
- 206010028980 Neoplasm Diseases 0.000 description 1
- 241000700159 Rattus Species 0.000 description 1
- QAOWNCQODCNURD-UHFFFAOYSA-L Sulfate Chemical compound [O-]S([O-])(=O)=O QAOWNCQODCNURD-UHFFFAOYSA-L 0.000 description 1
- 230000002411 adverse Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000008485 antagonism Effects 0.000 description 1
- 239000008277 atmospheric particulate matter Substances 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 125000005605 benzo group Chemical group 0.000 description 1
- 201000011510 cancer Diseases 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 239000003795 chemical substances by application Substances 0.000 description 1
- 238000007635 classification algorithm Methods 0.000 description 1
- 238000011109 contamination Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 201000010099 disease Diseases 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 231100000305 environmental exposure data Toxicity 0.000 description 1
- 239000003256 environmental substance Substances 0.000 description 1
- 231100000613 environmental toxicology Toxicity 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 239000010419 fine particle Substances 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 239000007789 gas Substances 0.000 description 1
- 231100000206 health hazard Toxicity 0.000 description 1
- 238000010438 heat treatment Methods 0.000 description 1
- 229910001385 heavy metal Inorganic materials 0.000 description 1
- 210000005260 human cell Anatomy 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000000099 in vitro assay Methods 0.000 description 1
- 238000001727 in vivo Methods 0.000 description 1
- 238000010801 machine learning Methods 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000003305 oral gavage Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000001105 regulatory effect Effects 0.000 description 1
- 241000894007 species Species 0.000 description 1
- 231100000440 toxicity profile Toxicity 0.000 description 1
- 238000009827 uniform distribution Methods 0.000 description 1
- 108700026215 vpr Genes Proteins 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H50/00—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
- G16H50/30—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for calculating health indices; for individual health risk assessment
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A50/00—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE in human health protection, e.g. against extreme weather
- Y02A50/20—Air quality improvement or preservation, e.g. vehicle emission control or emission reduction by using catalytic converters
Landscapes
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Strategic Management (AREA)
- Human Resources & Organizations (AREA)
- Economics (AREA)
- Medical Informatics (AREA)
- Health & Medical Sciences (AREA)
- Public Health (AREA)
- Development Economics (AREA)
- Entrepreneurship & Innovation (AREA)
- Primary Health Care (AREA)
- Epidemiology (AREA)
- Pathology (AREA)
- Databases & Information Systems (AREA)
- Game Theory and Decision Science (AREA)
- Data Mining & Analysis (AREA)
- Biomedical Technology (AREA)
- General Health & Medical Sciences (AREA)
- Marketing (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Tourism & Hospitality (AREA)
- Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Measuring Or Testing Involving Enzymes Or Micro-Organisms (AREA)
Abstract
一种大气污染二元混合物健康风险评价新方法,包括以下步骤:从文献或实际环境监测中获得大气中污染物浓度,选取近似浓度的污染物作为待评价的二元混合组分,并将实际环境暴露浓度转化为体外等价浓度;从公开数据库或文献中获得每种污染物的转录组数据和特定通路的体外生物测试信息,计算其混合效应数值、校正相对效能因子和联合毒性作用模式;对污染物体外等价浓度进行蒙特卡洛模拟并耦合校正相对效能因子,评价混合污染物的健康风险。本发明建立在污染物单体环境暴露浓度和生物测试数据的基础上,不依赖混合物毒性数据,便可评价其健康风险,节约人力物力成本,在实际环境混合暴露健康风险评价和环境安全性评估等领域具有广阔的应用前景。
Description
技术领域
本发明涉及污染混合暴露的健康风险评价技术领域,更具体地涉及一种大气污染二元混合物健康风险评价方法。
背景技术
环境污染问题从来都不仅仅是单一化学品的污染问题。真实环境中种类繁杂的污染物往往以混合形式存在,并暴露于人群产生联合毒性效应。由于混合污染联合作用产生复杂多变的健康效应,计算评价联合毒性及其健康风险一直以来都是科学界关注的难点和焦点。早在2013年美国国立环境卫生科学研究所(National Institute of EnvironmentalHealth Sciences,NIEHS)就指出应以阐明环境混合物的健康效应作为其未来优先研究重点之一。大气细颗粒物(Fine Particulate Matters,PM2.5)是天然混合体系。以大气颗粒物(Particulate Matters,PM)污染为例,灰霾通常涉及硫酸盐、氮氧化物、多环芳烃(Polycyclic Aromatic Hydrocarbons,PAHs)、多氯联苯(Polychlorinated Biphenyls,PCBs)和有毒重金属等各种化学成分的污染问题。即使主要有机污染组分多环芳烃,在大气中也是以多种同类物混合的方式出现,其中包括众所周知的人类致癌物苯并[a]芘(benzo[a]pyrene,BaP)。国际癌症研究机构(International Agency for Research on Cancer,IARC)在2013年将室外空气污染和大气细颗粒物归类为对人类致癌(第1组)(http://www.iarc.fr/en/media-center/iarcnews/pdf/pr221_E.pdf)。然而,PM2.5暴露对人类健康的影响不仅因特定地区污染组成不同而在不同地区之间存在显着差异,同一地区不同来源的大气污染物亦呈现迥异的健康影响。2005年Kan等人指出中国空气污染所产生的健康影响与发达国家不同。Laden等人发现机动车、燃煤来源的PM2.5每增高10μg/m3,人群每日死亡率分别增高3.4%、1.1%,机动车来源的细颗粒物健康危害更高。可见,准确预测大气污染混合物的联合毒性是科学评估其累积健康风险的前提条件。
遗憾的是,目前关于混合物联合毒性作用非实验评估的进展远不如单个污染物毒性预测方法的研究。例如已有研究应用机器学习分类算法对毒性数据库中对单一化学品特定毒性进行分类和预测。无论等效线图还是析因分析法,现有的联合毒性预测模型建立中,污染混合物的毒性实验数据仍是不可或缺的。当前大气污染健康风险评价在混合物层面仍依赖于简单加和假设,集中在多环芳烃混合物相对于BaP的致癌风险方面,而忽略了污染物可能产生协同作用而导致低估风险。污染混合物尤其是机制复杂的混合物的联合毒性由于其生物体内靶点和作用方式不同而更难分析预测。如何实现完全不依赖于混合物毒性实验的联合毒性计算评价是污染物毒性预测面临的巨大的挑战和机遇。
综上所述,单一化合物的毒性预测方法已渐趋成熟,并已在环境化学品风险评估中得以广泛应用,而混合化学品产生的复杂联合毒性效应是环境毒理学研究的难点问题,其非实验评估方法亟待发展。而受制于联合毒性研究手段的限制,大气污染混合物健康风险评价仍局限于污染组分毒性简单相加的工作假设,难以回避风险低估所带来的问题。亟需发展一种接近环境真实混合体系中组分联合作用模式的快速毒性预测评价方法,从而更为客观地评价污染物的健康风险和潜在威胁。
发明内容
有鉴于此,本发明的主要目的在于提供一种大气污染二元混合物健康风险评价方法,以期至少部分地解决上述技术问题中的至少之一。
作为本发明的一个方面,提供了一种二元混合物的联合毒性预测方法,所述二元混合物包含摩尔浓度相接近的两组分,包括以下步骤:
步骤1:从公开数据库或文献中获得两组分分别暴露的转录组数据和与特定靶点或通路相对应的测试终点的体外测试数据,选取两组分中毒性最大的一组分作为参照体,另一组分作为非参照体,得到参照体和非参照体相互的基因表达相似度、以及参照体和非参照体分别对应于不同靶点或通路的测试终点的单体剂量-效应曲线方程;
步骤2:根据所述单体剂量-效应曲线方程分别获取对应靶点或通路非参照体和参照体的体外测试半数有效浓度,进而得到非参照体相对于参照体的相对效能因子;
步骤3:根据所述基因表达相似度和单体剂量-效应曲线方程分别计算参照体或非参照体在先暴露时对在后暴露的非参照体或参照体产生的基因剩余效应,通过所述参照体和非参照体的基因剩余效应、单体效应得到二元混合物的混合效应数值,进而得到二元混合物相对于参照体的校正相对效能因子;
步骤4:利用所述非参照体相对于参照体的相对效能因子和二元混合物相对于参照体的校正相对效能因子之间变化的比率表征二元混合物的联合毒性效应的强弱和作用模式。
作为本发明的另一个方面,提供了一种大气污染二元混合物健康风险评价方法,包括以下步骤:
步骤A,通过文献或实际环境监测搜集大气污染物数据,获得一段时期内大气样本里若干种污染物浓度,选取摩尔浓度相接近的两种污染物作为待评价的二元混合物;
步骤B,利用步骤1~4所述的联合毒性预测方法,将二元混合物中毒性最大的污染物作为参照体,另一污染物作为非参照体,分别获取二元混合物相对于参照体的校正相对效能因子、以及非参照体相对于参照体的相对效能因子;
步骤C,将所述二元混合物中的非参照体在实际环境中的大气浓度作为环境暴露浓度,将所述环境暴露浓度换算为体外等价浓度;
步骤D,利用所述二元混合物相对于参照体的校正相对效能因子和非参照体相对于参照体的相对效能因子,结合非参照体的体外等价浓度和参照体的毒性起效浓度,针对相应的靶点或通路,进行二元混合物健康风险评估。
基于上述技术方案可知本发明大气污染二元混合物健康风险评价方法具有如下有益效果:
(1)本发明在混合物环境健康风险评价联合毒性评估方面,充分考虑污染物基因表达谱变化所反应的复杂毒性机制,通过将其映射到特定通路或靶点的体外测试数据,获得协同、拮抗等二元混合物的联合毒性作用模式,相比现有环境健康风险评价中污染组分毒性简单加和模型获得的混合效应估算值更接近真实情况;
(2)本发明首次提出基于联合毒性预测的环境污染健康风险评估的技术方案,充分结合现有的转录组和生物测试信息,通过基因剩余效应法估算等比例二元混合物联合毒性,耦合污染物体外等价浓度模拟数据进行混合暴露的健康风险评估,应用模拟抽样获得污染物体外等价浓度中不确定因素的分布范围,最终得到科学合理的健康风险估算值;
(3)本发明所述联合毒性评价方法直接建立在污染物单体的转录组数据和体外活性测试数据上,不需要提前获得混合体系的毒性实验数据,便可预测二元混合物联合毒性并进一步其评价健康风险,节约人力和物力成本,方法快速,效率高,该方法在实际环境混合物健康风险评价、环境安全性评估和混合物毒性预测等领域具有广阔的应用前景。
附图说明
图1为本发明实施例1采用大气污染物联合毒性预测进行健康风险评价示意图;
图2为本发明实施例1中多环芳烃苯并[a]蒽(BaA)和苯并[a]芘(BaP)剂量-效应曲线示意图;
图3为本发明实施例1中BaA单体和二元混合物的健康风险分布密度和累积概率示意图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚明白,以下结合具体实施例,并参照附图,对本发明作进一步的详细说明。
根据对现有大气污染健康风险评价的计算模型的文献调研,可见所涉及的方法或技术存在缺点。基于简单加和的毒性预测模型,建立在污染物之间独立作用的假设上,忽略了污染物之间的混合效应。往往低估毒性效应。本发明的目的是要提供一种基于结合转录组数据和特定通路体外生物测试信息的混合物联合毒性预测方法,不依赖先验混合毒性测试,通过耦合大气污染物暴露数据对大气污染健康进行评价。
本发明基本原理是,对于特定靶点或通路,两种污染物同时暴露引发的混合效应可认为是其中一种单体暴露产生的基因剩余效应与另一种单体暴露产生的综合影响;因此可结合二元混合物在基因表达相似度和二元混合物单体的在体外测试数据,估算二元混合物的混合效应。二元混合物的联合毒性模式和强弱通过污染物单体的相对效能因子和二元混合物的校正相对效能因子之间变化的比率表征。最后,通过二元混合物的校正相对效能因子和大气污染暴露数据便获得大气混合污染健康风险的分布密度和累积概率。
文献和专利检索结果表明,本发明首次提出基于联合毒性预测的环境污染健康风险评估的计算解决方案。在本发明完成之前,还未发现将基于转录组数据和生物测试信息,并耦合污染物体外等价浓度模拟数据用于大气污染混合物健康风险评价的报道。
首先,本发明提供了一种二元混合物的联合毒性预测方法,该二元混合物包含摩尔浓度相接近的两组分,例如摩尔浓度比例介于(1∶1.2)~(1.2∶1)之间的两组份,其包括以下步骤:
步骤1:从公开数据库或文献中获得两组分分别暴露的转录组数据和与特定靶点或通路相对应测试终点的单体体外测试数据,根据文献报道或体外测试数据选取两组分中毒性最大的一组分作为参照体ref,另一组分作为非参照体i,其中:
利用该转录组数据得到参照体ref和非参照体i相互的基因表达相似度:
其中Gi和Gref分别为非参照体i和参照体ref的差异基因集合,ri,ref表示非参照体i和参照体ref共同差异基因的数量与参照体ref总体差异基因数量的比值,rref,i表示非参照体i和参照体ref共同差异基因的数量与非参照体i总体差异基因数量的比值,且0≤ri,ref,rref,i≤1;
利用单体体外活性测试数据,即为不同给药剂量下靶点或通路响应的比例,得到参照体ref和非参照体i分别对应于不同靶点或通路的单体剂量-效应曲线方程。
步骤2:根据所述单体剂量-效应曲线方程分别获取对应靶点或通路k非参照体i和参照体ref的半数有效浓度,计算其比值进而得到非参照体i相对于参照体ref的相对效能因子。
具体的,本步骤包括:
由单体剂量-效应曲线方程分别获取引发靶点或通路k响应50%时非参照体i和参照体ref的单体浓度,即半数有效浓度和计算参照体ref和非参照体i的半数有效浓度的比值,即为非参照体i相对于参照体ref的相对效能因子
步骤3:根据基因表达相似度和单体剂量-效应曲线方程分别计算参照体ref或非参照体i在先暴露时对在后暴露的非参照体i或参照体产生的基因剩余效应,通过所述参照体ref和非参照体i的基因剩余效应、单体效应得到二元混合物的混合效应数值,进而得到二元混合物相对于参照体ref的校正相对效能因子。
具体地,本步骤包括以下子步骤:
其中ri,ref为非参照体i相对于参照体ref的基因相似度;为两倍量参照体ref暴露引发特定靶点或通路k响应的百分比;表示参照体ref和非参照体i生物学功能相同部分的效应;表示参照体ref相对于非参照体i而言生物学功能不同部分的效应;
其中rref,i为参照体ref相对于非参照体i的基因相似度;为两倍量非参照体i暴露引发特定靶点或通路k响应的百分比;表示非参照体i和参照体ref生物学功能相同部分的效应;(1-rref,i)fi c,k表示非参照体i相对于参照体ref而言生物学功能不同部分的效应;
子步骤34:综合考虑在浓度C下,非参照体i和参照体ref依次暴露引发特定靶点或通路k的混合效应数值与参照体ref和非参照体i依次暴露引发特定靶点或通路k的混合效应数值计算得到非参照体i和参照体ref同时等量暴露引发的混合效应数值
子步骤35:根据子步骤31~34计算在0~20倍浓度范围内等间隔的多个(例如10000个)浓度C,等量的非参照体i和参照体ref同时暴露下引发靶点或通路k响应的混合效应数值获得二元混合物的剂量-效应关系;
步骤4:二元混合物的联合毒性效应的强弱和作用模式可通过如式(8)以非参照体i相对于参照体ref的相对效能因子和二元混合物相对于参照体ref的校正相对效能因子之间变化的比率表征。即在有/无参照体ref存在下,二元混合物相对等量单体暴露引发特定靶点或通路k毒性的变化率:
其中若绝对值越大,表示在参照体ref存在情况下,参照体ref和非参照体i共同暴露产生联合毒性效应相比非参照体i毒性变化越显著;若为正数,说明在参照体ref存在下非参照体i毒性增强,产生协同效应;若为负数,说明参照体ref存在下非参照体i毒性减弱,产生拮抗效应;若为0,说明参照体ref存在下非参照体i毒性无变化,该二元体系为简单加和效应。
本发明基于前述二元混合物联合毒性预测方法,还提供了一种大气污染二元混合物健康风险评价方法,包括以下步骤:
步骤A,通过文献或实际环境监测搜集大气污染物数据,获得一段时期内大气样本中若干种污染物浓度,选取一组污染物作为待评价的混合组分,进一步地,选取摩尔浓度比例介于1∶1.2~1.2∶1之间的两种污染物作为待评价的二元混合物。
步骤B,利用前述的联合毒性预测方法,将二元混合物中毒性最大的污染物作为参照体ref,另一污染物作为非参照体i,分别获取二元混合物相对于参照体ref的校正相对效能因子以及非参照体i相对于参照体ref的相对效能因子
步骤C,将二元混合物中的非参照体i在实际环境中大气浓度作为环境暴露浓度,将环境暴露浓度换算为体外等价浓度;
具体地,待评价的二元混合组分中每种非参照体i的体外等价浓度ECi为:
其中QCi为一污染物的在实际环境中的大气浓度,即环境暴露浓度;AB为中国成年人呼吸率;f为污染物接触肺泡区的比例,若QCi为污染物气相浓度,则f=1,若QCi为污染物颗粒相浓度,则f为颗粒物肺泡沉积分数;SA为肺泡区域表面积;PSA为体外测试实验中细胞培养板的单孔表面积;V为体外测试实验中细胞培养板的单孔体积。
具体地,本步骤包括以下子步骤:
子步骤D1:由参照体ref的单体剂量-效应曲线方程获取在引发靶点或通路k响应的百分比时对应的ref单体浓度即参照体ref的毒性起效浓度;可通过参照体ref单体剂量-效应曲线方程推导得到,如公式(10):
子步骤D4:对于靶点或通路k,通过对非参照体i环境暴露浓度QCi和中国成年人呼吸率AB进行蒙特卡洛抽样模拟,获得和的分布密度和累积概率,表征非参照体i和参照体ref组成的二元混合物的健康风险。RA在累积概率的第99百分位数越大说明健康风险越大。
为了能够更清楚地理解本发明的技术内容,以下通过具体实施例结合附图进一步说明本发明的技术方案。需要说明的是,下述实施例仅是用于说明本发明,而不是对本发明的限制。
实施例1
多环芳烃作为一种普遍存在于大气颗粒物中的典型有毒成分,由于其已知的致癌特性和不良健康影响,三种多环芳烃同类物已被美国有毒物质和疾病登记机构(ATSDR,2015)列为优先污染物的前十名。本实施例通过整合转录组数据和体外测试数据,以大气颗粒物中多环芳烃混合污染为例说明大气污染二元混合物健康风险评价新方法,流程如图1所示,具体包括以下步骤:
(1)大气颗粒物污染数据搜集。
从文献中获得2015年8月至2016年3月北京大气颗粒相上多环芳烃的浓度(Sci.Total Environ.,2017,593-594,390-398)。其中两种多环芳烃苯并[a]芘(BaP)和苯并[a]蒽(BaA)的污染水平在非供暖季在1∶1附近,表1列出了4组浓度比例在(1∶1.2)~(1.2∶1)范围典型污染数据。选取BaP和BaA作为待评价二元混合物,其中BaP在环境中广泛存在并较稳定,是IARC中唯一归类为人类致癌物的多环芳烃(IARC,2010),被认为是多环芳烃中毒性最大的单体。已有文献报道BaP可调节PAH其他单体在人体细胞中形成DNA加合物的效率(Toxicology,2014,321,27-39.)。因此,选择BaP作为二元混合组分的参照体,BaA作为非参照体,二者的化学结构如下式所示。
(2)二元混合物联合毒性预测。
首先,从文献中获得BaP和BaA分别经口腔灌胃给药后,28天雄性大鼠肺部组织样本的基因芯片数据(Arch.Toxicol.,2016,90,2461-2480;Toxicol.Sci.,2012,129,213-224)。原始数据在GEO数据库(https://www.ncbi.nlm.nih.gov/geo/)的编号为GSE43438和GSE51321。使用R语言的affy包读取基因表达芯片数据并生成表达矩阵,RMA函数对基因数据进行标准化处理,limma包进行差异基因提取,保留校正p值(FDR)≤0.05的基因作为显著差异表达的基因。BaP和BaA三个化合物的差异基因数量分别为1736和169个。其次,根据其明显调节的基因,从ToxCast毒性数据库(https://actor.epa.gov/dashboard/)搜索获得BaA和BaP在AhR(AEID:806)、p53(AEID:1321)和NF-κB(AEID:1346)三种通路上的高通量体外测试数据。BaP和BaA的AC50(μM)值和单体剂量-效应曲线系数如表2所示,单体剂量-效应曲线如图2所示。
采用本发明所述方法,获得混合组分的高通量体外生物测试数据和基因表达数据后,根据公式(1a)和(1b)计算二元混合物的组分BaP和BaA基因表达相似度rBaA,BaP和rBap,BaA,分别为0.005和0.050。
根据公式(3)~(6),分别计算在0~20倍浓度范围内等间隔的10000个浓度点C等量暴露下,BaP单体和BaA单体在AhR、p53和NF-κB三种通路k上响应的百分比以及二元混合物所产生混合效应具体数值即二元混合物的剂量-效应关系。并反推得到时二元混合物相对于BaP单体的结果如表3所示。
根据公式(7)~(8),计算二元混合物相对于BaP单体的校正相对效能因子以及相对效能因子和校正相对效能因子之间的相对变化率得到三种通路上和的值,分别为0.346、0.15和0.14,且都为正。说明BaP存在情况下,BaA的毒性增强。其中,p53通路上的预测结果与基于实验方法的二元混合物毒性研究结果一致(Mutat Res.2007,620,123-134;Toxicology.2011,279,36-44)。
(3)二元混合物健康风险评估。
根据公式(9)将待评价的二元混合组分中污染物的环境暴露浓度QCBaA转换为体外等价浓度ECBaA。其中QCBaA为污染物BaA在实际环境中大气颗粒物上的浓度(nmol/m3)(表1);由于污染物浓度数据来源于大气颗粒物,因此f为预测的颗粒物肺泡沉积分数(0.2350),预测值来自于肺计量模型软件MPPD v3.04(https://www.ara.com/products/multiple-path-particle-dosimetry-model-mppd-v-304);SA为肺泡区域表面积(75m2);AB为中国成年人呼吸率(14.5~18m3/day);PSA为ToxCast体外测试中使用1536孔板的单孔表面积(0.023cm2);V为ToxCast体外测试中使用1536孔板的单孔体积(5μL)。
根据公式(11),体外等价浓度ECBaA联合相应通路的和计算得到二元混合物通路激活比率和对ECBaA中不确定变量QCBaA和AB通过蒙特卡洛抽样模拟获得RA(表征健康风险)的分布密度和累积概率,其累积概率第99百分位数分别为:0.23,0.10,0.041
其中蒙特卡洛抽样模拟获得10000个BaA的环境暴露浓度QCBaA和呼吸率AB的样本,QCBaA根据BaA环境暴露浓度的均匀分布抽样获得,AB从呼吸率的对数正态分布获得,用于描述中国人群呼吸率的变化范围。
二元混合物和BaA单体的健康风险的分布密度和累积概率结果如图3所示。由图3可知BaA和BaP混合物在三种通路上的健康风险相比BaA单体均增大。由图3可知,三种通路上混合物相比BaA单体的RA分布密度靠右,说明BaP存在情况下BaA的毒性增强,混合物可能引起更大的健康风险。可见通过本发明方法可以通过污染物单体的污染暴露数据和生物测试数据的进行混合毒性预测和健康风险评估,而不依赖混合物的生物测试数据。该方法在实际环境混合暴露健康风险评价和环境安全性评估等领域具有广阔的应用前景,有望为化学品混合污染健康风险评价提供技术支撑和指导。
综上所述,本发明通过充分利用已有的体外生物测试数据和转录组信息,将污染物联合毒性作用模式预测与实际环境暴露数据相结合,以模拟环境混合物的累积健康风险分布。特别是,本发明方法没有关于化学品混合物的作用模式的预知假设,混合效应由混合物中单个组分的体外剂量-效应曲线和转录组信息估算确定。
以上所述的具体实施例,对本发明的目的、技术方案和有益效果进行了进一步详细说明,应理解的是,以上所述仅为本发明的具体实施例而已,并不用于限制本发明,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
附表:
表1大气样本中多环芳烃单体的浓度(nmol/m3)
表2三种通路上多环芳烃单体的AC50(μM)值及剂量-效应曲线系数
表3三种通路上多环芳烃二元混合物的AC50(μM)及BaP的AC10(μM)
Claims (10)
1.一种二元混合物的联合毒性预测方法,所述二元混合物包含摩尔浓度相接近的两组分,其特征在于,包括以下步骤:
步骤1:从公开数据库或文献中获得两组分分别暴露的转录组数据和与特定靶点或通路相对应的测试终点的体外测试数据,选取两组分中毒性最大的一组分作为参照体,另一组分作为非参照体,得到参照体和非参照体相互的基因表达相似度、以及参照体和非参照体分别对应于不同靶点或通路的测试终点的单体剂量-效应曲线方程;
步骤2:根据所述单体剂量-效应曲线方程分别获取对应靶点或通路非参照体和参照体的体外测试半数有效浓度,计算其比值进而得到非参照体相对于参照体的相对效能因子;
步骤3:根据所述基因表达相似度和单体剂量-效应曲线方程分别计算参照体或非参照体在先暴露时对在后暴露的非参照体或参照体产生的基因剩余效应,通过所述参照体和非参照体的基因剩余效应、单体效应得到二元混合物的混合效应数值,进而得到二元混合物相对于参照体的校正相对效能因子;
步骤4:利用所述非参照体相对于参照体的相对效能因子和二元混合物相对于参照体的校正相对效能因子之间变化的比率表征二元混合物的联合毒性效应的强弱和作用模式。
2.根据权利要求1所述的联合毒性预测方法,其特征在于,所述二元混合物中参照体ref和非参照体i的摩尔浓度比例为(1∶1.2)~(1.2∶1)。
5.根据权利要求1所述的联合毒性预测方法,其特征在于,步骤3具体包括以下子步骤:
其中,ri,ref为非参照体i相对于参照体ref的基因相似度,rref,i为参照体ref相对于非参照体i的基因相似度;为两倍量参照体ref暴露引发特定靶点或通路k响应的百分比;为两倍量非参照体i暴露引发特定靶点或通路k响应的百分比;
7.一种大气污染二元混合物健康风险评价方法,其特征在于,包括以下步骤:
步骤A,通过文献或实际环境监测搜集大气污染物数据,获得一段时期内大气样本里若干种污染物浓度,选取摩尔浓度相接近的两种污染物作为待评价的二元混合物;
步骤B,利用如权利要求1至6任意一项所述的联合毒性预测方法,将二元混合物中毒性最大的污染物作为参照体,另一污染物作为非参照体,分别获取二元混合物相对于参照体的校正相对效能因子、以及非参照体相对于参照体的相对效能因子;
步骤C,将所述二元混合物中的非参照体在实际环境中的大气浓度作为环境暴露浓度,将所述环境暴露浓度换算为体外等价浓度;
步骤D,利用所述二元混合物相对于参照体的校正相对效能因子和非参照体相对于参照体的相对效能因子,结合非参照体的体外等价浓度和参照体的毒性起效浓度,针对相应的靶点或通路,进行二元混合物健康风险评估。
8.根据权利要求7所述的大气污染二元混合物健康风险评价方法,其特征在于,步骤A中,选取摩尔浓度的比例为1∶1.2~1.2∶1之间的两种污染物作为待评价的二元混合组分。
10.根据权利要求7所述的方法,其特征在于,步骤D具体包括以下子步骤:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811500291.5A CN109583662B (zh) | 2018-12-07 | 2018-12-07 | 大气污染二元混合物健康风险评价方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811500291.5A CN109583662B (zh) | 2018-12-07 | 2018-12-07 | 大气污染二元混合物健康风险评价方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109583662A CN109583662A (zh) | 2019-04-05 |
CN109583662B true CN109583662B (zh) | 2022-09-02 |
Family
ID=65929406
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811500291.5A Active CN109583662B (zh) | 2018-12-07 | 2018-12-07 | 大气污染二元混合物健康风险评价方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109583662B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113205859B (zh) * | 2021-03-15 | 2023-02-17 | 安徽建筑大学 | 一种用于定量表征复合污染物联合毒性相互作用的面积均值法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102867116A (zh) * | 2012-08-30 | 2013-01-09 | 常州大学 | 基于健康风险的区域水环境优先控制污染物筛选方法 |
CN104751373A (zh) * | 2015-03-23 | 2015-07-01 | 国网重庆市电力公司电力科学研究院 | 计及污染气体排放风险的环境经济调度方法 |
CN106153510A (zh) * | 2015-03-25 | 2016-11-23 | 日本电气株式会社 | 大气污染物浓度的估算方法和装置 |
CN107832579A (zh) * | 2017-11-08 | 2018-03-23 | 中科宇图科技股份有限公司 | 一种参数本地化的水环境健康经济损害评价方法及*** |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20140188495A1 (en) * | 2012-12-30 | 2014-07-03 | Nanjing University | Method for assessing and managing health risks associated with heavy metal pollution |
-
2018
- 2018-12-07 CN CN201811500291.5A patent/CN109583662B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102867116A (zh) * | 2012-08-30 | 2013-01-09 | 常州大学 | 基于健康风险的区域水环境优先控制污染物筛选方法 |
CN104751373A (zh) * | 2015-03-23 | 2015-07-01 | 国网重庆市电力公司电力科学研究院 | 计及污染气体排放风险的环境经济调度方法 |
CN106153510A (zh) * | 2015-03-25 | 2016-11-23 | 日本电气株式会社 | 大气污染物浓度的估算方法和装置 |
CN107832579A (zh) * | 2017-11-08 | 2018-03-23 | 中科宇图科技股份有限公司 | 一种参数本地化的水环境健康经济损害评价方法及*** |
Non-Patent Citations (2)
Title |
---|
城市空气质量健康风险评估研究进展;张莉等;《地理科学进展》;20060610(第03期);全文 * |
雾霾重污染期间北京居民对高浓度PM_(2.5)持续暴露的健康风险及其损害价值评估;谢元博等;《环境科学》;20140115(第01期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN109583662A (zh) | 2019-04-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Graney et al. | Concentrations and solubility of metals from indoor and personal exposure PM2. 5 samples | |
Shi et al. | Combined source apportionment, using positive matrix factorization–chemical mass balance and principal component analysis/multiple linear regression–chemical mass balance models | |
Romanazzi et al. | PM10 size distribution of metals and environmental-sanitary risk analysis in the city of Torino | |
Wickramasinghe et al. | PM10-bound polycyclic aromatic hydrocarbons: Biological indicators, lung cancer risk of realistic receptors and ‘source-exposure-effect relationship’under different source scenarios | |
Zavala et al. | The Gillings Sampler–An electrostatic air sampler as an alternative method for aerosol in vitro exposure studies | |
Tian et al. | Using an improved Source Directional Apportionment method to quantify the PM2. 5 source contributions from various directions in a megacity in China | |
CN109583662B (zh) | 大气污染二元混合物健康风险评价方法 | |
Santos et al. | Biomonitoring of metal contamination in a marine prosobranch snail (Nassarius reticulatus) by imaging laser ablation inductively coupled plasma mass spectrometry (LA-ICP-MS) | |
CN112304813B (zh) | 一种大气颗粒物健康风险来向源解析方法 | |
Adgate et al. | The use of isotope ratios to apportion sources of lead in Jersey City, NJ, house dust wipe samples | |
Sysalova et al. | Mobility assessment and validation of toxic elements in tunnel dust samples—Subway and road using sequential chemical extraction and ICP-OES/GF AAS measurements | |
Yajuan et al. | Mixed uncertainty analysis of polycyclic aromatic hydrocarbon inhalation and risk assessment in ambient air of Beijing | |
Shi et al. | Source apportionment of synchronously size segregated fine and coarse particulate matter, using an improved three-way factor analysis model | |
Liu et al. | Physically constrained source apportionment (PCSA) for polycyclic aromatic hydrocarbon using the Multilinear Engine 2-species ratios (ME2-SR) method | |
Chio et al. | Health risk assessment for residents exposed to atmospheric diesel exhaust particles in southern region of Taiwan | |
Moolla et al. | Occupational health risk assessment of benzene and toluene at a landfill site in Johannesburg, South Africa | |
Ebersviller et al. | Gaseous VOCs rapidly modify particulate matter and its biological effects–Part 2: Complex urban VOCs and model PM | |
Morris et al. | No excess risk of adverse birth outcomes in populations living near special waste landfill sites in Scotland | |
Teoldi et al. | Air quality in the Olona Valley and in vitro human health effects | |
Gadkari et al. | Source apportionment of personal exposure of fine particulates among school communities in India | |
Ulbrich et al. | Interpretation of organic components from positive matrix factorization of aerosol mass spectrometric data | |
Gaskill et al. | Binary mixtures of polycyclic aromatic hydrocarbons display nonadditive mixture interactions in an in vitro liver cell model | |
Konarski et al. | Comparison of urban and rural particulate air pollution characteristics obtained by SIMS and SSMS | |
Hopke | Applying multivariate curve resolution to source apportionment of the atmospheric aerosol | |
Dubé et al. | Representativeness of laboratory sampling procedures for the analysis of trace metals in soil |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |