Contrasting post-settlement selection results in many-to-one mapping of high performance phenotypes in the Hawaiian waterfall-climbing goby <Emphasis Type="Italic">Sicyopterus stimpsoni</Emphasis>

Natural selection drives adaptive evolution, but contrasting environmental pressures may lead to trade-offs between phenotypes that confer different performances. Such trade-offs may weaken the strength of selection and/or generate complex fitness surfaces with multiple local optima that correspond to different selection regimes. We evaluated how differences in patterns of phenotypic selection might promote morphological differences between subpopulations of the amphidromous Hawaiian waterfall-climbing goby, Sicyopterus stimpsoni. We conducted laboratory experiments on fish from the islands of Kaua‘i and Hawai‘i (the “Big Island”) to compare patterns of linear and nonlinear selection, and the opportunity for selection, that result from two contrasting pressures, predator evasion and waterfall climbing, which vary in intensity between islands. We found directional and nonlinear selection were strongest when individuals were exposed to their primary selective pressures (predator evasion on Kaua‘i, waterfall climbing on the Big Island). However, the opportunity for selection was greater for the non-primary pressure: climbing on Kaua‘i, predator evasion on the Big Island. Canonical rotation of the nonlinear gamma matrix demonstrated that individuals from Kaua‘i and the Big Island occupy regions near their local fitness peaks for some traits. Therefore, selection for predator evasion on Kaua‘i and climbing on the Big Island may be less effective in promoting morphological changes in this species, because variation of functionally important traits in their respective environments may have been reduced by directional or stabilizing selection. These results demonstrate that despite constraints on the opportunities for selection, population differences in phenotypic traits can arise due to differences in selective regimes. For S. stimpsoni, sufficient variation exists in other locomotor traits, allowing for necessary levels of performance in the contrasting selective regime (i.e., climbing on Kaua‘i and predator evasion on the Big Island) through many-to-one-mapping, which may be essential for the survival of local populations in an evanescent island environment.     

AN ADAPTIVE RADIATION OF FROGS IN A SOUTHEAST ASIAN ISLAND ARCHIPELAGO

Living amphibians exhibit a diversity of ecologies, life histories, and species‐rich lineages that offers opportunities for studies of adaptive radiation. We characterize a diverse clade of frogs (Kaloula, Microhylidae) in the Philippine island archipelago as an example of an adaptive radiation into three primary habitat specialists or ecotypes. We use a novel phylogenetic estimate for this clade to evaluate the tempo of lineage accumulation and morphological diversification. Because species‐level phylogenetic estimates for Philippine Kaloula are lacking, we employ dense population sampling to determine the appropriate evolutionary lineages for diversification analyses. We explicitly take phylogenetic uncertainty into account when calculating diversification and disparification statistics and fitting models of diversification. Following dispersal to the Philippines from Southeast Asia, Kaloula radiated rapidly into several well‐supported clades. Morphological variation within Kaloula is partly explained by ecotype and accumulated at high levels during this radiation, including within ecotypes. We pinpoint an axis of morphospace related directly to climbing and digging behaviors and find patterns of phenotypic evolution suggestive of ecological opportunity with partitioning into distinct habitat specialists. We conclude by discussing the components of phenotypic diversity that are likely important in amphibian adaptive radiations.     

Making pore choices: repeated regime shifts in stomatal ratio

Ecologically important traits do not evolve without limits. Instead, evolution is constrained by the set of available and viable phenotypes. In particular, natural selection may only favour a narrow range of adaptive optima constrained within selective regimes. Here, I integrate data with theory to test whether selection explains phenotypic constraint. A global database of 599 plant species from 94 families shows that stomatal ratio, a trait affecting photosynthesis and defence against pathogens, is highly constrained. Most plants have their stomata on the lower leaf surface (hypostomy), but species with half their stomata on each surface (amphistomy) form a distinct mode in the trait distribution. A model based on a trade-off between maximizing photosynthesis and a fitness cost of upper stomata predicts a limited number of adaptive solutions, leading to a multimodal trait distribution. Phylogenetic comparisons show that amphistomy is the most common among fast-growing species, supporting the view that CO₂ diffusion is under strong selection. These results indicate that selective optima stay within a relatively stable set of selective regimes over macroevolutionary time.     

