Heterozygous Trees Rebound the Fastest after Felling by Beavers to Positively Affect Arthropod Community Diversity

: Although genetic diversity within stands of trees is known to have community-level consequences, whether such effects are present at an even ﬁner genetic scale is unknown. We examined the hypothesis that genetic variability (heterozygosity) within an individual plant would affect its dependent community, which adds a new dimension to the importance of genetic diversity. Our study contrasted foliar arthropod community diversity and microsatellite marker-derived measures of genetic diversity of cottonwood ( Populus fremontii ) trees that had been felled by beavers ( Castor canadensis ) and were resprouting, relative to adjacent standing, unfelled trees. Three patterns emerged: 1. Productivity (speciﬁc leaf area), phytochemical defenses (salicortin), and arthropod community richness, abundance, and diversity were positively correlated with the heterozygosity of individual felled trees, but not with that of unfelled trees; 2. These relationships were not explained by population substructure, genetic relatedness of the trees, or hybridization; 3. The underlying mechanism appears to be that beaver herbivory stimulates increased productivity (i.e., 2 × increase from the most homozygous to the most heterozygous tree) that is the greatest in more heterozygous trees. Salicortin defenses in twigs were also expressed at higher concentrations in more heterozygous trees (i.e., 3 × increase from the most homozygous to the most heterozygous tree), which suggests that this compound may dissuade further herbivory by beavers, as has been found for other mammalian herbivores. We suggest that high stress to trees as a consequence of felling reveals a heterozygosity– productivity linkage, which in turn is attractive to arthropods. Although experiments are required to demonstrate causality, these results link the genetic diversity of individual trees to community diversity, supporting the hypothesis that interactions among foundation species (beavers and trees) have community-level effects, and underscores the importance of genetic


