Independent axes of genetic variation and parallel evolutionary divergence of opercle bone shape in threespine stickleback

Evolution of similar phenotypes in independent populations is often taken as evidence of adaptation to the same fitness optimum. However, the genetic architecture of traits might cause evolution to proceed more often toward particular phenotypes, and less often toward others, independently of the adaptive value of the traits. Freshwater populations of Alaskan threespine stickleback have repeatedly evolved the same distinctive opercle shape after divergence from an oceanic ancestor. Here we demonstrate that this pattern of parallel evolution is widespread, distinguishing oceanic and freshwater populations across the Pacific Coast of North America and Iceland. We test whether this parallel evolution reflects genetic bias by estimating the additive genetic variance-covariance matrix (G) of opercle shape in an Alaskan oceanic (putative ancestral) population. We find significant additive genetic variance for opercle shape and that G has the potential to be biasing, because of the existence of regions of phenotypic space with low additive genetic variation. However, evolution did not occur along major eigenvectors of G, rather it occurred repeatedly in the same directions of high evolvability. We conclude that the parallel opercle evolution is most likely due to selection during adaptation to freshwater habitats, rather than due to biasing effects of opercle genetic architecture.     

Adaptation of rainbow fish to lake and stream habitats

Fish occupy a range of hydrological habitats that exert different demands on locomotor performance. We examined replicate natural populations of the rainbow fishes Melanotaenia eachamensis and M. duboulayi to determine if colonization of low-velocity (lake) habitats by fish from high-velocity (stream) habitats resulted in adaptation of locomotor morphology and performance. Relative to stream conspecifics, lake fish had more posteriorly positioned first dorsal and pelvic fins, and shorter second dorsal fin bases. Habitat dimorphism observed between wild-caught fish was determined to be heritable as it was retained in M. eachamensis offspring raised in a common garden. Repeated evolution of the same heritable phenotype in independently derived populations indicated body shape divergence was a consequence of natural selection. Morphological divergence between hydrological habitats did not support a priori expectations of deeper bodies and caudal peduncles in lake fish. However, observed divergence in fin positioning was consistent with a family-wide association between habitat and morphology, and with empirical studies on other fish species. As predicted, decreased demand for sustained swimming in lakes resulted in a reduction in caudal red muscle area of lake fish relative to their stream counterparts. Melanotaenia duboulayi lake fish also had slower sustained swimming speeds (Ucrit) than stream conspecifics. In M. eachamensis, habitat affected Ucrit of males and females differently. Specifically, females exhibited the pattern observed in M. duboulayi (lake fish had faster Ucrit than stream fish), but the opposite association was observed in males (stream males had slower Ucrit than lake males). Stream M. eachamensis also exhibited a reversed pattern of sexual dimorphism in Ucrit (males slower than females) relative to all other groups (males faster than females). We suggest that M. eachamensis males from streams responded to factors other than water velocity. Although replication of muscle and Ucrit phenotypes across same habitat populations within and/or among species was suggestive of adaptation, the common garden experiment did not confirm a genetic basis to these associations. Kinematic studies should consider the effect of the position and base length of dorsal fins.     

Characterizing the evolution of genetic variance using genetic covariance tensors

Determining how genetic variance changes under selection in natural populations has proved to be a very resilient problem in evolutionary genetics. In the same way that understanding the availability of genetic variance within populations requires the simultaneous consideration of genetic variance in sets of functionally related traits, determining how genetic variance changes under selection in natural populations will require ascertaining how genetic variance–covariance (G) matrices evolve. Here, we develop a geometric framework using higher-order tensors, which enables the empirical characterization of how G matrices have diverged among populations. We then show how divergence among populations in genetic covariance structure can then be associated with divergence in selection acting on those traits using key equations from evolutionary theory. Using estimates of G matrices of eight male sexually selected traits from nine geographical populations of Drosophila serrata, we show that much of the divergence in genetic variance occurred in a single trait combination, a conclusion that could not have been reached by examining variation among the individual elements of the nine G matrices. Divergence in G was primarily in the direction of the major axes of genetic variance within populations, suggesting that genetic drift may be a major cause of divergence in genetic variance among these populations.     

Evidence for morphological and adaptive genetic divergence between lake and stream habitats in European minnows (Phoxinus phoxinus, Cyprinidae)