The effect of parity on morphological evolution among phrynosomatid lizards

The shift from egg laying to live‐bearing is one of the most well‐studied transitions in evolutionary biology. Few studies, however, have assessed the effect of this transition on morphological evolution. Here, we evaluated the effect of reproductive mode on the morphological evolution of 10 traits, among 108 species of phrynosomatid lizards. We assess whether the requirement for passing shelled eggs through the pelvic girdle has led to morphological constraints in oviparous species and whether long gestation times in viviparous species have led to constraints in locomotor morphology. We fit models to the data that vary both in their tempo (strength and rate of selection) and mode of evolution (Brownian or Ornstein‐Uhlenbeck) and estimates of trait optima. We found that most traits are best fit by a generalized multipeak OU model, suggesting differing trait optima for viviparous vs. oviparous species. Additionally, rates (σ²) of both pelvic girdle and forelimb trait evolution varied with parity; viviparous species had higher rates. Hindlimb traits, however, exhibited no difference in σ² between parity modes. In a functional context, our results suggest that the passage of shelled eggs constrains the morphology of the pelvic girdle, but we found no evidence of morphological constraint of the locomotor apparatus in viviparous species. Our results are consistent with recent lineage diversification analyses, leading to the conclusion that transitions to viviparity increase both lineage and morphological diversification.     

Adaptive evolution of reproductive and vegetative traits driven by breeding systems

The evolution of inflorescence size, a key trait in reproductive success, was studied in the genus Acer under a perspective of adaptive evolution. Breeding systems, hypothesized to indicate different levels of mating competition, were considered as the selective scenarios defining different optima of inflorescence size. Larger inflorescences, which increase male fitness by generating larger floral displays, were hypothesized to be selected under scenarios with higher competition with unisexuals. An identical approach was used to test if the same selective regimes could be driving the evolution of leaf size, a vegetative trait that was found to be correlated with inflorescence size. A Brownian motion model of inflorescence/leaf-size evolution (which cannot distinguish between changes caused by pure drift processes and changes caused by natural selection in rapidly and randomly changing environments) was compared with several adaptive Ornstein-Uhlenbeck (OU) models, which can quantify the effects of both stochasticity and natural selection. The best-fitting model for inflorescence/leaf-size evolution was an OU model with three optima that increased with the level of mating competition. Both traits evolved under the same selective regimes and in the same direction, confirming a pattern of correlated evolution. These results show that a selective regime hypothetically related to the evolution of a reproductive trait can also explain the evolution of a vegetative trait.     

The causes of variation in tree seedling traits: the roles of environmental selection versus chance

A key aspect of biodiversity is the great quantitative variation in functional traits observed among species. One perspective asserts that trait values should converge on a single optimum value in a particular selective environment, and consequently trait variation would reflect differences in selective environment, and evolutionary outcomes would be predictable. An alternative perspective asserts that there are likely multiple alternative optima within a particular selective environment, and consequently different lineages would evolve toward different optima due to chance. Because there is evidence for both of these perspectives, there is a long-standing controversy over the relative importance of convergence due to environmental selection versus divergence due to chance in shaping trait variation. Here, I use a model of tree seedling growth and survival to distinguish trait variation associated with multiple alternative optima from variation associated with environmental differences. I show that variation in whole plant traits is best explained by environmental differences, whereas in organ level traits variation is more affected by alternative optima. Consequently, I predict that in nature variation in organ level traits is most closely related to phylogeny, whereas variation in whole plant traits is most closely related to ecology.     

Optimal foraging, specialization, and a solution to Liem's paradox