Introduction
A growing number of reviews show that dominant or foundation species can have extended genetic effects that impact associated communities, ecosystem processes, and evolution [1][2][3][4][5][6][7][8]. These effects are a consequence of genetic variability, which is organized hierarchically. The broadest levels of genetic effects involve a number of unique individuals (genotypes) within a population, which can consist of genotypes within species or, more loosely, hybrid genotypes of closely related sympatric species. Such genetic variation is most accurately referred to as genotypic diversity [9]. Molecularly assessed genetic diversity, also usually a within-population metric, is measured as allelic diversity (the number of unique alleles in a sampled population) and as mean heterozygosity (the average heterozygosity in a population). The finest level of genetic diversity is measured within the genotype itself as individual heterozygosity (the proportion of heterozygous loci sampled in an individual).
At present, the within-population levels of genotypic and genetic diversity have been evaluated among pure parental species and their hybrids [10,11] and among genotypes within an individual species [4,9]. In both cases, greater genetic variation among plants within a plot or stand is associated with greater biodiversity in associated communities that range from soil microbes to arthropods [2,10,[12][13][14][15]. Intraspecific variation in the communities of individual plant genotypes is also associated with heritable differences in community stability, resistance to environmental change, and community network structure [9,16]. These results are largely driven by differences in plant traits that make each plant unique in functional traits (i.e., the multivariate phenotype; [17]), which in turn results in different genotypes supporting different communities of organisms and ecosystem processes such as decomposition and nutrient cycling [5,18].
Elegant examples in trees come from species such as cottonwoods (Populus spp.) and eucalypts (Eucalyptus globulus), for which associated arthropod community diversity was dictated by tree genotypic diversity [11,19,20]. In common garden experiments, Evans et al. [15] found that growth rate was a heritable trait and that faster growing genotypes of P. angustifolia supported greater arthropod species richness, abundance, and Shannon diversity H' than slower growing genotypes. Similarly, Barker et al. [21] demonstrated that P. tremuloides size, phenology, and foliar chemistry influenced arthropod abundance and diversity. These studies reveal that tree productivity is positively correlated with community diversity metrics, which was observed at individual tree, genotype mean, and population mean levels. The similarity of results in evolutionarily disparate tree species suggests that the genetic diversity-begets biodiversity relationship may commonly occur. However, although this effect has been demonstrated at two different levels, the even finer level of genetic diversity remains to be fully explored, i.e., the genetic variation within an individual plant that can be quantified as heterozygosity. Bringing community effects of genetic diversity down to their root level-that of the foundation individual-would demonstrate the fundamental importance of genetic diversity to biodiversity, and reinforce the relevance of the individual to ecology and restoration.
Individual multilocus heterozygosity (IMH) at microsatellite or SNP loci has received considerable attention due to heterozygosity-fitness correlations, which have been discussed for over 40 years [22]. Heterozygosity-fitness correlations are broadly distributed across taxa in natural populations [23,24], but the signal is usually weak and is contextspecific [23,25]. Strong correlations can occur, but only when particular population structures or histories promote them, and as such, they are a characteristic of populations and not of species per se. For heterozygosity-fitness correlations to arise, there must be a proportion of consanguineous matings [26], population bottlenecks or genetic drift in small populations [27], or admixture or immigration [28]. These processes generate variation in inbreeding and a correlation in heterozygosity across loci (identity disequilibrium). Hence, heterozygosity-fitness correlations are not expected in large panmictic populations at equilibrium. Instead, they are primarily present when inbreeding, as a result of population characteristics, is pronounced and when it explains the variance in fitness [22]. Inbreeding depression, the reduction of performance in fitness-related traits [29], likely generates this variation in inbreeding levels and the relationship between heterozygosity and fitness. It is caused by the expression of deleterious recessive or partially recessive alleles due to mating between relatives in small populations or self-pollination in a mixed mating system. Heterozygosity-fitness correlations are only expected when many loci controlling a trait have directional dominance, such as for growth or life-history traits. A heterozygote excess is present when alleles of a neutral molecular marker locus are in linkage disequilibrium with different partially recessive deleterious variants that are masked by dominant alleles on the homologous chromosome (associative overdominance; [27]). Homozygotes, on the other hand, have lowered fitness because they are homozygous for the associated detrimental variants.
In the present study, we focused on a small, semi-arid riparian system that is likely to have the population structure and/or history necessary for heterozygosity-fitness correlations. We targeted trees that had been felled by beavers (Castor canadensis) and were resprouting from their stumps ( Figure 1). Beaver felling represents a stressful event for these trees, as felled trees lose nearly all of their aboveground biomass and may die. However, many resprout from the stump using reserves stored in the roots, essentially converting a mature tree to a juvenile state. We sought a trait that is both indicative of tree productivity and would affect a diverse arthropod community. A few key attributes of leaf quality appear to predictably drive generalist herbivore preferences, regardless of the complexities of plant-herbivore interactions [44]. We examined one of these, i.e., specific leaf area (SLA), which is consistently positively correlated with relative growth rate in young trees [39,45,46]. We also examined the phytochemistry of resprout growth, since defensive chemistry has been found to deter browsing mammals [47][48][49], and the degree of inbreeding has been found to be negatively correlated with defensive chemistry [50][51][52].
With these data, we hypothesized that (1) individual Fremont cottonwood trees harboring greater genetic diversity (IMH) would affect the diversity of their dependent arthropod communities. Because individual multilocus heterozygosity is an individual tree metric and not a population metric, herein we often refer to IMH in the singular (i.e., a tree's IMH). If substantiated, this IMH/arthropod diversity relationship would likely result from other more proximal factors that affect arthropods such as productivity. To explore this potential mechanism, we hypothesized that (2) greater genetic diversity (IMH) would be associated with greater productivity (i.e., the heterozygosity-productivity hypothesis), which has been shown to be attractive to arthropods [15]. We also hypothesized that (3) the IMH/arthropod relationship would not be present in the juvenile growth of trees not felled by beavers because they would not have been stimulated by herbivory. Additionally, we hypothesized that (4) greater genetic diversity (IMH) would be associated with phytochemical defenses (i.e., the phenolic glycoside salicortin; the heterozygosity-phytochemistry hypothesis), as beaver felling of trees could induce phytochemical changes in resprout twigs that may act to deter further browsing by beavers. Answers to these hypotheses are important because our observations appear to reveal changes in tree and community phenotypes resulting from biotic stress caused by beaver herbivory. In a time of great environmental change in the American Southwest, where an ongoing 19-year-long drought is considered to be the second worst in 1200 years [53], multiple stressors are likely to affect plant survival and the communities they support. Our research enhances the understanding of (a) how biotic stressors such as beavers interact with genetics to affect plant performance and (b) whether genetic diversity at the individual tree level has the potential to impact associated communities. Recent studies have shown that selection acting on the tree is also acting on its associated community of herbivores and can lead to community evolution [5]; therefore, it is important to consider the cascading effects of interactions at higher levels.

Study Site
This study took place on Arizona's Verde River, a perennial stream, and on two of its intermittently flowing tributaries, Dry Beaver Creek and Red Tank Draw. The maximum straight-line distance between sites was 21.1 km. P. fremontii is a dominant tree and the preferred food of beavers in this riparian system [54] and others [55]. Fifty trees were selected for the study (Verde, N = 14; Dry Beaver Creek, N = 26; Red Tank Draw, N = 10). Half of the selected trees at each site had been felled by beavers and were resprouting from their stumps, and the other half were unfelled. These were a set of the study trees of Durben et al. [54], who randomly selected a tree that had been felled by beaver and had resprouted from its base and then selected the nearest tree of the same size and at the same distance from the river that had not been felled by beaver to serve as a control. This pairing allowed standardization for resource availability (i.e., distance from the river, within 30 m) and tree size, as beavers can be sensitive to both of these variables. All necessary permits were obtained for the described study, which complied with all relevant USA regulations.

Arthropod Surveys and Twig Chemistry
In August 2009, we performed non-destructive visual arthropod surveys. Because juvenile and mature foliage can support different arthropod communities [56], to standardize for plant ontogeny, only the juvenile foliage at the base of unfelled trees was used to contrast with the juvenile foliage of resprouting beaver-felled trees. We sampled the same amount of foliage on resprouting beaver-felled trees and unfelled trees by following Figure 1. A Fremont cottonwood tree that was felled by beaver and subsequently resprouted (a) and cottonwood trees that were felled by beaver during the winter, that may or may not resprout in the spring (b). For scale, the cooler in (a) was 25 cm wide, and the trees in (b) were 25-30 cm in diameter.

Study Site
This study took place on Arizona's Verde River, a perennial stream, and on two of its intermittently flowing tributaries, Dry Beaver Creek and Red Tank Draw. The maximum straight-line distance between sites was 21.1 km. P. fremontii is a dominant tree and the preferred food of beavers in this riparian system [54] and others [55]. Fifty trees were selected for the study (Verde, N = 14; Dry Beaver Creek, N = 26; Red Tank Draw, N = 10). Half of the selected trees at each site had been felled by beavers and were resprouting from their stumps, and the other half were unfelled. These were a set of the study trees of Durben et al. [54], who randomly selected a tree that had been felled by beaver and had resprouted from its base and then selected the nearest tree of the same size and at the same distance from the river that had not been felled by beaver to serve as a control. This pairing allowed standardization for resource availability (i.e., distance from the river, within 30 m) and tree size, as beavers can be sensitive to both of these variables. All necessary permits were obtained for the described study, which complied with all relevant USA regulations.

Arthropod Surveys and Twig Chemistry
In August 2009, we performed non-destructive visual arthropod surveys. Because juvenile and mature foliage can support different arthropod communities [56], to standardize for plant ontogeny, only the juvenile foliage at the base of unfelled trees was used to contrast with the juvenile foliage of resprouting beaver-felled trees. We sampled the same amount of foliage on resprouting beaver-felled trees and unfelled trees by following the methods of Wimp et al. [57,58], by which we measured each branch until we reached a total of 20 mm of branch diameter to standardize for less growth on resprout trees compared to unfelled trees. Thus, both groups of trees represented random samples, and our standardizations across groups allowed them to be compared relative to their phytochemistry, genetics, and arthropod communities. We observed each set of branches for 20 min and identified arthropods as species or morphospecies within a family or genus. Unknown arthropods were identified via a reference collection at the Colorado Plateau Biodiversity Center at Northern Arizona University. We found 24 species in six orders (Lepidoptera, Hemiptera, Thysanoptera, Coleoptera, Hymenoptera, and Araneae).
We analyzed phytochemistry from twigs instead of leaves because beavers forage on the cambium layer of trees. We wished to assess whether beaver felling of trees might induce chemical defenses against further herbivory and other chemical changes of resprout growth (reported in Durben et al. [54]) and whether these might be associated with tree heterozygosity. We chose phytochemical traits that were previously found important: salicortin and total nitrogen (Martinsen et al. [59]). We quantified salicortin and nitrogen in 45 cm of terminal plant growth of the juvenile region of each tree according to the methods of Lindroth et al. [60]. We used dry ice to flash-freeze twigs in the field. Twigs were ground in a Wiley mill to pass through a 40-mesh screen. We employed high-performance thinlayer chromatography (HPTLC) to quantify salicortin [60]. We used a Thermo Finnigan Flash 1112 elemental analyzer (Thermo Finnigan, San Jose, CA, USA) to assess nitrogen [61]. Results are reported as % dry weight; additional methods are provided by Durben et al. [54].

Microsatellite Genotyping
We collected fresh leaves from each tree and immediately placed them in Ziploc bags containing calcium sulfate (Drierite; [62]). After desiccation, seven~0.6 cm discs of tissue were removed from a leaf from each tree using a standard paper punch, total genomic DNA was extracted after leaf homogenization using the DNeasy plant miniprep kit (QI-AGEN, Valencia, CA, USA), and DNA concentrations were estimated using a NanoDrop fluorometer (Thermo Fisher Scientific Inc., Waltham, MA, USA).

Productivity Measurement
We examined specific leaf area (SLA, cm 2 g −1 ; ratio between leaf area and dry leaf mass), a correlate of tree productivity [46,65,66], of resprouting beaver-felled trees at Dry Beaver Creek in August 2010. This site was selected over the others because trees had similar canopy cover (0-5%), their resprout growth was of similar age (estimated 2-3 years based upon bud scale scars) and height (mean 2.5 m), and they occupied a small area within the river channel (30 × 70 m), minimizing potentially confounding factors [67]. At Dry Beaver Creek, we also examined SLA of juvenile foliage of mature trees that were still standing at the time of leaf collection (N = 4); although we had planned to sample all of the standing study trees, few were unfelled when we returned to the study site. SLA was not assessed at the other sites because of their heterogeneity in all of the above traits, each of which can dramatically drive variation, even among leaves of the same individual [67]. Twenty leaves were collected per tree, as recommended by Cornelissen et al. [46], and were standardized by the 5th plastocron position on each resprouting branch. SLA was computed using WinFOLIA software for total leaf area measurements, followed by oven-drying (60 • C for 48 h) for each tree's total collected leaf weight.

Statistical Approach
We used GENHET [68], an R function that estimates 5 commonly used measures of individual multilocus heterozygosity. Because these metrics are highly correlated and provided similar results, here we report the proportion of heterozygous loci (pHt) per individual tree, which was recommended by Chapman et al. [23]. Shannon's Diversity Index [69], arthropod richness, and arthropod abundance were used to quantify leaf arthropod communities. To examine the relationship between tree heterozygosity and the arthropod metrics, we employed multiple regression followed by simple linear regression after residual, and quantile plots indicated that the data were of constant variance and normally distributed. Analysis of covariance (ANCOVA) was used to determine whether the regression slopes of felled vs. unfelled trees differed significantly. To assess the association between tree productivity, individual multilocus heterozygosity, and arthropod community diversity, we regressed SLA on IMH and the above arthropod metrics. We performed partial correlation with SLA as the control variable to determine whether the relationship between IMH and arthropod community diversity was influenced by SLA. As a quality control exercise to assess whether some loci contributed to the results more than others, we examined the effect at each locus by classifying heterozygosity and performing student t-tests with Shannon's Diversity Index.

Examination of Potentially Confounding Factors
We assessed whether population substructure, underlying tree relatedness, or hybridization was present, which could bias our results. Two analytical approaches were employed to examine population substructure. First, we used a Bayesian clustering approach (STRUCTURE 2.3.3, [70][71][72]) to test whether the Verde River, Dry Beaver Creek, and Red Tank Draw sampling sites are more than one genetic population. Ten replicate runs were performed (admixture model, correlated allele frequencies) for 1-5 genetic clusters (K), for 50,000 burn-in, followed by 200,000 Markov chain Monte Carlo (MCMC) iterations. Second, we split tree IMH into categories, i.e., High, Intermediate, and Low, then ran Analysis of Molecular Variance (AMOVA) (program ARLEQUIN 3.5, [73]) with the three groups as subpopulations to arrive at R ST and its significance. To examine underlying patterns of relatedness in the aforementioned groups, we estimated the mean relatedness [74] of pairs of trees from genotypic data, then tested for among-and within-group differences. Although introgression from narrowleaf cottonwoods (Populus angustifolia, the most geographically proximate Populus spp.) to Fremont is not common [75,76], we nonetheless employed program STRUCTURE 2.3.3 to ascertain whether any trees were indeed hybrids, using the settings mentioned above. We tested our study population with two populations of narrowleaf: trees from the Blue River, Arizona, using 6 microsatellite loci in common (N = 17; [77]), and trees from the San Miguel River, Colorado, using 8 loci in common (N = 8; F. Walker, unpublished data).

IMH Correlations for Felled Trees
In agreement with our hypothesis that individual tree heterozygosity affects the arthropod community, we found that a felled tree's IMH was positively correlated with its arthropod community diversity (Shannon's Diversity Index; R 2 = 0.56, p = 0.005; Figure 2a),

IMH Correlations for Felled Trees
In agreement with our hypothesis that individual tree heterozygosity affects the arthropod community, we found that a felled tree's IMH was positively correlated with its arthropod community diversity (Shannon's Diversity Index; R 2 = 0.56, p = 0.005; Figure  2a), species richness (R 2 = 0.49, p = 0.01; Figure 2b), and abundance (R 2 = 0.36, p = 0.04; Figure 2c).   In examining potential hypotheses that could mechanistically account for the above IMH/arthropod community findings, we found support for the IMH/productivity hypothesis. At Dry Beaver Creek, where it was tractable to measure productivity, multiple regression analysis indicated that for beaver-felled trees, a tree's arthropod community diversity (Shannon's Diversity Index; response variable) could be predicted by its IMH and SLA (factor variables) (R 2 = 0.52; ANOVA for model: F(2, 11) = 4.28, p = 0.05). A felled tree's IMH was positively correlated with its SLA (simple linear regression, R 2 = 0.45, p = 0.02; Figure 3a). Partial correlation revealed that the relationship between IMH and arthropod community diversity was dependent, in part, on SLA (correlation coefficient between IMH and Shannon's Diversity Index decreased from 0.72 to 0.57 when SLA was partialed out; p = 0.08).
In examining potential hypotheses that could mechanistically account for the above IMH/arthropod community findings, we found support for the IMH/productivity hypothesis. At Dry Beaver Creek, where it was tractable to measure productivity, multiple regression analysis indicated that for beaver-felled trees, a tree's arthropod community diversity (Shannon's Diversity Index; response variable) could be predicted by its IMH and SLA (factor variables) (R 2 = 0.52; ANOVA for model: F(2, 11) = 4.28, p = 0.05). A felled tree's IMH was positively correlated with its SLA (simple linear regression, R 2 = 0.45, p = 0.02; Figure 3a). Partial correlation revealed that the relationship between IMH and arthropod community diversity was dependent, in part, on SLA (correlation coefficient between IMH and Shannon's Diversity Index decreased from 0.72 to 0.57 when SLA was partialed out; p = 0.08).
A felled tree's IMH was also positively correlated with twig salicortin concentration. This finding supports the hypothesis that beaver felling of trees induces phytochemical changes in resprout growth that may be important for deterring future browsing by beavers (i.e., the heterozygosity-phytochemistry hypothesis). Figure 3b shows the relationship between IMH and salicortin level (p = 0.03). There was no relationship between IMH and nitrogen concentration (R 2 = 0.28, p = 0.12). The other nearby drainages had similar relationships between individual multilocus heterozygosity and arthropod communities. We combined data from Verde and Red Tank Draw (due to their small sample sizes) and found a positive non-significant trend for Shannon's Diversity Index and IMH (N = 12, R 2 = 0.15, p = 0.17), as well as for salicortin level and IMH (N = 9, R 2 = 0.34, p = 0.10), in resprouting beaver-felled trees. For all sites combined, arthropod community diversity (Figure 4; R 2 = 0.33, p = 0.004), richness (R 2 = 0.24, p = 0.01), and abundance (R 2 = 0.18, p = 0.04) significantly increased with individual multilocus heterozygosity in resprouting beaver-felled trees (N = 25). A felled tree's IMH was also positively correlated with twig salicortin concentration. This finding supports the hypothesis that beaver felling of trees induces phytochemical changes in resprout growth that may be important for deterring future browsing by beavers (i.e., the heterozygosity-phytochemistry hypothesis). Figure 3b shows the relationship between IMH and salicortin level (p = 0.03). There was no relationship between IMH and nitrogen concentration (R 2 = 0.28, p = 0.12).

