High outcrossing rates in a self-compatible and highly aggregated host-generalist mistletoe.

Plants have evolved various strategies to avoid inbreeding, but the mass flowering displayed by many plants predisposes them to within-plant pollen movements and self-pollination. Mistletoes often aggregate at multiple spatial scales. Their bird pollinators often visit several flowers of the same individual and of others on the same host tree. We hypothesized that hermaphroditic mistletoes have self-incompatibility mechanisms that reduce or prevent selfing. Whether their spatial distribution, affected by host specificity, host distribution, and the behaviour of seed dispersers, influences their mating system and population genetic structure remains unclear. We studied how mating system and spatial distribution affect genetic structure in four populations of the host-generalist mistletoe Dendrophthoe pentandra in southwestern China using microsatellite markers and progeny arrays. We also characterized the fine-scale spatial genetic structure among 166 mistletoes from four host trees in one population. Prevalence and intensity of infection both appeared to vary among host species, strongly affecting the degree of aggregation. Host tree size had a strong effect on infection intensity. Surprisingly, manual pollination experiments indicated that D. pentandra is self-compatible, but genetic analyses revealed that outcrossing rates were higher than expected in all four populations (MLTR tm 0.83-1.20, Bayesian tm 0.772-0.952). Spatial genetic structure was associated with distance between host trees but not at shorter scales (within hosts). Our results demonstrate that the combination of bird pollination, bird-mediated seed dispersal, and post-dispersal processes result in outcrossing and maintain relatively high diversity in the presence of biparental inbreeding, despite very high local densities and possible self-compatibility.

Quantifying heritability and estimating evolutionary potential in the wild when individuals that share genes also share environments.

Accurate heritability estimates for fitness-related traits are required to predict an organism's ability to respond to global change. Heritability estimates are theoretically expected to be inflated if, due to limited dispersal, individuals that share genes are also likely to share similar environments. However, if relatives occupy similar environments due, at least partly, to genetic variation for habitat selection, then accounting for environmental similarity in quantitative genetic models may result in diminished heritability estimates in wild populations. This potential issue has been pointed out in the literature, but has not been evaluated by empirical studies. Here, we investigate whether environmental similarity among individuals can be partly explained by genetic variation for habitat selection, and how this link potentially blurs estimates for heritability in fitness-related traits. Using intensive GPS monitoring, we quantified home-range habitat composition for 293 roe deer inhabiting a heterogeneous landscape to assess environmental similarity. To investigate if environmental similarity might harbour genetic variation, we combined genome-wide data in a quantitative genetic framework to evaluate genetic variation for home-range habitat composition, which is partly the result of habitat selection at settlement. Finally, we explored how environmental similarity affects heritability estimates for behaviours related to the risk avoidance-resource acquisition trade-off (i.e. being in open habitat and distance to roads) and proxies of individual performance (i.e. body mass and hind foot length). We found substantial heritability for home-range habitat composition, with estimates ranging from 0.40 (proportion of meadows) to 0.85 (proportion of refuge habitat). Accounting for similarity in habitat composition between relatives decreased the heritability estimates for both behavioural and morphological traits (reduction ranging from 55% to 100% and from 22% to 41% respectively). As a consequence, only half of these heritability estimates remained significantly different from zero. Our results show that similar genotypes occupy similar environments, which could lead to heritable variation being incorrectly attributed to environmental effects. To accurately distinguish the sources of phenotypic variation and predict the ability of organisms to respond to global change, it is necessary to develop quantitative genetic studies investigating the mechanisms underpinning environmental similarity among relatives.

Non-reproducible signals of adaptation to elevation between open and understorey microhabitats in snapdragon plants.

