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. Supplementary Information The online version contains supplementary material available at 10.1007/s00248-023-02259-x.


Introduction
Documenting and explaining large-scale biodiversity patterns associated with environmental gradients (e.g.elevation or latitude) are major topics in biogeography and macroecology [37].Theories and methods of macroecology were historically developed using plants and animals, but microbial biogeography is now a dynamic research field [9,26,41], largely due to the development of high-throughput sequencing (HTS) of environmental DNA (eDNA) [35,46].Microbial macroecology originated from the debate over microbial biogeography, the initially dominant idea being that microbes did not have restricted geographical distribution [21,22].This view stemmed from the early microscopic studies that did not allow observation of fine morphological details [17].Methodological improvements and especially the development of molecular tools revealed the true diversity of microorganisms and more complex distribution patterns [56,57], leading to the moderate endemicity theory, with dispersal limitation being one of the key mechanisms explaining biogeographical patterns [4,35,36,68].Nevertheless, the distribution patterns and drivers of microbial diversity remain less well documented than for macroscopic plants and animals [54].One of the remaining challenges is that soils are characterized by a high degree of spatial heterogeneity, which represents a potential source of bias when exploring biodiversity patterns over large spatial scales (e.g., elevation or latitudinal gradients), especially when only few samples are available across large surfaces [60].
As for biogeography in general, studies of microbial diversity patterns along elevational gradients have become increasingly numerous in the last decade [63].This interest is motivated by their vast diversity, important ecological functions, relevance for climate change impact assessment (bioindicators and roles in biogeochemical cycles), and the complexity of the underlying driving mechanisms across global mountain systems [38,69].However, in order to compare elevational diversity patterns among different microbial groups and mountains, methodological aspects including the sampling strategy need to be optimized and standardized.Many studies followed the same study design as plant elevational diversity studies, which is five to ten samples on each elevational band along an elevation gradient (hereafter referred to as the replicate strategy) [9,39,55].As soil microbial community composition is known to being highly heterogeneous [1,44], it may be necessary to estimate this variability by replicating sampling plots within each elevation band.This may allow taking advantage of pre-existing information and making direct comparisons between soil microorganisms and above-ground macroorganisms possible.Alternatively, a single plot sample may be analyzed at each elevation band (hereafter referred to as the transect strategy) [19,62].A rationale for this choice is that the true replication unit is the elevation gradient and sampling multiple plots per elevation equates to pseudo-replication.Indeed, samples collected within elevational band are likely similar due to the fact that they are taken from the same or very similar environment in a limited area [60].This can lead to inflated sample size, false-positive results, and incorrect statistic inference [27].The reluctance in adopting the transect strategy could be due to its inability to capture the full range of community variation, resulting in insufficient sampling [43].
Microbial eDNA studies are often criticized for low or lack of replication [24,48,64].This limitation was at first mostly due to the financial cost of molecular analyses [48].Another important aspect to consider when discussing replication is the cost of field work, apart from the administrative difficulty in accessing the study sites.Researchers often obtain only a limited number of samples at high financial and time cost when working in remote areas.The transect strategy, requiring fewer resources, facilitates the study of biogeography and macroecology which favors sampling across diverse habitats and broader geographical regions [2,7].Nevertheless, replication is essential in ecological research aiming to reliably estimate the diversity variation of communities, thus allowing robust comparisons between multiple locations [20].Low replication can lead the null hypothesis to being erroneously accepted [32].Therefore, the degree of replication is a clear trade-off with sampling efforts and sound sampling design in ecological research.With chronically limited resources, it is crucial to optimize sampling strategies for estimating diversity patterns.Indeed, in the practical world, replication may be difficult to achieve due to logistic reasons and possible legal, cultural, and habitat conservation limitations, such as the amount of soil that can be taken from nature reserve.
Non-replicated studies do not necessarily prevent statistical robustness.Microbial biogeography usually involves sampling along environmental gradients, and the diversity patterns are estimated with regression.As regression does not require replication [32], there may not be such a thing as a "replication crisis" [20].Using simulated artificial data, Schweiger and colleagues (2016) demonstrated that having no replicate from many sites performed better than having some replicates from fewer sites, as true patterns might be obscured by lacking continuity along the focal gradient.Increasing the sampling intensity by reducing the elevational interval between adjacent samples also helps in detecting a strong elevational microbial diversity pattern in the Andes, highlighting the importance of continuous sampling in microbial biogeography [43].However, the effectiveness of the transect strategy compared to the replicate strategy has not been evaluated in empirical study.Such topic merits further research to explore how spatially sampling efforts should be allocated for investigating microbial diversity patterns efficiently.
This study focuses on a sampling trade-off that can be summarized by the general question: to what degree is replication within an elevation belt necessary to optimally assess alpha and beta diversity patterns along elevation gradients?To address this question, we collected soil samples and extracted eukaryotic microbial eDNA in natural vegetation along an elevational gradient using two different designs on the slope of Mt.Asahi, Hokkaido: five replicates at each of six elevational bands and one sample at each of 16 elevational bands.We hypothesized that both taxonomic richness and community composition would show moderate degrees of heterogeneity (i.e., similar alpha diversity and low beta diversity) among plots within a given elevation band but higher heterogeneity among elevation bands.If true, this would imply that (1) a single plot per elevation would be sufficient to obtain a good estimate of diversity patterns along an elevation gradient and (2) for a given total number of samples, gamma diversity estimates for the entire transect would be highest and beta diversity patterns along the elevation transect better assessed if a single sample was analyzed per elevation band and a higher number of elevation bands analyzed.This study would allow us to determine if little to no replicate could detect a pattern and how it may affect accessing the underlying mechanisms shaping biodiversity patterns.The goal of this study is therefore to optimize the sampling and analysis effort for the assessment of soil eukaryotic microbe diversity patterns along elevation gradients.

