一种基于绝迹稀释培养组的功能微生物高通量分析和选育方法

：
1.本发明属于微生物技术领域，具体涉及一种基于绝迹稀释培养组的功能微生物高通量分析和选育方法。


背景技术：
：
2.高通量测序技术和宏基因组学的出现和发展，使我们认识到目前已经分离培养的微生物仅仅是一小部分，环境中超过99％的微生物都仍然未实现纯培养。这些未培养的微生物虽然已经被发现，但是高通量测序获得的信息非常有限，我们很难通过核酸序列信息获知相应物种的功能和表型信息，以及它在环境中的生态作用。此外，生态系统中存在着许多低丰度的稀有物种，这些稀有物种是多样性的重要组成部分。在条件合适时，这些物种可能迅速响应，对推动元素循环和物质转化以及维持生态系统的稳定性发挥着重要的作用。但是，由于高通量测序技术自带的多个步骤的抽样过程，具有一定的检测下限，很难检测到这些稀有物种的存在，更难于研究它们的功能。因此，基于目前高通量测序等技术手段很难对环境中微生物的多样性和功能有全面的认识。大量功能未知的未培养微生物的培养难和低丰度稀有物种的测定难是限制微生物生态学发展的瓶颈。
3.近年来，培养组学概念的兴起为未培养微生物的挖掘和功能研究提供了新思路。培养组学是通过多种培养条件对目标微生物进行培养，并利用基质辅助激光解吸电离飞行时间质谱(maldi
‑
tof
‑
ms)和16s rrna基因测序技术对它们进行高通量鉴定的一种培养技术。通过培养组，许多研究在人体微生物群，特别是肠道微生物群中分离培养了大量的新型微生物，更新了我们对人体微生物组的认识。将培养组运用于微生物生态学研究中，将与高通量测序技术等非培养技术形成有效的互补，更加准确的阐明微生物生态学中的问题。培养组技术的特征在于对微生物的高通量鉴定上，而对于环境中特殊类群的微生物仍然需要结合不同的培养方法对不同的微生物种类进行分离。然而，目前对于环境中目标功能微生物类群的选育，仍然缺乏较好的方法。尤其是环境中相对丰度较低的功能物种的选育。