Experimental studies on local adaptation rarely investigate how different environmental variables might modify signals of adaptation or maladaptation. In plant common garden experiments, signals of adaptation or maladaptation to elevation are usually investigated in open habitats under full light. However, most plants inhabit heterogeneous habitats where environmental conditions differ. Understorey microhabitats are common and differ in terms of tree shade, temperature, water availability, microbiota, allelochemicals etc. Germination is a fitness-related trait of major importance for the adaptation of plants to contrasted climate conditions. It is affected by shade in snapdragon plants (Antirrhinum majus) and many other plant species. Here, we tested for the reproducibility of signals extrapolated from germination results between open and understorey microhabitats in two parapatric snapdragon plant subspecies (A. m. striatum and A. m. pseudomajus) characterized by a similar elevation range by using common garden experiments at different elevations. Signals observed under one microhabitat systematically differed in the other. Most scenarios could be inferred, with signals either shifting, appearing or disappearing between different environments. Our findings imply that caution should be taken when extrapolating the evolutionary significance of these types of experimental signals because they are not stable from one local environmental condition to the next. Forecasting the ability of plants to adapt to environmental changes based on common garden and reciprocal transplant experiments must account for the multivariate nature of the environment.

Natural selection fluctuates at an extremely fine spatial scale inside a wild population of snapdragon plants.

Spatial variation in natural selection is expected to shape phenotypic variation of wild populations and drive their evolution. Although evidence of phenotypic divergence across populations experiencing different selection regimes is abundant, investigations of intrapopulation variation in selection pressures remain rare. Fine-grained spatial environmental heterogeneity can be expected to influence selective forces within a wild population and thereby alter its fitness function by producing multiple fitness optima at a fine spatial scale. Here, we tested this hypothesis in a wild population of snapdragon plants living on an extremely small island in southern France (about 7500 m2 ). We estimated the spline-based fitness function linking individuals' fitness and five morphological traits in interaction with three spatially variable ecological drivers. We found that selection acting on several traits varied both in magnitude and direction in response to environmental variables at the scale of a meter. Our findings illustrate how different phenotypes can be selected at different locations within a population in response to environmental variation. Investigating spatial variation in selection within a population, in association with ecological conditions, represents an opportunity to identify putative ecological drivers of selection in the wild.

Thyroid hormones regulate the formation and environmental plasticity of white bars in clownfishes.

Determining how plasticity of developmental traits responds to environmental conditions is a challenge that must combine evolutionary sciences, ecology, and developmental biology. During metamorphosis, fish alter their morphology and color pattern according to environmental cues. We observed that juvenile clownfish (Amphiprion percula) modulate the developmental timing of their adult white bar formation during metamorphosis depending on the sea anemone species in which they are recruited. We observed an earlier formation of white bars when clownfish developed with Stichodactyla gigantea (Sg) than with Heteractis magnifica (Hm). As these bars, composed of iridophores, form during metamorphosis, we hypothesized that timing of their development may be thyroid hormone (TH) dependent. We treated clownfish larvae with TH and found that white bars developed earlier than in control fish. We further observed higher TH levels, associated with rapid white bar formation, in juveniles recruited in Sg than in Hm, explaining the faster white bar formation. Transcriptomic analysis of Sg recruits revealed higher expression of duox, a dual oxidase implicated in TH production as compared to Hm recruits. Finally, we showed that duox is an essential regulator of iridophore pattern timing in zebrafish. Taken together, our results suggest that TH controls the timing of adult color pattern formation and that shifts in duox expression and TH levels are associated with ecological differences resulting in divergent ontogenetic trajectories in color pattern development.

Phenotypic Response to Light Versus Shade Associated with DNA Methylation Changes in Snapdragon Plants (Antirrhinum majus).

