# mbctools: A user-friendly metabarcoding and cross-platform pipeline for analyzing multiple amplicon sequencing data across a large diversity of organisms

Christian Barnabé, Guilhem Sempéré, Vincent Manzanilla, Joel Moo-Millan, Antoine Amblard-Rambert, Etienne Waleckx

► To cite this version:

Christian Barnabé, Guilhem Sempéré, Vincent Manzanilla, Joel Moo-Millan, Antoine Amblard-Rambert, et al.. mbctools: A user-friendly metabarcoding and cross-platform pipeline for analyzing multiple amplicon sequencing data across a large diversity of organisms. Peer Community Journal, 2024, 4, pp.e114. 10.24072/pcjournal.501 . hal-05182581

HAL Id: hal-05182581

https://hal.science/hal-05182581v1

Submitted on 23 Jul 2025

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
Software tool

Published
2024-12-11

**Cite as**
Christian Barnabé, Guilhem Sempéré, Vincent Manzanilla, Joel Moo Millan, Antoine Amblard-Rambert and Etienne Waleckx (2024) mbctools: A User-Friendly Metabarcoding and Cross-Platform Pipeline for Analyzing Multiple Amplicon Sequencing Data across a Large Diversity of Organisms, Peer Community Journal, 4: e114.

**Correspondence**
guilhem.sempere@cirad.fr
vincent.manzanilla@ird.fr
etienne.waleckx@ird.fr

**Peer-review**
Peer reviewed and
recommended by
PCI Genomics,
https://doi.org/10.24072/pci-
genomics.100370

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

This article is licensed
under the Creative Commons
Attribution 4.0 License.

# mbctools: A User-Friendly Metabarcoding and Cross-Platform Pipeline for Analyzing Multiple Amplicon Sequencing Data across a Large Diversity of Organisms

Christian Barnabé<fnref n="1" />,#,<fnref n="1" />, Guilhem Sempéré<fnref n="1" />,#,²,³, Vincent Manzanilla<fnref n="1" />,<fnref n="1" />, Joel Moo Millan<fnref n="1" />,⁴, Antoine Amblard-Rambert<fnref n="1" />,⁴, and Etienne Waleckx<fnref n="1" />,⁵,<fnref n="1" />,⁴

Volume 4 (2024), article e114

https://doi.org/10.24072/pcjournal.501

## Abstract

We developed a python package called *mbctools*, designed to offer a cross-platform tool for processing amplicon data from various organisms in the context of metabarcoding studies. It can handle the most common tasks in metabarcoding pipelines such as paired-end merging, primer trimming, quality filtering, sequence denoising, zero-radius operational taxonomic unit (ZOTU) filtering, and has the capability to process multiple genetic markers simultaneously. *mbctools* is a menu-driven program that eliminates the need for expertise in command-line skills and ensures documentation of each analysis for reproducibility purposes. The software, designed to run in a console, offers an interactive experience, guided by keyboard inputs, assisting users along the way through data processing and hiding the complexity of command lines by letting them concentrate on selecting parameters to apply in each step of the process. In our workflow, VSEARCH is utilized for processing *fastq* files derived from amplicon-based Next-Generation Sequencing data. This software is a versatile open-source tool for processing amplicon sequences, offering advantages such as high speed, efficient memory usage, and the ability to handle large datasets. It provides functions for various tasks such as dereplication, clustering, chimera detection, and taxonomic assignment. VSEARCH is thus very efficient in retrieving the overall diversity of a sample. To adapt to the diversity of projects in metabarcoding, we facilitate the reprocessing of datasets with the possibility to adjust parameters. *mbctools* can also be launched in a headless mode, making it suited for integration into pipelines running on High-Performance Computing environments. *mbctools* is available at https://github.com/GuilhemSempere/mbctools, https://pypi.org/project/mbctools/.

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

<footnote n="1">¹Institut de Recherche pour le Développement, UMR INTERTRYP IRD, CIRAD, Université de Montpellier, Montpellier, France,
²Centre de Coopération Internationale en Recherche Agronomique pour le Développement, UMR INTERTRYP IRD, CIRAD, Université de Montpellier, Montpellier, France,
³South Green Bioinformatics Platform, Biodiversity, Montpellier, France,
⁴Laboratorio de Parasitología, Centro de Investigaciones Regionales "Dr Hideyo Noguchi", Universidad Autónoma de Yucatán, Mérida, Yucatán, México,
#Equal contribution</footnote>
Introduction

