水稻根系微生物组研究中的样本种植、取样和16S rRNA 基因扩增子文库制备方法
Planting, Sample Collection and Library Preparation of 16S rRNA Gene Amplicon Sequencing for Rice Root Microbiomes    

材料与试剂

1.  Parafilm封口膜 (Bemis, catalog number: PM-996)
2.  50 ml离心管 (BD Falcon, catalog number: 352070)
3.  滤纸 (SEP, catalog number: DXLZ11F)
4.  微量离心管 (1.5, 2, 5 ml; Eppendorf, catalog number: 0030125150, 30123344, 30119401)
5.  移液器枪头 (nuclease-free, 10, 200, 1000 µl; Axygen, catalog number: T-300, T-200-Y, T-1000-B)
6.  96孔PCR板 (Jet Keen Biotechnology, catalog number: PC-0200-9B)
7.  96孔PCR板封板膜 (Axygen, catalog number: PCR-TS)
8.  次氯酸钠 (Macklin, catalog number: S828471)
9.  乙醇 (SIGMA, catalog number: E7023-500ML)
10.  Murashige Skoog培养基 (含维生素) (MS, Caisson, catalog number: MSP09)
11.  蔗糖 (HuShi, catalog number: JC-SJ03064)
12.  氯化钠 (NaCl, SIGMA, catalog number: S7653-250G)
13.  磷酸氢二钠 (Na₂HPO₄·2H₂O, SIGMA, catalog number: 04272-1KG)
14.  磷酸二氢钠 (NaH₂PO₄·H₂O, SIGMA, catalog number: 71504-250G)
15.  琼脂粉 (SIGMA, catalog number: A7921-100G)
16.  土壤DNA提取试剂盒 (FastDNA^TM SPIN Kit for Soil, Mpbio, catalog number: 116560200)
17.  无核酸水 (Nuclease-Free Water, Qiagen, catalog number: 129115)
18.  1x TE缓冲液 (Coolaber, catalog number: SL2081)
19.  PicoGreen DNA定量试剂盒 (Quant-iT™ PicoGreen™ dsDNA Assay Kit, Invitorgen, catalog number: P11496)
20.  PrimeSTAR热稳定DNA聚合酶体系 (PrimeSTAR HS DNA Polymerase, TaKaRa, catalog number: R010A)
21.  琼脂糖 (Biowest, catalog number: G-10)
22.  TAE Buffer (Thermo, catalog number: 15558026)
23.  6x Loading Buffer (TaKaRa, catalog number: 9156)
24.  Beckmen Coulter Agencourt AMPure XP核酸纯化试剂盒 (Beckman, catalog number: A63880)
25.  琼脂糖电泳凝胶纯化回收试剂盒 (Promega Wizard^® SV Gel and PCR Clean-Up System, Promega, catalog number: A9282)
26.  核酸电泳染料 (Solarbio, catalog number: G8140)
    

仪器设备

