Impacts of Sampling Design on Estimates of Microbial Community Diversity and Composition in Agricultural Soils

Association between bacterial community and cadmium distribution across Colombian cacao crops

Assessing the bacterial community composition across cacao crops is important to understand its potential role as a modulator of cadmium (Cd) translocation to plant tissues under field conditions; Cd mobility between soil and plants is a complex and multifactorial problem that cannot be captured only by experimentation. Although microbes have been shown to metabolize and drive the speciation of Cd under controlled conditions, regardless of the link between soil bacterial community (SBC) dynamics and Cd mobilization in the rhizosphere, only a few studies have addressed the relationship between soil bacterial community composition (SBCC) and Cd content in cacao seeds (Cdseed). Therefore, this study aimed to explore the association between SBCC and different factors influencing the distribution of Cd across cacao crop systems. This study comprised 225 samples collected across five farms, where we used an amplicon sequencing approach to characterize the bacterial community composition. The soil Cd concentration alone (Cdsoil) was a poor predictor of Cdseed. Still, we found that this relationship was more apparent when the variation within farms was controlled, suggesting a role of heterogeneity within farms in modulating Cd translocation and, thus, seed Cd content. Our results provide evidence of the link between soil bacterial communities and the distribution of Cd across Colombian cacao crops, and highlight the importance of incorporating fine-spatial-scale studies to advance the understanding of factors driving Cd uptake and accumulation in cacao plants.

IMPORTANCE

Cadmium (Cd) content in cacao crops is an issue that generates interest due to the commercialization of chocolate for human consumption. Several studies provided evidence about the non-biological factors involved in its translocation into the cacao plant. However, factors related to this process, including soil bacterial community composition (SBCC), still need to be addressed. It is well known that soil microbiome could impact compounds' chemical transformation, including Cd, on the field. Here, we found the first evidence of the link between soil bacterial community composition and Cd concentration in cacao soils and seeds. It highlights the importance of including the variation of bacterial communities to assess the factors driving the Cd translocation into cacao seeds. Moreover, the results highlight the relevance of the spatial heterogeneity within and across cacao farms, influencing the variability of Cd concentrations.

Soil Micro-eukaryotic Diversity Patterns Along Elevation Gradient Are Best Estimated by Increasing the Number of Elevation Steps Rather than Within Elevation Band Replication

The development of high-throughput sequencing (HTS) of environmental DNA (eDNA) has stimulated the study of soil microbial diversity patterns and drivers at all scales. However, given the heterogeneity of soils, a challenge is to define effective and efficient sampling protocols that allow sound comparison with other records, especially vegetation. In studies of elevational diversity pattern, a trade-off is choosing between replication within elevation bands vs. sampling more elevation bands. We addressed this question for soil protists along an elevation gradient on Mt. Asahi, Hokkaido, Japan. We compared two sampling approaches: (1) the replicate strategy (five replicates at six elevational bands, total = 30) and (2) the transect strategy (one sample in each of 16 different elevational bands). Despite a nearly twofold lower sampling effort, the transect strategy yielded congruent results compared to the replicate strategy for the estimation of elevational alpha diversity pattern: the regression coefficients between diversity indices and elevation did not differ between the two options. Furthermore, for a given total number of samples, gamma diversity estimated across the entire transect was higher when sampling more elevational bands as compared to replication from fewer elevational bands. Beta diversity (community composition turnover) was lower within a given elevational band than between adjacent bands and increased with elevation distance. In redundancy analyses, soil organic matter-related variable (the first principal component of soil organic matter, water content, total organic carbon, and nitrogen by whom were highly correlated) and elevation best explained elevational beta diversity pattern for both sampling approaches. Taken together, our results suggest that sampling a single plot per elevation band will be sufficient to obtain a good estimate of soil micro-eukaryotic diversity patterns along elevation gradients. This study demonstrated the effectiveness of the transect strategy in estimating diversity patterns along elevation gradients which is instructive for future environmental or even experimental studies. While not advocating for completely replacing replication-based sampling practices, it is important to note that both replicate and transect strategies have their merits and can be employed based on specific research goals and resource limitations.

