Lichen holobionts show compositional structure along elevation.

Holobionts are dynamic ecosystems that may respond to abiotic drivers with compositional changes. Uncovering elevational diversity patterns within these microecosystems can further our understanding of community-environment interactions. Here, we assess how the major components of lichen holobionts-fungal hosts, green algal symbionts, and the bacterial community-collectively respond to an elevational gradient. We analyse populations of two lichen symbioses, Umbilicaria pustulata and U. hispanica, along an elevational gradient spanning 2100 altitudinal metres and covering three major biomes. Our study shows (i) discontinuous genomic variation in fungal hosts with one abrupt genomic differentiation within each of the two host species, (ii) altitudinally structured bacterial communities with pronounced turnover within and between hosts, and (iii) altitude-specific presence of algal symbionts. Alpha diversity of bacterial communities decreased with increasing elevation. A marked turnover in holobiont diversity occurred across two altitudinal belts: at 11°C-13°C average annual temperature (here: 800-1200 m a.s.l.), and at 7°C-9°C average annual temperature (here: 1500-1800 m a.s.l.). The two observed zones mark a clustering of distribution limits and community shifts. The three ensuing altitudinal classes, that is, the most frequent combinations of species in holobionts, approximately correspond to the Mediterranean, cool-temperate, and alpine climate zones. We conclude that multitrophic microecosystems, such as lichen holobionts, respond with concerted compositional changes to climatic factors that also structure communities of macroorganisms, for example, vascular plants.

Transposable Elements in the Genome of the Lichen-Forming Fungus Umbilicaria pustulata and Their Distribution in Different Climate Zones along Elevation.

Transposable elements (TEs) are an important source of genome plasticity across the tree of life. Drift and natural selection are important forces shaping TE distribution and accumulation. Fungi, with their multifaceted phenotypic diversity and relatively small genome size, are ideal models to study the role of TEs in genome evolution and their impact on the host's ecological and life history traits. Here we present an account of all TEs found in a high-quality reference genome of the lichen-forming fungus Umbilicaria pustulata, a macrolichen species comprising two climatic ecotypes: Mediterranean and cold temperate. We trace the occurrence of the newly identified TEs in populations along three elevation gradients using a Pool-Seq approach to identify TE insertions of potential adaptive significance. We found that TEs cover 21.26% of the 32.9 Mbp genome, with LTR Gypsy and Copia clades being the most common TEs. We identified 28 insertions displaying consistent insertion frequency differences between the two host ecotypes across the elevation gradients. Most of the highly differentiated insertions were located near genes, indicating a putative function. This pioneering study of the content and climate niche-specific distribution of TEs in a lichen-forming fungus contributes to understanding the roles of TEs in fungal evolution.

Expanding the mutualistic niche: parallel symbiont turnover along climatic gradients.

Keystone mutualisms, such as corals, lichens or mycorrhizae, sustain fundamental ecosystem functions. Range dynamics of these symbioses are, however, inherently difficult to predict because host species may switch between different symbiont partners in different environments, thereby altering the range of the mutualism as a functional unit. Biogeographic models of mutualisms thus have to consider both the ecological amplitudes of various symbiont partners and the abiotic conditions that trigger symbiont replacement. To address this challenge, we here investigate 'symbiont turnover zones'--defined as demarcated regions where symbiont replacement is most likely to occur, as indicated by overlapping abundances of symbiont ecotypes. Mapping the distribution of algal symbionts from two species of lichen-forming fungi along four independent altitudinal gradients, we detected an abrupt and consistent ß-diversity turnover suggesting parallel niche partitioning. Modelling contrasting environmental response functions obtained from latitudinal distributions of algal ecotypes consistently predicted a confined altitudinal turnover zone. In all gradients this symbiont turnover zone is characterized by approximately 12°C average annual temperature and approximately 5°C mean temperature of the coldest quarter, marking the transition from Mediterranean to cool temperate bioregions. Integrating the conditions of symbiont turnover into biogeographic models of mutualisms is an important step towards a comprehensive understanding of biodiversity dynamics under ongoing environmental change.

Evolutionary genomics can improve prediction of species' responses to climate change.

