Linking plant genes to insect communities: Identifying the genetic bases of plant traits and community composition

Genetic underpinnings of arthropod community distributions in Populus trichocarpa

Community genetics seeks to understand the mechanisms by which natural genetic variation in heritable host phenotypes can encompass assemblages of organisms such as bacteria, fungi, and many animals including arthropods. Prior studies that focused on plant genotypes have been unable to identify genes controlling community composition, a necessary step to predict ecosystem structure and function as underlying genes shift within plant populations. We surveyed arthropods within an association population of Populus trichocarpa in three common gardens to discover plant genes that contributed to arthropod community composition. We analyzed our surveys with traditional single-trait genome-wide association analysis (GWAS), multitrait GWAS, and functional networks built from a diverse set of plant phenotypes. Plant genotype was influential in structuring arthropod community composition among several garden sites. Candidate genes important for higher level organization of arthropod communities had broadly applicable functions, such as terpenoid biosynthesis and production of dsRNA binding proteins and protein kinases, which may be capable of targeting multiple arthropod species. We have demonstrated the ability to detect, in an uncontrolled environment, individual genes that are associated with the community assemblage of arthropods on a host plant, further enhancing our understanding of genetic mechanisms that impact ecosystem structure.

Genomic and transcriptomic analyses reveal polygenic architecture for ecologically important traits in aspen (Populus tremuloides Michx.)

Intraspecific genetic variation in foundation species such as aspen (Populus tremuloides Michx.) shapes their impact on forest structure and function. Identifying genes underlying ecologically important traits is key to understanding that impact. Previous studies, using single-locus genome-wide association (GWA) analyses to identify candidate genes, have identified fewer genes than anticipated for highly heritable quantitative traits. Mounting evidence suggests that polygenic control of quantitative traits is largely responsible for this "missing heritability" phenomenon. Our research characterized the genetic architecture of 30 ecologically important traits using a common garden of aspen through genomic and transcriptomic analyses. A multilocus association model revealed that most traits displayed a highly polygenic architecture, with most variation explained by loci with small effects (likely below the detection levels of single-locus GWA methods). Consistent with a polygenic architecture, our single-locus GWA analyses found only 38 significant SNPs in 22 genes across 15 traits. Next, we used differential expression analysis on a subset of aspen genets with divergent concentrations of salicinoid phenolic glycosides (key defense traits). This complementary method to traditional GWA discovered 1243 differentially expressed genes for a polygenic trait. Soft clustering analysis revealed three gene clusters (241 candidate genes) involved in secondary metabolite biosynthesis and regulation. Our work reveals that ecologically important traits governing higher-order community- and ecosystem-level attributes of a foundation forest tree species have complex underlying genetic structures and will require methods beyond traditional GWA analyses to unravel.

Genome-wide association study of aphid abundance highlights a locus affecting plant growth and flowering in Arabidopsis thaliana

Plant life-history traits, such as size and flowering, contribute to shaping variation in herbivore abundance. Although plant genes involved in physical and chemical traits have been well studied, less is known about the loci linking plant life-history traits and herbivore abundance. Here, we conducted a genome-wide association study (GWAS) of aphid abundance in a field population of Arabidopsis thaliana. This GWAS of aphid abundance detected a relatively rare but significant variant on the third chromosome of A. thaliana, which was also suggestively but non-significantly associated with the presence or absence of inflorescence. Out of candidate genes near this significant variant, a mutant of a ribosomal gene (AT3G13882) exhibited slower growth and later flowering than a wild type under laboratory conditions. A no-choice assay with the turnip aphid, Lipaphis erysimi, found that aphids were unable to successfully establish on the mutant. Our GWAS of aphid abundance unexpectedly found a locus affecting plant growth and flowering.

Environment and Genotype Influence Quantitative and Qualitative Variation in Condensed Tannins in Aspen