Uncovering the biogeographic pattern of the widespread nematode-trapping fungi Arthrobotrys oligospora: watershed is the key

Studies of biogeographic patterns of fungi have long been behind those of plants and animals. The presence of worldwide species, the lack of systematic sampling design and adequate sampling effort, and the lack of research units are responsible for this status. This study investigates the biogeographical patterns of Arthrobotrys oligospora, the most widespread globally distributed nematode-trapping fungi (NTF), by stratified collecting and analyzing 2,250 samples from 228 sites in Yunnan Province, China. The A. oligospora was isolated, and 149 strains were subjected to ITS, TUB, TEF and RPB2 gene sequencing and multi-gene association phylogeographic analysis. The results show that at population level A. oligospora is randomly distributed throughout Yunnan Province and has no biogeographical distribution pattern. At the genetic level, the phylogenetic tree of A. oligospora diverges into five major evolutionary clades, with a low degree of gene flow between the five clades. However, the correlation between the phylogenetic diversity of A. oligospora and geographical factors was low. There was no clear pattern in the phylogenetic clades distribution of A. oligospora either without dividing the study unit or when the grid was used as the study unit. When watersheds were used as the study unit, 67.4%, 63.3%, 65.9%, 83.3%, and 66.7% of clade 1-5 strains were distributed in the Jinsha river, Red river, Peal river, Lancang river, and Nujiang-Irawaddy river watersheds, respectively. The clades distribution of A. oligospora was highly consistent with the watersheds distribution. Training predictions of the clades distributions using randomly generated polygons were also less accurate than watersheds. These results suggest that watersheds are key to discovering the biogeographic distribution patterns of A. oligospora. The A. oligospora populations are blocked by mountains in the watershed, and gene flow barriers have occurred, which may have resulted in the formation of multiple cryptic species. Watersheds are also ideal for understanding such speciation processes, explaining factors affecting biodiversity distribution and coupling studies of plant and animal and microbial diversity.

Those Nematode-Trapping Fungi That are not Everywhere: Hints Towards Soil Microbial Biogeography

The existence of biogeography for microorganisms is a raising topic in ecology and researchers are employing better distinctions between single species, including the most rare ones, to reveal potential hidden patterns. An important volume of evidence supporting heterogeneous distributions for bacteria, archaea and protists is accumulating, and more recently a few efforts have targeted microscopic fungi. We propose an insight into this latter kingdom by looking at a group of soil nematode-trapping fungi whose species are well-known and easily recognizable. We chose a pure culture approach because of its reliable isolation procedures for this specific group. After morphologically and molecularly identifying all species collected from 2250 samples distributed in 228 locations across Yunnan province of China, we analyzed occurrence frequencies and mapped species, genera, and richness. Results showed an apparent cosmopolitan tendency for this group of fungi, including species richness among sites. However, only four species were widespread across the region, while non-random heterogeneous distributions were observed for the remaining 40 species, both in terms of statistical distribution of species richness reflected by a significant variance-to-mean ratio, as well as in terms of visually discernible spatial clusters of rare species and genera on the map. Moreover, several species were restricted to only one location, raising the question of whether endemicity exists for this microbial group. Finally, environmental heterogeneity showed a marginal contribution in explaining restricted distributions, suggesting that other factors such as geographical isolation and dispersal capabilities should be explored. These findings contribute to our understanding of the cryptic geographic distribution of microorganisms and encourage further research in this direction.

Effects of operational taxonomic unit inference methods on soil microeukaryote community analysis using long-read metabarcoding

