Study number PRJEB10570 at the European Nucleotide Archive contains the raw MiSeq reads from high-throughput sequencing of the 16S rRNA gene in hundreds of leaf, root, and soil samples. Note that some of these samples are associated with other experiments and were not analyzed for the Wagner et al. 2016 study.

Estimates of 16S rRNA gene copy number for OTUs analyzed in the Wagner et al. 2016 study. Created by following Kembel et al. (2012) Incorporating 16S Gene Copy Number Information Improves Estimates of Microbial Diversity and Abundance. PLoS Computational Biology 8 (10), e1002743.

Collection of glucosinolate data is described in Supplementary Note 1 of the article. HPLC peaks were called and area under each curve was integrated automatically by Agilent ChemStation software. Columns beginning with "area_" include the area under chromatogram peaks corresponding to the internal standard (Sinigrin) or one of the glucosinolate compounds (2OH1ME, 1ME, 6MSOH, 1MP, or I3M). "Weight" is the mass of air-dried samples, in mg. "Batch" is the identifier for groups of samples measured in the same HPLC run. "HPLC_ID" is a unique identifier for a glucosinolate sample, and "Plant_ID" is a unique identifier for individual plants in the field experiment.

Phylogeny of all OTUs (based on 97% 16S rRNA gene sequence identity). Generated from representative sequences using the midpoint method in QIIME.

Latitude and longitude of 5 sites where genotypes were collected from natural B. stricta populations / where experimental common gardens were located.

Data on soil & vegetation characteristics from field sites used in this experiment. "PlantDiv" = number of plant morphospecies present in each block; “Veg” = percent of each block covered by vegetation (aggregate of estimated percent cover for 50 10x10cm sub-plots per block). Units are ppm (except for pH, unitless; and conductivity, umho/cm).

Microbiome sample metadata. "SampleID" = unique identifier for a sequenced 16S rDNA amplicon pool corresponding to a single sample. "Name" = more descriptive unique sample name. "Plant_ID" = unique identifier for individual plants; may correspond to >1 SampleIDs when both roots and leaves were sampled from the same plant; "Experiment" = sub-experiment that a sample belongs to (e.g., "ecotypes" = main field experiment described in Wagner et al. 2016; "ecoGH" = greenhouse experiment described in supplementary information of Wagner et al. 2016; "rilGH" and "fieldBCMA" are not presented in Wagner et al. 2016); "Cohort" = year the individual (Plant_ID) was planted, if applicable; "Harvested" = year sample was harvested; "Treatment" = type of soil the individual was planted in (only applicable to greenhouse experiments); "oldPlate" = identifier of the 96-well plate containing the sample before full randomization; "newPlate" = MiSeq run, i.e., identifier of the amplicon library into which the sample was pooled, i.e., the 96-well plate containing the sample after full randomization; "Analysis" = broadest experimental grouping of samples; either "Ecotypes" (which contains samples described in Wagner et al. 2016) or "BCMA" (containing samples not analyzed or described in Wagner et al. 2016).

OTU table (for OTUs picked based on 97% 16S rRNA gene sequence identity) that forms the basis for most analyses in Wagner et al. 2016. Generated from raw MiSeq reads (see "Bioinformatics pipeline"). Columns are sample names (correspond to "SampleID" column in "Sample metadata" file); rows are OTUs. Decompress file using the Linux/UNIX command ' bunzip2 otuTable97.txt.bz2 '

List of OTUs identified as likely contaminants based on comparison to list of known contaminant 16S rRNA gene sequences (generated using custom scripts; see "Bioinformatics pipeline")

Taxonomic assignments of OTUs (binned at 97% 16S rRNA gene sequence identity), generated by comparing representative OTU sequences to the Greengenes database using the RDP classifier in Qiime (see "Bioinformatics pipeline").

Taxonomic assignments of OTUs (binned at 99% 16S rRNA gene sequence identity), generated by comparing representative OTU sequences to the Greengenes database using the RDP classifier in Qiime. * Note that these are OTUs based on 99% identity, NOT the OTUs that were used for most analyses in Wagner et al. 2016 *

Representative sequence for each OTU (based on 97% sequence identity), in FASTA format. Generated from raw MiSeq reads using FLASH and UPARSE software (see "Bioinformatics pipeline")

OTU table (for OTUs picked based on 99% 16S rRNA gene sequence identity) that forms the basis for most analyses in Wagner et al. 2016. Generated from raw MiSeq reads. Columns are sample names (correspond to "SampleID" column in "Sample metadata" file); rows are OTUs. Decompress file using the Linux/UNIX command ' bunzip2 otuTable99.txt.bz2 ' ** note that most analyses in this paper used OTUs picked based on 97% sequence similarity-- these are in a different file (otuTable97.txt) **

Master list of all individuals originally planted in the field experiment- columns as in the Sample Metadata file (SMD.txt)

Archive of R scripts used to produce results in Wagner et al. 2016. See README file for notes on contents and recommendations for use.

A description of the bioinformatics pipeline used to process raw MiSeq data, including parameter values and sample syntax.