中国预防医学杂志万维

刘晓端曾太文李奇房珂玮

（中国地质科学院生物环境地球化学研究中心，北京100037）

郭艳梅王峙王洪权吴涛刘晓芬

（北京市房山区预防医学研究中心，北京102400）

摘要本文根据氟元素迁移演化的时空特点，从煤－水－空气－农作物－人体完整的循环体系入手，基本查明病区内含氟煤层的分布、区域地质构造和变质作用对煤层物质成分的影响以及煤中主要微量元素成分，高氟煤是研究区致病的主要物源；系统地测定了环境介质中氟化合物，估算出氟在环境中的总体暴露量和人体摄入氟量，认为不科学的燃煤方式使居室空气中氟含量严重超标是导致燃煤型氟中毒的主要因素；燃煤过程中释放出的碳氢化合物及其他的无机化合物充满室内，经过复杂的化学反应产生了氟化氢和氟氯代烷烃化合物，包括CCl3F（F11）、CClF2（F112）和C2Cl3F3（F113），构成侵害人体的有害物质，特别是燃煤中形成的氟氯代烷烃化合物，是除工业生产的氟氯代烷烃化合物外的一个重要发现，对于人类的影响不亚于工业污染。初步掌握了氟元素生物地球化学环境特点及其转移过程和致病因素，厘定了地球化学模型和简单的回归方程。

关键词燃煤型氟中毒背景成因

燃煤型氟中毒是一种典型的地球化学疾病，这种疾病的流行与病区的地质、地球化学条件、气候变化、生活方式和习惯以及社会的经济发展等综合因素密切相关。在高氟环境中生活的人大量吸收氟与氟化物，从而破坏体内钙、磷的正常代谢，影响中枢神经的正常活动，造成肝、肾损伤，智商降低，引起氟斑牙、氟骨症等，以至于失去劳动能力。这种疾病在辽宁、河北、北京、山西、河南、陕西、湖北、湖南、云南、贵州、四川、广西、江西、浙江等14个省、市、自治区有较大范围出现，影响到约5000万人口的身心健康，严重地阻滞了这些地区农村经济的发展和农民群众的脱贫致富。为阻断氟元素对机体的侵害，提高人们的健康水平，更好地促进经济的持续发展和提前实现脱贫致富计划，摸索出有效的防治对策，选择北京市房山区堂上村、莲花庵村（山区）和路村（平原）三个重点村庄进行了燃煤型氟中毒的环境地球化学、流行病学调查及防治对策研究。

按照氟元素迁移演化的时空特点，从煤－水－空气－农作物－人体完整的循环体系入手，全面准确地把握氟源、携氟介质、氟中毒疾病程度三者之间的剂量－效应关系，基本查明病区内含氟煤层的分布、区域地质构造和变质作用对煤层物质成分的影响以及煤中主要微量元素成分，高氟煤是本区致病的主要物源；系统地测定了环境介质中氟化合物，估算出氟在环境中的总体暴露量和人体摄入氟量，认为不科学的燃煤方式使居室空气中氟含量严重超标是导致燃煤型氟中毒的主要因素。燃煤过程中释放出的碳氢化合物及其他的无机化合物充满室内，经过复杂的化学反应产生了氟化氢和氟氯代烷烃化合物，包括CCl3F（F11）、CClF2（F12）和C2Cl3F3（F113），构成侵害人体的有害物质，特别是燃煤中形成的氟氯代烷烃化合物，是除工业生产的氟氯代烷烃化合物外的一个重要发现，对于人类的影响不亚于工业污染。初步掌握了氟元素生物地球化学环境特点及其转移过程和致病因素，厘定了地球化学模型和简单的回归方程。

1研究工作方法

环境地球化学样品的采集和测定

在研究工作区内分别采集土壤、水、煤、拌煤土、炉灰、灰尘、粮食、室内空气粉尘和空气样品。除空气粉尘和空气以外的其他样品，经粉碎过200目筛后，测定总F、水溶性F、SiO2、Al2O3、Na2O、K2O、CaO、MgO、Fe、S、As、Sb、Bi、Hg、Se、V、Mn、Cu、Pb、Zn、Cd、Cr、Co、Ni、W、Mo、Ag、B、Be、Li、Ba和Ce。

空气样品由中国科学院大气物理所大气化学室用GC-MS测定其中的无机和有机成分。

粉尘样品由长春科技大学测定其中的氟含量。

流行病学调查采样和检测

以现状调查、不同地区病情对比和适当兼顾回顾性定群的调查方法，采集病区氟中毒病人和正常人的静脉血、晨尿、头发、指甲等生物材料，测定F，并同时登记相应的病况及体态特征，儿童最早氟斑牙出现的年龄以及病区的内外环境特征因素等。

2结果和讨论

区域地质背景

房山区处于祁吕—贺兰山字形构造东翼反射弧的南东侧与新华夏系交接部位，除缺失震旦系、上奥陶统、志留系、泥盆系、下石炭统及上白垩统外，其他地层均有出露，其中石炭—二叠系和侏罗系为本区的主要含煤地层。煤质以无烟煤为主，主要分布在百花山向斜和北岭向斜中。百花山向斜煤产地在大安山、金鸡台、史家营、莲花庵、秋林铺、碾子沟、堂上、宝儿水、东泥洼直至东村一带。北岭煤产地北起磁家务，南到周口店，西至南窖附近，东至房山花岗岩体的范围内。

煤中含有氟和砷等有害元素。郑宝山（1988）分析了北京市房山区卫生防疫站提供的12个样品，测定的氟含量为274μg/g。王云钊等（1990）对京西地区煤层分析结果表明，石炭—二叠纪煤层含氟为，侏罗纪煤层氟含量为。我们于1996年分别从石炭—二叠纪和侏罗纪煤层底板至顶板系统采集了38个样品，以厚度加权法，求得房山地区煤层中氟砷含量（表1）。

由上可见，房山地区石炭—二叠纪和侏罗纪煤层含氟、砷都较高。

此外，煤中S偏低，而Cu、Pb、Zn、Ni、Co等金属元素含量并不低。还含有Sb、Bi、Hg、Cd等有害元素。

表1房山区F、As元素含量统计表

注：世界均值F含量为80μg·g-1,As含量为5μg·g-1。

图1北京市房山区氟异常分布图