Global climate change (GCC) increasingly threatens biodiversity through the loss of species, and the transformation of entire ecosystems. Many species are challenged by the pace of GCC because they might not be able to respond fast enough to changing biotic and abiotic conditions. Species can respond either by shifting their range, or by persisting in their local habitat. If populations persist, they can tolerate climatic changes through phenotypic plasticity, or genetically adapt to changing conditions depending on their genetic variability and census population size to allow for de novo mutations. Otherwise, populations will experience demographic collapses and species may go extinct. Current approaches to predicting species responses to GCC begin to combine ecological and evolutionary information for species distribution modelling. Including an evolutionary dimension will substantially improve species distribution projections which have not accounted for key processes such as dispersal, adaptive genetic change, demography, or species interactions. However, eco-evolutionary models require new data and methods for the estimation of a species' adaptive potential, which have so far only been available for a small number of model species. To represent global biodiversity, we need to devise large-scale data collection strategies to define the ecology and evolutionary potential of a broad range of species, especially of keystone species of ecosystems. We also need standardized and replicable modelling approaches that integrate these new data to account for eco-evolutionary processes when predicting the impact of GCC on species' survival. Here, we discuss different genomic approaches that can be used to investigate and predict species responses to GCC. This can serve as guidance for researchers looking for the appropriate experimental setup for their particular system. We furthermore highlight future directions for moving forward in the field and allocating available resources more effectively, to implement mitigation measures before species go extinct and ecosystems lose important functions.

Understanding Maladaptation by Uniting Ecological and Evolutionary Perspectives.

Evolutionary biologists have long trained their sights on adaptation, focusing on the power of natural selection to produce relative fitness advantages while often ignoring changes in absolute fitness. Ecologists generally have taken a different tack, focusing on changes in abundance and ranges that reflect absolute fitness while often ignoring relative fitness. Uniting these perspectives, we articulate various causes of relative and absolute maladaptation and review numerous examples of their occurrence. This review indicates that maladaptation is reasonably common from both perspectives, yet often in contrasting ways. That is, maladaptation can appear strong from a relative fitness perspective, yet populations can be growing in abundance. Conversely, resident individuals can appear locally adapted (relative to nonresident individuals) yet be declining in abundance. Understanding and interpreting these disconnects between relative and absolute maladaptation, as well as the cases of agreement, is increasingly critical in the face of accelerating human-mediated environmental change. We therefore present a framework for studying maladaptation, focusing in particular on the relationship between absolute and relative fitness, thereby drawing together evolutionary and ecological perspectives. The unification of these ecological and evolutionary perspectives has the potential to bring together previously disjunct research areas while addressing key conceptual issues and specific practical problems.

Causes of maladaptation.

Evolutionary biologists tend to approach the study of the natural world within a framework of adaptation, inspired perhaps by the power of natural selection to produce fitness advantages that drive population persistence and biological diversity. In contrast, evolution has rarely been studied through the lens of adaptation's complement, maladaptation. This contrast is surprising because maladaptation is a prevalent feature of evolution: population trait values are rarely distributed optimally; local populations often have lower fitness than imported ones; populations decline; and local and global extinctions are common. Yet we lack a general framework for understanding maladaptation; for instance in terms of distribution, severity, and dynamics. Similar uncertainties apply to the causes of maladaptation. We suggest that incorporating maladaptation-based perspectives into evolutionary biology would facilitate better understanding of the natural world. Approaches within a maladaptation framework might be especially profitable in applied evolution contexts - where reductions in fitness are common. Toward advancing a more balanced study of evolution, here we present a conceptual framework describing causes of maladaptation. As the introductory article for a Special Feature on maladaptation, we also summarize the studies in this Issue, highlighting the causes of maladaptation in each study. We hope that our framework and the papers in this Special Issue will help catalyze the study of maladaptation in applied evolution, supporting greater understanding of evolutionary dynamics in our rapidly changing world.

Evolutionary Rates Standardized for Evolutionary Space: Perspectives on Trait Evolution.

Characterization of evolutionary radiations benefits from describing the temporal patterns of trait disparification. Comparative methods attempt this by evaluating the statistical fit of trait distributions to a phylogenetic hypothesis under assumed evolutionary models. However, it can be challenging to differentiate between models, with discriminatory power depending on the modes of evolution underlying trait distributions. We suggest rates of 'trait space saturation', standardized for limits to evolutionary change, as an additional tool to distinguish between modes of trait evolution. We evaluate this approach using simulations and show that trait space saturation can identify the true model of trait evolution in cases where traditional comparative methods can fail. We illustrate our approach using diverse empirical studies that represent contrasting scenarios of evolutionary radiation.

Environment and host identity structure communities of green algal symbionts in lichens.

