Fine-scale frequency differentiation along a herbivory gradient in the trichome dimorphism of a wild Arabidopsis

Altered stomatal patterning accompanies a trichome dimorphism in a natural population of Arabidopsis

Trichomes are large epidermal cells on the surface of leaves that are thought to deter herbivores, yet the presence of trichomes can also negatively impact plant growth and reproduction. Stomatal guard cells and trichomes have shared developmental origins, and experimental manipulation of trichome formation can lead to changes in stomatal density. The influence of trichome formation upon stomatal development in natural populations of plants is currently unknown. Here, we show that a natural population of Arabidopsis halleri that includes hairy (trichome-bearing) and glabrous (no trichomes) morphs has differences in stomatal density that are associated with this trichome dimorphism. We found that glabrous morphs had significantly greater stomatal density and stomatal index than hairy morphs. One interpretation is that this arises from a trade-off between the proportions of cells that have trichome and guard cell fates during leaf development. The differences in stomatal density between the two morphs might have impacts upon environmental adaptation, in addition to herbivory deterrence caused by trichome development.

Transcriptional Variation in Glucosinolate Biosynthetic Genes and Inducible Responses to Aphid Herbivory on Field-Grown Arabidopsis thaliana

Recently, increasing attempts have been made to understand how plant genes function in natura. In this context, transcriptional profiles represent plant physiological status in response to environmental stimuli. Herein, we combined high-throughput RNA-Seq with insect survey data on 19 accessions of Arabidopsis thaliana grown at a field site in Switzerland. We found that genes with the gene ontology (GO) annotations of "glucosinolate biosynthetic process" and "response to insects" were most significantly enriched, and the expression of these genes was highly variable among plant accessions. Nearly half of the total expression variation in the glucosinolate biosynthetic genes (AOPs, ESM1, ESP, and TGG1) was explained by among-accession variation. Of these genes, the expression level of AOP3 differed among Col-0 accession individuals depending on the abundance of the mustard aphid (Lipaphis erysimi). We also found that the expression of the major cis-jasmone activated gene CYP81D11 was positively correlated with the number of flea beetles (Phyllotreta striolata and Phyllotreta atra). Combined with the field RNA-Seq data, bioassays confirmed that AOP3 was up-regulated in response to attack by mustard aphids. The combined results from RNA-Seq and our ecological survey illustrate the feasibility of using field transcriptomics to detect an inducible defense, providing a first step towards an in natura understanding of biotic interactions involving phenotypic plasticity.

Evolutionary potential varies across populations and traits in the neotropical oak Quercus oleoides

Heritable variation in polygenic (quantitative) traits is critical for adaptive evolution and is especially important in this era of rapid climate change. In this study, we examined the levels of quantitative genetic variation of populations of the tropical tree Quercus oleoides Cham. and Schlect. for a suite of traits related to resource use and drought resistance. We tested whether quantitative genetic variation differed across traits, populations and watering treatments. We also tested potential evolutionary factors that might have shaped such a pattern: selection by climate and genetic drift. We measured 15 functional traits on 1322 1-year-old seedlings of 84 maternal half-sib families originating from five populations growing under two watering treatments in a greenhouse. We estimated the additive genetic variance, coefficient of additive genetic variation and narrow-sense heritability for each combination of traits, populations and treatments. In addition, we genotyped a total of 119 individuals (with at least 20 individuals per population) using nuclear microsatellites to estimate genetic diversity and population genetic structure. Our results showed that gas exchange traits and growth exhibited strikingly high quantitative genetic variation compared with traits related to leaf morphology, anatomy and photochemistry. Quantitative genetic variation differed between populations even at geographical scales as small as a few kilometers. Climate was associated with quantitative genetic variation, but only weakly. Genetic structure and diversity in neutral markers did not relate to coefficient of additive genetic variation. Our study demonstrates that quantitative genetic variation is not homogeneous across traits and populations of Q. oleoides. More importantly, our findings suggest that predictions about potential responses of species to climate change need to consider population-specific evolutionary characteristics.