Field Work
Mt. Asahi is a volcanic mountain of the Taisetsu Mountain Range located in central Hokkaido Island, Japan (43.66°N, 142.85°E).It is the tallest peak in Hokkaido which tops at 2290 m above sea level.Mt.Asahi lies in the hemiboreal zones characterized by severe winter conditions with snow cover potentially reaching 400 cm in extreme circumstances [45].The vegetation sequence from low elevation to the summit is deciduous broadleaf, mixed forest, conifer forest, thickets of dwarf pine at treeline around 1300 ~ 1400 m followed by alpine tundra, and patchy vegetation of herbaceous plants and dwarf shrubs at the top [45].In this study, we defined the area below the treeline, which is situated around 1400 m in our elevational transect, as forest habitat and the area above the treeline as alpine habitat.
To test how replicates may influence the estimates of microbial diversity patterns along the elevational gradient on Mt.Asahi, we collected soil samples using two sampling strategies: (1) The replicate strategy involved collecting five replicate samples at six elevational bands with ca.300-m elevation intervals, including three bands in forest (500, 900, and 1300 m.a.s.l.) and three in alpine habitat (1620, 1940, and 2260 m.a.s.l.).Within each of these elevational bands sampled for the replicate strategy, five samples were collected in plots that were distant by at least 50 m, resulting in a total of 30 samples.(2) The transect strategy involved collecting a single sample in each of 16 elevational bands (ten in forest and six in alpine habitats) with ca.100-m elevation intervals.Samples from the replicate strategy were utilized in both strategies, and, therefore, ten additional elevational bands were sampled for the transect strategy at 400, 600, 700, 800, 1000, 1100, 1200, 1460, 1780, and 2100 m.a.s.l.(Fig. 1).The total number of samples was therefore 40.The size of the sampling plots was 10 × 10 m in the forest and 5 × 10 m in the alpine.For each sample, ca. 100 g of leaf litter and upper 5 cm of soil were collected with a knife for 10 sub-samples, representative of the micro-habitats within the plot, and subsequently pooled into a 1-kg composite sample.Such composite samples are sufficient to capture the maximum diversity for a relatively small, locally heterogeneous area [58].