An understanding of how biotic interactions shape species' distributions is central to predicting host-symbiont responses under climate change. Switches to locally adapted algae have been proposed to be an adaptive strategy of lichen-forming fungi to cope with environmental change. However, it is unclear how lichen photobionts respond to environmental gradients, and whether they play a role in determining the fungal host's upper and lower elevational limits. Deep-coverage Illumina DNA metabarcoding was used to track changes in the community composition of Trebouxia algae associated with two phylogenetically closely related, but ecologically divergent fungal hosts along a steep altitudinal gradient in the Mediterranean region. We detected the presence of multiple Trebouxia species in the majority of thalli. Both altitude and host genetic identity were strong predictors of photobiont community assembly in these two species. The predominantly clonally dispersing fungus showed stronger altitudinal structuring of photobiont communities than the sexually reproducing host. Elevation ranges of the host were not limited by the lack of compatible photobionts. Our study sheds light on the processes guiding the formation and distribution of specific fungal-algal combinations in the lichen symbiosis. The effect of environmental filtering acting on both symbiotic partners appears to shape the distribution of lichens.

Season-Long Experimental Drought Alters Fungal Community Composition but Not Diversity in a Grassland Soil.

Using terrestrial model ecosystems (TMEs), we investigated how reduced moisture conditions impact soil fungal communities from a temperate grassland over the course of an entire season. Starting at about 65% of the soil's maximum water holding capacity (WHCmax), TME soils were adjusted to three moisture levels for 15 weeks: 70% WHCmax, approximating starting conditions, 50% WHCmax, and 30% WHCmax, representing reduced moisture conditions. Diversity and abundances of soil fungi at the start and at the end of the experiment were characterized using Illumina meta-barcoding. Community diversity at the end of the experiment did not differ between experimental moisture levels and was comparable to diversity measures from the field. However, fungal communities did change compositionally in both abundances and presence/absence of species. Analyzing class-level and individual contributions of fungi to these changes revealed that only a minor portion reacted significantly, indicating that most compositional change was likely driven by many consistent small-scale shifts in presence/absences or abundances. Together, our results show that prolonged reduction in soil moisture conditions will trigger compositional changes in soil fungal communities but not necessarily change overall diversity. We highlight the cumulative contribution of minor but consistent changes among community members, as opposed to significant responses of individual species. We also detected a strong general experimental effect on soil fungi that are moved from the field to experimental TMEs, suggesting the importance of acclimatization effects in these communities under laboratory conditions.

How Parallel Is Parallel Evolution? A Comparative Analysis in Fishes.

Evidence of phenotypic parallelism is often used to infer the deterministic role played by natural selection. However, variation in the extent or direction of divergence is often evident among independent evolutionary replicates, raising the following question: just how parallel, overall, is parallel evolution? We answer this question through a comparative analysis of studies of fishes, a taxon where parallel evolution has been much discussed. We first ask how much of the among-population variance in phenotypic traits can be explained by different "environment" types, such as high predation versus low predation or benthic versus limnetic. We then use phenotypic change vector analysis to quantify variation in the direction (vector angles) and magnitude (vector lengths) of environment-associated divergence. All analyses show high variation in the extent of parallelism-from very high to very low, along with everything in between-highlighting the importance of quantifying parallelism rather than just asserting its presence. Interestingly, instances of low extents of parallelism represent important components of divergence in many cases, promising considerable opportunities for inferences about the factors shaping phenotypic divergence.

Adaptive differentiation coincides with local bioclimatic conditions along an elevational cline in populations of a lichen-forming fungus.

BACKGROUND: Many fungal species occur across a variety of habitats. Particularly lichens, fungi forming symbioses with photosynthetic partners, have evolved remarkable tolerances for environmental extremes. Despite their ecological importance and ubiquity, little is known about the genetic basis of adaption in lichen populations. Here we studied patterns of genome-wide differentiation in the lichen-forming fungus Lasallia pustulata along an altitudinal gradient in the Mediterranean region. We resequenced six populations as pools and identified highly differentiated genomic regions. We then detected gene-environment correlations while controlling for shared population history and pooled sequencing bias, and performed ecophysiological experiments to assess fitness differences of individuals from different environments. RESULTS: We detected two strongly differentiated genetic clusters linked to Mediterranean and temperate-oceanic climate, and an admixture zone, which coincided with the transition between the two bioclimates. High altitude individuals showed ecophysiological adaptations to wetter and more shaded conditions. Highly differentiated genome regions contained a number of genes associated with stress response, local environmental adaptation, and sexual reproduction. CONCLUSIONS: Taken together our results provide evidence for a complex interplay between demographic history and spatially varying selection acting on a number of key biological processes, suggesting a scenario of ecological speciation.

Can Evolution Supply What Ecology Demands?