Metabarcoding consists in identifying a targeted subset of genomes within bulk samples by massive sequencing of amplicons of taxonomically informative genetic markers, also known as barcodes (Figure 1) (Taberlet et al., 2012). Over the past decade, this approach has quickly gained popularity since Hebert et al. (2003) first advocated the use of short variable DNA sequences, amplified using universal primers, for species identification, and discovery of new taxa. Its simplicity of implementation enables biologists, local governments, and NGOs to increase their understanding of spatial and temporal ecological networks (Holdaway et al., 2017). Using universal primers specific of kingdoms or subgroups of organisms of interest (e.g. "arthropods", "vertebrates", "reptiles", "amphibians", "mammals", "birds", ...), metabarcoding is highly effective for species-level identification and assessment of the whole diversity of a sample within the targeted kingdoms or subgroups of organisms of interest. For example, universal primers for vertebrates, targeting the 12S rRNA gene, have been proposed for metabarcoding studies focused on vertebrates (Riaz et al., 2011). In studies targeting animals, primers allowing the amplification of a portion of the mitochondrial marker Cytochrome oxidase 1 (COI) may be used (Hebert et al., 2016). In plants, standard DNA barcoding generally involves one to four plastid DNA regions (rbcL, matK, trnH-psbA, trnL), sometimes in combination with the internal transcribed spacers of nuclear ribosomal DNA (nrDNA, ITS) (CBOL Plant Working Group et al., 2009). For fungi, the internal transcribed spacer (ITS) regions ITS1 and ITS2 of the nuclear ribosomal DNA are the most commonly used barcodes (Schoch et al., 2012). For bacterial and archaeal communities, the 16S ribosomal RNA (16S rRNA) is the most widely used barcode and provides taxonomic resolution to the genus level (Caporaso et al., 2011). Beyond these markers, specific markers are used in various fields to genotype species. For instance, in the study of *Trypanosoma cruzi*, the causative agent of Chagas disease, the Glucose-6-Phosphate Isomerase (GPI) and the Cytochrome oxidase 2 (COII) genes are frequently used to distinguish between different strains (Lauthier et al., 2012; López-Domínguez et al., 2022; Barnabé et al., 2023). Although single markers are very informative, in many cases, a combination of markers is necessary to fully understand the spatial and temporal dynamics of ecological networks. A huge advantage of metabarcoding is that it is based on massive sequencing. As a consequence, all the amplicons amplified from different targeted markers can be pooled and sequenced simultaneously (Hernández-Andrade et al., 2019). This powerful approach has proven to be useful for diet profile analysis, air and water quality monitoring, biodiversity monitoring, and food quality control, e.g., beverage or ancient ecosystem composition (Leray & Knowlton, 2015; Thomsen & Willerslev, 2015; Raclariu et al., 2017). As another application example, our group has recently proposed to use this approach to untangle the transmission cycles of vector-borne pathogens and their dynamics, using the gut contents of vector insects as bulk samples, and simultaneously sequencing markers allowing for the molecular identification of vector species and/or genetic diversity, blood meal sources (using universal primers for vertebrates which may also serve as pathogen hosts), gut microbiome composition (which may modulate vectorial capacity), and pathogen diversity. Indeed, all the latter components interact together to shape transmission cycles (Dumonteil et al., 2018; Hernández-Andrade et al., 2019).

The increasing demand for efficient and versatile pipelines to analyze amplicon data has led to the development of over 30 programs and pipelines in recent years (Hakimzadeh et al., 2023). However, these applications often present some challenges in terms of reproducibility and user-friendliness for non-expert users in the field of bioinformatics. Indeed, the obtaining of pseudo quantitative taxonomic data from fastq files involves a series of complex processing steps that can be achieved with specific software implying the understanding of a large panel of parameters. To address this problem, mbctools was designed with a simplified interface, allowing scientists without advanced scripting skills to easily install, navigate, utilize the pipeline while focusing on key parameters. This accessibility ensures that non-expert users can effectively and autonomously
analyze their data without experiencing difficulties in reproducibility or requiring extensive bioinformatics knowledge.

Among pipelines and algorithms used in the field of metabarcoding data processing, substantial variations in sensitivity and specificity are observed (Bailet et al., 2020; Mathon et al., 2021). From these benchmarking studies, VSEARCH appears appropriate for our applications. mbctools' main purpose is to facilitate the use of VSEARCH, which has demonstrated strong performance in handling large-scale next-generation sequencing (NGS) data (Rognes et al., 2016), and it thus could be considered a "wrapper", hiding the complexity of command-line to the users. It thus allows them to easily process sequencing reads and obtain meaningful outputs, i.e, zero-radius operational taxonomic units (ZOTUs (Edgar, 2016), almost equivalent to ASVs (amplicon sequence variants)) along with abundance information.

mbctools processes demultiplexed metabarcoding data obtained from one or multiple genetic markers to distribute cleaned and filtered sequences into ZOTUs sorted per marker, which can be used for taxonomic assignment or further analyzes such as phylogeny (Figures 1, 2). While the original goal in developing mbctools was providing assistance to novice users through its interface, it can also prove useful to bioinformaticians willing to integrate it as part of a wider data processing pipeline, since its series of mandatory operations (initial analysis described below) can be launched in a headless mode by feeding the program with a configuration file. Additionally, addressing the challenge of processing large numbers of markers simultaneously across different kingdoms was one of the main motivations for designing mbctools, which allows for the simultaneous processing of multiple genetic markers and has a simplified graphical interface that guides the user through key parameters of the analysis.

## Prerequisites and dependencies

### Input data types - sequencing tech & lab

mbctools enables the analysis of amplicon data from multiple markers, obtained from short reads (Illumina) or long reads (Oxford Nanopore Technology, PacBio). It can handle both single-end reads and paired-end reads, provided that the latter can be merged. In the eventuality where the merging of paired-end reads is prevented by either large amplicon lengths (leading to insufficient overlap zone) or the presence of low complexity sequences, mbctools is able to utilize only forward reads for analysis. As it relies on the underlying VSEARCH software which is highly optimized in terms of efficiency and memory utilization, it is able to accommodate large datasets deriving from numerous samples and amplicons.

### Running environment

