Name: 使用微生物分子特征评估水力压裂对溪流的影响 (视频) | JoVE Journal
Uploaded: 2021-04-04T13:00:00.000+00:00
Duration: 9 min 11 s
Description: 了解利用微生物分子特征评估水力压裂对溪流生态系统影响的方法。本方案涵盖从样本采集、核酸提取到16S rRNA测序的完整流程，为环境微生物学研究提供标准化技术支持。

ai_agent_IP

ai_agent

maintenance

instantTrialDurationMinutes

isBiopharma

enableNewMyJoveSectionDesign

educationFeatures

seoMetaApi

ai_copilot_institution_blacklist

ai_chat_institution_blacklist

freeAccessWindow

ssoRoleCapture

scheduledMaintenanceBanner

joveLabsFeatures

disable_autoplay_for_bots

extended_trial_duration

5 minutes

作者希望感谢为开发这些方法的相关项目提供资助的机构，包括：通过本科前及本科科学教育项目提供资助的霍华德·休斯医学研究所（http://www.hhmi.org），以及通过资助编号为 DBI-1248096 和 CBET-1805549 的项目提供支持的&#160;国家科学基金会（http://www.nsf.gov）&#160;。

1. 核酸提取用沉积物样品的采集
<ol>
	<li>将一个无菌的50 mL锥形管浸入溪流水中。采样时应佩戴手套，以避免引入人为污染。此步骤可在岸边进行，或在水中时面向上游方向操作。</li>
	<li>保持锥形管浸没于水中时，取下管盖，用其从水底1至3 cm深度处舀取约3 mL的沉积物至锥形管中。</li>
	<li>将锥形管从水中取出，倒掉其中所有上清液，仅保留覆盖在沉积物表面的一薄层液体（约1 mL）。</li>
	<li>使用1000 µL移液器及配套的移液器吸头，向采集的样品中加入3 mL DNA/RNA保存液（保存液的具体规格见材料表）。移液器吸头应存放于无菌的吸头盒中，仅在使用前立即安装，并在使用后丢弃。盖紧管盖后，上下颠倒混匀10次，确保保存液与样品充分混合。 
	注意：步骤1.4并非必需，但如果后续需从沉积物中提取RNA，则强烈建议执行此步骤。</li>
	<li>在完成后续采样期间，将样品置于冰上。采样结束后，若样品将用于16S分析（DNA），则储存于-20 °C冰箱中；若用于宏转录组分析（RNA），则储存于-70 °C冰箱中。</li>
</ol>
2. 核酸提取的滤膜收集
<ol>
	<li>取下一个无菌的1 L瓶子的瓶盖。面对上游或从岸边开始，将溪水灌满瓶子至顶部，然后倒掉。重复此过程两次，以润洗瓶子。第四次将瓶子完全装满后盖上瓶盖。 
	注意：如果重复使用1 L瓶子，可通过以下步骤进行灭菌：用10%次氯酸钠溶液冲洗2分钟，接着用去离子水冲洗三次，再用70%乙醇冲洗一次，最后在121.1 °C条件下高压灭菌30分钟，并保持15分钟干燥时间。高压灭菌过程中，瓶盖应拧得非常松，以避免瓶子在过程中受压变形。</li>
	<li>将瓶子置于稳定台面后，使用无菌鲁尔锁式注射器吸取满体积液体。然后将注射器连接至一个无菌且无DNA/RNA的聚醚砜滤膜（直径1.7 cm，孔径0.22 µm），通过完全压下活塞将全部液体推过滤膜。重复此操作，直至瓶中收集的总体积（1 L）全部通过滤膜。 
	注意：注射器体积可灵活选择，只要准确记录通过滤膜的总水量即可。通常推荐使用60 mL注射器。虽然1 L为理想体积，但根据经验，至少200 mL的水量（假设每毫升约含20,000个细胞）仍可能足以用于后续DNA和RNA的提取。</li>
	<li>用注射器吸入约20 mL空气，并将其推过滤膜，以去除滤膜中多余的水分。 
	注意：此步骤有助于防止在执行步骤2.4时保存剂的损失。</li>
	<li>使用P1000移液器，将2 mL DNA/RNA保存剂通过滤膜较大的开口（即连接注射器的一端）注入，操作时保持滤膜水平放置。按下移液器时，移液枪头应位于滤膜筒内部，以确保保存剂进入滤膜内部。每次使用后更换枪头。 
	注意：与沉积物采集类似，此步骤并非必需，但强烈建议执行，以提高后续核酸得率，尤其是RNA。</li>
	<li>撕下一块石蜡膜，紧密包裹滤膜的每个开口/两端以密封。将包裹好的滤膜放入无菌样品袋中，并在采样期间将整个样品袋置于冰上保存。 
	注意：确保用于包裹滤膜的石蜡膜表面无菌，即此前未暴露于环境中。</li>
	<li>采样返回后，将滤膜在-20 °C保存用于16S分析，或在-70 °C保存用于宏转录组分析。</li>
</ol>
3. 核酸提取与定量
<ol>
	<li>在开始样品转移前，使用10%次氯酸钠溶液和70%乙醇清洁工作区域。</li>
	<li>对于来自步骤1.5的沉积物，通常使用约0.25 g样品。将金属工具浸入盛有70%乙醇的烧杯中，通过火焰灼烧乙醇实现工具的灼烧灭菌，每处理一个样品之间均需进行此操作。</li>
	<li>对于滤膜（来自步骤2.6），将滤纸转移至无菌管中用于提取。操作步骤如下。
	<ol>
		<li>通过折叠铝箔以确保其内层不暴露于外部环境，然后在121.1 °C条件下高压灭菌并设置5分钟干燥时间，从而制备无菌、无DNA和无RNA的表面。</li>
		<li>用70%乙醇和明火对尖嘴钳进行灭菌。然后在无菌表面上用尖嘴钳打开滤器外壳，并从外壳中取出滤芯。</li>
		<li>使用无菌手术刀在芯体的顶部和底部切开，然后沿接缝处切割，将滤纸从芯体上分离。用无菌镊子折叠滤纸，再用手术刀将滤纸切成小片。</li>
		<li>将滤纸片放入微量离心管中用于提取。确保滤纸不接触任何未灭菌或可能含有核酸的表面，以免导致样本发生非预期的污染。</li>
	</ol>
	</li>
	<li>按照先前所述方法进行DNA提取13 或使用市售的基于柱层析的试剂盒（参见 材料表）。以下简要描述所列商用试剂盒的操作步骤。
	<ol>
		<li>将样品转移至珠管中，使用细胞破碎仪高速处理至少5分钟以裂解细胞。离心后，将上清液转移至无菌微量离心管中。</li>
		<li>向上清液中加入裂解缓冲液（体积比1:1），转移至提供的滤膜（黄色）中，离心过滤。</li>
		<li>将滤膜转移至新的无菌微量离心管中，加入400 µL预处理缓冲液，离心后弃去滤过液。</li>
		<li>加入洗涤缓冲液（700 µL），离心，弃去流出液。然后加入洗涤缓冲液（400 µL），离心，再次弃去流出液。</li>
		<li>将滤膜转移至新的无菌微量离心管中，加入50 µL无DNase/RNase水，室温静置5分钟，然后离心。</li>
		<li>在孵育期间，将III-HRC滤芯放入收集管中，加入HRC预处理液（600 µL），随后以8,000 x g离心3 min。 g.</li>
		<li>将制备好的滤膜移至一支无菌的微量离心管中。将步骤 3.4.5 中洗脱的 DNA 转移至该滤膜上，于 16,000 x g 离心 g 3 分钟。流出液中含有提取的 DNA。</li>
	</ol>
	</li>
	<li>将沉积物和滤膜的DNA提取物在-20 °C下保存。 
	注意：假设温度稳定、避光且无有害污染物，DNA 提取物可在 -20 °C 下保存约 8 年14.</li>
	<li>按照制造商的说明书进行RNA提取。将RNA提取物保存在-80 °C。
	<ol>
		<li>将样品转移至珠管中，使用细胞破碎仪高速处理至少五分钟以裂解细胞。离心后，将上清液转移至无菌微量离心管中。</li>
		<li>向 supernatant 中加入裂解缓冲液（体积比1:1），转移至提供的离心柱（黄色）中。离心柱离心。</li>
		<li>加入等体积的95-100%乙醇至流穿液中，用移液器上下吹打混匀五次。</li>
		<li>将IICG柱（绿色）置于无菌微量离心管上。将混合溶液转移至柱中，离心。</li>
		<li>加入洗涤缓冲液（400 µL），离心后弃去滤过液。</li>
		<li>向柱中加入5 µL DNase I和75 µL DNA消化缓冲液，在室温下孵育15分钟。</li>
		<li>加入预平衡缓冲液（400 µL），离心，弃去流出液。</li>
		<li>加入洗涤缓冲液（700 µL），离心，弃去流出液。然后加入洗涤缓冲液（400 µL），离心，再次弃去流出液。</li>
		<li>将柱子转移至新的无菌微量离心管中，加入50 µL无DNase/RNase水，静置5分钟，然后离心洗脱。</li>
		<li>在孵育期间，将III-HRC滤芯放入收集管中，加入HRC预处理液（600 µL），随后以8,000 x g离心3分钟，完成滤芯的准备。 g.</li>
		<li>将制备好的滤器转移至一支无菌的微量离心管中。将步骤 3.6.9 中洗脱的 RNA 加入此滤器中，于 16,000 x g 离心 g 3 分钟。流出液中含有提取的 RNA。 
		注意：RNA 提取物在开始降解前最多只能保存一年15DNA和RNA提取物在反复冻融过程中会被降解。某些实验方案允许从同一样本中同时提取DNA和RNA。16,17.</li>
	</ol>
	</li>
	<li>使用荧光计或分光光度计对提取的DNA和RNA样本进行定量。参见 表1 例如荧光计测定的DNA浓度值。有关分光光度计定量方案的示例，请参见参考文献18使用所列试剂盒的沉积物DNA浓度值 材料表 通常在 1 至 40 ng/µL 之间，而滤膜 DNA 浓度值通常在 0.5 至 10 ng/µL 之间。使用清单中所列试剂盒测得的沉积物 RNA 浓度值 材料表 通常在1至20 ng/µL之间，而滤膜RNA的浓度值往往较低，通常在0.5至5 ng/µL之间。</li>