Contrasts of Felled and Unfelled Trees
In agreement with our hypothesis that trees not felled by beavers would lack a link between heterozygosity and arthropods, we found no relationship between a tree's IMH and its arthropod community diversity (Figure 3, N = 23; p = 0.62, R 2 = 0.01), richness (p = 0.25, R 2 = 0.06), or abundance (p = 0.42, R 2 = 0.03). There was also no relationship between IMH and salicortin (N = 25; p = 0.15, R 2 = 0.08) or IMH and specific leaf area (N = 4; p = 0.91, R 2 = 0.009), although the sample size for the latter was small due to beavers subsequently felling trees.

Alternative Hypotheses
Differences between felled and unfelled trees were not due to confounding factors. There was no evidence of genetic substructuring between sites (K = 1 resulted in the highest log-likelihood) or between trees with different levels of heterozygosity (AMOVA, RST = 0.001, p = 0.61). Mean pairwise relatedness between heterozygosity classes did not differ (p = 0.24), nor was there evidence of hybridization (K = 2 resulted in the highest log-likelihoods for both datasets). Sites did not differ in tree heterozygosity (see below) or arthropod community diversity or richness (p = 0.70 and 0.40); it was thus appropriate to pool site data.

Locus Behavior
The microsatellite loci employed in this study were highly discriminating with regard to individual tree identity: the probability of two randomly selected trees sharing a multilocus genotype was 1.66 × 10 −6 , and the probability of full siblings doing so was 0.003 [78]. All loci conformed to Hardy-Weinberg expectations, and there was no evidence of . Individual multilocus heterozygosity of beaver-felled P. fremontii was associated with leaf arthropod community diversity at all sites. In Arizona's Verde River area, the relationship was significantly positive for resprouting beaver-felled trees, but not for trees that were unfelled (beaver-felled: y = 0.69 + 1.44x, R 2 = 0.33, p = 0.004; unfelled: y = 1.19 − 0.25x, R 2 = 0.01, p = 0.62). This difference between felled and unfelled trees was significant (ANCOVA interaction term: F(1, 48) = 4.89, p = 0.03).