技术实现要素：
：
4.本发明的目的是提供一种基于绝迹稀释培养组的功能微生物高通量分析和选育方法，其能有效分离环境中不同丰度和不同功能特征的物种。
5.本发明的基于绝迹稀释培养组的功能微生物高通量分析和选育方法，包括以下步骤：
6.a、利用含底物的选择性培养基对环境样品进行多个系列的梯度稀释；
7.b、将梯度稀释的样品进行培养，并测定选择性培养基中底物的浓度变化情况；
8.c、选择降解底物最多的、活性的最高的若干个稀释度的样品组成绝迹稀释培养组，然后收集菌体，按照统一的生物量接种到新鲜的含底物的选择性培养基上进行培养，再测定底物的变化情况，利用底物的变化率反映绝迹稀释培养组不同样品的潜在活性；
9.d、收集步骤c测定潜在活性后的样品中的菌体，提取总dna，并扩增16s rrna基因v4
‑
v5区进行高通量测序，利用生物信息学方法分析绝迹稀释培养组样品微生物群落的物种组成，及其与潜在活性特征的关系；
10.e、根据绝迹稀释培养组样品的微生物群落组成特征，选取样品先用固体富集培养基对目标功能物种进行单菌落分离筛选，纯化，并进行物种鉴定和功能鉴定，验证是否分离到目标功能物种。
11.所述的含底物的选择性培养基以及测定选择性培养基中底物的浓度变化情况，是根据筛选目的微生物的功能来决定，例如可以以na2s作为底物，来筛选硝酸盐还原硫氧化细菌，硝酸盐还原硫氧化细菌介导的硝酸盐还原硫氧化作用是偶联氮循环和硫循环重要的过程，可以通过离子色谱法测定硝酸盐、亚硝酸盐和硫酸盐的浓度以及生物量，通过底物的变化确定哪些样品具有功能活性。也可以以nas2o3作为底物，来筛选氧化硫代硫酸盐的细菌，硫代硫酸盐是硫化物氧化过程和硫循环的一个重要的中间产物。
12.所述的环境样品可以是土壤、沉积物、水样等环境样品
13.优选，所述的梯度稀释的目的是稀释至目标功能微生物绝迹的稀释度，根据样品中总的微生物量，最高梯度的稀释度可以调整。梯度稀释过程不限于10倍稀释，稀释倍数可以根据需要进行调整。
14.优选，所述的梯度稀释是：取1g或者1ml环境样品到10ml离心管中，加入9ml培养基，定义为10
‑1稀释度样品；再从10
‑1稀释度取1ml混合液到第二个10ml离心管，加入9ml培养基，定义为10
‑2稀释度样品；再从10
‑2稀释度取1ml混合液到第三个10ml离心管，加入9ml培养基，定义为10
‑3稀释度样品；重复类似的稀释操作，直到10
‑8稀释度。
15.优选，所述的选择降解底物最多的、活性的最高的若干个稀释度的样品组成绝迹稀释培养组是选择降解底物最多的、活性的最高的3
‑
4个稀释度的样品组成绝迹稀释培养组。
16.优选，还可以根据绝迹稀释培养组中各个物种分别达到的最高稀释度，分别根据公式“物种个数/每克干沉积物＝(最高位数稀释度稀释倍数/沉积物干重所占质量分数)”估算这些物种在原环境样品中的绝对丰度。
17.优选，所述的步骤e的物种鉴定是利用菌落pcr扩增16s rrna基因进行初步物种鉴定，验证是否分离到目标物种，对于群落中无法获得的一些物种，根据物种的种属信息设计培养基进一步尝试分离，同时也进一步划线纯化和鉴定。
18.本发明能从环境样品中，特别是如土壤和沉积物等复杂的环境体系中，经过梯度稀释，在不同的稀释度下富集出不同功能活性特征的微生物群落；对这些群落进行高通量测序，分析其结构组成。再根据这些群落的组成特征，针对出现的不同种类微生物，设计分离培养基进行分离筛选。该方法通过稀释培养过程，能有效分离环境中不同丰度和不同功能特征的物种，不仅能获得相对丰度高的物种，也能利于分离相对丰度较低的稀有物种。
附图说明：
19.图1是基于绝迹稀释培养组的功能微生物高通量分析和选育方法的操作流程图；
20.图2是实施例1中硝酸盐还原硫化物氧化功能微生物的绝迹稀释培养组样品的功能活性特征，a为硝酸盐的减少量，b为亚硝酸盐的产生量，c为硫酸盐的生产量；
21.图3是实施例1中硝酸盐还原硫化物氧化功能微生物的绝迹稀释培养组样品的群落结构物种组成气泡图，气泡大小表示物种的相对丰度，s01和s02为原始沉积物样品物种组成；
22.图4是实施例2中硝酸盐还原硫代硫酸盐氧化功能微生物的绝迹稀释培养组样品的功能活性特征，a为硝酸盐的减少量，b为硫代硫酸盐的减少量，c为硫酸盐的生产量；
23.图5是实施例2中硝酸盐还原硫代硫酸盐氧化功能微生物的绝迹稀释培养组样品的群落结构物种组成气泡图。气泡大小表示物种的相对丰度。
具体实施方式：
24.下面结合附图，进一步通过实施例对本发明的基于绝迹稀释培养组对特定功能微生物进行高通量分析和选育的方法进行详细说明。本发明不限于以下实施例的描述。
25.实施例1：绝迹稀释培养组对河道沉积物中硝酸盐还原硫化物氧化功能微生物进行多样性研究和高通量选育
26.硝酸盐还原硫氧化细菌介导的硝酸盐还原硫氧化作用是偶联氮循环和硫循环重要的过程，同时是河道沉积物中硝酸盐和硫化物等污染物去除的主要驱动力。硝酸盐还原硫氧化细菌参与硝酸盐还原硫化物氧化过程的阐明，需要首先对该类细菌在该环境中的多样性和功能特征有充分的认识。但是，对于河道沉积物中的硝酸盐还原硫氧化细菌的报道很少。而复杂的沉积物环境必然存在非常多种类的硝酸盐还原硫氧化细菌，这些分类和功能特征未知的硝酸盐还原硫氧化细菌将阻碍我们对沉积物中硝酸盐还原硫氧化过程的认识。
27.本实施例具体通过绝迹稀释培养组对河道沉积物中能利用硫化物做底物的硝酸盐还原硫氧化细菌进行高通量选育，以研究其在河道沉积物中的多样性和功能特征。
28.首先，对河道沉积物进行采样和理化参数的测定，通过测定包括硝酸盐、硫化物和二价铁等浓度的测定表征污染的情况。绝迹稀释培养组的操作流程如图1所示。分别取1g的河道沉积物通过ns
‑
m培养基进行梯度稀释，每升ns
‑
m培养基包含：2g kno3、1g na2s
·
9h2o、1g nahco3、2g k2hpo4、0.1g mgcl2和1000ml h2o，混合均匀，灭菌备用。梯度稀释的过程是取约1g河道沉积物到10ml离心管中，加入9ml ns
‑
m培养基，定义为10
‑1稀释度样品；再从10
‑1稀释度取1ml混合液到第二个10ml离心管，加入9ml培养基，定义为10
‑2稀释度样品；再从10
‑2稀释度取1ml混合液到第三个10ml离心管，加入9ml培养基，定义为10
‑3稀释度样品；重复操作，直到获得10
‑8稀释度。共做了12次重复的梯度稀释，获得的96个样品。将样品放置于厌氧培养箱30℃培养一个月后，测定底物变化情况。具体是用离子色谱法测定硝酸盐、亚硝酸盐和硫酸盐的浓度以及生物量，通过底物的变化确定哪些样品具有功能活性。经过梯度稀释，高稀释度样品的群落之间往往出现功能分化。因此，选取具有功能活性的最高几个稀释度进行群落组成分析。本实施例中，出现功能分化的最高三个稀释度为10
‑4、10
‑5和10
‑6。分别对这三个稀释度样品的菌体进行收集，按照一致的生物量(最终浓度总蛋白量约1μg/ml)分别接种到新鲜的ns
‑
m培养基的中，放置于厌氧培养箱中30℃培养，分别在0天、3天、5天、7天和9天取样测定底物变化情况，以底物的浓度变化率来表示硝酸盐还原和硫化物氧化活性的大小。此外，分别提取这些样品的总dna，通用高通量测序引物515f：5
’‑
gtg cca gcm gcc gcg gta a
‑3’
(如seq id no.1所示)和806r：5
’‑
gga cta chv ggg twt cta at
‑3’
(如
seq id no.2所示)，对16s rrna基因的v5区进行pcr扩增和高通量测序，通过生物信息学方法分析这些样品的群落组成特征及其与功能活性的关系。最后根据不同样品的物种组成，选择样品利用多种琼脂糖固体培养基对关键微生物菌种进行平板单菌落分离筛选。获得的单菌落通过16s rrna基因通用引物27f：5
’‑
aga gtt tga tcm tgg ctc ag
‑3’
(如seq id no.3所示)和1492r：5
’‑
tac ggy tac ctt gtt acg act t
‑3’
(如seq id no.4所示)进行pcr扩增和测序，序列拼接后在ncbi数据库进行比对，初步鉴定其分类特征。并测定菌株的硝酸盐还原硫氧化功能活性特征。
29.结果表明，10
‑4、10
‑5和10
‑6稀释度下的不同系列的样品具有不同的功能活性特征，其中10
‑4稀释度的所有样品都具有较高的功能活性，10
‑5比10
‑4较低，10
‑6稀释度下只有几个样品具有显著的功能活性，而且活性较低(如图2所示)。对这些样品的细菌群落组成进行分析，发现没有明显硝酸盐还原硫化物氧化活性的样品，菌群依然非常复杂，没有产生明显的富集效果(如图3所示)。而在具有明显功能活性的样品中，几个物种的不同组合构成了不同样品细菌群落的主体，这些物种包括thiobacillus(otu11)、ciceribacter(otu11)、rhizobium(otu3)和azonexus(otu5)与pseudoxanthomonas(otu8)等(如图3所示)，这些都可能是具有硝酸盐还原硫化物氧化功能的物种。而这些物种中，只有thiobacillus被报道过具有硝酸盐还原硫化物氧化功能。
30.进一步对绝迹稀释培养组中显著富集的物种进行分离纯化，成功获得了thiobacillus、ciceribacter、azonexus和pseudoxanthomonas四个优势属的代表菌株，以及shinella、ferrovibrio、hydrogenophaga和paracoccus等非优势属的代表菌株。对它们进行硝酸盐还原硫化物氧化功能验证，虽然只有paracoccus和thiobacillus的菌株能进行硝酸盐还原硫化物氧化，但是其他属，包括azonexus、ciceribacter、ferrovibrio、hydrogenophaga、pseudoxanthomonas和shinella的菌株都具有硝酸盐还原硫代硫酸盐氧化功能。虽然thiobacillus的菌株有直接的硫化物氧化能力，但是其他的菌株可能通过与thiobacillus进行代谢互作，间接的参与到硫化物的氧化过程中。此外，ciceribacter、ferrovibrio、hydrogenophaga和thiobacillus的菌株与16s rrna基因序列与ncbi数据库中近缘物种模式菌株的相似度都比较低，可能是新的物种分类。因此，绝迹稀释培养组不仅能获得具有直接功能的菌株，还能获得存在互作共生的微生物。
31.根据这些物种在不同稀释度出现的频数，thiobacillus和pseudoxanthomonas菌属在10
‑6稀释度下依然被富集起来，说明它们在原沉积物中的丰度最高。而ciceribacter、rhizobium和azonexus等在较低的10
‑4和10
‑5稀释度下富集，说明它们在原沉积物中的丰度较低。参考最大可能数法(mpn)，用公式：物种个数/每克干沉积物＝(最高位数稀释度稀释倍数/沉积物干重所占质量分数)对这些物种在原沉积物中的绝对丰度进行估算，thiobacillus和pseudoxanthomonas的绝对丰度约为107个/克干泥，rhizobium、azonexus和ciceribacter的绝对丰度约为106个/克干泥。这些结果与原沉积物中总细菌16s rrna基因高通量测序结果一致，thiobacillus的相对丰度约为1.27％，而rhizobium和azonexus为0.2
‑
0.3％，ciceribacter则低于0.1％。rhizobium、azonexus和ciceribacter等这些低丰度物种很容易由于高通量测序过程中的抽样效应而无法测到。因此，绝迹稀释培养组能将环境中低丰度的物种富集起来，而进一步的利于分离培养。
32.实施例2：绝迹稀释培养组对沉积物中硝酸盐还原硫代硫酸盐氧化功能微生物进
行多样性研究和高通量选育
33.承接实施例1，硫代硫酸盐是硫化物氧化过程和硫循环的一个重要的中间产物。许多硝酸盐还原硫氧化细菌都能氧化硫代硫酸盐，它们进行的硫代硫酸盐氧化作用是沉积物中硝酸盐还原硫氧化过程的重要组成部分。对具有硫代硫酸盐氧化功能的硝酸盐还原硫氧化细菌进行多样性研究也是探索沉积物硝酸盐还原硫氧化过程的基础。
34.本实施例的沉积物来源与实施例1一致，绝迹稀释培养组的操作流程也与实施例1一致。其中培养基(nts
‑
m培养基)的成分是：2g kno3、2g nas2o3·
5h2o、1g nahco3、2g k2hpo4和0.1g mgcl2和1000ml h2o。后续的培养、功能活性测定和细菌群落组成分析与实施例1一致。以硝酸盐消耗速率，硫代硫酸盐消耗速率和硫酸盐产生速率表征功能活性特点。
35.结果发现，10
‑7和10
‑8稀释度的样品出现大量的功能绝迹，而且10
‑5到10
‑8的四个稀释度的样品群落的功能活性差异分化较大，群落的功能活性沿着稀释度呈线性下降趋势(图4)。对这几个稀释度的样品进行细菌群落分析，与实施例1一致的是，明显具有硝酸盐还原硫代硫酸盐氧化活性的样品，微生物的群落组成发生了简化，说明产生了富集效果。由achromobacter、afipia、aminobacter、aridibacter、azonexus、bosea、comamonas、hydrogenophaga、micromonospora、pseudomonas、rhizobium、rhodanobacter、sphingopyxis、thermomonas和thiobacillus等属的不同组合构成了多种多样的群落(图5)。这些属很可能是具有硫代硫酸盐氧化功能的硝酸盐还原硫氧化细菌。其中achromobacter、bosea、pseudomonas、rhodanobacter和thiobacillus都已经被报道能进行硝酸盐还原硫代硫酸盐氧化作用。而没有活性的样品，微生物群落多样性依然很高，结构组成复杂。
36.通过固体培养基对这些物种进行分离纯化，成功获得了包括achromobacter、afipia、aminobacter、azonexus、bosea、comamonas、hydrogenophaga、micromonospora、pseudomonas、rhodanobacter、sphingopyxis、thermomonas和thiobacillus等优势属对应的菌株。此外，还获得了包括bacillus、bradyrhizobium、brevibacterium、castellaniella、chelatococcus、ciceribacter、ferrovibrio、mesorhizobium、ochrobactrum、paracoccus、pseudolabrys、pseudoxanthomonas、pusillimonas、ralstonia、shinella和stenotrophomonas等在群落中非优势属对应的菌株。对分离获得的菌株进行硝酸盐还原硫代硫酸盐氧化能力测定，发现所有的菌株都具有硝酸盐还原硫代硫酸盐氧化功能(如表1所示)。其中，包括afipia、aminobacter、bradyrhizobium、brevibacterium、bosea、castellaniella、chelatococcus、ferrovibrio、hydrogenophaga、mesorhizobium、ochrobactrum、pseudolabrys、pusillimonas和thiobacillus等菌属的菌株的16s rrna基因序列与ncbi数据库中的模式菌株的相似度都比较低，它们都可能是新的物种分类。因此，通过绝迹稀释培养组，我们获得了大量以往仍未报道过的具有硫代硫酸盐氧化功能的物种，其中发现了多个可能的新物种。目前已报道具有硝酸盐还原硫氧化功能物种约有50个菌属，而能进行硫代硫酸盐氧化作用的也仅有约30个属。通过本案例，我们把硝酸盐还原硫氧化细菌扩展到了约70个菌属，而进行硫代硫酸盐氧化功能的物种扩展到了约50个，大大加深了我们对这类细菌多样性的认识。
37.根据这些物种在不同稀释度出现的频数，用公式：物种个数/每克干沉积物＝(最高位数稀释度稀释倍数/沉积物干重所占质量分数)对它们在原沉积物中的绝对丰度进行
估算。如表1所示，估算的结果多数与高通量测序的相对丰度结果一致，相对丰度低的或者无法检测的，估算得到的绝对丰度也较低。表明多数的这些物种在沉积物中的丰度很低，高通量测序技术带有的随机性和抽样效应往往会难以测到这些物种的存在。然而，通过绝迹稀释培养组，我们成功将这些物种富集起来，并通过分离纯化获得了纯培养。因此，通过绝迹稀释培养组，我们发现了沉积物中这些物种的多样性分布规律，大多数的物种丰度很低，它们可能在合适的条件下进行繁殖，发挥它们的生态功能。
38.表1通过绝迹稀释培养组分离获得的菌株及其功能和在原沉积物中的相对丰度
39.[0040][0041]
(+)表示阳性，(
‑
)表示阴性，nd表示未检测到。