The phenotypic plasticity of plants in response to change in their light environment, and in particularly, to shade is a schoolbook example of ecologically relevant phenotypic plasticity with evolutionary adaptive implications. Epigenetic variation is known to potentially underlie plant phenotypic plasticity. Yet, little is known about its role in ecologically and evolutionary relevant mechanisms shaping the diversity of plant populations in nature. Here we used a reference-free reduced representation bisulfite sequencing method for non-model organisms (epiGBS) to investigate changes in DNA methylation patterns across the genome in snapdragon plants (Antirrhinum majus L.). We exposed plants to sunlight versus artificially induced shade in four highly inbred lines to exclude genetic confounding effects. Our results showed that phenotypic plasticity in response to light versus shade shaped vegetative traits. They also showed that DNA methylation patterns were modified under light versus shade, with a trend towards global effects over the genome but with large effects found on a restricted portion. We also detected the existence of a correlation between phenotypic and epigenetic variation that neither supported nor rejected its potential role in plasticity. While our findings imply epigenetic changes in response to light versus shade environments in snapdragon plants, whether these changes are directly involved in the phenotypic plastic response of plants remains to be investigated. Our approach contributed to this new finding but illustrates the limits in terms of sample size and statistical power of population epigenetic approaches in non-model organisms. Pushing this boundary will be necessary before the relationship between environmentally induced epigenetic changes and phenotypic plasticity is clarified for ecologically relevant mechanisms with evolutionary implications.

Wild snapdragon plant pedigree sheds light on limited connectivity enhanced by higher migrant reproductive success in a fragmented landscape.

Background: In contrast with historical knowledge, a recent view posits that a non-negligible proportion of populations thrive in a fragmented landscape. One underlying mechanism is the maintenance of functional connectivity, i.e., the net flow of individuals or their genes moving among suitable habitat patches. Alternatively, functional connectivity might be typically limited but enhanced by a higher reproductive success of migrants. Methods: We tested for this hypothesis in wild snapdragon plants inhabiting six patches separated by seawater in a fragmented Mediterranean scrubland landscape. We reconstructed their pedigree by using a parentage assignment method based on microsatellite genetic markers. We then estimated functional connectivity and the reproductive success of plants resulting from between-patch dispersal events. Results: We found that wild snapdragon plants thrived in this fragmented landscape, although functional connectivity between habitat patches was low (i.e. 2.9%). The progeny resulting from between-patch dispersal events had a higher reproductive success than residents. Conclusion: Our findings imply that low functional connectivity in a fragmented landscapes may have been enhanced by higher reproductive success after migration. This original mechanisms might be partly compensating the negative impact of fragmentation.

Potential adaptive divergence between subspecies and populations of snapdragon plants inferred from QST -FST comparisons.

Phenotypic divergence among natural populations can be explained by natural selection or by neutral processes such as drift. Many examples in the literature compare putatively neutral (FST ) and quantitative genetic (QST ) differentiation in multiple populations to assess their evolutionary signature and identify candidate traits involved with local adaptation. Investigating these signatures in closely related or recently diversified species has the potential to shed light on the divergence processes acting at the interspecific level. Here, we conducted this comparison in two subspecies of snapdragon plants (eight populations of Antirrhinum majus pseudomajus and five populations of A. m. striatum) in a common garden experiment. We also tested whether altitude was involved with population phenotypic divergence. Our results identified candidate phenological and morphological traits involved with local adaptation. Most of these traits were identified in one subspecies but not the other. Phenotypic divergence increased with altitude for a few biomass-related traits, but only in A. m. striatum. These traits therefore potentially reflect A. m. striatum adaptation to altitude. Our findings imply that adaptive processes potentially differ at the scale of A. majus subspecies.

Intraspecific difference among herbivore lineages and their host-plant specialization drive the strength of trophic cascades.

Trophic cascades - the indirect effect of predators on non-adjacent lower trophic levels - are important drivers of the structure and dynamics of ecological communities. However, the influence of intraspecific trait variation on the strength of trophic cascade remains largely unexplored, which limits our understanding of the mechanisms underlying ecological networks. Here we experimentally investigated how intraspecific difference among herbivore lineages specialized on different host plants influences trophic cascade strength in a terrestrial tri-trophic system. We found that the occurrence and strength of the trophic cascade are strongly influenced by herbivores' lineage and host-plant specialization but are not associated with density-dependent effects mediated by the growth rate of herbivore populations. Our findings stress the importance of intraspecific heterogeneities and evolutionary specialization as drivers of trophic cascade strength and underline that intraspecific variation should not be overlooked to decipher the joint influence of evolutionary and ecological factors on the functioning of multi-trophic interactions.