德普昌达克冰川和中山站抗生素抗性基因及影响因素研究

张淑红, 何瑞, 侯书贵, 刘科, 李治国. 德普昌达克冰川和中山站抗生素抗性基因及影响因素研究[J]. 环境保护科学, 2021, 47(3): 151-158. doi: 10.16803/j.cnki.issn.1004-6216.2021.03.026
引用本文: 张淑红, 何瑞, 侯书贵, 刘科, 李治国. 德普昌达克冰川和中山站抗生素抗性基因及影响因素研究[J]. 环境保护科学, 2021, 47(3): 151-158. doi: 10.16803/j.cnki.issn.1004-6216.2021.03.026
ZHANG Shuhong, HE Rui, HOU Shugui, LIU Ke, LI Zhiguo. Antibiotic resistance genes and their influencing factors in Depuchangdake Glacier and Zhongshan Station[J]. Environmental Protection Science, 2021, 47(3): 151-158. doi: 10.16803/j.cnki.issn.1004-6216.2021.03.026
Citation: ZHANG Shuhong, HE Rui, HOU Shugui, LIU Ke, LI Zhiguo. Antibiotic resistance genes and their influencing factors in Depuchangdake Glacier and Zhongshan Station[J]. Environmental Protection Science, 2021, 47(3): 151-158. doi: 10.16803/j.cnki.issn.1004-6216.2021.03.026

德普昌达克冰川和中山站抗生素抗性基因及影响因素研究

    作者简介: 张淑红(1977-),女,博士、教授。研究方向:环境微生物。E-mail:shuhongzhang_2013@163.com
  • 基金项目:
    国家自然科学基金重点项目(41830644);国家自然科学基金面上基金项目(41671072);河南省自然科学基金资助项目(182300410137)
  • 中图分类号: X171