Laboratory Work and eDNA Sequencing
The composite soil samples were homogenized and sieved through 5-and 2-mm mesh.Ca. 0.5 g of this sieved material was fixed in LifeGuard DNA preservation buffer for DNA analyses.The remaining part of each sample was used to analyze basic edaphic physicochemical characteristics (Table 1).Five grams of soil was mixed with distilled water in a 1:2.5 (wt/vol) prepared for pH measurement with Metrohm 621 pH meter.Water content (WC) was measured Fig. 1 Sampling design on Mt.Asahi, Hokkaido, Japan.Blue dots represent locations where only one sample was collected, and red dots represent locations where 5 samples were taken from the same elevational bands.Digital elevation model data was downloaded from Geospatial Information Authority of Japan (www.gsi.go.jp).Elevation change was exaggerated three times relative to geographical distance by heating the soil at 105 °C for 24 h.The determination of soil organic matter (SOM) was conducted by performing loss-on-ignition at 450 °C using a muffle furnace (Nabertherm).Organic carbon (C) and nitrogen (N) were measured with the FLASH 2000 CHN analyzer (Thermo Fisher Scientific).Bioavailable phosphorous was measured by colorimetry [47].
Environmental DNA was extracted from 0.25 g of soil using the PowerSoil Pro DNA extraction kit (Qiagen, Venlo, Netherlands).The V4 region of the SSU rDNA gene was amplified by PCR with universal eukaryotic primers TAReuk454FWD1/ TAReukREV3 targeting the V4 hypervariable region [59].Pooled triplicate PCR products were sequenced with Illumina MiSeq (2*300 bp paired-end reads).Bioinformatic processing was conducted on a high-performance computing system using DeltaMP v0.5 [33].Raw reads were first demultiplexed with maximum of two mismatches on the barcodes and three mismatches on the primer sequences using Cutadapt v2.10 [40].End of reads was truncated to a minimum sequence length of 230 nt so that the expected error rate was below 4 using VSEARCH v2.13.6 [51].Reads not passing this quality filters were eliminated.Amplicon sequence variants (ASV) were called following the classical workflow of DADA2 [10] and pair-end joined with a minimum overlap of 10 nt and maximum of two mismatches in the overlapping region.Chimera were removed using the DADA2 function removeBimer-aDenovo.Taxonomic assignment of ASV was done against the PR2 database [23] using global pairwise alignment from VSEARCH.Tag jumps were controlled and removed using ASVs from positive controls from two cultivated aquatic algal species.Detailed information on soil chemical analyses, molecular analyses, and bioinformatic processing was given in Lentendu et al. [34].Raw sequence data are available on the SRA (bioproject PRJEB60066).

Bootstrapping
The two sampling strategies utilized the same elevational transect, allowing for the selection of any sample within an elevational band with replicates for the transect strategy.To obtain a more accurate estimation of diversity patterns and facilitate comparisons with the replicate strategy, we constructed 1000 potential sample combinations for the transect strategy.These combinations were generated with bootstrapping by randomly selecting one out of the five replicates from each of the six bands with replicates [16], along with samples from elevational bands without replicate (Fig S1).The diversity patterns estimated from 1000 potential combinations of the transect strategy were averaged.

Alpha Diversity
For alpha diversity, we employed the diversity estimation framework of the Hill number developed by Hill [25] and advanced by Chao et al. [13]: , q ≥ 0, q ≠ 1 where p i is the relative abundance or relative incidence frequency for i th species in an abundance sampling unit or sampling unit-based replicated occurrence data with S species.The parameter q determines the weight given to rare species.When q = 0, all species are weighted equally, so that 0 D is the species richness.Where q = 1, all individuals are given equal weight, and 1 D is the limit of q D as the order q tends to 1 which reduces to the exponential Shannon entropy.Where q = 2, the index inclines to the dominant species, and 2 D is the inverse of the Simpson index [13].q D is regarded as the effective number of equally abundant species [15].In this study, ASVs were treated as species and the number of sequences of each ASV was treated as abundance information.To remove the effects of sequencing depth on the estimation of alpha diversity, we calculated the Hill numbers when the sample size was rarefied/extrapolated to double the size of the minimum number of sequences which was 3834 (this was done after removing all the other sequences except eukaryotic microbes).Diversity estimation is reliable up to the twofold extrapolation of the sequencing depth [13].In this way, we maximized the use of the ASVs by discarding fewer sequences.We applied linear regression to estimate the diversity pattern between the Hill number (when q = 0, 1, and 2) and elevation with the replicate strategy and all the datasets of the transect strategy.Then, we tested if the regression coefficients (slopes) from the replicate strategy were different from each of the transect strategy using analysis of covariance (ANCOVA).Regression coefficients were averaged, and the p values were corrected for multiple tests using the procedure of Benjamini and Hochberg [3].