A simplistic view of the adaptive process pictures a hillside along which a population can climb: when ecological 'demands' change, evolution 'supplies' the variation needed for the population to climb to a new peak. Evolutionary ecologists point out that this simplistic view can be incomplete because the fitness landscape changes dynamically as the population evolves. Geneticists meanwhile have identified complexities relating to the nature of genetic variation and its architecture, and the importance of epigenetic variation is under debate. In this review, we highlight how complexity in both ecological 'demands' and the evolutionary 'supply' influences organisms' ability to climb fitness landscapes that themselves change dynamically as evolution proceeds, and encourage new synthetic effort across research disciplines towards ecologically realistic studies of adaptation.

Adaptive capabilities and fitness consequences associated with pollution exposure in fish.

Many fish populations are exposed to harmful levels of chemical pollution and selection pressures associated with these exposures have led to the evolution of tolerance. Our understanding of the physiological basis for these adaptations is limited, but they are likely to include processes involved with the absorption, distribution, metabolism and/or excretion of the target chemical. Other potential adaptive mechanisms include enhancements in antioxidant responses, an increased capacity for DNA and/or tissue repair and alterations to the life cycle of fish that enable earlier reproduction. Analysis of single-nucleotide polymorphism frequencies has shown that tolerance to hydrocarbon pollutants in both marine and estuarine fish species involves alteration in the expression of the xenobiotic metabolism enzyme CYP1A. In this review, we present novel data showing also that variants of the CYP1A gene have been under selection in guppies living in Trinidadian rivers heavily polluted with crude oil. Potential costs associated with these adaptations could reduce fitness in unpolluted water conditions. Integrating knowledge of local adaptation to pollution is an important future consideration in conservation practices such as for successful restocking, and improving connectivity within river systems.This article is part of the themed issue 'Human influences on evolution, and the ecological and societal consequences'.

Do stressful conditions make adaptation difficult? Guppies in the oil-polluted environments of southern Trinidad.

The ability of populations to rapidly adapt to new environments will determine their future in an increasingly human-modified world. Although meta-analyses do frequently uncover signatures of local adaptation, they also reveal many exceptions. We suggest that particular constraints on local adaptation might arise when organisms are exposed to novel stressors, such as anthropogenic pollution. To inform this possibility, we studied the extent to which guppies (Poecilia reticulata) show local adaptation to oil pollution in southern Trinidad. Neutral genetic markers revealed that paired populations in oil-polluted versus not-polluted habitats diverged independently in two different watersheds. Morphometrics revealed some divergence (particularly in head shape) between these environments, some of which was parallel between rivers. Reciprocal transplant experiments in nature, however, found little evidence of local adaptation based on survival and growth. Moreover, subsequent laboratory experiments showed that the two populations from oil-polluted sites showed only weak local adaptation even when compared to guppies from oil-free northern Trinidad. We conclude that guppies show little local adaptation to oil pollution, which might result from the challenges associated with adaptation to particularly stressful environments. It might also reflect genetic drift owing to small population sizes and/or high gene flow between environments.

When maladaptive gene flow does not increase selection.

Populations receiving high maladaptive gene flow are expected to experience strong directional selection-because gene flow pulls mean phenotypes away from local fitness peaks. We tested this prediction by means of a large and replicated mark-recapture study of threespine stickleback (Gasterosteus aculeatus) in two stream populations. One of the populations (outlet) experiences high gene flow from the lake population and its morphology is correspondingly poorly adapted. The other population (inlet) experiences very low gene flow from the lake population and its morphology is correspondingly well adapted. Contrary to the above prediction, selection was not stronger in the outlet than in the inlet, a result that forced us to consider potential reasons for why maladaptive gene flow might not increase selection. Of particular interest, we show by means of a simple population genetic model that maladaptive gene flow can-under reasonable conditions-decrease the strength of directional selection. This outcome occurs when immigrants decrease mean fitness in the resident population, which decreases the strength of selection against maladapted phenotypes. We argue that this previously unrecognized effect of gene flow deserves further attention in theoretical and empirical studies.

Host preference of an introduced 'generalist' parasite for a non-native host.

Parasites can invade new ecosystems if they are introduced with their native hosts or if they successfully infect and colonise new hosts upon arrival. Here, we ask to what extent an introduced parasite demonstrates specialisation among novel host species. Infection surveys across three field sites in Gatun Lake, Panama, revealed that the invasive peacock bass, Cichla monoculus, was more commonly infected by the introduced trematode parasite Centrocestus formosanus than were three other common cichlid fishes. Laboratory infection experiments were conducted to determine whether parasitism might be driven by differential encounter/exposure to parasites or by differential infection susceptibility/preference across different host species. These experiments were performed by controlling for parasite exposure in single host (compatibility) experiments and in mixed host (preference) experiments. In all cases, the peacock bass exhibited higher infection rates with viable metacercariae relative to the other potential fish hosts. Our experiments thus support that an introduced generalist parasite shows apparent specialisation on a specific novel host. Further studies are needed to determine whether these patterns of specialisation are the result of local adaptation following invasion by the parasite.