</ol>
4.&#160;16S rRNA 文库构建
<ol>
	<li>使用10%次氯酸钠和70%乙醇清洁工作区域。工作区域应为能够产生层流条件的封闭空间（层流罩）。</li>
	<li>使用DNA提取物（来自步骤3.5），在层流条件下，按照标准PCR方案制备用于16S rRNA扩增子测序的样品，例如Earth Microbiome网站所述的扩增16S rRNA V4高变区的方案19。</li>
	<li>按先前所述方法制备2%琼脂糖凝胶，并使其凝固17。将7 µL PCR产物与13 µL无DNase水混合，加入凝胶上样染料至终浓度为1×。待琼脂糖凝固后，将该PCR产物混合液上样至2%琼脂糖凝胶。 
	注意：也可使用预制胶，因其运行速度更快且已预先制备完成。</li>
	<li>以90 V电压电泳60–90分钟，通过检测约386 bp条带的存在判断16S rRNA V4扩增子是否成功扩增，具体参照Earth Microbiome的实验方案。</li>
</ol>
5. DNA 16S rRNA 文库纯化
<ol>
	<li>将产生明亮条带的样品 PCR 产物各取 10 µL，产生微弱条带的样品 PCR 产物各取 13 µL，汇集至适当大小的无菌微量离心管中。</li>
	<li>使用荧光计或分光光度计检测所得混合物的浓度，并按前述方法制备 2% 琼脂糖凝胶。理想情况下，混合物的浓度应至少达到 10 ng/µL，大多数样品的浓度应约为 25 ng/µL。</li>
	<li>在浓度和体积允许的情况下，将约 150–200 ng 的样品上样至 2% 琼脂糖凝胶的加样孔中。</li>
	<li>在 90 伏电压下电泳运行 60–90 分钟。</li>
	<li>通过运行 2% 琼脂糖凝胶对混合文库进行纯化。
	<ol>
		<li>从凝胶中切下 386 bp 的 DNA 条带，并按照先前所述方法使用 commercially available kit（商业化试剂盒）纯化混合文库20。用 30 µL 的 10 mM Tris-Cl（pH 8.5）洗脱纯化后的 DNA。此步骤应在与 DNA 或 RNA 提取不同的区域进行，以防止后续污染，因为切胶操作会将 PCR 扩增子扩散至实验人员及周围环境中。</li>
	</ol>
	</li>
	<li>使用荧光计或分光光度计检测纯化后混合物的浓度。若纯化效果良好，其浓度应至少为纯化前混合物的一半。通常，最终浓度应在 5 至 20 ng/µL 范围内。</li>
	<li>将纯化后的文库送交进行高通量测序（下一代测序）。确保运输过程中保持低温，需在运输容器中加入干冰。</li>
</ol>
6. RNA文库构建与纯化
<ol>
	<li>可使用多种商业试剂盒构建RNA文库。无论使用哪一种试剂盒，均需在无菌的层流环境中严格按照制造商提供的说明书进行操作。下文简要概述了材料表中所列试剂盒的操作流程21。
	<ol>
		<li>配制第一链cDNA合成预混液（8 µL无核酸酶水和2 µL第一链合成酶混合液），加入样本中。将样本放入PCR仪，按照说明书指定的条件进行反应。</li>
		<li>在冰上配制第二链cDNA合成预混液（8 µL第二链合成反应缓冲液、4 µL第二链合成酶混合液和48 µL无核酸酶水），加入样本中。将样本置于设定为16 °C的PCR仪中反应1小时。</li>
		<li>加入提供的磁珠（144 µL）纯化反应产物，并用80%乙醇洗涤两次（每次200 µL）。</li>
		<li>用提供的TE缓冲液（53 µL）洗脱，转移50 µL上清液至新的洁净PCR管中。将PCR管置于冰上。</li>
		<li>在冰上配制末端修复预混液（7 µL末端修复反应缓冲液和3 µL末端修复酶混合液），加入PCR管中。将PCR管放入PCR仪，按照说明书指定的条件进行反应。</li>
		<li>在冰上混合稀释的接头（2.5 µL）、连接预混液（30 µL）和连接增强剂（1 µL）。将混合液加入样本中，置于PCR仪中在20 °C反应15分钟。</li>
		<li>加入提供的磁珠（87 µL）纯化反应产物，用乙醇洗涤（200 µL）并洗脱，操作同前，但仅加入17 µL TE缓冲液。</li>
		<li>加入index引物（10 µL）和Q5预混液（25 µL），放入PCR仪，按照说明书指定的条件进行扩增。</li>
		<li>加入提供的磁珠（45 µL）纯化反应产物，再进行两次乙醇洗涤（每次200 µL），用23 µL TE缓冲液洗脱。转移20 µL至新的洁净PCR管中。</li>
	</ol>
	</li>
	<li>使用生物分析仪、荧光计或分光光度计检测文库中RNA的可检测浓度。</li>
	<li>将宏转录组文库按大致等摩尔比例混合。</li>
	<li>按照16S文库纯化的相同流程对文库进行纯化，但切取250至400 bp之间的片段。与16S文库具有明显条带不同，此处结果为弥散条带（smear）。</li>
	<li>如前所述，检测纯化后文库的浓度。</li>
	<li>将纯化后的文库用干冰运输至测序机构。 
	&#8203;注意：也可将RNA提取物送至大学或私营公司进行文库构建和测序。</li>
</ol>
7. 微生物群落分析
<ol>
	<li>测序完成后，获取样本数据，并将其下载到可用的计算机上。 
	注意：理想情况下，设备应至少配备 16 GB 的内存。有关计算需求（针对 Qiime2）的讨论，请参见 https://forum.qiime2.org/t/recommended-specifications-to-run-qiime2/9808。</li>
	<li>使用 mothur、QIIME2 和 R 等软件分析 16S rRNA 数据。示例 QIIME2 16S 分析教程请参见 https://docs.qiime2.org/2020.11/tutorials/moving-pictures/。</li>
	<li>对于宏转录组学（RNA）数据，使用 HUMAnN2 和 ATLAS 确定样本中存在的基因和通路。 
	注意：示例宏转录组学分析流程（最终进行多样性分析和随机森林分析）见于补充信息文件中。所有命令均通过命令行运行，例如 Mac 用户使用 Terminal。</li>
