Abstract:Background: High throughput sequencing of environmental DNA has applications in biodiversity monitoring, taxa abundance estimation, understanding the dynamics of community ecology, and marine species studies and conservation. Environmental DNA, especially, marine eDNA, has a fast degradation rate. Aside from the good quality reads, the data could have a significant number of reads that fall slightly below the default PHRED quality threshold of 30 on sequencing. For quality control, trimming methods are employed, which generally precede the merging of the read pairs. However, in the case of eDNA, a significant percentage of reads within the acceptable quality score range are also dropped. Methods: To infer the ideal merge tool that is sensitive to eDNA, two Hiseq paired-end eDNA datasets were utilized to study the merging by the tools - FLASH (Fast Length Adjustment of SHort reads), PANDAseq, COPE, BBMerge, and VSEARCH without preprocessing. We assessed these tools on the following parameters: Time taken to process, the quality, and the number of merged reads. Trimmomatic, a widely-used preprocessing tool, was also assessed by preprocessing the datasets at different parameters for the two approaches of preprocessing: Sliding Window and Maximum Information. The preprocessed read pairs were then merged using the ideal merge tool identified earlier. Results: FLASH is the most efficient merge tool balancing data conservation, quality of reads, and processing time. We compared Trimmomatic's two quality trimming options with increasing strictness with FLASH's direct merge. The raw reads processed with Trimmomatic then merged, yielded a significant drop in reads compared to the direct merge. An average of 29% of reads was dropped when directly merged with FLASH. Maximum Information option resulted in 30.7% to 68.05% read loss with lowest and highest stringency parameters, respectively. The Sliding Window approach conserves approximately 10% more reads at a PHRED score of 25 set as the threshold for a window of size 4. The lowered PHRED cut off conserves about 50% of the reads that could potentially be informative. We noted no significant reduction of data while optimizing the number of reads read in a window with the ideal quality (Q) score. Conclusions: Losing reads can negatively impact the downstream processing of the environmental data, especially for sequence alignment studies. The quality trim-first-merge-later approach can significantly decrease the number of reads conserved. However, direct merging of pair-end reads using FLASH conserved more than 60% of the reads. Therefore, direct merging of the paired-end reads can prevent potential removal of informative reads that do not comply by the trimming tool's strict checks. FLASH to be an efficient tool in conserving reads while carrying out quality trimming in moderation. Overall, our results show that merging paired-end reads of eDNA data before trimming can conserve more reads.

MeFiT: merging and filtering tool for illumina paired-end reads for 16S rRNA amplicon sequencing

MeFiT: Merging and Filtering Tool for Illumina Paired-End Reads for 16S rRNA Amplicon Sequencing

MEEPTOOLS: A maximum expected error based FASTQ read filtering and trimming toolkit

Fast and Simple Analysis of MiSeq Amplicon Sequencing Data with MetaAmp.

Merfin: improved variant filtering and polishing via k-mer validation

High-quality single amplicon sequencing method for illumina MiSeq platform using pool of 'N' (0-10) spacer-linked target specific primers without PhiX spike-in

PEAR: a fast and accurate Illumina Paired-End reAd mergeR

SimpleMetaPipeline: Breaking the bioinformatics bottleneck in metabarcoding

Benefits of merging paired-end reads before pre-processing environmental metagenomics data

A read-filtering algorithm for high-throughput marker-gene studies that greatly improves OTU accuracy

Merging and concatenation of sequencing reads: a bioinformatics workflow for the comprehensive profiling of microbiome from amplicon data

Mipe: A Metagenome-Based Community Structure Explorer and Ssu Primer Evaluation Tool

16S-FASAS: an integrated pipeline for synthetic full-length 16S rRNA gene sequencing data analysis

From reads to operational taxonomic units: an ensemble processing pipeline for MiSeq amplicon sequencing data

MeRIPseqPipe: an integrated analysis pipeline for MeRIP-seq data based on Nextflow

Microbe Finder (MiFi®): Implementation of an Interactive Pathogen Detection Tool in Metagenomic Sequence Data

FIGARO: An efficient and objective tool for optimizing microbiome rRNA gene trimming parameters

An improved dual-indexing approach for multiplexed 16S rRNA gene sequencing on the Illumina MiSeq platform

MitoFinder: Efficient automated large‐scale extraction of mitogenomic data in target enrichment phylogenomics

BiomiX, a User-Friendly Bioinformatic Tool for Automatized Multiomics Data Analysis and Integration

Merfin: improved variant filtering, assembly evaluation and polishing via k-mer validation