Long amplicon metabarcoding has opened the door for phylogenetic analysis of the largely unknown communities of microeukaryotes in soil. Here, we amplified and sequenced the ITS and LSU regions of the rDNA operon (around 1500 bp) from grassland soils using PacBio SMRT sequencing. We tested how three different methods for generation of operational taxonomic units (OTUs) effected estimated richness and identified taxa, and how well large-scale ecological patterns associated with shifting environmental conditions were recovered in data from the three methods. The field site at Kungsängen Nature Reserve has drawn frequent visitors since Linnaeus's time, and its species rich vegetation includes the largest population of Fritillaria meleagris in Sweden. To test the effect of different OTU generation methods, we sampled soils across an abrupt moisture transition that divides the meadow community into a Carex acuta dominated plant community with low species richness in the wetter part, which is visually distinct from the mesic-dry part that has a species rich grass-dominated plant community including a high frequency of F. meleagris. We used the moisture and plant community transition as a framework to investigate how detected belowground microeukaryotic community composition was influenced by OTU generation methods. Soil communities in both moisture regimes were dominated by protists, a large fraction of which were taxonomically assigned to Ciliophora (Alveolata) while 30%-40% of all reads were assigned to kingdom Fungi. Ecological patterns were consistently recovered irrespective of OTU generation method used. However, different methods strongly affect richness estimates and the taxonomic and phylogenetic resolution of the characterized community with implications for how well members of the microeukaryotic communities can be recognized in the data.

In situ aerobic composting eliminates the toxicity of Ageratina adenophora to maize and converts it into a plant- and soil-friendly organic fertilizer

Ageratina adenophora has invaded many subtropical and tropical countries and caused tremendous ecological and economic losses. This necessitates a new way to use the debris left after clearing this plant. Therefore, the allelochemicals in fresh and aerobically composted A. adenophora plants (FA and CA, respectively) were compared, and their allelopathy against maize was evaluated. The results showed that CA decreased the allelochemicals (6-hydroxy-5-isopropyl-3,8-dimethyl-4a,5,6,7,8,8a-hexahydronaphthalen-2(1H)-one and 4,7-dimethyl-1-(propan-2-ylidene)-1,4,4a,8a-tetrahydronaphthalene-2,6(1H, 7H)-dione) by over 95% compared to FA. In a seed germination test, CA aqueous solutions improved the seed germination and seedling growth, whereas FA solutions led to opposite results. Chemical fertilizers (CF) plus FA resulted in much lower plant biomass and nutrient uptake than CF in a greenhouse experiment. Compared with CF, CF+CA showed positive effects on maize, soil microbial biomass and diversity and enzyme activities in the field. However, the compositions of the predominant microbes were almost unaffected by the application of CA and CF+CA. These significant findings extended our knowledge regarding the elimination of A. adenophora toxicity against other plants and soil microbes through allelochemical degradation in the composting process. In situ aerobic composting provides a new, simple and economical method to convert A. adenophora into a plant- and soil-friendly organic fertilizer.

Imprinting statistically sound conclusions for gut microbiota in comparative animal studies: A case study with diet and teleost fishes

Despite the technical progress in high-throughput sequencing technologies, defining the sample size which is capable of yielding representative inferences in metabarcoding analysis still remains debatable. The present study addresses the influence of individual variability in assessing dietary effects on fish gut microbiota parameters and estimates the biological sample size that is sufficient to imprint a statistically secure outcome. European sea bass (Dicentrarchus labrax) and gilthead sea bream (Sparus aurata) were fed three alternative animal protein diets and a fishmeal control diet. Gut microbiota data from 12 individuals per diet, derived from Illumina sequencing of the V3-V4 region of the 16S rRNA gene, were randomized in all possible combinations of n-1 individuals. Results in this study showcased that increasing the sample size can limit the prevalence of individuals with high microbial load on the outcome and can ensure the statistical confidence required for an accurate validation of dietary-induced microbe shifts. Inter-individual variability was evident in the four dietary treatments where consequently misleading inferences arose from insufficient biological replication. These findings have critical implications for the design of future metabarcoding studies and highlight the urgency in selecting an adequate sample size able to safely elucidate the dietary effects on fish gut microbial communities.

Soil Bacterial Communities Exhibit Strong Biogeographic Patterns at Fine Taxonomic Resolution