</ol>

<ol><li>The process of unconventional natural gas production. US EPA. , Available from: 
 <a target="_blank" href="https://www.epa.gov/uog/process-unconventional-natural-gas-production">https://www.epa.gov/uog/process-unconventional-natural-gas-production</a> (2013).</li>
<li>Brittingham, M. C., Maloney, K. O., Farag, A. M., Harper, D. D., Bowen, Z. H. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Ecological+risks+of+shale+oil+and+gas+development+to+wildlife%2C+aquatic+resources%2C+and+their+habitats%5BTitle%5D)+AND+%22Environmental+Science+%26+Technology%22%5BJournal%5D)">Ecological risks of shale oil and gas development to wildlife, aquatic resources, and their habitats.</a> Environmental Science & Technology. 48 (19), 11034-11047 (2014).</li>
<li>McBroom, M., Thomas, T., Zhang, Y. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Soil+erosion+and+surface+water+quality+impacts+of+natural+gas+development+in+East+Texas%2C+USA%5BTitle%5D)+AND+%22Water%22%5BJournal%5D)">Soil erosion and surface water quality impacts of natural gas development in East Texas, USA.</a> Water. 4 (4), 944-958 (2012).</li>
<li>Maloney, K. O., Weller, D. E. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Anthropogenic+disturbance%2C+and+streams%3A+land+use+and+land-use+change+affect+stream+ecosystems+via+multiple+pathways%5BTitle%5D)+AND+%22Freshwater+Biology%22%5BJournal%5D)">Anthropogenic disturbance, and streams: land use and land-use change affect stream ecosystems via multiple pathways.</a> Freshwater Biology. 56 (3), 611-626 (2011).</li>
<li>Meyer, J. L., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((The+contribution+of+headwater+streams+to+biodiversity+in+river+networks1%5BTitle%5D)+AND+%22JAWRA+Journal+of+the+American+Water+Resources+Association%22%5BJournal%5D)">The contribution of headwater streams to biodiversity in river networks1.</a> JAWRA Journal of the American Water Resources Association. 43 (1), 86-103 (2007).</li>
<li>Alexander, R. B., Boyer, E. W., Smith, R. A., Schwarz, G. E., Moore, R. B. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((The+role+of+headwater+streams+in+downstream+water+quality%5BTitle%5D)+AND+%22Journal+of+the+American+Water+Resources+Association%22%5BJournal%5D)">The role of headwater streams in downstream water quality.</a> Journal of the American Water Resources Association. 43 (1), 41-59 (2007).</li>
<li>Ulrich, N., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Response+of+aquatic+bacterial+communities+to+hydraulic+fracturing+in+Northwestern+Pennsylvania%3A+A+five-year+study%5BTitle%5D)+AND+%22Scientific+Reports%22%5BJournal%5D)">Response of aquatic bacterial communities to hydraulic fracturing in Northwestern Pennsylvania: A five-year study.</a> Scientific Reports. 8 (1), 5683(2018).</li>
<li>Chen See, J. R., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Bacterial+biomarkers+of+Marcellus+shale+activity+in+Pennsylvania%5BTitle%5D)+AND+%22Frontiers+in+Microbiology%22%5BJournal%5D)">Bacterial biomarkers of Marcellus shale activity in Pennsylvania.</a> Frontiers in Microbiology. 9, 1697(2018).</li>
<li>Rausch, P., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Comparative+analysis+of+amplicon+and+metagenomic+sequencing+methods+reveals+key+features+in+the+evolution+of+animal+metaorganisms%5BTitle%5D)+AND+%22Microbiome%22%5BJournal%5D)">Comparative analysis of amplicon and metagenomic sequencing methods reveals key features in the evolution of animal metaorganisms.</a> Microbiome. 7 (1), 133(2019).</li>
<li>Louca, S., Doebeli, M., Parfrey, L. W. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Correcting+for+16S+rRNA+gene+copy+numbers+in+microbiome+surveys+remains+an+unsolved+problem%5BTitle%5D)+AND+%22Microbiome%22%5BJournal%5D)">Correcting for 16S rRNA gene copy numbers in microbiome surveys remains an unsolved problem.</a> Microbiome. 6 (1), 41(2018).</li>
<li>Trexler, R., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Assessing+impacts+of+unconventional+natural+gas+extraction+on+microbial+communities+in+headwater+stream+ecosystems+in+Northwestern+Pennsylvania%5BTitle%5D)+AND+%22Frontiers+in+Microbiology%22%5BJournal%5D)">Assessing impacts of unconventional natural gas extraction on microbial communities in headwater stream ecosystems in Northwestern Pennsylvania.</a> Frontiers in Microbiology. 5, 522(2014).</li>
<li>Mumford, A. C., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Shale+gas+development+has+limited+effects+on+stream+biology+and+geochemistry+in+a+gradient-based%2C+multiparameter+study+in+Pennsylvania%5BTitle%5D)+AND+%22Proceedings+of+the+National+Academy+of+Sciences%22%5BJournal%5D)">Shale gas development has limited effects on stream biology and geochemistry in a gradient-based, multiparameter study in Pennsylvania.</a> Proceedings of the National Academy of Sciences. 117 (7), 3670-3677 (2020).</li>
<li>JoVE Core Biology DNA Isolation. Journal of Visualized Experiments. , Cambridge, MA. Available from: 
 <a target="_blank" href="https://www.jove.com/cn/science-education/10814/dna-isolation">https://www.jove.com/cn/science-education/10814/dna-isolation</a> (2020).</li>
<li>Oxford Gene Technology DNA Storage and Quality. OGT. , Available from: 
 <a target="_blank" href="https://www.ogt.com/resources/literature/403_dna_storage_and_quality">https://www.ogt.com/resources/literature/403_dna_storage_and_quality</a> (2011).</li>
<li>ThermoFisher SCIENTIFIC Technical Bulletin #159: Working with RNA. Thermoscientific. , Available from: 
 <a target="_blank" href="https://www.thermofisher.com/us/en/home/references/ambion-tech-support/nuclease-enzymes/general-articles/working-with-rna.html">https://www.thermofisher.com/us/en/home/references/ambion-tech-support/nuclease-enzymes/general-articles/working-with-rna.html</a> (2020).</li>
<li>QIAGEN AllPrep DNA/RNA Mini Kit. Qiagen. , Available from: 
 <a target="_blank" href="https://www.qiagen.com/us/products/discovery-and-translational-research/dna-rna-purification/multianalyte-and-virus/allprep-dnarna-mini-kit/#orderinginformation">https://www.qiagen.com/us/products/discovery-and-translational-research/dna-rna-purification/multianalyte-and-virus/allprep-dnarna-mini-kit/#orderinginformation</a> (2020).</li>
<li>ZymoBIOMICS DNA/RNA Miniprep Kit. Zymo Research. , Available from: 
 <a target="_blank" href="https://www.zymoresearch.com/products/zymobiomics-dna-rna-miniprep-kit">https://www.zymoresearch.com/products/zymobiomics-dna-rna-miniprep-kit</a> (2020).</li>
<li>Desjardins, P., Conklin, D. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((NanoDrop+microvolume+quantitation+of+nucleic+acids%5BTitle%5D)+AND+%22Journal+of+Visualized+Experiments%22%5BJournal%5D)">NanoDrop microvolume quantitation of nucleic acids.</a> Journal of Visualized Experiments. (45), e2565(2010).</li>
<li>16S Illumina amplicon protocol: Earth microbiome project. Earth microbiome project. , Available from: 
 <a target="_blank" href="https://earthmicrobiome.org/protocols-and-standards/16s/">https://earthmicrobiome.org/protocols-and-standards/16s/</a> (2018).</li>
<li>Gel Purification: Binding, washing and eluting a sample | Protocol. Journal of Visualized Experiments. , Available from: 
 <a target="_blank" href="https://www.jove.com/v/5063/gel-purification">https://www.jove.com/v/5063/gel-purification</a> (2020).</li>
<li>New England Biolabs protocol for the use with NEBNext Poly(A) mRNA magnetic isolation module (E7490) and NEBNext Ultra II RNA library prep kit for Illumina (E7770, E7775). New England Biolabs. , Available from: 
 <a target="_blank" href="https://www.neb.com/protocols/2017/03/04/protocol-for-use-with-purified-mrna-or-rrna-depleted-rna-and-nebnext-ultra-ii-rna-library-prep-ki">https://www.neb.com/protocols/2017/03/04/protocol-for-use-with-purified-mrna-or-rrna-depleted-rna-and-nebnext-ultra-ii-rna-library-prep-ki</a> (2020).</li>