Gamma Diversity
To determine if the replicate strategy caught more overall eukaryotic microbe diversity for the whole transect (gamma diversity) than the transect strategy, we rarefied the ASV richness (q = 0) based on number of samples using the ASV occurrence data.All 1000 datasets of the transect strategy were rarefied to obtain an averaged rarefaction curve.ASV richness were extrapolated to twice the number of sampling units.

Beta Diversity
The patterns of soil eukaryotic microbe communities along the elevation gradient were assessed by principal coordinate analysis (PCoA) on the Bray-Curtis dissimilarity [65] calculated on the Hellinger-transformed percentage data.The explanatory power of elevation and habitat type was tested using permutational analysis of variance (PERMANOVA) on the replicate and transect strategies separately.p values were corrected and averaged together with F values for the PERMANOVA tests on 1000 datasets of the transect strategy.To assess beta diversity variation within and among elevational bands, we estimated the elevation-decay relationship with a linear model and tested, with ANOVA and Tukey's test, for each elevational interval if the beta diversity among different elevation bands was greater than the beta diversity within elevational band.We analyzed forest and alpine habitats separately due to their slightly different elevational sampling intervals.We also assessed the predictive value of environmental variables on the eukaryotic microbial community structure (beta diversity) using redundancy analysis (RDA) on both strategies.For addressing the collinearity of environmental variables, we first checked their variance inflation factors (VIF) for samples altogether.We found that SOM, TOC, TON, and WC had VIFs greater than ten which means they were highly correlated.Bivariate correlations between all pairs of these four variables also showed high collinearity (all Pearson's r > 0.94).Then, principal component analysis (PCA) was conducted on those four variables.The first axis of the principal component (PC1) which accounted for > 97% of variance was selected to represent the four organic matter related variables for further analyses as OM_pc1.Both forward and backward selections were used to determine the best explaining variables in the RDA model.Hierarchical partitioning was used to measure the individual and unique contributions of the selected variables explaining community structure in the parsimonious redundancy model [30].Overall constrained variance (adjusted R 2 ) of the RDA was assigned to each of the predicting variables to assess their relative importance.Similar to alpha diversity, we repeated the RDA and model selection process for all datasets of the transect strategy.Eventually, we checked how the best explaining variables changed between replicate and transect strategies.All statistics and figure production were conducted in R 4.21 [49]: dplyr 1.0.9 for data handling [67],iNEXT.3D1.0.1 for rarefaction and diversity estimation [14],vegan 2.6.2 for PCoA, PCA, PERMANOVA, and RDA [61],rdacca.hp1.0.8 for hierarchical partitioning [30],ggplot2 3.3.6for statistical graphics [66],and rayshader 0.34.6 for sampling location map [42].The full R code to replicate the analyses is available as supplement.

Alpha Diversity
The averaged elevational diversity pattern over 1000 datasets of the transect strategy was congruent with the pattern from the replicate strategy (Fig. 2).Both strategies showed a significant negative relationship between alpha diversity (Hill number of ASVs) and elevation when q = 0, but non-significant relationship when q = 2.For q = 1, a significant relationship was detected with the transect strategy but not with the replicate strategy.The slopes of the transect strategy were steeper than the one of the replicate strategy, but this difference was not significant (ANCOVA, Table 2).

Gamma Diversity
Rarefaction curves of ASV occurrence data showed that the replicate strategy yielded higher observed overall richness, but neither strategy reached saturation even when extrapolated to twice the sample size (Fig. 3).However, for the same given number of samples, ASV richness estimated from the transect strategy with samples from more elevational bands was higher than from the replicate strategy with samples from fewer elevational bands.

