CN112226524A - 判别土壤中参与硝酸盐依赖性锑氧化过程的菌种及其关键功能基因的方法 - Google Patents

判别土壤中参与硝酸盐依赖性锑氧化过程的菌种及其关键功能基因的方法 Download PDF

Info

Publication number
CN112226524A
CN112226524A CN202010938608.4A CN202010938608A CN112226524A CN 112226524 A CN112226524 A CN 112226524A CN 202010938608 A CN202010938608 A CN 202010938608A CN 112226524 A CN112226524 A CN 112226524A
Authority
CN
China
Prior art keywords
nitrate
antimony
soil
group
iii
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.)
Granted
Application number
CN202010938608.4A
Other languages
English (en)
Other versions
CN112226524B (zh
Inventor
孙蔚旻
张苗苗
李宝琴
宋本如
林瀚智
马克斯·科尔顿
高品
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Institute of Eco Environmental and Soil Sciences of Guangdong Academy of Sciens
Original Assignee
Institute of Eco Environmental and Soil Sciences of Guangdong Academy of Sciens
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Institute of Eco Environmental and Soil Sciences of Guangdong Academy of Sciens filed Critical Institute of Eco Environmental and Soil Sciences of Guangdong Academy of Sciens
Priority to CN202010938608.4A priority Critical patent/CN112226524B/zh
Publication of CN112226524A publication Critical patent/CN112226524A/zh
Application granted granted Critical
Publication of CN112226524B publication Critical patent/CN112226524B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • CCHEMISTRY; METALLURGY
    • C12BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
    • C12QMEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
    • C12Q1/00Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions
    • C12Q1/68Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions involving nucleic acids
    • C12Q1/6876Nucleic acid products used in the analysis of nucleic acids, e.g. primers or probes
    • C12Q1/6888Nucleic acid products used in the analysis of nucleic acids, e.g. primers or probes for detection or identification of organisms
    • C12Q1/689Nucleic acid products used in the analysis of nucleic acids, e.g. primers or probes for detection or identification of organisms for bacteria
    • CCHEMISTRY; METALLURGY
    • C12BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
    • C12QMEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
    • C12Q1/00Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions
    • C12Q1/68Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions involving nucleic acids
    • C12Q1/6844Nucleic acid amplification reactions
    • C12Q1/6851Quantitative amplification
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16BBIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
    • G16B25/00ICT specially adapted for hybridisation; ICT specially adapted for gene or protein expression
    • G16B25/20Polymerase chain reaction [PCR]; Primer or probe design; Probe optimisation
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16BBIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
    • G16B30/00ICT specially adapted for sequence analysis involving nucleotides or amino acids
    • G16B30/10Sequence alignment; Homology search
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16BBIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
    • G16B30/00ICT specially adapted for sequence analysis involving nucleotides or amino acids
    • G16B30/20Sequence assembly
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16BBIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
    • G16B40/00ICT specially adapted for biostatistics; ICT specially adapted for bioinformatics-related machine learning or data mining, e.g. knowledge discovery or pattern finding
    • CCHEMISTRY; METALLURGY
    • C12BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
    • C12QMEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
    • C12Q2600/00Oligonucleotides characterized by their use
    • C12Q2600/158Expression markers
    • CCHEMISTRY; METALLURGY
    • C12BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
    • C12QMEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
    • C12Q2600/00Oligonucleotides characterized by their use
    • C12Q2600/178Oligonucleotides characterized by their use miRNA, siRNA or ncRNA

Landscapes

  • Life Sciences & Earth Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Health & Medical Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Proteomics, Peptides & Aminoacids (AREA)
  • Biophysics (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Biotechnology (AREA)
  • Analytical Chemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • Organic Chemistry (AREA)
  • Medical Informatics (AREA)
  • Theoretical Computer Science (AREA)
  • Zoology (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Wood Science & Technology (AREA)
  • Evolutionary Biology (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Molecular Biology (AREA)
  • Genetics & Genomics (AREA)
  • Immunology (AREA)
  • General Engineering & Computer Science (AREA)
  • Biochemistry (AREA)
  • Microbiology (AREA)
  • Chemical Kinetics & Catalysis (AREA)
  • Artificial Intelligence (AREA)
  • Bioethics (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Data Mining & Analysis (AREA)
  • Databases & Information Systems (AREA)
  • Epidemiology (AREA)
  • Evolutionary Computation (AREA)
  • Public Health (AREA)
  • Software Systems (AREA)
  • Measuring Or Testing Involving Enzymes Or Micro-Organisms (AREA)

Abstract

本发明公开了一种判别土壤中参与硝酸盐依赖性锑氧化过程的菌种及其关键功能基因的方法,在消耗土壤自身的底物后,加入HCO3 作为代谢底物,提供Sb(III)作为电子供体,NO3 作为电子受体。微生物依赖硝酸盐进行代谢的同时耦联驱动了锑的氧化。本发明利用DNA‑SIP技术识别能够在培养体系中驱动Sb(III)氧化过程的微生物系统发育信息,再利用宏基因组‑分箱组装技术评估Sb(III)氧化微生物的代谢潜力,分析被DNA‑SIP确定为参与了锑氧化过程的微生物的关键功能基因。本发明方法对于廓清稻田土壤锑氧化功能微生物类群及关键功能微生物的代谢信息,对于了解锑氧化过程、锑氧化菌的关键功能基因有重大的意义。

Description

判别土壤中参与硝酸盐依赖性锑氧化过程的菌种及其关键功 能基因的方法
技术领域
本发明属于土壤微生物生态学领域,具体涉及一种利用稳定同位素示踪-宏基因组-分箱组装联用平台,判别土壤中参与硝酸盐依赖性锑氧化过程的菌种及其关键功能基因的方法。
背景技术
锑已被列为疑似致癌物,长期接触可对包括心、肺、肝和肾在内的各种器官系统造成损害。虽然锑在自然界中广泛存在,但采矿、冶炼和工业等人为活动大大增加了锑的释放,从而导致环境中的锑污染。
锑的两种氧化状态,即锑铁矿(Sb(III))和锑酸盐(Sb(V)),是自然界中的主要形态,强烈地影响了Sb的迁移率、生物利用度和毒性。例如,Sb(III)比Sb(V)毒性大。因此,Sb(III)氧化可以通过将毒性较大的Sb(III)氧化为毒性较小的Sb(V)来降低其毒性,这对于减少Sb污染的生态毒性具有重要意义。
与广泛的生物地球化学循环研究相比,微生物介导的Sb循环近年来成为关注热点。微生物对Sb(III)的氧化可能是一种有效的锑污染生物修复手段。因此,了解锑转化的微生物过程和相关功能基因具有重要的环境意义,但是目前对Sb(III)氧化微生物的了解还十分有限。尤其是与好氧Sb(III)氧化相比,厌氧Sb(III)氧化研究较少,因此其机理在很大程度上尚不清楚。在稻田、湿地等自然缺氧环境中,厌氧Sb(III)氧化比好氧Sb(III)氧化更为普遍。因此,对厌氧锑氧化菌的研究和鉴定对这一生物地球化学过程有重要意义。
中国专利申请CN 201910539841.2公开了一种判别土壤中参与锑还原过程的菌种及其关键功能基因的方法,该发明利用13C-已酸和Sb(V)对采集的水稻土样品进行微宇宙培育,利用超高速离心机对13C-已酸和Sb(V)培育土壤的微生物总DNA进行离心分层,对分层后多个浮力密度中的Sb(V)还原菌基因进行实时定量PCR分析和微生物群落分析,判断该水稻土中是否含有代谢活性的已酸利用型Sb(V)还原菌。
上述专利的缺陷在于:该发明利用12C-已酸和Sb(V)及水稻土样品建立微宇宙体系,并待体系内所有Sb(V)还原成Sb(III)后进行多次稀释富集培养,对富集培养体系中土壤的总DNA建立宏基因组文库及分箱组装分析,由于进行了多次传代富集后微生物群落与原位环境存在差异,对原位环境中的功能微生物及相应的代谢通路的揭示有限。
发明内容
考虑到稻田中普遍存在的缺氧条件和高浓度的硝酸盐,硝酸盐依赖性锑氧化菌可能对Sb在这些环境中的循环起着重要作用。DNA稳定同位素示踪(SIP)技术可以在不进行分离培养的情况下将微生物特性与环境样品中的功能联系起来,从而能够更好地了解稻田土壤中与硝酸盐依赖性锑氧化生物过程有关的功能微生物群落。另外,由于现有的功能基因数据库中暂无与锑氧化功能相关的功能基因,但由于锑和砷拥有相似的化学结构,也有多项研究提出细菌可以使用相似的代谢途径来转化锑和砷,因此可以以砷氧化基因作为参考分析对象。本发明的目的在于利用稳定同位素示踪-宏基因组-分箱组装联用平台,判别土壤中参与硝酸盐依赖性锑氧化过程的菌种,揭示驱动锑氧化过程发生的关键功能基因,同时也为将来鉴定参与其他重(类)金属转化微生物的研究提供了一种先进有效的方法。
本发明的目的通过下述技术方案实现:
一种判别土壤中参与硝酸盐依赖性锑氧化过程的菌种的方法,通过微宇宙培养体系,经过饥饿培养消耗土壤自身的底物后,在体系中加入HCO3 -作为代谢底物,以及向体系提供Sb(III)作为电子供体,NO3 -作为电子受体。使培养体系中只存在一个主导的电子交换过程,即以加入的Sb(III)作为电子供体,以NO3 -作为电子受体,微生物依赖硝酸盐进行代谢的同时耦联驱动了锑的氧化,使Sb(III)失去电子氧化成Sb(V)。
本发明对受Sb(III)胁迫下稻田土壤厌氧培养体系的Sb(III)氧化现象进行观察,利用DNA-SIP技术来识别能够在该培养体系中驱动Sb(III)氧化过程的微生物的系统发育信息,再利用宏基因组-分箱组装技术评估Sb(III)氧化微生物的代谢潜力,并重点分析被DNA-SIP确定为参与了锑氧化过程的微生物的关键功能基因。
具体地,本发明的方法包括以下步骤:
(1)采集存在锑污染的土壤,加入到无机盐培养基中,厌氧培养,直至消耗完土壤本底有机物;然后分成三组微宇宙体系,向第一组(13C+N+Sb)微宇宙体系加入13C-HCO3 -、硝酸盐和Sb(III)进行培养,以标记并观察具有硝酸盐还原耦联锑氧化能力的微生物;向第二组(12C+N+Sb)微宇宙体系加入12C-HCO3 -、硝酸盐和Sb(III)进行培养,以进一步确认体系中具有硝酸盐还原耦联锑氧化能力的微生物,排除第一组中因为其他原因(如ATCG比例)存在于重质组分的干扰微生物;向第三组(13C+N)微宇宙体系加入13C-HCO3 -和硝酸盐进行培养,以排除没有参与锑氧化代谢的微生物;三组微宇宙体系培养至少60天;
在第一组(13C+N+Sb)中利用13C-HCO3 -的硝酸盐依赖性锑氧化微生物的DNA会因为13C的标记而变重,而第二组(12C+N+Sb)中利用12C-HCO3 -的硝酸盐依赖性锑氧化微生物的DNA会因为无13C的标记存在于轻质组分(用来排除第一组中因为其他原因(如ATCG比例)存在于重质组分的干扰微生物),第三组(13C+N)是为了排除体系中(第一组中)能利用碳酸氢钠和硝酸盐但没有耦合锑氧化能力的菌(会有不少种类必须排除)。
所述的土壤是自然缺氧环境下的土壤,比如稻田或湿地土壤;
所述无机盐培养基的组成是:10.55g/L Na2HPO4·12H2O,1.5g/L KH2PO4,0.3g/LNH4Cl,0.1g/L MgCl2,0.01mg/L维生素H,0.02mg/L烟酸,0.1mg/L维生素B1,0.01mg/L对氨基苯甲酸,0.005mg/L维生素B5,0.05mg/L盐酸吡哆胺,0.01mg/L氰钴胺,10μL/L HCl(25%,w/w),1.5mg/L FeCl2·4H2O,0.19mg/L CoCl2·6H2O,0.1mg/L MnCl2·2H2O,0.07mg/L ZnCl2,0.024mg/L NiCl2·6H2O,0.036mg/L NaMoO4·2H2O,0.006mg/L H3BO3,0.002mg/L CuCl2·2H2O;
所述的厌氧培养,是培养过程中向培养体系做N2吹扫;所述的厌氧培养不少于一个月;
优选地,所述各组微宇宙体系中,13C-HCO3 -的终浓度为0.3-0.6g/L,12C-HCO3 -的终浓度为0.3-0.6g/L,硝酸盐的终浓度为0.1-0.3g/L,Sb(III)的终浓度为0.1-0.2g/L;培养过程中定期(优选每12天)添加以上物质,强化促进微生物代谢,提高后续提取的含13C的DNA产量,以便直接进行宏基因组测序;
所述的13C-HCO3 -12C-HCO3 -优选采用钠盐;所述的硝酸盐优选硝酸钠;所述的Sb(III)优选采用C8H4K2O12Sb2·3H2O;
(2)对培养后的三组微宇宙体系,提取其中的土壤微生物总DNA,将所得DNA提取物加入CsCl溶液做超高速离心,收集各层离心组分;
步骤(2)所述的超高速离心是以178,000×g离心48小时;
(3)测定各层离心组分的浮力密度(Buoyant Density,BD)值,以BD值将三组微宇宙体系的土壤微生物总DNA分成20-25个分层,去除各分层中的CsCl,并对各层离心组分做q-PCR扩增;
以各层离心组分的BD值为横坐标,以功能基因aioA的q-PCR扩增数值为纵坐标,作图;第一组(13C+N+Sb)微宇宙体系形成1个峰,BD值在1.70-1.76,为重质组分;第二组(12C+N+Sb)微宇宙体系形成1个峰,BD值在1.68-1.72,为轻质组分;第三组(13C+N)微宇宙体系形成1个峰,BD值在1.68-1.72,为轻质组分;
在三组微宇宙体系的重质或轻质组分中各选择1-3个q-PCR扩增数值较高的分层,进行16S rRNA基因V4-V5区域的高通量测序;
(4)将得到的16S rRNA V4-V5区域测序数据与已有的16S rRNA数据库进行比对,将具有97%以上相似度的划归为一个小组(OTU),分析各个OTU在三组微宇宙体系中的丰度;
每个小组(OTU)对应于一个不同的16S rRNA序列,也即对应于一种不同的细菌(微生物);
(5)在微生物群落测序结果中关注丰度较高的小组(OTUs),通过以下步骤确定代表硝酸盐依赖性锑氧化的微生物:
第一步,排除(13C+N)组中在重质组分明显富集的OTUs,由于(13C+N)组没有加入Sb(III),说明这些在重质组分富集的OTUs为其他同化13C-HCO3 -的微生物,而非参与锑氧化代谢的微生物;
第二步,在(12C+N+Sb)组的轻质组分富集,且在(13C+N+Sb)组重质组分富集的OTUs,就是具有硝酸盐依懒性锑氧化能力的微生物;这些OTUs所代表的微生物在微宇宙系统中发生了硝酸盐还原耦联锑氧化的过程,(13C+N+Sb)组中微生物同化了体系中的底物13C-HCO3 -13C使其DNA“变重”,由(12C+N+Sb)组的轻质组分“移动”到(13C+N+Sb)组的重质组分,具有以上特性的OTUs则被判断为具有硝酸盐依懒性锑氧化能力的微生物。
一种判别土壤中参与硝酸盐依赖性锑氧化过程的关键功能基因的方法,包括以下步骤:
(1)建立第一组(13C+N+Sb)微宇宙体系重质组分DNA的宏基因组文库,获得原始测序Reads,进行测序数据质控,过滤低质量数据,然后进行序列拼接,得到Contigs,再进行序列比对,把Contigs的独立数据集的Reads进行映射用来评估其丰度,再对Contigs进行Binning组装,取完整度>90%和冗余度<10%的分箱(bins)进行下游分析;
所述建立宏基因组文库,优选使用Illumina Hiseq 4000平台建立;
所述的数据质控,优选使用Trimmomatic-0.36进行分析;
所述的序列拼接,优选使用Megahit进行;
所述的序列比对,优选使用Bowtie 2进行;
所述Binning组装,优选使用默认设置的CONCOCT(版本0.4.0)进行;
(2)在宏基因组分箱(bins)中分析与锑氧化、氮循环功能相关的基因的丰度大小,由此获得与Sb(III)氧化代谢能力相关的基因,并获知在锑氧化过程中与N元素耦联作用相关的关键基因,由此更进一步了解稻田土壤中微生物驱动的硝酸盐耦联锑氧化过程的生物地球化学机制。
本发明相对于现有技术具有如下的优点及效果:
本发明利用的稳定同位素示踪-宏基因组-分箱组装联用平台,能够不依赖分离或富集培养准确判别土壤中参与了重(类)金属转化过程的微生物。利用DNA-SIP直接锚定群落中的关键功能微生物,避免了多次传代富集后微生物群落与原位环境存在差异,能直接对DNA-SIP锚定的功能微生物DNA进行宏基因组分析,进一步依靠分箱组装方法,廓清稻田土壤锑氧化功能微生物类群及关键功能微生物的代谢信息,对于了解微生物驱动的锑氧化过程、锑氧化菌以及关键功能基因的认知有着重大的意义。
附图说明
图1为SIP微宇宙培育体系的Sb(III)(a)、Sb(V)(b)、NO3 -(c)和NO2 -(d)浓度变化示意图。
图2为不同培养时间13/12C+N+Sb和13/12C+N处理的DNA在不同浮力密度组分中aioA基因的相对丰度。
图3为SIP微宇宙培育体系中具有不同的浮力密度DNA组分的OTU相对丰度气泡图(a)及在13/12C+N+Sb和13/12C+N处理中丰度线性判别分析(b)。
图4为具有锑氧化潜能微生物的基因组bins和read深度的体统发育从属关系示意图。
图5为锑氧化微宇宙体系的DNA宏基因组-分箱组装分析中,与锑氧化、氮循环基因相关的热图。
具体实施方式
下面结合实施例及附图对本发明作进一步详细的描述,但本发明的实施方式不限于此。
实施例
一种判别土壤中参与硝酸盐依赖性锑氧化过程的菌种的方法,包括以下步骤:
(1)样品采集和处理
土壤样本采集于湖南锡矿山锑矿附近,长期受锑污染的稻田土壤已经形成对微生物群选择性,可能富集了锑代谢微生物,在地表30-40厘米深的稻田土壤中采集样本,低温保存运送至实验室。
(2)建立DNA-SIP微宇宙培育体系
用采集回来的稻田土壤样品建立三组微宇宙体系,微宇宙体系用50mL的灭菌血清瓶建立,在瓶内加入约2g土壤以及40mL无机盐培养基(Mineral Salts Medium,MSM),在瓶内注入N2吹扫使微宇宙体系保持厌氧状态;
MSM溶液组成:10.55g/L Na2HPO4·12H2O,1.5g/L KH2PO4,0.3g/L NH4Cl,0.1g/LMgCl2,0.01mg/L维生素H,0.02mg/L烟酸,0.1mg/L维生素B1,0.01mg/L对氨基苯甲酸,0.005mg/L维生素B5,0.05mg/L盐酸吡哆胺,0.01mg/L氰钴胺,10μL/L HCl(25%,w/w),1.5mg/L FeCl2·4H2O,0.19mg/L CoCl2·6H2O,0.1mg/L MnCl2·2H2O,0.07mg/L ZnCl2,0.024mg/L NiCl2·6H2O,0.036mg/L NaMoO4·2H2O,0.006mg/L H3BO3,0.002mg/L CuCl2·2H2O;
经过一个月的饥饿培养,消耗土壤本底底物,然后分成三组:
向第一组(13C+N+Sb)微宇宙体系加入0.672g/L(终浓度)13C-NaHCO3、0.255g/L(终浓度)的NaNO3和0.632g/L(终浓度)C8H4K2O12Sb2·3H2O进行培育,向第二组(12C+N+Sb)微宇宙培育体系加入0.672g/L(终浓度)12C-NaHCO3、0.255g/L(终浓度)的NaNO3和0.632g/L(终浓度)C8H4K2O12Sb2·3H2O进行培育,向第三组(13C+N)微宇宙培育体系加入0.672g/L(终浓度)13C-NaHCO3和0.255g/L(终浓度)的NaNO3进行培育,并每12天对相应微宇宙体系原料物质进行再次补充。
在培育的第0天、30天和60天对三组微宇宙体系取样,用土壤DNA提取试剂盒进行土壤微生物总DNA的提取。
另外,利用连续流动分析仪(SAN++,Skalar,Breda,Holland)对三组微宇宙体系上清液中的NO3 -和NO2 -浓度进行定量,用高效液相色谱-氢化物发生-原子荧光分析仪(HPLC-HG-AFS)对体系整个培养过程中溶液的Sb(III)和Sb(V)浓度进行测定,观察到只有同时添加了Sb(III)和NO3 -的微宇宙体系里存在Sb(III)氧化成Sb(V)的现象,即Sb(III)浓度下降,Sb(V)浓度上升;同时也出现NO3 -还原为NO2 -现象,即NO3 -浓度下降,NO2 -浓度上升(图1)。而没有添加Sb(III)或NO3 -或灭菌的对照组中没有发生以上现象,说明同时添加了Sb(III)和NO3 -的活性微宇宙中发生了微生物驱动的硝酸盐依赖性锑氧化过程。
(3)利用超高速离心机对DNA-SIP培育体系的土壤微生物总DNA进行离心分层
取10μg DNA提取物放入5.1mL超速离心专用快封管中,加入CsCl溶液至接近装满离心管,并用Tris-EDTA(pH值8.0)和CsCl溶液调整离心管内BD值至1.73g/mL(BD值用折光仪测量),随后封管。把离心管放入超速离心机中以178,000×g的离心力在20℃下超速离心48小时,取出离心管放入组分回收装置中,用固定流速泵对组分回收装置中的离心管混合液进行分层收集回收组分,每层大概210μL共24层。
(4)对每个DNA分层组分进行BD值测定以及PCR扩增
测定各层回收组分的BD值,然后用核酸助沉剂和乙醇沉淀去除CsCl,得到净化后的回收组分。
由于锑和砷拥有相似的化学结构,细菌可以使用相似的代谢途径来转化锑和砷,目前还没有明确的锑氧化基因,因此以砷氧化基因作为参考分析对象。利用已知的砷氧化基因引物aioA F(5’-CCA CTT CTG CAT CGT GGG NTG YGG NTA-3’)和aioA R(5’-GGA GTTGTA GGC GGG CCK RTT RTG DAT-3’)对回收组分进行q-PCR扩增。结果如图2所示:
第一组(13C+N+Sb)微宇宙体系形成1个峰,BD值在1.70-1.76,为重质组分;第二组(12C+N+Sb)微宇宙体系形成1个峰,BD值在1.68-1.72,为轻质组分;第三组(13C+N)微宇宙体系形成1个峰,BD值在1.68-1.72,为轻质组分;
(5)挑选多个浮力密度组分DNA进行16S rRNA测序
根据BD值及功能基因q-PCR的作图确定DNA重质和轻质组分在全组分分层中的分布,在第一组(13C+N+Sb)和第三(13C+N)组微宇宙体系的重质组分中分别选择3个分层,第二(12C+N+Sb)组微宇宙体系轻质组分中选择3个分层,对以上分层进行16S rRNA基因V4-V5区域的高通量测序。
将得到的测序数据与已有的16S rRNA数据库进行比对分析以及归类操作,将序列按照彼此的相似性分归为许多小组,一个小组就是一个操作分类单元(OTU),通常按97%以上的相似度,每个OTU对应于一个不同的16S rRNA序列,也就是每个OTU对应于一个不同的细菌(微生物)种。
(6)通过DNA-SIP对稻田土壤中锑氧化活性细菌群落的分析
通过OTU分析,分析样品中的微生物群落多样性和不同微生物群种的丰度。图3为群落中丰度最高的前25个OTU。
在(13C+N+Sb)组微宇宙体系中13C-DNA富集的OTU可能具有硝酸盐还原耦联Sb(III)氧化的能力。然而,由于NaHCO3还可以被很多其他的微生物同化,因此需要对比(13C+N+Sb)组与(13C+N)组微宇宙体系,以及对比(13C+N+Sb)组与(12C+N+Sb)组微宇宙体系排除干扰,来判断参与Sb(III)氧化过程的细菌。
由于体系中只存在硝酸盐还原耦联锑氧化为主导的电子交换过程,锑氧化体系中的参与氧化过程的微生物利用了13C-NaHCO3作为底物进行代谢,使其DNA得到了13C的标记,在DNA超速离心过程中13C在重质组分中聚集,因此在第一组(13C+N+Sb)微宇宙体系重质组分中富集的OTUs系统发育信息所代表的微生物可能具有硝酸盐还原耦联锑氧化发生的能力;而在第二组(12C+N+Sb)微宇宙培育体系中,由于加入的是12C,因此硝酸盐还原耦联锑氧化发生的微生物OTUs则在轻质组分中富集;在第三组(13C+N)微宇宙培养体系重组分中富集的OTUs被判断能同化NaHCO3的微生物的系统发育信息,但由于体系中不含有锑,因此这些微生物不一定能发生硝酸盐还原耦联锑氧化反应。
判断方法为:
第一步,排除(13C+N)组在重质组分明显富集的OTUs,由于(13C+N)组没有加入Sb(III),这些在重质组分富集的OTUs为其他同化13C-NaHCO3的微生物,而非参与锑氧化代谢的微生物;
第二步,排除上述干扰后继续观察,OTUs在(12C+N+Sb)组的轻质组分富集,且在(13C+N+Sb)组的重质组分富集,说明这些OTUs所代表的微生物在微宇宙系统中发生了硝酸盐还原耦联锑氧化的过程,(13C+N+Sb)组中微生物同化了的NaHCO3因为13C使NaHCO3“变重”,所以由(12C+N+Sb)组的轻质组分“移动”到(13C+N+Sb)组的重质组分,具有以上特性的OTUs则被判断为具有硝酸盐还原耦联Sb(III)氧化能力的微生物。
如图3所示,根据上述判断方法,对比微生物在(13C+N+Sb)、(12C+N+Sb)、(13C+N)三组微宇宙体系的DNA-SIP群落OTU结果,判断出Azoarcus sp.、Azospira sp.和Chelativorans sp.具有上述特征,对这些锑氧化代谢菌在DNA浮力密度组分中的丰度分布进行分析,发现硝酸盐依赖性锑氧化体系中这3种锑氧化菌都同化了13C-NaHCO3,使得自身DNA“变重”,从而出现在浮力密度重层,因此在浮力密度重层中的这些锑氧化菌的分布和轻层明显不同。因此推断为参与锑氧化代谢的细菌。
一种判别土壤中参与硝酸盐依赖性锑氧化过程的关键功能基因的方法,包括以下步骤:
(1)重质组分DNA样品的宏基因组分析,并进行分箱组装
在Illumina Hiseq 4000平台建立含锑氧化功能微生物的(13C+N+Sb)组重质组分DNA的宏基因组文库,获得原始测序Reads,利用Trimmomatic-0.36进行测序数据质控,过滤低质量数据,然后使用Megahit进行序列拼接,得到Contigs,再使用Bowtie 2进行序列比对,把Contigs的独立数据集的Reads进行映射用来评估其丰度,再对Contigs使用默认设置的CONCOCT(版本0.4.0)进行Binning组装,取完整度>90%和冗余度<10%的bins进行下游分析。
(2)宏基因组-分箱组装揭示锑氧化代谢相关的功能基因的分析
利用宏基因组分析研究Sb(III)氧化微生物群落的代谢潜力,特别针对那些通过DNA-SIP鉴定的锑氧化菌。
经过分析,从宏基因组数据中获得了15个高质量的bins,这些bins属于Proteobacteria、Gemmatimonadetes和Verrucomicrobia三个门类(图4),在这些bins中也检测到在DNA-SIP中推断的锑氧化菌(完整度>95%,受污染程度<5%),如Azoarcus sp.(bin 9)、Azospira sp.(bin 10)和Chelativorans sp.(bin 14)(图4)。
在宏基因组bins中分析与氮循环、锑氧化等相关的功能注析基因的丰度大小(图5),另因现有数据库中没有与锑循环相关的功能基因,但由于锑和砷拥有相似的化学结构,已经有研究提出细菌可以使用相似的代谢途径来转化锑和砷,因此以砷循环基因作为参考分析对象。由此获得与Sb(III)氧化代谢能力相关的基因,并获知在锑氧化过程中与氮等元素耦联作用的生物地球化学循环相关关键基因,由此更进一步了解稻田土壤中微生物驱动的硝酸盐耦联锑氧化过程的生物地球化学机制。
锑氧化基因分析:由图5的分箱遗传注释表明,大多数分箱(除与Devosia相关的bin 6和与Gemmatimons相关的bin 8及bin 15外)含有aioA基因,而仅在三个分箱(与Thauera相关的bin 4、与Vulgatibacterae相关的bin 12和与Chelativorans相关的bin14)中观察到aioB基因(图5)。
锑氧化耦联的氮循环基因分析:对与反硝化作用有关的基因(包括负责NO3 -、NO2 -、NO-和N2O还原的基因)分析发现:总体来说硝酸盐还原基因(narG、narZ、narI、napA等)是所有分箱中检测到的最丰富的基因,其次是NO2 -还原基因(nirK、nirS等)和NO还原基因(norB、norC等),而编码N2O还原酶的基因(nosZ)在8个分箱中检测到(图5)。其中,Ramlibacter相关联的bin 3、Ramlibacter相关联的bin 9、Azospira相关联的bin 10和Anaeromyxobacter相关联的bin 13被观察到包含了从NO3 -到N2整个反硝化过程的驱动基因(图5)。
此结果揭示了稻田土壤中潜在的硝酸盐依赖性锑氧化菌及其功能基因,拓展了目前对稻田土壤锑氧化菌的生态功能以及微生物驱动的锑的地球化学循环过程的认知。
上述实施例为本发明较佳的实施方式,但本发明的实施方式并不受上述实施例的限制,其他的任何未背离本发明的精神实质与原理下所作的改变、修饰、替代、组合、简化,均应为等效的置换方式,都包含在本发明的保护范围之内。

Claims (10)

1.一种判别土壤中参与硝酸盐依赖性锑氧化过程的菌种的方法,其特征在于包括以下步骤:
(1)采集存在锑污染的土壤,加入到无机盐培养基中,厌氧培养,直至消耗完土壤本底有机物;然后分成三组微宇宙体系,向第一组(13C+N+Sb)微宇宙体系加入13C-HCO3 -、硝酸盐和Sb(III)进行培养,向第二组(12C+N+Sb)微宇宙体系加入12C-HCO3 -、硝酸盐和Sb(III)进行培养,向第三组(13C+N)微宇宙体系加入13C-HCO3 -和硝酸盐进行培养,三组微宇宙体系培养至少60天;
(2)对培养后的三组微宇宙体系,提取其中的土壤微生物总DNA,将所得DNA提取物加入CsCl溶液做超高速离心,收集各层离心组分;
(3)测定各层离心组分的BD值,以BD值将三组微宇宙体系的土壤微生物总DNA分成20-25个分层,去除各分层中的CsCl,并对各层离心组分做q-PCR扩增;
以各层离心组分的BD值为横坐标,以功能基因aioA的q-PCR扩增数值为纵坐标,作图;第一组(13C+N+Sb)微宇宙体系形成1个峰,BD值在1.70-1.76,为重质组分;第二组(12C+N+Sb)微宇宙体系形成1个峰,BD值在1.68-1.72,为轻质组分;第三组(13C+N)微宇宙体系形成1个峰,BD值在1.68-1.72,为轻质组分;
在三组微宇宙体系的重质或轻质组分中各选择1-3个q-PCR扩增数值较高的分层,进行16S rRNA基因V4-V5区域的高通量测序;
(4)将得到的16S rRNA V4-V5区域测序数据与已有的16S rRNA数据库进行比对,将具有97%以上相似度的划归为一个OTU,分析各个OTU在三组微宇宙体系中的丰度;
(5)在微生物群落测序结果中关注丰度较高的OTUs,通过以下步骤确定代表硝酸盐依赖性锑氧化的微生物:
第一步,排除(13C+N)组中在重质组分明显富集的OTUs;
第二步,在(12C+N+Sb)组的轻质组分富集,且在(13C+N+Sb)组重质组分富集的OTUs,就是具有硝酸盐依懒性锑氧化能力的微生物。
2.根据权利要求1所述的方法,其特征在于:步骤(1)所述各组微宇宙体系中,13C-HCO3 -的终浓度为0.3-0.6g/L,12C-HCO3 -的终浓度为0.3-0.6g/L,硝酸盐的终浓度为0.1-0.3g/L,Sb(III)的终浓度为0.1-0.2g/L;培养过程中定期添加以上物质。
3.根据权利要求1所述的方法,其特征在于:步骤(1)所述无机盐培养基的组成是:10.55g/L Na2HPO4·12H2O,1.5g/L KH2PO4,0.3g/L NH4Cl,0.1g/L MgCl2,0.01mg/L维生素H,0.02mg/L烟酸,0.1mg/L维生素B1,0.01mg/L对氨基苯甲酸,0.005mg/L维生素B5,0.05mg/L盐酸吡哆胺,0.01mg/L氰钴胺,10μL/L HCl(25%,w/w),1.5mg/L FeCl2·4H2O,0.19mg/LCoCl2·6H2O,0.1mg/L MnCl2·2H2O,0.07mg/L ZnCl2,0.024mg/L NiCl2·6H2O,0.036mg/LNaMoO4·2H2O,0.006mg/L H3BO3,0.002mg/L CuCl2·2H2O。
4.根据权利要求1所述的方法,其特征在于:步骤(1)所述的厌氧培养,是培养过程中向培养体系做N2吹扫;所述的厌氧培养不少于一个月。
5.根据权利要求1所述的方法,其特征在于:步骤(1)所述的土壤是自然缺氧环境下的土壤。
6.根据权利要求1所述的方法,其特征在于:步骤(1)所述的13C-HCO3 -12C-HCO3 -采用钠盐;所述的硝酸盐为硝酸钠;所述的Sb(III)采用C8H4K2O12Sb2·3H2O。
7.根据权利要求1所述的方法,其特征在于:步骤(2)所述的超高速离心是以178,000×g离心48小时。
8.一种判别土壤中参与硝酸盐依赖性锑氧化过程的关键功能基因的方法,其特征在于包括以下步骤:
(1)建立第一组(13C+N+Sb)微宇宙体系重质组分DNA的宏基因组文库,获得原始测序Reads,进行测序数据质控,过滤低质量数据,然后进行序列拼接,得到Contigs,再进行序列比对,把Contigs的独立数据集的Reads进行映射用来评估其丰度,再对Contigs进行Binning组装,取完整度>90%和冗余度<10%的分箱进行下游分析;
(2)在宏基因组分箱中分析与锑氧化、氮循环功能相关的基因的丰度大小,由此获得与Sb(III)氧化代谢能力相关的基因,并获知在锑氧化过程中与N元素耦联作用相关的关键基因。
9.根据权利要求8所述的方法,其特征在于:
步骤(1)所述建立宏基因组文库,使用Illumina Hiseq 4000平台建立;
步骤(1)所述的数据质控,使用Trimmomatic-0.36进行分析。
10.根据权利要求8所述的方法,其特征在于:
步骤(1)所述的序列拼接,使用Megahit进行;
步骤(1)所述的序列比对,使用Bowtie 2进行;
步骤(1)所述Binning组装,使用默认设置的CONCOCT进行。
CN202010938608.4A 2020-09-09 2020-09-09 判别土壤中参与硝酸盐依赖性锑氧化过程的菌种及其关键功能基因的方法 Active CN112226524B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010938608.4A CN112226524B (zh) 2020-09-09 2020-09-09 判别土壤中参与硝酸盐依赖性锑氧化过程的菌种及其关键功能基因的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010938608.4A CN112226524B (zh) 2020-09-09 2020-09-09 判别土壤中参与硝酸盐依赖性锑氧化过程的菌种及其关键功能基因的方法

Publications (2)

Publication Number Publication Date
CN112226524A true CN112226524A (zh) 2021-01-15
CN112226524B CN112226524B (zh) 2023-10-27

Family

ID=74116257

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010938608.4A Active CN112226524B (zh) 2020-09-09 2020-09-09 判别土壤中参与硝酸盐依赖性锑氧化过程的菌种及其关键功能基因的方法

Country Status (1)

Country Link
CN (1) CN112226524B (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113930477A (zh) * 2021-09-14 2022-01-14 广东省科学院生态环境与土壤研究所 一种检测植物根内生菌含有自养砷氧化功能微生物的方法
WO2022242427A1 (zh) * 2021-05-20 2022-11-24 北京工业大学 一种微宇宙培养装置及其在土壤碳扩散与微生物利用过程的定量分析中的应用

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102071471A (zh) * 2010-12-30 2011-05-25 南京大学 大豆根际土壤微生物群落宏基因组文库的构建及其应用
CN103966318A (zh) * 2014-04-04 2014-08-06 中国科学院南京土壤研究所 利用dna稳定性同位素探针原位揭示判别稻田甲酸利用型产甲烷古菌的方法
CN107557312A (zh) * 2017-07-20 2018-01-09 武汉纺织大学 一种硝酸盐依赖亚铁氧化菌群富集培养基及制备方法和应用
US20190034588A1 (en) * 2017-07-28 2019-01-31 Tata Consultancy Services Limited Method and system for identification and classification of operational taxonomic units in a metagenomic sample
US20190194740A1 (en) * 2016-05-19 2019-06-27 Perfect (China) Co., Ltd. Determination of microorganism operational taxonomic unit and sequence-assisted separation
CN110317863A (zh) * 2019-06-21 2019-10-11 广东省生态环境技术研究所 一种判别土壤中参与锑还原过程的菌种及其关键功能基因的方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102071471A (zh) * 2010-12-30 2011-05-25 南京大学 大豆根际土壤微生物群落宏基因组文库的构建及其应用
CN103966318A (zh) * 2014-04-04 2014-08-06 中国科学院南京土壤研究所 利用dna稳定性同位素探针原位揭示判别稻田甲酸利用型产甲烷古菌的方法
US20190194740A1 (en) * 2016-05-19 2019-06-27 Perfect (China) Co., Ltd. Determination of microorganism operational taxonomic unit and sequence-assisted separation
CN107557312A (zh) * 2017-07-20 2018-01-09 武汉纺织大学 一种硝酸盐依赖亚铁氧化菌群富集培养基及制备方法和应用
US20190034588A1 (en) * 2017-07-28 2019-01-31 Tata Consultancy Services Limited Method and system for identification and classification of operational taxonomic units in a metagenomic sample
CN110317863A (zh) * 2019-06-21 2019-10-11 广东省生态环境技术研究所 一种判别土壤中参与锑还原过程的菌种及其关键功能基因的方法

Non-Patent Citations (6)

* Cited by examiner, † Cited by third party
Title
MAX KOLTON等: "Bacteria responsible for antimonite oxidation in antimony-contaminated soil revealed by DNA-SIP coupled to metagenomics." *
MIAOMIAO ZHANG等: "Characterization of Nitrate-Dependent As(III)-Oxidizing Communities in Arsenic-Contaminated Soil and Investigation of Their Metabolic Potentials by the Combination of DNA-Stable Isotope Probing and Metagenomics" *
WEI XING等: "Identification of the autotrophic denitrifying community in nitrate removal reactors by DNA-stable isotope probing" *
XIAOMIN LI等: "Bacterial Communities and Functional Genes Stimulated During Anaerobic Arsenite Oxidation and Nitrate Reduction in a Paddy Soil" *
李爽;敖俊华;王庆;陈迪文;周文灵;吴启华;: "稻田NO3-还原耦合As(Ⅲ)氧化过程及微生物群落特征" *
王年;鲁小璐;邬梦晓俊;王猛;王红梅;: "微生物氧化As(Ⅲ)和Sb(Ⅲ)的研究进展" *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2022242427A1 (zh) * 2021-05-20 2022-11-24 北京工业大学 一种微宇宙培养装置及其在土壤碳扩散与微生物利用过程的定量分析中的应用
CN113930477A (zh) * 2021-09-14 2022-01-14 广东省科学院生态环境与土壤研究所 一种检测植物根内生菌含有自养砷氧化功能微生物的方法

Also Published As

Publication number Publication date
CN112226524B (zh) 2023-10-27

Similar Documents

Publication Publication Date Title
CN110317863B (zh) 一种判别土壤中参与锑还原过程的菌种及其关键功能基因的方法
Wang et al. Microbial insights into the biogeochemical features of thallium occurrence: a case study from polluted river sediments
Neufeld et al. Methodological considerations for the use of stable isotope probing in microbial ecology
Johnson et al. Techniques for detecting and identifying acidophilic mineral-oxidizing microorganisms
Ding et al. Investigation of microbial structure and composition involved in membrane fouling in the forward osmosis membrane bioreactor treating anaerobic bioreactor effluent
Dekas et al. Characterizing chemoautotrophy and heterotrophy in marine archaea and bacteria with single-cell multi-isotope NanoSIP
CN112226524B (zh) 判别土壤中参与硝酸盐依赖性锑氧化过程的菌种及其关键功能基因的方法
Zhang et al. Bacteria responsible for nitrate-dependent antimonite oxidation in antimony-contaminated paddy soil revealed by the combination of DNA-SIP and metagenomics
Kellermann et al. Microbial CO2 fixation potential in a tar-oil-contaminated porous aquifer
Zhao et al. Self-regulating microbiome networks ensure functional resilience of biofilms in sand biofilters during manganese load fluctuations
Tominski et al. Insights into carbon metabolism provided by fluorescence in situ hybridization-secondary ion mass spectrometry imaging of an autotrophic, nitrate-reducing, Fe (II)-oxidizing enrichment culture
Hayer et al. Identification of growing bacteria during litter decomposition in freshwater through quantitative stable isotope probing
WO2021212919A1 (zh) 一种基于绝迹稀释培养组的功能微生物高通量分析和选育方法
Abe et al. A rapid collection of yet unknown ammonia oxidizers in pure culture from activated sludge
CN112662727A (zh) 一种基于绝迹稀释培养组的功能微生物高通量分析和选育方法
White Is there anything else you need to understand about the microbiota that cannot be derived from analysis of nucleic acids?
Ito et al. Isolation, characterization, and in situ detection of a novel chemolithoautotrophic sulfur-oxidizing bacterium in wastewater biofilms growing under microaerophilic conditions
Du et al. Tracking de novo protein synthesis in the activated sludge microbiome using BONCAT-FACS
Jangir et al. In situ electrochemical studies of the terrestrial deep subsurface biosphere at the sanford underground research facility, South Dakota, USA
Zhao et al. Biofilm: A strategy for the dominance of comammox Nitrospira
Huang et al. Metagenomic analysis revealed the sulfur-and iron-oxidation capabilities of heterotrophic denitrifying sludge
CN111606371A (zh) 一种用于酸性矿水治理的新型可渗透反应墙装置和填料
Sieradzki et al. Metagenomics and stable isotope probing offer insights into metabolism of polycyclic aromatic hydrocarbons degraders in chronically polluted seawater
CN114289499A (zh) 一种基于宏基因组学的土著细菌修复石油污染土壤的方法
CN109593673B (zh) 一种黄杆菌jx-1及其在污水治理中的应用

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