Bacteria have been inferred to exhibit relatively weak biogeographic patterns. To what extent such findings reflect true biological phenomena or methodological artifacts remains unclear. Here, we addressed this question by analyzing the turnover of soil bacterial communities from three data sets. We applied three methodological innovations: (i) design of a hierarchical sampling scheme to disentangle environmental from spatial factors driving turnover; (ii) resolution of 16S rRNA gene amplicon sequence variants to enable higher-resolution community profiling; and (iii) application of the new metric zeta diversity to analyze multisite turnover and drivers. At fine taxonomic resolution, rapid compositional turnover was observed across multiple spatial scales. Turnover was overwhelmingly driven by deterministic processes and influenced by the rare biosphere. The communities also exhibited strong distance decay patterns and taxon-area relationships, with z values within the interquartile range reported for macroorganisms. These biogeographical patterns were weakened upon applying two standard approaches to process community sequencing data: clustering sequences at 97% identity threshold and/or filtering the rare biosphere (sequences lower than 0.05% relative abundance). Comparable findings were made across local, regional, and global data sets and when using shotgun metagenomic markers. Altogether, these findings suggest that bacteria exhibit strong biogeographic patterns, but these signals can be obscured by methodological limitations. We advocate various innovations, including using zeta diversity, to advance the study of microbial biogeography.IMPORTANCE It is commonly thought that bacterial distributions show lower spatial variation than for multicellular organisms. In this article, we present evidence that these inferences are artifacts caused by methodological limitations. Through leveraging innovations in sampling design, sequence processing, and diversity analysis, we provide multifaceted evidence that bacterial communities in fact exhibit strong distribution patterns. This is driven by selection due to factors such as local soil characteristics. Altogether, these findings suggest that the processes underpinning diversity patterns are more unified across all domains of life than previously thought, which has broad implications for the understanding and management of soil biodiversity.

Inoculum handling alters the strength and direction of plant-microbe interactions

The pooling of soil samples in plant-microbe interaction studies is commonly employed, but the impact of sample handling has rarely been explored experimentally. Concerns have been raised that sample pooling may reduce biological variation leading to inflated type I errors or may alter the magnitude of microbial effects observed, invalidating the results achieved. To assess the impact of inocula pooling on plant-microbe interactions, we examined the reciprocal influence of unpooled and pooled soil microbial inocula on growth of Solidago altissima and Schizachyrium scoparium, with and without inoculum sterilization. Soil pooling had no effect on the variance among replicates in either plant species. However, pooling dramatically altered the magnitude and direction of microbial impacts on plant performance. Pooling of Solidago altissima soil increased the antagonistic effects on growth of both target species. In contrast, pooling of Schizachyrium scoparium soil shifted impacts on Solidago altissima from effectively neutral to slightly positive. Pooling in this system altered both the strength and direction of plant-microbe interactions relative to unpooled soils. Therefore soil mixing should be avoided when the research goal is to determine naturally occurring interaction strengths, even within a single habitat.

Faecal pollution affects abundance and diversity of aquatic microbial community in anthropo-zoogenically influenced lotic ecosystems

The aquatic microbiota is known to be an important factor in the sustainability of the natural water ecosystems. However, the microbial community also might include pathogens, which result in very serious waterborne diseases in humans and animals. Faecal pollution is the major cause of these diseases. Therefore, it is of immense importance to assess the potential impact of faecal pollution, originating from both anthropogenic and zoogenic sources, on the profile of microbial communities in natural water environments. To this end, the microbial taxonomic diversity of lotic ecosystems in different regions of Norway, representing urban and rural areas, exposed to various levels of faecal pollution, was investigated over the course of a 1-year period. The highest microbial diversity was found in rural water that was the least faecally polluted, while the lowest was found in urban water with the highest faecal contamination. The overall diversity of the aquatic microbial community was significantly reduced in severely polluted water. In addition, the community compositions diverged between waters where the dominant pollution sources were of anthropogenic or zoogenic origin. The results provide new insight into the understanding of how faecal water contamination, specifically that of different origins, influences the microbial diversity of natural waters.