<li>Anderson, M. J. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Permutational+multivariate+analysis+of+variance+%28PERMANOVA%29%5BTitle%5D)+AND+%22Wiley+StatsRef%3A+Statistics+Reference+Online%22%5BJournal%5D)">Permutational multivariate analysis of variance (PERMANOVA).</a> Wiley StatsRef: Statistics Reference Online. , 1-15 (2017).</li>
<li>Schrader, C., Schielke, A., Ellerbroek, L., Johne, R. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((PCR+inhibitors+-+occurrence%2C+properties+and+removal%5BTitle%5D)+AND+%22Journal+of+Applied+Microbiology%22%5BJournal%5D)">PCR inhibitors - occurrence, properties and removal.</a> Journal of Applied Microbiology. 113 (5), 1014-1026 (2012).</li>
</ol>

作者无任何利益冲突需要披露。

本文所述方法是在2014年至2018年间由我研究团队发表的多项研究中逐步开发和完善而成的7,8,10，并已成功应用于一项为期三年的协作项目中，用于研究水力压裂对水生群落的影响，该项目即将提交一篇待发表的论文。这些方法将在项目剩余时间内继续使用。此外，目前其他研究水力压裂对溪流和生态系统影响的文献也描述了类似的样品采集、处理和分析方法7,8,10,11。然而，这些文献均未采用宏转录组学分析，因此本文是首篇描述如何利用此类分析方法揭示水力压裂对邻近溪流影响的研究。此外，本文所提出的数据采集方法更为详尽，且在防止污染方面所采取的步骤也更加细致。
本方案中最重要的步骤之一是初始样品的采集与保存。野外采样和收集过程面临某些挑战，因为在采集过程中维持无菌或无菌环境可能较为困难。在此步骤中，避免样品污染至关重要。为此，应佩戴手套，并且仅允许无菌容器和工具接触样品。采集后，样品应立即置于冰上，以减缓核酸降解。在采集时添加商业化的核酸保存剂，也可提高核酸得率，并使样品在采集后能够长期保存。无论何时进行核酸提取，使用适量的样品都非常重要：样品过多可能会堵塞用于提取的离心过滤柱（针对使用此类柱子的方案），而样品过少则可能导致得率偏低。请务必遵循所用试剂盒的说明书操作。
与野外采样类似，在核酸提取和样品制备过程中，避免或尽量减少污染也非常重要，尤其是在处理核酸产量较低的样品时，例如非最佳沉积物样品（含有大量碎石或岩石的样品）或水样。因此，与样品采集一样，在所有这些步骤中都应佩戴手套，以减少污染。此外，实验过程中使用的所有工作台面应在操作前用10%次氯酸钠溶液擦拭消毒，随后再用70%乙醇溶液擦拭。在移液步骤（3–6）中，应使用滤芯吸头，以避免移液器本身带来的污染，并且每次接触非无菌表面后都应更换吸头。所有用于实验操作的工具（包括移液器）均应在使用前后用次氯酸钠和乙醇溶液进行擦拭。为评估污染情况，每次进行核酸提取和PCR反应时都应包含提取空白对照和阴性对照（无菌液体）。如果提取后的定量结果显示阴性对照中可检测到DNA/RNA，则在样品量充足的情况下可重复提取。如果PCR的阴性样品出现扩增，应进行故障排查以确定污染来源，然后重新运行样品。为应对低水平污染，建议对提取空白对照和PCR阴性对照进行测序，以便在必要时通过计算分析识别并去除污染物。相反，PCR扩增也可能因多种原因失败。对于环境样品，PCR反应受到抑制通常是主要原因，这可能是由于多种物质干扰了Taq DNA聚合酶23的活性。若怀疑存在抑制作用，可使用PCR级水（见材料表）稀释DNA提取物。
本方案存在一些显著的局限性和潜在困难。水样和沉积物样品的采集均可能具有挑战性。为了获得足够的生物量，理想情况下需要将1 L的溪流水样通过滤膜过滤。滤膜的孔径需足够小以捕获微生物，但也可能同时截留沉积物颗粒。若因近期降雨导致水体中含有大量沉积物，滤膜容易堵塞，难以将全部水样体积过滤通过。在采集沉积物时，准确估计采样深度也可能较为困难。此外，需确保所采集的沉积物主要成分为土壤，因为砾石和岩石会降低核酸得率，并且可能无法准确代表微生物群落的真实情况。最后，样品采集后必须始终保持在冰上保存，尤其是在未使用保存剂的情况下，这一点至关重要。
尽管本实验方案同时涵盖了宏转录组学和16S的实验操作流程，但需要强调的是，这两种方法在操作过程及其所提供的数据类型方面存在显著差异。16S rRNA基因是细菌和古菌中一个高度保守的常用靶向区域，适用于对样本中的细菌群落进行表征。然而，作为一种靶向性特定方法，其通常难以实现物种水平的分辨，且对新近分化的物种或菌株的鉴定较为困难。相比之下，宏转录组学是一种更广泛的方法，能够捕获样本中所有活跃表达的基因和微生物。16S仅能提供物种鉴定相关数据，而宏转录组学则可提供功能层面的信息，例如基因表达情况和代谢通路。两种方法均具有重要价值，联合使用时可揭示样本中存在哪些细菌以及这些细菌正在表达哪些基因。
本文介绍了在研究水力压裂（fracking）背景下，进行16S rRNA和宏转录组分析的野外样品采集与处理方法。此外，还详细描述了从生物量极低的样品中提取高质量DNA/RNA以及长期保存样品的采集方法。文中所述方法总结了我们在样品采集与处理方面的实践经验，旨在通过分析微生物群落的结构与功能，揭示水力压裂对周边溪流生态系统的影响。微生物对环境扰动响应迅速，因此，存在的微生物种类及其表达的基因可提供有关水力压裂对生态系统影响的重要信息。总体而言，这些方法对于深入理解水力压裂如何影响这些重要生态系统具有重要价值。

水力压裂（HF），或称&#34;压裂&#34;，是一种天然气开采方法，随着化石燃料需求持续上升，该方法正变得日益普遍。该技术利用高压钻探设备，将水、沙子和化学物质的混合物注入富含甲烷的页岩层中，通常用于释放被封存的气体1。
由于这些非常规的开采技术相对较新，因此研究此类作业对附近水道的影响至关重要。水力压裂作业需要清理大片土地，用于设备运输和井场建设。每个井场需清理约1.2至1.7公顷的土地2，可能影响系统的径流和水质3。关于压裂液确切的化学成分，包括所使用的杀菌剂种类，目前缺乏透明度。此外，压裂废水通常具有高盐度2。进一步地，废水中可能含有金属和天然存在的放射性物质2。因此，由于人为错误或设备故障导致压裂液泄漏或溢出的可能性令人担忧。
溪流生态系统对外部环境变化极为敏感4 对于维持生物多样性具有重要意义5 以及适当的养分循环6 在整个流域内，微生物是淡水中数量最丰富的生物类群，因此在养分循环、生物降解和初级生产过程中发挥着关键作用。由于微生物对环境扰动具有高度敏感性，其群落组成与功能可作为获取生态系统信息的重要工具；近期研究表明，根据与水力压裂活动的距离远近，所观测到的细菌群落结构呈现出显著差异。7,8例如， 贝日阿托氏菌属, 伯克霍尔德菌属，以及 Methanobacterium 在压裂作业附近的溪流中被发现富集 Pseudonocardia, Nitrospira，以及 Rhodobacter 在远离水力压裂作业的溪流中富集7.
对16S核糖体RNA（rRNA）基因进行高通量测序是一种经济高效的方法，用于确定细菌群落组成，其速度更快、成本低于全基因组测序方法9。在分子生态学领域，常用16S rRNA基因中高度可变的V4区进行测序，通常可在较广范围内实现属水平的鉴定分辨率9，因此特别适用于不可预测的环境样本。该技术已在大量已发表的研究中广泛应用，并成功用于识别水力压裂作业对水生环境的影响7,8。然而需要注意的是，不同细菌的16S rRNA基因拷贝数存在差异，这会影响其检测丰度10。尽管已有少数工具用于校正这一偏差，但其有效性仍存疑10。另一种迅速普及且无此缺陷的方法是宏转录组测序，该方法对全部RNA进行测序，使研究人员能够同时鉴定出活跃的细菌及其基因表达情况。
因此，与先前已发表研究中的方法7,8,11,12相比，本实验方案还涵盖了用于研究微生物群落功能（宏转录组学）的样品采集、保存、处理和分析。本文详述的步骤可使研究人员了解水力压裂对溪流中微生物所表达的基因和代谢通路（包括抗微生物药物耐药性基因）是否产生影响以及影响程度。此外，本方案在样品采集方面的细节描述更为完善。&#160;尽管其中一些步骤和注意事项对经验丰富的研究人员而言可能显而易见，但对于刚开始从事研究的人员来说，这些内容可能极为宝贵。
本文基于我们实验室多年的实践经验，介绍了用于采集和处理样本以获取细菌遗传数据的方法，旨在研究水力压裂对附近溪流的影响。这些数据可用于后续分析，以识别与水力压裂状态相关的差异。