1990年曾太文等在“北京市山区地球化学图说明”中求得房山区土壤氟平均值为μg/g，以800μg/g为异常下限圈出京山90H5-75号氟异常分布在张坊、南尚乐及长操一带，面积约120km2，峰值2315μg/g；以1mg/L为下限圈出水氟异常约4km2，分布在良乡地区（图1）。

地氟病患病情况及流行病学调查统计

房山区燃煤型氟中毒影响到全区2/3的人口，约计40万人群，遍布山区和平原（图2）。堂上村、莲花庵村地处百花山南麓，处于背风的沟谷中，年平均气温约9℃，日平均气温在－5～24℃之间，无霜期仅150天左右，植被繁茂，气候寒冷。人均年烧煤～2t，最多可达5t，地炉敞烧，煤烟及粉尘排于室内，污染室内空气和粮食及蔬菜等，加之过山风在背风坡的堆积作用和稳定山谷风所形成的逆温层使当地氟污染气体聚集，浓度增高，两村氟中毒较为严重。路村经济条件略好，烧房山产的高氟、高砷煤，用量较少，仅在冬季取暖，炉灶与居室绝然分开的占多数。全区发现氟中毒症状有差异，以山区为重。主要表现为氟斑牙、牙缺损、头晕、头痛、疲倦、失眠、记忆减退、肢体麻木、肌肉及关节疼痛、骨骼变形等症状，严重的失去劳动能力，生活不能自理。

（1）氟斑牙发病率

氟斑牙在三个村庄较为普遍，根据《中国地方性氟中毒防治试行标准（1981年）》，分为白垩、着色、缺损三型，每型又分三度，统计结果见表2。

表2氟斑牙发病率统计

图2北京市房山区燃煤型氟中毒病区分布图

氟斑牙伴有牙龈萎缩以及黑斑等。发病年龄为3～72岁，15岁以前以白垩型和着色型为主，15岁以后缺损型多见。缺损率山区高于平原，其发病部位多见于近中切齿、侧切齿，磨齿较少见。

（2）氟骨症发病率

本区氟骨症症状主要表现在关节疼痛、功能性障碍和肢体变形等。三个村比较，堂上村较重，各村发病率和发病类型见表3和表4。

表3氟骨症发病率

表4氟骨症类型统计

（3）尿氟水平

分别在冬季和夏季收集患者晨尿，以氟离子选择电极法检测氟含量，见表5。

表5尿氟含量（mg/L）

1982年北京医科大学在北京地区正常人群中调查485人，获得正常人群尿氟值为（±）mg/L。可见，在燃煤污染严重的房山区三个村庄中，尿氟分别是正常值的、、倍。在置信水平下，与北京市正常值比较具有显著性差异。

（4）血氟水平

采取静脉血，以中子活化法检测血氟含量，结果见表6。

表6血氟含量

病区血氟值较北京市正常值高出两倍以上，且在置信水平下，有显著性差异。

总之，病区地氟病患者的尿氟、血氟都较正常值高，但年龄组之间及性别无明显的差异。统计结果还表明，尿氟和血氟具有显著相关性。

环境地球化学物质中氟水平统计

1991年曾太文等人对房山区西部蒲洼、十渡、张坊、南尚乐、长沟等地土壤、玉米中的氟含量进行了研究，结果如表7。

将上述地区土壤中氟含量与玉米中氟含量作相关分析，其相关系数r=，大于，具有显著相关性。由此推论，张坊土壤为μg/g，玉米含氟量仅μg/g，而堂上村、莲花庵村、路村土壤氟含量分别是606μg/g、491μg/g和570μg/g，所引致玉米氟含量最高也不过μg/g，不可能导致氟中毒。可是三个村庄环境介质中氟水平统计结果（表8）中，玉米含氟～μg/g，显然，除土壤因素外，与玉米收获后的空气氟污染有直接关系。

表7房山区西部土壤、粮食中氟含量（μglg）

表8环境介质中氟含量水平统计（μg/g）

1982年胡永增等对堂上村的调查发现，当地居民习惯将玉米、辣椒等挂在房梁上，使用敞开式炉灶烘烤，粮食氟含量变化如表9。显然随着烘烤时间的加长，粮食中的氟含量成百倍地增长，极大值可达460～8000μg/g，即使洗涤也保留下12～99倍的氟含量。在食物的烘干过程中，煤燃烧释放出大量的氟及其化合物沉积在上面，食用后造成氟中毒。这就清楚地表明，食用污染的粮食是摄入氟的途径，而粮食中的氟则直接来源于室内的空气。

表9粮食烘烤后氟含量变化统计表（μg/g）

随着改革开放的深入，农民的生活水平有了很大的提高，虽然现在大多数农民已不在屋内用煤烘烤粮食，但几年前由于食用高氟食物而造成的体内氟含量过高的影响依然存在。

煤燃烧产物的主要物质组分

（1）粉尘

房山地区燃煤方式是地炉敞烧，属低温燃烧，时刻有碳氢化合物裂解、热合形成煤烟，其他种类繁多的金属及非金属无机化合物也随之排入空气。煤烟在空气中呈气溶胶态，表面吸附了大量的有害气体和液体，这些粉尘通过呼吸进入人体，随着深度的增加，对人体健康的危害增大。病区室内空气中不同粒度粉尘年平均浓度见表10。

表10室内空气中粉尘年平均浓度

表11室内粉尘中有害物质平均浓度统计（mg/m3）

三个村庄室内粉尘平均浓度均明显超过的国家标准（中华人民共和国工业企业设计卫生标准），分别是国标的、、倍。小于5μm的粉尘浓度占总粉尘浓度的～～10μm的粉尘浓度仅占总粉尘浓度的～。现将室内粉尘有害元素含量换算成空气中体积浓度，结果见表11。

按照表中列出的参考标准，空气粉尘中有害元素超标的有氟、砷两元素，分别超标倍和2～倍。

（2）空气

①室内空气中F、S02、H2S、As

室内有害气体组分浓度是按采样时室内现场温度、压力换算成标准状态下的浓度，其结果见表12。

表12室内气态有害物质平均浓度统计（mg/m3）

四种有害气体都严重超标，构成以氟为主的多种有害物质并存的生活环境。据测定，含氟化合物以氟化氢为主。

②有机氟化合物