mbctools is a versatile cross-platform python package running seamlessly on Windows, Unix, or OSX systems. The installation of Python3.7 (or higher) and VSEARCH (version 2.19.0 or higher, available at https://github.com/torognes/vsearch) is required to run the pipeline, and those binaries must be added to the PATH environment variable. Additionally, for Windows users, Powershell script execution must be enabled using the Set-ExecutionPolicy command, as the default Windows command prompt does not provide enough flexibility to run the software. The mbctools package, available from the Python Package Index (PyPI) (https://pypi.org/project/mbctools/), can be installed and upgraded using the command line tool pip (pip3 install mbctools). Before starting the analyses, the package requires mandatory directories and files to be pre-arranged as illustrated in Figure 3:

A default directory, named 'fastq', should contain all the fastq files corresponding to the R1 (and R2, when working with paired-ends sequencing) reads. These files are expected to be demultiplexed, meaning that each fastq file should contain reads from a single sample.
<figure><img src="image_4.png" /><figcaption>Figure 1 - Metabarcoding workflow example and area of operation of mbctools.</figcaption></figure>
<figure><img src="image_5.png" /><figcaption>**Figure 2** - Overview of the mbctools workflow, where each station icon represents a distinct step in the process. In the initial analysis (menu 1), data processing begins with the merging for paired-end reads (step 1) and the conversion of reads from `fastq` to `fasta` (step 2). This is followed by dereplication (step 3), denoising (step 4, which can be re-run with various options via menu 1a, 1b, 1c, 1d), chimera removal (step 5), sequence affiliation to markers (step 6), and finally, sequence orientation harmonization (step 7). Following this, primer removal (step 8) and denoised sequence abundance filtering (step 9) are applied. The last step in the analysis involves obtaining per-marker `fasta` files containing retained ZOTUs (step 10). The menu labeled as "4" comes into play after taxonomic assignment, tuning outputs for visualization in `metaXplor` and ensuring long-term data accessibility.</figcaption></figure>

A default directory named 'refs' that contains one reference file per genetic marker in the multi-FASTA format. Each of them should contain one or more unaligned reference sequences, from the same strand and without gaps. Providing multiple reference sequences for a given marker improves diversity support and thus leads to a better assignment accuracy. In addition, the 'sequence orientation fixing' step runs faster when using several reference sequences. Filenames must end with the extension '.fas' and must be called by the name of the corresponding marker, referred to as locus in mbctools (Figure 3).

'lociPE.txt' provides the list of markers (one per line) corresponding to paired-end R1/R2 reads.

'lociSE.txt' provides the list of markers (one per line) corresponding to single-end reads and for which only R1 reads will be analyzed.

'samples.txt' includes a list of sample names, where each name appears on a separate line.

'fastqR1.txt' contains the names of the `fastq` R1 files corresponding to the samples.

'fastqR2.txt' contains the names of the `fastq` R2 files corresponding to the samples.

The contents of the last three files must be sorted in the same order, i.e., line (i) of 'samples.txt' shall correspond to line (i) of both 'fastqR1.txt' and 'fastqR2.txt'.

mbctools shall be launched from within the working directory in a shell console, typically by typing 'mbctools' if installed as a package.
<figure><img src="image_6.png" /><figcaption>**Figure 3** - Pre-required files and directories</figcaption></figure>

## Software features and interface overview

mbctools features a series of menus, each containing two to four options as illustrated in Figures 4 to 8. On startup, the program displays the main menu shown in Figure 4. This menu presents descriptive keywords that enable users to easily navigate through the available menus and sections. From each menu, the user is required to select one single option at a time.

<figure><img src="image_7.png" /><figcaption>**Figure 4** - Contents of the main menu, explaining navigation conventions and providing access to the four main sections.</figcaption></figure>

To analyze a new dataset, it is necessary to perform a complete initial analysis by selecting the first option from both the first and subsequent menu. The initial analysis processes the data through a succession of steps such as merging, sample-level dereplication, sequence denoising, chimera detection and removal, affiliation of sequences to markers (e.g., loci), and sequence re-orientation (menu 1, Figure 5). Once those are completed, primers are removed, and the denoising process may be refined by tuning the cluster size threshold (menu 2, Figure 6). Finally, centroid sequences associated with each marker can be exported to fasta files (menu 3, Figure 7). The above represents the complete cycle of data processing. Additionally (menu 4, Figure 8), mbctools software offers the possibility to export results in a format that is compatible with the metaXplor platform (Sempéré et al., 2021). The entire process is described as supplementary material.
If, at any point, users wish to exit the program and conduct further analyses of their choice, `mbctools` allows it since it generates intermediate output files at each step.

```
--- MENU 1: BASIC ANALYSIS - only option 1 is strictly mandatory ---
1 -> INITIAL ANALYSIS (mandatory): read merging, sample-level dereplication, sequence clustering, chimera detection, affiliation of sequences to loci, and sequence re-orientation
1a -> Re-analyze all loci, from the clustering step, modifying parameters
1b -> Re-analyze only one locus of paired-end amplicon (merged reads), modifying parameters
1c -> Re-analyze only one locus of single-end amplicon (R1 only), modifying parameters
1d -> Re-analyse a given sample, modifying parameters
1e -> Optional quality checking of fastq files (slow)
```

**Figure 5** - Contents of section 1 menu, in which only option 1, INITIAL ANALYSIS, is mandatory as it executes all steps required to proceed with menu 2. Options 1a to 1d may be used after option 1 in case the user wants to refine some parameters. Analysis 1e is optional and provides information about the reads' quality.

When option 1 (initial analysis) is selected from menu 1, the menu shown on Figure 5 allows to run the initial analysis (section 1), namely read merging, sample-level dereplication, sequence clustering (denoising), chimera detection, affiliation of sequences to markers, and sequence re-orientation. In the same menu, four other sections (1a, 1b, 1c, 1d) provide means to re-process some particular markers and/or samples with adjusted parameters. In most cases the default values for those parameters will work, but there might be situations where they need to be adapted, either globally (e.g., depending on the sequencing depth, the level of similarity between targeted genetic markers), at the marker (i.e., locus) level (e.g., according to its size), or at the sample level (e.g., to account for heterogeneous sample quality). Option 1e may be used to assess the read quality for each sample, and write results to files with suffix "*.quality.txt" in the outputs directory.

When the initial analysis is completed, menu 2 should be used to remove primer sequences, and to exclude denoised sequences with an abundance considered too low for being relevant (Figure 6). This operation may be launched either in batch mode, or sample by sample in case different thresholds need to be applied.

```
--- MENU 2: PRIMER REMOVAL, SELECTION OF MINIMUM SEQUENCE ABUNDANCE LEVELS ACCORDING TO USER-DEFINED THRESHOLDS ---
2a -> Apply the SAME size threshold for ALL SAMPLES for the loci based on PAIRED-END reads (R1/R2 merged)
    i.e. you want to keep only sequences whose abundance is greater than x% of the total number of
    sequences for a given sample. This threshold of x% can be chosen for each locus.
2b -> Apply the SAME size threshold for ALL SAMPLES for the loci based on SINGLE-END reads (R1 only)
    same as option 2a but only using the R1 reads instead of merged ones.
2c -> Apply a SPECIFIC size threshold for a given sample, for the loci based on PAIRED-END reads (R1/R2 merged)
    i.e. you want to modulate the threshold of x% by locus but also by sample within a particular locus.
2d -> Apply a SPECIFIC size threshold for a given sample, for the loci based on SINGLE-END reads (R1 only)
    same as option 2c but only using the R1 sequences instead of merged ones.
```

**Figure 6** - Menu 2 combines two procedures that need to be launched successively for each marker: primer removal and application of an abundance threshold for retaining centroid sequences.

In menu 3, the program groups sample files for each marker, creating a new file per locus with the suffix "\*\_allseq_select.fasta" (Figure 7). During this process, `mbctools` can, as an option, dereplicate sequences at the project level, ensuring that each entry in the fasta file is unique. In such cases, an additional tabulated file is generated to maintain a record of sequence abundance across the samples. Optionally generated files are then named with the following suffixes: "\*\_allseq_select_derep.fasta", "\*\_allseq_select_derep.tsv".
```
--- MENU 3: GENERATION OF A UNIQUE SEQUENCE FILE FOR EACH LOCUS (comprising all samples' data) ---
3a -> Process all loci at once
3b -> Process loci one by one
```

**Figure 7** - Contents of section 3 menu, allowing to obtain per-marker (per-locus) fasta files. Menu 3 provides means to generate fasta files containing the retained sequences, either for all markers (loci) at once, or for a given marker.

Finally, menu 4 converts *mbctools* outputs into a format suitable for importing into metaXplor (Figure 8), a web-oriented platform for storing, sharing, exploring and manipulating metagenomic pipeline outputs. metaXplor significantly contributes to enhancing the FAIR (Findable, Accessible, Interoperable, and Reusable) (Wilkinson et al., 2016) characteristics of this specific data type by serving as a robust solution for centralized data management, intuitive visualization, and long-term accessibility.

```
--- MENU 4: EXPORTING ANALYSIS RESULTS INTO metaXplor FORMAT ---
4a -> Generate sequence files
    Compiles all sequences selected for all loci into a single fasta
    Outputs a .tsv file indicating samples weights for each sequence
4b -> Generate assignment file
    Converts blastn results (obtained from blasting above-mentioned fasta file) from 'Hit table (text)'
    (format #7) into metaXplor format
4c -> Build metaXplor-format sample metadata file from provided tabulated file
4d -> Compresses all metaXplor files into a final, ready to import, zip archive
```

**Figure 8** - Contents of section 4 menu, guiding users in obtaining necessary files to create a metaXplor import archive (i.e., sample metadata, sequence fasta and abundance, assignment to reference taxa) in order to summarize all project results.

## Test data

The test data were generated in the scope of a study conducted on Chagas disease, a zoonotic parasitic illness caused by *Trypanosoma cruzi*. This disease is primarily transmitted to mammals through the feces and urine of hematophagous insects known as triatomines, commonly referred to as kissing bugs (Hemiptera: Reduviidae). *Triatoma dimidiata* is the main vector in the Yucatan Peninsula in Mexico. In 2019, the owner of an isolated rural dwelling collected thirty-three samples of *T. dimidiata* individuals from a conserved sylvatic site within the Balam Kú State Reserve, located in Campeche, Mexico (Figure 9). We used a metabarcoding approach to simultaneously characterize the genetic diversity of the vector and identify its blood meals, characterize its hindgut and midgut microbiome composition, and assess the genetic diversity of *T. cruzi* infecting *T. dimidiata* in this sylvatic ecotope. We extracted the total DNA from each *T. dimidiata* gut, and amplified different selected markers of interest for our question: (i) for *T. dimidiata* genotyping, a fragment of the vector's Internal Transcribed Spacer 2 (ITS-2) nuclear marker was amplified with primers ITS2_200F (5'-TCG YAT CTA GGC ATT GTC TG-3') and ITS2_200R (5'-CTC GCA GCT ACT AAG GGA ATC C-3') previously described in (Richards et al., 2013); (ii) for the identification of vertebrate blood meals, a fragment of the vertebrate 12S rRNA gene was amplified with primers L1085 (5'- CCC AAA CTG GGA TTA GAT ACC C-3') and H1259 (5'-GTT TGC TGA AGA TGG CGG TA-3') (Kitano et al., 2007; Moo-Millan et al., 2019); (iii) for the identification of hindgut and midgut bacterial microbiome composition, a fragment of the bacterial 16S rRNA gene was amplified with primers E786F (5'-GAT TAG ATA CCC TGG TAG-3') and U926R (5'-CCG TCA ATT CCT TTR AGT TT-3') (Baker et al. 2003); and (iv) to accurately genotype *T. cruzi*, we targeted markers COII, GPI, ND1 and the intergenic region of the mini-exon gene, using the primers described in (10). Then we pooled the different amplicons obtained per *T. dimidiata* specimen and sequenced them on an Illumina MiSeq. On average, we obtained 222,244 reads per sample. The reads were processed with the *mbctools* package, whose outputs were imported into metaXplor.
In *mbctools*, we retained the reads with a minimum length of 100bp, clusters with a minimum size of 8 and used an alpha value of 2. The taxonomic identification of the ZOTUs was performed using BLASTn (Johnson et al., 2008). The exported CSV file with the metadata served as input for metaXplor. The analysis results are directly available on the web application under the dataset 'mbctools_review' ("metaXplor: mbctools data set," 2024) (Figure 9).

<figure><img src="image_8.png" /><figcaption>**Figure 9** - Graphical output from the metaXplor web interface. On the left side the geographic origin of the thirty three samples. On the right side the taxonomic proportion of each identified organism</figcaption></figure>

We identified an average of 33 ZOTUs per *T. dimidiata* sample, and only 9.61% of the reads were discarded during the analysis (Table 1). After analysis, we found thirteen blood-meal sources (the most frequently identified being *Homo sapiens*, *Sus scrofa* and *Mus musculus*), six bacterial Orders (Bacillales, Oceanospirillales and Micrococcales being the most abundant), along with *T. cruzi* and *T. dimidiata* (Figure 9). The test data input and output files are freely accessible on the DataSuds repository (Manzanilla et al., 2024).

<table><caption>**Table 1** - Statistics concerning the processing of the test data set with the number of reads per sample and the attribution to each marker.</caption><thead><tr><th rowspan="2">Sample</th><th rowspan="2">Reads per sample</th><th rowspan="2">Merging</th><th rowspan="2">Number of clusters</th><th rowspan="2">clusters with no chimera</th><th colspan="2">Cluster per marker</th><th rowspan="2">COII</th><th rowspan="2">GPI</th><th rowspan="2">ITS2</th><th rowspan="2">ND1</th><th rowspan="2">MINEX</th><th rowspan="2">MINEX_SE</th></tr><tr><th>12S</th><th>16S</th></tr></thead><tbody><tr><td>CON001</td><td>215137</td><td>0.8696</td><td>11</td><td>11</td><td>0</td><td>15</td><td>0</td><td>0</td><td>2</td><td>3</td><td>6</td><td>0</td></tr><tr><td>CON002</td><td>268415</td><td>0.9229</td><td>17</td><td>13</td><td>18</td><td>12</td><td>0</td><td>0</td><td>2</td><td>0</td><td>0</td><td>0</td></tr><tr><td>CON003</td><td>242823</td><td>0.8636</td><td>87</td><td>87</td><td>0</td><td>9</td><td>0</td><td>1</td><td>2</td><td>0</td><td>12</td><td>0</td></tr><tr><td>CON004</td><td>205358</td><td>0.9344</td><td>3</td><td>3</td><td>0</td><td>11</td><td>0</td><td>0</td><td>3</td><td>0</td><td>0</td><td>0</td></tr><tr><td>CON005</td><td>149945</td><td>0.9408</td><td>7</td><td>7</td><td>1</td><td>6</td><td>0</td><td>0</td><td>2</td><td>0</td><td>0</td><td>0</td></tr><tr><td>CON006</td><td>180763</td><td>0.9311</td><td>2</td><td>2</td><td>0</td><td>11</td><td>0</td><td>0</td><td>3</td><td>0</td><td>0</td><td>0</td></tr><tr><td>CON007</td><td>7898</td><td>0.891</td><td>2</td><td>2</td><td>1</td><td>0</td><td>0</td><td>0</td><td>2</td><td>0</td><td>0</td><td>0</td></tr><tr><td>CON008</td><td>200589</td><td>0.924</td><td>3</td><td>3</td><td>0</td><td>19</td><td>0</td><td>0</td><td>2</td><td>0</td><td>0</td><td>0</td></tr><tr><td>CON009</td><td>261914</td><td>0.9226</td><td>2</td><td>2</td><td>0</td><td>15</td><td>0</td><td>0</td><td>2</td><td>0</td><td>0</td><td>0</td></tr><tr><td>CON010</td><td>283721</td><td>0.9325</td><td>139</td><td>126</td><td>16</td><td>16</td><td>0</td><td>0</td><td>2</td><td>0</td><td>0</td><td>0</td></tr><tr><td>CON011</td><td>186936</td><td>0.9112</td><td>1</td><td>1</td><td>0</td><td>16</td><td>0</td><td>0</td><td>2</td><td>0</td><td>0</td><td>0</td></tr><tr><td>CON012</td><td>194776</td><td>0.8326</td><td>61</td><td>60</td><td>14</td><td>10</td><td>2</td><td>1</td><td>2</td><td>0</td><td>24</td><td>4</td></tr><tr><td>CON013</td><td>189441</td><td>0.9303</td><td>14</td><td>14</td><td>4</td><td>8</td><td>0</td><td>0</td><td>4</td><td>0</td><td>5</td><td>0</td></tr><tr><td>CON014</td><td>151880</td><td>0.9086</td><td>25</td><td>25</td><td>0</td><td>8</td><td>0</td><td>0</td><td>30</td><td>0</td><td>0</td><td>0</td></tr><tr><td>CON015</td><td>221302</td><td>0.9174</td><td>1</td><td>1</td><td>0</td><td>14</td><td>0</td><td>0</td><td>2</td><td>0</td><td>0</td><td>0</td></tr><tr><td>CON016</td><td>139546</td><td>0.8699</td><td>46</td><td>46</td><td>2</td><td>6</td><td>2</td><td>1</td><td>10</td><td>0</td><td>14</td><td>3</td></tr><tr><td>CON017</td><td>228038</td><td>0.8812</td><td>3</td><td>3</td><td>0</td><td>11</td><td>0</td><td>0</td><td>2</td><td>0</td><td>0</td><td>0</td></tr><tr><td>CON018</td><td>258481</td><td>0.9012</td><td>28</td><td>27</td><td>6</td><td>12</td><td>0</td><td>0</td><td>15</td><td>0</td><td>0</td><td>0</td></tr><tr><td>CON019</td><td>251403</td><td>0.7904</td><td>77</td><td>77</td><td>1</td><td>14</td><td>4</td><td>2</td><td>10</td><td>0</td><td>3</td><td>0</td></tr><tr><td>CON020</td><td>223972</td><td>0.9283</td><td>24</td><td>24</td><td>2</td><td>1</td><td>0</td><td>0</td><td>4</td><td>0</td><td>0</td><td>0</td></tr><tr><td>CON021</td><td>309618</td><td>0.9346</td><td>8</td><td>8</td><td>0</td><td>15</td><td>0</td><td>0</td><td>14</td><td>0</td><td>0</td><td>0</td></tr><tr><td>CON022</td><td>197125</td><td>0.8479</td><td>83</td><td>83</td><td>11</td><td>8</td><td>2</td><td>2</td><td>4</td><td>0</td><td>5</td><td>0</td></tr><tr><td>CON023</td><td>193857</td><td>0.7558</td><td>53</td><td>53</td><td>1</td><td>7</td><td>2</td><td>1</td><td>5</td><td>5</td><td>15</td><td>0</td></tr><tr><td>CON024</td><td>190310</td><td>0.9132</td><td>62</td><td>62</td><td>13</td><td>20</td><td>0</td><td>0</td><td>2</td><td>0</td><td>0</td><td>0</td></tr><tr><td>CON025</td><td>273603</td><td>0.9104</td><td>44</td><td>37</td><td>29</td><td>15</td><td>0</td><td>0</td><td>2</td><td>0</td><td>0</td><td>0</td></tr><tr><td>CON026</td><td>284204</td><td>0.945</td><td>2</td><td>2</td><td>1</td><td>17</td><td>0</td><td>0</td><td>2</td><td>0</td><td>0</td><td>0</td></tr><tr><td>CON027</td><td>436822</td><td>0.9225</td><td>5</td><td>5</td><td>1</td><td>34</td><td>0</td><td>0</td><td>2</td><td>0</td><td>0</td><td>0</td></tr><tr><td>CON028</td><td>243235</td><td>0.9385</td><td>63</td><td>63</td><td>1</td><td>12</td><td>0</td><td>0</td><td>4</td><td>0</td><td>0</td><td>0</td></tr><tr><td>CON029</td><td>219907</td><td>0.9486</td><td>13</td><td>13</td><td>4</td><td>11</td><td>0</td><td>0</td><td>2</td><td>0</td><td>0</td><td>0</td></tr><tr><td>CON030</td><td>306738</td><td>0.9456</td><td>17</td><td>17</td><td>0</td><td>10</td><td>0</td><td>0</td><td>26</td><td>0</td><td>0</td><td>0</td></tr><tr><td>CON031</td><td>172815</td><td>0.9398</td><td>87</td><td>87</td><td>1</td><td>11</td><td>0</td><td>0</td><td>2</td><td>0</td><td>0</td><td>0</td></tr><tr><td>CON032</td><td>206078</td><td>0.9069</td><td>34</td><td>33</td><td>1</td><td>13</td><td>0</td><td>0</td><td>24</td><td>0</td><td>0</td><td>0</td></tr><tr><td>CON033</td><td>237429</td><td>0.9158</td><td>92</td><td>92</td><td>37</td><td>14</td><td>0</td><td>0</td><td>2</td><td>0</td><td>0</td><td>0</td></tr><tr><td>Average</td><td>222244.8</td><td>0.9</td><td>33.8</td><td>33</td><td>5</td><td>12.2</td><td>0.36</td><td>0.24</td><td>5.88</td><td>0.24</td><td>2.55</td><td>0.21</td></tr></tbody></table>
Conclusion

mbctools facilitates the processing of amplicon datasets for metabarcoding studies by offering menus that guide users through the various steps of the data analysis. It is, to our knowledge, the only tool eliminating the need for the computer literacy normally required to utilize command-line based software as complex as the underlying VSEARCH, while allowing to process multiple genetic markers simultaneously. mbctools offers the possibility to reprocess specific markers and samples with customized parameters, providing means to adapt to the diversity of metabarcoding datasets. In the end, it generates fasta files composed of the filtered ZOTUs ready for the taxonomic assignment step. Although BLASTn is a commonly used and straightforward solution for this task, given the diversity of tools and approaches (k-mer based Kraken2 (Wood et al., 2019), Kaiju (Menzel et al. 2016), phylogeny-based PhyloSift (Darling et al., 2014)), we leave to the user to proceed with taxonomic assignment, based on his own practice or the literature. To promote data management, visualization and long-term accessibility, mbctools offers a file conversion functionality compatible with the metaXplor web application which aims at centralizing online meta-omic data, while offering user-friendly means to interact with it. The metaXplor instance hosted at CIRAD (https://metaxplor.cirad.fr/) is indeed used by several teams to keep track of previous project data (be it private or public) and proves useful in helping scientists quickly recover precise information. As an example, metaXplor was utilized in the scope of 24 shotgun metagenomics projects based on VANA (Virion-Associated Nucleic Acids), conducted to uncover viral diversity (Moubset et al., 2022). It allowed for the incremental building of a data repository that efficiently manages and provides means to analyze the extensive sequence datasets thus generated. This platform facilitated similarity-based searches and phylogenetic analyses, significantly enhancing the retrieval and re-analysis of data, thereby promoting viral discovery and classification. Although we initially developed the mbctools pipeline and software for reproducibility purposes across our teams working on infectious diseases, the pipeline's strength resides in its capacity to handle diverse data and be tailored to different research needs. This tool plays a central role in training sessions provided by our teams in developing countries, and its functionalities and user-friendliness are recurrently extended according to the feedback they generate. Depending on the success of this first release, future versions may add support for new sequencing technologies, embedded taxonomic assignment solutions or phylogenetic tree building.

Authors' contributions

CB and GS designed and conceived the software based on EW's specifications and JMM's
feedbacks. VM did the Python packaging. CB, GS, VM, JMM, AAR, and EW wrote the article. AAR
provided the data to test mbctools. All authors read and approved the final manuscript.

Funding

This work received financial support from CONACYT (National Council of Science and Technology, Mexico) Basic Science (Project ID: CB2015-258752) and National Problems (Project ID: PN2015-893) grants awarded to EW, as well as from IRD (French National Research Institute for Sustainable Development).

Data, scripts, code, and supplementary information availability

mbctools is openly accessible on Github at https://github.com/GuilhemSempere/mbctools and on PyPI at https://pypi.org/project/mbctools/. The main pipeline description is available as supplementary material, as are higher-definition versions of Figures 1 to 3.
Acknowledgements

The authors would like to thank Frédéric Mahé for taking the time to introduce VSEARCH functionalities, as well as Dorian Grasset for testing the OSX version and creating a GitHub action to automatize PyPI package generation. Preprint version 2 of this article has been peer-reviewed and recommended by Peer Community In Genomics (https://doi.org/10.24072/pci.genomics.100370; Pollet, 2024).

Conflict of interest disclosure

The authors declare that they comply with the PCI rule of having no financial conflicts of interest in relation to the content of the article.

References

Bailet B, Apothéloz-Perret-Gentil L, Baričević A, Chonova T, Franc A, Frigerio J-M, Kelly M, Mora D, Pfannkuchen M, Proft S. 2020. Diatom DNA metabarcoding for ecological assessment: Comparison among bioinformatics pipelines used in six European countries reveals the need for standardization. *Science of the Total Environment* 745:140948.
https://doi.org/10.1016/j.scitotenv.2020.140948

Baker GC, Smith JJ, Cowan DA. 2003. Review and re-analysis of domain-specific 16S primers. *Journal of Microbiological Methods* 55:541-555. https://doi.org/10.1016/j.mimet.2003.08.009

Barnabé C, Brenière SF, Santillán-Guayasamín S, Douzery EJP, Waleckx E. 2023. Revisiting gene typing and phylogeny of *Trypanosoma cruzi* reference strains: Comparison of the relevance of mitochondrial DNA, single-copy nuclear DNA, and the intergenic region of mini-exon gene. *Infection, Genetics and Evolution* 115:105504.
https://doi.org/10.1016/j.meegid.2023.105504

Caporaso JG, Lauber CL, Walters WA, Berg-Lyons D, Lozupone CA, Turnbaugh PJ, Fierer N, Knight R. 2011. Global patterns of 16S rRNA diversity at a depth of millions of sequences per sample. *Proceedings of the National Academy of Sciences* 108:4516-4522.
https://doi.org/10.1073/pnas.1000080107

CBOL Plant Working Group, Hollingsworth PM, Forrest LL, Spouge JL, Hajibabaei M, Ratnasingham S, van der Bank M, Chase MW, Cowan RS, Erickson DL. 2009. A DNA barcode for land plants. *Proceedings of the National Academy of Sciences* 106:12794-12797.
https://doi.org/10.1073/pnas.0905845106

Darling AE, Jospin G, Lowe E, Matsen IV FA, Bik HM, Eisen JA. 2014. PhyloSift: phylogenetic analysis of genomes and metagenomes. *PeerJ* 2:e243. https://doi.org/10.7717/peerj.243

Dumontel E, Ramirez-Sierra M-J, Pérez-Carrillo S, Teh-Poot C, Herrera C, Gourbière S, Waleckx E. 2018. Detailed ecological associations of triatomines revealed by metabarcoding and next-generation sequencing: implications for triatomine behavior and *Trypanosoma cruzi* transmission cycles. *Scientific Reports* 8:4140. https://doi.org/10.1038/s41598-018-22455-x

Edgar RC. 2016. UNOISE2: improved error-correction for Illumina 16S and ITS amplicon sequencing. *bioRxiv*, 081257. https://doi.org/10.1101/081257

Hakimzadeh A, Asbun AA, Albanese D, Bernard M, Buchner D, Callahan B, Caporaso G, Curd E, Djemiel C, Durling MB et al. 2023. A pile of pipelines: an overview of the bioinformatics software for metabarcoding data analyses. *Molecular Ecology Resources* 24:e13847.
https://doi.org/10.1111/1755-0998.13847

Hebert PDN, Cywinska A, Ball SL, deWaard JR. 2003. Biological identifications through DNA barcodes. *Proceedings of the Royal Society of London. Series B: Biological Sciences* 270:313-321. https://doi.org/10.1098/rspb.2002.2218

Hebert PD, Ratnasingham S, Zakharov EV, Telfer AC, Levesque-Beaudin V, Milton MA, Pedersen S, Jannetta P, DeWaard JR. 2016. Counting animal species with DNA barcodes: Canadian
insects. *Philosophical Transactions of the Royal Society B: Biological Sciences* 371:20150333.
https://doi.org/10.1098/rstb.2015.0333

Hernández-Andrade A, Moo-Millan J, Cigarroa-Toledo N, Ramos-Ligonio A, Herrera C, Bucheton B, Bart J-M, Jamonneau V, Bañuls A-L, Paupy C. 2019. Metabarcoding: a powerful yet still underestimated approach for the comprehensive study of vector-borne pathogen transmission cycles and their dynamics. In Claborn D, Bhattacharya S, Roy S. (eds), Vector-Borne Diseases - Recent Developments in Epidemiology and Control. Intech Open, London, UK. ISBN: 978-1-83880-022-2. https://doi.org/10.5772/intechopen.89839

Holdaway RJ, Wood JR, Dickie IA, Orwin KH, Bellingham PJ, Richardson SJ, Lyver PO, Timoti P, Buckley TR. 2017. Using DNA metabarcoding to assess New Zealand's terrestrial biodiversity. *New Zealand Journal of Ecology* 41:251–262. https://doi.org/10.20417/nzjecol.41.28

López-Dominguez J, López-Monteón A, Ochoa-Martínez P, Dumonteil E, Barnabé C, Waleckx E, Hernández-Giles RG, Ramos-Ligonio A. 2022. Molecular characterization of four Mexican isolates of *Trypanosoma cruzi* and their profile susceptibility to nifurtimox. *Acta parasitologica* 67:1584–1593. https://doi.org/10.1007/s11686-022-00608-3

Johnson M, Zaretskaya I, Raytselis Y, Merezhuk Y, McGinnis S, Madden TL. 2008. NCBI BLAST: a better web interface. *Nucleic acids research* 36:W5–W9. https://doi.org/10.1093/nar/gkn201

Kitano T, Umetsu K, Tian W, Osawa M. 2007. Two universal primer sets for species identification among vertebrates. *International Journal of Legal Medicine* 121:423–427. https://doi.org/10.1007/s00414-006-0113-y

Lauthier JJ, Tomasini N, Barnabé C, Rumi MMM, D'Amato AMA, Ragone PG, Yeo M, Lewis MD, Llewellyn MS, Basombrío MA. 2012. Candidate targets for Multilocus Sequence Typing of *Trypanosoma cruzi*: validation using parasite stocks from the Chaco Region and a set of reference strains. *Infection, Genetics and Evolution* 12:350–358. https://doi.org/10.1016/j.meegid.2011.12.008

Leray M, Knowlton N. 2015. DNA barcoding and metabarcoding of standardized samples reveal patterns of marine benthic diversity. *Proceedings of the National Academy of Sciences* 112:2076–2081. https://doi.org/10.1073/pnas.1424997112

Manzanilla V, Barnabé C, Guilhem S, Millan JM, Amblard-Rambert A, Waleckx E. 2024. Test data for mbctools. Available at https://doi.org/10.23708/W3TODJ (accessed November 18th, 2024).

Mathon L, Valentini A, Guérin P, Normandeau E, Noel C, Lionnet C, Boulanger E, Thuiller W, Bernatchez L, Mouillot D. 2021. Benchmarking bioinformatic tools for fast and accurate eDNA metabarcoding species identification. *Molecular Ecology Resources* 21:2565–2579. https://doi.org/10.1111/1755-0998.13430

Menzel P, Ng KL, Krogh A. 2016. Fast and sensitive taxonomic classification for metagenomics with Kaiju. *Nature communications* 7:11257. https://doi.org/10.1038/ncomms11257

metaXplor: mbctools data set. 2024. Available at https://metaxplor.cirad.fr/metaXplor/main.jsp?module=mbctools_review (accessed February 5, 2024).

Moo-Millan JI, Arnal A, Pérez-Carrillo S, Hernandez-Andrade A, Ramírez-Sierra M-J, Rosado-Vallado M, Dumonteil E, Waleckx E. 2019. Disentangling *Trypanosoma cruzi* transmission cycle dynamics through the identification of blood meal sources of natural populations of *Triatoma dimidiata* in Yucatán, Mexico. *Parasites & Vectors* 12:572. https://doi.org/10.1186/s13071-019-3819-7

Moubset O, François S, Maclot F, Palanga E, Julian C, Claude L, Fernandez E, Rott P, Daugrois J-H, Antoine-Lorquin A, Bernardo P, Blouin AG, Temple C, Kraberger S, Fontenele RS, Harkins GW, Ma Y, Marais A, Candresse T, Chéhida SB, Lefeuvre P, Lett J-M, Varsani A, Massart S, Ogliastro M, Martin DP, Filloux D, Roumagnac P. 2022. Virion-Associated Nucleic Acid-Based Metagenomics: A Decade of Advances in Molecular Characterization of Plant Viruses. *Phytopathology* 112:2253–2272. https://doi.org/10.1094/PHYTO-03-22-0096-RVW

Pollet N. 2024. One tool to metabarcode them all. Peer Community Genomics 100370. https://doi.org/10.24072/pci.genomics.100370
Raclariu AC, Paltinean R, Vlase L, Labarre A, Manzanilla V, Ichim MC, Crisan G, Brysting AK, De Boer H. 2017. Comparative authentication of *Hypericum perforatum* herbal products using DNA metabarcoding, TLC and HPLC-MS. *Scientific Reports* 7:1291. https://doi.org/10.1038/s41598-017-01389-w

Riaz T, Shehzad W, Viari A, Pompanon F, Taberlet P, Coissac E. 2011. ecoPrimers: inference of new DNA barcode markers from whole genome sequence analysis. *Nucleic Acids Research* 39:e145–e145. https://doi.org/10.1093/nar/gkr732

Richards B, De La Rua NM, Monroy C, Stevens L, Dorn PL. 2013. Novel polymerase chain reaction-restriction fragment length polymorphism assay to determine internal transcribed spacer-2 group in the Chagas disease vector, *Triatoma dimidiata* (Latreille, 1811). *Memórias do Instituto Oswaldo Cruz* 108:395–398. https://doi.org/10.1590/s0074-0276108042013001

Rognes T, Flouri T, Nichols B, Quince C, Mahé F. 2016. VSEARCH: a versatile open source tool for metagenomics. *PeerJ* 4:e2584. https://doi.org/10.7717/peerj.2584

Schoch CL, Seifert KA, Huhndorf S, Robert V, Spouge JL, Levesque CA, Chen W, Fungal Barcoding Consortium, et al. 2012. Nuclear ribosomal internal transcribed spacer (ITS) region as a universal DNA barcode marker for Fungi. *Proceedings of the National Academy of Sciences* 109:6241–6246. https://doi.org/10.1073/pnas.1117018109

Sempéré G, Pétel A, Abbé M, Lefeuvre P, Roumagnac P, Mahé F, Baurens G, Filloux D. 2021. metaXplor: an interactive viral and microbial metagenomic data manager. *GigaScience* 10:giab001. https://doi.org/10.1093/gigascience/giab001

Taberlet P, Coissac E, Pompanon F, Brochmann C, Willerslev E. 2012. Towards next-generation biodiversity assessment using DNA metabarcoding. *Molecular Ecology* 21:2045–2050. https://doi.org/10.1111/j.1365-294x.2012.05470.x

Thomsen PF, Willerslev E. 2015. Environmental DNA–An emerging tool in conservation for monitoring past and present biodiversity. *Biological conservation* 183:4–18. https://doi.org/10.1016/j.biocon.2014.11.019

Wilkinson MD, Dumontier M, Aalbersberg IJJ, Appleton G, Axton M, Baak A, Blomberg N, Boiten J-W, da Silva Santos LB, Bourne PE. 2016. The FAIR Guiding Principles for scientific data management and stewardship. *Scientific data* 3:1–9. https://doi.org/10.1038/sdata.2016.18

Wood DE, Lu J, Langmead B. 2019. Improved metagenomic analysis with Kraken 2. *Genome biology* 20:1–13. https://doi.org/10.1186/s13059-019-1891-0