- Research article
- Open Access
Plant trichomes and a single gene GLABRA1 contribute to insect community composition on field-grown Arabidopsis thaliana
BMC Plant Biologyvolume 19, Article number: 163 (2019)
Genetic variation in plants alters insect abundance and community structure in the field; however, little is known about the importance of a single gene among diverse plant genotypes. In this context, Arabidopsis trichomes provide an excellent system to discern the roles of natural variation and a key gene, GLABRA1, in shaping insect communities. In this study, we transplanted two independent glabrous mutants (gl1–1 and gl1–2) and 17 natural accessions of Arabidopsis thaliana to two localities in Switzerland and Japan.
Fifteen insect species inhabited the plant accessions, with the insect community composition significantly attributed to variations among plant accessions. The total abundance of leaf-chewing herbivores was negatively correlated with trichome density at both field sites, while glucosinolates had variable effects on leaf chewers between the sites. Interestingly, there was a parallel tendency for the abundance of leaf chewers to be higher on gl1–1 and gl1–2 than on their different parental accessions, Ler-1 and Col-0, respectively. Furthermore, the loss of function in the GLABRA1 gene significantly decreased the resistance of plants to the two predominant chewers; flea beetles and turnip sawflies.
Overall, our results indicate that insect community composition significantly varies among A. thaliana accessions across two distant field sites, with GLABRA1 playing a key role in altering the abundance of leaf-chewing herbivores. Given that such a trichome variation is widely observed in Brassicaceae plants, the present study exemplifies the community-wide effect of a single plant gene on crucifer-feeding insects in the field.
Plants develop various resistance traits, such as spines and toxins, to deter herbivory . A growing number of studies on community genetics has revealed that genetic variation in plant resistance traits exerts cascading effects on insect abundance and community composition [2,3,4,5]. These insect indices projected on individual plants, called extended phenotype , can be explained by variation among plant genotypes [6,7,8]. Some researchers have reported the association of particular genetic polymorphisms with leaf damage [9, 10], insect abundance [3, 11], and community composition  in the field. In comparison, other studies have focused on how single genes affect the insect community using transformed plants [12, 13]. These lines of evidence from diverse plant species suggest that both quantitative genetic variation and single genes contribute to the community genetics of plant-insect interactions.
Arabidopsis thaliana (L.) Heynh. (Brassicaceae) is well-studied as a model system of Brassicaceae-insect interactions , within which intensive genomic and phenotypic information is available for the world-wide collection of natural accessions . In Arabidopsis-herbivore interactions, plant trichomes (epidermal hairs) function as a mechanical barrier against feeding and oviposition by insect herbivores [11, 16,17,18]. Glucosinolates (GSLs) are major secondary metabolites of Brassicales that act as toxic chemicals against generalists [19, 20], but can be detoxified by specialist herbivores [14, 20, 21]. For example, previous studies on A. thaliana focused on how these physical and chemical traits confer resistance against specific herbivore species, such as the cabbage butterfly Pieris rapae [22, 23], the diamond back moth Plutella xylostella [16, 20] and the green peach aphid Myzus persicae [24, 25]. However, knowledge remains limited about (i) how many insect species occupy A. thaliana in the field, (ii) whether plant defense traits contribute to insect abundance and community composition and (iii) the host genes that are responsible for community members and overall community composition.
Complex genetic mechanisms underlie anti-herbivore defenses when plants are exposed to multiple biotic and abiotic stresses [26, 27]; consequently, the phenotype under constant laboratory conditions might not be adequate for understanding how genes function in the field [28,29,30,31,32]. In this concept of in natura study [28, 30,31,32], Arabidopsis thaliana provides an ideal opportunity to test the effects of single plant genes on herbivores because its molecular mechanisms of defense have been studied using a variety of accessions, with respect to trichomes [33,34,35,36] and secondary metabolites [20, 36, 37]. For example, loss of function mutations in a few transcription factor genes (including GLABRA1 (GL1), GLABRA2, GLABRA3, and TRANSPARENT TESTA GLABRA1) result in glabrous phenotypes in A. thaliana [38,39,40,41,42] and related species [9, 43, 44]. While the loss of function of the latter three genes results in pleiotropic defects in root hairs, the loss of function of GL1 does not affect root hairs, due to the subfunctionalization of GL1 and its homolog WEREWOLF [40,41,42, 45]. Laboratory experiments have shown that loss of function in GL1 decreases resistance against leaf-chewing herbivores  and improves plant growth by saving on the cost of defense [17, 25]. However, these genetic effects remain unexplored in the field.
Common garden experiments using single-gene mutants provide a powerful tool to determine the causal links between a particular gene and its phenotypes [46,47,48,49]. In this study, we transplanted two glabrous mutants and 17 natural accessions of A. thaliana. In particular, we focused on gl1–1 and gl1–2 accessions, of which the former is a null trichome mutant derived from the Ler-1 accession and the latter is a hypomorphic mutant from the Col-0 accession [38, 45]. In addition, the natural accessions were selected to cover variations in trichome density and GSLs content [35, 37, 43]. Common garden experiments with these plants were conducted at two field sites in Switzerland and Japan, to identify common patterns between the two insect communities. Three specific questions were addressed: (i) is there significant variation in herbivore abundance and community composition among the A. thaliana accessions; (ii) which plant traits (physical, chemical or other life-history traits) influence herbivore abundance and community composition; (iii) does the loss of function of a single gene, GL1, alter insect abundance and community composition?
Plant materials and defense traits
Arabidopsis thaliana (L.) Heynh., commonly known as thale cress or mouse-ear cress, is an annual weed native to Eurasia and Africa and naturalized in North America and East Asia . This species is predominantly self-fertilizing  and when plants are collected from wild populations or when mutants are isolated by mutagenesis, selfed seeds can be maintained as an inbred line called an “accession”. Weak dormancy and early-flowering accessions, such as Col-0 and Ler-1 [51, 52], form both the spring and summer cohort owing to their rapid life-cycles . The spring cohort flowers and sets seeds in spring, and the summer cohort germinates in early summer and flowers in autumn . The accessions with strong dormancy, such as Cvi-0 and Shahdara , pass the summer as seeds. The accessions with a strong dormancy and late-flowering phenotype, such as Kas-2, are predominantly winter-annuals with only one generation within a calendar year . These different life-cycles of A. thaliana accessions depend on the level of seed dormancy, which can be attributed to the allelic status of the DELAY OF GERMINATION1 (DOG1) and DOG6 genes [51, 52, 54] and the duration to flower development, which is determined by FRIGIDA, FLOWERING LOCUS C and several other genes [54, 55]. In wild populations within Europe, generalist slugs and seed weevils feed on A. thaliana during late spring, while more diverse herbivores, such as Phyllotreta beetles, green peach aphids Myzus persicae and diamondback moths Plutella xylostella, occur during summer [56, 57]. In wild populations near our study locations, we discovered A. thaliana plants during early summer (Fig. 1), though the population size seemed smaller than those that overwintered. We also observed flowering and vegetative A. thaliana co-occurring during early summer and the plants displayed leaf damage from insect herbivores (Fig. 1b). Therefore, to investigate the diverse herbivores on A. thaliana we simulated the summer cohort using accessions with different life-cycles and defense traits.
To cover the wide variation in trichome density (physical defense) and GSL accumulation (chemical defense) with early- and late-flowering cycles, we selected 17 natural accessions and two glabrous mutants (Table 1). The natural accessions selected in this study should represent world-wide genetic variation, because the genome-wide pairwise genetic distance was 5.7% in median, which is comparable to that of all accessions analyzed by the 1001 Genome Consortium . These 17 accessions include both early- and late-flowering accessions (e.g., Col-0 and Kas-2 analyzed by Taylor et al. ), such that the flowering time under a long-day laboratory conditions ranged from 23 (Ws-2 accession) to 92 days (Br-0 accession) . To examine the effects of plant life-history traits on insect community composition, we measured and incorporated the plant size and presence of flowering stems (see ‘Common garden experiment’ and ‘Statistical analysis’ below).
To test the functional advantages of the GL1 gene in producing trichomes, we added two glabrous mutants, gl1–1 and gl1–2, to the set of natural accessions (Table 1). The former mutant, gl1–1, has the background of Ler accession with a null mutation due to a 6.5-kb deletion on GL1 and lacks leaf surface trichomes. The latter, gl1–2, has the background of Col accession with the deletion of 27 amino acids induced by X-ray radiation, showing a hypomorphic mutation with a lower density of trichomes on leaf surfaces [38, 45]. Out of the 17 natural accessions, Br-0 and C24 have no or few trichomes due to a frameshift mutation and one amino acid change in the myb DNA binding domain of GL1, respectively . We compiled the data on leaf trichome density (no./cm2) from the GWA-portal (https://gwas.gmi.oeaw.ac.at/: ).
All the natural accessions were included in previous quantitative genetic studies of GSL. Seven were used as the parental genotypes of recombinant inbred lines (e.g., Col × Ler and Cvi × Ler ; Bay × Sha ; Kas × Tsu ) and the other accessions were used in genome-wide association mapping . To test whether genetic potentials in GSL profiles explain herbivory rates, we used the data in Chan et al.  on 21 GSLs of 96 A. thaliana accessions. They used a mature leaf at 35-days post germination from a plant grown under short-day laboratory conditions without herbivory. As they performed two trials to quantify GSL, we used the average GSL contents (nmol/mg flesh weight). We focused on variation in aliphatic GSLs and their chain-length, because these parameters play a major role in preventing above-ground herbivory [48, 60, 61]. Regarding the data of Chan et al. , we applied a principal component analysis (PCA) to total C3-, C4-, C5-, C7-, and C8- aliphatic GSLs. The first and second principal components explained 44 and 33% variation in the GSL profiles among our 17 accessions, respectively (Additional file 1); therefore, these two components were used in our statistical analyses.
Common garden experiment
We used the experimental gardens of the University of Zurich at Irchel campus (Zurich, Switzerland: 47° 23′ N, 8° 33′ E, alt. ca. 500 m) and the Center for Ecological Research, Kyoto University (Otsu, Japan: 35° 06′ N, 134° 56′ E, alt. ca. 200 m) (Fig. 2). The Zurich site is close to a deciduous forest and the surroundings of the common garden are covered with concrete tiles to prevent weeds. The Otsu site is a suburb of cultivated fields and the ground of the study site is covered with short grasses. In the Otsu site, the grass weeds were mown, and the surroundings were covered with agricultural sheets before the experiment. No large Brassica cultivars were grown at either site during early summer. Average air temperature and total precipitation was 19 °C and 198 mm in Zurich (during July 2016; MeteoSwiss, http://www.meteoswiss.admin.ch/home.html) and 22 °C and 321 mm in Otsu (during June 2016; Japan Meteorological Agency, http://www.jma.go.jp/jma/index.html).
We prepared 10 replicates of 19 accessions (= 190 plants in total) for each experiment. To minimize unnecessary variation due to competition and consequent size bias, experimental plants were grown in separate pots under uniform conditions before they were exposed to the field environment in the common garden. We cultivated plants using mixed soils of agricultural composts (Profi Substrat Classic CL ED73, Einheitserde Co. in Switzerland; MetroMix 350, SunGro Co. in Japan) and perlites with a compost to perlite ratio of 3:1 litter volume. No additional fertilizers were supplied because the agricultural soils contain fertilizers. Seeds were sown on the soil and stratified under constant dark conditions at 4–5 °C air temperature for a week. Plants were then grown under short-day conditions (8 h:16 h light:dark [L:D], 20 °C air temperature, and 60% relative humidity) for 1 month to prevent flowering before the field experiment. The plant positions were rotated every week to minimize growth bias from light conditions. Each individual plant was moved to a plastic pot (7.5-cm diameter, with 6.5-cm depth in Japan; 6.0 × 6.0 × 6.0 cm in Switzerland) and acclimated for 3 days in a shaded outdoor place before the field experiments. The potted plants were randomly placed among three blocks in each common garden: 68, 69, and 53 plants were assigned to each block in Zurich; and 76, 76, and 38 plants were assigned to each block in Otsu. The potted plants were set on water-permeable plastic sheets in a checkered manner within the blocks, without being embedded in the ground. The blocks were more than 1.0 m apart. These experiments were conducted from June 18 to July 1, 2016 in Otsu and from July 13 to August 3, 2016 in Zurich. Plants were watered every 3 days in Otsu and every day in Zurich.
Insects and herbivorous collembola on individual plants were visually counted every 2–3 days. These species were identified ocularly with a magnifying glass. Dwelling traces and mummified aphids were also counted as a proxy of the number of leaf miners and parasitoid wasps, respectively. Eggs, larvae, and adults were counted for all species, as long as they could be observed by the naked eye. The abundance of each species was evaluated by the cumulative number of individuals over the experimental period to reflect the herbivory load on the plants . Small holes made by flea beetles were counted at the Zurich site and the maximum number throughout the experiment was used as an indicator of damage by flea beetles; however, this phenotyping was difficult in Japan, due to heavier and simultaneous infestation by sawflies. All counting was conducted by a single observer during the daytime (08:00–17:00). Early-flowering accessions of A. thaliana start reproduction under long-day field conditions and terminate their life-cycle within a month. Thus, to avoid the disturbance during analysis due to plant senescence, the duration of field observations was for 3 weeks, after the beginning of the field experiment. We rotated the order of census plants to avoid any bias of census time. The mortality of plants during insect monitoring was moderate, 4% at Zurich, and 9% at Otsu. The field experiment at Otsu ended on July 1, 2016 because one third of the plants died due to temporary high air temperatures on July 2 and 3, 2016.
We recorded the initial plant size and the presence/absence of flowering stems to incorporate the effects of plant life-history traits on insect abundance. Initial plant size was evaluated by the length of the largest rosette leaf (mm) at the beginning of the field experiment, because this parameter represents plant size at the growth stage. The presence/absence of flowering stems was recorded 2 weeks after transplanting the plants. The final leaf size could not be evaluated because some large leaves received herbivore damage or started drying due to senescence at the end of the experiment; therefore, only the initial size was incorporated into our statistical analyses.
Response variables - Community indices were examined at three levels (i.e., component species, guilds and entire communities) as response variables in the following analyses. At the species level, we analyzed the number of individuals of each herbivorous species. We analyzed species for which more than 20 individuals were observed at each site, because statistical tests were difficult to apply to rare species. For the Zurich data, we analyzed the number of leaf holes as an indicator of damage by flea beetles. At the guild level, we classified herbivorous species into those feeding on external leaf tissues (i.e., leaf chewers) and those feeding on internal plant tissues (including sap suckers and leaf miners). We also separated herbivorous species into specialists on Brassicaceae (e.g., cabbage butterflies, cabbage sawflies and turnip flea beetles) and generalists on multiple plant families (some species of aphids and thrips) (Table 2; Additional file 2). The total number of insect individuals in each category were analyzed as guild level statistics. At the entire community level, we calculated species richness (i.e., number of species), Shannon’s diversity index H′ and the total number of insect individuals on individual plants. All of the response variables were ln(x + 1)-transformed to improve normality before statistical analyses. All statistical analyses were conducted using R version 3.2.0 . We utilized the rda function (in the vegan package: ) to perform the redundancy analysis. We used the lme function (in the nlme package: ) to estimate heritability, as described below.
Variation in insects on plant accessions - To quantify variation in insect communities among plant accessions and study sites, we performed a redundancy analysis to partition sources of variation in community composition into the plant accession, study sites and accession-by-site effects. The accession-by-site interaction was first analyzed 999-times by permutation tests and then the main effects of accessions and sites were examined without the interaction term. Then, we estimated broad-sense heritability H2 in a focal response, as the proportion of variance attributable to plant accessions. We used linear mixed models, in which the accession ID was assigned as a random effect. This variance component of random effect was estimated by the restricted maximum likelihood method [66, 67]. The significance of heritability was examined by likelihood ratio tests by comparing the linear models with or without the random effect of accession ID. This estimation of heritability was separately performed for the data from Zurich and Otsu. P-values were corrected by the false discovery rate (FDR) of multiple testing . Another option of estimating heritability is to incorporate a genetic distance matrix among natural accessions, as used in genome-wide association studies [e.g., 26, 27]. However, it was difficult to apply the same approach to single-gene mutants and the limited number of accessions; thus, we adopted the linear mixed models without the distance matrix to estimate broad-sense heritability.
Effects of plant traits - To address whether particular plant traits contributed to community members and composition, we used multiple regressions that considered the trichome density, PC1 and PC2 of aliphatic GSLs, the presence/absence of flowering stems, and initial plant size (mm) as explanatory variables. No explanatory variables were heavily correlated with each other (|r| < 0.6 for all pairs). We considered the difference of experimental block as a covariate. Because trichome density had a highly skewed distribution due to completely glabrous phenotypes, this variable was ln(x + 1)-transformed before the analysis. First, we tested the effects of plant traits on each response variable without the dataset on glabrous mutants. When detecting significant effects of trichomes on a particular herbivore, we then compared the two glabrous mutants and their parental accessions to test how the GL1 genes affect guild and community indices encompassing the focal herbivore. Multiple regression with linear mixed models was used to analyze trichome production, initial plant size (mm), and the presence/absence of flowering stems as explanatory variables. The difference in parental background (i.e., Ler-1 or Col-0) was considered as a random effect. We used the lme function with the maximum likelihood method for these mixed models. All of the continuous response and explanatory variables were standardized following a normal distribution, with zero mean and one variance, to make coefficients comparable between the linear models. P-values were corrected by FDR .
Abundance and communities of insects among plant accessions and study sites
We observed 15 insect species including flea beetles, sawflies, butterflies, moths, aphids and thrips on A. thaliana in the two field experiments (Table 2; Fig. 2). Of these insects, five and three species were specific to the Otsu and Zurich sites, respectively. Redundancy analysis and permutation tests confirmed that the plant accession, study site and accession-by-site effects exhibited significant sources of variation in the community composition (Accession, Sum of Squares (SS) = 0.99, F = 1.57, P < 0.001; Site, SS = 1.31, F = 36.2, P < 0.001; Accession-by-site, SS = 0.96, F = 1.54, P < 0.001 with 999 permutations; Fig. 3). Based on the sum of squares, the redundancy analysis indicated that 13% of variation in the insect community composition can be explained by the plant accessions and accession-by-site effects. We also found significant broad-sense heritability in species richness, Shannon diversity and total abundance of insects on A. thaliana, and its magnitude varied between the two study sites (10–11% and 15–30% heritability in Zurich and Otsu, respectively: Additional file 3).
Plant traits underlying the abundance and communities of insects
We examined whether the species, guild, and community structure of insects were affected by trichome density, glucosinolates, and life-history traits among natural accessions (Fig. 4). Significant effects of trichomes at the species and guild levels were observed (Fig. 4; Additional file 4). Two predominant leaf chewers, the striped flea beetle in Zurich and the turnip sawfly in Otsu, occurred less on hairy accessions than on accessions that produced low quantities of trichomes (Figs. 4 and 5b; Additional file 4). The number of leaf holes was lower on hairy plants compared to glabrous plants at the Zurich site (Figs. 4 and 5a; Additional file 4), indicating that trichomes have a resistance function against flea beetles. At the Otsu site, the abundance of the cabbage butterfly Pieris rapae was also lower on hairy plants (Fig. 4; Additional file 4). At the guild level, trichomes had significant negative effects on the leaf chewers at both sites (Figs. 4 and 5c, d; Additional file 4). In contrast to trichomes, aliphatic GSLs did not have any consistent effects on herbivore abundance. The first principal component of GSLs had negative effects on leaf chewers, specialist herbivores, species richness and total abundance at the Zurich site but no significant effects on these indices at the Otsu site (Fig. 4). The second principal component of GSLs was positively correlated with the abundance of turnip sawflies at the Otsu site and western flower thrips at the Zurich site (Fig. 4: Additional file 4). These effects of trichomes and GSLs were variable between the two sites with respect to insect richness, Shannon diversity and total abundance (Figs. 4 and 5e, f; Additional file 4). Initial plant size or the presence of flowering stems significantly increased insect richness, diversity and total abundance at both sites (Fig. 4; Additional file 4). The result that Kas-2 in Otsu and C24 in Zurich were less likely to be occupied by leaf chewers (Fig. 5) was due to the small plant size of these accessions (Additional file 5).
Comparing glabrous mutants and parental hairy accessions
We examined the effects of a single gene GLABRA1 on herbivory, guild and community indices encompassing two predominant leaf chewers, flea beetles and sawflies. At the species level, compared to parental accessions, the two glabrous mutants had significantly more leaf holes made by the flea beetles and the larvae of the turnip sawfly (Fig. 5a, b; Additional file 6). At the guild level, leaf chewers tended to occur more often on the two glabrous mutants than on each of the parental accessions (Fig. 5c, d), although this difference was not statistically significant in Zurich (Additional file 6). Among community indices, total abundance at the Otsu site was significantly lower on the glabrous mutants than on the hairy parents (coef. ± SE = − 0.46 ± 0.16, Z = − 2.86, P < 0.01; Additional file 6).
Insect species diversity on field-grown Arabidopsis thaliana
A number of researchers have studied leaf herbivory [17, 48, 60, 69], herbivorous fauna [19, 56, 57, 70], and plant fitness [17, 48, 60, 69] in Arabidopsis thaliana under field conditions; however, quantitative evidence remains limited in relation to insect community composition on this plant species. In the present study, we found that the community composition of 15 arthropod species significantly varied among A. thaliana accessions in the two distant field sites. Consistent with our results in the Zurich site, Harvey et al.  found that the summer cohort of A. thaliana were heavily attacked by Phyllotreta beetles and also harbored diamondback moths and aphids in a wild population of Europe. While winter-annual is considered the major life-cycle of A. thaliana [56, 57], the short generation time and seed bank may allow them to form alternative life-cycles in the wild. Near the Zurich site, in addition to the overwintering cohort, flowering and vegetative plants co-occurred during July and more importantly, the wild A. thaliana had many holes in their leaves (Fig. 1). Our evidence of a summer cohort suggests that the present field experiments using a variety of accessions should represent a part of the Arabidopsis-herbivore interactions observed in nature, and thereby highlights a further need to study insect communities across the multiple life stages of A. thaliana.
The plant apparency hypothesis assumed the importance of plant size and architectural traits in anti-herbivore defense , and this hypothesis has been supported by a meta-analysis , comparative study , and genome-wide association mapping . In the context of community genetics, plant architectures are a key predictor of insect community composition on perennial herbs  and woody plants . In the present study, we also found that the presence of flowering stems and larger plant size increased insect species richness and diversity on A. thaliana (Fig. 4; Additional file 4). Even though similar accessions produced flowering stems between the two study sites (Additional file 5), only at the Otsu site did the presence of flowering stems contribute to insect species diversity. This result could be partly because the flowering stems attracted three herbivore species, the cabbage butterfly, green peach aphid and western flower thrip, at the Otsu site (Fig. 4).
Effects of trichomes and the single gene GLABRA1 on insect abundance
In reverse genetic analysis, multiple independent mutants with a consistent phenotype are required to prove the function of a particular gene. In addition, using multiple genetic backgrounds of parent can also give a strong indication of gene function. In previous studies, the roles of single genes in modulating herbivory were quantified using mutants derived from a single parental accession [12, 13, 48]. Our common garden experiments illustrate the function of the GL1 gene against herbivory in natura using two distinct lines, Ler (gl1–1) and Col (gl1–2). Notably, the Br-0 and C24 accessions were the most susceptible accessions to leaf chewers at each site (Fig. 5c, d), and these two accessions also have disruptive mutations on GL1 [34, 43]. Such associations between GL1 polymorphism and anti-herbivore functions have also been reported in field populations of A. lyrata [9, 10] and A. halleri [11, 18]. Plant trichomes also prevent herbivory by sawflies  and flea beetles [75, 76] on Brassica cultivars. Together with these results, our present results indicate that plant trichomes and a single gene, GL1, play a key role in physical defense against crucifer feeders.
Laboratory experiments on single-gene mutants and natural accessions of A. thaliana suggested that plants with high trichome density resisted infestation by aphids [25, 77]. Under the two tested field conditions, trichomes had no significant effect on the abundance of aphids, possibly because aphids primarily occurred on flowering stems, on which the trichome density is low. In fact, the presence of flowering stems was positively correlated with the abundance of aphids (Additional file 4). These results support the limited associations between aphid abundance and GL1 polymorphism detected in field-grown A. halleri . In addition, we could not detect any significant effects of trichomes and GL1 on the abundance of larval P. xylostella, even though trichomes prevent adult moths ovipositing on A. thaliana under laboratory conditions . Handley et al.  focused on several northern accessions of A. thaliana, whereas the current experimental setting covers a wider geographical range of natural accessions. A recent genome-wide association study using 350 natural accessions also found no significant association between GL1 polymorphism and herbivory by P. xylostella . Combined with the previous studies, our present results from field-grown Arabidopsis exemplify the importance of testing diverse accessions and environmental conditions.
Varying effects of glucosinolates on specialist herbivores
Some components of GSLs had negative effects whereas others had no or positive effects on the specialist herbivores. A possible explanation for these varying effects is that GSL data quantified under laboratory conditions might be insufficient to reflect the GSL accumulation in the field, due to phenotypic plasticity and the induced response of GSLs to herbivory [23, 24, 48, 78]. For example, the cabbage butterfly P. rapae can modify the expression level of the MAM1 gene, which involves a chain elongation of aliphatic GSL [23, 24]. It is also possible that some specialist herbivores overcome GSLs [20, 21, 61, 79, 80] and have thus obscured the effects of GSLs on the community composition. Specifically, the striped flea beetle P. striolata efficiently sequesters 4-methylthiobutyl from A. thaliana, a short-chain aliphatic GSL . The larvae of A. rosae sawflies also sequester GSLs , whereas adults utilize isothiocyanates, which are breakdown products of aliphatic GSLs, to find host plants . Unlike the case of trichomes, an on-site quantification of GSLs is needed to resolve their effects on herbivore abundance in the field.
Our field investigation showed a genetic basis in the insect community assemblage on A. thaliana, and the advantages of the functional allele of GL1 in avoiding leaf chewers. In Brassicaceae plants, evidence is accumulating to suggest that genetic variation within a plant species alters insect community composition and, in turn, exerts selection on plant defense [4, 81, 82]. Variation in the trichome density is also observed across Brassicaceae plants [74,75,76], where GL1 orthologs affect the trichome density . In the context of community genetics, the present study on GL1 provides evidence of a key gene affecting the community composition of crucifer-feeding insects. Future studies should assess the relative importance of single genes and quantitative genetic variation towards a comprehensive understanding of plant genetic effects on insect community assembly.
False discovery rate
- GL1 :
Principal component analysis
Schoonhoven LM, van Loon JJA, Dicke M. Insect-plant biology. 2nd ed. Oxford, UK: Oxford University Press; 2005.
Barbour MA, Rodriguez-Cabal MA, Wu ET, et al. Multiple plant traits shape the genetic basis of herbivore community assembly. Funct Ecol. 2015;29:995–1006.
Bernhardsson C, Robinson KM, Abreu IN, et al. Geographic structure in metabolome and herbivore community co-occurs with genetic structure in plant defence genes. Ecol Lett. 2013;16:791–8.
Poelman EH, van Dam NM, van Loon JJA, Vet LE, Dicke M. Chemical diversity in Brassica oleracea affects biodiversity of insect herbivores. Ecology. 2009;90:1863–77.
Whitham TG, Bailey JK, Schweitzer JA, et al. A framework for community and ecosystem genetics: from genes to ecosystems. Nat Rev Genet. 2006;7:510–23.
Johnson MTJ, Agrawal AA. Plant genotype and environment interact to shape a diverse arthropod community on evening primrose (Oenothera biennis). Ecology. 2005;86:874–85.
Keith AR, Bailey JK, Whitham TG. A genetic basis to community repeatability and stability. Ecology. 2010;91:3398–406.
Maddox GD, Root RB. Resistance to 16 diverse species of herbivorous insects within a population of goldenrod, Solidago altissima: genetic variation and heritability. Oecologia. 1987;72:8–14.
Kivimäki M, Kärkkäinen K, Gaudeul M, et al. Gene, phenotype and function: GLABROUS1 and resistance to herbivory in natural populations of Arabidopsis lyrata. Mol Ecol. 2007;16:453–62.
Løe G, Toräng P, Gaudeul M, Ågren J. Trichome production and spatiotemporal variation in herbivory in the perennial herb Arabidopsis lyrata. Oikos. 2007;116:134–42.
Sato Y, Kudoh H. Tests of associational defence provided by hairy plants for glabrous plants of Arabidopsis halleri subsp. gemmifera against insect herbivores. Ecol Entomol. 2015;40:269–79.
Kessler A, Halitschke R, Baldwin IT. Silencing the jasmonate cascade: induced plant defenses and insect populations. Science. 2004;305:665–8.
Kallenbach M, Bonaventure G, Gilardoni PA, et al. Empoasca leafhoppers attack wild tobacco plants in a jasmonate-dependent manner and identify jasmonate mutants in natural populations. Proc Natl Acad Sci U S A. 2012;109:E1548–57.
van Poecke RMP. Arabidopsis-insect interactions. Arabidopsis Book. 2007;5:e0107.
The 1001 Genomes Consortium. 1,135 genomes reveal the global pattern of polymorphism in Arabidopsis thaliana. Cell. 2016;166:481–91.
Handley R, Ekbom B, Ågren J. Variation in trichome density and resistance against a specialist insect herbivore in natural populations of Arabidopsis thaliana. Ecol Entomol. 2005;30:284–92.
Mauricio R. Costs of resistance to natural enemies in field populations of the annual plant Arabidopsis thaliana. Am Nat. 1998;151:20–8.
Sato Y, Kudoh H. Herbivore-mediated interaction promotes the maintenance of trichome dimorphism through negative frequency-dependent selection. Am Nat. 2017;190:E67–77.
Falk KL, Kästner J, Bodenhausen N, et al. The role of glucosinolates and the jasmonic acid pathway in resistance of Arabidopsis thaliana against molluscan herbivores. Mol Ecol. 2014;23:1188–203.
Kliebenstein D, Pedersen D, Barker B, Mitchell-Olds T. Comparative analysis of quantitative trait loci controlling glucosinolates, myrosinase and insect resistance in Arabidopsis thaliana. Genetics. 2002;161:325–32.
Müller C. Interactions between glucosinolate- and myrosinase-containing plants and the sawfly Athalia rosae. Phytochem Rev. 2009;8:121–34.
Reymond P, Bodenhausen N, van Poecke RMP, et al. A conserved transcript pattern in response to a specialist and a generalist herbivore. Plant Cell. 2004;16:3132–47.
Mewis I, Tokuhisa JG, Schultz JC, Appel HM, Ulrichs C, Gershenzon J. Gene expression and glucosinolate accumulation in Arabidopsis thaliana in response to generalist and specialist herbivores of different feeding guilds and the role of defense signaling pathways. Phytochemistry. 2006;67:2450–62.
Kempema LA, Cui X, Holzer FM, Walling LL. Arabidopsis transcriptome changes in response to phloem-feeding silverleaf whitefly nymphs: similarities and distinctions in responses to aphids. Plant Physiol. 2007;143:849–65.
Züst T, Joseph B, Shimizu KK, et al. Using knockout mutants to reveal the growth costs of defensive traits. Proc R Soc Lond B. 2011;278:2598–603.
Davila Olivas NH, Frago E, Thoen MPM, et al. Natural variation in life history strategy of Arabidopsis thaliana determines stress responses to drought and insects of different feeding guilds. Mol Ecol. 2017;26:2959–77.
Davila Olivas NH, Kruijer W, Gort G, et al. Genome-wide association analysis reveals distinct genetic architectures for single and combined stress responses in Arabidopsis thaliana. New Phytol. 2017;213:838–51.
Shimizu KK, Kudoh H, Kobayashi MJ. Plant sexual reproduction during climate change: gene function in natura studied by ecological and evolutionary systems biology. Ann Bot. 2011;108:777–87.
Nagano AJ, Sato Y, Mihara M, et al. Deciphering and prediction of transcriptome dynamics under fluctuating field conditions. Cell. 2012;151:1358–69.
Kudoh H. Molecular phenology in plants: in natura systems biology for the comprehensive understanding of seasonal responses under natural environments. New Phytol. 2016;210:399–412.
Yamasaki E, Altermatt F, Cavender-Bares J, et al. Genomics meets remote sensing in global change studies: monitoring and predicting phenology, evolution and biodiversity. Curr Opin Environ Sustain. 2017;29:177–86.
Zaidem M, Groen SC, Purugganan MD. Evolutionary and ecological functional genomics, from lab to the wild. Plant J. 2018; https://doi.org/10.1111/tpj.14167.
Symonds VV, Godoy AV, Alconada T, et al. Mapping quantitative trait loci in multiple populations of Arabidopsis thaliana identifies natural allelic variation for trichome density. Genetics. 2005;169:1649–58.
Bloomer RH, Juenger TE, Symonds VV. Natural variation in GL1 and its effects on trichome density in Arabidopsis thaliana. Mol Ecol. 2012;21:3501–15.
Atwell S, Huang YS, Vilhjálmsson BJ, et al. Genome-wide association study of 107 phenotypes in Arabidopsis thaliana inbred lines. Nature. 2010;465:627–31.
Roux F, Bergelson J. The genetics underlying natural variation in the biotic interactions of Arabidopsis thaliana: the challenges of linking evolutionary genetics and community ecology. Curr Top Dev Biol. 2016;119:111–56.
Chan EKF, Rowe HC, Kliebenstein DJ. Understanding the evolution of defense metabolites in Arabidopsis thaliana using genome-wide association mapping. Genetics. 2010;185:991–1007.
Oppenheimer DG, Herman PL, Sivakumaran S, Esch J, Marks MD. A myb gene required for leaf trichome differentiation in Arabidopsis is expressed in stipules. Cell. 1991;67:483–93.
Yoshida Y, Sano R, Wada T, et al. Jasmonic acid control of GLABRA3 links inducible defense and trichome patterning in Arabidopsis. Development. 2009;136:1039–48.
Hülskamp M. Plant trichomes: a model for cell differentiation. Nat Rev Mol Cell Biol. 2004;5:471–80.
Ishida T, Kurata T, Okada K, Wada T. A genetic regulatory network in the development of trichomes and root hairs. Annu Rev Plant Biol. 2008;59:365–86.
Grebe M. The patterning of epidermal hairs in Arabidopsis - updated. Cur Opin Plant Biol. 2012;15:31–7.
Hauser MT, Harr B, Schlötterer C. Trichome distribution in Arabidopsis thaliana and its close relative Arabidopsis lyrata: molecular analysis of the candidate gene GLABROUS1. Mol Biol Evol. 2001;18:1754–63.
Kawagoe T, Shimizu KK, Kakutani T, Kudoh H. Coexistence of trichome variation in a natural plant population: a combined study using ecological and candidate gene approaches. PLoS One. 2011;6:e22184.
Larkin JC, Walker JD, Bolognesi-Winfield AC, et al. Allele-specific interactions between ttg and gl1 during trichome development in Arabidopsis thaliana. Genetics. 1999;151:1591–604.
Shimizu KK. Ecology meets molecular genetics in Arabidopsis. Popul Ecol. 2002;44:221–33.
Wilczek AM, Roe JL, Knapp MC, et al. Effects of genetic perturbation on seasonal life history plasticity. Science. 2009;323:930–4.
Kerwin R, Feusier J, Corwin J, et al. Natural genetic variation in Arabidopsis thaliana defense metabolism genes modulates field fitness. eLife. 2015;4:e05604.
Heidel AJ, Clarke JD, Antonovics J, Dong X. Fitness costs of mutations affecting the systemic acquired resistance pathway in Arabidopsis thaliana. Genetics. 2004;168:2197–206.
Shimizu KK, Tsuchimatsu T. Evolution of selfing: recurrent patterns in molecular adaptation. Ann Rev Ecol Evol Syst. 2015;46:593–622.
Chiang GC, Barua D, Dittmar E, et al. Pleiotropy in the wild: the dormancy gene DOG1 exerts cascading control on life cycles. Evolution. 2013;67:883–93.
Bentsink L, Hanson J, Hanhart CJ, et al. Natural variation for seed dormancy in Arabidopsis is regulated by additive genetic and molecular pathways. Proc Natl Acad Sci U S A. 2010;107:4264–9.
Thompson L. The spatiotemporal effects of nitrogen and litter on the population dynamics of Arabidopsis thaliana. J Ecol. 1994;82:63–88.
Taylor MA, Cooper MD, Sellamuthu R, et al. Interacting effects of genetic variation for seed dormancy and flowering time on phenology, life history, and fitness of experimental Arabidopsis thaliana populations over multiple generations in the field. New Phytol. 2017;216:291–302.
Shindo C, Aranzana MJ, Lister C, et al. Role of FRIGIDA and FLOWERING LOCUS C in determining variation in flowering time of Arabidopsis. Plant Physiol. 2005;138:1163–73.
Arany AM, de Jong TJ, van der Meijden E. Herbivory and abiotic factors affect population dynamics of Arabidopsis thaliana in a sand dune area. Plant Biol. 2005;7:549–56.
Harvey JA, Witjes LMA, Benkirane M, Duyts H, Wagenaar R. Nutritional suitability and ecological relevance of Arabidopsis thaliana and Brassica oleracea as food plants for the cabbage butterfly, Pieris rapae. Plant Ecol. 2007;189:117–26.
Wentzell AM, Rowe HC, Hansen BG, Ticconi C, Halkier BA, Kliebenstein DJ. Linking metabolic QTLs with network and cis-eQTLs controlling biosynthetic pathways. PLoS Genet. 2007;3:e162.
Joseph B, Corwin JA, Zust T, et al. Hierarchical nuclear and cytoplasmic genetic architectures for plant growth and defense within Arabidopsis. Plant Cell. 2013;25:1929–45.
Brachi B, Meyer CG, Villoutreix R, et al. Coselected genes determine adaptive variation in herbivore resistance throughout the native range of Arabidopsis thaliana. Proc Natl Acad Sci U S A. 2015;112:4032–7.
Okamura Y, Sawada Y, Hirai MY, Murakami M. Effects of different secondary metabolite profiles in plant defense syndromes on specialist and generalist herbivores. Entomol Sci. 2016;19:97–103.
Ruppel RF. Cumulative insect-days as an index of crop protection. J Econ Entomol. 1983;76:375–7.
R Core Team. R: A language and environment for statistical computing. Vienna: R Foundation for Statistical Computing; 2015. http://www.R-project.org/;.
Oksanen J, Blanchet FG, Kindt R, et al. vegan: community ecology package. R package version 2.3–4. http://CRAN.R-project.org/package=vegan; 2015.
Pinheiro J, Bates D, DebRoy S, et al. nlme: linear and nonlinear mixed effects models. R package version 3.1–126. URL: http://CRAN.R-project.org/package=nlme; 2016.
Patterson HD, Thompson R. Maximum likelihood estimation of variances. Biometrika. 1971;58:545–54.
Falconer DS. Introduction to quantitative genetics. In: Pearson Education India. 3rd ed; 1989.
Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J Royal Stat Soc B. 1995;57:289–300.
Mauricio R, Rausher MD. Experimental manipulation of putative selective agents provides evidence for the role of natural enemies in the evolution of plant defense. Evolution. 1997;51:1435–44.
Bidart-Bouzat MG, Kliebenstein DJ. Differential levels of insect herbivory in the field associated with genotypic variation in glucosinolates in Arabidopsis thaliana. J Chem Ecol. 2008;34:1026–37.
Feeny P. Plant apparency and chemical defense. In: Wallace JW, Mansell RL, editors. Biochemical interaction between plants and insects. Boston, MA: Springer; 1976. p. 1–40.
Carmona D, Lajeunesse MJ, Johnson MTJ. Plant traits that predict resistance to herbivores. Funct Ecol. 2011;25:358–67.
Schlinkert H, Westphal C, Clough Y, et al. Feeding damage to plants increases with plant size across 21 Brassicaceae species. Oecologia. 2015;179:455–66.
Travers-Martin N, Müller C. Matching plant defence syndromes with performance and preference of a specialist herbivore. Funct Ecol. 2008;22:1033–43.
Agerbirk N, Ørgaard M, Nielsen JK. Glucosinolates, flea beetle resistance, and leaf pubescence as taxonomic characters in the genus Barbarea (Brassicaceae). Phytochemistry. 2003;63:69–80.
Soroka JJ, Holowachuk JM, Gruber MY, Grenkow LF. Feeding by flea beetles (Coleoptera: Chrysomelidae; Phyllotreta spp.) is decreased on canola (Brassica napus) seedlings with increased trichome density. J Econ Entomol. 2011;104:125–36.
Züst T, Heichinger C, Grossniklaus U, et al. Natural enemies drive geographic variation in plant defenses. Science. 2012;338:116–9.
Shirakawa M, Hara-Nishimura I. Specialized vacuoles of myrosin cells: chemical defense strategy in Brassicales plants. Plant Cell Physiol. 2018;59:1309–16.
Barker AM, Molotsane R, Müller C, et al. Chemosensory and behavioural responses of the turnip sawfly, Athalia rosae, to glucosinolates and isothiocyanates. Chemoecology. 2006;16:209–18.
Beran F, Pauchet Y, Kunert G, et al. Phyllotreta striolata flea beetles use host plant defense compounds to create their own glucosinolate-myrosinase system. Proc Natl Acad Sci U S A. 2014;111:7349–54.
Lankau R, Strauss S. Community complexity drives patterns of natural selection on a chemical defense of Brassica nigra. Am Nat. 2008;171:150–61.
Lankau RA, Kliebenstein DJ. Competition, herbivory and genetics interact to determine the accumulation and fitness consequences of a defence metabolite. J Ecol. 2009;97:78–88.
Li F, Kitashiba H, Nishio T. Association of sequence variation in brassica GLABRA1 orthologs with leaf hairiness. Mol Breeding. 2011;28:577–84.
The authors thank Jordi Bascompte and Matthew Barbour for their valuable discussions, Satoshi Oda and Aki Morishima for helping with plant cultivation and fieldwork, and three anonymous reviewers for their comments on the manuscript.
This study was supported by JSPS Postdoctoral Fellowship (Grant Number, 16 J30005) and JST PRESTO (JPMJPR17Q4) to Y.S., JST CREST (JPMJCR15O2) and JSPS KAKENHI (JP16H06171 and JP16H01473) to A.J.N., Swiss National Foundation, URPP Global Change and Biodiversity of the University of Zurich, JST CREST (JPMJCR16O3), Japan, MEXT KAKENHI (18H04785) to K.K.S. The field experiment in Japan was supported by the Joint Usage/Research Grant of Center for Ecological Research, Kyoto University. The funders have no role in study design, data analysis, or decision to publish.
Availability of data and materials
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
The first and second principal component (PC1 and PC2) summarizing the total amount (nmol/mg flesh weight) of C3-, C4-, C5-, C7- and C8-Aliphatic glucosinolates for 17 accessions of A. thaliana (compiled from Chan et al. ). Red arrows indicate contributions of each glucosinolate to PC1 and PC2. (PNG 145 kb)
Cumulative number of each insect species in Zurich, Switzerland (left chart) and Otsu, Japan (right chart) throughout the experiments. See Table 2 for the name of the arthropod species. Notes: 1Total number of parasitoid wasps and mummified aphids; 2This species is a non-insect arthropod; 3Only a dwelling trace was observed. (PNG 319 kb)
Likelihood ratio tests for estimating broad-sense heritability H2. Likelihood ratio, LR-χ2, was tested by comparing the models with and without a random effect of the accession ID. P-values were based on a χ2 distribution with one degree of freedom and corrected by false discovery rate (FDR) . Bold values indicate significant H2 at Pfdr < 0.05. Bars indicate no information available due to low abundance. (XLSX 11 kb)
Effects of trichome density, the first and second principal component (PC1 and PC2) of aliphatic glucosinolates (GSLs), presence of flowering stems and initial plant sizes on insect abundance and community composition for 17 natural accessions of Arabidopsis thaliana. Standardized coefficient (Coef.), standard error (SE), t-value, and FDR-corrected P-values are shown for each explanatory variable from Zurich, Switzerland and Otsu, Japan. Bold values highlight significant effects at Pfdr < 0.05. (XLSX 14 kb)
Insect abundance data. (XLS 137 kb)
Effects of trichome density, presence of flowering stem and initial plant size on insect abundance and community composition in a comparison between glabrous mutants and their parental accessions. Standardized coefficient (Coef.), standard error (SE), t-value, and FDR-corrected P-values are shown for each explanatory variable. Bold values highlight significant effects at Pfdr < 0.05. NA means no data available. The trichome density represents differences between the null and hypomorphic mutants; as the trichome density had a significant effect on flea beetles and sawflies (Additional file 4), we focused on species, guild and community indices comprising these two leaf chewers. (XLSX 11 kb)
R source code for statistical analyses (TXT 10 kb)