空气抽样检测结果发现，煤的低温燃烧在室内所形成的有机化合物达数十种，这些有机化合物排入空气后，在水、金属元素和温度的影响下，经过一系列的加成聚合反应、氧化反应、同分异构、水解、降解作用逐步形成氟代烷烃类的含氟有机化合物，其中有二氟二氯甲烷（F12）、一氟三氯甲烷（F11）、三氟乙烷（F113）等氟代烷烃类化合物；还有氟代烯烃类有机化合物，如二氟二氯乙烯等。现将几个主要问题讨论如下。

煤的低温燃烧煤燃烧过程是一种类似干馏和分解的过程，形成固态的焦炭，液态的焦油、氨水，气态的氢、甲烷、－氧化碳、二氧化碳、氮、烯烃等。温度在低于160℃时馏出苯、甲苯、二甲苯等单环的有机化合物：当温度达到160～230℃时馏出的主要是酚类和少量的萘，即多环芳烃等；230～270℃时则馏出萘，即稠环芳烃；在温度升至270～360℃时所馏出的则是杂环芳香族化合物。本次所测定的化合物除以上述及的外，尚有五节杂环类的二氢呋喃、吡咯及六节单环的吡咯。因此，煤在燃烧过程中，不但释放大量的碳氢化合物、硅氧化合物，而且有金属元素、有机化合物、卤族化合物。这些化合物在温度、压力和各种组分浓度变化的条件下，化学性质在逐渐改变，既有氟氯代烷烃类的化合物，也有卤化氢等化合物出现。

表13主要含氟有机化合物年平均含量统计表（mg/m3）

经测定，无论是堂上村、莲花庵村还是平原区的路村，室内空气中的有机化合物种类均繁多，现将主要的几种含氟有机化合物列于表13中。

除表中所列10种化合物以外，尚有二甲胺（C2H7N）、甲氧基丙胺（C4H11NO）、环丙基甲醇（C4H8O）、丙烷（C3H8）、溴化氢（HBr）、氯丙烯（C3H5CI）、氯甲基丙烷（C4H9Cl）、异氰基丁烷（C5H9N）、戊二猜（C5H6N2）、二甲基吡喃酮（C7H12O2）、羟基乙基肼（C2H8N2O）、二氢呋喃（C4H60）、甲基戊炔（C6H10）、二甲基环丙烷（C5H10）、氯丙烯（C3H5Cl）、氯代丁烯醇（C4H5Cl0）、吡咯（C4H5N）、丁烯醛（C4H60）、溴甲醚（C2H5BrO）、二氟二氯乙烯（C2Cl2F2）、三氯乙烯（C2HCl3）、甲基咪唑（C4H6N2）、甲苯（C7H8）、氯甲基戊烷（C>6H13Cl）、二氯三氟乙烷（C2HCl2F3）、四氯乙烯（C2Cl4）、乙苯（C8H10）、十一烷（C11H24）、十二烷（C12H26）、萘（C10H8）等多种有机化合物。另外，煤中的K、Na、Ca、Mg、Si和重金属元素等在燃烧中释放，参与各种化学反应，形成了一个充满化学变化危机的环境。

氟氯代烷烃的形成及毒性上述几种氟氯代烷烃是煤在燃烧过程中所释放出的烷烃和氟化氢气体相互作用而形成的。

煤燃烧释放出甲烷，在其他卤素的作用下，烷烃的氢原子为卤素所取代而形成四氯化碳。

CH4+Cl2→CH3Cl+HCl

CH3Cl+Cl2→CH2Cl2+HCl（二氯甲烷）

CH2Cl+Cl2→CHCl3+HCl（三氯甲烷）

CHCl3+Cl2→CCl4+HCl（四氯化碳）

四氯化碳在金属元素催化作用下，与无水氟化氢反应得到氟里昂：

地球化学环境：农业·健康

或：

地球化学环境：农业·健康

地球化学环境：农业·健康

由于空气中存在大量的HF和CCl4，使得上述反应反复进行，室内F12的浓度不断增加。同时还产生CCl3F、CClF2、CClF3等氟氯代烷烃化合物。

氟氯代烷烃类化合物具有沸点低、常温下易挥发等特点，易从呼吸道进入机体。经实验14C—F11或14C—F12被动物吸入后，一分钟即可在血液中检出。肺泡对氟氯代烷烃的吸收率取决于他们的化学活动性质，并与他们在肺泡的分压和脂溶性呈正比。试验还表明，在给家兔吸入384mg/m3的F22一小时后，约有80%被吸收，各肺泡内血液中F22分压迅速达到平衡状态。在堂上村、莲花庵村、路村三个村检出～，年平均;～，年平均。而前苏联1974年颁布的标准，F11、F12均为小于10mg/m3。显然本区F11、F12是较高的。氟氯代烷烃能随血液分布至全身组织，并透过细胞膜进入细胞。进入机体的氟氯烷烃主要分布在肺组织，并能通过脑脊液到达脑实质。根据研究获得的成果，氟氯代烷烃类随着氯原子数量增加而毒性增大。另外，研究中还发现溴化氢存在，其浓度达到70mg/m3。由于氟溴代烷稳定性较差，因此不能排除氟溴代烷烃的存在。高浓度的氟溴烷烃可以引起呼吸道刺激、震颤、痉挛，在致死浓度动物产生急性肺水肿而死亡。

至此，我们清楚地看到，高氟煤燃烧过程中释放出的氟及其化合物不仅污染了室内空气，也使粮食中的氟含量增加，大量的氟由呼吸道和消化道进入人体。

总摄氟量估算及危害程度的评价

不论是从什么渠道进入人体的过量氟，导致体内各种生化指标的改变和生理机能的变异程度，与环境要素氟水平是有密切关系的。这里所指环境要素系与人体有关的空气、粉尘、水、粮食、蔬菜、水果等。

（1）摄氟量计算

按照ICRP肺模型，空气进入肺泡沉积份额按计算；室内粉尘中包括100～＜5μm的总体沉积份额按计算。

粮食摄入量计算为原粮×出粉率×元素保存率×元素浓度值。

食品结构模型按1982年膳食模型计算，三个村庄的特点略有改变（表14）。

表14膳食结构模型（kg/d）

计算中的主要参数见表8。

根据表8和表14参数分别计算从呼吸道和消化道进入人体的氟量，见表15。

堂上村摄氟总量最高达到，路村最少也有，分别是标准摄入量的～倍。经消化道和呼吸道进入人体氟量占总氟量的比例，见表16。