Condensed tannins (CTs) are abundant, ecologically-relevant secondary metabolites in many plants, which respond to variables associated with anthropogenic environmental change. While many studies have reported how genetic and environmental factors affect CT concentrations, few have explored how they influence CT molecular structure. Here, using trembling aspen (Populus tremuloides) as a model organism, we report how foliar CT concentrations, polymer sizes, representation of procyanidins and prodelphinidins, and stereochemistry vary in response to changes in air temperature (warming and freeze damage), air composition (elevated CO2 and O3), soil quality (nutrients and microbiome), and herbivory (mammal and lepidopteran). Use of multiple aspen genotypes enabled assessment of genetic influences on aspen CTs. CT concentration and composition were analyzed by thiolysis-ultra high performance liquid chromatography/mass spectrometry in archived leaf samples from prior experiments. All environmental variables explored except for soil microbiome influenced both CT quantity and quality, with climate factors appearing to have larger effect magnitudes than herbivory. Climate, soil, and herbivory effects varied among genotypes, while air composition effects were consistent across genotypes. Considering that CT properties (concentrations and molecular structures) mediate functions at the organismal through ecosystem scales, intraspecific variation in responses of CT properties to environmental factors could provide a pathway through which environmental change exerts selective pressure on Populus populations. Future studies are needed to identify the molecular-level mechanisms by which environmental factors influence CT concentrations and structures, and to establish their ecological and evolutionary significance.

Effects of condensed tannins on behavior and performance of a specialist aphid on aspen