Species that appear highly specialized on the basis of their phenotype (e.g., morphology, behavior, and physiology) also sometimes act as ecological generalists. This apparent paradox has been used to argue against the importance of competition as a diversifying evolutionary force. We provide an alternative explanation based on optimal foraging theory. Some resources are intrinsically easy to use and are widely preferred, while others require specialized phenotypic traits on the part of the consumer. This asymmetry allows optimally foraging consumers to evolve phenotypic specializations on nonpreferred resources without greatly compromising their ability to use preferred resources. The evolution of phenotypic specialization on nonpreferred resources can be driven by competition, but the specialists act as ecological generalists whenever their preferred resources are available. Our model identifies at least three different concepts of specialization that need to be distinguished, based on diet, prey utilization efficiencies, and phenotypic adaptations. The relationships among these concepts are complex and often counterintuitive. Specialists should often reject the very resources that they have evolved traits to use. The most extreme phenotypic specializations should occur in the absence of a trade-off between using preferred and nonpreferred resources. Our model may explain why extreme phenotypic-specializations evolve more often in fish communities than in terrestrial vertebrate communities and provides a mechanism whereby species can coexist in stable communities despite common preferences for some resources.     

Testing limits to adaptation along altitudinal gradients in rainforest Drosophila

Given that evolution can generate rapid and dramatic shifts in the ecological tolerance of a species, what prevents populations adapting to expand into new habitat at the edge of their distributions? Recent population genetic models have focused on the relative costs and benefits of migration between populations. On the one hand, migration may limit adaptive divergence by preventing local populations from matching their local selective optima. On the other hand, migration may also contribute to the genetic variance necessary to allow populations to track these changing optima. Empirical evidence for these contrasting effects of gene flow in natural situations are lacking, largely because it remains difficult to acquire. Here, we develop a way to explore theoretical models by estimating genetic divergence in traits that confer stress resistance along similar ecological gradients in rainforest Drosophila. This approach allows testing for the coupling of clinal divergence with local density, and the effects of genetic variance and the rate of change of the optimum on the response to selection. In support of a swamping effect of migration on phenotypic divergence, our data show no evidence for a cline in stress-related traits where the altitudinal gradient is steep, but significant clinal divergence where it is shallow. However, where clinal divergence is detected, sites showing trait means closer to the presumed local optimum have more genetic variation than sites with trait means distant from their local optimum. This pattern suggests that gene flow also aids a sustained response to selection.     

Fixed‐effect variance and the estimation of repeatabilities and heritabilities: issues and solutions

Linear mixed‐effects models are frequently used for estimating quantitative genetic parameters, including the heritability, as well as the repeatability, of traits. Heritability acts as a filter that determines how efficiently phenotypic selection translates into evolutionary change, whereas repeatability informs us about the individual consistency of phenotypic traits. As quantities of biological interest, it is important that the denominator, the phenotypic variance in both cases, reflects the amount of phenotypic variance in the relevant ecological setting. The current practice of quantifying heritabilities and repeatabilities from mixed‐effects models frequently deprives their denominator of variance explained by fixed effects (often leading to upward bias of heritabilities and repeatabilities), and it has been suggested to omit fixed effects when estimating heritabilities in particular. We advocate an alternative option of fitting models incorporating all relevant effects, while including the variance explained by fixed effects into the estimation of the phenotypic variance. The approach is easily implemented and allows optimizing the estimation of phenotypic variance, for example by the exclusion of variance arising from experimental design effects while still including all biologically relevant sources of variation. We address the estimation and interpretation of heritabilities in situations in which potential covariates are themselves heritable traits of the organism. Furthermore, we discuss complications that arise in generalized and nonlinear mixed models with fixed effects. In these cases, the variance parameters on the data scale depend on the location of the intercept and hence on the scaling of the fixed effects. Integration over the biologically relevant range of fixed effects offers a preferred solution in those situations.     

Divergent ontogenies of trophic morphology in two closely related haplochromine cichlids