表15摄入氟量统计（mg/d）

表16总摄氟量的构成比（%）

各村构成比例虽然不同，但经两种途径的摄氟量都超过人体所需氟水平。前已述及，煤烟污染的室内空气含氟量超过国家标准的几十至百倍，而室内的粮食及蔬菜被污染现象也相当严重，使玉米、辣椒等含氟量升至几十毫克至几百毫克每千克。通过室内空气和食物进入人体的氟量严重超标，因此，房山地区的燃煤型氟中毒的病因是食用被氟污染的粮食和由呼吸道摄人氟量过多共同作用的结果，都是由燃煤污染造成的。由此可见，改变燃煤方式或改用低氟煤是控制氟中毒疾病流行的主要方法。

（2）摄氟量与疾病相关性

为探求燃煤型氟中毒与各项生化指标及总摄氟量之间的关系。将有关指标与总摄氟量作相关分析，结果见表17。

表17总摄氟量与生化指标相关矩阵（γ＜）

①显著性相关。

氟斑牙、牙缺损、氟骨症与总摄氟量呈明显正相关，氟骨症与牙缺损程度呈直线相关。血氟与尿氟密切相关，但与疾病发生率相关不明显，提示血氟、尿氟两者仅能作为判断疾病发生的指标，难以准确量化疾病的程度。已有研究成果表明，人体内血氟与骨氟处于动态平衡中，尿氟与血氟密切相关则有可能通过增加尿的排氟量而降低血氟，逐步地将骨氟转移到血液中经尿排出体外使之降到正常水平，从而解除氟中毒对人体的危害。

以总摄氟量为自变量，牙缺损、氟骨症发生率可以用下列回归方程表示：

牙缺损率（%）=（r=）

氟骨症患病率（%）=（r=）

由以上回归方程可以看出，总摄氟量如降到6mg/d以下，在置信水平下，牙缺损和氟骨症患病率可降至零，这与已有研究成果人体每日摄入4～6mg/d的结论是一致的。

病区划分及危害程度评估

由于生活燃煤引起人体摄氟量过度而导致中毒，从发现本病到现在已20余年，但病区的划分及危害程度的评估，目前仍无正式标准。吴德良等（1990）在研究燃煤型氟病区总摄氟量与氟斑牙的相关性时提出“5～15岁儿童缺损氟斑牙与总摄氟量、尿氟呈直线相关，较氟斑牙更为可靠，更能客观反映氟危害的严重程度。”孙玉富等（1993）认为“燃煤型氟中毒人均日摄氟量超过时氟斑牙总病率即达90%以上，因此，不能作为不同病区程度的划分指标，而仅能作为病区与非病区的划分标准。”江远福（1994年）提出燃煤型氟中毒病区轻、中重、特重病区的划分标准是以8～15岁缺损型氟斑牙患病率依次是5～、20～、40～、60～。还有其他的研究者也做了不同的划分方案。

本区调查结果表明，尽管8～15岁儿童氟斑牙患病率高达98%，而牙缺损率为零，因此，以这个年龄段作为划分病区标准是不现实的。但当不同村庄总摄氟量有差异时，总体人群发病率也明显不同，特别是牙缺损和氟骨症发病率较明显。

至于氟骨症患病率，白学信（1990）认为：“在轻、中病区大量调查氟骨症是不合算的，实际调查时拍片也有困难”。流行病学调查中常使用尿氟含量，但尿氟含量与氟斑牙无相关关系，且受气候、饮水等因素影响。室内空气氟含量只能说明室内污染程度，但不是决定氟中毒程度的指标。决定氟中毒程度的是总摄氟量，总摄氟量的多少则明显反映在发病程度上，即：总摄氟量与其生物效应具有高度的一致性，可以作为病区相互比较的基础。

根据以上作者的意见，结合本次调查结果，可以得出以下几点简单的结论：

（1）综合的环境因素包括地质、地球化学、气象条件、生活方式、氟源分布等构成燃煤型氟中毒症流行的前提条件；

（2）血氟、尿氟及氟斑牙患病率只能作为划分病区与非病区的客观指标；

（3）缺损型氟斑牙患病率则是病区发病程度的指标；

（4）氟骨症是使用困难较多但较有决定意义的划分发病程度的客观指标；

（5）人均每日总的摄氟量则是反映燃煤型氟中毒的客观指标，能全面反映出人群总体摄氟水平，也能作为不同病区之间对比的客观基础。

根据以上几点，我们拟将区域地球化学条件、氟斑牙发病率、氟斑牙缺损率、氟骨症发病率、总摄氟量等为指标，将燃煤型氟中毒流行范围划分为疑似病区、轻病区、中病区、重病区、特重病区等五种不同程度的病区，具体划分标准如下。

疑似病区：具有明显的高氟和氟积累的地球化学背景，有利于氟污染的气象条件，落后的燃煤方式等。

轻病区：每年四个月以下燃烧高氟煤，氟斑牙患病率在10%以上，成年人氟斑牙患病率在50%以上，缺损率在10%以下，氟骨症患病率在5%以下，且多为I，每日人均总摄氟量在10mg/以下。

中病区：每年6个月燃烧高氟煤，通风条件差，氟斑牙患病率为50%～70%，氟斑牙缺损率为10%～25%，氟骨症患病率为5%～12%，多为Ⅱ，人均每日摄氟量为10～15mg。

重病区：每年有6～8个月燃烧高氟煤，通风条件差，粮食蔬菜在室内储藏，氟斑牙患病率在70%以上，氟斑牙缺损率为25%～60%，氟骨症患病率为13%～35%，多为Ⅲ，人均日摄氟量为15～30mg。

特重病区：长年燃烧高氟煤，通风条件极差，粮食蔬菜室内存放，氟斑牙患病率在90%以上，成年人氟斑牙缺损率在60%以上，氟骨症患病率在35%以上，为Ⅲ。

以上划分是初步设想，尚需要更多的资料加以验证和补充。

3主要结论

（1）由于区域地质地球化学特点以及区域性变质作用，使煤层中的有机、无机化合物聚集、分散、相互反应，形成一个十分庞杂的资源库，这是地氟病发生的最重要的物源。地球化学环境是发生地氟病的根本原因，而导致地方性氟中毒的决定性因素是人体摄入氟的途径和氟存在的形态及溶解度。