1.  超净工作台/生物安全柜 (Thermo Fisher Scientific, catalog number: 51025411)
2.  镊子 (Jinzhong, catalog number: JD5020)
3.  三角瓶 (aladdin, catalog number: T4227-500ml-1EA)
4.  130 mm方形培养皿 (MaiSiNuo, catalog number: HZX064-1)
5.  1 L试剂瓶 (Cleman, catalog number: CN-600-1000)
6.  电子天平 (Mettler Toledo, catalog number: AL104)
7.  高压蒸汽灭菌锅 (PHCbi, catalog number: MLS-3781-PC)
8.  剪刀 (Jinzhong, catalog number: J21130)
9.  高速离心机 (Eppendorf, model no. 5810R)
10.  平板摇床 (Kylin-Bell Lab Instruments, model no. TS-2)
11.  单道移液器 (10, 20, 100, 200, 1000 µl; Eppendorf, catalog number: 3120000020, 3120000038, 3120000046, 3120000054, 3120000062)
12.  12通道移液器 (10, 100, 300 µl; Eppendorf, catalog number: 3122000027, 3122000043, 3122000060)
13. 旋涡混合器 (ZangHan, catalog number: VM100)
14.  金属水浴锅 (Dry Baths, Thermo, catalog number: 88870011)
15.  破碎仪 (Mpbio FastPrep-24TM, Mpbio, catalog number: MPFastPrep-24 5G)
16.  高速微型离心机 (Thermo Fisher Scientific, model no. Heraeus Pico 17)
17.  紫外分光光度计 (Thermo Fisher Scientific, model no. NanoDrop ND-3300)
18.  96孔酶标板 (Costar, catalog number: 3590)
19. 酶标仪 (Molecular Devices^®, model no: PARADIGM)
20.  PCR仪 (T100TM; BIO-RAD, model no. 1861096)
21.  电泳仪 (JUNYI-DONGFANG, catalog number: JY300HC)
22.  凝胶成像仪 (BIO-RAD, model no. Universal Hood II)
23.  切胶刀片 (Jinzhong, catalog number: J11010)
24.  磁力架 (Invitrongen DynaMag™-96 Side Magnet, Invitrongen DynaMag™-2 Magne, Invitrongen, catalog number: 12027, 12321D)
25.  -80 °C超低温冰箱 (Thermo Fisher Scientific, model no. 907)
26.  -20 °C冰箱 (Haier, catalog number: DW-40L92)
    

实验步骤

一、种子消毒、萌发及移栽
注：以下实验需在超净台内完成。

1. 挑选饱满的水稻种子。去除水稻种子颖壳，弃去干瘪的种子，不要损伤胚，置于无菌三角瓶内。
2.  使用酒精灭菌。加入70%酒精，液面没过种子，消毒30 s，期间不停晃动瓶身，保证每粒种子与酒精充分接触。弃去酒精。
3.  使用次氯酸钠灭菌。加入2.5%有效氯浓度的次氯酸钠溶液，液面没过种子，消毒15 min，期间不停晃动瓶身。弃去次氯酸钠。此步重复3次。
4.  使用无菌去离子水清洗。加入无菌去离子水，液面没过种子，清洗10 min，期间不停晃动瓶身。弃去无菌水。此步重复3次。
5.  接种至培养基。使用无菌镊子将种子整齐均匀的平铺于1/2 MS固体培养基表面，每皿可容纳10粒种子，胚向上放置于平板底部1/3处 (图 1A)。使用Parafilm封口膜封口。
6.  培养水稻种子。将培养皿竖直放置 (图 1B)，于25 °C，16 h光照，21%湿度的培养间内培养5-7 d。
   
   
   图1. 水稻种子培养方法. A. 灭菌种子铺板位置； B. 竖直培养方法。
   
   
7.  为防止培养基中残留的营养在移栽后对植物生长发育产生影响，需去除根系残留培养基。挑选长势一致的植物无菌苗，使用无菌镊子将无菌苗从培养基中拔出，用无菌去离子水洗去根上残留的培养基。
8.  将无菌苗移栽至大田培养。将清洗干净的幼苗以6 × 6的方阵移栽至田间，每株水稻间隔20 cm。不同品种间隔至少30 cm。在正常的水肥条件下生长8周。
   
   

二、样本收取

