Data and code from: Environmental stress strengthens plant-mycorrhizal associations: A novel extension of the stress-gradient hypothesis
Data files
May 20, 2026 version files 22.22 MB
-
HJA_StressGradientElevationMicrobiome.zip
22.20 MB
-
README.md
15.45 KB
May 27, 2026 version files 10.38 MB
-
HJA_StressGradientElevationMicrobiome.zip
10.38 MB
-
README.md
4.13 KB
Abstract
This dataset supports Doolittle et al. (2026), an investigation of how abiotic environmental stress shapes plant–fungal community relationships across an elevational gradient in the central Cascade Range of Oregon, USA. Sampling was conducted in 2021 at seven forest reference stands (RS02, RS21, RS22, RS27, RS28, RS31, RS38) within the H.J. Andrews Experimental Forest, spanning approximately 500 to 1,290 m elevation a.s.l. Three integrated data types were collected at the quadrat level (25 x 25-m grid) within each stand. Forest inventory data include tree species identity, diameter at breast height (DBH), stem position, and quadrat coordinates from large, high-resolution censuses. Soil chemistry was characterized for each quadrat. Fungal community data come from ITS amplicon sequencing of paired forest-floor (litter) and mineral-soil samples, processed into amplicon sequence variants (ASVs) and annotated with ecological guild assignments (ectomycorrhizal, pathogenic, saprotrophic, and unassigned) grouped from FUNGuild output. The repository contains the sample metadata, processed soil chemistry tables, FUNGuild ASV taxonomic assignments, ASV abundance tables, and tree census files needed to reconstruct the phyloseq objects used in analysis. R scripts reproducing the full analytical pipeline, data preparation, alpha- and beta-diversity calculations, PERMANOVA, PCoA on both raw and ZIPPCA-nb-denoised ASV counts, and variation partitioning of plant–soil–fungal community relationships are included alongside the data so the published results can be re-executed end-to-end from the raw inputs.
This repository contains the data and R scripts that reproduce all analyses and figures in:
Doolittle et al. (2026). Environmental stress strengthens plant–mycorrhizal associations: A novel extension of the stress-gradient hypothesis.
A description of the dataset is provided in dryad_abstract.md. Repository is shipped as a zipped file, to view layout please open HJA_StressGradientElevationMicrobiome.zip.
Repository layout
.
├── HJA_StressGradientElevationMicrobiome.Rproj
├── README.md
├── dryad_abstract.md # Dataset abstract (Dryad)
├── data/ # Raw input data
│ ├── FGEO_SoilSamples_Meta2021.csv # sample metadata
│ ├── oregon_funguild_sttax_table.guilds.rds # FUNGuild output
│ ├── soildata.csv # per-quadrat soil chemistry
│ ├── soildata_quadlocs.csv # quadrat coordinates (gx, gy)
│ └── tree_data/
│ ├── (HJA_refstand_data.RData) # NOT distributed here, and NOT used in main analyses, only suppelementary — see "Required external data" below
│ ├── RS02_main_plot.RData # tree census, large main-plot at RS02
│ ├── RS21_main_plot.RData # tree census, RS21
│ ├── RS28_main_plot.RData # tree census, RS28
│ ├── RS31_main_plot.RData # tree census, RS31
│ └── RS38_main_plot.RData # tree census, RS38
├── scripts/ # Numbered analytical pipeline
│ ├── 00_set_workspace_prep_phyloseq.R
│ ├── 01_subset_ps_div_metrics.R
│ ├── 02_hja_soil_data.R
│ ├── 03_hja_tree_prep.R
│ ├── 04_beta_analyses.R
│ ├── 05_beta_deviation_models.R
│ └── 06_dispersion_tests.R
├── output/ # Created by scripts (RDS, tables)
└── figures/ # Created by scripts (TIFF, PNG, PDF)
Required external data
scripts/03_hja_tree_prep.R reads one file that is not redistributed in this repository. This dataset is not used in the main analyses, only supplementary tests.
-
data/tree_data/HJA_refstand_data.RData— the H.J. Andrews reference-stand survival and growth dataset originally published by:LaManna, J. A., F. A. Jones, D. M. Bell, R. J. Pabst, and D. C. Shaw. 2022. "Tree species diversity increases with conspecific negative density dependence across an elevation gradient." Ecology Letters 25(5): 1237–1249. https://doi.org/10.1111/ele.13996
Download the original
.RDatafile from LaManna et al.'s deposit and place it atdata/tree_data/HJA_refstand_data.RData. The file should contain the objecthja.survival.growth.data, a nested list with one entry per reference stand. Script 03 indexes into this list withhja.survival.growth.data$RS28[7]and…$RS38[5].
This dataset is not used in the main analyses, only supplementary tests. The remaining analyses (scripts 01, 02, 05, 06) do not depend on the reference-stand data and will run from the files included here.
Notes: Site codes RS02–RS38 are reference stands within the H.J. Andrews Experimental Forest, central Cascades, Oregon. Approximate elevations: RS38 = 500 m, RS02 = 520 m, RS27 = 790 m, RS31 = 900 m, RS28 = 1,060 m, RS21 = 1,190 m, RS22 = 1,290 m. Site RS22 is dropped from most analyses due to limited sample recovery.
Raw FASTQ files are deposited at NCBI under BioProject PRJNA1210581 (https://www.ncbi.nlm.nih.gov/bioproject/PRJNA1210581). The processed ASV table and FUNGuild guild annotations distributed here as data/oregon_funguild_sttax_table.guilds.rds were derived from those raw reads.
Changes after May 20, 2026: Improved readme, fixed figure generation, harmonized repository structure.
