Distinct prokaryotic and eukaryotic communities and networks in two agricultural fields of central Japan with different histories of maize-cabbage rotation.

Crop rotation is an important agricultural practice for homeostatic crop cultivation. Here, we applied high-throughput sequencing of ribosomal RNA gene amplicons to investigate soil biota in two fields of central Japan with different histories of maize-cabbage rotation. We identified 3086 eukaryotic and 17,069 prokaryotic sequence variants (SVs) from soil samples from two fields rotating two crops at three different growth stages. The eukaryotic and prokaryotic communities in the four sample groups of two crops and two fields were clearly distinguished using ß-diversity analysis. Redundancy analysis showed the relationships of the communities in the fields to pH and nutrient, humus, and/or water content. The complexity of eukaryotic and prokaryotic networks was apparently higher in the cabbage-cultivated soils than those in the maize-cultivated soils. The node SVs (nSVs) of the networks were mainly derived from two eukaryotic phyla: Ascomycota and Cercozoa, and four prokaryotic phyla: Pseudomonadota, Acidobacteriota, Actinomycetota, and Gemmatimonadota. The networks were complexed by cropping from maize to cabbage, suggesting the formation of a flexible network under crop rotation. Ten out of the 16 eukaryotic nSVs were specifically found in the cabbage-cultivated soils were derived from protists, indicating the potential contribution of protists to the formation of complex eukaryotic networks.

Distinct community structures of soil nematodes from three ecologically different sites revealed by high-throughput amplicon sequencing of four 18S ribosomal RNA gene regions.

Quantitative taxonomic compositions of nematode communities help to assess soil environments due to their rich abundance and various feeding habitats. DNA metabarcoding by the 18S ribosomal RNA gene (SSU) regions were preferentially used for analyses of soil nematode communities, but the optimal regions for high-throughput amplicon sequencing have not previously been well investigated. In this work, we performed Illumina-based amplicon sequencing of four SSU regions (regions 1-4) to identify suitable regions for nematode metabarcoding using the taxonomic structures of nematodes from uncultivated field, copse, and cultivated house garden soils. The fewest nematode-derived sequence variants (SVs) were detected in region 3, and the total nematode-derived SVs were comparable in regions 1 and 4. The relative abundances of reads in regions 1 and 4 were consistent in both orders and feeding groups with prior studies, thus suggesting that region 4 is a suitable target for the DNA barcoding of nematode communities. Distinct community structures of nematodes were detected in the taxon, feeding habitat, and life-history strategy of each sample; i.e., Dorylamida- and Rhabditida-derived plant feeders were most abundant in the copse soil, Rhabditida-derived bacteria feeders in the house garden soil, and Mononchida- and Dorylamida-derived omnivores and predators and Rhabditida-derived bacteria feeders in the field soil. Additionally, low- and high-colonizer-persister (cp) groups of nematodes dominated in the house garden and copse soils, respectively, whereas both groups were found in the field soil, suggesting bacteria-rich garden soil, undisturbed and plant-rich copse soil, and a transient status of nematode communities in the field soil. These results were also supported by the maturity indices of the three sampling sites. Finally, the influence of the primer tail sequences was demonstrated to be insignificant on amplification. These findings will be useful for DNA metabarcoding of soil nematode communities by amplicon sequencing.

Use of universal primers for the 18S ribosomal RNA gene and whole soil DNAs to reveal the taxonomic structures of soil nematodes by high-throughput amplicon sequencing.

Nematodes are abundant metazoans that play crucial roles in nutrient recycle in the pedosphere. Although high-throughput amplicon sequencing is a powerful tool for the taxonomic profiling of soil nematodes, polymerase chain reaction (PCR) primers for amplification of the 18S ribosomal RNA (SSU) gene and preparation of template DNAs have not been sufficiently evaluated. We investigated nematode community structure in copse soil using four nematode-specific (regions 1-4) and two universal (regions U1 and U2) primer sets for the SSU gene regions with two DNAs prepared from copse-derived mixed nematodes and whole soil. The major nematode-derived sequence variants (SVs) identified in each region was detected in both template DNAs. Order level taxonomy and feeding type of identified nematode-derived SVs were distantly related between the two DNA preparations, and the region U2 was closely related to region 4 in the non-metric multidimensional scaling (NMDS) based on Bray-Curtis dissimilarity. Thus, the universal primers for region U2 could be used to analyze soil nematode communities. We further applied this method to analyze the nematodes living in two sampling sites of a sweet potato-cultivated field, where the plants were differently growing. The structure of nematode-derived SVs from the two sites was distantly related in the principal coordinate analysis (PCoA) with weighted unifrac distances, suggesting their distinct soil environments. The resultant ecophysiological status of the nematode communities in the copse and field on the basis of feeding behavior and maturity indices was fairly consistent with those of the copse- and the cultivated house garden-derived nematodes in prior studies. These findings will be useful for the DNA metabarcoding of soil eukaryotes, including nematodes, using soil DNAs.

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.