1.  土体土壤 (Bulk soil) 样本的收取 (图 2)。于田间未种植植物的地块上使用铲子挖取10 cm × 10 cm × 10 cm的土块。于切面均匀收取各深度土壤约10 g，此样本即土体土壤样本。
2.  根际土样本的收取 (图 2)。于水稻方阵中央4 × 4的区域中选取长势一致的水稻植株，使用铲子在其四周挖出15 cm × 15 cm × 10 cm的土块，将水稻拔离土壤。佩戴无菌手套将水稻从茎秆基部均匀的劈分成两半。使用无菌镊子和剪刀在断层中挑取2至3个分蘖的全根剪下，于30 ml无菌去离子水中大力摇晃至根系表面附着的土壤颗粒脱落，使用无菌镊子挑出根系弃去。剩余液体于3,000 rpm条件下离心15 min。弃去部分上清至保留约5 ml土壤沉淀和上清液，混合均匀。此样本即根际土样本。
3.  根系样本的收取 (图 2)。为了获得和水稻根系紧密联系的微生物数据，需要获得尽可能干净的根系样本。首先佩戴无菌手套将植株根系上的大块土壤去除，保留3-4个分蘖的全根；之后使用无菌去离子水冲洗根系至无可见的土壤颗粒；最后置于含有35 ml 1x PBS的50 ml管中，使用平板摇床于180 rpm条件下清洗15 min，更换新管重复清洗3次。使用无菌滤纸吸干根系残留水分。使用无菌剪刀将根系破碎成2 mm长的小段，此样本即为根系样本。
   注：所有样本应置于无菌的离心管内于-80 °C环境保存。
   
   
   图2. 根、根际土、土体土壤样本的收取方法
   
   

三、PCR模板DNA制备
注：以下实验应在超净台内完成。

1.  确定不同类型样本的上样量。将样本转移至土壤DNA提取试剂盒的Lysing Matrix E管内。样本量为根系样本0.2 g，根际土样本450 µl，土体土壤样本0.3 g。将样本和Lysing Matrix E内研磨珠混合均匀。
2.  研磨样本。将Lysing Matrix E置于破碎仪内，7,200 rpm破碎30 s，间隔30 s，再次破碎30 s，重复2次。根系样本因韧性较大，可在2次重复间将Lysing Matrix E置于液氮中冷冻5 s使研磨更加充分。
3.  溶解细胞膜蛋白。向Lysing Matrix E中添加978 µl的Sodium Phosphate Buffer以及122 µl的MT Buffer。置于破碎仪内，7,200 rpm破碎30 s，间隔30 s，再次破碎30 s，重复2次，使细胞充分裂解。
4.  去除不可溶杂质。将Lysing Matrix E置于离心机中，14,000 rcf离心15 min。将上清约900 µl转移至新2 ml离心管。
5.  沉淀样本中蛋白。向样本中添加250 µl PPS，轻柔的上下翻动20次，室温孵化10 min，14,000 rcf离心5 min。将上清转移至新2 ml离心管内。
6.  使用Binding Matrix吸附样品中DNA。向2 ml管内加入900 µl的Binding Matrix。于旋转混合器上慢速混合3至4 min。静置8 min。去除上清液1100 µl。剩余样品混合均匀后转移至SPIN Modules内。14,000 rcf离心1 min，去除滤液。
7.  再次溶解并去除样品中蛋白质。向SPIN Modules内添加500 µl Concentrated SEWS-M (需按试剂盒要求添加乙醇后使用)。使用移液器吹吸混匀过滤柱内液体。14,000 rcf离心1 min，去除滤液。为保证滤液去除彻底，再次空管14,000 rcf离心2 min，将过滤柱转移至新的2 ml离心管内。
8.  从Binding Matrix中溶解并分离DNA。向过滤柱内添加100 µl DES。置于金属水浴锅内，55 °C孵化5 min。14,000 rcf离心1 min。弃去过滤柱。置于-80 °C冻存。
9.  粗测DNA浓度以确定稀释倍数。因后续流程中的DNA精确定量对低浓度DNA敏感，需提前稀释样本DNA，定量DNA上样的适宜浓度为5-20 ng/µl。使用分光光度计测量DNA浓度以确定各样本合适的稀释倍数。
10.  稀释DNA。将2 µl样本转移至96孔PCR板，按照计算好的稀释倍数添加Nuclease-Free Water稀释并混合均匀。
11.  制备定量DNA的标准曲线。使用1x TE稀释PicoGreen DNA定量试剂盒中的λ DNA 50倍至2 μg/ml。按表 1所示将稀释后的λ DNA和1x TE加入96孔酶标板中，混合均匀，每个浓度设置两个重复。
    
    表 1. 定量DNA的标准曲线的浓度及配置方法
    

    λ DNA浓度 (ng/µl)	1x TE体积 (µl)	2 μg/ml λ DNA体积 (µl)
    0.75	12.5	37.5
    0.5	25	25
    0.4	30	20
    0.25	37.5	12.5
    0.1	45	5
    0.05	47.5	2.5
    0.01	49.5	0.5
    0	50	0

    
    