<table><tbody><tr><td>200 证明乙醇</td><td>Thermo Fisher Scientific</td><td>A4094</td><td>需向 Buffer PE 中加入 400 mL（参见 Qiagen QIAQuck 凝胶提取试剂盒说明书），并向 DNA/RNA 洗涤缓冲液中加入 96 mL（参见 ZymoBIOMICS DNA/RNA Miniprep 试剂盒说明书）。&lt;br/&gt; ZymoBIOMICS DNA/RNA Miniprep 试剂盒和 NEBNext&amp;reg; Ultra&amp;trade; II RNA 文库构建试剂盒（含样本纯化磁珠）还需额外乙醇。</td></tr><tr><td>琼脂糖</td><td>Thermo Fisher Scientific</td><td>BP1356-100</td><td>每瓶 100 g。制备一块 30 mL 的 2% 凝胶需使用 0.6 g 琼脂糖。</td></tr><tr><td>消毒漂白剂</td><td>Walmart (Clorox)</td><td>无目录号</td><td>实验前后需使用 10% 漂白剂溶液清洁工作区域</td></tr><tr><td>DNA 凝胶上样染料</td><td>Thermo Fisher Scientific</td><td>R0611</td><td>每个自制（非 e-gel）样品均需加入上样染料，比例为每 5 &amp;micro;L 样品加入 1 &amp;micro;L 染料</td></tr><tr><td>DNA Marker</td><td>MilliporeSigma</td><td>D3937-1VL</td><td>每次电泳（包括 e-gel）均应运行 DNA Marker</td></tr><tr><td>DNA/RNA Shield (2x)</td><td>Zymo Research</td><td>R1200-125</td><td>每个沉积物样品（50 mL 离心管）需 3 mL，每个水样（滤膜）需 2 mL</td></tr><tr><td>溴化乙锭</td><td>Thermo Fisher Scientific</td><td>BP1302-10</td><td>用于染色自制 e-gel</td></tr><tr><td>正向引物</td><td>Integrated DNA Technologies (IDT)</td><td>51-01-19-06</td><td>每轮 PCR 反应需 0.5 &amp;micro;L</td></tr><tr><td>异丙醇</td><td>MilliporeSigma</td><td>563935-1L</td><td>通常每文库用量少于 2 mL。具体体积取决于切胶片段的质量（参见 Qiagen QIAQuick 凝胶提取试剂盒说明书）。&amp;nbsp;</td></tr><tr><td>PCR 级水</td><td>MilliporeSigma</td><td>3315932001</td><td>每轮 PCR 反应需 13 &amp;micro;L（假设使用 1 &amp;micro;L 样本 DNA 模板）</td></tr><tr><td>Platinum Hot Start PCR Master Mix (2x)</td><td>Thermo Fisher Scientific</td><td>13000012</td><td>每轮 PCR 反应需 10 &amp;micro;L</td></tr><tr><td>反向引物</td><td>Integrated DNA Technologies (IDT)</td><td>51-01-19-07</td><td>每轮 PCR 反应需 0.5 &amp;micro;L</td></tr><tr><td>TBE 缓冲液 (Tris-硼酸-EDTA)</td><td>Thermo Fisher Scientific</td><td>B52</td><td>需 1 L 10x TBE 缓冲液（制备一块 30 mL 凝胶需 30 mL 1x TBE 缓冲液）</td></tr><tr><td>1 L 瓶</td><td>Thermo Fisher Scientific</td><td>02-893-4E</td><td>每条溪流需一个（若在使用间进行灭菌，同一瓶子可用于多条溪流）</td></tr><tr><td>1.5 mL 微量离心管</td><td>MilliporeSigma</td><td>BR780400-450EA</td><td>每次 DNA 提取需 5 个微量离心管，额外需 3 个用于 RNA 纯化（参见 ZymoBIOMICS DNA/RNA Miniprep 试剂盒说明书）</td></tr><tr><td>2% 琼脂糖 e-gel</td><td>Thermo Fisher Scientific</td><td>G401002</td><td>每块凝胶可运行 10 个样品（因此若包含一个 PCR 阴性对照，则为 9 个；若同时运行提取阴性对照，则为 8 个）</td></tr><tr><td>50 mL 离心管</td><td>CellTreat</td><td>229421</td><td>每个沉积物样品需 1 个 50 mL 离心管</td></tr><tr><td>500 mL 烧杯</td><td>MilliporeSigma</td><td>Z740580</td><td>仅需 1 个（用于火焰灭菌）</td></tr><tr><td>铝箔</td><td>Walmart (Reynolds KITCHEN)</td><td>无编号</td><td>铝箔可折叠后高压灭菌。未暴露于环境的部分可用作处理滤膜时的无菌、无 DNA 和 RNA 污染的表面&lt;br/&gt;（每个滤膜使用一块折叠铝箔，以避免交叉污染）</td></tr><tr><td>高压灭菌器</td><td>Gettinge</td><td>LSS 130</td><td>仅需 1 台</td></tr><tr><td>离心机</td><td>MilliporeSigma</td><td>EP5404000138-1EA</td><td>仅需 1 台</td></tr><tr><td>保温箱</td><td>ULINE</td><td>S-22567</td><td>几乎任何保温箱均可使用。此处列出的型号因由泡沫制成而更轻，便于野外采样时携带。</td></tr><tr><td>Disruptor Genie</td><td>Bio-Rad</td><td>3591456</td><td>仅需 1 台</td></tr><tr><td>电泳槽</td><td>Bio-Rad</td><td>1664000EDU</td><td>仅需 1 台</td></tr><tr><td>电泳电源</td><td>Bio-Rad</td><td>1645050</td><td>仅需 1 台</td></tr><tr><td>冷冻柜 (-20 °C)</td><td>K2 SCIENTIFIC</td><td>K204SDF</td><td>需 1 台用于储存 DNA 提取物</td></tr><tr><td>超低温冰箱 (-80 °C)</td><td>K2 SCIENTIFIC</td><td>K205ULT</td><td>需 1 台用于储存 RNA 提取物</td></tr><tr><td>手套</td><td>Thermo Fisher Scientific</td><td>19-020-352</td><td>目录号对应中号手套。</td></tr><tr><td>加热块</td><td>MilliporeSigma</td><td>Z741333-1EA</td><td>仅需 1 台</td></tr><tr><td>实验室喷灯</td><td>Sterlitech</td><td>177200-00</td><td>仅需 1 台</td></tr><tr><td>层流罩</td><td>AirClean Systems</td><td>AC624LFUV</td><td>仅需 1 台</td></tr><tr><td>文库纯化试剂盒</td><td>Qiagen</td><td>28704</td><td>每套试剂盒可进行 50 次反应</td></tr><tr><td>磁力架</td><td>Alpaqua</td><td>A001219</td><td>仅需 1 个</td></tr><tr><td>微量离心机</td><td>Thermo Fisher Scientific</td><td>75004061</td><td>仅需 1 台</td></tr><tr><td>移液器 (1000 &amp;micro;L 量程)</td><td>Pipette.com</td><td>L-1000</td><td>仅需 1 支</td></tr><tr><td>移液器 (2 &amp;micro;L 量程)</td><td>Pipette.com</td><td>L-2</td><td>仅需 1 支</td></tr><tr><td>移液器 (20 &amp;micro;L 量程)</td><td>Pipette.com</td><td>L-20</td><td>仅需 1 支</td></tr><tr><td>移液器 (200 &amp;micro;L 量程)</td><td>Pipette.com</td><td>L-200R</td><td>仅需 1 支</td></tr><tr><td>NEBNext Ultra II RNA 文库构建试剂盒（含样本纯化磁珠）</td><td>New England BioLabs Inc.</td><td>E7775S</td><td>每套试剂盒包含 24 个样本所需的试剂。</td></tr><tr><td>Parafilm 封口膜</td><td>MilliporeSigma</td><td>P7793-1EA</td><td>每个滤膜需 2 块 1&quot; × 1&quot; 的正方形封口膜</td></tr><tr><td>PCR 管</td><td>Thermo Fisher Scientific</td><td>AM12230</td><td>每轮反应需 1 支</td></tr><tr><td>移液器吸头 (适用于 1000 &amp;micro;L 量程)</td><td>Pipette.com</td><td>LF-1000</td><td>每包 576 个吸头</td></tr><tr><td>移液器吸头 (适用于 20 &amp;micro;L 量程)</td><td>Pipette.com</td><td>LF-20</td><td>每包 960 个吸头</td></tr><tr><td>移液器吸头 (适用于 200 &amp;micro;L 量程)</td><td>Pipette.com</td><td>LF-250</td><td>每包 960 个吸头</td></tr><tr><td>PowerWulf ZXR1+ 计算机集群</td><td>PSSC Labs</td><td>无编号</td><td>此为一台足够强大的超级计算机示例，可及时完成宏转录组学分析。仅需 1 台。&amp;nbsp;</td></tr><tr><td>Qubit 荧光计启动试剂盒</td><td>Thermo Fisher Scientific</td><td>Q33239</td><td>包含 Qubit 4 荧光计、100 次 DNA 检测所需的试剂以及 500 个 Qubit 管</td></tr><tr><td>药勺</td><td>Thermo Fisher Scientific</td><td>14-357Q</td><td>仅需 1 个</td></tr><tr><td>无菌刀片</td><td>AD Surgical</td><td>A600-P10-0</td><td>每个滤膜需 1 片</td></tr><tr><td>Sterivex-GP 压力过滤装置</td><td>MilliporeSigma</td><td>SVGP01050</td><td>每个水样需 1 个滤器</td></tr><tr><td>PCR 仪</td><td>Bio-Rad</td><td>1861096</td><td>仅需 1 台</td></tr><tr><td>老虎钳</td><td>Irwin</td><td>2078500</td><td>仅需 1 个（用于打开滤器）</td></tr><tr><td>Vortex-Genie 2 涡旋振荡器</td><td>MilliporeSigma</td><td>Z258415-1EA</td><td>仅需 1 台</td></tr><tr><td>WHIRL-PAK 袋</td><td>ULINE</td><td>S-22729</td><td>每个滤膜需 1 个</td></tr><tr><td>ZymoBIOMICS DNA/RNA Miniprep 试剂盒</td><td>Zymo Research</td><td>R2002</td><td>每套试剂盒包含 50 个样本所需的试剂。&amp;nbsp;</td></tr></tbody></table>