Heritable plant phenotypes track light and herbivory levels at fine spatial scales

Organismal phenotypes often co-vary with environmental variables across broad geographic ranges. Less is known about the extent to which phenotypes match local conditions when multiple biotic and abiotic stressors vary at fine spatial scales. Bittercress (Brassicaceae: Cardamine cordifolia), a perennial forb, grows across a microgeographic mosaic of two contrasting herbivory regimes: high herbivory in meadows (sun habitats) and low herbivory in deeply shaded forest understories (shade habitats). We tested for local phenotypic differentiation in plant size, leaf morphology, and anti-herbivore defense (realized resistance and defensive chemicals, i.e., glucosinolates) across this habitat mosaic through reciprocal transplant-common garden experiments with clonally propagated rhizomes. We found habitat-specific divergence in morphological and defensive phenotypes that manifested as contrasting responses to growth in shade common gardens: weak petiole elongation and attenuated defenses in populations from shade habitats, and strong petiole elongation and elevated defenses in populations from sun habitats. These divergent phenotypes are generally consistent with reciprocal local adaptation: plants from shade habitats that naturally experience low herbivory show reduced investment in defense and an attenuated shade avoidance response, owing to its ineffectiveness within forest understories. By contrast, plants from sun habitats with high herbivory show shade-induced elongation, but no evidence of attenuated defenses canonically associated with elongation in shade-intolerant plant species. Finally, we observed differences in flowering phenology between habitat types that could potentially contribute to inter-habitat divergence by reducing gene flow. This study illuminates how clonally heritable plant phenotypes track a fine-grained mosaic of herbivore pressure and light availability in a native plant.

The Long-Term "In Natura" Study Sites of Arabidopsis halleri for Plant Transcription and Epigenetic Modification Analyses in Natural Environments

The majority of organismal phenomena show functional significance in the context of natural environments. However, we know little about how dynamic gene expression is controlled under natural complex conditions. One of the most attractive challenges in current biology is to understand organismal functions in natural environments. We established and have developed long-term "in natura" study sites of Arabidopsis halleri to evaluate precise control of gene expression in natural environments. At the sites, we monitored meteorological factors, recorded plant growth and phenology, and collected RNA and chromatin samples to investigate dynamics of transcription and epigenetic modifications. Here, we introduce the in natura study sites, especially with the emphasis on methodologies for setting up study sites in natural plant populations and collecting samples used in transcriptomics and epigenetics in natural environments. Although the methods introduced here need to be modified depending on situations of one's study systems, our case can be a model for planning new in natura studies.

Fine-scale frequency differentiation along a herbivory gradient in the trichome dimorphism of a wild Arabidopsis

Geographic variation is commonly observed in plant resistance traits, where plant species might experience different selection pressure across a heterogeneous landscape. Arabidopsis halleri subsp. gemmifera is dimorphic for trichome production, generating two morphs, trichome-producing (hairy) and trichomeless (glabrous) plants. Trichomes of A. halleri are known to confer resistance against the white butterfly, cabbage sawfly, and brassica leaf beetle, but not against flea beetles. We combined leaf damage, microclimate, and microsatellite loci data of 26 A. halleri populations in central Japan, to explore factors responsible for fine-scale geographic variation in the morph frequency. We found that hairy plants were less damaged than glabrous plants within populations, but the among-site variation was the most significant source of variation in the individual-level damage. Fixation index (Gst″) of a putative trichome locus exhibited a significant divergence along population-level damage with an exception of an outlier population, inferring the local adaptation to herbivory. Notably, this outlier was a population wherein our previous study reported a balancing role of the brassica leaf beetle Phaedon brassicae on the morph frequency. This differentiation of the trichome locus was unrelated to neutral genetic differentiation (evaluated by Gst″ of microsatellite loci) and meteorological factors (including temperature and solar radiation). The present findings, combined with those of our previous work, provide suggestive evidence that herbivore-driven divergence and occasional outbreak of a specific herbivore have jointly contributed to the ecogeographic pattern in the frequency of two morphs.