12.  使用1x TE稀释样本DNA。将2 µl稀释好的样本DNA和48 µl 1x TE加入96孔酶标板的空孔中，混合均匀。
13.  添加荧光染料PicoGreen。于避光的环境中，将试剂盒中的PicoGreen荧光染料使用1x TE稀释200倍，并添加50 µl至标准曲线样品孔及待测样品孔内，混合均匀，孵育5 min。
14.  测量样品的荧光。使用酶标仪于480 nm的激发光和520 nm的吸收光下读取荧光值，每个样品读取3次，读数取平均值。
15.  计算样本DNA的浓度。使用标准曲线的0 ng/µl的2个样品孔的读数平均值作为空白对照值校准标准曲线样品孔和待测样品孔的读数。根据标准曲线样品的DNA浓度和荧光值绘制标准曲线。将样品孔的荧光值带入标准曲线以求得DNA浓度。
16.  稀释DNA模板至3.5 ng/µl。根据荧光定量的DNA浓度计算稀释量，将DNA样本于96孔PCR板内使用无核酸水稀释至浓度为3.5 ng/µl。
    
    

四、两步法扩增样品16S rRNA基因
注：以下实验应在超净台内完成。

1.  引物设计及稀释。两步法扩增16S rRNA基因使用两对引物逐步添加用以区分同一文库中不同样本的标签 (Barcode) 序列、用以区分不同文库的索引 (index) 序列及Illumina测序所需序列，引物结构如图 3所示，引物序列见附件 1，引物合成于英潍捷公司，引物纯化方式为HPLC纯化。使用前需使用无核酸水将引物稀释至100 nmol/ml。
   
   
   图 3. 16S rRNA基因两步扩增所用引物的结构
   
   
2.  第一轮PCR扩增。第一轮PCR扩增通过引物在16S rDNA的799序列位置前引入了不同的Barcode序列用以区分同一文库中的不同样本，每个文库可承载60个不同的样本对应60个不同的Barcode序列。使用PrimeSTAR热稳定DNA聚合酶体系扩增，扩增体系如表 2所示。所有样本平行扩增3份，同时每个样本设立一个使用无核酸水替代模板DNA的阴性对照，用于指示使用的实验试剂及实验过程是否存在污染。扩增条件如表 3所示。扩增产物置于-20 °C环境保存。
   
   表2. 第一轮PCR扩增体系
   

   试剂名称	添加体积 (µl)
   Nuclease-Free Water	16.8
   5x PrimeSTAR Buffer (Mg2+ Plus)	6
   dNTP Mixture	2.4
   正反向引物	0.75
   3.5 ng/µl的DNA模板	3
   PrimeSTAR HS DNA Polymerase	0.3

   
   表3. 第一轮PCR扩增条件
   

   步骤	温度 (°C)	时间
   1	98	30 s
   2	98	10 s
   3	55	15 s
   4	72	1 min
   重复步骤2至4共25个循环
   5	72	5 min
   6	12	保温

   
   