利用微生物分子特征评估水力压裂对溪流的影响

Evaluating the Impact of Hydraulic Fracturing on Streams using Microbial Molecular Signatures

水力压裂（HF），通常称为&#34;压裂&#34;，是利用高压水、砂和化学物质的混合物来压裂岩石，从而释放石油和天然气。这一技术彻底改变了美国的能源产业，使其能够获取以往无法利用的资源，目前美国三分之二的天然气产量均来自该技术。尽管压裂对美国经济产生了积极影响，但多项研究已指出其对环境造成的有害影响。尤其令人关注的是压裂对源头溪流的影响，因为这些溪流对整个流域的健康具有不成比例的重要作用。这些溪流中的细菌可作为溪流健康状况的指示指标，因为在受干扰溪流中存在的细菌种类及其丰度，预期会与条件相似但未受干扰的溪流有所不同。因此，本实验方案旨在利用细菌群落来判断溪流是否受到压裂活动的影响。为此，必须采集靠近压裂作业区域（可能受影响）以及位于压裂活动上游或不同流域（未受影响）的溪流中的沉积物和水样。随后对这些样品进行核酸提取、文库构建和测序，以分析微生物群落组成。接着可采用相关性分析和机器学习模型，识别出能够解释群落差异的关键特征，并鉴定可用于预测压裂影响的生物标志物。这些方法能够揭示不同源头溪流中微生物群落在压裂距离远近上的多种差异，并为今后进一步研究压裂活动的环境影响奠定基础。

DNA 和 RNA 提取的成功与否可以通过多种仪器和方案进行评估。通常，只要能够检测到任意浓度的 DNA 或 RNA，即可认为提取成功。因此，查看表 1 可知，除一个样本外，其余所有提取均被视为成功。此步骤失败的原因通常包括初始生物量过低、样品保存不当或提取过程中的人为操作失误。对于滤膜样品而言，即使浓度低于检测限，提取仍可能实际成功。如果这些提取物在 PCR 扩增（如进行 16S 扩增）中未产生条带，或在文库构建后（宏转录组学）仍无法检测到浓度，则很可能确实提取失败。
如果遵循16S实验方案，PCR扩增后出现明亮的条带（如图1中第4和第6孔所示）表明实验成功，而未出现条带（如第一行其他孔所示）则表明实验失败。此外，若含有阴性PCR对照的凝胶泳道中出现明亮条带，也表明实验失败，因为此时很难排除污染对阴性对照的影响未波及样本。
对于16S和宏转录组测序，均可通过获得的序列数量来评估测序结果的成功与否（图2）。16S样品的序列数应至少达到1,000条，理想情况下应不少于5,000条（图2A）。同样，宏转录组样品的序列数应至少达到500,000条，理想情况下应不少于2,000,000条（图2B）。序列数低于最低要求的样品不应用于后续分析，因其可能无法准确代表其细菌群落结构。然而，序列数介于最低值与理想值之间的样品仍可使用，但若大量样品均处于该范围，则结果解释需更加谨慎。
后续下游分析的成功与否，可简单根据是否获得预期的输出文件来判断。无论如何，QIIME2 和 R 等程序应能够基于水力压裂活动对细菌群落之间可能存在的显著差异进行评估（图3）。图3 的数据来自在13条不同溪流的21个不同地点采集的沉积物样品，用于16S和宏转录组学分析。在这21个采样点中，有12个位于水力压裂活动的下游，归类为HF+；另外9个位于水力压裂活动的上游，或处于未进行水力压裂的流域，归类为HF-。除水力压裂活动的存在与否外，这些溪流在其他方面具有可比性。
这些差异可能表现为基于水力压裂状态的群落组成持续性变化。如果情况如此，HF+ 和 HF- 样本在 PCoA 图中应彼此分开聚类，如图 3A 和 图 3B 所示。 为了确认这些明显的差异并非排序方法的人为假象，还需进行进一步的统计分析。例如，对生成图 3A 和 图 3B 所依据的距离矩阵进行 PERMANOVA22 检验，结果显示基于水力压裂状态存在显著聚类，表明图中观察到的分离现象反映了样本间细菌群落的真实差异，而非排序方法所致的人为假象。显著的 PERMANOVA 或 ANOSIM 结果强烈提示 HF+ 与 HF- 样本之间存在一致性的差异，这表明 HF+ 样本受到了水力压裂的影响；而较高的 p 值则表明样本未受水力压裂影响。宏转录组数据同样可以使用相同的方法进行可视化和评估。
分析差异性特征（微生物或功能）可以揭示样本可能受到干扰的证据。确定差异性特征的一种方法是构建随机森林模型。该随机森林模型可用于评估样本的水力压裂状态能否被准确分类。如果模型的分类效果优于随机预期，则进一步表明样本间的差异与水力压裂状态相关。此外，模型中最重要的预测因子将揭示哪些特征在正确区分样本时最为关键（图3C）。这些特征的数值也会根据水力压裂状态表现出持续性的差异。一旦确定了这些差异性特征，便可查阅文献，判断它们是否曾被报道与水力压裂相关。然而，寻找涉及差异性功能的研究可能具有挑战性，因为大多数研究仅使用了16S rRNA组成数据。因此，在评估差异性功能的意义时，一种可行的方法是判断这些功能是否曾被报道与压裂液中常用生物杀灭剂的潜在抗性相关，或是否可能有助于耐受高盐环境。此外，分析特定分类单元的功能谱也可为水力压裂的影响提供证据（图3D）。例如，若某一分类单元被随机森林模型识别为差异性类群，则可将其在HF+样本中的抗微生物耐药谱与HF-样本中的谱型进行比较；若两者存在显著差异，则提示含有生物杀灭剂的压裂液可能已进入溪流。
<table border="1" fo:keep-together.within-page="1" fo:keep-with-next.within-page="always">
	<tbody>
		<tr>
			<td>样本编号</td>
			<td>浓度 (ng/µL)</td>
		</tr>
		<tr>
			<td>1</td>
			<td>1.5</td>
		</tr>
		<tr>
			<td>2</td>
			<td>1.55</td>
		</tr>
		<tr>
			<td>3</td>
			<td>0.745</td>
		</tr>
		<tr>
			<td>4</td>
			<td>0.805</td>
		</tr>
		<tr>
			<td>5</td>
			<td>7.82</td>
		</tr>
		<tr>
			<td>6</td>
			<td>0.053</td>
		</tr>
		<tr>
			<td>7</td>
			<td>0.248</td>
		</tr>
		<tr>
			<td>8</td>
			<td>0.945</td>
		</tr>
		<tr>
			<td>9</td>
			<td>1.82</td>
		</tr>
		<tr>
			<td>10</td>
			<td>0.804</td>
		</tr>
		<tr>
			<td>11</td>
			<td>0.551</td>
		</tr>
		<tr>
			<td>12</td>
			<td>1.69</td>
		</tr>
		<tr>
			<td>13</td>
			<td>4.08</td>
		</tr>
		<tr>
			<td>14</td>
			<td>Below_Detection</td>
		</tr>
		<tr>
			<td>15</td>
			<td>7.87</td>
		</tr>
		<tr>
			<td>16</td>
			<td>0.346</td>
		</tr>
		<tr>
			<td>17</td>
			<td>2.64</td>
		</tr>
		<tr>
			<td>18</td>
			<td>1.15</td>
		</tr>
		<tr>
			<td>19</td>
			<td>0.951</td>
		</tr>
	</tbody>