Linking macrotrends and microrates: Re-evaluating microevolutionary support for Cope's rule.

Cope's rule, wherein a lineage increases in body size through time, was originally motivated by macroevolutionary patterns observed in the fossil record. More recently, some authors have argued that evidence exists for generally positive selection on individual body size in contemporary populations, providing a microevolutionary mechanism for Cope's rule. If larger body size confers individual fitness advantages as the selection estimates suggest, thereby explaining Cope's rule, then body size should increase over microevolutionary time scales. We test this corollary by assembling a large database of studies reporting changes in phenotypic body size through time in contemporary populations, as well as studies reporting average breeding values for body size through time. Trends in body size were quite variable with an absence of any general trend, and many populations trended toward smaller body sizes. Although selection estimates can be interpreted to support Cope's rule, our results suggest that actual rates of phenotypic change for body size cannot. We discuss potential reasons for this discrepancy and its implications for the understanding of Cope's rule.

Individual differences in migratory behavior shape population genetic structure and microhabitat choice in sympatric blackcaps (Sylvia atricapilla).

In migratory birds, traits such as orientation and distance are known to have a strong genetic background, and they often exhibit considerable within-population variation. How this variation relates to evolutionary responses to ongoing selection is unknown because the underlying mechanisms that translate environmental changes into population genetic changes are unclear. We show that within-population genetic structure in southern German blackcaps (Sylvia atricapilla) is related to individual differences in migratory behavior. Our 3-year study revealed a positive correlation between individual migratory origins, denoted via isotope (δ (2)H) values, and genetic distances. Genetic diversity and admixture differed not only across a recently established migratory polymorphism with NW- and SW-migrating birds but also across δ (2)H clusters within the same migratory route. Our results suggest assortment based on individual migratory origins which would facilitate evolutionary responses. We scrutinized arrival times and microhabitat choice as potential mechanisms mediating between individual variation in migratory behavior and assortment. We found significant support that microhabitat choice, rather than timing of arrival, is associated with individual variation in migratory origins. Moreover, examining genetic diversity across the migratory divide, we found migrants following the NW route to be genetically more distinct from each other compared with migrants following the traditional SW route. Our study suggests that migratory behavior shapes population genetic structure in blackcaps not only across the migratory divide but also on an individual level independent of the divide. Thus, within-population variation in migratory behavior might play an important role in translating environmental change into genetic change.

Contrasting patterns of genetic differentiation among Blackcaps (Sylvia atricapilla) with divergent migratory orientations in Europe.

Migratory divides are thought to facilitate behavioral, ecological, and genetic divergence among populations with different migratory routes. However, it is currently contentious how much genetic divergence is needed to maintain distinct migratory behavior across migratory divides. Here we investigate patterns of neutral genetic differentiation among Blackcap (Sylvia atricapilla) populations with different migratory strategies across Europe. We compare the level of genetic divergence of populations migrating to southwestern (SW) or southeastern (SE) wintering areas with birds wintering in the British Isles following a recently established northwesterly (NW) migration route. The migratory divide between SW and SE wintering areas can be interpreted as a result of a re-colonization process after the last glaciation. Thus we predicted greater levels of genetic differentiation among the SW/SE populations. However, a lack of genetic differentiation was found between SW and SE populations, suggesting that interbreeding likely occurs among Blackcaps with different migratory orientations across a large area; therefore the SW/SE migratory divide can be seen as diffuse, broad band and is, at best, a weak isolating barrier. Conversely, weak, albeit significant genetic differentiation was evident between NW and SW migrants breeding sympatrically in southern Germany, suggesting a stronger isolating mechanism may be acting in this population. Populations located within/near the SW/SE contact zone were the least genetically divergent from NW migrants, confirming NW migrants likely originated from within the contact zone. Significant isolation-by-distance was found among eastern Blackcap populations (i.e. SE migrants), but not among western populations (i.e. NW and SW migrants), revealing different patterns of genetic divergence among Blackcap populations in Europe. We discuss possible explanations for the genetic structure of European Blackcaps and how gene flow influences the persistence of divergent migratory behaviors.