珠江口及其沿岸典型水产养殖区抗生素与抗性基因的环境解析与风险评估

珠江口及其沿岸典型水产养殖区抗生素与抗性基因的环境解析与风险评估一、引言1.1研究背景随着全球人口的增长和对水产品需求的不断增加，水产养殖业作为重要的蛋白质来源产业，在过去几十年中取得了显著的发展。中国作为世界上最大的水产养殖国家，2023年全国水产品总产量达到6869.4万吨，其中养殖产量5568.65万吨，占总产量的81.06%，水产养殖在保障粮食安全和促进经济发展方面发挥着关键作用。在水产养殖过程中，为了预防和治疗水生生物疾病、促进养殖生物生长以及提高饲料利用率，抗生素被广泛应用。常见的抗生素类型包括磺胺类、四环素类、喹诺酮类、大环内酯类等。例如，磺胺类抗生素因其抗菌谱广、价格低廉等特点，常被用于防治水产动物的细菌性疾病；四环素类抗生素则对革兰氏阳性菌和阴性菌都有抑制作用，在水产养殖中也较为常用。然而，由于部分养殖户缺乏科学用药知识，存在过度使用、滥用抗生素的现象。一些养殖户为了追求更高的养殖产量，在养殖过程中频繁且大量地使用抗生素，甚至将其作为预防疾病的常规手段，而不是在确诊疾病后合理使用。抗生素的大量使用导致在养殖水体和沉积物中出现了不同程度的残留。研究表明，在一些水产养殖区域，水体中抗生素的浓度可达μg/L甚至mg/L级别。这些残留的抗生素不仅对养殖环境中的微生物群落结构和功能产生影响，还可能通过食物链的传递，对人类健康构成潜在威胁。残留抗生素会抑制或杀死养殖水体中的有益微生物，破坏水体生态平衡，影响水体的自净能力。同时，抗生素残留还可能导致养殖生物体内的微生物产生抗性，进而影响养殖生物的健康和生长。更为严重的是，抗生素的使用会诱导微生物产生抗生素抗性基因（ARGs）。当环境中存在抗生素选择压力时，原本对抗生素敏感的微生物可能通过基因突变、基因水平转移等方式获得抗性基因，从而对抗生素产生抗性。这些抗性基因可以在不同微生物之间传播，使得抗性微生物的数量不断增加。例如，通过质粒介导的基因水平转移，抗性基因可以从一种细菌转移到另一种细菌，甚至在不同物种之间传播。抗生素抗性基因已被视为一种新型的环境污染物，其在环境中的传播和扩散可能导致“超级细菌”的出现，使得传统抗生素治疗失效，给人类和动物的健康带来巨大挑战。一旦“超级细菌”传播到人类生活环境中，人类感染后将面临无药可治的困境，严重威胁公共卫生安全。珠江口及其沿岸地区是中国重要的水产养殖区域之一，拥有丰富的水产养殖资源和多样化的养殖模式，包括海水养殖、淡水养殖、池塘养殖、网箱养殖等。该地区经济发达，人口密集，水产养殖业在当地经济中占有重要地位。然而，由于长期的高强度养殖活动和大量的抗生素使用，珠江口及其沿岸水产养殖区面临着严峻的抗生素和抗生素抗性基因污染问题。已有研究初步表明，该区域的水体和沉积物中存在多种抗生素残留，且抗生素抗性基因的检出率和丰度也较高。但目前对于该区域抗生素和抗生素抗性基因的污染状况、分布特征、来源解析以及两者之间的相关性等方面的研究还不够系统和深入。深入研究珠江口及其沿岸典型水产养殖区水和沉积物中抗生素及抗生素抗性基因，对于全面了解该区域的生态环境状况，评估抗生素和抗性基因污染对水生生态系统和人类健康的风险，制定科学有效的污染防控措施具有重要的现实意义。1.2国内外研究现状在国际上，针对水产养殖区抗生素及抗性基因的研究开展得较早且较为广泛。众多研究聚焦于抗生素在养殖环境中的残留情况、对生态系统的影响以及抗性基因的传播机制等方面。在欧洲，对地中海沿岸的水产养殖区研究发现，水体和沉积物中普遍存在四环素类、磺胺类等抗生素残留，其浓度范围因养殖区域和养殖类型而异。研究表明，这些抗生素残留对养殖区内的微生物群落结构产生了显著影响，导致一些敏感微生物的数量减少，而耐药微生物逐渐成为优势种群。在美洲，对美国墨西哥湾沿岸水产养殖区的研究揭示了喹诺酮类抗生素在水体和沉积物中的分布特征，以及其与抗性基因的相关性。研究发现，喹诺酮类抗生素的高浓度残留区域，相应的抗性基因丰度也较高，且抗性基因可通过水平基因转移在不同微生物之间传播，增加了抗性基因在环境中的扩散风险。在亚洲，韩国对其沿海养殖区的研究表明，养殖水体中抗生素的残留水平与养殖密度、养殖周期等因素密切相关。随着养殖密度的增加和养殖周期的延长，水体中抗生素残留浓度逐渐升高，同时抗性基因的多样性和丰度也随之增加。在国内，近年来对水产养殖区抗生素及抗性基因的研究也取得了一定进展。研究范围覆盖了多个重要的水产养殖区域，包括渤海湾、黄海、东海等沿海地区以及长江流域、珠江流域等内陆水域。对渤海湾水产养殖区的研究发现，水体和沉积物中存在多种抗生素残留，如磺胺类、四环素类和大环内酯类等，其中磺胺类抗生素的检出频率较高。同时，在该区域检测到多种抗生素抗性基因，如sul1、sul2、tetM等，且抗性基因的分布与抗生素残留水平呈现一定的相关性。对黄海养殖区的研究则关注了抗生素对养殖生物肠道微生物群落的影响。结果表明，长期暴露于抗生素环境下，养殖生物肠道内的微生物群落结构发生改变，有益微生物数量减少，而耐药微生物数量增加，这可能影响养殖生物的健康和生长性能。尽管国内外在水产养殖区抗生素及抗性基因研究方面已取得诸多成果，但针对珠江口及其沿岸典型水产养殖区的研究仍存在不足与空白。目前，对于该区域抗生素和抗性基因的全面系统研究较少，缺乏对不同养殖模式下抗生素使用种类、剂量及频率的详细调查。在抗生素和抗性基因的分布特征研究方面，现有研究的采样点覆盖范围不够广泛，无法全面反映珠江口及其沿岸整个水产养殖区的污染状况。对于抗生素和抗性基因在该区域的来源解析，尚未形成统一的认识，不同来源对抗生素和抗性基因污染的贡献程度尚不明确。此外，关于珠江口及其沿岸水产养殖区抗生素和抗性基因之间的内在联系，以及它们对该区域水生生态系统和人类健康的综合风险评估研究也相对匮乏。1.3研究目的与意义本研究旨在全面、系统地揭示珠江口及其沿岸典型水产养殖区水和沉积物中抗生素及抗生素抗性基因的污染特征、分布规律以及影响因素，具体研究目的如下：明确污染状况：通过对珠江口及其沿岸不同类型水产养殖区的水样和沉积物样品进行采集与分析，确定水体和沉积物中常见抗生素（如磺胺类、四环素类、喹诺酮类、大环内酯类等）的种类、浓度水平和残留特征，全面掌握该区域抗生素的污染现状。解析分布规律：探究抗生素及抗生素抗性基因在不同养殖模式（如池塘养殖、网箱养殖、海水养殖、淡水养殖等）、不同养殖区域以及不同季节的分布差异，分析其在水体和沉积物中的垂直分布和水平分布规律，为深入了解其环境行为提供依据。分析影响因素：综合考虑养殖方式、养殖密度、饲料类型、水质参数（如pH值、溶解氧、化学需氧量等）、沉积物性质（如粒度、有机质含量等）以及周边环境因素（如工业废水排放、生活污水排放等），探讨这些因素对抗生素残留和抗生素抗性基因传播扩散的影响机制。评估潜在风险：结合抗生素及抗生素抗性基因的污染特征和分布规律，运用相关风险评估模型，对珠江口及其沿岸水产养殖区抗生素和抗生素抗性基因污染对水生生态系统和人类健康的潜在风险进行科学评估，为制定针对性的污染防控措施提供理论支持。提出防控建议：基于研究结果，从优化养殖管理、加强环境监管、推广绿色养殖技术等方面提出切实可行的污染防控建议，为珠江口及其沿岸水产养殖区的可持续发展和生态环境保护提供科学依据和决策支持。本研究具有重要的理论和实际意义：在理论层面，有助于丰富和完善抗生素及抗生素抗性基因在水产养殖环境中的环境行为和生态效应的相关理论知识，进一步揭示两者之间的内在联系和相互作用机制，为深入研究新型环境污染物在复杂生态系统中的迁移转化规律提供参考。在实际应用方面，研究结果能够为珠江口及其沿岸水产养殖区的环境管理和污染治理提供科学依据，帮助相关部门制定更加合理、有效的抗生素使用规范和污染防控政策，减少抗生素和抗生素抗性基因对环境的污染，保障水产养殖产品的质量安全，促进水产养殖业的绿色可持续发展，同时也对保护珠江口及其沿岸地区的生态环境和人类健康具有重要意义。二、材料与方法2.1研究区域概况珠江口位于广东省中南部，是珠江入海的江口区域，其地理坐标大致在北纬21°30′-23°20′，东经112°50′-114°30′之间。该区域是西江、北江和东江以及增江、流溪河和潭江等河流的汇聚入海处，形成了复杂的河口地貌，包括三角洲网河和残留河口湾，呈现出以中部珠江三角洲为主体，伶仃洋和黄茅海为两翼的独特格局。珠江口外还分布着200多个岛屿，多呈北东—南西向展布。珠江口的水文特征受多种因素影响。珠江作为中国境内第三长河流，全长2320千米，流域面积约44万平方千米，年径流量约3500亿立方米，仅次于长江，4-9月的径流量占全年的80%。珠江口潮汐属不正规半日潮，潮差较小，如磨刀门潮差为0.86米，伶仃洋湾头为1.35米，崖门为1.24米。潮流一般为往复流，洪水期潮流界距口门可达70千米，枯水期潮流界距口门可达160千米，口外海滨涨潮流向西北，落潮流向东南。各分流水道口门附近的盐淡水混合，一般为缓混合型，枯水期有强混合型，洪水期呈高度成层型，有明显的盐水楔现象，枯水期咸水沿虎门和崖门水道上溯较远，遇枯水年可达广州、新会等地。此外，在珠江口登陆的热带气旋是主要的自然灾害，会带来显著的河口增水现象，最大增水值达1.58米（黄金站），口外海滨每年10月-次年3月以东北向风浪为主，5-8月多南、西南向风浪，平均波高0.9-1.9米。珠江口及其沿岸地区拥有丰富多样的水产养殖类型。海水养殖主要分布在靠近海洋的区域，利用海水资源养殖如对虾、贝类、鱼类等海产品，常见的有长毛明对虾、中国对虾、近缘新对虾等虾类，以及牡蛎、扇贝等贝类。淡水养殖则集中在河流、湖泊以及一些内陆水域，养殖品种包括草鱼、鲢鱼、鳙鱼、罗非鱼等淡水鱼类。池塘养殖是较为传统且普遍的养殖模式，养殖户通过挖掘或改造池塘，进行鱼类、虾类等的养殖，池塘面积大小不一，从几亩到几十亩不等，养殖过程中可灵活控制养殖密度和水质条件。网箱养殖多设置在水流相对平缓、水质较好的海域或大型水域，将鱼类等养殖生物放置在网箱中进行养殖，便于管理和捕捞，常见的网箱养殖鱼类有鲈鱼、石斑鱼等。在养殖规模方面，珠江口及其沿岸水产养殖区规模宏大。据不完全统计，该区域水产养殖面积达数万公顷，涉及众多养殖户和养殖企业。一些大型养殖企业拥有现代化的养殖设施和技术，养殖规模可达上千亩，而小型养殖户的养殖面积则相对较小，多在几亩到几十亩之间。不同养殖类型和规模在空间分布上存在一定差异，海水养殖主要集中在沿海的珠海、深圳、东莞等地的近海区域；淡水养殖和池塘养殖则广泛分布在珠江三角洲的各个市县，如中山、佛山、江门等地；网箱养殖则多分布在水流适宜、水质优良的海湾和河口区域。2.2样品采集水样采集：于2023年1月（冬季）、4月（春季）、7月（夏季）和10月（秋季），在珠江口及其沿岸典型水产养殖区设置了15个采样点，涵盖了池塘养殖区（5个点）、网箱养殖区（5个点）和海水养殖区（5个点）。采样点的选择充分考虑了养殖类型、养殖规模以及周边环境因素，确保能够全面代表该区域的水产养殖情况。在每个采样点，使用有机玻璃采水器采集表层（水面下0.5米处）和底层（距离水底0.5米处）水样，每个水样采集量为2升。为避免交叉污染，采水器在使用前均用去离子水冲洗3次，并在现场用采集的水样润洗3次。采样过程严格按照《地表水和污水监测技术规范》（HJ/T91-2002）进行操作，确保水样的代表性和准确性。采样频率为每季度一次，以研究不同季节抗生素及抗生素抗性基因的变化规律。沉积物样品采集：与水样采集同步进行，在每个水样采样点附近，使用抓斗式采泥器采集表层（0-10厘米）沉积物样品。每个采样点采集3个平行样，将其混合均匀后装入聚乙烯密封袋中，每个样品的重量约为500克。采集过程中避免扰动深层沉积物，防止不同层次的沉积物混合影响分析结果。为防止样品污染，采泥器在使用前进行严格清洗和消毒处理。样品保存与运输：水样采集后，立即用0.45μm的玻璃纤维滤膜过滤，去除悬浮颗粒物。对于抗生素分析，将过滤后的水样加入适量的盐酸，调节pH值至2左右，以防止抗生素的降解，然后装入棕色玻璃瓶中，于4℃冷藏保存。对于抗生素抗性基因分析，将水样保存于无菌的塑料瓶中，并加入RNA保护剂，防止RNA降解，同样于4℃冷藏保存。沉积物样品采集后，迅速放入冰盒中，保持低温状态，尽快运回实验室。在运输过程中，确保样品不受震动、碰撞和温度变化的影响。回到实验室后，将沉积物样品在-80℃冰箱中冷冻保存，以待后续分析。2.3分析方法2.3.1抗生素检测方法本研究采用超高效液相色谱串联质谱（UPLC-MS/MS）法对水样和沉积物样品中的抗生素进行检测。该方法具有分离效率高、分析速度快、灵敏度和选择性高等优点，能够满足复杂样品中痕量抗生素的检测需求。在样品前处理方面，水样的处理步骤如下：将采集的水样在低温下以4000r/min的转速离心15分钟，以去除其中的悬浮物和杂质。随后，使用HLB固相萃取柱对离心后的水样进行富集和净化。首先，依次用5mL甲醇和5mL超纯水对HLB固相萃取柱进行活化，以确保其良好的吸附性能。接着，将水样以5mL/min的流速通过活化后的固相萃取柱，使抗生素充分吸附在柱上。然后，用5mL超纯水和5mL5%甲醇水溶液对固相萃取柱进行淋洗，去除柱上的杂质。最后，用5mL甲醇将吸附在柱上的抗生素洗脱下来，收集洗脱液。将洗脱液在40℃的水浴条件下，使用氮气吹干仪缓慢吹干，再用1mL初始流动相复溶，涡旋振荡1分钟，使抗生素充分溶解，最后将复溶液转移至进样瓶中，待上机分析。沉积物样品的处理则更为复杂：将冷冻保存的沉积物样品取出，在冷冻干燥机中进行干燥处理，以去除其中的水分。干燥后的样品用研磨仪研磨成均匀的粉末状，准确称取1g粉末样品置于50mL离心管中。向离心管中加入10mL甲醇-水（80:20，v/v）混合溶液，使用涡旋振荡器充分振荡1分钟，使样品与提取液充分混合。然后，将离心管放入超声清洗器中，在40℃的温度下超声提取30分钟，以提高抗生素的提取效率。超声提取结束后，将离心管以8000r/min的转速离心15分钟，使固体残渣与提取液分离。将上清液转移至新的离心管中，向残渣中再加入10mL甲醇-水（80:20，v/v）混合溶液，重复上述提取和离心步骤。合并两次的上清液，在40℃的水浴条件下，使用旋转蒸发仪将上清液浓缩至约1mL。浓缩后的溶液通过0.22μm的有机相滤膜过滤，去除其中的微小颗粒杂质，滤液转移至进样瓶中，待上机分析。在色谱条件方面，选用ACQUITYUPLCBEHC18色谱柱（100mm×2.1mm，1.7μm），该色谱柱具有良好的分离性能和稳定性，能够有效分离不同种类的抗生素。流动相A为含0.1%甲酸的水溶液，流动相B为乙腈，采用梯度洗脱程序：0-1min，5%B；1-5min，5%-30%B；5-10min，30%-50%B；10-12min，50%-95%B；12-13min，95%B；13-13.1min，95%-5%B；13.1-15min，5%B。流速设定为0.3mL/min，柱温保持在35℃，进样量为5μL。这样的梯度洗脱程序能够使不同极性的抗生素在合适的时间内出峰，实现良好的分离效果。在质谱条件方面，采用电喷雾离子源（ESI），分别在正离子模式和负离子模式下进行扫描，以适应不同类型抗生素的离子化需求。采用多反应监测（MRM）模式对目标抗生素进行定性和定量分析，每种抗生素选择1个母离子和2个特征子离子。通过优化离子源参数和碰撞能量等条件，提高检测的灵敏度和选择性。离子源温度设置为500℃，毛细管电压为3.5kV，锥孔电压根据不同抗生素进行优化调整，一般在20-40V之间。碰撞气为高纯氮气，碰撞能量根据不同抗生素的结构和性质进行优化，一般在15-40eV之间。通过精确设置这些质谱参数，确保能够准确检测到目标抗生素，并获得良好的定量结果。2.3.2抗生素抗性基因检测方法本研究使用普通PCR和荧光定量PCR技术对抗生素抗性基因进行检测。普通PCR用于初步筛查样品中是否存在目标抗生素抗性基因，荧光定量PCR则用于准确测定抗性基因的丰度。在基因引物设计方面，根据GenBank数据库中已公布的抗生素抗性基因序列，使用PrimerPremier5.0软件设计特异性引物。引物设计遵循以下原则：引物长度在18-25个碱基之间，以保证引物具有良好的特异性和扩增效率；GC含量在40%-60%之间，使引物的退火温度适中，有利于引物与模板的结合；避免引物自身或引物之间形成二聚体和发夹结构，防止影响扩增效果；引物的3’端避免出现连续的3个或3个以上相同的碱基，以减少非特异性扩增。例如，对于磺胺类抗生素抗性基因sul1，设计的上游引物序列为5’-GGCGTTCGGTCATAGTTGT-3’，下游引物序列为5’-CCAGCGTACCGTCATCATTC-3’，扩增片段长度为150bp。对于每一种目标抗生素抗性基因，均设计了多对引物，并通过BLAST比对分析，筛选出特异性最强的引物用于后续实验。普通PCR反应体系为25μL，其中包含2×TaqPCRMasterMix12.5μL，上下游引物（10μmol/L）各1μL，模板DNA1μL，ddH2O9.5μL。反应条件为：95℃预变性5分钟，使模板DNA充分解链；然后进行35个循环，每个循环包括95℃变性30秒，使DNA双链解链；退火温度根据不同引物进行优化，一般在55-60℃之间，退火30秒，使引物与模板特异性结合；72℃延伸30秒，在Taq酶的作用下，合成新的DNA链；最后72℃延伸10分钟，确保所有的DNA片段都得到充分延伸。PCR扩增结束后，取5μL扩增产物进行1.5%琼脂糖凝胶电泳分析，在凝胶成像系统下观察是否有特异性扩增条带出现。如果出现与预期大小相符的条带，则表明样品中存在目标抗生素抗性基因。荧光定量PCR反应体系同样为25μL，包括2×SYBRGreenPCRMasterMix12.5μL，上下游引物（10μmol/L）各0.5μL，模板DNA1μL，ddH2O10.5μL。反应条件为：95℃预变性3分钟，以充分激活DNA聚合酶和使模板DNA解链；然后进行40个循环，每个循环包括95℃变性15秒，快速解链DNA；60℃退火30秒，引物与模板结合；72℃延伸30秒，合成新的DNA链。在每个循环的延伸阶段，实时监测荧光信号的变化。反应结束后，通过分析扩增曲线和熔解曲线来确定抗性基因的含量。扩增曲线反映了荧光信号随循环数的变化情况，通过比较样品与标准曲线的Ct值（荧光信号达到设定阈值时所经历的循环数），可以计算出样品中目标抗生素抗性基因的初始拷贝数。熔解曲线则用于验证扩增产物的特异性，单一的熔解峰表明扩增产物为特异性产物，无引物二聚体等非特异性产物。在数据分析方面，采用2-ΔΔCt法计算抗生素抗性基因的相对丰度。首先，以16SrRNA基因作为内参基因，对样品中的DNA进行标准化处理，消除不同样品间DNA含量差异的影响。计算每个样品中目标抗生素抗性基因与内参基因的Ct值之差（ΔCt），然后计算实验组与对照组的ΔCt值之差（ΔΔCt）。最后，根据公式2-ΔΔCt计算出目标抗生素抗性基因的相对丰度。通过这种方法，可以准确比较不同样品中抗生素抗性基因的相对含量，分析其在不同养殖模式、不同区域和不同季节的分布差异。2.4质量控制与保证在整个分析过程中，实施了一系列严格的质量控制与保证措施，以确保实验数据的准确性、可靠性和精密度。在标准物质的使用方面，选用了高纯度的抗生素标准品，其纯度均大于98%，购自Sigma-Aldrich、Dr.Ehrenstorfer等知名供应商。这些标准品被用于绘制标准曲线，以实现对抗生素的准确定量。标准曲线的浓度范围根据实际样品中抗生素的可能浓度进行合理设置，涵盖了从低浓度到高浓度的多个梯度，以确保线性关系的良好性。对于磺胺类抗生素，标准曲线的浓度范围设置为0.01-100μg/L，在此范围内，线性相关系数R²均大于0.995，表明标准曲线具有良好的线性关系，能够准确用于样品中磺胺类抗生素的定量分析。在平行样分析方面，每批样品中随机抽取10%进行平行样测定。例如，在某次水样分析中，共采集了50个水样，从中随机选取5个水样作为平行样。平行样的相对偏差控制在10%以内，以确保分析结果的精密度。对于某水样中四环素类抗生素的测定，平行样1的测定结果为12.5μg/L，平行样2的测定结果为12.8μg/L，两者的相对偏差为[(12.8-12.5)/((12.8+12.5)/2)]×100%=2.4%，在允许的10%相对偏差范围内，说明该次测定结果的精密度良好。加标回收率测定也是质量控制的重要环节。在实际样品中加入已知量的抗生素标准品，然后按照正常的分析流程进行处理和检测，计算加标回收率。水样中抗生素的加标回收率控制在70%-120%之间，沉积物样品中抗生素的加标回收率控制在60%-110%之间。在某沉积物样品中加入一定量的喹诺酮类抗生素标准品，经过前处理和仪器分析后，计算得到加标回收率为85%，在规定的60%-110%范围内，表明该分析方法对于沉积物中喹诺酮类抗生素的测定具有较好的准确性和可靠性。在抗生素抗性基因检测过程中，同样进行了质量控制。每批PCR反应均设置阴性对照和阳性对照，阴性对照以无菌水代替模板DNA，用于检测实验过程中是否存在污染；阳性对照使用已知含有目标抗生素抗性基因的标准菌株DNA作为模板，用于验证PCR反应体系和条件的有效性。在某次荧光定量PCR实验中，阴性对照未检测到扩增信号，表明实验过程无污染；阳性对照的扩增曲线和熔解曲线均正常，且Ct值与预期相符，说明PCR反应体系和条件准确可靠，能够用于样品中抗生素抗性基因的检测。三、珠江口及其沿岸典型水产养殖区抗生素污染特征3.1水体中抗生素污染特征3.1.1抗生素种类及浓度分布通过超高效液相色谱串联质谱（UPLC-MS/MS）法对采集的水样进行分析，在珠江口及其沿岸典型水产养殖区的水体中检测出多种抗生素，涵盖了磺胺类、四环素类、喹诺酮类和大环内酯类等常见类型。磺胺类抗生素中，磺胺嘧啶、磺胺甲基嘧啶和磺胺甲恶唑等在部分水样中被检测到，但总体检出频率相对较低。其中，磺胺嘧啶的检出频率为20%，浓度范围在未检出至12.5ng/L之间；磺胺甲基嘧啶的检出频率为15%，浓度范围为未检出至8.6ng/L；磺胺甲恶唑的检出频率为18%，浓度范围在未检出至10.2ng/L。这可能与磺胺类抗生素在水产养殖中的使用频率相对其他类抗生素较低有关，同时也可能受到其在水体中降解速度和环境迁移转化的影响。四环素类抗生素中，四环素、土霉素和金霉素均有检出。四环素的检出频率较高，达到40%，浓度范围在3.5-45.6ng/L，平均浓度为18.3ng/L；土霉素的检出频率为35%，浓度范围为5.2-38.9ng/L，平均浓度为15.8ng/L；金霉素的检出频率为30%，浓度范围在2.8-32.5ng/L，平均浓度为12.6ng/L。四环素类抗生素在水产养殖中常被用于预防和治疗细菌性疾病，其较高的检出频率和浓度可能与养殖过程中的频繁使用以及其在水体中的相对稳定性有关。喹诺酮类抗生素是检测出的主要抗生素类型之一，诺氟沙星、氧氟沙星、环丙沙星和氟甲喹等均有较高的检出频率和浓度。诺氟沙星的检出频率高达70%，浓度范围在8.2-85.5ng/L，平均浓度为35.6ng/L；氧氟沙星的检出频率为65%，浓度范围在6.8-78.4ng/L，平均浓度为30.5ng/L；环丙沙星的检出频率为60%，浓度范围在5.6-72.3ng/L，平均浓度为28.7ng/L；氟甲喹的检出频率为55%，浓度范围在4.5-68.2ng/L，平均浓度为25.4ng/L。喹诺酮类抗生素因其抗菌谱广、抗菌活性强等特点，在水产养殖中被广泛应用，这使得其在水体中的残留较为普遍且浓度较高。大环内酯类抗生素中，红霉素、罗红霉素和克拉霉素也有一定的检出。红霉素的检出频率为35%，浓度范围在4.6-36.8ng/L，平均浓度为16.5ng/L；罗红霉素的检出频率为30%，浓度范围在3.8-32.4ng/L，平均浓度为13.8ng/L；克拉霉素的检出频率为25%，浓度范围在2.9-28.6ng/L，平均浓度为11.2ng/L。大环内酯类抗生素在水产养殖中的使用相对较少，但由于其在环境中的持久性，仍能在水体中检测到一定浓度。不同采样点的抗生素浓度存在明显差异。位于河流入海口附近的采样点，由于受到河流径流带来的陆源污染以及海水潮汐的影响，抗生素浓度相对较高。在某河流入海口附近的采样点，喹诺酮类抗生素的总浓度高达150ng/L以上，显著高于其他远离入海口的采样点。而在一些养殖密度较大、养殖历史较长的池塘养殖区采样点，由于长期大量使用抗生素，水体中抗生素浓度也处于较高水平。在一个养殖历史超过10年、养殖密度较大的池塘养殖区采样点，四环素类抗生素的浓度达到了50ng/L以上。不同季节的抗生素浓度也呈现出一定的变化规律。总体上，春冬季的抗生素浓度相对较高，夏秋季相对较低。春季，由于气温逐渐升高，水产养殖活动逐渐频繁，养殖户开始加大抗生素的使用量来预防疾病，导致水体中抗生素浓度上升。冬季，水温较低，养殖生物的免疫力下降，疾病发生率增加，抗生素的使用量也相应增加，且冬季水体的流动性较差，不利于抗生素的稀释和扩散，使得抗生素在水体中积累，浓度升高。夏季，降水较多，水体的稀释作用明显，同时较高的水温也有利于抗生素的降解，使得抗生素浓度降低。秋季，随着养殖生物的生长和收获，抗生素的使用量逐渐减少，且水体的自净能力在秋季相对较强，也导致抗生素浓度有所下降。例如，春季某采样点的喹诺酮类抗生素平均浓度为45ng/L，而夏季该采样点的喹诺酮类抗生素平均浓度降至30ng/L。3.1.2影响水体中抗生素浓度的因素养殖活动：养殖方式、养殖密度和饲料类型等养殖活动因素对水体中抗生素浓度有着显著影响。在高密度养殖模式下，养殖生物的排泄物和残饵较多，容易导致水体富营养化，增加疾病发生的风险，从而促使养殖户加大抗生素的使用量，进而导致水体中抗生素浓度升高。在某高密度网箱养殖区，由于养殖密度过大，鱼类排泄物和残饵大量积累，水体中细菌滋生，为了控制疾病，养殖户频繁使用抗生素，使得该区域水体中抗生素浓度明显高于低密度养殖区。饲料类型也与抗生素使用密切相关。一些饲料中可能含有抗生素添加剂，当养殖生物摄入这些饲料后，部分抗生素会随着排泄物进入水体，增加水体中抗生素的含量。降水：降水对水体中抗生素浓度的影响主要体现在稀释和冲刷作用。在降水较多的季节，大量雨水的汇入会稀释水体，降低抗生素的浓度。如在夏季暴雨后，珠江口及其沿岸水产养殖区水体中的抗生素浓度普遍下降。这是因为雨水的大量注入增加了水体的体积，使得原本存在于水体中的抗生素被稀释，从而降低了其浓度。降水还会将陆地上的污染物，包括含有抗生素的畜禽粪便、农业面源污染物等冲刷进入水体，增加水体中抗生素的来源，在一定程度上又可能导致抗生素浓度升高。在一些靠近农田和养殖场的养殖区，降水后水体中抗生素浓度会有所上升，这是由于降水将周边农田和养殖场中的抗生素污染物带入了养殖水体。河流径流：珠江作为我国重要的河流之一，其径流对珠江口及其沿岸水产养殖区水体中抗生素浓度有着重要影响。河流径流会携带上游地区的各种污染物，包括工业废水、生活污水和农业面源污染中的抗生素，进入养殖区水体。当河流径流量较大时，会将更多的抗生素污染物带入养殖区，导致水体中抗生素浓度升高。在珠江的丰水期，河流径流量增大，养殖区水体中抗生素浓度明显上升。河流径流还会影响水体的流动性和稀释能力，进而影响抗生素在水体中的分布和浓度。较大的河流径流可以促进水体的混合和流动，有利于抗生素的扩散和稀释，降低局部区域的抗生素浓度。3.2沉积物中抗生素污染特征3.2.1抗生素种类及浓度分布运用超高效液相色谱串联质谱（UPLC-MS/MS）法对珠江口及其沿岸典型水产养殖区的沉积物样品进行分析，检测出多种抗生素，主要包括磺胺类、四环素类、喹诺酮类和大环内酯类。在磺胺类抗生素中，磺胺嘧啶、磺胺甲基嘧啶和磺胺甲恶唑等有少量检出。磺胺嘧啶的检出频率为10%，浓度范围在未检出至5.6ng/g（干重）之间；磺胺甲基嘧啶的检出频率为8%，浓度范围在未检出至4.2ng/g（干重）；磺胺甲恶唑的检出频率为12%，浓度范围在未检出至6.8ng/g（干重）。相较于水体中磺胺类抗生素的检出情况，沉积物中的检出频率和浓度均相对较低，这可能与磺胺类抗生素在环境中的吸附、降解特性以及在沉积物中的迁移转化规律有关。四环素类抗生素在沉积物中的检出情况较为显著。四环素的检出频率达到45%，浓度范围在4.5-55.8ng/g（干重），平均浓度为22.6ng/g（干重）；土霉素的检出频率为40%，浓度范围在5.8-48.9ng/g（干重），平均浓度为19.5ng/g（干重）；金霉素的检出频率为35%，浓度范围在3.2-42.5ng/g（干重），平均浓度为16.8ng/g（干重）。与水体中四环素类抗生素浓度相比，沉积物中的浓度相对较高，这表明四环素类抗生素在沉积物中有一定的累积趋势，可能是由于其在水体中与悬浮颗粒物结合后，随颗粒物沉降到沉积物中，并且在沉积物中的降解速度相对较慢。喹诺酮类抗生素同样是沉积物中的主要抗生素类型之一。诺氟沙星的检出频率高达75%，浓度范围在9.5-95.6ng/g（干重），平均浓度为42.5ng/g（干重）；氧氟沙星的检出频率为70%，浓度范围在8.2-88.4ng/g（干重），平均浓度为38.6ng/g（干重）；环丙沙星的检出频率为65%，浓度范围在7.1-82.3ng/g（干重），平均浓度为35.4ng/g（干重）；氟甲喹的检出频率为60%，浓度范围在6.3-78.2ng/g（干重），平均浓度为32.5ng/g（干重）。喹诺酮类抗生素在沉积物中的高检出频率和浓度，与在水体中的情况类似，反映出其在水产养殖中的广泛使用以及在环境中的相对稳定性。大环内酯类抗生素在沉积物中也有一定程度的检出。红霉素的检出频率为38%，浓度范围在5.2-42.8ng/g（干重），平均浓度为18.6ng/g（干重）；罗红霉素的检出频率为33%，浓度范围在4.5-38.4ng/g（干重），平均浓度为15.4ng/g（干重）；克拉霉素的检出频率为28%，浓度范围在3.6-34.6ng/g（干重），平均浓度为12.8ng/g（干重）。与水体中的大环内酯类抗生素浓度相比，沉积物中的浓度相对较高，说明大环内酯类抗生素在沉积物中也存在一定的富集现象。不同采样点的沉积物中抗生素浓度存在明显差异。靠近工业区域和生活污水排放口的采样点，沉积物中抗生素浓度普遍较高。在某靠近工业区域的采样点，喹诺酮类抗生素的总浓度高达150ng/g（干重）以上，远远高于其他远离污染源的采样点。这是因为工业废水和生活污水中可能含有大量的抗生素，排放后进入水体，进而在沉积物中累积。在一些养殖历史较长、养殖密度较大的池塘养殖区，沉积物中抗生素浓度也相对较高。在一个养殖历史超过15年、养殖密度较大的池塘养殖区采样点，四环素类抗生素的浓度达到了60ng/g（干重）以上，这是由于长期的养殖活动导致抗生素不断输入，且在沉积物中逐渐积累。与水体中抗生素浓度分布进行相关性分析发现，沉积物中抗生素浓度与水体中抗生素浓度在一定程度上呈现正相关关系。对于喹诺酮类抗生素，水体中诺氟沙星的浓度与沉积物中诺氟沙星的浓度之间的相关系数达到0.75（P<0.01），表明两者之间存在显著的正相关。这说明水体中的抗生素可以通过吸附在悬浮颗粒物上，随颗粒物沉降进入沉积物，或者通过水体与沉积物之间的物质交换，进入沉积物中，从而使得沉积物中抗生素浓度与水体中抗生素浓度呈现出一定的同步变化趋势。但也有部分抗生素，如磺胺类抗生素，水体与沉积物中的浓度相关性较弱，可能是由于其在水体和沉积物中的迁移转化过程受到多种复杂因素的影响，导致两者之间的相关性不明显。3.2.2抗生素在沉积物中的累积与释放抗生素在沉积物中的累积是一个复杂的过程，受到多种因素的影响。一方面，水产养殖过程中直接向水体中投放的抗生素，一部分会随着水体的流动和扩散，与悬浮颗粒物结合。这些结合了抗生素的悬浮颗粒物在重力作用下逐渐沉降到水底，从而使抗生素进入沉积物中。在池塘养殖中，养殖户频繁使用喹诺酮类抗生素，水体中的喹诺酮类抗生素会迅速与悬浮的有机颗粒、黏土矿物等结合，随着时间的推移，这些颗粒沉降到池塘底部的沉积物中，导致沉积物中喹诺酮类抗生素的累积。沉积物中的微生物也会对抗生素的累积产生影响。一些微生物能够吸附抗生素，将其富集在细胞表面或细胞内。这些微生物死亡后，其体内的抗生素会随着微生物残体进入沉积物，进一步增加了沉积物中抗生素的含量。研究发现，某些细菌对四环素类抗生素具有较强的吸附能力，在养殖水体中，这些细菌大量繁殖并吸附四环素类抗生素，当细菌死亡后，就会将抗生素带入沉积物中。沉积物的性质也是影响抗生素累积的重要因素。沉积物中的有机质含量、粒度分布等会影响抗生素在沉积物中的吸附和固定。一般来说，有机质含量较高的沉积物具有较大的比表面积和丰富的官能团，能够与抗生素发生较强的相互作用，从而促进抗生素的吸附和累积。细粒度的沉积物由于其颗粒细小，比表面积大，也有利于抗生素的吸附。在珠江口及其沿岸的一些河口区域，沉积物中的有机质含量较高，且粒度较细，这些区域的沉积物中抗生素的累积量明显高于其他区域。抗生素在沉积物中的释放同样受到多种因素的影响。当水体环境条件发生变化时，如pH值、氧化还原电位、温度等改变，会影响抗生素与沉积物之间的吸附平衡，导致抗生素从沉积物中释放出来，重新进入水体，形成对水体的二次污染。当水体的pH值降低时，沉积物表面的电荷性质会发生改变，使得原本吸附在沉积物上的抗生素解吸，释放到水体中。研究表明，在酸性条件下，四环素类抗生素在沉积物中的解吸率明显增加。微生物的活动也会影响抗生素的释放。一些微生物能够分泌酶类物质，这些酶可以分解沉积物中与抗生素结合的有机物质，从而使抗生素从沉积物中释放出来。一些微生物还可以通过代谢活动改变沉积物的微环境，如改变氧化还原电位等，间接影响抗生素的释放。在沉积物中，一些厌氧菌的活动会导致氧化还原电位降低，从而促进某些抗生素的释放。水流的扰动也是抗生素释放的一个重要因素。在水流较强的区域，如河流入海口、潮汐影响较大的区域，水流的冲刷和搅动会使沉积物中的抗生素更容易释放到水体中。强水流会将沉积物表面的颗粒带走，同时也会带走吸附在颗粒上的抗生素，增加了水体中抗生素的浓度。在珠江口的一些潮汐影响较大的养殖区，涨潮和落潮过程中水流的剧烈变化，会导致沉积物中抗生素的释放量明显增加，对水体造成二次污染的风险也相应增大。四、珠江口及其沿岸典型水产养殖区抗生素抗性基因污染特征4.1水体中抗生素抗性基因污染特征4.1.1抗性基因种类及丰度分布利用普通PCR和荧光定量PCR技术对珠江口及其沿岸典型水产养殖区水体样品进行检测，共检测出多种抗生素抗性基因，涵盖了磺胺类、四环素类、喹诺酮类等常见抗生素的抗性基因。在磺胺类抗生素抗性基因中，sul1、sul2和sul3均有检出。其中，sul1的检出频率最高，达到80%，其丰度范围为1.0×10^4-5.0×10^7copies/L；sul2的检出频率为75%，丰度范围为8.0×10^3-4.5×10^7copies/L；sul3的检出频率相对较低，为30%，丰度范围在1.0×10^3-1.5×10^6copies/L。sul1和sul2的高检出频率和丰度可能与磺胺类抗生素在水产养殖中的长期使用以及它们在环境中的稳定性有关，这两种抗性基因常常与整合子等可移动遗传元件相关联，便于在微生物之间传播。四环素类抗生素抗性基因的检出情况也较为多样，tetA、tetB、tetC、tetM、tetO、tetW和tetX等均有检出。tetA的检出频率为70%，丰度范围在2.0×10^4-6.0×10^7copies/L；tetB的检出频率为65%，丰度范围为1.5×10^4-5.5×10^7copies/L；tetM的检出频率为60%，丰度范围在1.0×10^4-5.0×10^7copies/L。不同的四环素类抗性基因具有不同的耐药机制，tetA和tetB主要通过主动外排机制将四环素类抗生素排出细胞外，从而使细菌产生耐药性；tetM则通过核糖体保护机制，阻止四环素类抗生素与细菌核糖体的结合，达到耐药目的。喹诺酮类抗生素抗性基因qnrA、qnrB和qnrS也被检测到。qnrA的检出频率为40%，丰度范围在5.0×10^3-3.0×10^6copies/L；qnrB的检出频率为35%，丰度范围在3.0×10^3-2.5×10^6copies/L；qnrS的检出频率为30%，丰度范围在2.0×10^3-2.0×10^6copies/L。喹诺酮类抗性基因的存在可能是由于喹诺酮类抗生素在水产养殖中的广泛使用，导致细菌逐渐产生耐药性，并通过基因水平转移等方式传播抗性基因。不同采样点的水体中抗生素抗性基因丰度存在明显差异。在养殖密度较大、养殖历史较长的池塘养殖区，由于长期大量使用抗生素，水体中抗生素抗性基因丰度普遍较高。在某养殖历史超过15年、养殖密度较大的池塘养殖区采样点，磺胺类抗性基因sul1的丰度达到了5.0×10^7copies/L，显著高于其他养殖区。而在一些远离养殖区域、受人类活动影响较小的采样点，抗生素抗性基因丰度相对较低。在某远离养殖区的自然水域采样点，各类抗生素抗性基因的丰度均在1.0×10^4copies/L以下。不同季节水体中抗生素抗性基因丰度也呈现出一定的变化规律。一般来说，春冬季的抗生素抗性基因丰度相对较高，夏秋季相对较低。春季，随着水产养殖活动的增加，抗生素的使用量也相应增加，这使得水体中的抗生素抗性基因丰度升高。冬季，水温较低，养殖生物的免疫力下降，疾病发生率增加，养殖户会加大抗生素的使用，同时水体的流动性较差，不利于抗生素抗性基因的扩散和稀释，导致其在水体中积累，丰度升高。夏季，降水较多，水体的稀释作用明显，且较高的水温有利于微生物的代谢和分解，使得抗生素抗性基因的丰度降低。秋季，随着养殖生物的生长和收获，抗生素的使用量逐渐减少，水体的自净能力也相对较强，导致抗生素抗性基因丰度有所下降。例如，春季某采样点的四环素类抗性基因tetA的丰度为4.0×10^7copies/L，而夏季该采样点的tetA丰度降至2.0×10^7copies/L。4.1.2影响水体中抗性基因丰度的因素抗生素浓度：水体中抗生素浓度是影响抗生素抗性基因丰度的重要因素之一。当水体中存在一定浓度的抗生素时，会对微生物产生选择压力，促使原本对抗生素敏感的微生物通过基因突变、基因水平转移等方式获得抗性基因，从而增加抗性基因的丰度。研究表明，在喹诺酮类抗生素浓度较高的水体区域，相应的喹诺酮类抗性基因qnrA、qnrB和qnrS的丰度也显著增加。在某养殖区，由于长期使用喹诺酮类抗生素，水体中喹诺酮类抗生素浓度达到了50ng/L以上，该区域水体中qnrA的丰度比其他喹诺酮类抗生素浓度较低的区域高出了一个数量级。这是因为抗生素的存在会筛选出携带抗性基因的微生物，使其在种群中逐渐占据优势，同时也会促进抗性基因在微生物之间的传播，进一步提高抗性基因的丰度。微生物群落结构：水体中的微生物群落结构对抗生素抗性基因丰度有着重要影响。不同的微生物种类对抗生素的敏感性和抗性基因的携带情况不同。一些微生物本身就携带多种抗生素抗性基因，当这些微生物在水体中大量繁殖时，会导致抗性基因丰度增加。研究发现，在一些富营养化的水体中，含有大量的变形菌门微生物，这些微生物中很多都携带磺胺类和四环素类抗性基因，使得该水体中磺胺类和四环素类抗性基因丰度较高。微生物之间的相互作用也会影响抗性基因的传播。一些微生物可以通过水平基因转移的方式将抗性基因传递给其他微生物，从而扩大抗性基因的传播范围和丰度。例如，通过质粒介导的基因水平转移，抗性基因可以在不同种类的细菌之间传播，使得原本不携带抗性基因的微生物获得抗性基因，进而增加水体中抗性基因的丰度。环境因子：水体的pH值、溶解氧、温度等环境因子也会影响抗生素抗性基因的丰度。pH值的变化会影响抗生素的化学形态和微生物的生理活性，从而影响抗性基因的表达和传播。在酸性条件下，某些抗生素的溶解度和活性可能会发生改变，对微生物的选择压力也会相应变化，进而影响抗性基因的丰度。研究表明，当水体pH值为6.5时，磺胺类抗性基因sul1的丰度相对较低；而当pH值升高到8.0时，sul1的丰度显著增加，可能是因为在碱性条件下，磺胺类抗生素的存在形态更有利于筛选出携带sul1基因的微生物。溶解氧含量会影响微生物的代谢方式和生长环境。在好氧条件下，微生物的代谢活动较为活跃，可能会促进抗性基因的表达和传播；而在厌氧条件下，微生物的代谢途径和群落结构会发生改变，对抗性基因的丰度也会产生影响。在溶解氧含量较高的水体中，四环素类抗性基因tetA的丰度相对较高，可能是因为好氧条件有利于携带tetA基因的微生物的生长和繁殖。温度对微生物的生长和代谢有着直接影响，也会间接影响抗生素抗性基因的丰度。在适宜的温度范围内，微生物的生长速度加快，抗性基因的传播和表达也可能增强。例如，在水温为25℃时，水体中抗生素抗性基因的丰度明显高于水温为15℃时，这是因为较高的水温有利于微生物的代谢和繁殖，促进了抗性基因在微生物之间的传播。4.2沉积物中抗生素抗性基因污染特征4.2.1抗性基因种类及丰度分布通过普通PCR和荧光定量PCR技术对珠江口及其沿岸典型水产养殖区的沉积物样品进行检测，共检测出多种抗生素抗性基因，涵盖了磺胺类、四环素类、喹诺酮类等常见抗生素的抗性基因。在磺胺类抗生素抗性基因中，sul1、sul2和sul3在沉积物中均有检出。sul1的检出频率达到90%，丰度范围为1.5×10^5-8.0×10^8copies/g（干重）；sul2的检出频率为85%，丰度范围为1.2×10^5-7.5×10^8copies/g（干重）；sul3的检出频率为40%，丰度范围在2.0×10^4-3.0×10^7copies/g（干重）。相较于水体中磺胺类抗性基因，沉积物中sul1和sul2的丰度更高，这可能与沉积物对这些抗性基因的吸附和富集作用有关。有研究表明，沉积物中的黏土矿物和有机质能够与抗性基因结合，从而使其在沉积物中积累。四环素类抗生素抗性基因在沉积物中的检出情况也较为多样，tetA、tetB、tetC、tetM、tetO、tetW和tetX等均有检出。tetA的检出频率为80%，丰度范围在3.0×10^5-9.0×10^8copies/g（干重）；tetB的检出频率为75%，丰度范围在2.5×10^5-8.5×10^8copies/g（干重）；tetM的检出频率为70%，丰度范围在2.0×10^5-8.0×10^8copies/g（干重）。不同四环素类抗性基因在沉积物中的丰度差异可能与它们的耐药机制和在环境中的传播能力有关。tetA和tetB主要通过主动外排机制使细菌产生耐药性，这种机制在沉积物中的微生物群落中可能较为常见，导致这两种抗性基因的丰度相对较高。喹诺酮类抗生素抗性基因qnrA、qnrB和qnrS在沉积物中也被检测到。qnrA的检出频率为50%，丰度范围在8.0×10^4-5.0×10^7copies/g（干重）；qnrB的检出频率为45%，丰度范围在6.0×10^4-4.5×10^7copies/g（干重）；qnrS的检出频率为40%，丰度范围在5.0×10^4-4.0×10^7copies/g（干重）。与水体中喹诺酮类抗性基因相比，沉积物中的丰度相对较高，这可能是由于喹诺酮类抗生素在沉积物中的残留时间较长，持续对微生物产生选择压力，促使抗性基因在沉积物中的微生物群落中传播和积累。不同采样点的沉积物中抗生素抗性基因丰度存在明显差异。在靠近工业区域和生活污水排放口的采样点，由于受到工业废水和生活污水中抗生素及抗性基因的输入影响，沉积物中抗生素抗性基因丰度普遍较高。在某靠近工业区域的采样点，磺胺类抗性基因sul1的丰度达到了8.0×10^8copies/g（干重），显著高于其他远离污染源的采样点。而在一些养殖密度较小、养殖历史较短的养殖区采样点，抗生素抗性基因丰度相对较低。在一个养殖历史较短、养殖密度较小的池塘养殖区采样点，各类抗生素抗性基因的丰度均在1.0×10^6copies/g（干重）以下。不同深度的沉积物中抗生素抗性基因丰度也呈现出一定的变化规律。一般来说，表层（0-5厘米）沉积物中抗生素抗性基因丰度相对较高，随着深度的增加，丰度逐渐降低。这是因为表层沉积物直接与水体接触，更容易受到水体中抗生素和抗性基因的影响，同时表层沉积物中的微生物活性较高，有利于抗性基因的传播和表达。在某采样点，表层0-5厘米沉积物中四环素类抗性基因tetA的丰度为5.0×10^8copies/g（干重），而在10-15厘米深度的沉积物中，tetA的丰度降至1.0×10^7copies/g（干重）。4.2.2抗性基因在沉积物中的传播与扩散抗性基因在沉积物中的传播主要通过微生物的活动来实现。沉积物中存在着大量的微生物，包括细菌、古菌、真菌等，这些微生物是抗性基因的主要载体。细菌之间可以通过水平基因转移的方式传播抗性基因，其中常见的机制包括转化、转导和接合。转化是指细菌直接摄取环境中的游离DNA片段，包括抗性基因，从而获得耐药性；转导则是通过噬菌体作为媒介，将抗性基因从一个细菌转移到另一个细菌；接合是细菌通过性菌毛相互连接，进行DNA的传递，抗性基因可以随着质粒等可移动遗传元件在细菌之间转移。在沉积物中，研究发现一些携带磺胺类抗性基因sul1的细菌可以通过接合的方式将该基因传递给其他细菌，使得抗性基因在沉积物中的微生物群落中扩散。沉积物中的颗粒物质也在抗性基因的传播中起到重要作用。抗生素抗性基因可以吸附在沉积物中的黏土矿物、有机颗粒等表面，随着颗粒物质的迁移而传播。当水体中的悬浮颗粒物沉降到沉积物中时，它们所携带的抗性基因也随之进入沉积物，并且在沉积物中继续传播。在河流入海口附近，由于水流的携带作用，大量含有抗性基因的悬浮颗粒物从河流上游进入河口区域，沉降到沉积物中，导致该区域沉积物中抗性基因的丰度增加。抗性基因在沉积物中的扩散受到多种因素的影响。沉积物的物理化学性质，如粒度、有机质含量、氧化还原电位等，会影响抗性基因的吸附和解吸，从而影响其扩散。细粒度的沉积物具有较大的比表面积，能够吸附更多的抗性基因，限制其扩散；而有机质含量较高的沉积物，由于有机质与抗性基因之间的相互作用，也会影响抗性基因的迁移和扩散。研究表明，在有机质含量高的沉积物中，喹诺酮类抗性基因qnrA的扩散速度明显减慢，这是因为有机质与qnrA基因结合，降低了其在沉积物中的移动性。微生物群落结构的变化也会影响抗性基因的扩散。当沉积物中的微生物群落结构发生改变时，如优势菌种的更替、微生物多样性的变化等，会影响抗性基因在微生物之间的传播和扩散。在富营养化的沉积物中，一些适应富营养环境的微生物大量繁殖，这些微生物可能携带更多的抗性基因，并且它们之间的相互作用可能促进抗性基因的传播，导致抗性基因在沉积物中的扩散范围扩大。抗性基因在沉积物中的传播与扩散对周边环境具有潜在的影响。一方面，抗性基因可能通过沉积物与水体之间的物质交换，重新进入水体，导致水体中抗性基因的浓度增加，进一步威胁水生生态系统的健康。当沉积物中的抗性基因随着水流的扰动或生物扰动重新释放到水体中时，会使水体中的抗性基因污染范围扩大，影响水体中的微生物群落结构和生态功能。另一方面，抗性基因还可能通过食物链的传递，从沉积物中的微生物转移到水生生物体内，进而对人类健康构成潜在威胁。如果水生生物摄取了含有抗性基因的沉积物颗粒或微生物，抗性基因可能在其体内积累，当人类食用这些水生生物时，就有可能摄入抗性基因，增加人类感染耐药菌的风险。五、抗生素与抗生素抗性基因的相关性及影响因素5.1抗生素与抗生素抗性基因的相关性分析为了深入探究抗生素残留与抗生素抗性基因之间的内在联系，本研究运用Pearson相关性分析方法，对珠江口及其沿岸典型水产养殖区水体和沉积物中抗生素浓度与抗性基因丰度数据进行了详细分析。在水体中，结果显示部分抗生素与相应抗性基因之间存在显著的正相关关系。喹诺酮类抗生素中的诺氟沙星浓度与喹诺酮类抗性基因qnrA丰度的相关系数达到了0.78（P<0.01），呈现出极显著的正相关。这表明随着水体中诺氟沙星浓度的升高，qnrA基因的丰度也会显著增加，诺氟沙星的存在对qnrA基因的产生和传播具有明显的促进作用。这种相关性的产生可能是由于诺氟沙星作为一种广泛使用的喹诺酮类抗生素，长期存在于水体中，对微生物产生了强烈的选择压力。在这种选择压力下，原本对诺氟沙星敏感的微生物通过基因突变、基因水平转移等方式获得qnrA抗性基因，从而能够在含有诺氟沙星的环境中生存和繁殖，导致qnrA基因的丰度增加。四环素类抗生素中的四环素浓度与四环素类抗性基因tetA丰度之间也存在显著的正相关，相关系数为0.65（P<0.05）。四环素在水产养殖中的广泛使用，使得水体中存在一定浓度的四环素，这促使微生物产生对四环素的耐药性，tetA基因作为一种常见的四环素抗性基因，其丰度也随之增加。tetA基因主要通过主动外排机制，将四环素类抗生素排出细胞外，使细菌获得耐药性。当水体中四环素浓度升高时，携带tetA基因的细菌在竞争中更具优势，其数量和丰度也会相应增加。然而，并非所有抗生素与抗性基因之间都呈现出明显的相关性。磺胺类抗生素中的磺胺嘧啶浓度与磺胺类抗性基因sul1丰度之间的相关性较弱，相关系数仅为0.25（P>0.05），无统计学意义。这可能是因为磺胺嘧啶在水体中的降解速度较快，或者其在环境中的迁移转化过程较为复杂，导致其与sul1基因之间的关联不明显。磺胺类抗生素的作用机制与其他类抗生素有所不同，其耐药机制也更为复杂，除了sul1基因介导的耐药机制外，还可能存在其他耐药途径，这也可能导致磺胺嘧啶与sul1基因之间的相关性不显著。在沉积物中，同样观察到部分抗生素与抗性基因之间的显著相关性。四环素类抗生素中的土霉素浓度与四环素类抗性基因tetM丰度的相关系数为0.72（P<0.01），呈现出极显著的正相关。土霉素在沉积物中的残留可能持续对微生物产生选择压力，促使携带tetM基因的微生物在沉积物中大量繁殖，从而导致tetM基因丰度增加。tetM基因通过核糖体保护机制，阻止四环素类抗生素与细菌核糖体的结合，使细菌产生耐药性。当沉积物中土霉素浓度较高时，携带tetM基因的细菌能够更好地生存和繁殖，tetM基因的丰度也会相应提高。喹诺酮类抗生素中的氧氟沙星浓度与喹诺酮类抗性基因qnrB丰度之间也存在显著的正相关，相关系数为0.68（P<0.05）。氧氟沙星在沉积物中的积累会诱导微生物产生对其的耐药性，qnrB基因作为一种重要的喹诺酮类抗性基因，其丰度也会随之增加。qnrB基因可以通过改变细菌细胞膜的通透性，降低喹诺酮类抗生素进入细胞的量，从而使细菌产生耐药性。当沉积物中氧氟沙星浓度升高时，携带qnrB基因的细菌在生存竞争中更具优势，其丰度也会相应上升。与水体情况类似，沉积物中也有部分抗生素与抗性基因之间相关性不明显。磺胺类抗生素中的磺胺甲恶唑浓度与磺胺类抗性基因sul2丰度之间的相关系数为0.30（P>0.05），无统计学意义。这可能是由于磺胺甲恶唑在沉积物中的吸附、解吸过程以及与其他物质的相互作用较为复杂，影响了其与sul2基因之间的关联。沉积物中的微生物群落结构、物理化学性质等因素也可能对磺胺甲恶唑与sul2基因之间的相关性产生影响，导致两者之间的关系不显著。5.2环境因子对抗生素及抗生素抗性基因的影响环境因子在抗生素降解和抗性基因传播过程中扮演着关键角色，它们通过多种复杂的机制影响着抗生素和抗生素抗性基因在珠江口及其沿岸典型水产养殖区的环境行为。温度是影响抗生素降解和抗性基因传播的重要环境因子之一。温度对微生物的生长、代谢和酶活性有着直接的影响，进而影响抗生素的降解和抗性基因的传播。在一定温度范围内，随着温度的升高，微生物的代谢活动增强，酶的活性提高，这有利于抗生素的生物降解。对于四环素类抗生素，在温度为30℃时，其在水体中的降解速率明显高于20℃时，这是因为较高的温度促进了微生物的生长和代谢，使得能够降解四环素类抗生素的微生物数量增加，酶活性增强，从而加快了抗生素的降解速度。然而，当温度过高时，可能会对微生物产生不利影响，导致微生物的生长受到抑制，酶活性降低，反而不利于抗生素的降解。当温度超过40℃时，部分微生物的蛋白质和核酸结构可能会受到破坏，影响其正常的生理功能，使得抗生素的降解速率下降。温度还会影响抗性基因的传播。较高的温度可以增加微生物的活性和繁殖速度，促进抗性基因在微生物之间的水平转移。在夏季高温时，水体中微生物的代谢活动旺盛，细菌之间的相互作用频繁，这使得抗性基因更容易通过转化、转导和接合等方式在微生物之间传播，导致抗性基因的丰度增加。温度的变化还可能影响微生物细胞膜的流动性和通透性，从而影响抗性基因的摄取和表达。在低温条件下，微生物细胞膜的流动性降低，可能会阻碍抗性基因的进入，减少抗性基因的传播。pH值同样对抗生素降解和抗性基因传播具有重要影响。pH值的变化会影响抗生素的化学形态和微生物的生理活性，进而影响抗生素的降解和抗性基因的传播。不同抗生素在不同pH值条件下的稳定性和降解途径存在差异。磺胺类抗生素在酸性条件下相对稳定，而在碱性条件下容易发生水解反应，从而促进其降解。当水体pH值为8.0时，磺胺类抗生素的降解速率明显加快，这是因为碱性条件促进了磺胺类抗生素的水解反应，使其分解为小分子物质，从而实现降解。而对于喹诺酮类抗生素，其在酸性和中性条件下相对稳定，在碱性条件下可能会发生结构变化，影响其抗菌活性和降解特性。pH值还会影响微生物的生长和代谢，从而影响抗性基因的传播。不同微生物对pH值的适应范围不同，当环境pH值适宜时，微生物的生长和繁殖速度加快，抗性基因的传播也会相应增加。在pH值为7.5的水体中，携带四环素类抗性基因的微生物生长良好，其数量增加，抗性基因的传播也更为活跃。pH值的变化还可能影响微生物表面的电荷性质，进而影响抗性基因的吸附和解吸，以及微生物之间的相互作用，最终影响抗性基因的传播。溶解氧是水产养殖环境中的关键参数之一，对抗生素降解和抗性基因传播有着显著影响。溶解氧含量会影响微生物的代谢方式和生长环境，进而影响抗生素的降解和抗性基因的传播。在好氧条件下，微生物可以利用氧气进行有氧呼吸，代谢活动较为活跃，这有利于抗生素的生物降解。一些好氧微生物能够分泌特定的酶，将抗生素分解为无害物质。在溶解氧充足的水体中，四环素类抗生素可以被好氧微生物通过氧化、水解等方式降解。而在厌氧条件下，微生物则进行无氧呼吸，其代谢途径和产物与好氧条件下不同，可能会导致抗生素的降解速率降低，或者产生一些难以降解的中间产物。溶解氧含量还会影响抗性基因的传播。在好氧条件下，微生物的活性较高，抗性基因的表达和传播可能会增强。一些携带抗性基因的微生物在好氧环境中更容易生长和繁殖，从而增加抗性基因的传播机会。而在厌氧条件下，微生物的群落结构和代谢方式发生改变，可能会抑制抗性基因的传播。在厌氧的沉积物中，由于微生物的代谢活动受到限制，抗性基因在微生物之间的转移频率相对较低，从而减少了抗性基因的传播范围。营养盐（如氮、磷等）也是影响抗生素降解和抗性基因传播的重要环境因子。营养盐是微生物生长和代谢所必需的物质，其含量的变化会影响微生物的生长和群落结构，进而影响抗生素的降解和抗性基因的传播。当水体中营养盐充足时，微生物的生长速度加快，数量增加，这有利于抗生素的生物降解。充足的氮、磷营养盐可以为微生物提供合成细胞物质和酶的原料，促进微生物的代谢活动，使其能够更好地降解抗生素。在富营养化的水体中，由于营养盐丰富，微生物数量增多，对四环素类抗生素的降解能力增强。然而，当营养盐过量时，可能会导致水体富营养化，引发藻类大量繁殖，从而改变水体的生态环境，对抗生素的降解和抗性基因的传播产生间接影响。藻类的大量繁殖可能会消耗水中的溶解氧，导致水体缺氧，影响微生物的生长和代谢，进而影响抗生素的降解和抗性基因的传播。营养盐含量还会影响抗性基因的传播。营养盐的变化会改变微生物的群落结构，不同微生物对抗生素的敏感性和抗性基因的携带情况不同，因此营养盐的变化可能会导致抗性基因在微生物群落中的分布和传播发生改变。在营养盐丰富的水体中，一些携带抗性基因的微生物可能会成为优势种群，从而增加抗性基因的传播范围和丰度。营养盐还可能通过影响微生物之间的相互作用，如竞争、共生等关系，来影响抗性基因的传播。5.3养殖模式对抗生素及抗生素抗性基因的影响不同养殖模式下，抗生素的使用量、残留水平以及抗性基因的污染特征存在显著差异。在池塘养殖模式中，由于水体相对封闭，自净能力较弱，养殖户为了预防和控制疾病，往往会频繁且大量地使用抗生素。据调查，池塘养殖中抗生素的平均使用量达到了每立方米水体5-10克，显著高于其他养殖模式。在某高密度池塘养殖区，为了控制草鱼的细菌性疾病，养殖户每月使用抗生素的频率高达3-4次，每次使用量为每立方米水体8克左右。这种大量使用抗生素的情况导致池塘水体和沉积物中抗生素残留水平较高。在该池塘养殖区，水体中喹诺酮类抗生素的浓度可达到50-100ng/L，沉积物中喹诺酮类抗生素的浓度则可达到40-80ng/g（干重）。长期的高浓度抗生素残留使得池塘养殖环境中微生物受到强烈的选择压力，从而导致抗生素抗性基因的丰度较高。在该池塘养殖区，磺胺类抗性基因sul1在水体中的丰度可达5.0×10^7copies/L，在沉积物中的丰度可达8.0×10^8copies/g（干重）。网箱养殖模式设置在水流相对平缓、水质较好的海域或大型水域，水体交换相对较为频繁，一定程度上有利于抗生素的稀释和扩散。然而，由于网箱养殖密度通常较大，养殖生物的排泄物和残饵容易在网箱周围积累，增加了疾病发生的风险，促使养殖户使用抗生素。网箱养殖中抗生素的平均使用量为每立方米水体3-5克，相对池塘养殖有所降低。在某网箱养殖区，为了预防鲈鱼的疾病，养殖户每2-3周使用一次抗生素，每次使用量为每立方米水体4克左右。由于水体的流动性，网箱养殖水体中抗生素残留水平相对较低，一般在10-30ng/L之间，但沉积物中抗生素残留水平仍较高，可达30-60ng/g（干重）。这是因为虽然水体中的抗生素能够被稀释，但随着水流的作用，抗生素会逐渐沉积到网箱底部的沉积物中。在抗生素的选择压力下，网箱养殖环境中也检测到了较高丰度的抗生素抗性基因，如四环素类抗性基因tetA在水体中的丰度可达3.0×10^7copies/L，在沉积物中的丰度可达6.0×10^8copies/g（干重）。海水养殖模式利用海水资源，水体交换更为频繁，自净能力相对较强。在海水养殖中，抗生素的使用量相对较少，平均使用量为每立方米水体1-3克。在某海水养殖区，养殖户主要依靠良好的水质管理和生态防控措施来预防疾病，每月使用抗生素的次数不超过1次，每次使用量为每立方米水体2克左右。因此，海水养殖水体中抗生素残留水平较低，一般在5-15ng/L之间，沉积物中抗生素残留水平也相对较低，在10-30ng/g（干重）之间。由于抗生素使用量较少，海水养殖环境中抗生素抗性基因的丰度相对较低，如喹诺酮类抗性基因qnrA在水体中的丰度一般在1.0×10^6copies/L以下，在沉积物中的丰度在2.0×10^7copies/g（干重）以下。不同养殖模式下抗生素抗性基因的种类和分布也存在差异。池塘养殖中，由于长期大量使用多种抗生素，导致多种抗生素抗性基因共存，且分布较为广泛。除了常见的磺胺类、四环素类和喹诺酮类抗性基因外，还检测到一些少见的抗性基因，如氯霉素类抗性基因cat等。在某池塘养殖区，同时检测到了sul1、sul2、tetA、tetB、qnrA、qnrB以及cat等多种抗性基因，这些抗性基因在水体和沉积物中均有较高的丰度。网箱养殖中，抗性基因的种类相对较少，但某些抗性基因的丰度较高。由于网箱养殖中喹诺酮类抗生素的使用相对较多，喹诺酮类抗性基因qnrA、qnrB等在网箱养殖环境中较为常见，且丰度较高。在某网箱养殖区，qnrA基因在水体中的丰度可达3.0×10^7copies/L，在沉积物中的丰度可达6.0×10^8copies/g（干重），而其他类抗性基因的丰度相对较低。海水养殖中，由于抗生素使用量少，抗性基因的种类和丰度都相对较低，主要检测到一些常见的抗性基因，如sul1、tetA等，且丰度明显低于池塘养殖和网箱养殖。在某海水养殖区，sul1基因在水体中的丰度为5.0×10^5copies/L，在沉积物中的丰度为1.0×10^6copies/g（干重），tetA基因在水体中的丰度为3.0×10^5copies/L，在沉积物中的丰度为8.0×10^5copies/g（干重）。六、抗生素及抗生素抗性基因的环境风险评估6.1风险评估方法本研究采用风险商值法（RiskQuotient，RQ）对珠江口及其沿岸典型水产养殖区抗生素及抗生素抗性基因的环境风险进行评估。风险商值法是一种广泛应用于环境风险评估的方法，通过计算预测环境浓度（PredictedEnvironmentalConcentration，PEC）与预测无效应浓度（PredictedNo-EffectConcentration，PNEC）的比值，来判断污染物对环境的风险程度。对于抗生素的风险评估，预测环境浓度（PEC）直接采用本研究中在水体和沉积物样品中检测到的抗生素实测浓度。在计算水体中喹诺酮类抗生素诺氟沙星的PEC时，直接使用在各采样点水体中检测到的诺氟沙星浓度数据。预测无效应浓度（PNEC）则通过查阅相关文献资料获取，不同类型抗生素的PNEC值因作用机制、毒性等因素而异。磺胺类抗生素的PNEC值通常在0.01-0.1μg/L之间，这是基于其对水生生物的毒性研究以及环境安全阈值的设定得出的。对于缺乏相关PNEC数据的抗生素，采用物种敏感度分布（SSD）法进行推导。通过收集该抗生素对多种水生生物的毒性数据，构建物种敏感度分布曲线，从而计算出其PNEC值。在评估抗生素抗性基因的环境风险时，由于目前对抗生素抗性基因的风险评估尚缺乏统一的标准和方法，本研究参考相关研究，以抗性基因的丰度作