Beta Diversity
The PCoA plot demonstrated that samples from a given elevational band tend to cluster together (Fig. 4).PER-MANOVA tests showed that soil micro-eukaryotic community composition was significantly explained by elevation (R 2 = 0.208, p = 0.001), habitat (R 2 = 0.076, p = 0.002), and their interaction (R 2 = 0.116, p = 0.001) using the replicate strategy.PERMANOVA tests on the transect strategy delivered the similar results, but the effect of habitat was not significant (elevation: R 2 = 0.212, p correct = 0.002; habitat: R 2 = 0.075, p correct = 0.1; interaction: R 2 = 0.131, p correct = 0.006).Elevation-decay relationships were significant both in the forest (R 2 = 0.47, p < 0.001) and alpine habitat (R 2 = 0.54, p < 0.001) (Fig. 5).Tukey's test showed beta diversity between samples from the same elevational bands was significantly lower than for samples with elevational interval of 200 m or more, but there was no statistical difference for samples with elevational interval of 100 m in forest habitat.In alpine habitat, beta diversity between samples with 160-m elevation interval or more was significantly higher than samples from the same elevational bands.
RDA results showed that the best explaining variables were OM_pc1, elevation, pH, and C/N ratio which altogether explained 25% of the variance (Fig. 6, Table 3) with OM_pc1 and elevation already accounted for 20% (Fig. 7).Among the 1000 datasets of the transect strategy, the combination that includes OM_pc1 and elevation emerged as the best explaining variables up to 840 times.They were usually the strongest predicting variables for beta diversity variations with high R 2 value (Fig. 7), especially OM_pc1 which was always selected and had the highest explaining power.

Discussion
In this study, based on commonly used diversity indices (species richness, the exponential Shannon entropy, and the inverse of the Simpson index), we found no difference in elevational alpha diversity patterns estimated from the transect strategy and the replicate strategy.Community composition of soil eukaryotic microbes showed moderate heterogeneity, i.e. low beta diversity between plots within same elevational bands or plots with short elevation difference (less than 100-m elevation).The estimated total species richness was higher for the transect strategy than for the replicate strategy for an equal number of samples.In the RDA, the same set of environmental variables explained the community structure of soil eukaryotic microbes for the transect and replicate strategies.Although we demonstrated the reliability of the transect strategy in estimating microbial elevational diversity patterns, the current study does not intend to question the importance and effectiveness of replication in ecological studies.Soil microbial communities are highly heterogeneous even in homogeneous environment at the centimeter scale [1,44].Increased sample size could help us detect rare soil microorganisms and assess the role of the rare biosphere in ecosystem functions [12,24].Sample replicates at fine scale also facilitate the estimation of soil microbial community variation and the underlying mechanisms [50].However, the critical lack of large-scale biogeographic knowledge, in particular, of soil microorganisms motivates the employment of the transect strategy with manageable efforts and limited resources [11].
The trade-off between single vs. replication is always a key issue when designing empirical data collection.Simulation experiments show that the balance between random and systematic errors determine if replication is the better strategy or not [52]: Replication is preferred only in situations where systematic errors (such as spatial autocorrelation and deterministic factors) are well controlled and the data variance is shaped mostly by random error.Most ecological modeling, however, lacks such fine controls of various factors.Our finding was in line with this, i.e., soil eukaryotic microbial communities supported the transect strategy over the replicate strategy.Continuous sampling by reducing replication but increasing factor levels could mitigate the effects of high systemic error and improve the accuracy of diversity estimation [29,52].Of course, regression analysis only requires a few samples, but the true pattern may not be revealed by few factor levels with rather large interval [28].Fierer et al. [19] found no bacterial alpha diversity gradient along a 3450-m elevation   3 transect in the Andes with one sample collected on each of the only six elevational bands.However, by increasing the number of elevation band to 14, they found that bacteria exhibited a decreasing diversity with increasing elevation, similar to plants [43].Although they had three technical replicates for each sample, the finding of significant microbial diversity pattern was not due to the replicates as they had detected weak or no pattern when using the only six elevational bands in their previous study with new replicate dataset [19].The biological replicates in this study caught more variability of microbial diversity, but it is the addition of more elevational bands which allowed for detecting stronger, and likely more realistic, biodiversity patterns for both cases.Our findings confirmed the above-mentioned simulation and empirical studies which highlighted the importance of continuous sampling in estimation of alpha diversity pattern along environmental gradients.More importantly, this study demonstrated that the transect strategy without replication was sufficient to estimate alpha diversity pattern along elevational gradient.Low beta diversity within elevational band has been observed numerous times in soil microbial elevation studies [9,39,55].Microbial communities from the same elevational band are filtered by similar climatic and edaphic conditions and thus exhibit low composition difference [53].On Mt.Asahi, soil eukaryotic microbe communities were significantly different for elevational intervals larger than 100 m, which means that adding more elevational bands facilitated compositional variation estimations.This is also shown by the fact that, for the same total number of samples, single samples from more elevational bands yield higher gamma diversity than replicate samples from fewer elevational bands.Low beta diversity within elevational band might also be caused by the composite sampling method, as it homogenized microbial communities among different micro-habitats [12,18].However, our goal was to catch the maximum diversity in the plot and assess elevational diversity patterns in general rather than to quantify the influence of micro-habitats.The transect strategy only required half the efforts of the replicate strategy, and we were still able to identify the major environmental drivers (soil organic matter related variable and elevation) in shaping soil eukaryotic microbe community.Although some less significant variables may be missed, for example pH and C/N ratio were not frequently identified by the transect strategy but were significant variables in the replicate strategy.We should be cautious with low sample size in multivariate analysis as it inevitably reduces the sensitivity in finding the best predictors [6].We also need to limit the number of variables that we are interested in when we adopt the transect strategy as it cannot surpass the number of samples [6].
The 40 samples in this study took five researchers 5 days to collect in the field and process (mixing and sieving) for fixing DNA at the field station.Increasing the number of elevational bands with fewer replicates on the same elevation is a practical strategy to reduce efforts and costs in the field and laboratory.While still yielding reliable estimation of biodiversity patterns, this approach allows relocating resources to include more transects and mountain regions to generalize the study at a broader scale [34], or to assess the comparative magnitude of elevation vs. seasonal patterns in mountain soil microbial communities [8,31], in the longer term, how they may respond to climate change [5].