（2）煤中的化学物质经低温燃烧后，释放出大量的有机和无机含氟化合物及其他有毒物质，使室内空间形成一个复杂的化学环境，这些物质经呼吸道进入人体，同时也污染了粮食和蔬菜，并被人体所吸收，造成氟中毒。这是当前房山区地氟病发病的首要病因。

（3）不科学的燃煤方式和落后的生活方式是造成室内空气污染的直接原因，也是增加人体摄氟量的重要因素。因此，改变燃煤方式，降低室内空气污染程度是防治燃煤型氟中毒的关键措施。

（4）煤在燃烧过程中经复杂的化学反应，形成大量的氟氯代烷烃化合物，其中包括CFCl3（F11）、CF2Cl（F12）和C2Cl3F3（F113），这是这次调查的重大发现。由于F11、F12较为稳定，极可能对大气臭氧层造成破坏，使得臭氧耗竭、人群皮肤癌患者急剧增加，这个影响人类生存环境的问题应引起充分的重视，而这些氟氯代烷烃化合物的形成机理有待进一步深入研究。

（5）研究了氟斑牙发病率、氟斑牙缺损率、氟骨症发病率与人均每日总摄氟量之间的相互关系，建立了用总摄氟量预报氟中毒发病率的回归方程。以地球化学条件、地氟病发病率和总摄氟量为综合指标，提出了将病区划分为疑似病区和轻、中、重、特重病区。对不同程度的病区，应当采取有针对性的防治燃煤型氟中毒的措施。

参考文献

耿精忠等.环境与建康.北京：华夏出版社，1993

吴鹏鸣等.环境空气监测质量保证手册.北京：中国环境科学出版社，1989

曹守仁.煤烟污染与健康.北京：中国环境科学出版社，1992

北京市地矿局.北京市区域地质志.北京：地质出版社，1991

北京市地矿局.北京市区域矿产总结，北京：地质出版社，1992

郑宝山.中国煤炭含氟量的研究.中国地方病防治杂志1988，第3卷，第2期

郭英廷等.贵州西部晚二叠世煤层中的有害微量元素及其分布.贵州工学院学报，1993

孙景信.煤中微量元素及其在燃烧过程中的分布特征.中国科学（A）,1998,12

国家环保局.化学品毒性法规环境数据手册.北京：中国环境科学出版社，1992

苏企洵.有机化学.北京：高等教育出版社，1988

霍亚贞.北京自然地理.北京：北京师范学院出版社，1989

戴国均.地方性氟中毒.贵阳：贵州人民出版社，1965

王云钊.氟骨症X射线诊断学图析.北京：中国环境科学出版社，1990

戴瑞庭.北京门头沟燃煤污染型氟关节病.中国地方病防治杂志1988，第3卷，第2期

吴德良.燃煤污染病区总摄氟量与氟斑牙患病的相关研究.中国预防医学杂志1990，第24卷，第1期

孙玉富.摄氟量与燃煤污染型氟中毒关系的流行病学研究.中国地方病学杂志1993，第12卷，第6期

江远福.生活燃煤污染型地方性氟中毒研究进展.中国地方病防治杂志1994，第9卷，第2期

白学信.燃煤型地氟病病区划分标准问题.中国地方病防治杂志1990.第5卷.第4期

谢妮.砷对氟中毒影响的实验研究.中国地方病学杂志1990，第13卷，第1期

周振荣.中医药对氟骨症的治疗研究.见：中国地方病学杂志1996，第15卷.第五届全国地氟病地砷病学术会议论文集

胡乔木等.中国大百科全书（化学卷I）.北京：中国大百科全书出版社

中国预防医学杂志，其实都还行，望采纳

呵呵，中华预防医学杂志是北大核心，中国预防是啥？大家都是核心而已

多吃些瘦肉、鱼，蛋、牛奶、豆及豆制品。

