Nested PCR to optimize rpoB metabarcoding for low-concentration and host-associated bacterial DNA

Nested PCR to optimize rpoB metabarcoding for low-concentration and host-associated bacterial DNA | bioRxiv /* */ /* */ <!-- <!-- /*! * yepnope1.5.4 * (c) WTFPL, GPLv2 */ (function(a,b,c){function d(a){return"[object Function]"==o.call(a)}function e(a){return"string"==typeof a}function f(){}function g(a){return!a||"loaded"==a||"complete"==a||"uninitialized"==a}function h(){var a=p.shift();q=1,a?a.t?m(function(){("c"==a.t?B.injectCss:B.injectJs)(a.s,0,a.a,a.x,a.e,1)},0):(a(),h()):q=0}function i(a,c,d,e,f,i,j){function k(b){if(!o&&g(l.readyState)&&(u.r=o=1,!q&&h(),l.onload=l.onreadystatechange=null,b)){"img"!=a&&m(function(){t.removeChild(l)},50);for(var d in y[c])y[c].hasOwnProperty(d)&&y[c][d].onload()}}var j=j||B.errorTimeout,l=b.createElement(a),o=0,r=0,u={t:d,s:c,e:f,a:i,x:j};1===y[c]&&(r=1,y[c]=[]),"object"==a?l.data=c:(l.src=c,l.type=a),l.width=l.height="0",l.onerror=l.onload=l.onreadystatechange=function(){k.call(this,r)},p.splice(e,0,u),"img"!=a&&(r||2===y[c]?(t.insertBefore(l,s?null:n),m(k,j)):y[c].push(l))}function j(a,b,c,d,f){return q=0,b=b||"j",e(a)?i("c"==b?v:u,a,b,this.i++,c,d,f):(p.splice(this.i++,0,a),1==p.length&&h()),this}function k(){var a=B;return a.loader={load:j,i:0},a}var l=b.documentElement,m=a.setTimeout,n=b.getElementsByTagName("script")[0],o={}.toString,p=[],q=0,r="MozAppearance"in l.style,s=r&&!!b.createRange().compareNode,t=s?l:n.parentNode,l=a.opera&&"[object Opera]"==o.call(a.opera),l=!!b.attachEvent&&!l,u=r?"object":l?"script":"img",v=l?"script":u,w=Array.isArray||function(a){return"[object Array]"==o.call(a)},x=[],y={},z={timeout:function(a,b){return b.length&&(a.timeout=b[0]),a}},A,B;B=function(a){function b(a){var a=a.split("!"),b=x.length,c=a.pop(),d=a.length,c={url:c,origUrl:c,prefixes:a},e,f,g;for(f=0;f<d;f++)g=a[f].split("="),(e=z[g.shift()])&&(c=e(c,g));for(f=0;f<b;f++)c=x[f](c);return c}function g(a,e,f,g,h){var i=b(a),j=i.autoCallback;i.url.split(".").pop().split("?").shift(),i.bypass||(e&&(e=d(e)?e:e[a]||e[g]||e[a.split("/").pop().split("?")[0]]),i.instead?i.instead(a,e,f,g,h):(y[i.url]?i.noexec=!0:y[i.url]=1,f.load(i.url,i.forceCSS||!i.forceJS&&"css"==i.url.split(".").pop().split("?").shift()?"c":c,i.noexec,i.attrs,i.timeout),(d(e)||d(j))&&f.load(function(){k(),e&&e(i.origUrl,h,g),j&&j(i.origUrl,h,g),y[i.url]=2})))}function h(a,b){function c(a,c){if(a){if(e(a))c||(j=function(){var a=[].slice.call(arguments);k.apply(this,a),l()}),g(a,j,b,0,h);else if(Object(a)===a)for(n in m=function(){var b=0,c;for(c in a)a.hasOwnProperty(c)&&b++;return b}(),a)a.hasOwnProperty(n)&&(!c&&!--m&&(d(j)?j=function(){var a=[].slice.call(arguments);k.apply(this,a),l()}:j[n]=function(a){return function(){var b=[].slice.call(arguments);a&&a.apply(this,b),l()}}(k[n])),g(a[n],j,b,n,h))}else!c&&l()}var h=!!a.test,i=a.load||a.both,j=a.callback||f,k=j,l=a.complete||f,m,n;c(h?a.yep:a.nope,!!i),i&&c(i)}var i,j,l=this.yepnope.loader;if(e(a))g(a,0,l,0);else if(w(a))for(i=0;i (function(w,d,s,l,i){w[l]=w[l]||[];w[l].push({'gtm.start':new Date().getTime(),event:'gtm.js'});var f=d.getElementsByTagName(s)[0];var j=d.createElement(s);var dl=l!='dataLayer'?'&l='+l:'';j.src='//www.googletagmanager.com/gtm.js?id='+i+dl;j.type='text/javascript';j.async=true;f.parentNode.insertBefore(j,f);})(window,document,'script','dataLayer','GTM-M677548'); Skip to main content Home About Submit ALERTS / RSS Search for this keyword Advanced Search New Results Nested PCR to optimize rpoB metabarcoding for low-concentration and host-associated bacterial DNA Laëtitia Leclerc , Gabryelle Agoutin , Thierry Brévault , View ORCID Profile Antony Champion , Jean Mainguy , Nicolas Nègre , Sudeeptha Yainna , Géraldine Pascal , Sophie Gaudriault , Jean-Claude Ogier doi: https://doi.org/10.1101/2025.03.17.643641 Laëtitia Leclerc 1 DGIMI, Univ Montpellier , INRAE, Montpellier, France Find this author on Google Scholar Find this author on PubMed Search for this author on this site Gabryelle Agoutin 2 GenPhySE, Université de Toulouse , INRAE, ENVT, 31326, Castanet Tolosan, France Find this author on Google Scholar Find this author on PubMed Search for this author on this site Thierry Brévault 3 AIDA, Univ Montpellier , CIRAD, Montpellier, France Find this author on Google Scholar Find this author on PubMed Search for this author on this site Antony Champion 4 DIADE, Univ Montpellier , IRD, Montpellier, France Find this author on Google Scholar Find this author on PubMed Search for this author on this site ORCID record for Antony Champion Jean Mainguy 5 Université de Toulouse , INRAE, BioinfOmics, GenoToul Bioinformatics facility, 31326, Castanet-Tolosan, France 6 Université de Toulouse , INRAE, UR 875 MIAT, 31326, Castanet-Tolosan, France Find this author on Google Scholar Find this author on PubMed Search for this author on this site Nicolas Nègre 1 DGIMI, Univ Montpellier , INRAE, Montpellier, France Find this author on Google Scholar Find this author on PubMed Search for this author on this site Sudeeptha Yainna 1 DGIMI, Univ Montpellier , INRAE, Montpellier, France Find this author on Google Scholar Find this author on PubMed Search for this author on this site Géraldine Pascal 2 GenPhySE, Université de Toulouse , INRAE, ENVT, 31326, Castanet Tolosan, France Find this author on Google Scholar Find this author on PubMed Search for this author on this site Sophie Gaudriault 1 DGIMI, Univ Montpellier , INRAE, Montpellier, France Find this author on Google Scholar Find this author on PubMed Search for this author on this site Jean-Claude Ogier 1 DGIMI, Univ Montpellier , INRAE, Montpellier, France Find this author on Google Scholar Find this author on PubMed Search for this author on this site For correspondence: jean-claude.ogier{at}umontpellier.fr Abstract Full Text Info/History Metrics Supplementary material Preview PDF Abstract Background Housekeeping genes have proved effective taxonomic markers for characterizing bacterial microbiota in short-read amplicon metabarcoding studies. A region of the rpoB gene, in particular, has been shown to minimize OTU overestimation bias with a high degree of accuracy, providing better species-level taxonomic resolution (Ogier et al. BMC Microbiol. 19:171). However, the primers for rpoB are highly degenerate, leading to potential problems in the amplification of bacterial DNA present at low concentration in the sample or embedded within eukaryotic matrices, as for the host-associated microbiota. We addressed these limitations by using a two-step PCR approach to optimize the rpoB procedure. The first PCR amplifies a 906-nucleotide region of the rpoB gene with the classical primers, referred to here as outer primers, and the second PCR then uses primers incorporating Illumina adapters, referred to here as inner primers, to amplify a 435-nucleotide subregion, the taxonomic marker for metabarcoding. Results We first used in silico approaches to evaluate the universality of the outer and inner rpoB primers. We then tested the nested rpoB PCR method on commercial mock samples of known composition. The nested PCR approach increased amplification efficiency for dilute samples without biasing the bacterial composition of the mock sample revealed by metabarcoding relative to single-step PCR. We also tested the nested rpoB PCR method on field-collected samples of the lepidopteran Spodoptera frugiperda . The nested PCR outperformed single-step PCR, increasing amplification efficiency for bacterial DNA present at low concentrations (oral secretions from S. frugiperda ) or embedded in eukaryotic DNA matrices ( S. frugiperda larvae). Conclusions This method provides a promising new strategy for characterizing insect-associated microbiota that can also be applied to other host microbiomes. Background Microbial communities are frequently characterized by polymerase chain reaction (PCR) amplification followed by next-generation sequencing (NGS), particularly through short-read metabarcoding approaches, typically targeting variable regions of the 16S rRNA gene, such as the V3-V4 regions. However, short-amplicon sequencing studies based on the 16S rRNA gene are subject to certain limitations, particularly in terms of taxonomic resolution, which is often restricted to the genus and family levels [ 1 , 2 ]. Full-length 16S rRNA gene sequencing with PacBio technology has been shown to improve taxonomic resolution in human microbiome samples [ 3 ], but its use is currently limited due to the higher cost required to obtain an equivalent number of reads per sample. Housekeeping genes encoding conserved proteins ubiquitous in the bacterial kingdom provide a possible alternative to the short amplicon of the 16S rRNA gene. These genes evolve much more rapidly than the 16S rRNA gene and can therefore be used to differentiate between bacterial lineages [ 4 ]. Moreover, unlike the 16S rRNA gene, housekeeping genes are generally present as a single copy in genomes, minimizing biases in diversity assessment and providing a more accurate representation of bacterial community composition. Housekeeping genes have proved effective for the characterization of bacterial diversity, particularly at the species level, across various ecosystems [ 4 – 8 ]. The rpoB marker, in particular, has been successfully used in metabarcoding analyses of samples rich in Pseudomonadota (formerly Proteobacteria), such as entomopathogenic nematodes, and has been shown to have a finer taxonomic resolution than the 16S (V3V4) marker [ 5 , 9 ]. In particular, phylogenies based on the same rpoB subregion are highly concordant with core genome phylogenies [ 10 ], indicating that this subregion is a robust marker for precise bacterial identification. One key challenge in the metabarcoding of microbial communities is detecting prokaryotic DNA in samples with very low bacterial DNA concentrations, or in samples in which eukaryotic DNA predominates, as in the host-associated microbiota. In most cases, physically separating microbes from host tissues is difficult or even impossible, making global DNA extraction the only option. However, global DNA extraction typically results in a low microbial-to-host DNA ratio, limiting the efficiency of amplification with bacterial primers. This issue is further exacerbated by the use of primers targeting housekeeping genes, such as rpoB , which has a highly variable sequence, necessitating the use of degenerate bases to maximize species coverage. This may increase cross-reactivity with eukaryotic DNA, further hampering target amplification. Finally, in metabarcoding protocols, the long sequence and GC-rich regions of barcoding primers (about 60 nucleotides, including the target region, index sequences, and Illumina adapters) reduce the specificity of hybridization, making it even more challenging to amplify bacterial DNA from a host DNA matrix. Nested PCR can be used to increase amplification sensitivity. In nested PCR approaches, there are two successive reactions, the first to amplify a larger region, and the second to amplify a smaller region nested within the first amplicon. Nested PCR is particularly useful for samples with low target DNA concentrations [ 11 ], but this method has only rarely been used in microbial ecology [ 12 ], and its potential for metabarcoding has yet to be explored. In this study, we evaluated the potential of a nested PCR approach combined with a metabarcoding analysis with rpoB primers for characterizing bacterial communities in samples with low bacterial DNA concentrations or in which the bacterial DNA is embedded in a eukaryotic matrix. For validation of the efficacy of this strategy, we compared the performance of the single-step and nested rpoB PCR methods on a set of mock samples. We then compared the efficiency of the nested rpoB PCR approach to that of single-step PCR on biological samples consisting of oral secretions (OS) from the lepidopteran Spodoptera frugiperda , which contain low concentrations of bacterial DNA, and S. frugiperda larvae in which the bacterial DNA is diluted in large amounts of eukaryotic DNA. Materials and Methods Biological material The two commercial mock samples contained DNA from eight bacterial strains in equal proportions (Mock_8sp; ZymoBIOMICS Microbial Community DNA Standard; D6305), and various proportions DNA from the same eukaryotic species (Mock_8sp log; ZymoBIOMICS Microbial Community DNA Standard II_log distribution; D6311). The composition of the mock communities is described in Table 1 . View this table: View inline View popup Table 1. Composition of the commercial DNA mock samples Wild S. frugiperda larvae and their oral secretions (OS) were collected from maize ( Zea mays , Poaeceae) fields near Velingara in the Anambé Valley in Southern Senegal (13.163268/-14.039137) in September 2021. We collected OS from approximately 100 larvae with flexible forceps that had been disinfected with 70% ethanol, and sterile plastic Pasteur pipettes. The samples were pooled to obtain a final volume of 200 µL and stored at approximately 4°C for a few hours before being frozen at −20°C. Extraction of bacterial DNA from insect samples Total DNA was extracted from both the OS pool (about 200 µL) and six individual larvae (posterior half body): three biological replicates at stage L3 and three at stage L5. OS and larval samples were lysed by freezing at −80°C for 15 min and then heating at 80°C for 20 min. We then added 100 μL (500 µL for the larval samples) Quick Extract lysis solution (Bacterial DNA extraction kit from Epi-centre, USA), 1 μL (2 µL for the larval samples) Ready-Lyse Lysozyme Solution (Epi-centre, USA) and 20 µL (50µL for the larval samples) EDTA (0.5 M, pH 8). The samples were mechanically lysed with six 2 mm glass beads and three cycles of 7 m/s for 40 s in a FastPrep device (MP Biomedicals, Illkirch-Graffenstaden, France). The complete lysis of insect and prokaryotic cells was achieved by incubation at room temperature for 90 min and then at 55°C with 10 µL (20 µL for the larval samples) proteinase K (20 mg/mL), with shaking at 1200 rpm until the solution cleared. RNA was eliminated by incubation with 10 μL (20 µL for the larval samples) RNase A at 20 mg/mL (Invitrogen PureLinkTM RNaseA, France) for 15 min at 37°C. After a final incubation with 100 µL protein precipitation solution (Wizard kit, Promega, France) for 5 min on ice, the samples were centrifuged for 5 min at 17,000 x g and 4°C and the supernatant was recovered. OS samples were subjected to phenol and chloroform treatments due to the presence of polysaccharides. The DNA present in the OS and larval samples was precipitated with absolute ethanol and washed with 70% ethanol. The extracted DNA was resuspended in 50 μL ultrapure water and stored at −20°C. We assessed the levels of contaminating bacterial DNA at each of the various steps in DNA sample preparation, by including several negative control extractions with sterile ultrapure water. Library preparation and Illumina MiSeq sequencing We used two PCR strategies to generate rpoB amplicon libraries for MiSeq sequencing: single-step PCR, involving direct amplification, and nested PCR, consisting of two successive PCRs ( Fig. 1 ). Each PCR was performed in duplicate (mock samples) or triplicate (OS samples). Our protocol also included RNase free water as a negative control and positive E. mundtii genomic DNA as a positive control. Download figure Open in new tab Figure 1. Protocol for the amplification of a region of the rpoB gene by (A) single-step PCR and (B) nested PCR strategies. The positions of the primers are indicated by green arrows (Uni_ rpoB _deg_F/R primers) and blue arrows ( rpoB _F/R primers). nt , nucleotides. For the single-step rpoB PCR, we used the existing Univ_ rpoB _deg_F (5’ - TCG TCG GCA GCG TCA GAT GTG TAT AAG AGA CAG N GGYTWYGAAGTNCGHGACGTDCA - 3’) and Univ_ rpoB _deg_R ( 5’ - GTC TCG TGG GCT CGG AGA TGT GTA TAA GAG ACA GN T GACGYTGCATGTTBGMRCCCATMA - 3’) primers, which generate a 435 bp rpoB amplicon. These primers include a 5’ nucleotide overhang with adapters for Illumina library indexing, indicated by the bases in italics in the primer sequences above. The conditions for the single-step PCR rpoB were as previously described [ 5 ]. Briefly, 35 cycles of amplification were performed in a Bio-Rad thermocycler with 1 to 50 ng genomic DNA, the high-fidelity iProof™ DNA Polymerase (Bio-Rad), and an annealing temperature of 57⍰°C. For the nested rpoB PCR, we designed degenerate consensus outer rpoB primers — rpoB _F (forward primer 5’ - CAGYTDTCNCARTTYATGGAYCA -3’) and rpoB _R (reverse primer 5’ - AGTTRTARCCDTYCCANGKCAT - 3’) — based on ClustalW alignments ( http://multalin.toulouse.inra.fr/multalin/ ) of rpoB gene sequences from bacterial reference genomes representing a broad taxonomic diversity of eubacteria. The binding sites of the selected primers correspond to Escherichia coli K12 nucleotide positions 1528 to 2434, resulting in the amplification of 906-nucleotide portion of the rpoB gene. The specificity of the rpoB primers was checked in silico as previously described [ 5 ]. The first round of amplification in the nested PCR was performed with the rpoB_F/R primers and the GoTaq DNA Polymerase (Promega) in a Bio-Rad thermocycler, with 1 to 50 ng genomic DNA, in accordance with the enzyme manufacturer’s instructions. Twenty-five amplification cycles were performed (annealing temperature of 50 °C and annealing time of 45 s), generating an amplicon of 906 bp. The second round of amplification in the nested PCR was performed with 1 µL of the product of the first round of PCR, with the Univ_ rpoB _deg_F/R primers incorporating Illumina adapters (see above) and the high-fidelity iProof™ DNA Polymerase (Bio-Rad). The same protocol was followed as for the single-step PCR (annealing temperature of 57 °C and annealing time of 45 s), except that the number of amplification cycles was reduced to 15. The presence of contaminating DNA was assessed in each PCR run by including negative controls with sterile ultrapure water as the template. The amounts of amplicon DNA and the size of the amplicons were systematically analyzed by agarose gel electrophoresis. All Illumina-indexed amplicons of the appropriate size on agarose gels were purified, multiplexed, and sequenced at the Genseq platform (University of Montpellier, France). The amplicons were purified with magnetic beads (Clean PCR, Proteigene, France) and a new round of PCR was performed in a total volume of 18 µL (5 µL of product from the first round of PCR, 9 µL Phusion® High-Fidelity PCR Master Mix, NEB, France, 2 µL index adapter I5, 2 µL index adapter I7). The cycling conditions were as follows: 95°C for 3 minutes, then 10 cycles of 95°C for 30 seconds, 55°C for 30 seconds, 72°C for 30 seconds, and a final elongation phase for 5 minutes at 72°C. We used a set of 384 index pairs based on the IDT’s unique dual index (UD) set to mark all samples and multiplex them in a single MiSeq run. The final PCR products were purified with magnetic beads, multiplexed and sequenced in pairs on an Illumina MiSeq sequencer with the MiSeq v3 reagent kit (600 cycles; Illumina) according to the manufacturer’s instructions. FastQ files were generated with bcltofastq at the end of the run. The quality of the raw data was assessed with a module developed by the MBB platform (Montpellier Bioinformatics Biodiversity platform supported by the LabEx CeMEB an ANR “Investissements d’avenir” program (ANR-10-LABX-04-01)) running the multiqc program [ 13 ]. The raw sequence data can be downloaded from http://www.ebi.ac.uk/ena/ (accession numbers: PRJEB81812) Bioinformatic processing of sequence data The sequenced reads were analyzed with FROGS v4.1.0 software [ 14 ]. We successively used the pre-processing, clustering_swarm, remove_chimera, cluster_filter, taxonomic_affiliation tools with the parameters shown in Table S1. Update of the rpoB database We created an updated rpoB sequence database by downloading the genomes with the genome_updater tool ( https://github.com/pirovc/genome_updater ). With this tool, we were able to create a RefSeq database [ 15 ] containing only complete genomes, complete chromosomes and protein sequences. We obtained protein sequences from 46,711 bacteria. We then used FetchMGS ( https://github.com/motu-tool/fetchMGs.pl ) [ 16 ] to extract DNA and protein sequence files for rpoB . We obtained 47,069 rpoB nucleic acid sequences from 46,298 bacteria. The rpoB database is available from the FROGS website: https://genoweb.toulouse.inrae.fr/frogs_databanks/assignation/rpoB/ . In silico PCR analysis for evaluating the universality of rpoB primers We used OBITools suite v1.2.11 [ 17 ] to evaluate the universality of the primers. The results were first formatted manually to make the identifiers unique and compatible with EcoPCR v1.0.1 [ 18 ], the in silico PCR tool. We then used obiconvert to construct a database that could be used by EcoPCR. We obtained three different types of amplicons targeting rpoB , by performing different runs of EcoPCR using: (i) the Univ_ rpoB _deg_F and Univ_ rpoB _deg_R primers (see above for the primer sequences) to obtain typical metabarcoding amplicons [ 5 ], (ii) the nested PCR primers to perform two successive PCR, one with rpoB _F and rpoB _R (see above) and the other based on previous output sequences with the Univ_ rpoB _deg_F and Univ_ rpoB _deg_R primers. For comparisons of taxonomic specificity, we compared only species to which these markers (amp_direct_ rpoB , amp_nested_primary_ rpoB , amp_nested_ rpoB ) were common. In total, 47,069 genomes were sufficiently complete to contain the rpoB gene. The bacterial diversity of the amplicons obtained is displayed with the krona tool [ 19 ]. All these steps and scripts are detailed on github: https://github.com/geraldinepascal/RPOB_paper/tree/main/pcr_in_silico_script . Bacterial community and statistical analyses ASV diversity and statistical analyses were performed with the R packages Phyloseq [ 20 ], Vegan [ 21 ], Microbiota Process [ 22 ] and Ampvis2 [ 23 ], as previously described [ 9 ]. Two custom-written scripts were used to generate plots to display diversity by taxonomy rank: add_multiaffi_to_abd_table.py ( https://github.com/geraldinepascal/RPOB_paper/blob/main/FROGS_analysis_results/scripts/add_multiaffi_to_abd_table.py ) and plot_taxo_ranks.py ( https://github.com/geraldinepascal/RPOB_paper/blob/main/FROGS_analysis_results/scripts/plot_taxo_ranks.py ). Results Design of single-step and nested rpoB PCR strategies for metabarcoding The single-step and nested PCR strategies for the amplification of the rpoB gene are illustrated in Fig. 1 . For the single-step PCR strategy ( Fig. 1A ), the region used for metabarcoding is directly amplified (35 cycles) with the Uni_ rpoB _deg_F/R inner primers containing Illumina adapters [ 5 ] to yield a 435 bp amplicon. For the nested PCR strategy ( Fig. 1B ), we first designed the “ rpoB _F/R” outer primers to bind to conserved regions of the rpoB gene encompassing the metabarcoding target region and generating a 906 bp amplicon. This amplicon was then subjected to a second amplification reaction with the inner primers incorporating Illumina adapters (Uni_ rpoB _deg_F/R primers), generating a 435 bp amplicon. The first PCR should enrich the substrate for the second reaction in the rpoB target, thereby increasing the efficiency of the second round of PCR. We optimized the number of amplification cycles, using genomic DNA from biological samples available in the laboratory: 25 cycles for the first PCR and 15 for the second (data not shown). This optimization minimized the total number of cycles while ensuring a robust signal for Illumina sequencing and no detectable signal in the PCR negative control. In silico evaluation of the taxonomic coverage of rpoB primers We checked the universality of the outer and inner rpoB primers in silico against an rpoB database with 47,069 entries, this database having been updated in 2024 ( https://web-genobioinfo.toulouse.inrae.fr/frogs_databanks/assignation/rpoB/ ). The single-step and nested strategies had the potential to amplify 72% and 68%, respectively, of the 47,069 sequences stored in the rpoB database ( Fig. 2 ). A Krona representation identified the bacterial taxa potentially amplified by the rpoB primers in the single-step PCR strategy ( Fig. 2A ) and after the first (Fig. S1) and second rounds ( Fig. 2B ) of the nested PCR strategy. The outer primers used in the first PCR of the nested approach had a high degree of taxonomic coverage, encompassing over 88% of all rpoB sequences, accounting for the similar levels of taxonomic coverage achieved with the single-step and nested strategies. The Pseudomonadota and most of the Bacillota, excluding the Staphylococcaceae , were well covered by the successive use of outer rpoB primers in the first PCR and inner rpoB primers in the second PCR. By contrast, the Actinomycetota and Bacteroidota were poorly covered by both the single-step and nested PCR strategies. For more details on taxonomic coverage, an interactive link is available from https://github.com/geraldinepascal/RPOB_paper/tree/main/krona_figure . Download figure Open in new tab Figure 2. Coverage of bacterial taxonomic diversity, evaluated with the EcoPCR v1.0.1 in silico tool for (A) the single-step rpoB PCR strategy and (B) the nested rpoB PCR strategy. The in silico PCR was performed on the 47,069 rpoB sequences stored in the rpoB database. The pie charts illustrate the number of sequences amplified with two mismatches allowed (green) and the number of unretrieved sequences (pink) from the 47,069 sequences. Below, the Krona representations display the major taxa potentially amplified by the rpoB marker (green) in the single-step and nested PCR strategies. Taxa that are not amplified are shown in pink. For more details on the taxa covered by the two PCR strategies, all Krona visualizations are accessible via the interactive link: https://github.com/geraldinepascal/RPOB_paper/tree/main/krona_figure . Experimental comparison of the single-step and nested rpoB PCR strategies with different mock samples We compared the sensitivity and reliability of nested and single-step PCR strategies on two commercial mock samples — Mock_8sp and Mock_8sp_log (see Table 1 for composition) — containing DNA from eight bacterial species, non-diluted (ND) or diluted 1/10 or 1/100. After amplification by PCR, gel electrophoresis of the rpoB amplicons showed that the nested PCR strategy resulted in successful amplification, even at dilutions as high as 1:100, for both mock samples ( Fig. 3A and 3B ). By contrast, with the single-step PCR strategy, amplification was possible only up to a dilution of 1/10 dilution for Mock_8sp ( Fig. 3A ) and only with the undiluted sample for the Mock_8sp_log ( Fig. 3B ). The nested PCR strategy was, therefore, more sensitive than the direct amplification approach. No amplification was observed from the negative control, regardless of the strategy used. Download figure Open in new tab Figure 3. Comparison of single-step and nested rpoB PCR strategies with commercial mock communities used as DNA matrices to assess amplification sensitivity ( A and B ) and the ability to determine taxonomic composition by metabarcoding ( B and C ). Panels A and B show agarose gel electrophoresis of rpoB amplicons from mock_8sp and mock_8sp_log samples, respectively, amplified by single-step or nested PCR strategies. The dilutions of the mock samples are indicated on the agarose gel, below the wells. For each dilution, two technical replicates are shown. ND: not diluted. Panels C and D show heatmaps for the bacterial taxonomic composition obtained by rpoB metabarcoding of the mock_8sp and mock_8sp_log samples, respectively (only samples with successful amplification were sequenced). The columns ‘Single’ and ‘Nested’ correspond to the single-step and nested rpoB PCR strategies, and the column ‘T’ indicates the expected composition of the mocks. The dilutions of the mock samples are indicated at the bottom of the heatmap. The 20 most abundant ASVs across samples, assigned at species level, are listed on the left. Each species name is followed by an RDP bootstrap confidence score, indicating the confidence level of the taxonomic assignment. The relative read abundance is represented by a gradient of red hues. We then sequenced the amplicons visible on agarose gels with Illumina MiSeq technology. The number of reads per sample and the ASV composition are detailed in Table S2; rarefaction curves are shown in Fig. S2. The relative abundance of ASVs by taxonomic rank clearly confirmed that the rpoB primers provided a high degree of taxonomic resolution, with 100% of taxa assigned to species level in the mock samples (Fig. S3). The bacterial compositions obtained for Mock_8sp after the single-step and nested amplification strategies closely resembled the theoretical composition of the mock sample ( Fig. 3C and 3D ). However, several additional sequences were assigned to the species Shigella flexneri in Mock_8sp and Listeria innocua in both mock samples, regardless of the strategy used. These additional assignments probably correspond to artifactual sequence variants, attributable to the close phylogenetic relationship between the genera Shigella and Escherichia , or between L. innocua and L. monocytogenes . However, neither E. faecalis nor S. aureus was detected in Mock 8sp_log by either single-step or nested PCR, consistent with the low proportion of genomes attributable to these species in Mock 8sp_log, below 0.1%, the threshold generally required for signal detection by an amplicon metabarcoding approach. We then compared the relative abundance of taxa at species level obtained for the Mock_8sp samples with the single-step and nested rpoB PCR strategies (Fig. S4). Both PCR strategies overestimated the abundances of E. coli, Enterococcus faecalis , and Listeria monocytogenes , while significantly underestimating the abundances of Limosilactobacillus fermentum and Staphylococcus aureus . This result is consistent with the in silico analysis of the taxonomic diversity targeted by the rpoB primers, which revealed limited coverage for certain bacterial groups within the Bacillota (see above). The single-step PCR strategy significantly overestimated the abundance of Salmonella enterica, whereas the nested PCR strategy did not (Fig. S4). Overall, our results demonstrate that the nested PCR strategy provides taxonomic assignment profiles similar to those generated by the single-step method, but with a greater sensitivity and more efficient amplification for samples containing a low concentration of the target DNA. Application of the nested rpoB strategy to characterization of the microbiota of insect samples We experimentally validated the nested PCR strategy on an insect, Spodoptera frugiperda , sampled in Senegal (near Velingara). We tested our methodology on two types of biological matrix: (i) DNA extracted from the oral secretion (OS) of larvae, a biological fluid composed of saliva, the gut microbiota, and a food bolus [ 24 ], which contains only very low concentrations (< 4 µg/mL) of bacterial DNA, and (ii) DNA extracted from demi-larvae (caterpillars), in which the bacterial DNA is embedded in a eukaryotic DNA matrix. The rpoB region was targeted by both the single-step and nested PCR strategies. PCR amplifications were performed on the six replicate samples of larval DNA (three stage L3, three stage L5), on DNA from the pooled OS sample, for which three PCR replicates were processed, and on the extraction control samples. The single-step rpoB PCR strategy generated no visible amplicon for either the OS or larval samples, which it was therefore impossible to sequence for metabarcoding. By contrast, the nested rpoB strategy successfully amplified all OS and larval samples, yielding sufficient DNA for metabarcoding. Importantly, no amplification was observed for the extraction or PCR negative controls following the two successive PCR rounds of the nested rpoB strategy. All amplicons generated by the rpoB nested strategy were further subjected to Illumina MiSeq sequencing and statistical analysis. Details of the read count per sample and ASV composition are provided in Table S2, and rarefaction curves are provided in Fig. S5. The Pseudomonadota largely dominated in the pooled OS sample, accounting for three quarters of the observed diversity, with Bacillota accounting for the remaining quarter ( Fig. 4A ). Two species predominated ( Fig. 4B ): a species belonging to the genus Ochrobactrum , but with an insufficient RDP identification confidence threshold for reliable assignment to a species, and a species confidently identified (RDP bootstrap confidence of 1) as E. casseliflavus . There were slight variations between technical PCR replicates, highlighting the importance of such replicates for ensuring reliable and reproducible results for the metabarcoding of low-concentration DNA samples. Download figure Open in new tab Figure 4. Composition of the bacterial community according to metabarcoding of oral secretion (OS) samples from S. frugiperda amplified by the nested rpoB PCR strategy. ( A ) At phylum level, the relative abundances are represented on bar plots, with sample names indicated at the top of the plots. ( B ) At species level, the 30 most abundant ASVs across samples are listed on the left of the heatmap. Each species name is followed by an RDP bootstrap confidence score, indicating the confidence level of the taxonomic assignment. The relative abundance (percentage) is represented by a gradient of red hues. Sample names, corresponding to the three PCR replicates, are indicated at the top of the heatmap columns. For larval samples, the nested rpoB strategy revealed a predominance of Pseudomonadota in five of the six biological replicates, with only one L5-stage sample containing a majority of Bacillota ( Fig. 5A ). The nested rpoB strategy performed well for identification to species level (i.e., RDP assignment threshold generally > 0.8), making it possible to identify the species associated with S. frugiperda larvae with confidence ( Fig. 5B ). In particular, the metabarcoding analysis revealed the prevalence of Klebsiella quasipneumoniae, Acinetobacter soli, E. casseliflavus, K. pneumoniae, and A. baumannii in both L3-stage and L5-stage caterpillars. Interestingly, a few taxa, such as E. casseliflavus and K. quasipneumoniae, were common to both the larval and OS samples. Download figure Open in new tab Figure 5. Composition of the bacterial community revealed by metabarcoding of S. frugiperda demi-larva samples subjected to amplification with the nested rpoB PCR strategy. ( A ) At phylum level, relative abundances are represented on bar plots, with sample names indicated at the top of the plots. At the bottom of the graph, L3 corresponds to the third developmental stage, and L5 corresponds to the fifth developmental stage of the caterpillar ( B ) At species level, the 30 most abundant ASVs across samples are listed on the left of the heatmap. Each species name is followed by an RDP bootstrap confidence score, indicating the confidence level of the taxonomic assignment. The relative abundance (percentage) is represented by a gradient of red hues. Sample names are indicated at the top of the heatmap columns. Discussion Metabarcoding is a powerful tool for simultaneously identifying multiple taxa within a DNA sample. However, short-read metabarcoding studies, which typically rely on amplicons of variable regions of the 16S rRNA gene to characterize bacterial microbiota [ 25 – 27 ], lack the precision required for fine-scale taxonomic resolution [ 28 , 29 ]. In this study, we focused on the rpoB marker due to its ability to provide species-level taxonomic assignments [ 5 ], providing a promising alternative to 16S rRNA-based approaches. However, the primers for rpoB have several weaknesses. Due to the lack of highly conserved regions in the rpoB gene, the primers contain multiple degenerate bases, which can make specific amplification difficult, particularly in matrices with low concentrations of bacterial gDNA and in host-associated samples, which generally contain much less prokaryotic than eukaryotic DNA. For example, even with 35 amplification cycles, PCR with rpoB primers fails to amplify the target efficiently in some complex host-associated microbiota, such as insect samples. Bivand and coll. [ 30 ] tried to improve the specificity of the primers for target binding by designing rpoB primers according to the dual priming oligonucleotide (DPO) principle, and coupling the amplification with Sanger sequencing to detect pathogenic taxa in clinical samples dominated by human DNA. This approach was effective for targeted detection but is not suitable for NGS-based metabarcoding. In this study, we validated a nested PCR strategy to improve the amplification efficiency and sensitivity of rpoB primers for metabarcoding analysis. Unlike the standard single-step rpoB PCR, this nested PCR involves a first PCR to enrich the sample in the target DNA and a second PCR to amplify the rpoB fragments for metabarcoding. Nested PCR has long been known to increase PCR sensitivity [ 31 ] but, to our knowledge, this is the first time that this method has been used in a metabarcoding approach. We first demonstrate the reliability of nested PCR with the rpoB marker for the accurate identification of bacterial species in different mock samples of known composition and variable complexity, and its ability to increase amplification sensitivity relative to the single-step PCR strategy. Some studies have reported that nested PCR has the drawback of potentially increasing bias through preferential amplification in the two successive PCR steps [ 32 ]. By optimizing the number of cycles between the first and second PCRs (40 cycles in total), we minimized the biases associated with nested PCR, ensuring robust results comparable to those obtained by single-step PCR (35 cycles in total). We then compared the performance of the single-step and nested rpoB PCR strategies for describing the host-associated bacterial microbiota in two types of lepidopteran insect samples: (i) S. frugiperda larval oral secretions (OS), corresponding to a matrix containing low concentrations of bacterial DNA and, (ii) S. frugiperda larvae, corresponding to a eukaryotic matrix with low proportions of bacterial DNA. The single-step rpoB PCR failed to amplify the target, whereas the rpoB nested PCR strategy successfully amplified the target from all insect samples. The rpoB nested PCR strategy is, therefore, considerably more sensitive than the single-step rpoB PCR, for samples with low concentrations of DNA and samples in which eukaryotic DNA predominates over bacterial DNA. However, our findings also reveal several limitations of the rpoB primers, such as their limited coverage of certain bacterial phyla, contrasting with the universality of 16S primers. Indeed, the rpoB primers cover only about 70% of bacterial diversity (essentially the phyla Pseudomonodota and Bacillota). For example, rpoB lacks sensitivity for the detection of Bacteroidota, members of which may be present in host-associated microbiota samples, such as insect samples. As a means of addressing this bias, we propose the design of a multi-marker approach using phylum-specific rpoB primer pairs to achieve an optimal coverage of bacterial diversity. Conclusion In metabarcoding studies, combining a nested PCR strategy with the use of the rpoB marker appears to be a promising approach for describing the bacterial complexity and richness of insect-associated microbiota directly at species level. Our study highlights the benefits of nested PCR for efficiently amplifying prokaryotic DNA from eukaryotic hosts. The nested PCR rpoB metabarcoding method, with its species-level resolution and enhanced sensitivity, offers new opportunities for investigating the ecological roles of specific bacterial taxa in insect physiology, health, and environmental interactions. Finally, it would be of interest to extend this strategy to other host-associated microbiota samples. This study provides valuable insights for the growing scientific community exploring bacterial communities through metabarcoding across diverse ecosystems. List of abbreviations DNA deoxyribonucleic acid OS oral secretion ASV amplicon sequence variants PCR polymerase chain reaction RDP Classifier Ribosomal Database Project classifier Declarations Ethics approval and consent to participate Not applicable Consent for publication Not applicable Availability of data and material The sequence data supporting the findings of this study have been deposited in the European Nucleotide Archive under primary accession code PRJEB81812. Competing interests The authors have no competing interests to declare. Funding This work was funded by the SPE ‘MicroSpodoRice’ grant from INRAE. Jean Mainguy’s work was supported by “La Région Occitanie” and the European Union via the “Regional Platforms of Research and Innovation” call for proposals of the Occitanie Region as part of the Operational Program FEDER-FSE MIDI-PYRENEES ET GARONNE 2014–2020. Authors’ contributions L.L., S.G., J.C.O. conceived and designed the experiments. L.L, J.C.O. conducted the experiments. L.L., S.G., J.C.O. analyzed the data. J.C.O., G.A., J.M. and G.P. conducted the bioinformatics analysis and analyzed the data. T.B. organized and led the Senegal mission. LL and SY collected FAW larvae in Senegal. LL collected the OS in the field. A.C., N.N, S.G and J.C.O. supervised the research and contributed to the discussion of the results. L.L., S.G., J.C.O. wrote the manuscript. L.L., J.C.O., G.A., G.P. prepared the figures and supplementary files. All authors reviewed the manuscript. Supplemental Materials Table S1. FROGS parameters Table S2. Detailed amplicon sequence variants (ASV) table for the samples analyzed Fig. S1. Krona representations of the bacterial taxa potentially amplified by the outer primers used in the first round of PCR in the nested rpoB PCR strategy. The in silico EcoPCR v1.0.1 tool was used on the 47,069 rpoB sequences stored in the rpoB database, with two mismatches allowed. The Krona graph displays the major taxa potentially amplified by the outer rpoB primers (in green). Taxa that are not amplified are shown in pink. For more details on the taxa covered by the two PCR strategies, all Krona visualizations are accessible via the interactive link: https://github.com/geraldinepascal/RPOB_paper/tree/main/krona_figure . Fig. S2. Rarefaction curves for the metabarcoding sequences of the different commercial mock samples amplified by the single-step or nested rpoB PCR strategy. For each sample, species richness is shown on the y -axis and the number of sequences is shown on the x -axis. The ggrare function was used to generate the rarefaction curves ( https://rdrr.io/github/gauravsk/ranacapa/man/ggrare.html ). Fig. S3. Plot showing the relative abundance of ASVs at different taxonomic ranks for the various commercial mock samples following amplification by the single-step or nested rpoB PCR strategy. Two custom scripts were used to generate the plot: add_multiaffi_to_abd_table.py ( https://github.com/geraldinepascal/RPOB_paper/blob/main/FROGS_analysis_results/scripts/add_multiaffi_to_abd_table.py ) and plot_taxo_ranks.py ( https://github.com/geraldinepascal/RPOB_paper/blob/main/FROGS_analysis_results/scripts/plot_taxo_ranks.py ). Fig. S4. Histogram plot of bacterial composition showing the relative abundances of ASVs at species level in the non-diluted commercial mock_8sp sample, amplified by the single-step or nested rpoB PCR strategy. Two technical replicates (rep1 and rep2) were performed, with the last bar indicating the expected composition. Fig. S5. Rarefaction curves for the metabarcoding sequences of S. frugiperda OS samples ( panel A ) and S. frugiperda demi-larva samples ( panel B ) subjected to amplification by the single-step or nested rpoB PCR strategy. The OS samples correspond to three independent PCR replicates. The larval samples correspond to several biological replicates: L3 is the third developmental stage, and L5 is the fifth developmental stage of the caterpillar. For each sample, species richness is displayed on the y -axis and the number of sequences on the x -axis. The ggrare function was used to generate the rarefaction curves ( https://rdrr.io/github/gauravsk/ranacapa/man/ggrare.html ). Acknowledgments We would like to thank the Insect Quarantine Platform (PIQ), a member of the Vectopole Sud network, for providing the infrastructure required for experiments on insect pests. We would also like to thank Gaëtan Clabots, Raphaël Bousquet and Dylan Valenza for maintaining the FAW colonies at the DGIMI laboratory in Montpellier. The metabarcoding data were produced with the GenSeq technical facilities at the “Institut des Sciences de l’Evolution” in Montpellier, with support from LabEx CeMEB, an ANR “Investissements d’avenir” program (ANR-10-LABX-04-01). We would like to thank the Genotoul bioinformatics platform Toulouse Occitanie (Bioinfo Genotoul, https://doi.org/10.15454/1.5572369328961167E12 ) for providing assistance, computing and storage resources. We would like to thank Maïmouna Cissoko and Nils Poulicard for their help during field sampling in Senegal and Victoria Mariotti for her help with bioinformatics analyses. References 1. ↵ Schirmer M , Ijaz UZ , D’Amore R , Hall N , Sloan WT , Quince C . Insight into biases and sequencing errors for amplicon sequencing with the Illumina MiSeq platform . Nucleic Acids Research . 2015 ; 43 : e37 . OpenUrl CrossRef PubMed 2. ↵ Klindworth A , Pruesse E , Schweer T , Peplies J , Quast C , Horn M , et al. Evaluation of general 16S ribosomal RNA gene PCR primers for classical and next-generation sequencing-based diversity studies . Nucleic Acids Res . 2013 ; 41 : e1 . OpenUrl CrossRef PubMed 3. ↵ Buetas E , Jordán-López M , López-Roldán A , D’Auria G , Martínez-Priego L , De Marco G , et al. Full-length 16S rRNA gene sequencing by PacBio improves taxonomic resolution in human microbiome samples . BMC Genomics . 2024 ; 25 : 310 . 4. ↵ Poirier S , Rué O , Peguilhan R , Coeuret G , Zagorec M , Champomier-Vergès M-C , et al. Deciphering intra-species bacterial diversity of meat and seafood spoilage microbiota using gyrB amplicon sequencing: A comparative analysis with 16S rDNA V3-V4 amplicon sequencing . PLoS One . 2018 ; 13 : e0204629 . OpenUrl CrossRef PubMed 5. ↵ Ogier J-C , Pagès S , Galan M , Barret M , Gaudriault S . rpoB , a promising marker for analyzing the diversity of bacterial communities by amplicon sequencing . BMC Microbiol . 2019 ; 19 : 171 . 6. Roux S , Enault F , le Bronner G , Debroas D. Comparison of 16S rRNA and protein-coding genes as molecular markers for assessing microbial diversity (Bacteria and Archaea) in ecosystems . FEMS Microbiology Ecology . 2011 ; 78 : 617 – 28 . OpenUrl CrossRef PubMed 7. Vos M , Quince C , Pijl AS , de Hollander M , Kowalchuk GA . A comparison of rpoB and 16S rRNA as markers in pyrosequencing studies of bacterial diversity . PLoS One . 2012 ; 7 : e30600 . OpenUrl CrossRef PubMed 8. ↵ Barret M , Briand M , Bonneau S , Préveaux A , Valière S , Bouchez O , et al. Emergence Shapes the Structure of the Seed Microbiota . Appl Environ Microbiol . 2015 ; 81 : 1257 – 66 . OpenUrl Abstract / FREE Full Text 9. ↵ Ogier J-C , Pagès S , Frayssinet M , Gaudriault S . Entomopathogenic nematode-associated microbiota: from monoxenic paradigm to pathobiome . Microbiome . 2020 ; 8 : 25 . 10. ↵ Hassler HB , Probert B , Moore C , Lawson E , Jackson RW , Russell BT , et al. Phylogenies of the 16S rRNA gene and its hypervariable regions lack concordance with core genome phylogenies . Microbiome . 2022 ; 10 : 104 . 11. ↵ Thompson A , Bench S , Carter B , Zehr J. Chapter Three Coupling FACS and Genomic Methods for the Characterization of Uncultivated Symbionts . 2013 . doi: 10.1016/b978-0-12-407863-5.00003-4 . OpenUrl CrossRef 12. ↵ Ogier J-C , Son O , Gruss A , Tailliez P , Delacroix-Buchet A . Identification of the Bacterial Microflora in Dairy Products by Temporal Temperature Gradient Gel Electrophoresis . AEM . 2002 ; 68 : 3691 – 701 . OpenUrl CrossRef 13. ↵ Ewels P , Magnusson M , Lundin S , Käller M . MultiQC: summarize analysis results for multiple tools and samples in a single report . Bioinformatics . 2016 ; 32 : 3047 – 8 . OpenUrl CrossRef PubMed 14. ↵ Escudié F , Auer L , Bernard M , Mariadassou M , Cauquil L , Vidal K , et al. FROGS1: Find, Rapidly, OTUs with Galaxy Solution . Bioinformatics . 2018 ; 34 : 1287 – 94 . OpenUrl CrossRef PubMed 15. ↵ O’Leary NA , Wright MW , Brister JR , Ciufo S , Haddad D , McVeigh R , et al. Reference sequence (RefSeq) database at NCBI: current status, taxonomic expansion, and functional annotation . Nucleic Acids Res . 2016 ; 44 : D733 – 745 . OpenUrl CrossRef PubMed 16. ↵ Milanese A , Mende DR , Paoli L , Salazar G , Ruscheweyh H-J , Cuenca M , et al. Microbial abundance, activity and population genomic profiling with mOTUs2 . Nat Commun . 2019 ; 10 : 1014 . OpenUrl CrossRef PubMed 17. ↵ Boyer F , Mercier C , Bonin A , Le Bras Y , Taberlet P , Coissac E. obitools: a unix-inspired software package for DNA metabarcoding . Molecular Ecology Resources . 2016 ; 16 : 176 – 82 . OpenUrl CrossRef PubMed 18. ↵ Ficetola GF , Coissac E , Zundel S , Riaz T , Shehzad W , Bessière J , et al. An In silico approach for the evaluation of DNA barcodes . BMC Genomics . 2010 ; 11 : 434 . 19. ↵ Ondov BD , Bergman NH , Phillippy AM . Interactive metagenomic visualization in a Web browser . BMC Bioinformatics . 2011 ; 12 : 385 . 20. ↵ McMurdie PJ , Holmes S. phyloseq: an R package for reproducible interactive analysis and graphics of microbiome census data . PLoS One . 2013 ; 8 : e61217 . OpenUrl CrossRef PubMed 21. ↵ Dixon P . VEGAN, a package of R functions for community ecology . Journal of Vegetation Science . 2003 ; 14 : 927 – 30 . OpenUrl CrossRef Web of Science 22. ↵ Xu S , Zhan L , Tang W , Wang Q , Dai Z , Zhou L , et al. MicrobiotaProcess: A comprehensive R package for deep mining microbiome . Innovation (Camb ). 2023 ; 4 : 100388 . OpenUrl PubMed 23. ↵ Andersen KS , Kirkegaard RH , Karst SM , Albertsen M . ampvis2: an R package to analyse and visualise 16S rRNA amplicon data. preprint . Bioinformatics ; 2018 . 24. ↵ Peiffer M , Felton GW . Do Caterpillars Secrete “Oral Secretions”? J Chem Ecol . 2009 ; 35 : 326 – 35 . OpenUrl CrossRef PubMed Web of Science 25. ↵ Bukin YS , Galachyants YP , Morozov IV , Bukin SV , Zakharenko AS , Zemskaya TI . The effect of 16S rRNA region choice on bacterial community metabarcoding results . Sci Data . 2019 ; 6 : 190007 . OpenUrl CrossRef PubMed 26. Zemb O , Achard CS , Hamelin J , De Almeida M , Gabinaud B , Cauquil L , et al. Absolute quantitation of microbes using 16S rRNA gene metabarcoding: A rapid normalization of relative abundances by quantitative PCR targeting a 16S rRNA gene spike-in standard . Microbiologyopen . 2020 ; 9 : e977 . OpenUrl CrossRef PubMed 27. ↵ Roume H , Mondot S , Saliou A , Le Fresne-Languille S , Doré J . Multicenter evaluation of gut microbiome profiling by next-generation sequencing reveals major biases in partial-length metabarcoding approach . Sci Rep . 2023 ; 13 : 22593 . 28. ↵ Pan P , Gu Y , Sun D-L , Wu QL , Zhou N-Y . Microbial Diversity Biased Estimation Caused by Intragenomic Heterogeneity and Interspecific Conservation of 16S rRNA Genes . Appl Environ Microbiol . 89 : e02108 – 22 . 29. ↵ Hrovat K , Dutilh BE , Medema MH , Melkonian C . Taxonomic resolution of different 16S rRNA variable regions varies strongly across plant-associated bacteria . ISME Communications . 2024 ; 4 : ycae034 . OpenUrl 30. ↵ Bivand JM , Dyrhovden R , Sivertsen A , Tellevik MG , Patel R , Kommedal Ø . Broad-range amplification and sequencing of the rpoB gene: a novel assay for bacterial identification in clinical microbiology . Journal of Clinical Microbiology . 2024 ; 62 : e00266 – 24 . OpenUrl PubMed 31. ↵ Song Y . PCR-based diagnostics for anaerobic infections . Anaerobe . 2005 ; 11 : 79 – 91 . OpenUrl CrossRef PubMed Web of Science 32. ↵ Fan Z-Y , Li X-R , Mao D-P , Zhu G-F , Wang S-Y , Quan Z-X . Could nested PCR be applicable for the study of microbial diversity? World J Microbiol Biotechnol . 2009 ; 25 : 1447 – 52 . OpenUrl CrossRef View the discussion thread. Back to top Previous Next Posted March 17, 2025. Download PDF Supplementary Material Email Thank you for your interest in spreading the word about bioRxiv. NOTE: Your email address is requested solely to identify you as the sender of this article. Your Email * Your Name * Send To * Enter multiple addresses on separate lines or separate them with commas. You are going to email the following Nested PCR to optimize rpoB metabarcoding for low-concentration and host-associated bacterial DNA Message Subject (Your Name) has forwarded a page to you from bioRxiv Message Body (Your Name) thought you would like to see this page from the bioRxiv website. Your Personal Message CAPTCHA This question is for testing whether or not you are a human visitor and to prevent automated spam submissions. Share Nested PCR to optimize rpoB metabarcoding for low-concentration and host-associated bacterial DNA Laëtitia Leclerc , Gabryelle Agoutin , Thierry Brévault , Antony Champion , Jean Mainguy , Nicolas Nègre , Sudeeptha Yainna , Géraldine Pascal , Sophie Gaudriault , Jean-Claude Ogier bioRxiv 2025.03.17.643641; doi: https://doi.org/10.1101/2025.03.17.643641 Share This Article: Copy Citation Tools Nested PCR to optimize rpoB metabarcoding for low-concentration and host-associated bacterial DNA Laëtitia Leclerc , Gabryelle Agoutin , Thierry Brévault , Antony Champion , Jean Mainguy , Nicolas Nègre , Sudeeptha Yainna , Géraldine Pascal , Sophie Gaudriault , Jean-Claude Ogier bioRxiv 2025.03.17.643641; doi: https://doi.org/10.1101/2025.03.17.643641 Citation Manager Formats BibTeX Bookends EasyBib EndNote (tagged) EndNote 8 (xml) Medlars Mendeley Papers RefWorks Tagged Ref Manager RIS Zotero Tweet Widget Facebook Like Google Plus One Subject Area Microbiology Subject Areas All Articles Animal Behavior and Cognition (7642) Biochemistry (17715) Bioengineering (13907) Bioinformatics (42003) Biophysics (21470) Cancer Biology (18624) Cell Biology (25533) Clinical Trials (138) Developmental Biology (13390) Ecology (19935) Epidemiology (2067) Evolutionary Biology (24356) Genetics (15617) Genomics (22529) Immunology (17753) Microbiology (40432) Molecular Biology (17200) Neuroscience (88681) Paleontology (667) Pathology (2840) Pharmacology and Toxicology (4828) Physiology (7653) Plant Biology (15161) Scientific Communication and Education (2046) Synthetic Biology (4304) Systems Biology (9826) Zoology (2271)