Sampling cores and sequencing depths affected the measurement of microbial diversity in soil quadrats

Due to the massive quantity and broad phylogeny, an accurate measurement of microbial diversity is highly challenging in soil ecosystems. Initially, the deviation caused by sampling should be adequately considered. Here, we attempted to uncover the effect of different sampling strategies on α diversity measurement of soil prokaryotes. Four 1 m² sampling quadrats in a typical grassland were thoroughly surveyed through deep 16S rRNA gene sequencing (over 11 million reads per quadrat) with numerous replicates (33 soil sampling cores with total 141 replicates per quadrat). We found the difference in diversity was relatively small when pooling soil cores before and after DNA extraction and sequencing, but they were both superior to a non-pooling strategy. Pooling a small number of soil cores (i.e., 5 or 9) combined with several technical replicates is sufficient to estimate diversities for soil prokaryotes, and there is great flexibility in pooling original samples or data at different experimental steps. Additionally, the distribution of local α diversity varies with sampling core number, sequencing depth, and abundance distribution of the community, especially for high orders of Hill diversity index (i.e., Shannon entropy and inverse Simpson index). For each grassland soil quadrat (1 m²), retaining 100,000 reads after taxonomic clustering might be a realistic option, as these number of reads can efficiently cover the majority of common species in this area. Our findings provide important guidance for soil sampling strategy, and the general results can serve as a basis for further studies.