3.  通过琼脂糖电泳进行第一轮PCR的质量控制。将平行扩增的样本混合均匀后吸取5 µl与1 µl 6x Loading Buffer混合加入1%琼脂糖凝胶中进行电泳。若阴性对照无条带 (错误范例见图 4B)，各样品间条带亮度一致 (错误范例见图 4C)，位置位于400 bp左右，则扩增质量合格 (图 4A)。
   
   
   图 4. 利用琼脂糖电泳对PCR过程进行质量控制的范例
   图中单数泳道为样品PCR扩增产物，双数泳道为前一泳道样本的阴性对照； A. 质量合格的PCR扩增产物的电泳条带； B. 阴性对照出现扩增条带的错误范例； C. PCR产物DNA浓度低导致条带亮度不一致的错误范例。
   
   
4.  使用磁珠结合PCR产物中的DNA。将Beckmen Coulter Agencourt AMPure XP核酸纯化试剂盒中的磁珠放置于室温环境孵育30 min。取50 µl PCR产物和50 µl磁珠于96孔PCR板内混匀，室温放置10 min。
5. 去除PCR产物中的杂质以纯化DNA。将混合了磁珠的PCR产物置于磁力架上，放置5 min。使用移液器小心吸出96孔PCR板内的液体，注意不要触碰磁珠。向每个样品孔中加入80%乙醇100 µl洗涤样品30 s，吸出乙醇，重复洗涤2次。将96孔PCR板取下磁力架，置于操作台上晾干5 min至酒精挥发完全。
6.  洗脱纯化的DNA。向每个样品孔中加入50 µl 1x TE缓冲液，吹吸混匀至管壁上无黏附的磁珠。室温放置5 min。将96孔PCR板置于磁力架上5 min使磁珠吸附在管壁。将96孔PCR板内液体吸出转移至新板，吹吸混匀。纯化后的DNA置于-20 °C环境保存。
7.  第二轮PCR模板DNA的稀释。使用分光光度计对磁珠纯化后的DNA进行浓度测定，并使用无核酸水将DNA稀释至10 ng/µl。
8.  第二轮PCR扩增。第二轮PCR扩增在反向引物1193R后端添加了不同的index序列用以区分不同文库。使用PrimeSTAR热稳定DNA聚合酶体系扩增，扩增体系如表 4所示。所有样本平行扩增3份，同时每个样本设立一个使用第一轮PCR阴性对照的PCR产物替代模板DNA的阴性对照，用于指示使用的实验试剂及实验过程是否存在污染。扩增条件如表 5所示。扩增产物置于-20 °C环境保存。
   
   表4.  第二轮PCR扩增体系
   

   试剂名称	添加体积 (µl)
   Nuclease-Free Water	13.8
   5x PrimeSTAR Buffer (Mg2+Plus)	6
   dNTP Mixture	2.4
   正反向引物	0.75
   10 ng/µl的DNA模板	6
   PrimeSTAR HS DNA Polymerase	0.3

   
   表5. 第二轮PCR扩增条件
   

   步骤	温度 (°C)	时间
   1	98	30 s
   2	98	10 s
   3	55	15 s
   4	72	1 min
   重复步骤2至4共8个循环
   5	72	5 min
   6	12	保温

   
   
9.  通过琼脂糖电泳进行第二轮PCR的质量控制。将平行扩增的样本混合均匀后吸取5 µl样本与1 µl 6x Loading Buffer混合加入1%琼脂糖凝胶中进行电泳。质量检测标准同第一轮PCR。
   
   

五、PCR产物的纯化、混合及Illumina测序

