DOI: 10.5061/dryad.b1296
title: Data from: Multiple plant traits shape the genetic basis of herbivore community assembly

description: 1. Community genetics research has posited a genetic basis to the assembly
 of ecological communities. For arthropod herbivores in particular, there
 is strong support that genetic variation in host plants is a key factor
 shaping their diversity and composition. However, the specific plant
 phenotypes underlying herbivore responses remain poorly explored for most
 systems. 2. We address this knowledge gap by examining the influence of
 both genetic and phenotypic variation in a dominant host-plant species,
 Salix hookeriana, on its associated arthropod herbivore community in a
 common garden experiment. Specifically, we surveyed herbivore responses
 among five different arthropod feeding guilds to 26 distinct S. hookeriana
 genotypes. Moreover, we quantified the heritability of a suite of plant
 traits that determine leaf quality (e.g. phenolic compounds, trichomes,
 specific leaf area, C : N) and whole-plant architecture, to identify which
 traits best accounted for herbivore community responses to S. hookeriana
 genotype. 3. We found that total herbivore abundance and community
 composition differed considerably among S. hookeriana genotypes, with
 strong and independent responses of several species and feeding guilds
 driving these patterns. We also found that leaf phenolic chemistry
 displayed extensive heritable variation, whereas leaf physiology and plant
 architecture tended to be less heritable. Of these traits, herbivore
 responses were primarily associated with leaf phenolics and plant
 architecture; however, different herbivore species and feeding guilds were
 associated with different sets of traits. Despite our thorough trait
 survey, plant genotype remained a significant predictor of herbivore
 responses in most trait association analyses, suggesting that unmeasured
 host-plant characteristics and/or interspecific interactions were also
 contributing factors. 4. Taken together, our results support that the
 genetic basis of herbivore community assembly occurs through a suite of
 plant traits for different herbivore species and feeding guilds. Still,
 identifying these phenotypic mechanisms requires measuring a broad range
 of plant traits and likely further consideration of how these traits
 affect interspecific interactions.
description: Table S7: Willow trait
 REML modelsThis R script reproduces
 Table S7 of the manuscript entitled, "Multiple plant traits shape the
 genetic basis of herbivore community assembly." This script contains
 all of the restricted maximum likelihood (REML) models and broad-sense
 heritability estimates of the 40 willow (Salix hookeriana) traits measured
 in this common garden experiment.trait_REML_table.RTable S4: Heatmap of willow phenotypic correlationsThis R script reproduces Table S4 of the manuscript entitled, "Multiple plant traits shape the genetic basis of herbivore community assembly." Specifically, this script calculates willow (Salix hookeriana) phenotypic correlations (Pearson's r) among the 40 traits measured in this common garden experiment. This script also generates a heat map for visualizing these correlations.trait.corr_plot.RFig. 3: Multipanel figure of genetic variation in willow traitsThis script reproduces Fig. 3 of the manuscript entitled, "Multiple plant traits shape the genetic basis of herbivore community assembly." Specifically, it calculates and plots the mean and standard error of several plant traits measured among the 26 willow genotypes used in this common garden experiment.plant.trait_plots.RFig. 2: Multipanel figure of herbivore guild responses to willow genotypeThis R script reproduces Fig. 2 in the manuscript entitled, "Multiple plant traits shape the genetic basis of herbivore community assembly." Specifically, this script calculates and plots the mean and standard error of the abundance of each feeding guild occurring on the 26 willow genotypes measured in this common garden experiment.herb.guild_responses_genotype.RFig. 1 & S6: Herbivore community responses to willow genotypeThis R script reproduces Fig. 1 & S6 in the manuscript entitled, "Multiple plant traits shape the genetic basis of herbivore community assembly."herb.com_responses_genotype.RSummary of results for herbivore community responses to willow genotype.This R script runs one-way ANOVA, GLMs, and RDA models of herbivore community responses to willow genotype. These results are reported in the first 2 paragraphs of the Results section of the manuscript entitled, "Multiple plant traits shape the genetic basis of herbivore community assembly."herb.com_anova_table.RTable 1 & 2: Phenotypic and genetic correlations of herbivore guild and species responsesThis R script reproduces the results reported in Table 1 & 2 of the manuscript entitled, "Multiple plant traits shape the genetic basis of herbivore community assembly. Specifically, this script calculates Pearson correlations for phenotypic and genotypic correlations of herbivore community.herb.com_corr_table.RTable 3: willow traits that best predict herbivore community responsesThis R script runs the forward model selection approach advocated by Blanchet et al. (2008, Ecology) to identify which willow traits best predict herbivore community responses. After identifying these traits, I then determined whether willow genotype still had a significant effect. These results are reported in Table 3 of the manuscript entitled, "Multiple plant traits shape the genetic basis of herbivore community assembly."herb.trait_regressions.RTable S5: Principal Component Analysis of willow phenolic compoundsThis R script replicates the principal component analyses we performed on different classes of phenolic compounds (Table S5) in the manuscript entitled, "Multiple plant traits shape the genetic basis of herbivore community assembly." We explored different transformations (sqrt, log, log(x+1)) to linearize the relationships among leaf phenolic compounds. We then used scree plots and visual inspection of phenolic compound loadings to choose interpretable principal componenents. All PCs were performed on the correlation matrix to give equal weight to each phenolic compound.phenolic_pcas.RWillow herbivore community and plant trait management scriptThis R script manages the herbivore community and plant trait data for the manuscript entitled, "Multiple plant traits shape the genetic basis of herbivore community assembly". This file is sourced in for all of the other R scripts associated with this manuscript so that all of the figures, tables, and analyses are reproducible.herb_trait_manage.rWillow genotype, herbivore community, and plant trait dataContains all of the data that is analyzed and reported in the manuscript entitled, "Multiple plant traits shape the genetic basis of herbivore community assembly."herb.trait.data.csv