Contrasts of Felled and Unfelled Trees
In agreement with our hypothesis that trees not felled by beavers would lack a link between heterozygosity and arthropods, we found no relationship between a tree's IMH and its arthropod community diversity (Figure 3, N = 23; p = 0.62, R 2 = 0.01), richness (p = 0.25, R 2 = 0.06), or abundance (p = 0.42, R 2 = 0.03). There was also no relationship between IMH and salicortin (N = 25; p = 0.15, R 2 = 0.08) or IMH and specific leaf area (N = 4; p = 0.91, R 2 = 0.009), although the sample size for the latter was small due to beavers subsequently felling trees.

Alternative Hypotheses
Differences between felled and unfelled trees were not due to confounding factors. There was no evidence of genetic substructuring between sites (K = 1 resulted in the highest log-likelihood) or between trees with different levels of heterozygosity (AMOVA, R ST = 0.001, p = 0.61). Mean pairwise relatedness between heterozygosity classes did not differ (p = 0.24), nor was there evidence of hybridization (K = 2 resulted in the highest log-likelihoods for both datasets). Sites did not differ in tree heterozygosity (see below) or arthropod community diversity or richness (p = 0.70 and 0.40); it was thus appropriate to pool site data.

Locus Behavior
The microsatellite loci employed in this study were highly discriminating with regard to individual tree identity: the probability of two randomly selected trees sharing a multilocus genotype was 1.66 × 10 −6 , and the probability of full siblings doing so was 0.003 [78]. All loci conformed to Hardy-Weinberg expectations, and there was no evidence of linkage disequilibrium between locus pairs (program GENEPOP 4.0.10; [79]). The sig-nificant results for beaver-felled trees were not driven by a strong effect at particular loci. The per locus relationship between heterozygosity and arthropod community diversity was examined because it yielded the strongest overall result. Although one locus was significant (p = 0.03), it was an order of magnitude less than the multilocus result; five other loci, although not significant, followed in the same direction, and none was in the opposite direction.