【综述】几种用于发现未知病毒核酸序列的技术及其应用翁康生 病毒是引发人类传染性疾病的主要病原体之一, 它们极大地威胁着人类健康。目前还存在人类尚未认知或新出现的病毒, 随时可能严重危害人类健康安全〔1 - 3〕。及早地发现,鉴别未知的或新出现的病毒, 是有效的预防和控制的先决条件之一。因此, 建立、储备、改良、发展、乃至创新应用于发现、鉴别未知或新出现病毒的技术方法是十分必要的。近20 多年来, 常采用传统的微生物学技术方法和现代分子生物学技术方法相结合的途径, 发现和鉴别未知病毒。通过细胞培养的方法分离病毒、电镜观察、用已知病毒的抗血清建立的免疫学方法作排他性检测、用已知病毒核酸序列建立的PCR、杂交等方法, 作特异核酸序列的检测、用分子生物学技术获得未知病毒核酸序列, 查询基因数据库, 检出并确定未知病毒基因组序列, 最终发现鉴别出未知病毒。对于无法用细胞分离培养的未知病毒, 有的采用免疫学与分子生物学技术相结合, 筛选获取病毒特异抗原编码基因的克隆, 进而发现鉴别出该病毒。更多的则是采用相应的分子生物学技术, 从被检样品中发现获取未知病毒的核酸序列,进而发现鉴别未知病毒。无论未知病毒是否可以用细胞培养分离, 最终对其基因组序列的测定分析, 是鉴别和判断的决定性依据之一, 而获取未知病毒的核酸序列是前提条件。从少量样品中, 从高度复杂的宿主细胞核酸物质中, 分离、扩增、获取足够量的无基因序列资料的未知病毒的核酸片段, 供进一步克隆、测序、生物信息学分析, 是用分子生物学技术发现、鉴别未知和新出现病毒的关键之一〔4〕, 也是最终测定分析, 拼接出未知和新出现病毒基因组序列的瓶颈步骤。病毒所携核酸物质有DNA 和RNA 之分, 可采用的技术方法也有所不同, 现将有关技术与其应用作一简介, 以供参考。1 代表性差异分析法代表性差异分析法是为寻找分析两个生物样品复杂的基因组间有何差异而发展建立起来的分子生物学技术方法, 并不断得到演化, 发展和应用。病毒感染宿主细胞后, 与未感染的同类细胞相比, 二者核酸物质间的差异主要在于是否存在病毒核酸。消减去二者核酸间相同序列的背景部分, 扩增、比较、选取余下可能存在差异的部分, 进一步分析以发现未知病毒的核酸序列。病毒的核酸结构各有不同, 可选用相应的代表性差异分析法, 见表1 。111 DNA 代表性差异分析法(DNA Representation differenceanalysis , DNA RDA)此方法是Lisitsyn 等〔5〕利用核酸消减杂交技术〔6〕、PCR 方法和双链DNA 热变性后互补链退火复性的二级动力学原理〔7 - 8〕作者单位:上海市疾病预防控制中心 200336表1 病毒核酸类型与各代表性差异分析法的选用病毒核酸类型DNA RDA c DNA RDA非rRNA 序列6 核苷酸引导c DNA RDAds DNA 线状√ds DNA 环状√ss RNA polyA( + ) √ss RNA polyA( - ) 3 √ds RNA polyA( - ) √ 3 负链ss RNA polyA( - ) 视病毒在宿主细胞的转录机制而定。而建立的。方法中将需分析的样品DNA(Test DNA ,T- DNA) 和对照DNA(Driver DNA ,D - DNA) 设为二组,分别用同一种限制性内切酶酶切处理,并接上5′端去磷酸化的人工接头,补齐接头后,加入与接头序列互补的引物作PCR 扩增。切除扩增产物上的人工接头后,切出的T - DNA 连上第二种人工接头,变性后与过量的变性D - DNA 杂交。通过杂交,消减去T- DNA 中与D - DNA 中同源的核酸序列,而只存在于T - DNA 中的靶序列DNA(Target DNA) 则自我退火复性,其两端连有第二种人工接头。加入与第二种接头互补的引物作PCR ,只有靶序列DNA呈指数扩增,因而得到进一步富集。进过如此重复的几个轮回后,以电泳检测比较T- DNA 和D - DNA ,将T- DNA 中呈现的差异部分作分离,克隆,序列分析。Lisitsyn 等以10μg 人淋巴细胞基因组DNA 作为D - DNA ,在相同的人DNA 中加入相当于单拷贝量的120 pg 腺病毒DNA作T- DNA。以此作为实验模型,用DNA 代表性差异分析法成功的寻找、鉴定出外加入的腺病毒DNA 序列。应用此技术,Chang 等〔9〕在艾滋病相关的卡波西肉瘤(Kaposis Sarcoma) 中发现一段类似人类疱疹病毒的基因, 并由此发现一种新的病毒HPV8。以后人们又以此技术发现鉴定了HPV6、TTV 病毒、黄热病毒样基因组、MDV 等〔10 - 13〕DNA 病毒。112 cDNA 代表性差异分析法(cDNA Representation differenceanalysis , cDNA RDA)Hubank 等〔14〕针对mRNA 所含序列相对简单的特点,提出了cDNA 代表性差异分析法。它的基本原理与DNA RDA 相同,主要不同在于,采用识别4 核苷酸序列的限制性内切酶,它的识别位点在mRNA 反转录成的cDNA 中出现的频率更高,平均酶切片段长度约256 bp ,保证了cDNA 序列群中绝大多数序列,至少被切出一个片段可扩增,供差异分析,分离鉴定。cDNA RDA 技术相对经济,可高效灵敏地用于非常少的起始材料而获得结果〔15〕。具有polyA( + ) - RNA 病毒,其核酸可类似于mRNA 分离纯化,因此可应用此技术。利用cDNA RDA技术,发现鉴定了TiV、MenV ,等〔16 ,17〕RNA 病毒。113 非rRNA 序列6 聚核苷酸引导反转录的cDNA RDA中国预防医学杂志2007 年6 月第8 卷第3 期 Chin Prev Med , June 2007 , Vol18 No13 ·317 ·© 1994-2008 China Academic Journal Electronic Publishing House. All rights reserved. ( - ) - RNA 病毒,其核酸物质不似于mRNA ,需和宿主细胞总RNA 同时分离,并用随机引物作反转录。因为宿主细胞总RNA 中rRNA 约占80 % ,由于竞争反应、靶序列信号被湮灭等原因,从这样的总RNA 抽提物中,用随机6 聚核苷酸引物引导反转录的cDNA RDA 技术,发现鉴别polyA( - ) - RNA 病毒核酸序列是困难的。Endoh 等〔18〕罗列了6 聚核苷酸所有可能的排列组合,共计4 096 个序列模式,以大鼠18S、518S、28S 等rRNA、微卫星重复序列、SARS - CoV、BI - 3 病毒等的序列数据为模型,筛选出在rRNA 序列中出现频率极低或不出现的6 聚核苷酸序列模式共96 种。将这些序列分别合成并混合后,称之为非rRNA 序列6 聚核苷酸引物。生物信息学分析96 种序列模式在哺乳动物病毒科具代表性的1 791 个病毒基因组序列中出现的频率,数据表明,非rRNA 序列6 核苷酸引物可引导绝大多数病毒的cDNA 合成。分别用非rRNA 序列6 聚核苷酸引物和随机6 核苷酸引物作cDNA 反转录效率、cDNA RDA 试验,结果表明,二类引物对人工合成的RNA (二类引物在其序列中出现的频率相似) 反转录效率几乎相等,而前者对细胞总RNA反转录效率远低于随机引物。用二类引物作cDNA RDA ,检测人工合成的RNA ,前者灵敏度是用随机引物的30 倍。在模拟实验中用非rRNA 序列6 聚核苷酸引物引导反转录,串联cDNARDA 技术,检测鉴别出感染细胞的BI - 3 和SRAS - Cov 核酸序列片段。此方法能从1μg 总RNA 中检测出3 ng 的外来RNA ,其检测灵敏度不及普通的PCR 检测方法,但对于检测鉴别在宿主细胞中复制,但不知其基因序列的poly A( - ) - RNA 病毒而言,也是一个可选择的方法。114 抑制消减杂交cDNA RDA 技术结合消减杂交和PCR 抑制作用〔19〕的技术原理,Diatchenks 等〔20〕等发展出了抑制消减杂交技术( suppressionsubtractive hybridization ,SSH) 。与前两种RDA 技术不同点在于,SSH 技术将内切酶酶切处理的T - cDNA 分为两份,分别接上不同序列的去磷酸化的接头1 和2 ,分别于过量的D - cDNA作第一轮杂交。杂交过程中两组中的单链T - cDNA 浓度趋同,T- cDNA 中的非靶序列单链cDNA 与D - cDNA 中相应序列形成杂交双链而被消减, T - cDNA 中差异表达的单链cDNA被显著富集。合并一轮杂交物,加入过量变性D - cDNA ,作第二轮杂交。合并的二组份一轮杂交物中剩下的趋同化、经消减杂交后的单链T - cDNA 能互补杂交, 可以形成: 原组内T- cDNA 单链间的杂交、T - cDNA 与D - cDNA 单链间的杂交、二组间T- cDNA 单链间的杂交。补齐杂交反应后双链cD2NA 末端,用分别与接头1 和2 的外侧部分序列互补的寡核苷酸为引物,作PCR 扩增。二组份间T- cDNA 互补单链杂交物,因两端分别具有接头1 和2 ,可被指数扩增;T - cDNA 与D -cDNA 杂交物和剩余单链T- cDNA ,因一端具接头序列,被线性扩增;而同组间T- cDNA 杂交物两端具反转重复长序列,因抑制性PCR 效应,在PCR 反应循环中分子内退火形成稳定的“锅柄结构”〔19〕而不被扩增。因此,SSH 技术通过二轮消减杂交和抑制性PCR 特异扩增,使假阳性大大降低,提高了检出低丰度靶mRNA 的灵敏度。Hu 等〔21〕应用SSH 技术,结合反转录酶的模板切换(tem2plate - switching) 功能, 以HCV RNA 阳性血清体外感染的人MOLT- 4 急性淋巴母细胞白血病T 细胞系为模型,通过反转录合成全长cDNA、抑制性消减杂交、消减的cDNA 文库构建、反相斑点杂交筛选,在被筛的96 个克隆里,T- cDNA 探针杂交呈特异阳性的16 克隆中,序列分析后得到4 个插入HCV 序列的克隆。2 非特异多重引导滚环式扩增法乳头瘤病毒、痘病毒等,其基因物质为环状DNA 分子。在事前未知基因序列的情况下,发现和鉴别这类病毒核酸序列还可选择非特异多重引导滚环式扩增法(multiply primed rolling -circle amplification ,RCA) ,扩增、分离、获取其基因片段供进一步分析。自然状况下,环状DNA 经常以滚环方式进行复制。Dean等〔22〕应用随机6 聚核苷酸作引物,加入φ29 DNA 聚合酶,以质粒DNA 和噬菌体DNA 为模型,建立了多重引物引导的滚环式扩增法。φ29 DNA 聚合酶可长距离( > 70 000 nt) 地结合于DNA模板,进行链置换DNA 合成。而随机6 聚核苷酸引物可多位点的与单链环状DNA 互补复性。在φ29 DNA 聚合酶作用下,以随机引物引导,合成与模板互补的DNA 链。当合成链延伸到与模板结合的随机引物5′端时,在φ29 DNA 聚合酶的链置换活性作用下,下游被延伸的随机引物链被“甩”出模板。而上游的延伸链继续在环状模板上复制合成。同时,被从单链环状模板上“甩”出的互补链,又成为新的模板,随机引物与之结合,在φ29 DNA 聚合酶作用下,继续以枝杈的形式进行链延伸和链置换,最后以双链DNA 串联体形式释放。用此法可使1 ng 纯pCU18 环状DNA 模板延展式地扩增至107倍。Rector 等〔23〕以此原理建立了不依赖已知的特定基因序列(非序列依赖性) 的多重引导滚环式扩增环状DNA 病毒基因组方法,并应用其扩增获取了HPV 16 的基因组DNA。在接近实样的试验样品中,由于稀释倍数和环状DNA 分子较大等原因,将HPV 16 基因组DNA 扩增了214 ×104 倍。3 病毒颗粒相关核酸的非序列依赖性PCR 扩增病毒核酸可包裹于病毒外壳内,病毒的蛋白外壳或脂膜对病毒核酸具有保护作用。而病毒颗粒具有不同于细菌或其他真核细胞的理化特性。利用这样的特点Allender 等〔24〕和Stang等〔25〕各自建立了病毒颗粒相关核酸的非序列依赖性PCR 扩增方法(sequence - independent amplification) 。两种方法的共同点在于,依据病毒颗粒小、具一定密度,用0122μm 滤器过滤、或再串上超速密度梯度离心,从样品中分离出病毒颗粒,DNA 酶酶解游离的DNA ,裂解病毒颗粒,抽提获取较纯的病毒颗粒相关核酸。Allender 等〔24〕借鉴RDA 原理,对病毒颗粒相关核酸用限制性内切酶酶切后,作非序列依赖性单引物PCR 扩增( sequence -independent single primer amplification ,SISPA) :将抽提获取的DNA或RNA 分别补齐,合成第二链DNA ,或反转录,合成双链cDNA。限制性内切酶酶切后,酶切片段两端连接一种接头,并以与接头同序列的单一寡核苷酸为引物,作PCR 扩增。扩增产物进一步克隆与序列分析。用此法检验HBV 阳性血清和GBV - B 阳性血清样品,结果在相当于106/ ml 个基因组拷贝浓度的50μl样品中,可重复试验检出相应的病毒基因片段。Stang 等〔25〕则在得到双链DNA 或双链cDNA 后加入k - 随机引物,此种引物5′端含有20 个固定序列的核苷酸,3′端则有·318 · 中国预防医学杂志2007 年6 月第8 卷第3 期 Chin Prev Med , June 2007 , Vol18 No13© 1994-2008 China Academic Journal Electronic Publishing House. All rights reserved. 核苷酸随机简并序列。与变性模板退火时,6 核苷酸随机简并序列随机地与模板相应序列互补退火,在T4 DNA聚合酶作用下作链延伸。然后在延伸产物中加入k - 随机引物中固定序列部分的20 寡核苷酸作引物,进行PCR 扩增。扩增产物电泳分析、克隆、测序。用此方法检验由实时- PCR 定量,包括病毒颗粒相关核酸及游离核酸在内的病毒基因组拷贝数为109/ ml 的Cox - 3 和MAV - 1 培养物。前者的12 克隆中,9 个克隆插入有肠道病毒的四个不同区域的同源基因片段;而6 个MAV - 1 中分离的克隆内,5 个含有99 % 同源MAV - 1 基因片段。基于病毒颗粒分离纯化、DNase 处理、病毒颗粒相关核酸的非序列依赖性PCR 扩增,获取、鉴定未知病毒的基因片段,尽管灵敏度不够高,但其实验时间较短,步骤相对简单,对于病毒拷贝数高,时间紧急的样品鉴别,是较适宜的一套方法。病毒的种类、结构、特性多种多样,感染病毒后需要检验的样品又各不相同,因此用于发现鉴别未知病毒的核酸序列的技术,也不是固定不变和完全通用的。以上的技术方法各有优缺点和适用范围。而针对扩增获取未知病毒基因组序列片断,这一发现鉴定未知病毒的分子生物学技术的要点或瓶颈,必然还会有新的改进、创新技术出现,将会更快、更灵敏、更简便、更准确的发现鉴别未知病毒。参 考 文 献〔1〕 Drosten C , Gunther S , Preiser W, et al1 Identification of novel corona2virus in patients with severe acute respiratory syndrome1 N Engl J Med ,2003 , 348 : 1967 - 19761〔2〕 ven den Hoogen BG, de Jong JC , Groen J , et al , A newly discoveredhuman pneumovirus isolated from young children with respiratory tractdisease1 Nat Med , 2001 , 7 : 719 - 7241〔3〕 Fouchier RA , Hartwig NG, Bestebroer TM, et al1 A previously unde2scribed coronavirus associated with respiratory disease in human1 ProcNatl Acad Sci U S A , 2004 , 101 : 6212 - 62161〔4〕 Muerhoff AS , Leary TP , Desai SM, et al1 Amplification and subtractionmethods and their application to the discovery at novel human viruses1 JMed Virol , 1997 , 53 : 96 - 1031〔5〕 Lisitsyn N , Lisitsyn N , Wigler M1 Cloning the differences between twocomplex genomes1 Science , 1993 , 259 : 946 - 9511〔6〕 Lamar EE , Palmer E1 Y- encoded1 Species - specific DNA in mice :evidence that the Y chromosome exists in two polymorphic forms in in2bred strains1 Cell , 1984 , 37 : 171 - 1771〔7〕 Wieland I , Bolger G, Asouline G, et al1 A method for differencecloning : geng amplification following subtractive hybridization1 Proc NatlAcad Sci USA , 1990 , 87 : 2720 - 27241〔8〕 Milner JJ , Cecchini E , Doming PD1 A kinetic model for subtractive hy2bridizationg1 Nucleic Acids Res , 1995 , 23 : 176 - 1871〔9〕 Chang Y, Cesarman E , Pessin MS , et al1 Identification of herpesvirus- like DNA sequence in AIDS - Associated Kaposi’s Sarcoma1 Scie2nce , 1994 , 266 : 1865 - 18691〔10〕 Challoner PB , Smith KT , Parker JD , et al1 Plaque - associated expres2sion of human herpesvirus 6 in multiple selerosis1 Proc Natl Acad SciUSA , 1995 , 92 : 7440 - 74441〔11〕 Nishizawa T , Okamoto H , Konishi K, et al1 A novel DNA virus (TTV)associated with elevated transaminase levels in pasttransfusion hepatitisof unknown etiology1 Biochem Biophy Res Commun , 1997 , 24 : 92 -971〔12〕 Simons JN , Pilot - Matios TJ , Leary TP , et al1 Identification of two fla2vivirus - like genomes in the GB hepatitis agent1 Proc Natl Acad SciUSA , 1995 , 92 : 3401 - 34051〔13〕 Endoh D , Cho KO , Tsukamoto K, et al1 Application of representationaldifference analysis to genomic fragments of Mark’s disease virus1 J ClinMicrobiol , 2000 , 38 : 4310 - 43141〔14〕 Hubank M, Schatz DG1 Identifying differences in mRNA - expression byrepresentational difference analysis of cDNA1 Nucleic Acids Res ,1994 , 22 : 5640 - 56481〔15〕 Bowler LD1 Representational difference analysis of cDNA1 Methods MolMed , 2004 , 94 : 49 - 661〔16〕 Chua KB , Wang LF , Lam SK, et al1 Tioman virus , a novel paramyxo2virus isolated fromfruit bats in Malaysia1Virology , 2001 , 283 : 215 -2291〔17〕 Bowden TR , Westenberg M, Wang LF , et al1 Molecular characteriza2tion of Menangle virus , a novel paramyxovirus which infects pigs , frutbats , and humans1 Virology , 2001 , 283 : 358 - 373〔18〕 Endoh D , Mizatanil T , Kirisawa R , et al1 Species - independent detec2tion of RNA virus by representational difference analysis using non - ri2bosomal hexanncleotides for reverse transcription1 Nucleic Acids Res ,2005 , 33 : e651〔19〕 Siebert PD , Chenchik A , Kellogg DE , et al1 An improved PCR methodfor walking in uncloned genomic DNA1 Nucleic Acids Res , 1995 , 23 :1087 - 10881〔20〕 Diatchenko L , Lau YF , Campbell AP1 Suppression subtractive hy2bridization : a method for generating differentially regulated or tissue -specific cDNA probes and libraries1 Proc Natl Acad Sci U S A , 1996 ,93 : 6025 - 60301〔21〕 Hu Y, Hirshfield I1 Rapid approach to identify an unrecognized viral a2gent1 J Virol Methods , 2005 , 127 : 80 - 861〔22〕 Dean FB , Nelson JR , Giesler TL , et al1 Rapid amlification of plasmidand phage DNA using phi 29 DNA polymerase and multiply - primedrolling circle amplificationg1 Genome Res , 2001 , 11 : 1095 - 10991〔23〕 Rector A , Tachezy R , Ranst MV1 A sequence - independent strategyfor detection and cloning of circular DNA virus genomes by using multi2ply primed rolling - circle amplification1 J Virol , 2004 , 78 : 4993 -49981〔24〕 Allander T , Emerson SU , Engle RE , et al1 A virus discovery methodincorporating DNase treatment and its applicationg to the identificationgof two bovine parvovirus species1 Proc Natl Aced Sci USA , 2001 , 98 :11609 - 116141〔25〕 Stang A , Korn K, Wildner O , et al1 Characterization of virus isolates byparticle - associated nucleic acid PCR1 J Clin Microbiol , 2005 , 43 :716 - 7201(收稿日期: 2006 - 05 - 15)中国预防医学杂志2007 年6 月第8 卷第3 期 Chin Prev Med , June 2007 , Vol18 No13 ·319 ·© 1994-2008 China Academic Journal Electronic Publishing House. All rights reserved.