# Sequencing Depth Has a Stronger Effect than DNA Extraction on Soil Bacterial Richness Discovery

Concepcion Sanchez-Cid, Romie Tignat-Perrier, Franqueville Laure, Laurence

Delaurière, Trista Schagat, Timothy M. Vogel

► To cite this version:

Concepcion Sanchez-Cid, Romie Tignat-Perrier, Franqueville Laure, Laurence Delaurière, Trista Schagat, et al.. Sequencing Depth Has a Stronger Effect than DNA Extraction on Soil Bacterial Richness Discovery. Biomolecules, 2022, 10.3390/biom12030364. [hal-03662087](https://hal.science/hal-03662087)

HAL Id: hal-03662087

https://hal.science/hal-03662087

Submitted on 8 Jun 2022

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
Article

# Sequencing Depth Has a Stronger Effect than DNA Extraction on Soil Bacterial Richness Discovery

Concepcion Sanchez-Cid <fnref n="1" />,<fnref n="2" />,* , Romie Tignat-Perrier <fnref n="1" />,<fnref n="3" /> , Laure Franqueville <fnref n="1" /> , Laurence Delaurière <fnref n="2" />, Trista Schagat <fnref n="4" /> and Timothy M. Vogel <fnref n="1" />

**Abstract:** Although Next-Generation Sequencing techniques have increased our access to the soil microbiome, each step of soil metagenomics presents inherent biases that prevent the accurate definition of the soil microbiome and its ecosystem function. In this study, we compared the effects of DNA extraction and sequencing depth on bacterial richness discovery from two soil samples. Four DNA extraction methods were used, and sequencing duplicates were generated for each DNA sample. The V3–V4 region of the 16S rRNA gene was sequenced to determine the taxonomical richness measured by each method at the amplicon sequence variant (ASV) level. Both the overall functional richness and antibiotic resistance gene (ARG) richness were evaluated by metagenomics sequencing. Despite variable DNA extraction methods, sequencing depth had a greater influence on bacterial richness discovery at both the taxonomical and functional levels. Sequencing duplicates from the same sample provided access to different portions of bacterial richness, and this was related to differences in the sequencing depth. Thus, the sequencing depth introduced biases in the comparison of DNA extraction methods. An optimisation of the soil metagenomics workflow is needed in order to sequence at a sufficient and equal depth. This would improve the accuracy of metagenomic comparisons and soil microbiome profiles.

**Keywords:** sequencing depth; DNA extraction; soil microbiome; bacterial richness; metagenomics

<figure><img src="image_2.png" /></figure>

check for
updates

**Citation:** Sanchez-Cid, C.; Tignat-Perrier, R.; Franqueville, L.; Delaurière, L.; Schagat, T.; Vogel, T.M. Sequencing Depth Has a Stronger Effect than DNA Extraction on Soil Bacterial Richness Discovery.
*Biomolecules* **2022**, *12*, 364. https://doi.org/10.3390/biom12030364

Academic Editors: Michael Van Dyke and Alexander A. Bolshoy

Received: 11 January 2022

Accepted: 20 February 2022

Published: 25 February 2022

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_3.png" /></figure>

**Copyright:** © 2022 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

## 1. Introduction

The soil ecosystem arguably harbours the highest diversity of microorganisms of any ecosystem [1]. Unravelling the composition and function of the soil microbiome is critical to better understanding the role of these microbial communities in soil function and ecosystem services. The use of Next Generation-Sequencing (NGS) techniques has increased our access to the microbial communities present in soil, especially to the large proportion of uncultured microorganisms [2,3]. Nevertheless, each methodological step from soil sampling to sequence annotation presents inherent biases that limit the depth and reliability of soil microbiome analyses [4–7]. Of all these biases, the ones associated with DNA extraction have been particularly highlighted for their effects [8,9]. DNA can be adsorbed by soil compounds such as clay [10–12], which, when combined with the presence of lysis-recalcitrant bacteria [13], reduces DNA extraction efficiency. Moreover, organic matter and humic acids that are known to potentially inhibit enzymatic reactions [14,15] are often coextracted with DNA. Many studies have compared DNA extraction methods

<footnote n="1">¹ Environmental Microbial Genomics, Laboratoire Ampère, CNRS UMR 5005, Ecole Centrale de Lyon, Université de Lyon, 69134 Ecully, France; rom26.p@hotmail.fr (R.T.-P.); laure.franqueville@ec-lyon.fr (L.F.); vogel@univ-lyon1.fr (T.M.V.)</footnote>

<footnote n="2">² Scientific Applications and Training, Promega France, 69100 Charbonnières-les-Bains, France; laurence.delauriere@promega.com</footnote>

<footnote n="3">³ Institut des Géosciences de l'Environnement, Université Grenoble Alpes, CNRS, IRD, Grenoble INP, Grenoble, 38400 Saint-Martin-d'Hères, France</footnote>

<footnote n="4">⁴ Scientific Applications and Training, Promega Corporation, Madison, WI 53711, USA; trista.schagat@promega.com</footnote>

<footnote>* Correspondence: concepcion.sanchezcid-torres@ec-lyon.fr</footnote>
and documented the biases imposed by lysis procedures [16–19]. This has provoked
proposals of different methods for DNA extraction and purification from soil over the past
few decades [16,20–22] in an attempt to obtain an unbiased picture of soil microbiome
biodiversity [10]. Nevertheless, no method has been shown to overcome all the biases
described above, and the debate on the choice of DNA extraction methods is still ongoing.

The criteria used to define the performance of a DNA extraction method vary between studies and range from nucleic acid yield to phylogenetic diversity. However, higher yields, purity, and integrity do not always imply an improvement in bacterial diversity discovery [5]. Analyses of the relative abundance of taxonomic groups have a limited potential for selecting nucleic acid extraction methods, since the biases associated with soil metagenomics prevent us from determining the actual distribution of soil microbial populations within a community [23]. Furthermore, nucleic acid extraction methods may modify the relative abundance of detected communities without affecting the bacterial richness discovery [24]. In other words, the use of different DNA extraction methods might detect different proportions of the same communities [25] rather than different taxa or functions. The goal of this paper was not to resolve the debate concerning the use of relative abundance versus richness. We consider that bacterial richness measurements provide a more objective comparison of the performance of nucleic acid extraction methods, since biodiversity calculations with a limited set of sequences are strongly biased by evenness, which depends on the number of sequences. In addition, the actual relative abundances of the soil microbiome remain unknown.

The sequencing depth has also shown an important effect on bacterial discovery. Several studies have observed that low sequencing depths may bias the evaluation of the composition and function of microbial communities, both in terms of richness and diversity [26–28]. However, to the best of our knowledge, this is the first study to compare simultaneously the relative contribution of both sequencing depth and DNA extraction to soil bacterial richness discovery. DNA was extracted from two soil samples using four DNA extraction methods: two novel semiautomated methods for DNA extraction, a commercial kit, and the phenol/chloroform method as described by Griffiths et al. [29]. Sequencing duplicates were generated for each DNA sample to account for the effects of sequencing depth on bacterial richness discovery, which was evaluated at three levels. First, taxonomic richness obtained from the sequencing of the V3–V4 hypervariable region of the 16S rRNA gene was measured at the amplicon sequence variant (ASV) level. Second, the overall functional richness was determined from metagenomic samples annotated using the second level of classification of the SEED hierarchical subsystems [30]. Finally, the richness discovery of antibiotic resistance genes (ARGs) annotated from metagenomic reads using the CARD database [31] was also assessed. These genes are not abundant and are often associated with Mobile Genetic Elements (MGE) [32–34], which affects their identification using metagenomics and their association to bacterial hosts [35]. Given ARG's relatively low abundance in the environment, the sequencing depth could have a stronger impact on their richness discovery than on the overall functional richness discovery.

**2. Materials and Methods**

2.1. *Soil Sampling*

Two soils were selected for this study (Table 1). Soils were sampled at an experimental farm (Scottish Agricultural College, Craibstone, Scotland, grid reference NJ872104) and at a field planted with corn at La Côte Saint Andre, France. All samples were kept at 4 °C before DNA extraction. Details about the Scottish Agricultural College soil composition were provided by Kemp et al. [36]. Physical characterisation of the soil from La Côte de Saint André was performed by CESAR (Centre Scientifique Agricole Régional) using the standard methods (NFX 31-107, NFX 31-117, ISO 10694, and ISO 13878). According to the national legislation in both Scotland and France, where the samples were collected, no ethical application was needed for the collection of soil samples.
<table><caption>**Table 1.** Physical characterisation of the samples selected for this study.</caption><thead><tr><th></th><th>Scottish Agricultural College</th><th>La Côte de Saint André</th></tr></thead><tbody><tr><td>Sand</td><td>73.85%</td><td>42.9%</td></tr><tr><td>Silt</td><td>20.04%</td><td>43.6%</td></tr><tr><td>Clay</td><td>6.11%</td><td>13.5%</td></tr><tr><td>pH</td><td>4.5</td><td>7.24</td></tr><tr><td>Organic matter</td><td>5.97%</td><td>2.92%</td></tr><tr><td>Organic C</td><td>3.79%</td><td>1.7%</td></tr><tr><td>Total N</td><td>0.45%</td><td>0.17%</td></tr></tbody></table>

## 2.2. DNA Extraction and Purification

DNA was extracted from 250 mg of sample using the DNeasy PowerSoil Kit (QIAGEN, Hilden, Germany) and the phenol/chloroform extraction method described by Griffiths et al. [29], as well as a new semiautomated protocol in which the Maxwell RSC Instrument (Promega, Madison, WI, USA) and a prototype version of the Maxwell RSC Fecal Microbiome DNA Kit (Promega) are used for DNA purification. Two modifications of this protocol were tested and are referred to as the Maxwell 1 and Maxwell 2 methods. All DNA extractions were performed in triplicate. In the Maxwell 1 method, 250 mg of sample were diluted in 1 mL of Lysis Buffer (Promega) and heated for 5 min at 95 °C. Samples underwent bead-beating twice at 5.5 m/s for 30 s in Lysis Matrix E tubes (MP Biomedicals, Irvine, CA, USA) and centrifuged at 10,600 × g for 5 min. Then, 300 µL of supernatant were added to 300 µL of Binding Buffer (Promega) and loaded into a Maxwell RSC cartridge containing magnetic beads for DNA purification on the Maxwell RSC Instrument, according to Technical Manual TM640. A second purification using the ProNex Size-Selective Purification System (Promega) was carried out to reduce the humic acid carryover. In the Maxwell 2 method, two variants were introduced in the previously described protocol: 500 µL of Lysis Buffer were mixed with 500 µL of 0.5-M Sodium Phosphate Buffer (pH 7.0) and added to 250 mg of sample, and cells were lysed without bead-beating.

## 2.3. Quantitative PCR (qPCR) Assay

DNA was quantified using the Qubit Fluorometer and the Qubit dsDNA HS Assay Kits (Thermo Fisher Scientific, Waltham, MA, USA) and diluted to 2.5 ng/µL. Total bacterial abundance was estimated by quantifying the V3 region of the 16S rRNA gene by qPCR using the "universal" primers 341F (5'-CCT ACG GGA GGC AGC AG-3') and 534R (5'-ATT ACC GCG GCT GCT GGC A-3') [37,38]. A qPCR assay was carried out using the Corbett Rotor-Gene 6000 (QIAGEN) in a 20-µL reaction volume containing GoTaq qPCR Master Mix (Promega), 0.75 µM of each primer, and 2 µL of DNA at 2.5 ng/µL (5 ng of DNA). Two non-template controls were also included. The standard curve was obtained using 10-fold serial dilutions (10⁷ to 10² copies) of a linearised plasmid pGEM-T Easy Vector (Promega) containing the 16S rRNA gene of *Pseudomonas aeruginosa* PAO1. Cycling conditions for qPCR amplification were 95 °C for 2 min, followed by 30 cycles of 95 °C for 15 s, 60 °C for 30 s, and 72 °C for 30 s. Melting curves were generated after amplification by increasing the temperature from 60 °C to 95 °C.

## 2.4. 16S rRNA Gene Amplicon Sequencing and Analysis

All DNA samples were diluted to 2.5 ng/µL before 16S rRNA gene PCR amplification and sequencing to reduce PCR inhibition and normalise the amount of starting material for library preparation. Then, the V3–V4 hypervariable region of the bacterial 16S rRNA gene was amplified using Titanium Taq DNA Polymerase (Takara Clontech, Kyoto, Japan) and forward 341F with Illumina overhang (5'-TCG TCG GCA GCG TCA GAT GTG TAT AAG AGA CAG TCG TCG GCA GCG TCA GAT GTG TAT AAG AGA CAG CCT ACG GGN GGC WGC AG-3') and reverse 785F with Illumina overhang (5'-GTC TCG TGG GCT CGG AGA TGT GTA TAA GAG ACA GGT CTC GTG GGC TCG GAG ATG TGT ATA AGA GAC AGG ACT ACH VGG GTA TCT AAT CC-3') primers [6] to identify the
Amplicon Sequence Variants (ASVs). Cycling conditions for PCR amplification were 95 °C for 3 min, followed by 25 cycles of 95 °C for 30 s, 55 °C for 30 s, and 72 °C for 30 s and a final extension step at 72 °C for 5 min. DNA libraries were prepared based on Illumina's "16S Metagenomics Library Prep Guide" (15044223 Rev. B) using the Platinum Taq DNA Polymerase (Invitrogen, Waltham, MA, USA) and the Nextera XT Index Kit V2 (Illumina, San Diego, CA, USA). DNA sequencing with a 15% PhiX spike-in was performed using the MiSeq System and the MiSeq Reagent Kit v2 (Illumina). All samples were sequenced twice (in two different runs following the same protocol), except two samples extracted from La Côte de Saint André—one using the Maxwell 1 method and another one using the phenol/chloroform method—for which there was no DNA left for a second sequencing run. Sequences were trimmed, filtered, and denoised, ASVs were inferred, and chimeras were removed using DADA2 [39]. Since the amplicon size was too big to ensure the proper merging of forward and reverse reads, only forward reads were used in this analysis. Nineteen nucleotides and 10 nucleotides were trimmed from the left and right ends of the forward reads, respectively. Reads were filtered to have 0 Ns and a maximum of 2 expected errors. After following all the steps of the DADA2 pipeline, singletons were removed from inferred ASVs. In order to evaluate the effect of the sequencing depth in the taxonomic richness assessment, the ASV richness detected in each sample was determined using the vegan package (version 2.5-6) [40] in R (version 3.5.1) and plotted as a function of the sequencing depth. Finally, Venn diagrams comparing the taxonomic richness discovery of different DNA extraction methods, sequencing runs, individual samples sequenced at different depths, and samples extracted using different methods and sequenced at similar depths (a maximum difference of 1000 sequences between samples) were obtained using the VennDiagram package (version 1.6.20) [41] in R.

## 2.5. Metagenomics Sequencing and Analysis

Metagenomics libraries were prepared from <1 ng of DNA using the Nextera XT Library Prep Kit and Indexes (Illumina), as detailed in Illumina's "Nextera XT DNA Library Prep Kit" reference guide (15031942 v03). DNA sequencing with a 1% PhiX spike-in was performed using the MiSeq System and the MiSeq Reagent Kit v2 (Illumina). All samples were sequenced twice (in two different runs following the same protocol), except two samples extracted from La Côte de Saint André—one using the Maxwell 1 method and another one using the phenol/chloroform method—for which there was no DNA left for a second sequencing run. Metagenomics reads were trimmed using the Fastq Quality Trimmer tool of the FASTX Toolkit. Nucleotides that did not meet a minimum quality score of Q20 were trimmed from the sequences, and sequences shorter than 100 nucleotides after trimming were removed. Two samples extracted using the DNeasy PowerSoil Kit—one from La Côte de Saint André and one from the Scottish Agricultural College—were removed from the analysis, since they did not meet a sequencing depth of 10,000 sequences. Reads were blasted against the nr database using Diamond [42] and filtered at an e-value of 10<sup>-5</sup>, an identity of 60%, and a query cover of 70%. The best hit was selected. Then, the sequences were functionally annotated using MEGAN6 [43] and the SEED hierarchical subsystems [30]. The second level of hierarchical functional subsystems classification was selected for functional class richness discovery analysis, and singletons were removed. In parallel, reads were blasted against the CARD database [31] using Diamond and filtered at an identity of 60%, a length of 33 amino acids, and an e-value of 10<sup>-5</sup>. The best hit was chosen, and singletons were removed. To evaluate the effect of the sequencing depth in both the overall function and ARG richness discovery, the functional class richness and the ARG richness detected in each sample were determined using the vegan package in R and plotted as a function of the sequencing depth. Finally, Venn diagrams comparing the taxonomic richness discovery of different DNA extraction methods, sequencing runs, individual samples sequenced at different depths, and samples extracted using different methods and sequenced at similar depths (a maximum difference of 1000 sequences between samples) were obtained using the VennDiagram package (version 1.6.20) in R.
## 2.6. Statistical Analysis

Correlations between sequencing depth and measured ASV and functional and ARG richness were determined using the Pearson coefficient, and the significance of the correlation was evaluated by a two-tailed Student's *t*-distribution. Samples from each of the two soils were then grouped by a sequencing run or by DNA extraction method to determine the effect of these factors on ASV and functional and ARG richness discovery (Tables S1–S6 in the Supporting Information). Data normality was checked using the Shapiro-Wilk test. Statistical differences in richness between the conditions (sequencing run or DNA extraction method) were evaluated using ANOVA and pairwise two-tailed *t*-Student tests when data showed a normal distribution and using Kruskal-Wallis and pairwise Wilcoxon signed-rank tests when the data showed a non-normal distribution. Only significant differences are shown for pairwise comparisons.

# 3. Results

## 3.1. Total Bacterial Abundance Extracted by Each Method

The total abundance of the bacterial community extracted from two different soils, which was estimated by the number of copies of the 16S rRNA gene, varied between DNA extraction methods. Higher numbers of 16S rRNA gene copies were extracted from both soils using the phenol/chloroform method (Table 2). On the other hand, lower numbers of 16S rRNA gene copies were obtained from the Scottish Agricultural College using the Maxwell 2 modification of the prototype Maxwell RSC Fecal Microbiome DNA Kit (Promega) and from La Côte de Saint André soil using the DNeasy PowerSoil Kit (QIAGEN) compared to using any other method on the same soil.

<table><caption>**Table 2.** 16S rRNA gene copy numbers from the Scottish Agricultural College (Craibstone, Scotland) and La Côte de Saint André (France) soils (250 mg) using different DNA extraction and purification strategies. Maxwell 1 and 2: modifications of the prototype Maxwell RSC Fecal Microbiome DNA Kit (Promega). All data (averages and standard deviations) are based on three separate soil samples for each method.</caption><thead><tr><th>Soil</th><th>Method</th><th>16S rRNA Gene Copies/µL</th></tr></thead><tbody><tr><td rowspan="4">Scottish Agricultural College soil</td><td>Maxwell 1</td><td>153.96 ± 52.41</td></tr><tr><td>Maxwell 2</td><td>64.28 ± 9.13</td></tr><tr><td>Phenol/Chloroform method</td><td>189.21 ± 15.04</td></tr><tr><td>DNeasy PowerSoil Kit</td><td>165.31 ± 65.89</td></tr><tr><td rowspan="4">La Côte de Saint André soil</td><td>Maxwell 1</td><td>160.61 ± 14.7</td></tr><tr><td>Maxwell 2</td><td>170.39 ± 16.32</td></tr><tr><td>Phenol/Chloroform method</td><td>211.95 ± 45.41</td></tr><tr><td>DNeasy PowerSoil Kit</td><td>88.9 ± 25.64</td></tr></tbody></table>

## 3.2. Sequencing Depth and DNA Extraction Effect on Bacterial Richness Discovery

The sequencing depth was significantly correlated to ASV and functional and ARG richness discovery in both soils (Figures 1–3). A greater sequencing depth increased the taxonomic and functional richness discovery from the two soils included in this study, regardless of the method used for DNA extraction or the sequencing run (Figures 1–3). In other words, at a cursory level, all methods detected a similar taxonomical, functional, and ARG richness at equal sequencing depths. Furthermore, DNA extraction triplicates that were sequenced at different depths had access to different proportions of bacterial DNA richness (Figures 1–3), and similar results were observed when different depths were obtained between sequencing runs from individual samples (Figures 1–3). On the other hand, at similar sequencing depths, slight differences in the bacterial richness discovery were sometimes found between samples, regardless of whether they were extracted using
the same or different methods (Figures 1–3). These differences were also found between sequencing duplicates from the same sample sequenced at similar depths: some examples of this are the ASV richness measured from samples 3 and 2 extracted from La Côte de Saint André soil using the Maxwell 1 and the Maxwell 2 modifications of the prototype Maxwell RSC Fecal Microbiome DNA Kit (Promega), respectively (Figure 1B).

<figure><img src="image_4.png" /><figcaption>**Figure 1.** Effect of sequencing depth on the taxonomic richness discovery from (A) the Scottish Agricultural College soil (Craibstone, Scotland) and (B) La Côte de Saint André soil (France). Green: Maxwell 1 modification of the prototype Maxwell RSC Fecal Microbiome DNA Kit (Promega). Red: Maxwell 2 modification of the prototype Maxwell RSC Fecal Microbiome DNA Kit (Promega). Blue: Phenol/Chloroform. Yellow: DNeasy PowerSoil Kit (QIAGEN). The taxonomical richness measured in each sample was determined and plotted as a function of the sequencing depth. For each method/soil pair, triplicates (labelled 1, 2, and 3) were sequenced twice and plotted in the graph. Lighter colours represent the first sequencing run, and darker colours represent the second sequencing run. Correlation *p*-value: 5.15 × 10⁻¹⁰ (Scottish Agricultural College soil) and 7.38 × 10⁻¹⁰ (La Côte de Saint André soil).</figcaption></figure>
<figure><img src="image_5.png" /><figcaption>**Figure 2.** Effect of sequencing depth on the functional richness discovery from (A) the Scottish Agricultural College soil (Craibstone, Scotland) and (B) La Côte de Saint André soil (France). Green: Maxwell 1 modification of the prototype Maxwell RSC Fecal Microbiome DNA Kit (Promega). Red: Maxwell 2 modification of the prototype Maxwell RSC Fecal Microbiome DNA Kit (Promega). Blue: Phenol/Chloroform. Yellow: DNeasy PowerSoil Kit (QIAGEN). The functional richness (second SEED sublevel of classification) measured in each sample was determined and plotted as a function of the sequencing depth. For each method/soil pair, triplicates (labelled 1, 2, and 3) were sequenced twice and plotted in the graph. Lighter colours represent the first sequencing run, and darker colours represent the second sequencing run. Correlation *p*-value: 4.55 × 10<sup>-5</sup> (Scottish Agricultural College soil) and 1.72 × 10<sup>-7</sup> (La Côte de Saint André soil).</figcaption></figure>

The resequencing of the samples in two different runs did not have a significant impact on the ASV and ARG richness discovery, whereas significant differences were found between the functional richness measured from the La Côte de Saint André samples in the two sequencing runs (Tables S1–S3 in the Supporting Information). On the other hand, significant differences were observed in the ASV richness measured by some of the DNA extraction methods in La Côte de Saint André soil (Table S4 in the Supporting Information), as well as in functional and ARG richness in the Scottish Agricultural College soil (Tables S5 and S6 in the Supporting Information). Regarding the ASV richness discovery from La Côte de Saint André, significantly higher measurements were obtained from samples extracted using the DNeasy PowerSoil Kit (Table S4 in the Supporting Information), i.e., the most deeply sequenced method, than from the Maxwell 1 and 2 methods, which showed the lowest sequencing depths (Figure 1B). In addition, the metagenomic sequencing of samples extracted from the Scottish Agricultural College soil using the DNeasy PowerSoil Kit provided significantly higher functional richness than the Maxwell
1 method (Table S5 in the Supporting Information) and significantly higher ARG richness than both the Maxwell 1 and the Maxwell 2 methods (Table S6 in the Supporting Information). DNeasy PowerSoil Kit triplicates from the first metagenomic sequencing run showed a higher sequencing depth than any other sample, whereas all samples extracted using the Maxwell methods showed low sequencing depths (Figures 2A and 3A). In the comparisons where the distribution of sequencing depths between DNA extraction methods was more even (ASV richness in the Scottish Agricultural College soil and functional and ARG richness in La Côte de Saint André soil), the DNA extraction methods did not show a significant effect on the richness discovery (Tables S4–S6 in the Supporting Information).

<figure><img src="image_6.png" /><figcaption>**Figure 3.** Effect of sequencing depth on the ARG richness discovery from (A) the Scottish Agricultural College soil (Craibstone, Scotland) and (B) La Côte de Saint André soil (France). Green: Maxwell 1 modification of the prototype Maxwell RSC Fecal Microbiome DNA Kit (Promega). Red: Maxwell 2 modification of the prototype Maxwell RSC Fecal Microbiome DNA Kit (Promega). Blue: Phenol/Chloroform. Yellow: DNeasy PowerSoil Kit (QIAGEN). The ARG richness measured in each sample was determined and plotted as a function of the sequencing depth. For each method/soil pair, triplicates (labelled 1, 2, and 3) were sequenced twice and plotted in the graph. Lighter colours represent the first sequencing run, and darker colours represent the second sequencing run. Correlation *p*-value: 3.37 × 10⁻¹⁰ (Scottish Agricultural College soil) and 2.35 × 10⁻¹² (La Côte de Saint André soil).</figcaption></figure>

The effects of sequencing depth and DNA extraction were compared at a more exhaustive level using Venn diagrams. At an overall level, the number of unique ASVs, functional classes and ARGs detected by a DNA extraction method increased with the sequencing depth (Figures S4–S6 in the Supporting Information). In addition, higher sequenced depths also increased the proportion of ASVs (Figure 4A–D), functional classes (Figure 5A–D), and ARGs (Figure 6A–D) discovered by the triplicates extracted using same method. This trend
to detect more unique ASVs, functional classes, and ARGs at higher sequencing depths was observed for all the DNA extraction methods tested in this study. The sole exceptions to this were the similar ASV discoveries from triplicates extracted using the Maxwell 1 modification of the prototype Maxwell RSC Fecal Microbiome DNA Kit (Promega) and sequenced at different depths (Figure 4A) and the higher functional class discoveries obtained at lower sequencing depths from triplicates extracted using the same method (Figure 5A). When the same sample was sequenced twice and differences in the sequencing depth were observed, the richness discovery was also influenced. More unique ASVs (Figure 4E–H), functional classes (Figure 5E–H), and ARGs (Figure 6E–H) were identified when the sample was sequenced at a higher depth. These differences in richness discovery between the same sample increased with increasing differences in the sequencing depth and were more evident in the metagenomic sequencing, where the sequencing depth was more uneven (Figures 4E–H and 5E–H). Finally, even though differences in the proportion of the soil ASV (Figure 4I–L) and overall functional (Figure 5I–L) and ARG (Figure 6I–L) richness unravelled by samples extracted using different methods and sequenced at similar depths were also found, the distribution of unique elements was more even than that observed when the same sample was sequenced at different depths. Regarding metagenomic sequencing, the number of functional classes and ARGs discovered by the two samples extracted using different methods and sequenced at a similar depth increased along with the sequencing depth (Figures 5I–L and 6I–L). However, whereas the number of functional classes that were detected by only one sample decreased along with the sequencing depth (Figure 4I–L), at higher sequencing depths, both the number of core and unique discovered ARGs increased (Figure 6I–L).

<figure><img src="image_7.png" /><figcaption>LA CÔTE DE SAINT ANDRE</figcaption></figure>

**Figure 4.** Venn Diagrams representing shared and unique ASVs from triplicates extracted using the same method sequenced in two different runs (A–D): individual samples sequenced at different depths (E–H) and samples extracted using different methods and sequenced at similar depths (I–L). Maxwell 1 and 2: modifications of the prototype Maxwell RSC Fecal Microbiome DNA Kit (Promega). PowerSoil: DNeasy PowerSoil Kit (QIAGEN).
<figure><img src="image_8.png" /><figcaption>**Figure 5.** Venn Diagrams representing shared and unique functional classes from triplicates extracted using the same method sequenced in two different runs (A–D): individual samples sequenced at different depths (E–H) and samples extracted using different methods and sequenced at similar depths (I–L). Maxwell 1 and 2: modifications of the prototype Maxwell RSC Fecal Microbiome DNA Kit (Promega). PowerSoil: DNeasy PowerSoil Kit (QIAGEN).</figcaption></figure>

<figure><img src="image_9.png" /><figcaption>**Figure 6.** Venn Diagrams representing shared and unique ARGs from triplicates extracted using the same method sequenced in two different runs (A–D): individual samples sequenced at different depths (E–H) and samples extracted using different methods and sequenced at similar depths (I–L). Maxwell 1 and 2: modifications of the prototype Maxwell RSC Fecal Microbiome DNA Kit (Promega). PowerSoil: DNeasy PowerSoil Kit (QIAGEN).</figcaption></figure>
4. Discussion

Although the development of soil metagenomic techniques implied a significant breakthrough in the understanding of soil microbial ecology, biases associated with these approaches prevent the accurate definition of the soil microbiome [44] and its associated functions in the soil ecosystem [45]. The aim of this study was to compare two main sources of biases in the metagenomic workflow: sequencing depth and DNA extraction methods and determine their relative contribution to bacterial richness discovery. Whereas sequencing depth had a significant impact on bacterial richness discovery from both soils analysed in the study at all the analysed levels, different sequencing runs and DNA extraction methods only showed significant differences when the sequencing depth was highly unequal between samples. In other words, these significant differences seem to be influenced by a higher sequencing depth of some samples rather than by the sequencing run or the DNA extraction method. Yet, clearly, when the sequencing depths were more evenly distributed between methods (i.e., when all methods showed a similar sequencing depth or all methods had at least one deeply sequenced sample), these differences became insignificant. This supports our hypothesis that the sequencing depth has a greater impact than the choice of the DNA extraction method on richness discovery at the ASV, functional class, and ARG level.

The slight differences in the bacterial richness measured in different samples at similar sequencing depths could be partially due to DNA extraction biases if they were only observed between samples extracted using different methods. However, DNA extraction triplicates (Figures 1–3) and sequencing duplicates (Figure 1B) sequenced at similar depths also detected different levels of bacterial richness. These differences could also be related to the distribution of ASVs, functional classes, and ARGs in soil ecosystems (Figures S1–S3 in the Supporting Information). Only a small proportion of annotated ASVs, functional classes, and ARGs are relatively abundant in the soils included in this study, while most of them are present at low abundances in extracted environmental DNA. This pattern of taxonomical distribution has already been observed in a study comparing soils from 237 different locations, where only 2% of the ensemble of bacterial phylotypes were found to be dominant [46]. Regarding functional classes, this distribution is consistent with a few functions being shared between different taxa and implicated in common bacterial ecology processes [47], while a large pool of low-abundant genes confers functions specific to single taxonomical groups. Finally, the few abundant ARGs in soil could have been selected, along with their carrying bacteria, in response to indigenous or exogenous selective pressures [48], whereas the majority of the natural reservoir of ARGs would represent a low-abundance pool. Since metagenomic techniques do not provide access to the whole bacterial diversity present in soil [49,50] and sequencing is a random subsampling from pools containing unevenly distributed taxa and genes, the proportion of high- and low-abundant elements sequenced from each sample could vary between samples, regardless of DNA extraction. This would lead to biased observations of bacterial richness that are likely more related to an insufficient sequencing depth than to DNA extraction limitations (Figures S4–S6 in the Supplementary Information).

The improved bacterial richness discovery from method triplicates (Figures 4A–D, 5A–D and 6A–D) and from individual samples (Figures 4E–H, 5E–H and 6E–H) at higher sequencing depths is consistent with a better detection of low-abundant elements at increasing sequencing depths. However, two exceptions were observed with this trend (Figures 4A and 5A). These inconsistencies are an illustration of the random subsampling of pools leading to variations in the richness discovery. In addition, the more even distribution of unique elements between samples extracted using different methods and sequenced at similar depths (Figures 4I–L, 5I–L and 6I–L) than between sequencing duplicates from the same sample at different sequencing depths (Figures 4E–H, 5E–H and 6E–H) supports the idea that the sequencing depth has a stronger effect on the richness discovery than DNA extraction. Finally, whereas more core and less unique functional classes were detected by two extraction methods sequenced at equal depths when this depth increased (Figure 5I–L),
and the number of both core and unique ARGs detected from both soils increased with the sequencing depth (Figure 6I–L). This supports the concern that short-read sequencing such as the one performed in this study provides an incomplete picture of the diversity of soil ARGs even at high sequencing depths and that the combination of long and short reads should be used to improve the ARG richness discovery [51,52].

Overall, our results confirm that the sequencing depth is a major source of biases in metagenomic studies and suggest that it has a stronger impact than DNA extraction on the richness discovery. Previous studies have shown that the sequencing depth has an impact on the bacterial richness and diversity discovery [26–28,53], but unfortunately, they were not coupled to the extraction method variants, since these have also been considered as critical for bacterial richness and diversity [16–18,54]. Thus, this is the first study (to the best of our knowledge) evaluating both effects simultaneously and demonstrating that the impact of sequencing depth on the soil bacterial richness discovery is greater than that of DNA extraction methods. This comparison can help determine which investments are more urgently needed to improve the metagenomics workflow. Whereas DNA extraction methods may affect the richness discovery, other factors, such as DNA amplification [55], library preparation, sequencing techniques [8], and sequence annotation [56], may also generate biases and lead to inaccurate comparisons between methods. We are currently unable to determine which of the differences observed between the DNA extraction methods are related to DNA extraction technique variability and which are a consequence of insufficient sequencing depths. Therefore, efforts should be made to optimise each of the steps in the metagenomic workflow in order to sequence representative samples of extracted DNA at a sufficient and even depth. This would not only facilitate an accurate comparison between DNA extraction methods but would also help define standard methods for soil metagenomics that would improve metagenomic comparisons and eventually lead to accurate profiles of soil microbiomes.

**Supplementary Materials:** The following supporting information can be downloaded at https://www.mdpi.com/article/10.3390/biom12030364/s1: Figure S1: Abundance of the ensemble of ASVs from DNA extracted from (A) the Scottish Agricultural College soil and (B) La Côte de Saint André soil. Figure S2: Total abundance of the ensemble of functional classes classified using SEED from DNA extracted from (A) the Scottish Agricultural College soil and (B) La Côte de Saint André soil. Figure S3: Total abundance of the ensemble of ARGs from DNA extracted from (A) the Scottish Agricultural College soil and (B) La Côte de Saint André soil. Figure S4: Venn Diagrams representing shared and unique ASV from (A) the Scottish Agricultural College soil, both sequencing runs, (B) La Côte de Saint André soil, both sequencing runs, (C) the Scottish Agricultural College soil, first sequencing run, and (D) La Côte de Saint André soil, first sequencing run. Figure S5: Venn Diagrams representing shared and unique functional classes from (A) the Scottish Agricultural College soil, both sequencing runs, (B) La Côte de Saint André soil, both sequencing runs, (C) the Scottish Agricultural College soil, first sequencing run, and (D) La Côte de Saint André soil, first sequencing run. Figure S6: Venn Diagrams representing shared and unique ARGs from (A) the Scottish Agricultural College soil, both sequencing runs, (B) La Côte de Saint André soil, both sequencing runs, (C) the Scottish Agricultural College soil, first sequencing run, and (D) La Côte de Saint André soil, first sequencing run. Table S1: Statistics of the comparison between sequencing run effects on ASV richness discovery. Table S2: Statistics of the comparison between sequencing run effects on functional richness discovery. Table S3: Statistics of the comparison between sequencing run effects on ARG richness discovery. Table S4: Statistics of the comparison between DNA extraction method effects on ASV richness discovery. Table S5: Statistics of the comparison between DNA extraction method effects on functional richness discovery. Table S6: Statistics of the comparison between DNA extraction method effects on ARG richness discovery.

**Author Contributions:** Conceptualisation, C.S.-C., L.F., L.D., T.S. and T.M.V.; Methodology, C.S.-C. and R.T.-P.; Formal analysis, C.S.-C., R.T.-P. and T.M.V.; Original draft preparation, C.S.-C., R.T.-P. and T.M.V.; Writing—Review and Editing, R.T.-P., L.F., L.D. and T.S.; and Supervision, T.M.V. All authors have read and agreed to the published version of the manuscript.
**Funding:** This research was funded by the ANRT (Association Nationale de la Recherche et de la Technologie) and Promega Corporation as part of the CIFRE PhD grant 2017/0486.

**Institutional Review Board Statement:** Not applicable.

**Informed Consent Statement:** Not applicable.

**Data Availability Statement:** The datasets generated and analysed during the current study are publicly available in the University of Lyon Environmental Microbial Genomics Group's repository: https://www.mmnt.net/db/0/0/ftp.ec-lyon.fr/pub/ADN/. Last Accessed date 10 January 2022.

**Acknowledgments:** We would like to thank Robin Walker for the soil samples from SRUC, Craibstone, Aberdeen. The Woodlands Field research plots are funded through the Scottish Government RESAS 2016-21 programme. We would also like to thank Christina Hazard for soil sampling from La Côte de Saint André (France).

**Conflicts of Interest:** This study was partially funded by Promega Corporation. C.S.-C. was employed by Promega Corporation during this study. L.D. and T.S. are still employed by Promega Corporation. The sponsors had no role in the design, execution, interpretation, or writing of the study. R.T.-P., L.F. and T.M.V. declare no conflicts of interest.

## References

1. Roesch, L.F.W.; Fulthorpe, R.R.; Riva, A.; Casella, G.; Hadwin, A.K.M.; Kent, A.D.; Daroub, S.H.; Camargo, F.A.O.; Farmerie, W.G.; Triplett, E.W. Pyrosequencing enumerates contrasts soil microbial diversity. *ISME J.* **2007**, *1*, 283–290. [CrossRef] [PubMed]

2. Rappé, M.S.; Giovannoni, S.J. The Uncultured Microbial Majority. *Annu. Rev. Microbiol.* **2003**, *57*, 369–394. [CrossRef] [PubMed]

3. Hugenholtz, P.; Goebel, B.M.; Pace, N.R. Impact of culture-independent studies on the emerging phylogenetic view of bacterial diversity. *J. Bacteriol.* **1998**, *180*, 4765–4774. [CrossRef] [PubMed]

4. Cruaud, P.; Vigneron, A.; Lucchetti-Miganeh, C.; Ciron, P.E.; Godfroy, A.; Cambon-Bonavita, M.A. Influence of DNA extraction method, 16S rRNA targeted hypervariable regions, and sample origin on microbial diversity detected by 454 pyrosequencing in marine chemosynthetic ecosystems. *Appl. Environ. Microbiol.* **2014**, *80*, 4626–4639. [CrossRef] [PubMed]

5. Klindworth, A.; Pruesse, E.; Schweer, T.; Peplies, J.; Quast, C.; Horn, M.; Glöckner, F.O. Evaluation of general 16S ribosomal RNA gene PCR primers for classical and next-generation sequencing-based diversity studies. *Nucleic Acids Res.* **2013**, *41*, e1. [CrossRef]

6. Hugerth, L.W.; Andersson, A.F. Analysing microbial community composition through amplicon sequencing: From sampling to hypothesis testing. *Front. Microbiol.* **2017**, *8*, 1561. [CrossRef] [PubMed]

7. Tremblay, J.; Singh, K.; Fern, A.; Kirton, E.S.; He, S.; Woyke, T.; Lee, J.; Chen, F.; Dangl, J.L.; Tringe, S.G. Primer and platform effects on 16S rRNA tag sequencing. *Front. Microbiol.* **2015**, *6*, 771. [CrossRef] [PubMed]

8. Lombard, N.; Prestat, E.; van Elsas, J.D.; Simonet, P. Soil-specific limitations for access and analysis of soil microbial communities by metagenomics. *FEMS Microbiol. Ecol.* **2011**, *78*, 31–49. [CrossRef]

9. Robe, P.; Nalin, R.; Capellano, C.; Vogel, T.M.; Simonet, P. Extraction of DNA from soil. *Eur. J. Soil Biol.* **2003**, *39*, 183–190. [CrossRef]

10. Levy-Booth, D.J.; Campbell, R.G.; Gulden, R.H.; Hart, M.M.; Powell, J.R.; Klironomos, J.N.; Pauls, K.P.; Swanton, C.J.; Trevors, J.T.; Dunfield, K.E. Cycling of extracellular DNA in the soil environment. *Soil Biol. Biochem.* **2007**, *39*, 2977–2991. [CrossRef]

11. Paget, E.; Monrozier, L.J.; Simonet, P. Adsorption of DNA on Clay-Minerals—Protection against DNaseI and Influence on Gene-Transfer. *FEMS Microbiol. Lett.* **1992**, *97*, 31–39. [CrossRef]

12. Pathan, S.I.; Arfaioli, P.; Ceccherini, M.T.; Ascher-Jenull, J.; Nannipieri, P.; Pietramellara, G.; D'Acqui, L.P. Physical protection of extracellular and intracellular DNA in soil aggregates against simulated natural oxidative processes. *Appl. Soil Ecol.* **2021**, *165*, 104002. [CrossRef]

13. Frostegård, Å.; Courtois, S.; Ramisse, V.; Clerc, S.; Bernillon, D.; Le Gall, F.; Jeannin, P.; Nesme, X.; Simonet, P. Quantification of bias related to the extraction of DNA directly from soils. *Appl. Environ. Microbiol.* **1999**, *65*, 5409–5420. [CrossRef] [PubMed]

14. Tebbe, C.C.; Vahjen, W. Interference of humic acids and DNA extracted directly from soil in detection and transformation of recombinant DNA from bacteria and a yeast. *Appl. Environ. Microbiol.* **1993**, *59*, 2657–2665. [CrossRef] [PubMed]

15. Wnuk, E.; Waśko, A.; Walkiewicz, A.; Bartminski, P.; Bejger, R.; Mielnik, L.; Bieganowski, A. The effects of humic substances on DNA isolation from soils. *PeerJ* **2020**, *8*, e9378. [CrossRef] [PubMed]

16. Kauffmann, I.M.; Schmitt, J.; Schmid, R.D. DNA isolation from soil samples for cloning in different hosts. *Appl. Microbiol. Biotechnol.* **2004**, *64*, 665–670. [CrossRef] [PubMed]

17. Miller, D.N.; Bryant, J.E.; Madsen, E.L.; Ghiorse, W.C. Evaluation and optimization of DNA extraction and purification procedures for soil and sediment samples. *Appl. Environ. Microbiol.* **1999**, *65*, 4715–4724. [CrossRef]

18. Töwe, S.; Wallisch, S.; Bannert, A.; Fischer, D.; Hai, B.; Haesler, F.; Kleineidam, K.; Schloter, M. Improved protocol for the simultaneous extraction and column-based separation of DNA and RNA from different soils. *J. Microbiol. Methods* **2011**, *84*, 406–412. [CrossRef]
19. Zielińska, S.; Radkowski, P.; Blendowska, A.; Ludwig-Gałężowska, A.; Łos, J.M.; Ł., M. The choice of the DNA extraction method may influence the outcome of the soil microbial community structure analysis. *Microbiologyopen* **2017**, *6*, e00453. [CrossRef]

20. Roosevelt-Amsaleg, C.L.; Garnier-Sillam, E.; Harry, M. Extraction and Purification of Microbial DNA from Soil and Sediment Samples. *Appl. Soil Ecol.* **2001**, *18*, 47–60. [CrossRef]

21. Sagova-Mareckova, M.; Cermak, L.; Novotna, J.; Plhackova, K.; Jana Forstova, J.K. Innovative Methods for Soil DNA Purification Tested in Soils with Widely Differing Characteristics. *Appl. Environ. Microbiol.* **2008**, *74*, 2902–2907. [CrossRef] [PubMed]

22. Bollmann-Giolai, A.; Giolai, M.; Heavens, D.; Macaulay, L.; Malone, J.; Clark, M.D. A low-cost pipeline for soil microbiome profiling. *Microbiologyopen* **2020**, *9*, e1133. [CrossRef]

23. McLaren, M.R.; Willis, A.D.; Callahan, B.J. Consistent and correctable bias in metagenomic sequencing experiments. *Elife* **2019**, *8*, e46923. [CrossRef] [PubMed]

24. Xie, K.; Deng, Y.; Zhang, X.; Wang, X.; Kang, G.; Bai, L.; Huang, H. Biases in Prokaryotic community amplicon sequencing affected by DNA extraction methods in both saline and non-saline soil. *Front. Microbiol.* **2018**, *9*, 1796. [CrossRef]

25. Costea, P.I.; Zeller, G.; Sunagawa, S.; Pelletier, E.; Alberti, A.; Levenez, F.; Tramontano, M.; Driessen, M.; Hercog, R.; Jung, F.E.; et al. Towards standards for human fecal sample processing in metagenomic studies. *Nat. Biotechnol.* **2017**, *35*, 1069–1076. [CrossRef]

26. Zaheer, R.; Noyes, N.; Polo, R.O.; Cook, S.R.; Marinier, E.; Doomselaar, G.V.; Belk, K.E.; Morley, P.S.; McAllister, T.A. Impact of sequencing depth on the characterization of the microbiome and resistome. *Sci. Rep.* **2018**, *8*, 5890. [CrossRef]

27. Smith, D.P.; Peay, K.G. Sequence depth, not PCR replication, improves ecological inference from next generation DNA sequencing. *PLoS ONE* **2014**, *9*, e90234. [CrossRef] [PubMed]

28. Lundin, D.; Severin, I.; Logue, J.B.; Östman, O.; Andersson, A.F.; Lindström, E.S. Which sequencing depth is sufficient to describe patterns in bacterial α- and β-diversity? *Environ. Microbiol. Rep.* **2012**, *4*, 367–372. [CrossRef] [PubMed]

29. Griffiths, R.I.; Whiteley, A.S.; O'Donnell, A.G.; Bailey, M.J. Rapid method for coextraction of DNA and RNA from natural environments for analysis of ribosomal DNA- and rRNA-based microbial community composition. *Appl. Environ. Microbiol.* **2000**, *66*, 5488–5491. [CrossRef]

30. Overbeek, R.; Begley, T.; Butler, R.M.; Choudhuri, J.V.; Chuang, H.Y.; Cohoon, M.; Crécy-Lagard, V.; Diaz, N.; Disz, T.; Edward, R.; et al. The subsystems approach to genome annotation and its use in the project to annotate 1000 genomes. *Nucleic Acids Res.* **2005**, *33*, 5691–5702. [CrossRef] [PubMed]

31. Alcock, B.P.; Raphenya, A.R.; Lau, T.T.Y.; Tsang, K.K.; Bouchard, M.; Edalatmand, A.; Huynh, W.; Nguyen, A.-L.V.; Cheng, A.A.; Liu, S.; et al. CARD 2020: Antibiotic resistome surveillance with the comprehensive antibiotic resistance database. *Nucleic Acids Res.* **2020**, *48*, D517–D525. [CrossRef] [PubMed]

32. Garbisu, C.; Garaiyurrebaso, O.; Lanzén, A.; Álvarez-Rodríguez, I.; Arana, L.; Blanco, F.; Smalla, K.; Grohmann, E.; Alkorta, I. Mobile genetic elements and antibiotic resistance in mine soil amended with organic wastes. *Sci. Total Environ.* **2018**, *621*, 725–733. [CrossRef] [PubMed]

33. Jechalke, S.; Schreiter, S.; Wolters, B.; Dealtry, S.; Heuer, H.; Smalla, K. Widespread dissemination of class 1 integron components in soils and related ecosystems as revealed by cultivation-independent analysis. *Front. Microbiol.* **2013**, *4*, 420. [CrossRef]

34. Kopmann, C.; Jechalke, S.; Rosendahl, I.; Groeneweg, J.; Krügerrecklenfort, E.; Zimmerling, U.; Weichelt, V.; Siemens, J.; Amelung, W.; Heuer, H.; et al. Abundance and transferability of antibiotic resistance as related to the fate of sulfadiazine in maize rhizosphere and bulk soil. *FEMS Microbiol. Ecol.* **2013**, *83*, 125–134. [CrossRef]

35. Stalder, T.; Press, M.O.; Sullivan, S.; Liachko, I.; Top, E.M. Linking the resistome and plasmidome to the microbiome. *ISME J.* **2019**, *13*, 2437–2446. [CrossRef] [PubMed]

36. Kemp, J.S.; Paterson, E.; Gammack, S.M.; Cresser, M.S.; Killham, K. Leaching of genetically modified Pseudomonas fluorescens through organic soils: Influence of temperature, soil pH, and roots. *Biol. Fertil. Soils* **1992**, *13*, 218–224. [CrossRef]

37. Muyzer, G.; Hottentrager, S.; Teske, A.; Wawer, C. Denaturing gradient gel electrophoresis of PCR-amplified 16S rDNA—A new molecular approach to analyse the genetic diversity of mixed microbial communities. In *Molecular Microbial Ecology Manual*; Akkermans, A., van Elsas, J., de Bruijn, F., Eds.; Kluwer Academic: Alphen aan den Rijn, The Netherlands, 1995; pp. 1–23.

38. Watanabe, K.; Kodama, Y.; Harayama, S. Design and evaluation of PCR primers to amplify bacterial 16S ribosomal DNA fragments used for community fingerprinting. *J. Microbiol. Methods* **2001**, *44*, 253–262. [CrossRef]

39. Callahan, B.J.; McMurdie, P.J.; Rosen, M.J.; Han, A.W.; Johnson, A.J.A.; Holmes, S.P. DADA2: High resolution sample inference from Illumina amplicon data. *Nat. Methods* **2016**, *13*, 581–583. [CrossRef]

40. Oksanen, J.; Blanchet, F.G.; Friendly, M.; Kindt, R.; Legendre, P.; McGlinn, D.; Minchin, P.R.; O'Hara, R.B.; Simpson, G.L.; Solymos, P.; et al. Vegan: Community Ecology Package, R Package Version 2.4-2. 2019. Available online: https://CRAN.R-project.org/package=vegan (accessed on 11 January 2022).

41. Chenn, H. Generate High-Resolution Venn and Euler Plots. In *VennDiagram Package*, 2018; Volume 33, p. 21. Available online: https://cran.r-project.org/web/packages/VennDiagram/index.html (accessed on 11 January 2022).

42. Benjamin Buchfink Xie, C.; Huson, D.H. Fast and sensitive protein alignment using DIAMOND. *Nat. Methods* **2015**, *12*, 59–60.

43. Huson, D.H.; Beier, S.; Flade, I.; Górska, A.; El-Hadidi, M.; Mitra, S.; Ruscheweyh, H.-J.; Tappu, R. MEGAN Community Edition-Interactive Exploration and Analysis of Large-Scale Microbiome Sequencing Data. *PLoS Comput. Biol.* **2016**, *12*, e1004957. [CrossRef] [PubMed]

44. Delmont, T.O.; Robe, P.; Cecillon, S.; Clark, I.M.; Constancias, F.; Simonet, P.; Hirsch, P.R.; Vogel, T.M. Accessing the soil metagenome for studies of microbial diversity. *Appl. Environ. Microbiol.* **2011**, *77*, 1315–1324. [CrossRef] [PubMed]
45. Nesme, J.; Achouak, W.; Agathos, S.N.; Bailey, M.; Baldrian, P.; Brunel, D.; Frostegård, Å.; Heulin, T.; Jansson, J.K.; Jurkevitch, E.; et al. Back to the future of soil metagenomics. *Front. Microbiol.* **2016**, *7*, 73. [CrossRef] [PubMed]

46. Delgado-Baquerizo, M.; Oliverio, A.M.; Brewer, T.E.; Benavent-González, M.; Eldridge, D.J.; Bardgett, R.D.; Maestre, F.T.; Singh, B.K.; Fierer, N. A global atlas of the dominant bacteria found in soil. *Science* **2018**, *359*, 320–325. [CrossRef]

47. Delmont, T.O.; Robe, P.; Clark, I.; Simonet, P.; Vogel, T.M. Metagenomic comparison of direct and indirect soil DNA extraction approaches. *J. Microbiol. Methods* **2011**, *86*, 397–400. [CrossRef] [PubMed]

48. Martínez, J.L. Bottlenecks in the transferability of antibiotic resistance from natural ecosystems to human bacterial pathogens. *Front. Microbiol.* **2012**, *2*, 256. [CrossRef] [PubMed]

49. Gweon, H.S.; Shaw, L.P.; Swann, J.; Maio, N.D.; AbuOun, M.; Niehus, R.; Hubbard, A.T.M.; Bowes, M.J.; Bailey, M.J.; Peto, T.E.A.; et al. The impact of sequencing depth on the inferred taxonomic composition and AMR gene content of metagenomic samples. *Environ. Microbiomes* **2019**, *14*, 7. [CrossRef]

50. Delmont, T.O.; Simonet, P.; Vogel, T.M. Describing microbial communities and performing global comparisons in the omic era. *ISME J.* **2012**, *6*, 1625–1628. [CrossRef]

51. Slizovskiy, I.B.; Mukherjee, K.; Dean, C.J.; Boucher, C.; Noyes, N.R. Mobilization of Antibiotic Resistance: Are Current Approaches for Colocalizing Resistomes and Mobilomes Useful? *Front. Microbiol.* **2020**, *11*, 1376. [CrossRef]

52. Peter, S.; Bosio, M.; Gross, G.; Bezdan, D.; Gutierrez, J.; Oberhettinger, P.; Liese, J.; Vogel, W.; Dörfel, D.; Berger, L.; et al. Tracking of Antibiotic Resistance Transfer and Rapid Plasmid Evolution in a Hospital Setting by Nanopore Sequencing. *mSphere* **2020**, *5*, e00525-20. [CrossRef] [PubMed]

53. Pereira-Marques, J.; Hout, A.; Ferreira, R.M.; Weber, M.; Pinto-Ribeiro, I.; Doorn, L.-J.V.; Knetsch, C.W.; Figueiredo, C. Impact of host DNA and sequencing depth on the taxonomic resolution of whole metagenome sequencing for microbiome analysis. *Front. Microbiol.* **2019**, *10*, 1277. [CrossRef] [PubMed]

54. Wang, Y.; Morimoto, S.; Ogawa, N.; Oomori, T.; Fujii, T. An improved method to extract RNA from soil with efficient removal of humic acids. *J. Appl. Microbiol.* **2009**, *107*, 1168–1177. [CrossRef] [PubMed]

55. Gohl, D.M.; Vangay, P.; Garbe, J.; MacLean, A.; Hauge, A.; Becker, A.; Gould, T.J.; Clayton, J.B.; Johnson, T.J.; Hunter, R.; et al. Systematic improvement of amplicon marker gene methods for increased accuracy in microbiome studies. *Nat. Biotechnol.* **2016**, *34*, 942–949. [CrossRef] [PubMed]

56. Kusnezowa, A.; Leichert, L.I. In silico approach to designing rational metagenomic libraries for functional studies. *BMC Bioinform.* **2017**, *18*, 267. [CrossRef]