Fish develop morphological specializations in their trophic and locomotor systems as a result of varying functional demands in response to environmental pressures at different life stages. These specializations should maximize particular performances in specialists, adapting them to their trophic and habitat niches at each ontogenetic stage. Because differential growth rates of the structural components comprised in the head are likely to be linked to the diet of a fish throughout its development, we investigated the ontogenetic development of two haplochromine cichlid species belonging to different trophic guilds. We employed geometric morphometric techniques to evaluate whether starting from morphologically similar fry they diverge into phenotypes that characterize trophic guilds and locomotor types. Our examination of overall body shape shows that certain specialized morphological features are already present in fry, whereas other traits diverge through ontogeny due to differences in species‐specific allometric variation. Allometric shape variation was found to be more relevant for the biter specialist than for the sucker morphotype. Our results confirm that phenotypic changes during ontogeny can be linked to dietary and habitat shifts in these fish. Furthermore, evidence for an integrated development of trophic and locomotor specializations in morphology was observed. J. Morphol. 276:860–871, 2015. © 2015 Wiley Periodicals, Inc.     

Testing the stages model in the adaptive radiation of cichlid fishes in East African Lake Tanganyika

Adaptive radiation (AR) is a key process in the origin of organismal diversity. However, the evolution of trait disparity in connection with ecological specialization is still poorly understood. Available models for vertebrate ARs predict that diversification occurs in the form of temporal stages driven by different selective forces. Here, we investigate the AR of cichlid fishes in East African Lake Tanganyika and use macroevolutionary model fitting to evaluate whether diversification happened in temporal stages. Six trait complexes, for which we also provide evidence of their adaptiveness, are analysed with comparative methods: body shape, pharyngeal jaw shape, gill raker traits, gut length, brain weight and body coloration. Overall, we do not find strong evidence for the ‘stages model’ of AR. However, our results suggest that trophic traits diversify earlier than traits implicated in macrohabitat adaptation and that sexual communication traits (i.e. coloration) diversify late in the radiation.     

Phylogenetic patterns of trait and trait plasticity evolution: Insights from amphibian embryos

Environmental variation favors the evolution of phenotypic plasticity. For many species, we understand the costs and benefits of different phenotypes, but we lack a broad understanding of how plastic traits evolve across large clades. Using identical experiments conducted across North America, we examined prey responses to predator cues. We quantified five life‐history traits and the magnitude of their plasticity for 23 amphibian species/populations (spanning three families and five genera) when exposed to no cues, crushed‐egg cues, and predatory crayfish cues. Embryonic responses varied considerably among species and phylogenetic signal was common among the traits, whereas phylogenetic signal was rare for trait plasticities. Among trait‐evolution models, the Ornstein–Uhlenbeck (OU) model provided the best fit or was essentially tied with Brownian motion. Using the best fitting model, evolutionary rates for plasticities were higher than traits for three life‐history traits and lower for two. These data suggest that the evolution of life‐history traits in amphibian embryos is more constrained by a species’ position in the phylogeny than is the evolution of life history plasticities. The fact that an OU model of trait evolution was often a good fit to patterns of trait variation may indicate adaptive optima for traits and their plasticities.     

Intraspecific trait variation of woody species reduced in a savanna community,southwest China

Plants deploy various ecological strategies in response to environmental heterogeneity. In many forest ecosystems, plants have been reported to have notable inter- and intra-specific trait variation, as well as clear phylogenetic signals, indicating that these species possess a degree of phenotypic plasticity to cope with habitat variation in the community. Savanna communities, however, grow in an open canopy structure and exhibit little species diversification, likely as a result of strong environmental stress. In this study, we hypothesized that the phylogenetic signals of savanna species would be weak, the intraspecific trait variation (ITV) would be low, and the contribution of intraspecific variation to total trait variance would be reduced, owing to low species richness, multiple stresses and relatively homogenous community structure. To test these hypotheses, we sampled dominant woody species in a dry-hot savanna in southwestern China, focusing on leaf traits related to adaptability of plants to harsh conditions (year-round intense radiation, low soil fertility and seasonal droughts). We found weak phylogenetic signals in leaf traits and low ITV (at both individual and canopy-layer levels). Intraspecific variation (including leaf-, layer- and individual-scales) contributed little to the total trait variance, whereas interspecific variation and variation in leaf phenology explained substantial variance. Our study suggests that intraspecific trait variation is reduced in savanna community. Furthermore, our findings indicate that classifying species by leaf phenology may help better understand how species coexist under similar habitats with strong stresses.     