</table>
表1：基于Fluorometer 1x DS DNA高灵敏度检测法测定的DNA浓度示例. 除样品14外，其余所有样品的提取均被视为成功，因其含有可检测到的DNA量。
<img alt="显示DNA分离的琼脂糖凝胶电泳结果；可见分子量标准。" class="xfigimg" src="/files/ftp_upload/61904/61904fig01.jpg" data-altupdatedat="1788304762000" data-alt="Figure 1"> 
图1：PCR产物的电泳凝胶示例。 凝胶预先染色，并在紫外光下观察，使凝胶上存在的任何DNA发出荧光。第一行第4和第6孔的样品PCR成功，因为二者均出现一条亮度高且大小符合预期的单一条带（根据标准分子量标记判断）。其余六个孔的样品PCR失败，因为未产生任何条带。阳性对照（第二行第一孔）出现明亮条带，表明PCR操作正确；阴性对照（第二行第6和第7孔）无条带，表明样品未受污染。如果阴性对照出现与样品亮度相当的条带，则认为PCR失败，因为无法排除样品中的扩增产物仅来源于污染的可能性。<a href="https://www.jove.com/files/ftp_upload/61904/61904fig01large.jpg" target="_blank">请点击此处查看该图的放大版本。</a>
<img alt="测序数据分布；密度图；比例与序列；比较分析图" class="xfigimg" src="/files/ftp_upload/61904/61904fig02.jpg" data-altupdatedat="1788304762000" data-alt="Figure 2"> 
图2：序列数示例。（A）16S序列数示例。几乎所有这些16S样本的序列数均超过1,000。极少数序列数少于1,000的样本应从后续分析中剔除，因其序列数量不足以准确代表其细菌群落。部分样本序列数介于1,000至5,000之间；虽然并非理想，但由于超过了最低阈值，仍可使用，且大多数样本也超过了5,000这一理想最低值。（B）宏转录组学序列数示例。所有样本的序列数均超过了最低值（500,000）和理想最低值（2,000,000）。因此，所有样本的测序均成功，均可用于后续分析。<a href="https://www.jove.com/files/ftp_upload/61904/61904fig02large.jpg" target="_blank">请点击此处查看该图的放大版本。</a>
<img alt="水力压裂研究中变量的PCA与统计图表对比；包含饼图和散点图。" class="xfigimg" src="/files/ftp_upload/61904/61904fig03.jpg" data-altupdatedat="1788304762000" data-alt="图3"> 
图3：示例分析。（A）基于通过QIIME2创建并可视化的加权Unifrac距离矩阵计算出的坐标生成的PCoA图。（B）基于从QIIME2导出的加权Unifrac距离矩阵计算出的坐标生成的PCoA图。该坐标的可视化使用R语言中的Phyloseq和ggplot2软件包完成。利用Vegan软件包将元数据向量拟合到图中。每个点代表一个样本的细菌群落，点之间的距离越近，表明群落组成越相似。在这些16S沉积物样本中观察到基于水力压裂状态的聚类现象（PERMANOVA，p=0.001）。此外，向量显示HF+样本通常具有较高水平的钡（Barium）、溴化物（Bromide）、镍（Nickel）和锌（Zinc），其细菌群落组成与HF-样本不同。（C）随机森林模型中最佳预测因子的图表，用于检验细菌丰度是否可用于预测样本的水力压裂状态。该随机森林模型通过R语言中的randomForest软件包构建。图中展示了前20个预测因子，以及使用这些因子分离样本时导致的杂质减少程度（即HF+和HF-样本被错误归为同一组的数量度量），以平均Gini指数下降值（Mean Decrease in Gini Index）表示。（D）基于宏转录组数据的伯克氏菌目（Burkholderiales）抗微生物耐药性谱的饼图。首先使用Kraken2对序列进行注释，以确定其所属分类单元；随后将这些已注释的序列与MEGARes 2.0数据库进行BLAST比对，以鉴定正在活跃表达的抗微生物耐药基因（以“MEG_#”形式表示）。然后提取伯克氏菌目成员所表达的抗微生物耐药基因，以确定在该分类群中最常见的耐药基因类型。尽管宏转录组学成本更高且耗时更长，但它能够实现功能层面的分析，例如本例中的分析，而这是16S数据无法实现的。值得注意的是，本示例分析使用了Kraken2而非HUMAnN2。Kraken2比HUMAnN2更快，但它仅输出组成信息，而无法像HUMAnN2那样同时提供组成、贡献度、基因功能及代谢通路信息。<a href="https://www.jove.com/files/ftp_upload/61904/61904fig03large.jpg" target="_blank">请点击此处查看该图的放大版本。</a>
补充文件：一个宏转录组学分析流程的示例。 <a href="https://www.jove.com/files/ftp_upload/61904/Supplemental%20File_61904_R2_RE.docx" target="_blank">请点击此处下载该文件。</a>

本文介绍了一种通过分析水体和沉积物中的微生物群落来研究水力压裂对附近溪流影响的实验方案。

水力压裂（hydraulic fracturing，简称 HF），通常称为“压裂”（fracking），利用高压水、砂和化学物质的混合物来压裂岩石，从而释放出石油……