1.  利用琼脂糖电泳纯化PCR产物。因不同的样本已添加了不同的barcode及index序列，可以合并样本进行纯化回收。将同一index的文库中的样本每3个各取25 µl第二轮PCR产物混合，添加15 µl 6x Loading Buffer，吹吸混匀。置于1.2%的琼脂糖电泳凝胶中，125 V电压下电泳约40 min，至条带距离胶孔约5 cm。
2.  切胶并溶解。将400 bp位置的条带用干净的手术刀切下，置于2 ml离心管内，称量胶块的重量，以每10 mg凝胶添加10 µl的琼脂糖电泳凝胶纯化回收试剂盒中的Embrane Binding Solution的比例向离心管内添加Embrane Binding Solution。置于50-65 °C金属水浴锅中10 min至胶块溶解，期间每间隔2 min取出离心管颠倒混匀以加速胶块溶解。
3.  分离胶液并洗涤DNA。将融化的胶液冷却至室温，添加至带有SV Minicolumn过滤柱的Collection Tube中，孵育1 min。16,000 rcf离心1 min。弃去Collection Tube中滤液。加入700 µl Membrane Wash Solution (需按试剂盒要求添加乙醇后使用)。16,000 rcf离心1 min。弃去Collection Tube中滤液。再次加入500 µl Membrane Wash Solution。16,000 rcf离心5 min。弃去Collection Tube，将SV Minicolumn置入新的1.5 ml离心管，开盖干燥1 min以挥发酒精。
4.  洗脱纯化的DNA。向离心管内添加50 µl 无核酸水，静置3 min。16,000 rcf离心1 min。弃去SV Minicolumn。得到的纯化DNA置于-80 °C环境保存。
5.  切胶回收DNA的定量。使用PicoGreen DNA定量试剂盒对胶回收产物进行定量，方法参照步骤20-26。
6.  混合同一文库样本。带有同一index的样本按照每个胶回收产物200 ng DNA的量混合入1.5 ml离心管内。
7.  文库样本的纯化。使用Beckmen Coulter Agencourt AMPure XP核酸纯化试剂盒对文库DNA进行纯化以获得纯净的文库DNA样本，方法参照第一轮PCR产物的纯化。注意体系更换为1.5 ml离心管，应更换对应磁力架；DNA样品和磁珠的添加比例从1:1调整为1:0.9；使用80%酒精洗涤时体积调整为200 µl以没过磁珠；本次纯化进行两次以得到尽可能纯净的文库样本。纯化后的文库样本置于-80 °C环境保存。
8. 测量各文库样本的DNA浓度。使用分光光度计测量混合好并纯化了的的文库样本的DNA浓度。
9.  测序样本的混合。按照每个文库1200 ng DNA的量混合所有的待测序文库至1.5 ml离心管。若部分文库样品不满60个，应按比例减少其混样量。单次测序至少需要1500 ng DNA，测序由测序供应商处理。每个文库在Illumina HiSeq2500或NovaSeq6000平台采用双端250 bp模式测序，数据量 > 3 GB，即平均每个样本10万条读长 (10万条/样 × 60个样 × 250 bp × 2端 = 3 GB)。剩余文库样本和混合测序样本置于-80 °C环境保存。

溶液配方

1.  10x PBS (Phosphate Buffer Saline) 存储液：1.3 M NaCl、70 mM Na₂HPO₄和30 mM NaH₂PO₄ (pH = 7)。121 °C高压灭菌15 min，可在室温 (22-25 °C) 存储2个月。
2.  1x PBS 工作液：用无菌去离子水将10x PBS 存储液稀释10倍。可在室温下存储1个月。
3.  无菌去离子水：1 L去离子水装入1 L试剂瓶内，121 °C高压灭菌15 min。在室温条件下可储存1个月。
4.  2.5%有效氯浓度的次氯酸钠溶液：将200 ml 次氯酸钠溶液溶解于120 ml无菌去离子水中。现用现配。
5.  70%乙醇溶液：将70 ml无水乙醇溶解至30 ml无菌去离子水中。现用现配。
6.  80%乙醇溶液：将80 ml无水乙醇溶解至20 ml无菌去离子水中。现用现配。
7.  1x TAE：将20 ml 50x TAE溶解至980 ml无菌去离子水中。在室温条件下可储存1个月。
8.  1x TE Buffer：将5 ml 20x TE Buffer溶解至95 ml无菌去离子水中。在-4 °C条件下可储存1个月。
   

