Taxonomic profiling of individual nematodes isolated from copse soils using deep amplicon sequencing of four distinct regions of the 18S ribosomal RNA gene.

Nematodes are representative soil metazoans with diverged species that play crucial roles in nutrient recycling in the pedosphere. Qualitative and quantitative information on nematode communities is useful for assessing soil quality, and DNA barcode-mediated taxonomic analysis is a powerful tool to investigate taxonomic compositions and changes in nematode communities. Here, we investigated four regions (regions 1-4) of the 18S small subunit ribosomal RNA (SSU) gene as PCR targets of deep amplicon sequencing for the taxonomic profiling of individual soil nematodes. We determined the sequence variants (SVs) of 4 SSU regions for 96 nematodes (total 384 amplicons) isolated from copse soils and assigned their taxonomy using the QIIME2 software with dada2 or deblur algorithm and the SILVA database. Dada2 detected approximately 2-fold more nematode-derived SVs than deblur, and a larger number of SVs were obtained in regions 1 and 4 than those in other regions. These results and sufficient reference sequence coverage in region 4 indicated that DNA barcoding using a primer set for region 4 followed by dada2-based analysis would be most suitable for soil nematode taxonomic analysis. Eighteen SSU-derived operational taxonomic units (rOTUs) were obtained from 68 isolates, and their orders were determined based on the phylogenetic trees built by four regional sequences of rOTUs and 116 nematode reference species as well as the BLASTN search. The majority of the isolates were derived from three major orders Dorylaimida (6 rOTUs, 51.5% in 68 isolates), Rhabditida (4 rOTUs, 29.4%), and Triplonchida (7 rOTUs, 17.6%). The predicted feeding types of the isolates were fungivores (38.2% in total nematodes), plant feeders (32.4%), and 14.7% for both bacterivores and omnivores/predators. Additionally, we attempted to improve the branch structure of phylogenetic trees by using long nucleotide sequences artificially prepared by connecting regional sequences, but the effect was limited.