Natural selection drives local adaptation, potentially even at small temporal and spatial scales. As a result, adaptive genetic and phenotypic divergence can occur among populations living in different habitats. We investigated patterns of differentiation between contrasting lake and stream habitats in the cyprinid fish European minnow (Phoxinus phoxinus) at both the morphological and genomic levels using geometric morphometrics and AFLP markers, respectively. We also used a spatial correlative approach to identify AFLP loci associated with environmental variables representing potential selective forces responsible for adaptation to divergent habitats. Our results identified different morphologies between lakes and streams, with lake fish presenting a deeper body and caudal peduncle compared to stream fish. Body shape variation conformed to a priori predictions concerning biomechanics and swimming performance in lakes vs. streams. Moreover, morphological differentiation was found to be associated with several environmental variables, which could impose selection on body and caudal peduncle shape. We found adaptive genetic divergence between these contrasting habitats in the form of 'outlier' loci (2.9%) whose genetic divergence exceeded neutral expectations. We also detected additional loci (6.6%) not associated with habitat type (lake vs. stream), but contributing to genetic divergence between populations. Specific environmental variables related to trophic dynamics, landscape topography and geography were associated with several neutral and outlier loci. These results provide new insights into the morphological divergence and genetic basis of adaptation to differentiated habitats.     

Genetics of body shape and armour variation in threespine sticklebacks

Patterns of genetic variation and covariation can influence the rate and direction of phenotypic evolution. We explored the possibility that the parallel morphological evolution seen in threespine stickleback (Gasterosteus aculeatus) populations colonizing freshwater environments is facilitated by patterns of genetic variation and covariation in the ancestral (marine) population. We estimated the genetic (G) and phenotypic (P) covariance matrices and directions of maximum additive genetic (g(max) ) and phenotypic (p(max) ) covariances of body shape and armour traits. Our results suggest a role for the ancestral G in explaining parallel morphological evolution in freshwater populations. We also found evidence of genetic constraints owing to the lack of variance in the ancestral G. Furthermore, strong genetic covariances and correlations among traits revealed that selective factors responsible for threespine stickleback body shape and armour divergence may be difficult to disentangle. The directions of g(max) and p(max) were correlated, but the correlations were not high enough to imply that phenotypic patterns of trait variation and covariation within populations are very informative of underlying genetic patterns.     

Evolution of the G-matrix in life history traits in the common frog during a recent colonisation of an island system

Studies of genetic correlations between traits that ostensibly channel the path of evolution away from the direction of natural selection require information on key aspects such as ancestral phenotypes, the duration of adaptive evolution, the direction of natural selection, and genetic covariances. In this study we provide such information in a frog population system. We studied adaptation in life history traits to pool drying in frog populations on islands of known age, which have been colonized from a mainland population. The island populations show strong local adaptation in development time and size. We found that the first eigenvector of the variance–covariance matrix (g _max) had changed between ancestral mainland populations and newly established island populations. Interestingly, there was no divergence in g _max among island populations that differed in their local adaptation in development time and size. Thus, a major change in the genetic covariance of life-history traits occurred in the colonization of the island system, but subsequent local adaptation in development time took place despite the constraints imposed by the genetic covariance structure.     

ADAPTIVE RADIATION ALONG GENETIC LINES OF LEAST RESISTANCE

Are measurements of quantitative genetic variation useful for predicting long-term adaptive evolution? To answer this question, I focus on g_max, the multivariate direction of greatest additive genetic variance within populations. Original data on threespine sticklebacks, together with published genetic measurements from other vertebrates, show that morphological differentiation between species has been biased in the direction of g_max for at least four million years, despite evidence that natural selection is the cause of differentiation. This bias toward the direction of evolution tends to decay with time. Rate of morphological divergence between species is inversely proportional to θ, the angle between the direction of divergence and the direction of greatest genetic variation. The direction of greatest phenotypic variance is not identical with g_max, but for these data is nearly as successful at predicting the direction of species divergence. I interpret the findings to mean that genetic variances and covariances constrain adaptive change in quantitative traits for reasonably long spans of time. An alternative hypothesis, however, cannot be ruled out: that morphological differentiation is biased in the direction g_max because divergence and g_max are both shaped by the same natural selection pressures. Either way, the results reveal that adaptive differentiation occurs principally along “genetic lines of least resistance.”     

Genetic covariances promote climatic adaptation in Australian Drosophila*

Evolutionary potential for adaptation hinges upon the orientation of genetic variation for traits under selection, captured by the additive genetic variance-covariance matrix (G), as well as the evolutionary stability of G. Yet studies that assess both the stability of G and its alignment with selection are extraordinarily rare. We evaluated the stability of G in three Drosophila melanogaster populations that have adapted to local climatic conditions along a latitudinal cline. We estimated population- and sex-specific G matrices for wing size and three climatic stress-resistance traits that diverge adaptively along the cline. To determine how G affects evolutionary potential within these populations, we used simulations to quantify how well G aligns with the direction of trait divergence along the cline (as a proxy for the direction of local selection) and how genetic covariances between traits and sexes influence this alignment. We found that G was stable across the cline, showing no significant divergence overall, or in sex-specific subcomponents, among populations. G also aligned well with the direction of clinal divergence, with genetic covariances strongly elevating evolutionary potential for adaptation to climatic extremes. These results suggest that genetic covariances between both traits and sexes should significantly boost evolutionary responses to environmental change.     