Fig. 2 Table 2 Fig. 3
Fig. 2 Elevational alpha diversity patterns of soil micro-eukaryotic diversity along an elevation gradient on Mt.Asahi, Hokkaido, Japan, with different sampling strategies.Hill numbers 0, 1, and 2 represent the diversity metrics of ASV richness, exponential Shannon, and

Fig. 4
Fig.4 Principal coordinate analysis (PCoA) plot of soil micro-eukaryotic community structure along an elevation gradient on Mt.Asahi, Hokkaido, Japan, based on the Bray-Curtis dissimilarity.Ellipses show the 95% confidence interval around replicate samples from each elevational band

Fig. 6 Fig. 5
Fig. 6 RDA of soil microeukaryotic communities (ASVs) along an elevation gradient on Mt.Asahi, Hokkaido, Japan, showing replicate samples with environmental variables.Both forward selection and backward elimination were conducted in model selection of RDA to find the parsimonious variables in explaining the micro-eukaryotic community composition.Ellipses show the 95% confidence interval around replicate samples from each elevational band

Fig. 7
Fig.7 Relative contribution of environmental variables to soil microeukaryotic beta diversity patterns along an elevation gradient on Mt.Asahi, Hokkaido, Japan.Light blue bars represent unique contribution of the selected variable, and gray bars represent the individual contributions which include the unique and shared effect of the particular variable.In some case, the unique effect can exceed the indi-

Table 1
Soil physicochemical variable data of each elevational band.Mean value and standardized deviation are shown for the elevational band with replicate samples Abbreviations: WC, water content; SOM, soil organic matter; TON, total organic carbon; TOC, total organic nitrogen; P, bioavailable phosphorus; C/N, carbon to nitrogen ratio; C/P, nitrogen to phosphorus ratio

Table 3
Parsimonious variables and their explaining power (adjusted R 2 ) of RDA on soil eukaryotic microbe community composition along an elevation gradient on Mt.Asahi, Hokkaido, Japan, using different sampling strategies For the transect strategy, the results of 1000 datasets are shown in descending order based of the frequency of the combination of parsimonious variables