The IMH-Arthropod Correlation as It Relates to Productivity and Plant Defense
This work links genetic diversity (IMH) of the individual tree to community diversity. Arthropod diversity, as measured by Shannon's diversity index, and arthropod richness and abundance provided significant signals of the positive relationship with IMH, with the former being the strongest. In order to examine causality, the correlational patterns found in this study require experimental testing in a nearby common garden. A possible mechanism that may promote the arthropod-heterozygosity correlation that we examined is increased growth rate.
Since resprouting beaver-felled trees already have significant root systems, it is likely that they will allocate biomass to leaves, and thus any heterozygosity-driven differences in relative growth rate, such as SLA, will be readily apparent. A tradeoff often exists between root and leaf growth, with trees investing in one at the expense of the other [66]. In our study, higher SLA appeared to be attractive to arthropods. In many plant species, high SLA is correlated with high foliar nitrogen concentration, and both traits make leaves more attractive to herbivores [44,80,81]. Indeed, herbivores have been found to prefer high-SLA leaves in the field [44] and in the lab [67], and Martinsen et al. [59] found that beaver-felled P. angustifolia contained higher foliar nitrogen than unfelled control clones of the same genotypes. That the partial correlation results, showing a positive relationship between IMH and arthropod community diversity after accounting for joint covariance with SLA, approached significance (p = 0.08) may indicate that we have captured only one of multiple productivity measures that may be driving the relationship. We suggest that future studies include additional productivity estimates.
The significant correlation between IMH and salicortin may suggest that this phenolic glycoside increases following felling as a deterrent to further herbivory by beavers. It also suggests that heterozygosity facilitates this induction, although controlled experimental tests are required to draw this conclusion. Although no studies have examined beaver choice based on pre-felling salicortin levels, indirect evidence for the role of this phytochemical as a repellant was suggested by Durben et al. [54], who found higher concentrations of salicortin in resprout growth of beaver-felled P. fremontii than in juvenile growth of unfelled trees. Other mammalian herbivores avoid eating trembling aspen (Populus tremuloides) with leaves that have high concentrations of salicortin (North American porcupines, Erethizon dorsatum, [48]; elk, Cervus canadensis, [49]). Because we tested the phytochemistry of twigs that beavers feed upon and not of leaves (for which we quantified foliar arthropods), we do not infer responses of foliar arthropods to this IMH-salicortin relationship. However, Martinsen et al. [59] found that salicortin was the most abundant phenolic glycoside and that concentrations in the leaves were nearly two times greater in beaver-felled, resprouting trees (mostly complex backcross hybrids resembling P. angustifola with a few genotypes of P. fremontii and P. angustifolia) than in unfelled control trees of the same replicated genotypes. Similarly, they found that total nitrogen was nearly 20% greater in the leaves of resprout growth versus those of control trees. Further, Durben et al. [54] found these same salicortin and nitrogen relationships in the twigs of beaver-felled P. fremontii relative to those of unfelled trees. These findings suggest that twig and leaf salicortin could be positively correlated, although additional research is needed to confirm or reject this conclusion. We did not find a relationship between IMH and total nitrogen concentrations.
In plants, heterozygosity correlations have been found that can have clear relationships with productivity, fitness, environmental stress, and/or pathogens/herbivory. For example, the link between heterozygosity and fitness is believed to accrue with age in sexual tree species [82]; this may be the reason why long-lived species are most often observed to have an excess of heterozygotes [83]. For instance, in the American chestnut (Castanea dentata), heterozygotes had slightly higher growth rates, which over time produce a major excess of heterozygotes, and a clear heterozygote advantage as a consequence of the chestnut blight pandemic [43]. This pandemic has nearly eradicated sexual reproduction, allowing clonal heterozygous genotypes to dominate populations due to their growth advantage over the past 70 years. In aspen (Populus tremuloides), the growth rate increased with individual heterozygosity [42,83,84], and did so for males more than for females [41]. Cole et al. [39,41] found a tradeoff in which growth in aspens is positively related to heterozygosity and negatively related to condensed tannins. They suggested that episodes of pathogen or herbivore attacks alternated the selection for trees with high heterozygosity and high growth with that for trees with low heterozygosity and defensive tannin production. In mountain pine (Pinus uncinata Ram.), juveniles showed higher heterozygosity and growth rates, a relationship not seen for adult trees [85].