Quantitative genetic variation for thermal performance curves within and among natural populations of Drosophila serrata

Thermal performance curves (TPCs) provide a powerful framework for studying the evolution of continuous reaction norms and for testing hypotheses of thermal adaptation. Although featured heavily in comparative studies, the framework has been comparatively underutilized for quantitative genetic tests of thermal adaptation. We assayed the distribution of genetic (co)variance for TPC (locomotor activity) within and among three natural populations of Drosophila serrata and performed replicated tests of two hypotheses of thermal adaptation--that 'hotter is better' and that a generalist-specialist trade-off underpins the evolution of thermal sensitivity. We detected significant genetic variance within, and divergence among, populations. The 'hotter is better' hypothesis was not supported as the genetic correlations between optimal temperature (T(opt)) and maximum performance (z(max)) were consistently negative. A pattern of variation consistent with a generalist-specialist trade-off was detected within populations and divergence among populations indicated that performance curves were narrower and had higher optimal temperatures in the warmer, but less variable tropical population.     

The role of geography and ecology in shaping repeated patterns of morphological and genetic differentiation between European minnows (Phoxinus phoxinus) from the Pyrenees and the Alps

Neutral and selective processes can drive repeated patterns of evolution in different groups of populations experiencing similar ecological gradients. In this paper, we used a combination of nuclear and mitochondrial DNA markers, as well as geometric morphometrics, to investigate repeated patterns of morphological and genetic divergence of European minnows in two mountain ranges: the Pyrenees and the Alps. European minnows (Phoxinus phoxinus) are cyprinid fish inhabiting most freshwater bodies in Europe, including those in different mountain ranges that could act as major geographical barriers to gene flow. We explored patterns of P. phoxinus phenotypic and genetic diversification along a gradient of altitude common to the two mountain ranges, and tested for isolation by distance (IBD), isolation by environment (IBE) and isolation by adaptation (IBA). The results indicated that populations from the Pyrenees and the Alps belong to two well differentiated, reciprocally monophyletic mtDNA lineages. Substantial genetic differentiation due to geographical isolation within and between populations from the Pyrenees and the Alps was also found using rapidly evolving AFLPs markers (isolation by distance or IBD), as well as morphological differences between mountain ranges. Also, morphology varied strongly with elevation and so did genetic differentiation to a lower extent. Despite moderate evidence for IBE and IBA, and therefore of repeated evolution, substantial population heterogeneity was found at the genetic level, suggesting that selection and population specific genetic drift act in concert to affect genetic divergence.     

Limited genetic differentiation between acoustically divergent populations of urban and rural silvereyes (Zosterops lateralis)

The bioacoustic attributes of vocalisations made by birds in urban environments often differ markedly from those of rural conspecifics. Whether such differences are result from genetic divergence between urban and rural populations, or from plasticity or cultural evolution of song remains poorly understood. Silvereyes (Zosterops lateralis) show evidence of acoustic adaptation to urban noise, modifying both their songs and calls in cities when compared to rural areas. We investigated whether these differences were associated with corresponding morphological and neutral genetic differences. Across six pairs of geographically separate urban and rural populations, all morphological traits measured were similar. Furthermore, genetic analyses of variation at nine microsatellite loci revealed high levels of genetic connectivity between populations, and similar levels of heterozygosity in both habitat types. Consistent directional shifts in song attributes of city birds across large geographic areas thus do not appear to be accompanied by associated morphological or neutral genetic divergence.     

Sequence evolution and phylogenetic signal in control-region and cytochrome b sequences of rainbow fishes (Melanotaeniidae)

