¹³C-DNA-SIP Method for Tracing Rhizosphere and Hyphosphere Phosphate Solubilizing Bacteria

材料与试剂

1. 各种型号的移液枪和枪头
2. 2 ml离心管
3. PVC板
4. 30 μm孔径尼龙膜
5. 亚克力板
6. 凡士林
7. 水管
8. 温度计
   
9. 2 mm，80 μm和30 μm孔径筛子
10. 50 ml注射器
11. 橡皮塞（∅=1 cm，海佳硅橡胶制品有限公司，浙江，中国）
12. 直径1 mm玻璃珠
13. 中速定量滤纸
14. 1 L烧杯
15. 玻璃棒
16. 8.9 ml离心管 (Beckman)
17. 10 ml注射器及其针头
18. 土壤DNA提取试剂盒 (MP，美国)
19. 氯化钙
20. ¹³CO₂和¹²CO₂ (atm 99 %，明尼克科技有限公司，北京，中国)
21. 30 % H₂O₂
22. 三氟乙酸铯 (CsTFA)（GE, Healthcare Limited Compony, The UK）
23. KCl
24. EDTA-Na₂·2H₂O
25. Tris base
26. 37 % HCl
27. q-PCR SYBR Mix (Biomad，北京，中国)
28. 异丙醇
29. 无水乙醇
30. ddH₂O
    

仪器设备

1. 半导体制冷片
2. 风扇
3. 水泵
4. 散热器
5. CO₂同位素检测仪 (Los Gatos Research，美国)
6. Deltaplusxp质谱仪 (Thermo Scientific，Bremen，德国)
7. 元素分析仪 (Flashea 1112，CE Instruments，Wigan，英国)
8. Nano drop核酸蛋白检测仪
9. 折光仪 (Reichert AR2000 Digital Refractometer)
10. 蠕动泵 (master flex C/L)
11. -80摄氏度冰箱
12. 荧光定量PCR仪
    

软件和数据库