本方案可用于回答水力压裂是否以及如何影响附近溪流中的细菌，进而影响溪流本身的问题。该技术的主要优势在于其整体性，因为它将引导研究人员完成从样品采集到数据分析的全部流程。演示该实验流程的将是本实验室的生物信息学家Jeremy Chansee、朱尼亚塔学院的本科生Sydney Regal，以及本实验室的实验室主管Gillian Leister。为进行核酸提取而采集沉积物样品时，应佩戴手套，将一个带盖的无菌50毫升锥形管从岸边浸入溪流水中。在锥形管浸没于水中的状态下，打开管盖，并用管盖从水下1至3厘米深度处舀取约3毫升沉积物放入管中。采集样品后，将管中除约1毫升以外的水倒出，然后使用1,000微升移液器向样品中加入3毫升DNA/RNA保存液。将带盖的锥形管涡旋混匀5秒，充分混合保存液与样品，并将样品置于冰上保存。返回实验室后，将样品在-20°C条件下保存用于16S DNA分析，或在-70°C条件下保存用于宏转录组RNA分析。对于滤膜采集样品，先用溪水完全充满并排空一个无菌1升瓶三次，以润洗容器，最后再将瓶子完全充满一次。在稳定的表面上，用无菌鲁尔锁注射器吸取足量的溪水，并将该注射器连接至一个无菌且无DNA/RNA污染、直径为1.7厘米、孔径为0.22微米的聚醚砜滤膜。将溪水全部体积通过滤膜进行过滤。当整个样品体积完成同样过滤后，向注射器中吸入约20毫升空气，并推动空气通过滤膜，以去除滤膜上残留的多余水分。接下来，使用 P1000 微量移液器吸取 2 毫升 DNA/RNA 防腐剂，在水平持握滤器的同时，将移液器枪头置于滤器筒腔内，对准滤器较大开口处加样，确保防腐剂进入滤器内部。然后用石蜡膜将滤器的每个开口紧密包裹密封，并将滤器放入盛有冰块的无菌样品袋中。采样结束后，按照所示方法将滤器储存，用于后续 16S 或宏转录组分析。开始样品转移前，使用10%次氯酸钠和70%乙醇清洁工作区域。从沉积物样品中提取核酸时，使用经火焰和乙醇灭菌的金属工具，将约250毫克样品转移至微量离心管中。从滤膜样品中提取核酸时，使用经70%乙醇和火焰灭菌的老虎钳，在无菌表面上打开滤膜外壳，并将滤芯从外壳中取出。使用无菌手术刀在核心样本的顶部、底部以及接缝处进行切割，并使用无菌镊子将滤纸折叠后，再用手术刀将其切成小块。然后小心地将滤纸碎片放入微量离心管中，避免接触任何未灭菌的表面。为了裂解上述任一样本中的细胞，将离心管置于细胞破碎仪中高速处理。至少静置五分钟后，离心样品并将上清液转移至新的无菌微量离心管中。按1:1比例向上清液中加入裂解缓冲液，然后将溶液转移至提供的滤膜中。将滤膜放入新的微量离心管中，并向管中加入400 µL的预处理缓冲液。离心后，向离心管中加入 700 微升洗涤缓冲液，再次离心滤器。弃去穿流液后，向管中加入 400 微升洗涤缓冲液进行额外一次离心，并将滤器转移至新的无菌微量离心管中。为洗脱 DNA，在室温下用 50 微升无 DNase/RNase 的水处理滤器 5 分钟。同时，将一个 HRC 滤膜放入收集管中，并加入 600 µL HRC 提取液。离心后，将滤膜转移至新的无菌微量离心管中，并将洗脱的 DNA 加入滤膜中进行离心。流出液即为提取的 DNA。为了构建DNA 16S RNA文库，首先使用新提取的DNA产物，采用标准PCR方案对16S核糖体RNA进行扩增。取7微升PCR产物与13微升无DNase的水混合，将PCR溶液上样至2%琼脂糖凝胶。然后以90伏电压电泳60至90分钟，检测是否存在大小为386个碱基对的条带，以确认扩增是否成功。为了纯化DNA 16S核糖体RNA文库，将PCR产物中产生明亮条带的10微升和产生微弱条带的13微升样品汇集至一个无菌微量离心管中，并将每份混合样品约150至200纳克上样至新的2%琼脂糖凝胶的各个加样孔中。按照示范条件电泳后，从凝胶中切下386碱基对的条带，使用商业化试剂盒纯化DNA。随后用30微升10毫摩尔Tris盐酸溶液洗脱纯化的DNA，并将纯化后的文库置于干冰中保存，随后运送进行高通量测序。在此代表性分析中，除一个样本外，其余所有提取均被视为成功。PCR 扩增后观察到的明亮条带表明 16S 方案成功。16S 样本应至少包含 1,000 条序列，理想情况下应达到至少 5,000 条序列；而宏转录组样本应至少包含 500,000 条序列，理想情况下应达到至少 200 万条序列。展示了来自13条不同溪流的21个不同采样点的沉积物样本，用于16S和宏转录组学分析。在这些采样点中，有12个位于水力压裂活动的下游，被归类为水力压裂阳性；另外9个位于水力压裂活动的上游，或位于无水力压裂活动的流域，被归类为水力压裂阴性。通过PERMANOVA分析评估，水力压裂阳性和阴性样本数据之间的显著分离表明，水力压裂阳性样本受到了水力压裂活动的影响。最重要的随机森林预测因子将揭示哪些特征对于正确区分样本最为关键。如果随机森林模型识别出某一分类单元具有重要性，则可将其在水力压裂阳性样本中的抗微生物耐药性谱型与水力压裂阴性样本中的谱型进行比较。若两者存在显著差异，则可能表明含有杀菌剂的压裂液已进入溪流。微生物无处不在，因此这类工作存在显著的污染风险。最初的样品转移步骤尤其容易受到污染。如果研究者关注微生物的生物降解或其他代谢过程，可采用RNA宏基因组测序来研究微生物活跃的基因表达。该技术可实现对原位细菌群落研究的分子技术标准化，以及对细菌序列数据的生物信息学分析。

evaluating-impact-hydraulic-fracturing-on-streams-using-microbial

Research

科研

JoVE Journal

Environment

环境

使用微生物分子特征评估水力压裂对溪流的影响

本文描述了基于微孔板的胞外酶活性比色法或荧光法分析流程。这些方法可在合理的时间范围内，对大量环境样品中的酶活性进行快速测定。...

利用高通量微孔板法测定水体、土壤和沉积物中微生物胞外酶活性

表征微生物群落一直是环境微生物学中的长期目标。新一代测序技术如今能够以极低的成本和工作量，以前所未有的深度对微生物群落进行表征。本文详细介绍了我们使用台式测序仪对细菌16S核糖体RNA基因进行测序的方法。...

16S核糖体RNA基因扩增子的高通量测序

本文介绍了一种全面的实验方案，用于评估偏远海洋环境中有机和无机营养物质的可利用性，以及微生物和病毒群落的丰度与结构。...

揭开海洋中看不见的参与者——水化学与海洋微生物学野外指南

本论文描述了一种双重数字PCR检测方法，可用于同时定量Enterococcus...

一种用于同时定量水中粪肠球菌属（Enterococcus spp.）和人源粪便相关标记物HF183的双重数字PCR检测方法

本研究介绍了一种用于研究土壤渗漏计中地下水分、地球化学及微生物异质性的挖掘方法。该渗漏计模拟了一个最初处于均质条件的人工坡面，在18个月期间经过八次灌溉，累计接受了约5,000...

用于水文、地球化学与微生物学耦合研究的土壤溶滤器采样

流式细胞术结合可视化聚类分析为研究水生生物膜提供了一种简便快速的方法。该方法可用于生物膜的表征、生物膜群落结构变化的检测，以及嵌入生物膜中非生物颗粒的检测。...

利用流式细胞术表征水生生物膜

我们介绍了从水生环境中分离、扩增和表征降解烃类细菌的过程。该方案概述了细菌的分离方法、基于16S...

从水生环境分离、扩增与鉴定具有碳氢化合物代谢特性的细菌物种

该方案可检测美国南部德克萨斯州沿海湿地中与甲烷循环相关的关键基因，并可视化其空间分布，从而加深对这些动态生态系统中甲烷调控机制及其环境影响的理解。...

可视化滨海湿地中甲烷循环微生物的动态变化

许多重要的活性矿物在过程中形成 原位...

矿物采样器的部署与回收

可靠且准确的现场可部署DNA分析方法对许多行业至关重要。本文提供了两种用于野外环境微生物分析的方法：1）利用生物提取野外试剂盒从环境样本中分离微生物DNA；2）通过野外qPCR仪进行DNA分析。...

使用生物提取野外试剂盒和野外qPCR仪进行现场微生物DNA分析

了解利用微生物分子特征评估水力压裂对溪流生态系统影响的方法。本方案涵盖从样本采集、核酸提取到16S rRNA测序的完整流程，为环境微生物学研究提供标准化技术支持。

利用微生物分子特征评估水力压裂对溪流的影响

探索更多视频