The nucleotide sequences of segments of the cytochrome b gene (351 bp), the tRNA(Pro) gene (49 bp), and the control region (approximately 313 bp) of mitochondrial DNA were obtained from 26 fish representing different populations and species of Melanotaenia and one species of Glossolepis, freshwater rainbow fishes confined to Australia and New Guinea. The purpose was to investigate relative rates and patterns of sequence evolution. Overall levels of divergence were similar for the cytochrome b and tRNA control-region sequences, both ranging from < 1% within subspecies to 15%-19% between genera. However, the patterns of sequence evolution differed. For the cytochrome b gene, transitions consistently exceeded transversions, the bias ranging from 4.2:1 to 2:1, depending on the level of sequence divergence. However, in the control-region sequence, a bias toward transitions (2:1) was observed only in comparisons between very similar sequences, and transversions outnumbered transitions in comparisons of divergent sequences. Graphic comparisons suggested that the control region was saturated for transitions at relatively low levels of sequence divergence but accumulated transversions at a greater rate than did the cytochrome b sequence. These distinct patterns of base substitution are associated with differences in A+T content, which is 70% for the tRNA control- region segment versus 50% for cytochrome b. A test for skewness in the distribution of lengths of random trees indicated that both segments contained phylogenetic signal. Parsimony analyses of the data from the two regions, with or without weighting schemes appropriate to the respective patterns of sequence evolution, identified the same five groupings of sequences, but the relationships among the groups differed. However, in most cases the branches uniting different combinations of groups were poorly supported, and the differences among topologies were insignificant. Considering the observed patterns of base substitution and the results of the phylogenetic analyses, we deduce that both the control region and cytochrome b are appropriate for population genetic studies but that the control region is less effective than cytochrome b for resolving relationships among divergent lineages of rainbow fishes.      

Linking stream ecology with morphological variability in a native freshwater fish from semi‐arid Australia

Environmental variation is a potent force affecting phenotypic expression. While freshwater fishes have provided a compelling example of the link between the environment and phenotypic diversity, few studies have been conducted with arid‐zone fishes, particularly those that occur in geographically isolated regions where species typically inhabit intermittent and ephemeral creeks. We investigated morphological variation of a freshwater fish (the western rainbowfish, Melanotaenia australis) inhabiting creeks in the Pilbara region of northwest Australia to determine whether body shape variation correlated with local environmental characteristics, including water velocity, habitat complexity, predator presence, and food availability. We expected that the geographic isolation of creeks within this arid region would result in habitat‐specific morphological specializations. We used landmark‐based geometric morphometrics to quantify the level of morphological variability in fish captured from 14 locations within three distinct subcatchments of a major river system. Western rainbowfish exhibited a range of morphologies, with variation in body depth accounting for a significant proportion (>42%) of the total variance in shape. Sexual dimorphism was also apparent, with males displaying deeper bodies than females. While the measured local habitat characteristics explained little of the observed morphological variation, fish displayed significant morphological differentiation at the level of the subcatchment. Local adaptation may partly explain the geographic patterns of body shape variation, but fine‐scale genetic studies are required to disentangle the effects of genetic differentiation from environmentally determined phenotypic plasticity in body shape. Developing a better understanding of environment–phenotype relationships in species from arid regions will provide important insights into ecological and evolutionary processes in these unique and understudied habitats.     

GENETIC ARCHITECTURE AND POSTZYGOTIC REPRODUCTIVE ISOLATION: EVOLUTION OF BATESON–DOBZHANSKY–MULLER INCOMPATIBILITIES IN A POLYGENIC MODEL

The Bateson–Dobzhansky–Muller model predicts that postzygotic isolation evolves due to the accumulation of incompatible epistatic interactions, but few studies have quantified the relationship between genetic architecture and patterns of reproductive divergence. We examined how the direction and magnitude of epistatic interactions in a polygenic trait under stabilizing selection influenced the evolution of hybrid incompatibilities. We found that populations evolving independently under stabilizing selection experienced suites of compensatory allelic changes that resulted in genetic divergence between populations despite the maintenance of a stable, high‐fitness phenotype. A small number of loci were then incompatible with multiple alleles in the genetic background of the hybrid and the identity of these incompatibility loci changed over the evolution of the populations. For F₁ hybrids, reduced fitness evolved in a window of intermediate strengths of epistatic interactions, but F₂ and backcross hybrids evolved reduced fitness across weak and moderate strengths of epistasis due to segregation variance. Strong epistatic interactions constrained the allelic divergence of parental populations and prevented the development of reproductive isolation. Because many traits with varying genetic architectures must be under stabilizing selection, our results indicate that polygenetic drift is a plausible hypothesis for the evolution of postzygotic reproductive isolation.     

Natural selection drives patterns of lake-stream divergence in stickleback foraging morphology