Quantitative trait loci affecting survival and fertility-related traits in Caenorhabditis elegans show genotype-environment interactions, pleiotropy and epistasis.

We have identified, using composite interval mapping, quantitative trait loci (QTL) affecting a variety of life history traits (LHTs) in the nematode Caenorhabditis elegans. Using recombinant inbred strains assayed on the surface of agar plates, we found QTL for survival, early fertility, age of onset of sexual maturity, and population growth rate. There was no overall correlation between survival on solid media and previous measures of survival in liquid media. Of the four survival QTL found in these two environments, two have genotype-environment interactions (GEIs). Epistatic interactions between markers were detected for four traits. A multiple regression approach was used to determine which single markers and epistatic interactions best explained the phenotypic variance for each trait. The amount of phenotypic variance accounted for by genetic effects ranged from 13% (for internal hatching) to 46% (for population growth). Epistatic effects accounted for 9-11% of the phenotypic variance for three traits. Two regions containing QTL that affected more than one fertility-related trait were found. This study serves as an example of the power of QTL mapping for dissecting the genetic architecture of a suite of LHTs and indicates the potential importance of environment and GEIs in the evolution of this architecture.     

Phenotypic plasticity is not affected by experimental evolution in constant,predictable or unpredictable fluctuating thermal environments

The selective past of populations is presumed to affect the levels of phenotypic plasticity. Experimental evolution at constant temperatures is generally expected to lead to a decreased level of plasticity due to presumed costs associated with phenotypic plasticity when not needed. In this study, we investigated the effect of experimental evolution in constant, predictable and unpredictable daily fluctuating temperature regimes on the levels of phenotype plasticity in several life history and stress resistance traits in Drosophila simulans. Contrary to the expectation, evolution in the different regimes did not affect the levels of plasticity in any of the traits investigated even though the populations from the different thermal regimes had evolved different stress resistance and fitness trait means. Although costs associated with phenotypic plasticity are known, our results suggest that the maintenance of phenotypic plasticity might come at low and negligible costs, and thus, the potential of phenotypic plasticity to evolve in populations exposed to different environmental conditions might be limited.     

The phenotypic and genetic covariance structure of drosphilid wings

Evolutionary constraint results from the interaction between the distribution of available genetic variation and the position of selective optima. The availability of genetic variance in multitrait systems, as described by the additive genetic variance-covariance matrix (G), has been the subject of recent attempts to assess the prevalence of genetic constraints. However, evolutionary constraints have not yet been considered from the perspective of the phenotypes available to multivariate selection, and whether genetic variance is present in all phenotypes potentially under selection. Determining the rank of the phenotypic variance-covariance matrix (P) to characterize the phenotypes available to selection, and contrasting it with the rank of G, may provide a general approach to determining the prevalence of genetic constraints. In a study of a laboratory population of Drosophila bunnanda from northern Australia we applied factor-analytic modeling to repeated measures of individual wing phenotypes to determine the dimensionality of the phenotypic space described by P. The phenotypic space spanned by the 10 wing traits had 10 statistically supported dimensions. In contrast, factor-analytic modeling of G estimated for the same 10 traits from a paternal half-sibling breeding design suggested G had fewer dimensions than traits. Statistical support was found for only five and two genetic dimensions, describing a total of 99% and 72% of genetic variance in wing morphology in females and males, respectively. The observed mismatch in dimensionality between P and G suggests that although selection might act to shift the intragenerational population mean toward any trait combination, evolution may be restricted to fewer dimensions.     

FORAGING TRAIT (CO)VARIANCES IN STICKLEBACK EVOLVE DETERMINISTICALLY AND DO NOT PREDICT TRAJECTORIES OF ADAPTIVE DIVERSIFICATION