1. Greengene数据库 (http://greengenes.secondgenome.com/)

实验步骤

一、试验设计

1. 不接种AM真菌的对照NM；
2. 两侧根室分别接种Rizophagus intraradices（R.i）和Gigaspora margarita （G.m）。后续试验中会进行¹³CO₂脉冲标记，用于活性微生物组群落结构和功能的测定。分别进行¹³CO₂标记和¹²CO₂标记 (对照处理)，每个处理设3次重复，完全随机排列。
   

二、供试装置
        为了达到在同一植物根系上接种两种不同AM真菌的目的，试验采用分根培养的根盒。另外，为了能够得到不受根系影响的菌丝，我们将根室和菌丝室之间用30 μm的尼龙网分隔，使根系限制在根室中生长，而AM真菌的菌丝可以穿过该尼龙网，在菌丝室中生长，从菌丝室的土壤中吸收养分。试验装置用PVC板制作，包括根室 (中间两个) 和菌丝室 (两侧) 两部分 (图1)。收获前一周，对棉花植株进行¹³CO₂同位素脉冲标记，为了防止¹³CO₂本身对土壤¹³C丰度和土壤微生物的影响，使用如图1所示的装置将地上部和土壤分离，脉冲标记室中只有地上部。另外为防止温度过高影响植物光合作用，需要添加降温装置，降温装置由半导体制冷片、风扇、水泵及散热器组成，使得标记过程中标记室内的温度不超过35 °C。


 图1. 分根培养系统和二氧化碳同位素脉冲标记装置
注：培养体系由4个分室组成，中间有两个根室，另外两个为菌丝室。每个分室为5 cm × 5 cm × 12 cm，由PVC板组合而成。使用30 μm尼龙网将根室和菌丝室隔开，防止根系生长进入菌丝室，但允许菌丝生长到菌丝室中。脉冲标记设备包括支架、标记室和冷却系统。支架高度为15 cm，由透明亚克力板制成。支架中间有一个5厘米直径的孔，上端有一个5厘米高的水槽。罩子也用透明亚克力板制成，并用水槽中的水密封。冷却系统由半导体、风机、水泵、散热器和两根水管组成。利用散热器用水来交换热量以便于控制温度。

三、供试土壤
        土壤风干并过筛 (2 mm)。此外，土壤的一部分过30 μm筛，用于后续菌丝的收集。具体方法为：将土壤放到足量的水中 (土：水 = 1:5) 充分搅拌后，过30 μm筛子，保留过筛的水土混合物，丢弃掉不能过筛的部分，沉淀后小心倒掉上层水分，剩下的土壤于烘箱中70 °C烘干，用粉碎机粉碎待用。另外，将直径1 mm的玻璃珠用1%盐酸酸化，然后用去离子水清洗，直到其pH值变为中性。土壤和玻璃珠在使用前经过γ射线辐照以杀死土壤中的土著微生物群落 (25 kGy，60 Co γ-射线，中科院原子能研究所，房山，北京)。辐照后的材料在使用前存放7天，使土壤中各元素恢复到稳定状态。

四、土壤原始细菌群落
        为保持土壤原始细菌群落的一致性，在移栽过程中要向根室和菌丝室灭菌土壤中分别加入5 ml土壤滤液。土壤滤液是通过将30 g未经灭菌的新鲜土壤放入300 ml无菌水中，搅拌均匀，六层定量滤纸 (Whatman中速滤纸) 过滤后获得，这样既能允许普通土壤微生物通过，又能有效过滤掉AM真菌的孢子和菌丝 (Wang等， 2016)。另外，为了消除菌种本身所携带微生物对根室土壤原始细菌群落的影响，同样需要制备菌种滤液，两种供试菌剂R.i和G.m各30 g，分别放入300 ml无菌水中搅拌，六层中速滤纸过滤，得到菌种滤液。移栽过程中，在接种R.i的根室加入G.m菌种滤液5 ml，在接种G.m的根室加入R.i菌种滤液5 ml，不接种真菌的对照处理加入两种AM真菌菌种滤液各5 ml。接种R.i的根室需要加入5 ml土壤滤液，另外还需要加入5 ml G.m的菌种滤液。接种G.m的根室需要加入5 ml土壤滤液和5 ml R.i的菌种滤液。对于不接种的对照处理，两个根室需要加入5 ml土壤滤液以及两种菌种滤液各5 ml。

五、种子萌发和育苗
        供试植物选用棉花 (Gossypium herbaceum L.)，品种为“新陆早32”。挑选重量一致的棉花种子，将种子用10% H₂O₂表面消毒10 min，用去离子水彻底清洗8次以洗净种子表面H₂O₂，然后将种子放在潮湿的滤纸上 (40×20 cm)，26 °C黑暗中发芽2天。然后将种子放在新的潮湿滤纸顶部5 cm处，并将胚根朝下。将滤纸卷成圆柱形，并置于水培系统中 (12 h光照/12 h黑暗，26 °C) 培养17天，以使根系伸长。

六、移栽及培养过程
        育苗结束后，挑选长势一致，且有7条根系的棉花幼苗用于后续的试验。移栽时首先在两侧根室分别加入580 g过2 mm筛的供试土壤，剪掉棉花幼苗的主根后将侧根平均分配到两侧根室，加入20 g相应的AM真菌接种剂，约600个孢子 (对照处理接种20 g灭菌的菌剂)，5 ml土壤滤液，5 ml相应的菌种滤液，然后将水补齐到108 ml使土壤含水量达到18% (约为田间持水量的70%)，再将剩余的100 g土壤铺到最上层。需要注意的是在接种一种AM真菌时，需要严格密封另一侧根室和菌丝室，防止交叉污染。在两侧菌丝室分别加入600 g土壤与玻璃珠的混合物 (300 g过30 μm筛子的土壤和300 g直径1 mm的玻璃珠) 以及5 ml土壤滤液和103 ml的无菌去离子水。试验在中国农业大学资源与环境学院温室中进行。每天上午浇水至含水量的18 %，每周随机变换根盒位置，7周后进行¹³CO₂标记。

七、 ¹³CO₂标记
        播种7周后，在如图1密闭的标记室中进行¹³CO₂ (atm> 99%) 脉冲标记处理并以¹²CO₂ (atm> 99%) 脉冲标记作为对照。标记室 (体积约127000 cm3) 包含一个透明亚克力罩和一个底部托盘。底部托盘的外缘设置一个水槽用于标记室密封，植物地上部穿过托盘中央的孔进入标记室，并用凡士林密封，以防止标记期间标记室内气体泄漏或环境中CO₂的影响。在脉冲标记过程中，植物密封在标记室内，并在其中放置一个风扇，以循环空气。采用半导体散热片、风机、水泵、水管和散热器组合的冷却设备，将标记室温度控制在35 °C以下 (在35 °C时棉花的光合作用处于较高水平)。标记时间从9:00到17:00，这是白天光合作用最强的时间段，利于标记气体的同化。每小时使用注射器向标记室内注入100 ml气体，并通过二氧化碳同位素检测仪 (Los Gatos Research，美国) 实时测定标记室内的CO₂浓度，每天16:00 最后一次注入二氧化碳，待标记室内的二氧化碳消耗完后，打开标记室，第二天标记时再重新密封标记室，并检查标记装置是否密封完好。持续标记7 d。为了避免二氧化碳标记期间植物蒸发产生的蒸汽在标记室内壁形成水珠影响标记装置的透光性，将三袋氯化钙 (每袋100 g) 放在标记室内。每天在晚上标记结束后，取出氯化钙在105 °C的烘箱中干燥2 h以重复使用。

八、样品收获
        标记结束后的第二天，拆除标记装置，破坏性取样，地上部分用去离子水洗净，吸水纸擦干，地下部分首先用抖根法收集菌根际土壤，一部分存储于室温用于δ¹³C的测定，另一部分于-20°C存储用于测定细菌16S rDNA和解磷功能基因的群落多样性。菌丝际基质一部分保存于室温条件下，用于δ¹³C的测定。另一部分储存于-20°C，用于菌丝的收集和DNA的提取。具体方法为室温融化后，菌丝室土壤和玻璃珠的混合物用无菌水（121°C，30 min灭菌）经30 μm筛网冲洗，过滤掉土壤，菌丝会缠绕在玻璃珠上，然后将菌丝体和玻璃珠混合物全部转移到烧杯中，用玻璃棒搅拌，将漂浮在水面的菌丝挑出并收集到无菌的离心管中待用。通过这种方法，只有紧密地附着在菌丝上的细菌才能留存下来。

九、δ¹³C的测定
        菌根际土壤和菌丝际土壤¹³C丰度用δ¹³C‰表征。土壤样品风干，研磨，并过80 μm筛。在中国农业大学资源与环境学院稳定同位素实验室，使用Deltaplusxp质谱仪 (Thermo Scientific，Bremen，德国) 和元素分析仪 (Flashea 1112，CE Instruments，Wigan，英国) 以连续流动模式测定土壤样品的碳同位素丰度。元素分析仪的燃烧温度为1020 °C。δ13C‰以V-PDB (维也纳-佩迪-贝勒姆尼) (Deniro等，1978) 为标准计算，方程式如下：
δ¹³C ‰=(R_样品/R_标准品-1)×1000
式中，δ¹³C为样品碳同位素比值，单位为千分之一 (‰)，R样品和R标准品分别为样品和标准品的¹³C/¹²C比值。δ¹³C测量值的SD应该小于0.15‰。

十、DNA提取、密度梯度离心及q-PCR分析
        根据生产商的说明书，使用FastDNA Spin Kit (MP Biomedicals LLC，Santa Ana，CA，美国) 试剂盒提取¹³CO₂标记和¹²CO₂标记的菌根际土壤和菌丝样品的DNA，并通过三氟乙酸铯 (CsTFA) 密度梯度离心分离不同密度梯度的DNA，离心梯度从下到上分成16层。具体方法为：

1. 将提取纯化后的DNA用nano-drop测定其浓度，保证用Beckman超速离心机离心时的DNA质量 (3-5 µg)；
2. 在已灭菌的10 ml离心管中配置混合液：5 ml三氟乙酸铯 (CsTFA)，DNA样品和GB缓冲液的体积为3.9 ml，将混合液混匀；
3. 先用折光仪 (Reichert AR2000 Digital Refractometer) 测定灭菌超纯水的折光率 (nD-TC模式)，加入约150 µl水，经验数值为1.3333 ± 0.0002，将水吸干，测定混合液的折光率，加入混合液约75 µl，若超过1.3610 ± 0.0002的范围，则用CsTFA或GB调整；
4. 将调整好的混合液加入到8.9 ml的离心管 (Beckman) 中配平 (无误差)；
5. 45400 × g，转子型号为Ti90，离心36-40 h后关闭离心机，待转速降为0后打开离心机盖，去除转子轻放于水平实验台上；
6. 将蠕动泵 (master flex C/L) 的流量调节好，流速调节为550 μl/30 s，泵的一端连接灭菌超纯水，另一端用针头接到离心管上方接近管口处；在离心管底部扎孔，打开泵，开始计时，底部离心液用2 ml离心管接住，30 s换一个离心管，8.9 ml的离心液大约分成16个片段；
7. 用折光仪先测定ddH2O折光率，再从低密度到高密度测定各个片段的折光率，以估算浮力密度，浮力密度计算公式为：y = 384.44060*x*x-1031.00836*x + 692.65494。其中y为浮力密度，x为折光率 (nD-TC)；
8. 在各个片段中加入500 µl异丙醇，混匀配平，-80 °C放置1 h后以14,000 × g，4 °C离心1 h (放置离心管的方向便于抽取溶液)；将上清液转移至新的2 ml Tube中，-20 °C保存 (此步骤是为了防止DNA未沉淀好，可留作再次离心使用)；在有DNA沉淀的2 ml离心管中加入150 µl 70%乙醇，用移液枪吹打冲壁洗涤，14,000 × g，4 °C离心30 min；弃上清液，吹干乙醇；加入30 µl ddH2O溶解DNA，分装后在-20 °C保存。
   

十一、定量PCR检测
        将制备好的标准物与样品DNA同时进行定量PCR检测。具体步骤为：首先配制足量的混合液，每个样品混合液由12.5 μl Biomed q-PCR SYBR Mix，1.0 μl BSA，1.0 μl引物 (Ba519f/Ba907r)，7.5 μl ddH2O组成。引物序列为：Ba519f (5’-CAGCMGCCGCGGTAANWC-3’)-Ba907r (5’-CCGTCAATTCMTTTRAGTT -3’) (Stubner, 2002)。每个样品及标准物均准备3个技术平行。另外为了保证混合液足够，多准备5个样品的混合液。然后将混合液离心混匀后，向定量PCR板的每个孔中准确加入23 μl 的混合液，然后再加入2 μl的标准物或样品DNA。封膜后放入定量PCR仪q-Tower (Jena，德国) 中检测，反应条件为94 °C下预变性5 min；94 °C变性30 s，60 °C退火45 s，72 °C延伸1 min，从第二步开始循环40次，每次循环后检测荧光，收集溶解曲线。以标准物拷贝数的Log值及Ct值绘制标准曲线，并根据标曲计算每个样品中的DNA拷贝数 (Tellmann, 2008)。密度梯度离心后DNA拷贝数在不同分层的分布如图2所示。


图2. 密度梯度离心后R.i和G.m菌丝际DNA在不同浮力密度分层中的分布

十二、高通量测序
        根据密度梯度离心后每一层q-PCR结果，每个处理挑选浮力密度约为1.58的那一层DNA样品进行高通量测序，将DNA样品送往生物技术公司，利用Illumina MiSeq PE-300平台，进行16S rDNA V3-V4区测序和解磷功能基因 (以编码碱性磷酸酶alp的基因phoD为例) 的测序，所用引物为338f (5’-ACTCCTACGGGAGGCAGCA-3’) 和806r (50’-GGACTACHVGGTWTCTAAT-3’) (Xu et al., 2016); phoD-F733 (5’-TGGGAYGATCAYGARGT-3’)和phoD-R1083 (5’-CTGSGCSAKSACRTTCCA -3’) (Sato et al., 2005)。得到测序原始数据后，用QIIME (v1.8.0) 对所得数据进行分析。简而言之，根据Barcode序列区分样品序号，剔除掉低质量的序列，用FLASH软件合并两端序列。去除嵌合体后，按照 97 %的标准划分可操作分类单元 (OTUs)，使用UCLUST默认参数从每个OTU中选择一个代表性序列 (Edgar, 2010)。根据Greengene数据库对所选取的代表性序列进行BLAST分类，匹配度最高的即为该OTU所属的分类群 (DeSantis et al., 2006)。

十三、解磷功能基因片段拷贝数分析
        根据密度梯度离心后每一层q-PCR结果，每个处理挑选浮力密度约为1.58的那一层DNA样品进行功能基因拷贝数定量PCR检测，以编码碱性磷酸酶alp的phoD基因为例：定量PCR体系和程序同步骤11，引物为phoD-F733 (5’-TGGGAYGATCAYGARGT-3’)- phoD-R1083 (5’-CTGSGCSAKSACRTTCCA-3’) (Sato et al., 2005)。

溶液配方

1. GB缓冲液配方：
           配1 L GB所需：KCl: 7.455 g，EDTA-Na₂·2H₂O: 0.372 g，Tris base: 12.114 g，37% HCl: 2.43 ml
   

致谢

        该方法是在王菲 (河南科技学院，副教授)和石宁(山东省农业科学院，助理研究员) 两位老师的博士论文基础上改进的，在此一并感谢。同时，本工作受到了国家自然科学基金联合基金重点支持项目“灰漠土磷养分高效利用相关的生物肥力性状演变过程及培育机制” (U1703232) 的资助。

参考文献

1. Wang, F., Shi, N., Jiang, R. F., Zhang, F. S., Feng, G. (2016). In situ stable isotope probing of phosphate-solubilizing bacteria in the hyphosphere. J Exp Bot. 67(6): 1689-1701.
2. Deniro, M. J., Epstein, S. (1978). Carbon isotopic evidence for different feeding patterns in two hyrax species occupying the same habitat. Science. 201: 906-908.
3. Stubner, S. (2002). Enumeration of 16S rDNA of Desulfotomaculum lineage 1 in rice field soil by real-time PCR with SybrGreen™ detection. J Microbiol Methods. 50: 155-164.
4. Tellmann, G. (2008). The E-Method: a highly accurate technique for gene-expression analysis. Nat Methods. 3:1–2.
5. Xu, N., Tan, G., Wang, H., Gai, X. P. (2016). Effect of biochar additions to soil on nitrogen leaching, microbial biomass and bacterial community structure. Eur J Soil Bio. 74: 1-8.
6. Sato, K., Suyama, Y., Saito, M., Sugawara, K. (2005). A new primer for discrimination of arbuscular mycorrhizal fungi with polymerase chain reaction-denature gradient gel electrophoresis. Grassl Sci. 51: 179-181.
7. Edgar, R.C. (2010). Search and clustering orders of magnitude faster than BLAST. Bioinformatic. 26: 2460-2461.
8. DeSantis, T. Z., Hugenholtz, P., Larsen, N., Rojas, M., Brodie, E. L., Keller, K., Huber, T., Dalevi, D., Hu, P., Andersen, G.L. (2006). Greengenes, a chimera-checked 16S rRNA gene database and workbench compatible with ARB. Appl Environ Microbiol. 72:5069-5072.