Heterozygosity, Stress, and a Semi-Arid System
In a variety of taxa, both microsatellites and allozymes show that the heterozygosityfitness correlations can be hidden in favorable conditions but become apparent when the organism experiences stress (eelgrass: [38]; common frog [86]; Perez's frog [87]). In pinyon pines (Pinus edulis), trees in high-stress abiotic environments (cinder fields) were more heterozygous than in reduced-stress environments (sandy-loam soils), and in the high-stress environment, trees genetically resistant to insect attack were more heterozygous than trees susceptible to insect attack [88]. Thus, it appears that abiotic and biotic stressors can reveal fitness-associated traits that also correlate with increased heterozygosity. Further, plant stress associated with climate change and reduced productivity has been demonstrated to have dramatic negative effects on a diverse arthropod community associated with pinyon pine in the Southwest. Stone et al. [89,90] found that for both juvenile and mature trees, arthropod richness and abundance on individual pinyon pine were positively correlated with the tree's radial trunk growth and that highly stressed trees with low productivity supported a very depauperate community relative to vigorously growing trees. This relationship was found to occur among trees within a stand, as well as across stands [91], and these community findings are in general agreement with the vigor hypothesis [92], according to which most herbivores prefer vigorously growing plants to feed upon because of their increased nutritional content. The present study extends this relationship to the underlying genetics (IMH) of individual plants.
In aspen, the effect of individual multilocus heterozygosity, as assessed by microsatellite markers, on growth rate was a simple linear relationship [42]. Because the association between heterozygosity and growth rate was additionally found in aspen studies employing isozymes [83,84], Cole et al. [42] postulated that both marker types yielded a signal of genome-wide inbreeding, in line with the general effects hypothesis for heterozygosityfitness correlations [93]. As with aspen, Fremont cottonwoods are dioecious obligate outcrossers. Despite this, some trees would be expected to express effects of inbreeding if sufficiently localized pollen and seed dispersal caused breeding between relatives (biparental inbreeding). This scenario is particularly plausible in riparian systems of the desert Southwest, such as intermittently flowing tributaries of the Verde River like Dry Beaver Creek, where relatively few trees inhabit a riparian corridor surrounded by expansive desert and where we found the strongest pattern. We propose that in this cottonwood system, the heterozygosity-productivity relationship may be a surrogate for a heterozygosity-fitness correlation and that this correlation is a consequence of variation in inbreeding, in accordance with the general effects hypothesis for heterozygosity-fitness correlations. Although inbreeding depression has been found to result in post-dispersal selection against inbred progeny and a corresponding increase in heterozygosity from juvenile to mature forest trees (e.g., [94]), in this study system it appears that mature trees with some degree of inbreeding were present. The new generation of juvenile growth following the new biotic stress of felling by beavers provides another opportunity for selection against inbred individuals.

Conclusions
In agreement with our hypothesis that interactions among foundation species have community-level effects, we found that genetic diversity (IMH) of individual resprouting beaver-felled trees was positively correlated with productivity, a phenolic glycoside (salicortin), and arthropod community richness, diversity, and abundance, whereas these correlations were absent for trees not felled by beavers. Thus, it appears that the tree heterozygosity/arthropod community/productivity/salicortin relationship only reveals itself in the presence of beavers that severely stress trees by removing nearly all of their aboveground biomass so that trees often die. This work provides a quartet of implications for heterozygosity-productivity (fitness) studies and conservation. First, IMH-fitness investigations should compare groups of organisms associated with populations experiencing different degrees of abiotic and biotic stress. Second, genetic diversity at the individual level can scale up to effect communities. Third, this effect is likely to be strongest where two or more foundation species interact, as with cottonwoods and beavers. Fourth, where conservation biologists have a choice of which source populations and individual genotypes to use in restoration, our results suggest that genetic diversity within individuals can positively affect biodiversity. Because these relationships are a product of biotic stress that occurred during a 19-year period of record drought years in the American Southwest that are projected to continue, it seems critical to be cognizant of how plant genetics including heterozygosity may interact with multiple stressors to affect plant survival and the communities they support.