How does natural selection shape the structure of variance and covariance among multiple traits, and how do (co)variances influence trajectories of adaptive diversification? We investigate these pivotal but open questions by comparing phenotypic (co)variances among multiple morphological traits across 18 derived lake‐dwelling populations of threespine stickleback, and their marine ancestor. Divergence in (co)variance structure among populations is striking and primarily attributable to shifts in the variance of a single key foraging trait (gill raker length). We then relate this divergence to an ecological selection proxy, to population divergence in trait means, and to the magnitude of sexual dimorphism within populations. This allows us to infer that evolution in (co)variances is linked to variation among habitats in the strength of resource‐mediated disruptive selection. We further find that adaptive diversification in trait means among populations has primarily involved shifts in gill raker length. The direction of evolutionary trajectories is unrelated to the major axes of ancestral trait (co)variance. Our study demonstrates that natural selection drives both means and (co)variances deterministically in stickleback, and strongly challenges the view that the (co)variance structure biases the direction of adaptive diversification predictably even over moderate time spans.     

Evolutionary transition between bee pollination and hummingbird pollination in Salvia: Comparing means,variances and covariances of corolla traits

Covariation among traits can modify the evolutionary trajectory of complex structures. This process is thought to operate at a microevolutionary scale, but its long‐term effects remain controversial because trait covariation can itself evolve. Flower morphology, and particularly floral trait (co)variation, has been envisioned as the product of pollinator‐mediated selection. Available evidence suggests that major changes in pollinator assemblages may affect the joint expression of floral traits and their phenotypic integration. We expect species within a monophyletic lineage sharing the same pollinator type will show not only similarity in trait means but also similar phenotypic variance‐covariance structures. Here, we tested this expectation using eighteen Salvia species pollinated either by bees or by hummingbirds. Our findings indicated a nonsignificant multivariate phylogenetic signal and a decoupling between means and variance‐covariance phenotypic matrices of floral traits during the evolution to hummingbird pollination. Mean trait value analyses revealed significant differences between bee‐ and hummingbird‐pollinated Salvia species although fewer differences were detected in the covariance structure between groups. Variance‐covariance matrices were much more similar among bee‐ than hummingbird‐pollinated species. This pattern is consistent with the expectation that, unlike hummingbirds, bees physically manipulate the flower, presumably exerting stronger selection pressures favouring morphological convergence among species. Overall, we conclude that the evolution of hummingbird pollination proceeded through different independent transitions. Thus, although the evolution of hummingbird pollination led to a new phenotypic optimum, the process involved the diversification of the covariance structure.     

Trait diversity, heritability and genetic advance in selected germplasm lines of tef

The genetic improvement of the Ethiopian cereal, tef, Eragrostis tef (Zucc.) Trotter, depends upon the variability in the indigenous germplasm. A bi-replicated randomized complete block field experiment was, therefore, carried out at Debre Zeit and Alem Tena in Ethiopia during the 1996 main season to study the pheno-morphic and agronomic trait diversity in 320 tef germplasm lines. All of the 17 traits assessed showed substantial (p < or = 0.001) variation among the lines. Genotypes and locations interacted significantly (p < or = 0.05) on 11 of the traits. At about 50% similarity level, the tef lines grouped into six major clusters of nine to 243 lines. Five principal components (PCs) extracted about 71% of the entire variation of the lines. About 28% of the total variance explained by the first PC was due chiefly to variation in main shoot culm length, diameters of the two basal culm internodes, panicle length and grain yield/panicle. About 16% of the whole variance explained by the second PC originated mainly from variation in the length of the first and second basal culm internodes, grain yield/plant, and peduncle length. The third PC accounting for about 12% of the entire variance resulted largely from variation in harvest index and shoot phytomass yield/plant. Across traits, the phenotypic and genotypic coefficients of variation varied in that order from about 2% for grain yield/panicle to 58% for number of fertile tillers/plant, and from less than 1% for diameters of the two basal culm internodes and grain yield/panicle to 35% for panicle length. Estimates of broad sense heritability and genetic advance (as ratio of the mean) were highest for panicle length (71%) and number of fertile tillers/plant (21%), respectively. But both of these were lowest for the second basal culm internode diameter (< 1%). Overall, the study confirmed that tef is a highly versatile crop with broad trait diversity in the germplasm, and this offers ample opportunities for improvement through breeding.