Antibiotic resistance genes and their influencing factors in Depuchangdake Glacier and Zhongshan Station

  • 摘要: 为了解冰川环境抗生素抗性基因(Antibiotic Resistance Genes,ARGs)的存在状况及其与重金属和I型整合子(class I integron,intI)的关系,通过荧光定量PCR(Real-time Quantitative PCR,qPCR),研究青藏高原德普昌达克冰川和南极中山站雪样ARGs、重金属抗性基因(heavy metaI resistance genes,HMRGs)和intI的相对丰度及其相关性。结果表明,大部分基因的相对丰度在德普昌达克冰川小于中山站。德普昌达克冰川tetA、czcD、pcoA、znuA、intI1和intI2的相对丰度随海拔高度的增加而降低;南极的merA、zntA、znuA、sul2、tetA和intI2距离中山站建筑区域越近,相对丰度越高。2个研究位点ARGs均受Zn抗性基因znuA的影响,且znuA与intI之间显著相关。冰川环境虽然远离人类活动,但其ARGs同时受intI和与人类活动有关的重金属Zn的影响。
  • 加载中
  • 图 1  德普昌达克冰川和中山站雪样抗性基因及intI丰度分布

    表 1  qPCR引物名称、序列、片段扩增长度和退火温度

    目标基因目的基因引物序列(5′→3′)片段长度/bp退火温度/℃文献
    HMRGsmerAF: GTGCCGTCCAAGATCATG
    R: GGTGGAAGTCCAGTAGGGTGA
    32160[21]
    zntAF: GGTCGGGTCTGGCATTGAAG
    R: TTGCAGCATCGGCGCGCAGGGTA
    26360[22]
    czcDF: TCATCG CCG GTGCGATCATCAT
    R:TGT CAT TCA CGA CAT GAA CC
    27255[19]
    pcoAF: GCTGCAGATGGCCAGTATGTAAA
    R: CCCTCGAGCGTAACCGGTCC
    14760[23]
    znuAF: TCATCAGTAGCGGTTTCACA
    R: CACAGGTTGCTGAACTCGCC
    21358[24]
    terDF:AGTAAAGCAGCTCCGTCAAT
    R: CCGAACAGCATGGCAGTCT
    43455[25]
    chrAF: TCACGCCGGAATATAACTAC
    R: CGTACCCTGATCAATCACTT
    22955[26]
    czcAF:TCGACGGBGCCGTGGTSMTBGTCGAGAA
    R: GTVAWSGCCAKCGGVBGGAACA
    23263[22]
    mntHF: TCACCGTGGCATACAGTGACACAC
    R: TGAAATTGTTTTAGCGCACGACCT
    18763[27]
    ARGsaac(6')-Ib-crF: TTGCGATGCTCTATGAGTGGCTA
    R: CTCGAATGCCTGGCGTGTTT
    48255[28]
    qnrAF: AGAGGATTTCTCACGCCAGG
    R: GCAGCACTATKACTCCCAAGG
    61958[29]
    tetYF:GCTGATATTTGCGGGTTTCTA
    R:CGTCAAGCCTGTTAAAGTTCC
    17758[30]
    tetMFF:ACAGAAAGCTTATTATATAAC
    R:TGGCGTGTCTATGATGTTCAC
    17155[31]
    tetAF:GCTACATCCTGCTTGCCTTC
    R: CATAGATCGCCGTGAAGAGG
    21064[32]
    sul1F:CGCACCGGAAACATCGCTGCAC
    R: TGAAGTTCCGCCGCAAGGCTCG
    16362[33]
    sul2F:TCCGGTGGAGGCCGGTATCTGG
    R:CGGGAATGCCATCTGCCTTGAG
    19162[33]
    ermBF: TAAAGGGCATTTAACGACGAAACT
    R:TTTATACCTCTGTTTGTTAGGGAATTGAA
    17258[34]
    intIintI1F: GGCTTCGTGATGCCTGCTT
    R: CATTCCTGGCCGTGGTTCT
    14859[35]
    intI2F: GTTATTTTATTGCTGGGATTAGGC
    R: TTTTACGCTGCTGTATGGTGC
    16655[36]
    intI3F: GGATGTCTGTGCCTGCTTG
    R: GCCACCACTTGTTTGAGGA
    10060[37]
    16S rRNA338F: ACTCCTACGGGAGGCAGCAG
    518R: ATTACCGCGGCTGCTGG
    20055[38]
    目标基因目的基因引物序列(5′→3′)片段长度/bp退火温度/℃文献
    HMRGsmerAF: GTGCCGTCCAAGATCATG
    R: GGTGGAAGTCCAGTAGGGTGA
    32160[21]
    zntAF: GGTCGGGTCTGGCATTGAAG
    R: TTGCAGCATCGGCGCGCAGGGTA
    26360[22]
    czcDF: TCATCG CCG GTGCGATCATCAT
    R:TGT CAT TCA CGA CAT GAA CC
    27255[19]
    pcoAF: GCTGCAGATGGCCAGTATGTAAA
    R: CCCTCGAGCGTAACCGGTCC
    14760[23]
    znuAF: TCATCAGTAGCGGTTTCACA
    R: CACAGGTTGCTGAACTCGCC
    21358[24]
    terDF:AGTAAAGCAGCTCCGTCAAT
    R: CCGAACAGCATGGCAGTCT
    43455[25]
    chrAF: TCACGCCGGAATATAACTAC
    R: CGTACCCTGATCAATCACTT
    22955[26]
    czcAF:TCGACGGBGCCGTGGTSMTBGTCGAGAA
    R: GTVAWSGCCAKCGGVBGGAACA
    23263[22]
    mntHF: TCACCGTGGCATACAGTGACACAC
    R: TGAAATTGTTTTAGCGCACGACCT
    18763[27]
    ARGsaac(6')-Ib-crF: TTGCGATGCTCTATGAGTGGCTA
    R: CTCGAATGCCTGGCGTGTTT
    48255[28]
    qnrAF: AGAGGATTTCTCACGCCAGG
    R: GCAGCACTATKACTCCCAAGG
    61958[29]
    tetYF:GCTGATATTTGCGGGTTTCTA
    R:CGTCAAGCCTGTTAAAGTTCC
    17758[30]
    tetMFF:ACAGAAAGCTTATTATATAAC
    R:TGGCGTGTCTATGATGTTCAC
    17155[31]
    tetAF:GCTACATCCTGCTTGCCTTC
    R: CATAGATCGCCGTGAAGAGG
    21064[32]
    sul1F:CGCACCGGAAACATCGCTGCAC
    R: TGAAGTTCCGCCGCAAGGCTCG
    16362[33]
    sul2F:TCCGGTGGAGGCCGGTATCTGG
    R:CGGGAATGCCATCTGCCTTGAG
    19162[33]
    ermBF: TAAAGGGCATTTAACGACGAAACT
    R:TTTATACCTCTGTTTGTTAGGGAATTGAA
    17258[34]
    intIintI1F: GGCTTCGTGATGCCTGCTT
    R: CATTCCTGGCCGTGGTTCT
    14859[35]
    intI2F: GTTATTTTATTGCTGGGATTAGGC
    R: TTTTACGCTGCTGTATGGTGC
    16655[36]
    intI3F: GGATGTCTGTGCCTGCTTG
    R: GCCACCACTTGTTTGAGGA
    10060[37]
    16S rRNA338F: ACTCCTACGGGAGGCAGCAG
    518R: ATTACCGCGGCTGCTGG
    20055[38]
    下载: 导出CSV

    表 2  ARGs、HMRGs和intI之间斯皮尔曼相关性分析(n=3)

    基因名称德普昌达克冰川中山站
    sul2tetAintI1intI2intI3sul2tetAintI1intI2intI3
    merA−0.50−1.00**−0.50−1.00**−0.501.00**1.00**−1.00**1.00**0.50
    zntA−0.500.50−0.500.50−0.501.00**1.00**−1.00**1.00**0.50
    czcD0.501.00**0.501.00**0.50−0.50−0.500.50−0.50−1.00**
    pcoA0.501.00**0.501.00**0.500.500.50−0.500.501.00**
    znuA0.501.00**0.501.00**0.501.00**1.00**−1.00**1.00**0.50
    terD1.00**0.501.00**0.501.00**0.500.50−0.500.501.00**
    chrA−0.50−1.00**−0.50−1.00**−0.500.500.50−0.500.501.00**
    czcA1.00**0.501.00**0.501.00**0.500.50−0.500.501.00**
    mntH−0.50−1.00**−0.50−1.00**−0.500.500.50−0.500.501.00**
    sul21.000.501.00**0.501.00**1.001.00**−1.00**1.00**0.50
    tetA0.501.000.501.00**0.501.00**1.00−1.00**1.00**0.50
       注:**表示极显著相关(P<0.01,双边检验)。
    基因名称德普昌达克冰川中山站
    sul2tetAintI1intI2intI3sul2tetAintI1intI2intI3
    merA−0.50−1.00**−0.50−1.00**−0.501.00**1.00**−1.00**1.00**0.50
    zntA−0.500.50−0.500.50−0.501.00**1.00**−1.00**1.00**0.50
    czcD0.501.00**0.501.00**0.50−0.50−0.500.50−0.50−1.00**
    pcoA0.501.00**0.501.00**0.500.500.50−0.500.501.00**
    znuA0.501.00**0.501.00**0.501.00**1.00**−1.00**1.00**0.50
    terD1.00**0.501.00**0.501.00**0.500.50−0.500.501.00**
    chrA−0.50−1.00**−0.50−1.00**−0.500.500.50−0.500.501.00**
    czcA1.00**0.501.00**0.501.00**0.500.50−0.500.501.00**
    mntH−0.50−1.00**−0.50−1.00**−0.500.500.50−0.500.501.00**
    sul21.000.501.00**0.501.00**1.001.00**−1.00**1.00**0.50
    tetA0.501.000.501.00**0.501.00**1.00−1.00**1.00**0.50
       注:**表示极显著相关(P<0.01,双边检验)。
    下载: 导出CSV
  • [1] HU H W, WANG J T, LI J, et al. Field-based evidence for copper contamination induced changes of antibiotic resistance in agricultural soils[J]. Environmental Microbiology, 2016, 18(11): 3896 − 3909. doi: 10.1111/1462-2920.13370
    [2] LAFFITE A, KILUNGA P I, KAYEMBE J M, et al. Hospital effluents are one of several sources of metal, antibiotic resistance genes, and bacterial markers disseminated in sub-Saharan urban rivers[J]. Frontiers in Microbiology, 2016, 7: 1128.
    [3] MOSKVITCH K. UK prize lets public decide on world’s biggest science problem[J]. Nature, 2014.
    [4] HSU J T, CHEN C Y, YOUNG C W, et al. Prevalence of sulfonamide resistant bacteria, resistance genes and integron-associated horizontal gene transfer in natural water bodies and soils adjacent to a swine feedlot in northern Taiwan[J]. Journal of Hazardous Materials, 2014, 227: 34 − 43.
    [5] BENGTSSON-PALME J, BOULUND F, FICK J, et al. Shotgun metagenomics reveals a wide array of antibiotic resistance genes and mobile elements in a polluted lake in India[J]. Frontiers in Microbiology, 2014, 5: e648.
    [6] USHIDA K, SEGAWA T, KOHSHIMA S, et al. Application of real-time PCR array to the multiple detection of antibiotic resistant genes in glacier ice samples[J]. Journal of General and Applied Microbiology, 2010, 56(1): 43 − 52. doi: 10.2323/jgam.56.43
    [7] SEGAWA T, TAKEUCHI N, RIVERA A, et al. Distribution of antibiotic resistance genes in glacier environments[J]. Environmental Microbiology Reports, 2013, 5(1): 127 − 134. doi: 10.1111/1758-2229.12011
    [8] VAN GOETHEM MW, PIERNEEF R, BEZUIDT O K I, et al. A reservoir of ‘historical’ antibiotic resistance genes in remote pristine Antarctic soils[J]. Microbiome, 2018, 6: 40. doi: 10.1186/s40168-018-0424-5
    [9] WANG F, STEDTFELD R D, KIM O S, et al. Influence of soil characteristics and proximity to Antarctic research stations on abundance of antibiotic resistance genes in soils[J]. Environmental Science & Technology, 2016, 50(23): 12621 − 12629.
    [10] OKUBO T, AE R, NODA J, et al. Detection of the sul2-strA-strB gene cluster in an ice core from Dome Fuji Station, East Antarctica[J]. Journal of Global Antimicrobial Resistance, 2019, 17: 72 − 78. doi: 10.1016/j.jgar.2018.11.005
    [11] CZEKALSKI N, GASCÓN DÍEZ E, BÜRGMANN H. Wastewater as a point source of antibiotic-resistance genes in the sediment of a freshwater lake[J]. The ISME Journal, 2014, 8: 1381 − 1390. doi: 10.1038/ismej.2014.8
    [12] ZHOU T, LU J J, TONG Y B, et al. Distribution of antibiotic resistance genes in Bosten Lake, Xinjiang, China[J]. Water Science & Technology, 2014, 70(5): 925 − 931.
    [13] BERG J, THORSEN M K, HOLM P E, et al. Cu exposure under field conditions coselects for antibiotic resistance as determined by a novel cultivation-independent bacterial community tolerance assay[J]. Environmental Science & Technology, 2010, 44(22): 8724 − 8728.
    [14] CESARE A D, ECKERT E M, DURSO S, et al. Co-occurrence of integrase 1, antibiotic and heavy metal resistance genes in municipal wastewater treatment plants[J]. Water Research, 2016, 94: 208 − 214. doi: 10.1016/j.watres.2016.02.049
    [15] JANG H M, LEE J, CHOI S, et al. Response of antibiotic and heavy metal resistance genes to two different temperature sequences in anaerobic digestion of waste activated sludge[J]. Bioresource Technology, 2018, 267: 303 − 310. doi: 10.1016/j.biortech.2018.07.051
    [16] CUI E P, WU Y, JIAO Y N, et al. The behavior of antibiotic resistance genes and arsenic influenced by biochar during different manure composting[J]. Environmental Science and Pollution Research, 2017, 24(16): 14484 − 14490. doi: 10.1007/s11356-017-9028-z
    [17] SANDBERG K D, ISHII S, LAPARA T M. A microfluidic quantitative polymerase chain reaction method for the simultaneous analysis of dozens of antibiotic resistance and heavy metal resistance genes[J]. Environmental Science & Technology Letters, 2018, 5(1): 20 − 25.
    [18] YANG Y Y, XU C, CAO X H, et al. Antibiotic resistance genes in surface water of eutrophic urban lakes are related to heavy metals, antibiotics, lake morphology and anthropic impact[J]. Ecotoxicology, 2017, 26(6): 831 − 840. doi: 10.1007/s10646-017-1814-3
    [19] 刘亚平. 喜马拉雅山中段和东天山雪冰痕量元素环境记录研究[D]. 北京: 中国科学院大学, 2010.
    [20] ROOSA S, WATTIEZ R, PRYGIEL E, et al. Bacterial metal resistance genes and metal bioavailability in contaminated sediments[J]. Environmental Pollution, 2014, 189: 143 − 151. doi: 10.1016/j.envpol.2014.02.031
    [21] PITKÄNEN L K, TAMMINEN M, HYNNINEN A, et al. Fish farming affects the abundance and diversity of the mercury resistance merA in marine sediments[J]. Microbes and Environments, 2001, 26(3): 205 − 211.
    [22] ZHOU Q, WANG M Z, ZHONG X X, et al. Dissemination of resistance genes in duck/fish polyculture ponds in Guangdong Province: correlations between Cu and Zn and antibiotic resistance genes[J]. Environmental Science and Pollution Research, 2019, 26(8): 8182 − 8193. doi: 10.1007/s11356-018-04065-2
    [23] XIONG W G, ZENG Z L, ZHANG Y M, et al. Fate of metal resistance genes in arable soil after manure application in a microcosm study[J]. Ecotoxicology and Environmental Safety, 2015, 113: 59 − 63. doi: 10.1016/j.ecoenv.2014.11.026
    [24] HENSLEY M P, GUNASEKERA T S, EASTON J A, et al. Characterization of Zn(II)-responsive ribosomal proteins YkgM and L31 in E. coli[J]. Journal of Inorganic Biochemistry, 2012, 111: 164 − 172. doi: 10.1016/j.jinorgbio.2011.11.022
    [25] BIELASZEWSKA M, MELLMANN A, ZHANG W, et al. Characterisation of the Escherichia coli strain associated with an outbreak of haemolyticuraemic syndrome in Germany, 2011: a microbiological study[J]. The Lancet Infectious Diseases, 2011, 11(9): 671 − 676. doi: 10.1016/S1473-3099(11)70165-7
    [26] PATRA RC, MALIK S, BEER M. Molecular characterization of chromium (VI) reducing potential in Gram positive bacteria isolated from contaminated sites[J]. Soil Biology & Biochemistry, 2010, 42(10): 1857 − 1863.
    [27] TONG Y, ZHAI Q, WANG G, et al. System-wide analysis of manganese starvation induced metabolism in key elements of Lactobacillus plantarum[J]. RSC Advances, 2017, 7: 12959 − 12968. doi: 10.1039/C7RA00072C
    [28] XIA L N, LI L, WU C M, et al. A survey of plasmid-mediated fluoroquinolone resistance genes from Escherichia coli isolates and their dissemination in Shandong, China[J]. Foodborne Pathogens and Disease, 2010, 7(2): 207 − 215. doi: 10.1089/fpd.2009.0378
    [29] CHEN X, ZHANG W Q, PAN W J, et al. Prevalence of qnr, aac (6′)-Ib-cr, qepA, and oqxAB in Escherichia coli isolates from humans, animals, and the environment[J]. Antimicrobial Agents and Chemotherapy, 2012, 56(6): 3423 − 3427. doi: 10.1128/AAC.06191-11
    [30] He X L, XU Y B, CHEN J L, et al. Evolution of corresponding resistance genes in the water of fish tanks with multiple stresses of antibiotics and heavy metals[J]. Water Research, 2017, 124: 39 − 48. doi: 10.1016/j.watres.2017.07.048
    [31] 刘锐. 规模化猪场废水典型抗生素抗性基因的调查研究[J]. 家畜生态学报, 2017, 38(11): 68 − 71. doi: 10.3969/j.issn.1673-1182.2017.11.013
    [32] NG LK, MARTIN I, ALFA M, et al. Multiplex PCR for the detection of tetracycline resistant genes[J]. Molecular and Cellular Probes, 2001, 15(4): 209 − 215. doi: 10.1006/mcpr.2001.0363
    [33] SU H C, YING G G, TAO R, et al. Occurrence of antibiotic resistance and characterization of resistance genes and integrons in Enterobacteriaceae isolated from integrated fish farms in South China[J]. Journal of Environmental Monitoring, 2011, 13(11): 3229 − 3236. doi: 10.1039/c1em10634a
    [34] 任省涛, 郭夏丽, 芦阿虔, 等. 林可霉素菌渣堆肥抗生素抗性基因变化分析[J]. 中国环境科学, 2018, 38(11): 4276 − 4283. doi: 10.3969/j.issn.1000-6923.2018.11.037
    [35] LUO Y, MAO D, RYSZ M, et al. Trends in antibiotic resistance genes occurrence in the Haihe River, China[J]. Environmental Science & Technology, 2010, 44(19): 7220 − 7225.
    [36] HE L Y, LIU Y S, SU H C, et al. Dissemination of antibiotic resistance genes in representative broiler feedlots environments: identification of indicator ARGs and correlations with environmental variables[J]. Environmental Science & Technology, 2014, 48(22): 13120 − 13129.
    [37] BARRAUD O, BACLET M C, DENIS F. Quantitative multiplex real-time PCR for detecting class 1, 2 and 3 integrons[J]. Journal of Antimicrobial Chemotherapy, 2010, 65(8): 1642 − 1645. doi: 10.1093/jac/dkq167
    [38] SONG Y, KONG Y L, WANG J C, et al. Identification of the produced volatile organic compounds and the involved soil bacteria during decomposition of watermelon plant residues in a Fusarium-infested soil[J]. Geoderma, 2018, 315: 178 − 187. doi: 10.1016/j.geoderma.2017.11.021
    [39] YANG Y Y, LIU G H, YE C, et al. Bacterial community and climate change implication affected the diversity and abundance of antibiotic resistance genes in wetlands on the Qinghai-Tibetan Plateau[J]. Journal of Hazardous Materials, 2019, 5: 283 − 293.
    [40] 方精云. 北极冻土的化学元素背景及其分布特征[J]. 环境科学学报, 2000, 20(1): 69 − 75. doi: 10.3321/j.issn:0253-2468.2000.01.014
    [41] SEILER C, BERENDONK T U. Heavy metal driven co-selection of antibiotic resistance in soil and water bodies impacted by agriculture and aquaculture[J]. Frontiers in Microbiology, 2012, 3: 399.
    [42] JI X, SHEN Q, LIU F, et al. Antibiotic resistance gene abundances associated with antibiotics and heavy metals in animal manures and agricultural soils adjacent to feedlots in Shanghai, China[J]. Journal of Hazardous Materials, 2012, 235–236: 178 − 185.
    [43] HUANG J, KANG S C, ZHANG Q G, et al. Spatial distribution and magnification processes of mercury in snow from high-elevation glaciers in the Tibetan Plateau[J]. Atmospheric Environment, 2012, 46: 140 − 146. doi: 10.1016/j.atmosenv.2011.10.008
    [44] VANDE VELDE K P, FERRARI C P, BARBANTE C, et al. A 200 year record of atmospheric cobalt, chromium, molybdenum and antimony in high altitude alpine firn and ice[J]. Environmental Science & Technology, 1999, 33: 3495 − 3501.
    [45] 汪大立, 康建成, 孙波, 等. 南极中山站至内陆冰盖330 km剖面表面雪样微粒分布特征及其意义[J]. 冰川冻土, 2000, 22(2): 128 − 134. doi: 10.3969/j.issn.1000-0240.2000.02.005
    [46] ZHANG Y, KANG S, CHEN P, et al. Records of anthropogenic antimony in the glacial snow from the southeastern Tibetan Plateau[J]. Journal of Asian Earth Science, 2016, 131: 62 − 71. doi: 10.1016/j.jseaes.2016.09.007
    [47] CHRISTNER BC. Recovery and identification of viable bacteria immured in glacial ice[J]. Icarus, 2000, 144: 479 − 485. doi: 10.1006/icar.1999.6288
    [48] SUN M, YE M, SCHWAB A P, et al. Human migration activities drive the fluctuation of ARGs: Case study of landfills in Nanjing, eastern China[J]. Journal of Hazardous Materials, 2016, 315: 93 − 101. doi: 10.1016/j.jhazmat.2016.04.077
    [49] CHEN B, LLIANG X, HUANG X, et al. Differentiating anthropogenic impacts on ARGs in the Pearl River Estuary by using suitable gene indicators[J]. Water Research, 2013, 47(8): 2811 − 2820. doi: 10.1016/j.watres.2013.02.042
    [50] BACCOLO G, DI MAURO B, MASSABÒD, et al. Cryoconite as a temporary sink for anthropogenic species stored in glaciers[J]. Scientific Reports, 2017, 7: 9623. doi: 10.1038/s41598-017-10220-5
  • 加载中
图( 1) 表( 2)
计量
  • 文章访问数:  1465
  • HTML全文浏览数:  1465
  • PDF下载数:  17
  • 施引文献:  0
出版历程
  • 收稿日期:  2020-07-20
  • 刊出日期:  2021-06-20

德普昌达克冰川和中山站抗生素抗性基因及影响因素研究

    作者简介: 张淑红(1977-),女,博士、教授。研究方向:环境微生物。E-mail:shuhongzhang_2013@163.com
  • 1. 商丘师范学院生物与食品学院,河南 商丘 476000
  • 2. 南京大学地理与海洋学院,江苏 南京 210093
  • 3. 商丘师范学院测绘与规划学院,河南 商丘 476000
基金项目:
国家自然科学基金重点项目(41830644);国家自然科学基金面上基金项目(41671072);河南省自然科学基金资助项目(182300410137)

摘要: 为了解冰川环境抗生素抗性基因(Antibiotic Resistance Genes,ARGs)的存在状况及其与重金属和I型整合子(class I integron,intI)的关系,通过荧光定量PCR(Real-time Quantitative PCR,qPCR),研究青藏高原德普昌达克冰川和南极中山站雪样ARGs、重金属抗性基因(heavy metaI resistance genes,HMRGs)和intI的相对丰度及其相关性。结果表明,大部分基因的相对丰度在德普昌达克冰川小于中山站。德普昌达克冰川tetA、czcD、pcoA、znuA、intI1和intI2的相对丰度随海拔高度的增加而降低;南极的merA、zntA、znuA、sul2、tetA和intI2距离中山站建筑区域越近,相对丰度越高。2个研究位点ARGs均受Zn抗性基因znuA的影响,且znuA与intI之间显著相关。冰川环境虽然远离人类活动,但其ARGs同时受intI和与人类活动有关的重金属Zn的影响。

English Abstract

  • 由于抗生素在农业[1]和医疗[2]环境的频繁使用,使得抗生素抗性成为21世纪6大难题之一[3]。抗生素的持续性选择压力,使得抗生素抗性基因(Antibiotic resistance genes,ARGs)不仅出现在人类活动频繁的环境[1-2, 4-5],在远离人类的地区也有报道[6-7],因此人类对ARGs的不断扩散非常担忧。冰川环境因抗生素选择压力比较小,能反应微生物对抗生素的固有抗性,因此是抗生素抗性风险评估的重要环境。该环境ARGs的研究有2种不同的结果:VAN GOETHEM et al[8]使用宏基因组测序技术,在南极Mackay冰川表层土壤中检测到的ARGs序列与其他环境具有不同的聚类分枝,这些ARGs是其进化的结果。然而,WANG et al[9]对南极冈瓦纳研究站的研究结果表明,人类的长期和短期活动都会影响ARGs在南极的传播,但长期影响和短期影响的结果差不多。OKUBO et al[10]在1200~1400年的南极冰芯检测到sul2,虽然赋予人工合成的磺胺甲恶唑抗性的基因在抗生素时代之前就已经存在,但sul2-strA-strB基因簇与当今细菌具有非常高的序列相似性,该结果与VAN GOETHEM et al[8]的报道相反。

    ARGs不仅会受到抗生素的选择压力,重金属同样也会对其产生影响[1112],如Cu会选择临床上重要的抗生素—万古霉素的耐药性[13]。重金属污染具有持久性,能在环境中存在数百年,其对ARGs的选择压力甚至比抗生素自身还要严重。重金属同样也会选择重金属抗性基因(Heavy metal resistance genes,HMRGs),因此ARGs和HMRGs可同时出现,使得微生物产生多重抗性。这类研究多集中在人类活动频繁的环境,如污水处理厂[14]、活性污泥[15]、农业土壤[1]、医院废水[2]、肥料[16]、城市污水和饮用水[17]、富营养化的地表水中[18],但尚未见到冰川环境HMRGs的研究。而且上述人类活动频繁的环境HMRGs的研究多集中在Cu、Zn、Hg、Ag、As、Ni和Cd等的抗性基因,关于Te和Mn等的抗性基因尚未见报道。

    冰川雪冰中重金属的浓度非常低,有时候甚至低于pg/g(10−12 g/g)水平[19],给其研究带来诸多困难。由于HMRGs的丰度与重金属浓度的密切相关性[20],以及荧光定量PCR(Real-time Quantitative PCR,qPCR)的敏感性,HMRGs可代替重金属,从侧面反应冰川中重金属的存在状况。作者之前通过宏基因组测序技术,在青藏高原西昆仑山的崇测冰帽冰川沉积物检测出Zn(9.13%)、Cu(22.64%)、Fe(14.76%)、As(7.35%)、Ni(0.61%)、Mo(9.86%)、Mn(9.16%)、Co(5.38%)、Hg(2.40%)、Cr(6.36%)、Cd(5.50%)、W(0.75%)和Te(6.09%)这13种HMRGs,其中Te和Mn所占百分比并不低,位于第4和第8位。本研究以距离崇测冰帽比较近的青藏高原喀喇昆仑山脉德普昌达克冰川,以及南极中山站为研究位点,通过qPCR,检测代表Hg、Cu、Zn、Cr、Te、Mn和Cd等金属的抗性基因,人工合成的氟喹诺酮类、磺胺类的抗性基因,天然的四环素和红霉素类抗性基因,以及在耐药性传递方面起重要作用的I型整合子(class I integron,intI),一方面初步判断ARGs与HMRGs在不同研究位点存在与否及相对丰度,另一方面了解ARGs、HMRGs与intI之间的相关性。该研究有助于了解冰川环境ARGs的存在状况,以及重金属和intI对其影响。

  • 2016年7月,采集德普昌达克冰川不同海拔高度的3个雪样:德普昌1(34°16'15.79″N,81°10'55.76″E,5433.78 ma.s.l)、德普昌2(34°16'16.42″N,81°10'55.76″E,5422.40 ma.s.l)和德普昌3(34°16'17.81″N,81°10'55.80″E,5411.75 ma.s.l),运输至商丘师范学院河南省特色微生物资源开发与应用工程研究中心,在洁净室进行雪样融化与过滤及DNA提取。

    2016年2月,中国第32次南极科考队在中山站采集3个雪样:中山1(69°22′22.93″S,76°22′05.94″E,11 m a.s.l,站区莫愁湖畔)、中山2(69°22′17.74″S,76°22′11.18″E,20 m a.s.l,正北方六角楼下)和中山3(69°22′24.03″S,76°21′30.72″E,25 m a.s.l,站区西边丘陵顶部)。所采样品经雪龙号运输至上海。由上海运输至南京大学海岸与海岛开发教育部重点实验室的低温室进行保存。2016年7月,在该重点实验室万级洁净室进行雪样融化与过滤及DNA提取。再将提取的DNA用冰袋和塑料泡沫盒运输至商丘师范学院,进行后续的qPCR分析。

    为保证所采雪样不受人为污染,采样者在采样过程中一律穿洁净服、戴口罩和PE手套。用无菌不锈钢铲去掉表层雪样1 cm,5点梅花取样法采取1 cm以下的表层样品。所采雪样装入无菌采样瓶,将每个采样瓶装入克林莱密实袋,以避免外部污染。所有样品置于冰柜中,以低温(≤–5 ℃)方式运回实验室。

  • 雪样置于4 ℃缓慢融化后,在万级洁净室内的超净工作台上,将约400~700 mL融水用小型切向流超滤系统(Labscale,Millipore)配套0.22 μm无菌Millipore微孔滤膜进行过滤。用灭菌剪刀将收集雪样的滤膜剪碎,置于2 mL离心管,按照DNA提取试剂盒(FastDNA Spin Kit for Soil,MP)说明书进行DNA提取。提取的DNA用微量紫外-可见光分光光度计(Thermo Scientific™ NanoDrop™ One)测OD值。所有样品DNA的A260/A280均介于1.8~2.0。

  • 所涉及目标基因包括8个ARGs、和3个intI,16S rDNA作为分析时的内参基因。引物名称、序列、片段长度和退火温度等,见表1

    qPCR反应体系为:SYBRR Premix Ex TaqTMⅡ(TliRNaseH Plus, Takara)10 µL,上下游引物(10 µmol/L)各0.8 µL,10 ng/µL DNA模板1.0 μL,ROX Reference Dye(50×)0.4 µL,无菌ddH2O7 µL。反应程序:95 ℃预变性1 min;95 ℃变性10 s,不同引物退火温度下40 s,72 ℃延伸30 s,共40个循环。用无菌ddH2O代替DNA模板作阴性对照,每个样品3个重复。

  • 试验数据以3次重复的平均值表示,目标基因的相对丰度为其浓度与16S rDNA浓度的比值。用SPSS 22.0软件的Spearman对ARGs、HMRGs和intI的相对丰度进行相关性分析。

  • 在8种ARGs中,aac(6')-Ib-cr仅在德普昌达克冰川被检出,sul1ermB仅在中山站被检出,而qnrAtetYtetM在2个位点均未被检出。仅sul2tetA在德普昌达克冰川和中山站2个位点都能被检出。这2种ARGs在2个位点的分布,见图1

    图1可见,sul2的相对丰度在德普昌达克冰川不随海拔增加而递增或递减,tetA的相对丰度在德普昌达克冰川随海拔升高而降低,sul2tetA的相对丰度在中山站随海拔升高而增加。

    aac(6')-Ib-cr的相对丰度在德普昌达克冰川的分布与sul2相同,不随海拔增加而递增或递减。sul1的相对丰度在中山站也未因海拔的增加而递增或递减。ermB的相对丰度在中山站的分布与sul2tetA相同,随海拔升高而增加。

  • 9个HMRGs在德普昌达克冰川和中山站均能被检测到,但在2个研究位点随海拔高度的变化,相对丰度的分布趋势有所差异。图1可见,znuAczcDpcoA的相对丰度在德普昌达克冰川随海拔升高而降低。然而,merAchrAmntH的相对丰度则与海拔高度呈正比。在南极,仅merAznuAzntA的相对丰度随海拔变化呈现出规律性变化,即随海拔的升高而增加。因此,merA在德普昌达克冰川随海拔的变化趋势与中山站相一致,而znuA随海拔的变化趋势却与中山站相反。

    HMRGs的相对丰度在德普昌达克冰川和中山站这2个位点之间也有差异。如,terDmntHpcoAznuAchrAzntAczcA在德普昌达克冰川的相对丰度平均值低于中山站,而merAczcD则相反。

  • 3个intI在2个研究位点均被检出。图1可见,3个intI在德普昌达克冰川均低于中山站。intI3在2个位点不因海拔的增加而递增或递减,intI1在德普昌达克冰川不随海拔的增加而递增或递减,但在中山站随海拔升高而递减。intI2在德普昌达克冰川随海拔升高而降低,但在中山站则相反,即随海拔升高而递增。

    对2种ARGs(sul2tetA)、9种HMRGs和3种intI进行相关性分析,结果见表2

    表2可知,由于每个位点只有3个样品,因此相关系数为1的2个因数才具有显著相关性。在德普昌达克冰川,sul2tetA分别与不同的HMRGs和intI相关联。sul2terDczcAintI1intI3显著相关;tetAczcDpcoAznuAintI2显著相关。在中山站,sul2tetA正相关,与二者密切相关的HMRGs和intI也相同。如merAzntAznuAintI2均与sul2tetA显著相关,intI2还同时与merAzntAznuA显著相关。

  • 自然环境中重金属对ARGs影响的研究较少,尤其是冰川这种环境,重金属元素的自然本底值很低,给其研究带来一定难度。由于HMRGs与重金属之间密切相关性,可以借助HMRGs来反应冰川环境重金属的水平。本文通过对青藏高原德普昌达克冰川和南极中山站雪样ARGs、HMRGs和intI的相对丰度进行分析,来判断远离人类的冰川环境重金属和intI如何影响ARGs。

    德普昌达克冰川的2个ARGs中,仅tetA随海拔的变化呈现规律性变化,与之变化相一致的是czcDpcoAznuAintI1intI2,这些基因均随海拔高度的增加而降低。YANG et al[39]对青藏高原湿地中ARGs的研究结果也证明了海拔越高,ARGs的总丰度越低,sul等基因的丰度越低。这些结果说明高海拔的地区因人类活动强度低,因而这些ARGs、HMRGs和intI含量通常较低,因此,ARGs、HMRGs和intI可以看作是人类活动强度或人为污染的代表。

    南极中山站的merAzntAznuAsul2tetAintI2均随海拔升高而升高。由于3个样品海拔高度在11~25 m之间,差别不大,因此其差异并非海拔所引起。但从样品1到样品3,距离中山站建筑区域的距离降低,人类活动就会越来越频繁。因此,上述6个基因的相对丰度就会随人类活动的增加而提升。方精云[40]对北极的研究也说明人类活动较为密集的区域周边的Mn含量会上升。

    虽然2个研究位点ARGs、HMRGs和intI的分布有差异,但也有共性。比如,ARGs均受Zn抗性基因znuA的影响,而且znuAintI之间显著相关,说明代表Zn的抗性基因与intI一起会影响ARGs在冰川上的分布。Zn与农业及水产养殖关系密切[4142],这说明影响ARGs分布的是与人类活动密切相关的重金属。znuAtetAsul2的影响大于pcoA,这与Zn对ARGs的影响大于Cu的报道[30]相一致。

    但是有些HMRGs在冰川上的分布可能与其他因素有关,因此不能仅依赖相关性来判断重金属对ARGs的影响。如merA在德普昌达克冰川和中山站均随海拔升高而增加,这与Hg具有“海拔效应”的报道[43]相一致。即Hg和其相对应的抗性基因merA在冰川雪样中会因海拔升高引起的温度降低而提高其沉积效率。Cr分布在较大的微粒上[44],而微粒的含量在冰川表层雪中随海拔升高,总体上呈明显的递减趋势[45]。因此,chrA在德普昌达克冰川随海拔的增加而降低。merAchrA在冰川上的分布或许可以反映其他环境因子的变化。

    ARGs、HMRGs和intI不仅具有海拔上的差异,还有空间上的区别。sul2tetA、3个intIzntApcoAznuAterDchrAczcAmntH在德普昌达克冰川的相对丰度平均值低于中山站。这与青藏高原冰川表层雪的Sb浓度[46]、可培养细菌菌落形成单位[47]比南极高的报道均相反。这可能是由于中山站采样位点距离科考站比较近,科考队员的频繁活动增加了ARGs、HMRGs和intI的丰度。人为影响程度越深的地方,ARGs浓度和多样性就越高,这种现象在垃圾填埋场[48]、河流、河口和海岸[49]等地均有报道。

    在远离人类的环境,HMRGs与ARGs的相关性说明了重金属对ARGs的影响。冰川环境的选择压力包括永久性的低温、高剂量的紫外线照射、每日冻结以及通常累积并存储在冰晶洞中的重金属和放射性核素[50],所有这些都可能增强ARGs。因此,要想降低ARGs的产生和扩散,必须从抗生素和重金属等多方面去考虑。tetAsul2的抗性机制分别为外排泵和靶标替代,二者在中山站均与intI2相关。intI不仅可以影响tetAsul2的扩散,还可以使2种抗性机制不相同的ARGs共同出现在同一种细菌,从而给抗生素药物的治疗带来一定难度。

  • 本研究通过qPCR,分析了青藏高原德普昌达克冰川和南极中山站不同海拔及不同研究位点ARGs、HMRGs和intI相对丰度的差异。结果表明,德普昌达克冰川ARGs、HMRGs和intI的相对丰度并不比南极中山站高。2个冰川的研究都说明,虽然冰川环境远离人类活动,但ARGs受intI的影响,同时还与人类活动有关的重金属Zn有关。然而,上述研究仅基于qPCR及相关性分析,还需用高通量qPCR以及重金属分析的结果,判断重金属对ARGs的影响,以及从序列相似性角度判断冰川环境的ARGs、HMRGs和intI是原有的还是人类活动输入的。

参考文献 (50)

目录

/

返回文章
返回