失败经验

1.  控制环境微生物污染。微生物在环境中普遍存在的特性导致实验过程中样本容易受到环境微生物的污染。因此所有操作应尽量在超净工作台内完成，收样过程等不能在超净工作台内完成的步骤应使用70%酒精对实验环境进行消毒；流程中所使用的所有耗材应进行高温高压灭菌并烘干后使用，实验过程中应全程佩戴手套并使用70%酒精对手套消毒；所有流程使用的试剂应保证无菌。
2. 控制样本的交叉污染。收样过程中不同处理的样本易产生交叉污染，应在更换不同处理的样本时，使用新的或消毒所用耗材 (如手套、剪刀、镊子、铲子等)。
3.  防止DNA提取过程引入差异。样本破碎质量会影响DNA提取的整体性，进而在后续数据中引入差异。应保证所有样本破碎至匀浆状态、无可见大颗粒，保证DNA提取的完整。
4.  防止文库构建过程引入差异。为防止PCR过程中引物偏好性引入差异，PCR过程进行3组平行扩增并混合；为防止扩增效率不同引入差异，同类样本的PCR产物浓度、电泳条带亮度应在一定范围内保持一致，避免出现2倍以上的差异；为防止切胶回收过程引入差异，应保证所有切下的胶块重量差异在10 mg内，并且切下的胶块在凝胶上位于同一水平线上。
   

致谢

本项目由中国科学院战略先导专项 (编号：XDA24020104)、中国科学院前沿科学重点研究项目 (编号：QYZDB-SSW-SMC021)、国家自然科学基金项目 (编号：31772400, 31761143017, 31801945, 31701997) 和中国科学院青年创新促进会 (编号：2020101, 2021092) [Supported by the Strategic Priority Research Program of the Chinese Academy of Sciences (Precision Seed Design and Breeding, No. XDA24020104), the Key Research Program of Frontier Sciences of the Chinese Academy of Science (No. QYZDB-SSW-SMC021), the National Natural Science Foundation of China (No. 31772400, 31761143017, 31801945, 31701997), the Youth Innovation Promotion Association CAS (No. 2020101, 2021092)] 支持。此实验方法已经在Nature Biotechnology (Zhang et al., 2019) 等国际顶级期刊发表文章中使用。

参考文献

1. Zhang, J. Y., Liu, Y. X., Zhang, N., Hu, B., Jin, T., Xu, H. R., Qin, Y., Yan, P. X., Zhang, X. N., Guo, X. X., Hui, J., Cao, S. Y., Wang, X., Wang, C., Wang, H., Qu, B. Y., Fan, G. Y., Yuan, L. X., Garrido-Oter, R., Chu, C. C. and Bai, Y. (2019). NRT1.1B is associated with root microbiota composition and nitrogen use in field-grown rice. Nat Biotechnol 6: 676-689.
2. Paloma, D., Thorsten, T., Ruben, G., Matthew, A., Eric, K., Paul S. L. and Stephane, H. (2018). Microbial interkingdom interactions in roots promote Arabidopsis survival. Cell 175: 973-983.

Copyright: © 2021 The Authors; exclusive licensee Bio-protocol LLC.

引用格式：徐浩然, 张婧赢, 郭晓璇, 曲宝原, 刘永鑫, 白洋. (2021). 水稻根系微生物组研究中的样本种植、取样和16S rRNA 基因扩增子文库制备方法. // 微生物组实验手册. Bio-101: e2003697. DOI: 10.21769/BioProtoc.2003697.
How to cite: Xu, H. R., Zhang, J. Y., Guo, X. X., Qu, B. Y., Liu, Y. X. and Bai, Y.  (2021). Planting, Sample Collection and Library Preparation of 16S rRNA Gene Amplicon Sequencing for Rice Root Microbiomes . // Microbiome Protocols eBook. Bio-101: e2003697. DOI: 10.21769/BioProtoc.2003697.