Genes involved in plant defences against herbivores and pathogens are often highly polymorphic. This is a putative sign that balancing selection may have operated reciprocally on the hosts and their herbivores. Spatial and temporal variations (for example, in soil nutrients and the plants' ontogenetic development) may also modulate resistance traits, and thus selection pressures, but have been largely overlooked in theories of plant defences. Important elements of defences in Populus tremula (hereafter aspen) are phenolic compounds, including condensed tannins (CTs). Concentrations of CTs vary considerably with both variations in external factors and time, but they are also believed to provide genotype-dependent resistance, mainly against chewing herbivores and pathogens. However, evidence of their contributions to resistance is sparse. Detailed studies of co-evolved plant-herbivore associations could provide valuable insights into these contributions. Therefore, we examined correlations between CT levels in aspen leaves and both the feeding behavior and reproduction of the specialist aspen leaf aphid (Chaitophorus tremulae) in varied conditions. We found that xylem sap intake and probing difficulties were higher on genotypes with high-CT concentrations. However, aphids engaged in more nonprobing activities on low-CT genotypes, indicating that CTs were not the only defence traits involved. Thus, high-CT genotypes were not necessarily more resistant than low-CT genotypes, but aphid reproduction was generally negatively correlated with local CT accumulation. Genotype-specific resistance ranking also depended on the experimental conditions. These results support the hypothesis that growth conditions may affect selection pressures mediated by aphids in accordance with balancing selection theory.

Genetic divergence along a climate gradient shapes chemical plasticity of a foundation tree species to both changing climate and herbivore damage

Climate change is threatening the persistence of many tree species via independent and interactive effects on abiotic and biotic conditions. In addition, changes in temperature, precipitation, and insect attacks can alter the traits of these trees, disrupting communities and ecosystems. For foundation species such as Populus, phytochemical traits are key mechanisms linking trees with their environment and are likely jointly determined by interactive effects of genetic divergence and variable environments throughout their geographic range. Using reciprocal Fremont cottonwood (Populus fremontii) common gardens along a steep climatic gradient, we explored how environment (garden climate and simulated herbivore damage) and genetics (tree provenance and genotype) affect both foliar chemical traits and the plasticity of these traits. We found that (1) Constitutive and plastic chemical responses to changes in garden climate and damage varied among defense compounds, structural compounds, and leaf nitrogen. (2) For both defense and structural compounds, plastic responses to different garden climates depended on the climate in which a population or genotype originated. Specifically, trees originating from cool provenances showed higher defense plasticity in response to climate changes than trees from warmer provenances. (3) Trees from cool provenances growing in cool garden conditions expressed the lowest constitutive defense levels but the strongest induced (plastic) defenses in response to damage. (4) The combination of hot garden conditions and simulated herbivory switched the strategy used by these genotypes, increasing constitutive defenses but erasing the capacity for induction after damage. Because Fremont cottonwood chemistry plays a major role in shaping riparian communities and ecosystems, the effects of changes in phytochemical traits can be wide reaching. As the southwestern US is confronted with warming temperatures and insect outbreaks, these results improve our capacity to predict ecosystem consequences of climate change and inform selection of tree genotypes for conservation and restoration purposes.

Intraspecific variation in plant economic traits predicts trembling aspen resistance to a generalist insect herbivore

Patterns of trait expression within some plant species have recently been shown to align with the leaf economics spectrum paradigm. Resistance to herbivores is also expected to covary with leaf economics traits. We selected 36 mature Populus tremuloides genotypes in a common garden to assess whether aspen leaf economics patterns follow those observed among species globally. We also evaluated leaf economics strategies in the context of insect resistance by conducting bioassays to determine the effects of plant traits on preference and performance of Lymantria dispar. We found that: (1) intraspecific trait patterns of P. tremuloides parallel those exhibited by the interspecific leaf economics spectrum, (2) herbivores preferred leaves from genotypes with resource-acquisitive strategies, and (3) herbivores also performed best on genotypes with resource-acquisitive strategies. We conclude that a leaf economics spectrum that incorporates defense traits is a useful tool for explaining intraspecific patterns of variation in plant strategies, including resistance to herbivores.

Polyploidy and growth-defense tradeoffs in natural populations of western quaking Aspen

Polyploidy, the expression of more than two sets of chromosomes, is common in plants, and is thought to influence plant trait expression and drive plant species evolution. The degree to which polyploidy influences interactions among physiological processes such as growth and defense in natural populations through its effect on phenotypic variability is poorly understood. We link broad plant genotypic features (including polyploidy) to phenotypic expression of growth and chemical defense in natural populations of quaking aspen (Populus tremuloides) to examine patterns in resource allocation that might drive growth-defense tradeoffs. Quaking aspen are capable of rapid growth, and are also a primary food plant for a large range of herbivores, including insects and ungulates. While often diploid, aspen can exhibit polyploidy as triploid clones. We tested for the effect of genotype, cytotype (ploidy level, divided between diploids and triploids), and ramet age on relationships between growth and leaf chemistry across natural aspen clones in northern Utah. Substantial genotype variability in growth and leaf chemistry occurred across both cytotypes. Phenolic glycosides, but not condensed tannins, were negatively related to growth. Ramet age was also negatively related to growth. Phenolic glycosides were negatively related to condensed tannins, but only for the diploid clones. Triploid clones exhibited ~ 20% higher levels of phenolic glycosides than diploids. Growth in quaking aspen was likely sacrificed for the production of phenolic glycosides. Our study underscores the importance of considering polyploidy, genetic variability, and ramet age in understanding growth-defense tradeoffs in natural populations of clonal organisms, such as quaking aspen.

What Drives Caterpillar Guilds on a Tree: Enemy Pressure, Leaf or Tree Growth, Genetic Traits, or Phylogenetic Neighbourhood?

Communities of herbivorous insects on individual host trees may be driven by processes ranging from ongoing development via recent microevolution to ancient phylogeny, but the relative importance of these processes and whether they operate via trophic interactions or herbivore movement remains unknown. We determined the leaf phenology, trunk diameter, genotype, and neighbourhood of sessile oak trees (Quercus petraea), and sampled their caterpillar communities. We found that leaf development across a time period of days related to free-living caterpillars, which disappeared with leaf age. Tree growth across decades is related to increased parasitism rate and diversity of herbivores. The microevolution of oak trees across millennia is related to the abundance of leaf-mining casebearers, which is higher on more homozygous oaks. However, oak genome size was not important for any guild. In contrast to most previous studies, the phylogenetic distance of oaks from their neighbours measured in millions of years was associated with higher abundances of entire caterpillar guilds. Furthermore, on trees surrounded by only distantly related tree species, parasitism tended to be lower. Lower parasitism, in turn, was associated with higher abundances of codominant caterpillar species. Neighbourhoods and traits of trees were also related to community composition and diversity, but not to the average wingspans or specialization of species, consistent with the assembly of herbivore communities being driven by leaf traits and parasitism pressure on trees rather than by insect movement among trees. However, movement in rarer species may be responsible for concentration effects in more phylogenetically distant neighbourhoods. Overall, we suggest that the assembly of insects on a tree is mostly driven by trophic interactions controlled by a mosaic of processes playing out over very different time scales. Comparisons with the literature further suggest that, for oak trees, the consequences of growing amongst distantly related tree species may depend on factors such as geographic region and tree age.

A keystone gene underlies the persistence of an experimental food web

Genes encode information that determines an organism's fitness. Yet we know little about whether genes of one species influence the persistence of interacting species in an ecological community. Here, we experimentally tested the effect of three plant defense genes on the persistence of an insect food web and found that a single allele at a single gene promoted coexistence by increasing plant growth rate, which in turn increased the intrinsic growth rates of species across multiple trophic levels. Our discovery of a "keystone gene" illustrates the need to bridge between biological scales, from genes to ecosystems, to understand community persistence.

Growth-defense trade-offs shape population genetic composition in an iconic forest tree species

All organisms experience fundamental conflicts between divergent metabolic processes. In plants, a pivotal conflict occurs between allocation to growth, which accelerates resource acquisition, and to defense, which protects existing tissue against herbivory. Trade-offs between growth and defense traits are not universally observed, and a central prediction of plant evolutionary ecology is that context-dependence of these trade-offs contributes to the maintenance of intraspecific variation in defense [Züst and Agrawal, Annu. Rev. Plant Biol., 68, 513-534 (2017)]. This prediction has rarely been tested, however, and the evolutionary consequences of growth-defense trade-offs in different environments are poorly understood, especially in long-lived species [Cipollini et al., Annual Plant Reviews (Wiley, 2014), pp. 263-307]. Here we show that intraspecific trait trade-offs, even when fixed across divergent environments, interact with competition to drive natural selection of tree genotypes corresponding to their growth-defense phenotypes. Our results show that a functional trait trade-off, when coupled with environmental variation, causes real-time divergence in the genetic architecture of tree populations in an experimental setting. Specifically, competitive selection for faster growth resulted in dominance by fast-growing tree genotypes that were poorly defended against natural enemies. This outcome is a signature example of eco-evolutionary dynamics: Competitive interactions affected microevolutionary trajectories on a timescale relevant to subsequent ecological interactions [Brunner et al., Funct. Ecol. 33, 7-12 (2019)]. Eco-evolutionary drivers of tree growth and defense are thus critical to stand-level trait variation, which structures communities and ecosystems over expansive spatiotemporal scales.

Genetic and plastic rewiring of food webs under climate change

Climate change is altering ecological and evolutionary processes across biological scales. These simultaneous effects of climate change pose a major challenge for predicting the future state of populations, communities and ecosystems. This challenge is further exacerbated by the current lack of integration of research focused on these different scales.

We propose that integrating the fields of quantitative genetics and food web ecology will reveal new insights on how climate change may reorganize biodiversity across levels of organization. This is because quantitative genetics links the genotypes of individuals to population‐level phenotypic variation due to genetic (G), environmental (E) and gene‐by‐environment (G × E) factors. Food web ecology, on the other hand, links population‐level phenotypes to the structure and dynamics of communities and ecosystems.

We synthesize data and theory across these fields and find evidence that genetic (G) and plastic (E and G × E) phenotypic variation within populations will change in magnitude under new climates in predictable ways. We then show how changes in these sources of phenotypic variation can rewire food webs by altering the number and strength of species interactions, with consequences for ecosystem resilience. We also find evidence suggesting there are predictable asymmetries in genetic and plastic trait variation across trophic levels, which set the pace for phenotypic change and food web responses to climate change. Advances in genomics now make it possible to partition G, E and G × E phenotypic variation in natural populations, allowing tests of the hypotheses we propose.

By synthesizing advances in quantitative genetics and food web ecology, we provide testable predictions for how the structure and dynamics of biodiversity will respond to climate change.

Growing up aspen: ontogeny and trade-offs shape growth, defence and reproduction in a foundation species

BACKGROUND AND AIMS

Intraspecific variation in foundation species of forest ecosystems can shape community and ecosystem properties, particularly when that variation has a genetic basis. Traits mediating interactions with other species are predicted by simple allocation models to follow ontogenetic patterns that are rarely studied in trees. The aim of this research was to identify the roles of genotype, ontogeny and genotypic trade-offs shaping growth, defence and reproduction in aspen.

METHODS

We established a common garden replicating >500 aspen genets in Wisconsin, USA. Trees were measured through the juvenile period into the onset of reproduction, for growth, defence chemistry (phenolic glycosides and condensed tannins), nitrogen, extrafloral nectaries, leaf morphology (specific leaf area), flower production and foliar herbivory and disease. We also assayed the TOZ19 sex marker and heterozygosity at ten microsatellite loci.

KEY RESULTS

We found high levels of genotypic variation for all traits, and high heritabilities for both the traits and their ontogenetic trajectories. Ontogeny strongly shaped intraspecific variation, and trade-offs among growth, defence and reproduction supported some predictions while contradicting others. Both direct resistance (chemical defence) and indirect defence (extrafloral nectaries) declined during the juvenile stage, prior to the onset of reproduction. Reproduction was higher in trees that were larger, male and had higher individual heterozygosity. Growth was diminished by genotypic allocation to both direct and indirect defence as well as to reproduction, but we found no evidence of trade-offs between defence and reproduction.

CONCLUSIONS

Key traits affecting the ecological communities of aspen have high levels of genotypic variation and heritability, strong patterns of ontogeny and clear trade-offs among growth, defence and reproduction. The architecture of aspen's community genetics - its ontogeny, trade-offs and especially its great variability - is shaped by both its broad range and the diverse community of associates, and in turn further fosters that diversity.

Consistent community genetic effects in the context of strong environmental and temporal variation in Eucalyptus

Provenance translocations of tree species are promoted in forestry, conservation, and restoration in response to global climate change. While this option is driven by adaptive considerations, less is known of the effects translocations can have on dependent communities. We investigated the relative importance and consistency of extended genetic effects in Eucalyptus using two species-E. globulus and E. pauciflora. In E. globulus, the dependent arthropod and pathogen canopy communities were quantified based on the abundance of 49 symptoms from 722 progeny from 13 geographic sub-races across 2 common gardens. For E. pauciflora, 6 symptoms were quantified over 2 years from 238 progeny from 16 provenances across 2 common gardens. Genetic effects significantly influenced communities in both species. However, site and year effects outweighed genetic effects with site explaining approximately 3 times the variation in community traits in E. globulus and site and year explaining approximately 6 times the variation in E. pauciflora. While the genetic effect interaction terms were significant in some community traits, broad trends in community traits associated with variation in home-site latitude for E. globulus and home-site altitude for E. pauciflora were evident. These broad-scale trends were consistent with patterns of adaptive differentiation within each species, suggesting there may be extended consequences of local adaptation. While small in comparison to site and year, the consistency of genetic effects highlights the importance of provenance choice in tree species, such as Eucalyptus, as adaptive divergence among provenances may have significant long-term effects on biotic communities.

Proteome changes in pepper (Capsicum annuum L.) leaves induced by the green peach aphid (Myzus persicae Sulzer)

BACKGROUND

Aphid attack induces defense responses in plants activating several signaling cascades that led to the production of toxic, repellent or antinutritive compounds and the consequent reorganization of the plant primary metabolism. Pepper (Capsicum annuum L.) leaf proteomic response against Myzus persicae (Sulzer) has been investigated and analyzed by LC-MS/MS coupled with bioinformatics tools.

RESULTS

Infestation with an initially low density (20 aphids/plant) of aphids restricted to a single leaf taking advantage of clip cages resulted in 6 differentially expressed proteins relative to control leaves (3 proteins at 2 days post-infestation and 3 proteins at 4 days post-infestation). Conversely, when plants were infested with a high density of infestation (200 aphids/plant) 140 proteins resulted differentially expressed relative to control leaves (97 proteins at 2 days post-infestation, 112 proteins at 4 days post-infestation and 105 proteins at 7 days post-infestation). The majority of proteins altered by aphid attack were involved in photosynthesis and photorespiration, oxidative stress, translation, protein folding and degradation and amino acid metabolism. Other proteins identified were involved in lipid, carbohydrate and hormone metabolism, transcription, transport, energy production and cell organization. However proteins directly involved in defense were scarce and were mostly downregulated in response to aphids.

CONCLUSIONS

The unexpectedly very low number of regulated proteins found in the experiment with a low aphid density suggests an active mitigation of plant defensive response by aphids or alternatively an aphid strategy to remain undetected by the plant. Under a high density of aphids, pepper leaf proteome however changed significantly revealing nearly all routes of plant primary metabolism being altered. Photosynthesis was so far the process with the highest number of proteins being regulated by the presence of aphids. In general, at short times of infestation (2 days) most of the altered proteins were upregulated. However, at longer times of infestation (7 days) the protein downregulation prevailed. Proteins involved in plant defense and in hormone signaling were scarce and mostly downregulated.

Putting the genes into community genetics

As part of the long-term fusion of evolutionary biology and ecology (Ford, 1964), the field of community genetics has made tremendous progress in describing the impacts of plant genetic variation on community and ecosystem processes. In the "genes-to-ecosystems" framework (Whitham et al., 2003), genetically based traits of plant species have ecological consequences, but previous studies have not identified specific plant genes responsible for community phenotypes. The study by Barker et al. (2019) in this issue of Molecular Ecology uses an impressive common garden experiment of trembling aspen (Figure 1) to test for the genetic basis of tree traits that shape the insect community composition. Using a Genome-Wide Association Study (GWAS), they found that genomic regions associated with phytochemical traits best explain variation in herbivore community composition, and identified specific genes associated with different types of leaf-modifying herbivores and ants. This is one of the first studies to identify candidate genes underlying the heritable plant traits that explain patterns of insect biodiversity.