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

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.

Effects of source sample amount on biodiversity surveys of bacteria, fungi, and nematodes in soil ecosystems

Bacteria, fungi, and nematodes are major components of soil ecosystems, playing pivotal roles in belowground material cycles and biological community processes. A number of studies have recently uncovered the diversity and community structure of those organisms in various types of soil ecosystems based on DNA metabarcoding (amplicon sequencing). However, because most previous studies examined only one or two of the three organismal groups, it remains an important challenge to reveal the entire picture of soil community structure. We examined how we could standardize DNA extraction protocols for simultaneous DNA metabarcoding of bacteria, fungi, and nematodes. Specifically, in an Illumina sequencing analysis of forest and farmland soil samples, we performed DNA extraction at five levels of soil-amount (0.5, 2, 5, 10, and 20 g). We then found that DNA extraction with the 0.5 g soil setting, which had been applied as default in many commercial DNA extraction kits, could lead to underestimation of α-diversity in nematode community. We also found that dissimilarity (β-diversity) estimates of community structure among replicate samples could be affected by soil sample amount. Based on the assays, we conclude that DNA extraction from at least 20 g of soil is a standard for comparing biodiversity patterns among bacteria, fungi and nematodes.

Sensitive and accurate DNA metabarcoding of parasitic helminth mock communities using the mitochondrial rRNA genes

Next-generation sequencing technologies have accelerated the pace of helminth DNA metabarcoding research, enabling species detection in bulk community samples. However, finding suitable genetic markers with robust species-level resolution and primers targeting a broad species range among parasitic helminths are some of the challenges faced. This study aimed to demonstrate the potential use of the mitochondrial 12S and 16S rRNA genes for parasitic helminth (nematodes, trematodes, cestodes) DNA metabarcoding. To demonstrate the robustness of the 12S and 16S rRNA genes for DNA metabarcoding, we determined the proportion of species successfully recovered using mock helminth communities without environment matrix and mock helminth communities artificially spiked with environmental matrices. The environmental matrices are human fecal material, garden soil, tissue, and pond water. Our results revealed the robustness of the mitochondrial rRNA genes, through the high sensitivity of the 12S rRNA gene, and the effectiveness of the 12S and 16S primers targeting platyhelminths. With the mitochondrial rRNA genes, a broad range of parasitc helminths were successfully detected to the species level. The potential of the mitochondrial rRNA genes for helminth DNA metabarcoding was demonstrated, providing a valuable gateway for future helminth DNA metabarcoding applications like helminth detection and biodiversity studies.

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.