To what extent are patterns of biological diversification determined by natural selection? We addressed this question by exploring divergence in foraging morphology of threespine stickleback fish inhabiting lake and stream habitats within eight independent watersheds. We found that lake fish generally displayed more developed gill structures and had more streamlined bodies than did stream fish. Diet analysis revealed that these morphological differences were associated with limnetic vs. benthic foraging modes, and that the extent of morphological divergence within watersheds reflected differences in prey resources utilized by lake and stream fish. We also found that patterns of divergence were unrelated to patterns of phenotypic trait (co)variance within populations (i.e. the ‘line of least resistance’). Instead, phenotypic (co)variances were more likely to have been shaped by adaptation to lake vs. stream habitats. Our study thus implicates natural selection as a strong deterministic force driving morphological diversification in lake–stream stickleback. The strength of this inference was obtained by complementing a standard analysis of parallel divergence in means between discrete habitat categories (lake vs. stream) with quantitative estimates of selective forces and information on trait (co)variances.     

Adaptive life-history evolution in the livebearing fish Brachyrhaphis rhabdophora: genetic basis for parallel divergence in age and size at maturity and a test of predator-induced plasticity

I document a genetic basis for parallel evolution of life-history phenotypes in the livebearing fish Brachyrhaphis rhabdophora from northwestern Costa Rica. In previous work, I showed that populations of B. rhabdophora that co-occur with predators attain maturity at smaller sizes than populations that live in predator-free environments. I also demonstrated that this pattern of phenotypic divergence in life histories was independently repeated in at least five isolated drainages. However, life-history phenotypes measured from wild-caught fish could be attributed to environmental effects rather than to genetic differences among populations. In the present study, I reared male fish from four populations (two that co-occur with predators and two from predator-free environments) under four sets of environmental conditions. The pattern of phenotypic divergence in maturation size documented in the field between populations collected from different predation environments persisted after two generations in the laboratory. I also found a genetic basis for differences between populations in the age at which males attain maturity and in growth rates. By rearing fish in four different common environments, I tested for phenotypic plasticity in male life-history traits in response to nonlethal exposure to predators. There was a significant delay in the onset of sexual maturity in fish exposed to predators relative to those in the control, but no differences among treatments in size at maturity or growth rates. These results, coupled with previous work on B. rhabdophora, demonstrate a repeated pattern of parallel evolutionary divergence among genetically isolated populations that is strongly associated with predation.     

Multiple speciation events in an arthropod with divergent evolution in sexual morphology

Sexual selection can facilitate divergent evolution of traits related to mating and consequently promote speciation. Theoretically, independent operation of sexual selection in different populations can lead to divergence of sexual traits among populations and result in allopatric speciation. Here, we show that divergent evolution in sexual morphology affecting mating compatibility (body size and genital morphologies) and speciation have occurred in a lineage of millipedes, the Parafontaria tonominea species complex. In this millipede group, male and female body and genital sizes exhibit marked, correlated divergence among populations, and the diverged morphologies result in mechanical reproductive isolation between sympatric species. The morphological divergence occurred among populations independently and without any correlation with climatic variables, although matching between sexes has been maintained, suggesting that morphological divergence was not a by-product of climatic adaptation. The diverged populations underwent restricted dispersal and secondary contact without hybridization. The extent of morphological difference between sympatric species is variable, as is diversity among allopatric populations; consequently, the species complex appears to contain many species. This millipede case suggests that sexual selection does contribute to species richness via morphological diversification when a lineage of organisms consists of highly divided populations owing to limited dispersal.     

The effects of diadromy and its loss on genomic divergence: The case of amphidromous Galaxias maculatus populations

Understanding the evolutionary mechanisms that affect the genetic divergence between diadromous and resident populations across heterogeneous environments is a challenging task. While diadromy may promote gene flow leading to a lack of genetic differentiation among populations, resident populations tend to be affected by local adaptation and/or plasticity. Studies on these effects on genomic divergence in nonmodel amphidromous species are scarce. Galaxias maculatus, one of the most widespread fish species in the Southern Hemisphere, exhibits two life histories, an ancestral diadromous, specifically, amphidromous form, and a derived freshwater resident form. We examined the genetic diversity and divergence among 20 estuarine and resident populations across the Chilean distribution of G. maculatus and assessed the extent to which selection is involved in the differentiation among resident populations. We obtained nearly 4,400 SNP markers using a RADcap approach for 224 individuals. As expected, collections from estuarine locations typically consist of diadromous individuals. Diadromous populations are highly differentiated from their resident counterparts by both neutral and putative adaptive markers. While diadromous populations exhibit high gene flow and lack site fidelity, resident populations appear to be the product of different colonization events with relatively low genetic diversity and varying levels of gene flow. In particular, the northernmost resident populations were clearly genetically distinct and reproductively isolated from each other suggesting local adaptation. Our study provides insights into the role of life history differences in the maintenance of genetic diversity and the importance of genetic divergence in species evolution.