Abstract:Introduction: High-throughput sequencing (HTS) provides an efficient and cost-effective way to generate large amounts of sequence data, providing a very powerful tool to analyze biodiversity of soil organisms. However, marker-based methods and the resulting datasets come with a range of challenges and disputes, including incomplete reference databases, controversial sequence similarity thresholds for delimitating taxa, and downstream compositional data analysis. Methods: Here, we use HTS data from a soil nematode biodiversity experiment to explore standardized HTS data processing procedures. We compared the taxonomic assignment performance of two main rDNA reference databases (SILVA and PR2). We tested whether the same ecological patterns are detected with Amplicon Sequence Variants (ASV; 100% similarity) versus classical Operational Taxonomic Units (OTU; 97% similarity). Further, we tested how different HTS data normalization methods affect the recovery of beta diversity patterns and the identification of differentially abundant taxa. Results: At this time, the SILVA 138 eukaryotic database performed better than the PR2 4.12 database, assigning more reads to family level and providing higher phylogenetic resolution. ASV- and OTU-based alpha and beta diversity of nematodes correlated closely, indicating that OTU-based studies represent useful reference points. For downstream data analyses, our results indicate that loss of data during subsampling under rarefaction-based methods might reduce the sensitivity of the method, e.g. underestimate the differences between nematode communities under different treatments, while the clr-transformation-based methods may overestimate effects. The Analysis of Compositions of Microbiome with Bias Correction approach (ANCOM-BC) retains all data and accounts for uneven sampling fractions for each sample, suggesting that this is currently the optimal method to analyze compositional data. Discussion: Overall, our study highlights the importance of comparing and selecting taxonomic reference databases before data analyses, and provides solid evidence for the similarity and comparability between OTU- and ASV-based nematode studies. Further, the results highlight the potential weakness of rarefaction-based and clr-transformation-based methods. We recommend future studies use ASV and that both the taxonomic reference databases and normalization strategies are carefully tested and selected before analyzing the data.

De novo genomic analyses for non-model organisms: an evaluation of methods across a multi-species data set

De novo assembly of transcriptomes and differential gene expression analysis using short-read data from emerging model organisms – a brief guide

ExpressAnalyst: A unified platform for RNA-sequencing analysis in non-model species

Defining Loci in Restriction-Based Reduced Representation Genomic Data from Nonmodel Species: Sources of Bias and Diagnostics for Optimal Clustering

Using Illumina Next Generation Sequencing Technologies to Sequence Multigene Families in De Novo Species

Ocean to Tree: Leveraging Single-Molecule RNA-Seq to Repair Genome Gene Models and Improve Phylogenomic Analysis of Gene and Species Evolution

De novo clustering of long reads by gene from transcriptomics data

Sequencing breakthroughs for genomic ecology and evolutionary biology

Combining independent de novo assemblies optimizes the coding transcriptome for nonconventional model eukaryotic organisms

High-Throughput Sequencing and Advanced Computational Methods: Revolutionizing Genomic Data Analysis

A composite genome approach to identify phylogenetically informative data from next-generation sequencing

Optimising high-throughput sequencing data analysis, from gene database selection to the analysis of compositional data: a case study on tropical soil nematodes

Disentangling Cobionts and Contamination in Long-Read Genomic Data using Sequence Composition

Whole Animal Genome Sequencing: user-friendly, rapid, containerized pipelines for processing, variant discovery, and annotation of short-read whole genome sequencing data

Quest for Orthologs in the Era of Biodiversity Genomics

A De-Novo Genome Analysis Pipeline (DeNoGAP) for large-scale comparative prokaryotic genomics studies

Comparative analysis of amphibian genomes: An emerging resource for basic and applied research

Building better genome annotations across the tree of life

XSAnno: a framework for building ortholog models in cross-species transcriptome comparisons

OrthologAL: A Shiny application for quality-aware humanization of non-human pre-clinical high-dimensional gene expression data