Parallel evolution of Pitx1 underlies pelvic reduction in Scottish threespine stickleback (Gasterosteus aculeatus)

Inferring the evolution of reproductive isolation in a lineage of fossil threespine stickleback, Gasterosteus doryssus

Darwin attributed the absence of species transitions in the fossil record to his hypothesis that speciation occurs within isolated habitat patches too geographically restricted to be captured by fossil sequences. Mayr's peripatric speciation model added that such speciation would be rapid, further explaining missing evidence of diversification. Indeed, Eldredge and Gould's original punctuated equilibrium model combined Darwin's conjecture, Mayr's model and 124 years of unsuccessfully sampling the fossil record for transitions. Observing such divergence, however, could illustrate the tempo and mode of evolution during early speciation. Here, we investigate peripatric divergence in a Miocene stickleback fish, Gasterosteus doryssus. This lineage appeared and, over approximately 8000 generations, evolved significant reduction of 12 of 16 traits related to armour, swimming and diet, relative to its ancestral population. This was greater morphological divergence than we observed between reproductively isolated, benthic-limnetic ecotypes of extant Gasterosteus aculeatus. Therefore, we infer that reproductive isolation was evolving. However, local extinction of G. doryssus lineages shows how young, isolated, speciating populations often disappear, supporting Darwin's explanation for missing evidence and revealing a mechanism behind morphological stasis. Extinction may also account for limited sustained divergence within the stickleback species complex and help reconcile speciation rate variation observed across time scales.

Dissecting the genetic architecture of quantitative traits using genome-wide identity-by-descent sharing

Additive and dominance genetic variances underlying the expression of quantitative traits are important quantities for predicting short-term responses to selection, but they are notoriously challenging to estimate in most non-model wild populations. Specifically, large-sized or panmictic populations may be characterized by low variance in genetic relatedness among individuals which, in turn, can prevent accurate estimation of quantitative genetic parameters. We used estimates of genome-wide identity-by-descent (IBD) sharing from autosomal SNP loci to estimate quantitative genetic parameters for ecologically important traits in nine-spined sticklebacks (Pungitius pungitius) from a large, outbred population. Using empirical and simulated datasets, with varying sample sizes and pedigree complexity, we assessed the performance of different crossing schemes in estimating additive genetic variance and heritability for all traits. We found that low variance in relatedness characteristic of wild outbred populations with high migration rate can impair the estimation of quantitative genetic parameters and bias heritability estimates downwards. On the other hand, the use of a half-sib/full-sib design allowed precise estimation of genetic variance components and revealed significant additive variance and heritability for all measured traits, with negligible dominance contributions. Genome-partitioning and QTL mapping analyses revealed that most traits had a polygenic basis and were controlled by genes at multiple chromosomes. Furthermore, different QTL contributed to variation in the same traits in different populations suggesting heterogeneous underpinnings of parallel evolution at the phenotypic level. Our results provide important guidelines for future studies aimed at estimating adaptive potential in the wild, particularly for those conducted in outbred large-sized populations.

Morphological description of spontaneous pelvic fin loss in a neotropical cichlid fish

Pelvic fins are a characteristic structure of the vertebrate Bauplan. Yet, pelvic fin loss has occurred repeatedly across a wide diversity of other lineages of tetrapods and at least 48 times in teleost fishes. This pelvic finless condition is often associated with other morphological features such as body elongation, loss of additional structures, and bilateral asymmetry. However, despite the remarkable diversity in the several thousand cichlid fish species, none of them are characterized by the complete absence of pelvic fins. Here, we examined the musculoskeletal structure and associated bilateral asymmetry in Midas cichlids (Amphilophus cf. citrinellus) that lost their pelvic fins spontaneously in the laboratory. Due to this apparent mutational loss of the pelvic girdle and fins, the external and internal anatomy are described in a series of "normal" Midas individuals and their pelvic finless sibling tankmates. First, other traits associated with teleost pelvic fin loss, the genetic basis of pelvic fin loss, and the potential for pleiotropic effects of these genes on other traits in teleosts were all reviewed. Using these traits as a guide, we investigated whether other morphological differences were associated with the pelvic girdle/fin loss. The mean values of the masses of muscle of the pectoral fin, fin ray numbers in the unpaired fins, and oral jaw tooth numbers did not differ between the two pelvic fin morphotypes. However, significant differences in meristic values of the paired traits assessed were observed for the same side of the body between morphotypes. Notably, bilateral asymmetry was found exclusively for the posterior lateral line scales. Finally, we found limited evidence of pleiotropic effects, such as lateral line scale numbers and fluctuating asymmetry between the Midas pelvic fin morphotypes. The fast and relatively isolated changes in the Midas cichlids suggest minor but interesting pleiotropic effects could accompany loss of cichlid pelvic fins.

Spatial heterogeneity in pH, body size and habitat size generates ecological opportunity in an evolutionary radiation

Much of the biological diversity we see today is thought to be the product of evolutionary radiation, the rapid proliferation of species from a single ancestor into multiple discrete forms. Spatial heterogeneity in environmental variables has been proposed as creating the necessary ecological opportunity to stimulate evolutionary radiation. Nonetheless, the ecological mechanisms generating and maintaining diversity in spatially heterogeneous environments are not fully understood. The authors investigated the role of strong spatial heterogeneity in generating ecological opportunity in an evolutionary radiation of freshwater populations of the three-spined stickleback (Gasterosteus aculeatus L.) on the island of North Uist using a spatially explicit Bayesian model. The authors identified pH, loch surface area and body size as predictors of variance in the number of lateral plates that comprise anti-predator armour in G. aculeatus. An East-West gradient of pH, a product of the distinctive environment of North Uist, generates a robust selective environment facilitating G. aculeatus evolutionary radiation. Larger lochs were associated with atypical phenotypes, possibly related to larger population sizes and greater selection efficiency. An association between pH and lateral plate number is likely an effect of body size, with a positive relationship between body size and lateral plate number that is mediated by swimming efficiency in G. aculeatus.

Attempting genetic inference from directional asymmetry during convergent hindlimb reduction in squamates

Loss and reduction in paired appendages are common in vertebrate evolution. How often does such convergent evolution depend on similar developmental and genetic pathways? For example, many populations of the threespine stickleback and ninespine stickleback (Gasterosteidae) have independently evolved pelvic reduction, usually based on independent mutations that caused reduced Pitx1 expression. Reduced Pitx1 expression has also been implicated in pelvic reduction in manatees. Thus, hindlimb reduction stemming from reduced Pitx1 expression has arisen independently in groups that diverged tens to hundreds of millions of years ago, suggesting a potential for repeated use of Pitx1 across vertebrates. Notably, hindlimb reduction based on the reduction in Pitx1 expression produces left-larger directional asymmetry in the vestiges. We used this phenotypic signature as a genetic proxy, testing for hindlimb directional asymmetry in six genera of squamate reptiles that independently evolved hindlimb reduction and for which genetic and developmental tools are not yet developed: Agamodon anguliceps, Bachia intermedia, Chalcides sepsoides, Indotyphlops braminus, Ophisaurus attenuatuas and O. ventralis, and Teius teyou. Significant asymmetry occurred in one taxon, Chalcides sepsoides, whose left-side pelvis and femur vestiges were 18% and 64% larger than right-side vestiges, respectively, suggesting modification in Pitx1 expression in that species. However, there was either right-larger asymmetry or no directional asymmetry in the other five taxa, suggesting multiple developmental genetic pathways to hindlimb reduction in squamates and the vertebrates more generally.

The genetic basis of a novel reproductive strategy in Sulawesi ricefishes: How modularity and a low number of loci shape pelvic brooding

The evolution of complex phenotypes like reproductive strategies is challenging to understand, as they often depend on multiple adaptations that only jointly result in a specific functionality. Sulawesi ricefishes (Adrianichthyidae) evolved a reproductive strategy termed as pelvic brooding. In contrast to the more common transfer brooding, female pelvic brooders carry an egg bundle connected to their body for weeks until the fry hatches. To examine the genetic architecture of pelvic brooding, we crossed the pelvic brooding Oryzias eversi and the transfer brooding Oryzias nigrimas (species divergence time: ∼3.6 my). We hypothesize, that a low number of loci and modularity have facilitated the rapid evolution of pelvic brooding. Traits associated to pelvic brooding, like rib length, pelvic fin length, and morphology of the genital papilla, were correlated in the parental species but correlations were reduced or lost in their F1 and F2 hybrids. Using the Castle-Wright estimator, we found that generally few loci underlie the studied traits. Further, both parental species showed modularity in their body plans. In conclusion, morphological traits related to pelvic brooding were based on a few loci and the mid-body region likely could evolve independently from the remaining body parts. Both factors presumably facilitated the evolution of pelvic brooding.

(Non)Parallel developmental mechanisms in vertebrate appendage reduction and loss

Appendages have been reduced or lost hundreds of times during vertebrate evolution. This phenotypic convergence may be underlain by shared or different molecular mechanisms in distantly related vertebrate clades. To investigate, we reviewed the developmental and evolutionary literature of appendage reduction and loss in more than a dozen vertebrate genera from fish to mammals. We found that appendage reduction and loss was nearly always driven by modified gene expression as opposed to changes in coding sequences. Moreover, expression of the same genes was repeatedly modified across vertebrate taxa. However, the specific mechanisms by which expression was modified were rarely shared. The multiple routes to appendage reduction and loss suggest that adaptive loss of function phenotypes might arise routinely through changes in expression of key developmental genes.

Molecular signatures of selection on the human GLI3 associated central nervous system specific enhancers

The zinc finger-containing transcription factor Gli3 is a key mediator of Hedgehog (Hh) signaling pathway. In vertebrates, Gli3 has widespread expression pattern during early embryonic development. Along the anteroposterior axes of the central nervous system (CNS), dorsoventral neural pattern elaboration is achieved through Hh mediated spatio-temporal deployment of Gli3 transcripts. Previously, we and others uncovered a set of enhancers that mediate many of the known aspects of Gli3 expression during neurogenesis. However, the potential role of Gli3 associated enhancers in trait evolution has not yet received any significant attention. Here, we investigate the evolutionary patterns of Gli3 associated CNS-specific enhancers that have been reported so far. A subset of these enhancers has undergone an accelerated rate of molecular evolution in the human lineage in comparison to other primates/mammals. These fast-evolving enhancers have acquired human-specific changes in transcription factor binding sites (TFBSs). These human-unique changes within subset of Gli3 associated CNS-specific enhancers were further validated as single nucleotide polymorphisms through 1000 Genome Project Phase 3 data. This work not only infers the molecular evolutionary patterns of Gli3 associated enhancers but also provides clues for putative genetic basis of the population-specificity of gene expression regulation.

Genetic population structure constrains local adaptation in sticklebacks

Repeated and independent adaptation to specific environmental conditions from standing genetic variation is common. However, if genetic variation is limited, the evolution of similar locally adapted traits may be restricted to genetically different and potentially less optimal solutions or prevented from happening altogether. Using a quantitative trait locus (QTL) mapping approach, we identified the genomic regions responsible for the repeated pelvic reduction (PR) in three crosses between nine-spined stickleback populations expressing full and reduced pelvic structures. In one cross, PR mapped to linkage group 7 (LG7) containing the gene Pitx1, known to control pelvic reduction also in the three-spined stickleback. In the two other crosses, PR was polygenic and attributed to 10 novel QTL, of which 90% were unique to specific crosses. When screening the genomes from 27 different populations for deletions in the Pitx1 regulatory element, these were only found in the population in which PR mapped to LG7, even though the morphological data indicated large-effect QTL for PR in several other populations as well. Consistent with the available theory and simulations parameterized on empirical data, we hypothesize that the observed variability in genetic architecture of PR is due to heterogeneity in the spatial distribution of standing genetic variation caused by >2× stronger population structuring among freshwater populations and >10× stronger genetic isolation by distance in the sea in nine-spined sticklebacks as compared to three-spined sticklebacks.

Persistence of pelvic spine polymorphism in a panmictic population of brook stickleback (Culaeainconstans) in Alberta, Canada

In the stickleback family (Gasterosteidae), divergence in pelvic morphology is typically associated with reproductive isolation and ecological niche specialization. In this study, we assessed whether reproductive isolation may be involved in the persistence of within-lake polymorphism in brook stickleback (Culaea inconstans (Kirtland, 1840)). We studied populations of brook stickleback in central Alberta, Canada, that are polymorphic for pelvic reduction (i.e., presence–absence of pelvic spines). We investigated the persistence of the pelvic polymorphisms over 40 generations in five lakes. We found that, of the five populations examined, the frequencies of pelvic phenotypes have remained unchanged in three populations, the polymorphism has persisted with significantly changed frequencies of pelvic phenotypes in one population, and the polymorphism has been lost in one population. Using 10 microsatellite loci, we observed low genetic divergence between individuals with and without pelvic spines in one lake, suggesting a lack of assortative mating among individuals with divergent pelvic phenotypes. This observation suggests that the pelvic polymorphism in brook stickleback persists due to balancing selection in the absence of reproductive isolation, although the selective forces were not evaluated. The significant changes in phenotypic frequencies in two populations are likely associated with environmental disturbances such as species introductions.

Inferred genetic architecture underlying evolution in a fossil stickleback lineage

Inferring the genetic architecture of evolution in the fossil record is difficult because genetic crosses are impossible, the acquisition of DNA is usually impossible and phenotype-genotype maps are rarely obvious. However, such inference is valuable because it reveals the genetic basis of microevolutionary change across many more generations than is possible in studies of extant taxa, thereby integrating microevolutionary process and macroevolutionary pattern. Here, we infer the genetic basis of pelvic skeleton reduction in Gasterosteus doryssus, a Miocene stickleback fish from a finely resolved stratigraphic sequence that spans nearly 17,000 years. Reduction in pelvic score, a categorical measure of pelvic structure, resulted primarily from reciprocal frequency changes of two discrete phenotypic classes. Pelvic vestiges also showed left-side larger asymmetry. These patterns implicate Pitx1, a large-effect gene whose deletion generates left-side larger asymmetry of pelvic vestiges in extant, closely related Gasterosteus aculeatus. In contrast, reductions in the length of the pelvic girdle and pelvic spines resulted from directional shifts of unimodal, continuous trait distributions, suggesting an additional suite of genes with minor, additive pelvic effects, again like G. aculeatus. Similar genetic architectures explain shared but phyletically independent patterns across 10 million years of stickleback evolution.

A novel enhancer near the Pitx1 gene influences development and evolution of pelvic appendages in vertebrates

Vertebrate pelvic reduction is a classic example of repeated evolution. Recurrent loss of pelvic appendages in sticklebacks has previously been linked to natural mutations in a pelvic enhancer that maps upstream of Pitx1. The sequence of this upstream PelA enhancer is not conserved to mammals, so we have surveyed a large region surrounding the mouse Pitx1 gene for other possible hind limb control sequences. Here we identify a new pelvic enhancer, PelB, that maps downstream rather than upstream of Pitx1. PelB drives expression in the posterior portion of the developing hind limb, and deleting the sequence from mice alters the size of several hind limb structures. PelB sequences are broadly conserved from fish to mammals. A wild stickleback population lacking the pelvis has an insertion/deletion mutation that disrupts the structure and function of PelB, suggesting that changes in this ancient enhancer contribute to evolutionary modification of pelvic appendages in nature.

(Non)Parallel Evolution

Parallel evolution across replicate populations has provided evolutionary biologists with iconic examples of adaptation. When multiple populations colonize seemingly similar habitats, they may evolve similar genes, traits, or functions. Yet, replicated evolution in nature or in the laboratory often yields inconsistent outcomes: Some replicate populations evolve along highly similar trajectories, whereas other replicate populations evolve to different extents or in distinct directions. To understand these heterogeneous outcomes, biologists are increasingly treating parallel evolution not as a binary phenomenon but rather as a quantitative continuum ranging from parallel to nonparallel. By measuring replicate populations’ positions along this (non)parallel continuum, we can test hypotheses about evolutionary and ecological factors that influence the extent of repeatable evolution. We review evidence regarding the manifestation of (non)parallel evolution in the laboratory, in natural populations, and in applied contexts such as cancer. We enumerate the many genetic, ecological, and evolutionary processes that contribute to variation in the extent of parallel evolution.

The genetic and molecular architecture of phenotypic diversity in sticklebacks

A major goal of evolutionary biology is to identify the genotypes and phenotypes that underlie adaptation to divergent environments. Stickleback fish, including the threespine stickleback (Gasterosteus aculeatus) and the ninespine stickleback (Pungitius pungitius), have been at the forefront of research to uncover the genetic and molecular architecture that underlies phenotypic diversity and adaptation. A wealth of quantitative trait locus (QTL) mapping studies in sticklebacks have provided insight into long-standing questions about the distribution of effect sizes during adaptation as well as the role of genetic linkage in facilitating adaptation. These QTL mapping studies have also provided a basis for the identification of the genes that underlie phenotypic diversity. These data have revealed that mutations in regulatory elements play an important role in the evolution of phenotypic diversity in sticklebacks. Genetic and molecular studies in sticklebacks have also led to new insights on the genetic basis of repeated evolution and suggest that the same loci are involved about half of the time when the same phenotypes evolve independently. When the same locus is involved, selection on standing variation and repeated mutation of the same genes have both contributed to the evolution of similar phenotypes in independent populations.This article is part of the themed issue 'Evo-devo in the genomics era, and the origins of morphological diversity'.

Developmental timing differences underlie armor loss across threespine stickleback populations

Comparing ontogenetic patterns within a well-described evolutionary context aids in inferring mechanisms of change, including heterochronies or deletion of developmental pathways. Because selection acts on phenotypes throughout ontogeny, any within-taxon developmental variation has implications for evolvability. We compare ontogenetic order and timing of locomotion and defensive traits in three populations of threespine stickleback that have evolutionarily divergent adult forms. This analysis adds to the growing understanding of developmental genetic mechanisms of adaptive change in this evolutionary model species by delineating when chondrogenesis and osteogenesis in two derived populations begin to deviate from the developmental pattern in their immediate ancestors. We found that differences in adult defensive morphologies arise through abolished or delayed initiation of these traits rather than via an overall heterochronic shift, that intra-population ontogenetic variation is increased for some derived traits, and that altered armor developmental timing differentiates the derived populations from each other despite parallels in adult lateral plate armor phenotypes. We found that changes in ossified elements of the pelvic armor are linked to delayed and incomplete development of an early-forming pelvic cartilage, and that this disruption likely presages the variable pelvic vestiges documented in many derived populations.

Pigeonetics takes flight: Evolution, development, and genetics of intraspecific variation

Intensive artificial selection over thousands of years has produced hundreds of varieties of domestic pigeon. As Charles Darwin observed, the morphological differences among breeds can rise to the magnitude of variation typically observed among different species. Nevertheless, different pigeon varieties are interfertile, thereby enabling forward genetic and genomic approaches to identify genes that underlie derived traits. Building on classical genetic studies of pigeon variation, recent molecular investigations find a spectrum of coding and regulatory alleles controlling derived traits, including plumage color, feather growth polarity, and limb identity. Developmental and genetic analyses of pigeons are revealing the molecular basis of variation in a classic example of extreme intraspecific diversity, and have the potential to nominate genes that control variation among other birds and vertebrates in general.

Molecular shifts in limb identity underlie development of feathered feet in two domestic avian species

Birds display remarkable diversity in the distribution and morphology of scales and feathers on their feet, yet the genetic and developmental mechanisms governing this diversity remain unknown. Domestic pigeons have striking variation in foot feathering within a single species, providing a tractable model to investigate the molecular basis of skin appendage differences. We found that feathered feet in pigeons result from a partial transformation from hindlimb to forelimb identity mediated by cis-regulatory changes in the genes encoding the hindlimb-specific transcription factor Pitx1 and forelimb-specific transcription factor Tbx5. We also found that ectopic expression of Tbx5 is associated with foot feathers in chickens, suggesting similar molecular pathways underlie phenotypic convergence between these two species. These results show how changes in expression of regional patterning genes can generate localized changes in organ fate and morphology, and provide viable molecular mechanisms for diversity in hindlimb scale and feather distribution.

Genomics of Rapid Incipient Speciation in Sympatric Threespine Stickleback

Ecological speciation is the process by which reproductively isolated populations emerge as a consequence of divergent natural or ecologically-mediated sexual selection. Most genomic studies of ecological speciation have investigated allopatric populations, making it difficult to infer reproductive isolation. The few studies on sympatric ecotypes have focused on advanced stages of the speciation process after thousands of generations of divergence. As a consequence, we still do not know what genomic signatures of the early onset of ecological speciation look like. Here, we examined genomic differentiation among migratory lake and resident stream ecotypes of threespine stickleback reproducing in sympatry in one stream, and in parapatry in another stream. Importantly, these ecotypes started diverging less than 150 years ago. We obtained 34,756 SNPs with restriction-site associated DNA sequencing and identified genomic islands of differentiation using a Hidden Markov Model approach. Consistent with incipient ecological speciation, we found significant genomic differentiation between ecotypes both in sympatry and parapatry. Of 19 islands of differentiation resisting gene flow in sympatry, all were also differentiated in parapatry and were thus likely driven by divergent selection among habitats. These islands clustered in quantitative trait loci controlling divergent traits among the ecotypes, many of them concentrated in one region with low to intermediate recombination. Our findings suggest that adaptive genomic differentiation at many genetic loci can arise and persist in sympatry at the very early stage of ecotype divergence, and that the genomic architecture of adaptation may facilitate this.

Genetic Architecture of Conspicuous Red Ornaments in Female Threespine Stickleback

Explaining the presence of conspicuous female ornaments that take the form of male-typical traits has been a longstanding challenge in evolutionary biology. Such female ornaments have been proposed to evolve via both adaptive and nonadaptive evolutionary processes. Determining the genetic underpinnings of female ornaments is important for elucidating the mechanisms by which such female traits arise and persist in natural populations, but detailed information about their genetic basis is still scarce. In this study, we investigated the genetic architecture of two ornaments, the orange-red throat and pelvic spine, in the threespine stickleback (Gasterosteus aculeatus). Throat coloration is male-specific in ancestral marine populations but has evolved in females in some derived stream populations, whereas sexual dimorphism in pelvic spine coloration is variable among populations. We find that ornaments share a common genetic architecture between the sexes. At least three independent genomic regions contribute to red throat coloration, and harbor candidate genes related to pigment production and pigment cell differentiation. One of these regions is also associated with spine coloration, indicating that both ornaments might be mediated partly via pleiotropic genetic mechanisms.

Iterative development and the scope for plasticity: contrasts among trait categories in an adaptive radiation

Phenotypic plasticity can influence evolutionary change in a lineage, ranging from facilitation of population persistence in a novel environment to directing the patterns of evolutionary change. As the specific nature of plasticity can impact evolutionary consequences, it is essential to consider how plasticity is manifested if we are to understand the contribution of plasticity to phenotypic evolution. Most morphological traits are developmentally plastic, irreversible, and generally considered to be costly, at least when the resultant phenotype is mis-matched to the environment. At the other extreme, behavioral phenotypes are typically activational (modifiable on very short time scales), and not immediately costly as they are produced by constitutive neural networks. Although patterns of morphological and behavioral plasticity are often compared, patterns of plasticity of life history phenotypes are rarely considered. Here we review patterns of plasticity in these trait categories within and among populations, comprising the adaptive radiation of the threespine stickleback fish Gasterosteus aculeatus. We immediately found it necessary to consider the possibility of iterated development, the concept that behavioral and life history trajectories can be repeatedly reset on activational (usually behavior) or developmental (usually life history) time frames, offering fine tuning of the response to environmental context. Morphology in stickleback is primarily reset only in that developmental trajectories can be altered as environments change over the course of development. As anticipated, the boundaries between the trait categories are not clear and are likely to be linked by shared, underlying physiological and genetic systems.

Using the phenoscape knowledgebase to relate genetic perturbations to phenotypic evolution

The abundance of phenotypic diversity among species can enrich our knowledge of development and genetics beyond the limits of variation that can be observed in model organisms. The Phenoscape Knowledgebase (KB) is designed to enable exploration and discovery of phenotypic variation among species. Because phenotypes in the KB are annotated using standard ontologies, evolutionary phenotypes can be compared with phenotypes from genetic perturbations in model organisms. To illustrate the power of this approach, we review the use of the KB to find taxa showing evolutionary variation similar to that of a query gene. Matches are made between the full set of phenotypes described for a gene and an evolutionary profile, the latter of which is defined as the set of phenotypes that are variable among the daughters of any node on the taxonomic tree. Phenoscape's semantic similarity interface allows the user to assess the statistical significance of each match and flags matches that may only result from differences in annotation coverage between genetic and evolutionary studies. Tools such as this will help meet the challenge of relating the growing volume of genetic knowledge in model organisms to the diversity of phenotypes in nature. The Phenoscape KB is available at http://kb.phenoscape.org.

The diversification of Heliconius butterflies: what have we learned in 150 years?

Research into Heliconius butterflies has made a significant contribution to evolutionary biology. Here, we review our understanding of the diversification of these butterflies, covering recent advances and a vast foundation of earlier work. Whereas no single group of organisms can be sufficient for understanding life's diversity, after years of intensive study, research into Heliconius has addressed a wide variety of evolutionary questions. We first discuss evidence for widespread gene flow between Heliconius species and what this reveals about the nature of species. We then address the evolution and diversity of warning patterns, both as the target of selection and with respect to their underlying genetic basis. The identification of major genes involved in mimetic shifts, and homology at these loci between distantly related taxa, has revealed a surprising predictability in the genetic basis of evolution. In the final sections, we consider the evolution of warning patterns, and Heliconius diversity more generally, within a broader context of ecological and sexual selection. We consider how different traits and modes of selection can interact and influence the evolution of reproductive isolation.

Parallel developmental genetic features underlie stickleback gill raker evolution

Background

Convergent evolution, the repeated evolution of similar phenotypes in independent lineages, provides natural replicates to study mechanisms of evolution. Cases of convergent evolution might have the same underlying developmental and genetic bases, implying that some evolutionary trajectories might be predictable. In a classic example of convergent evolution, most freshwater populations of threespine stickleback fish have independently evolved a reduction of gill raker number to adapt to novel diets. Gill rakers are a segmentally reiterated set of dermal bones important for fish feeding. A previous large quantitative trait locus (QTL) mapping study using a marine × freshwater F2 cross identified QTL on chromosomes 4 and 20 with large effects on evolved gill raker reduction.

Results

By examining skeletal morphology in adult and developing sticklebacks, we find heritable marine/freshwater differences in gill raker number and spacing that are specified early in development. Using the expression of the Ectodysplasin receptor (Edar) gene as a marker of raker primordia, we find that the differences are present before the budding of gill rakers occurs, suggesting an early change to a lateral inhibition process controlling raker primordia spacing. Through linkage mapping in F2 fish from crosses with three independently derived freshwater populations, we find in all three crosses QTL overlapping both previously identified QTL on chromosomes 4 and 20 that control raker number. These two QTL affect the early spacing of gill raker buds.

Conclusions

Collectively, these data demonstrate that parallel developmental genetic features underlie the convergent evolution of gill raker reduction in freshwater sticklebacks, suggesting that even highly polygenic adaptive traits can have a predictable developmental genetic basis.

Modular skeletal evolution in sticklebacks is controlled by additive and clustered quantitative trait Loci

Understanding the genetic architecture of evolutionary change remains a long-standing goal in biology. In vertebrates, skeletal evolution has contributed greatly to adaptation in body form and function in response to changing ecological variables like diet and predation. Here we use genome-wide linkage mapping in threespine stickleback fish to investigate the genetic architecture of evolved changes in many armor and trophic traits. We identify >100 quantitative trait loci (QTL) controlling the pattern of serially repeating skeletal elements, including gill rakers, teeth, branchial bones, jaws, median fin spines, and vertebrae. We use this large collection of QTL to address long-standing questions about the anatomical specificity, genetic dominance, and genomic clustering of loci controlling skeletal differences in evolving populations. We find that most QTL (76%) that influence serially repeating skeletal elements have anatomically regional effects. In addition, most QTL (71%) have at least partially additive effects, regardless of whether the QTL controls evolved loss or gain of skeletal elements. Finally, many QTL with high LOD scores cluster on chromosomes 4, 20, and 21. These results identify a modular system that can control highly specific aspects of skeletal form. Because of the general additivity and genomic clustering of major QTL, concerted changes in both protective armor and trophic traits may occur when sticklebacks inherit either marine or freshwater alleles at linked or possible "supergene" regions of the stickleback genome. Further study of these regions will help identify the molecular basis of both modular and coordinated changes in the vertebrate skeleton.

Identification of major and minor QTL for ecologically important morphological traits in three-spined sticklebacks (Gasterosteus aculeatus)

Quantitative trait locus (QTL) mapping studies of Pacific three-spined sticklebacks (Gasterosteus aculeatus) have uncovered several genomic regions controlling variability in different morphological traits, but QTL studies of Atlantic sticklebacks are lacking. We mapped QTL for 40 morphological traits, including body size, body shape, and body armor, in a F₂ full-sib cross between northern European marine and freshwater three-spined sticklebacks. A total of 52 significant QTL were identified at the 5% genome-wide level. One major QTL explaining 74.4% of the total variance in lateral plate number was detected on LG4, whereas several major QTL for centroid size (a proxy for body size), and the lengths of two dorsal spines, pelvic spine, and pelvic girdle were mapped on LG21 with the explained variance ranging from 27.9% to 57.6%. Major QTL for landmark coordinates defining body shape variation also were identified on LG21, with each explaining ≥15% of variance in body shape. Multiple QTL for different traits mapped on LG21 overlapped each other, implying pleiotropy and/or tight linkage. Thus, apart from providing confirmatory data to support conclusions born out of earlier QTL studies of Pacific sticklebacks, this study also describes several novel QTL of both major and smaller effect for ecologically important traits. The finding that many major QTL mapped on LG21 suggests that this linkage group might be a hotspot for genetic determinants of ecologically important morphological traits in three-spined sticklebacks.

Challenges and prospects in genome-wide quantitative trait loci mapping of standing genetic variation in natural populations

A considerable challenge in evolutionary genetics is to understand the genetic mechanisms that facilitate or impede evolutionary adaptation in natural populations. For this, we must understand the genetic loci contributing to trait variation and the selective forces acting on them. The decreased costs and increased feasibility of obtaining genotypic data on a large number of individuals have greatly facilitated gene mapping in natural populations, particularly because organisms whose genetics have been historically difficult to study are now within reach. Here we review the methods available to evolutionary ecologists interested in dissecting the genetic basis of traits in natural populations. Our focus lies on standing genetic variation in outbred populations. We present an overview of the current state of research in the field, covering studies on both plants and animals. We also draw attention to particular challenges associated with the discovery of quantitative trait loci and discuss parallels to studies on crops, livestock, and humans. Finally, we point to some likely future developments in genetic mapping studies.

QTL analysis of behavior in nine-spined sticklebacks (Pungitius pungitius)

The genetic architecture of behavioral traits is yet relatively poorly understood in most non-model organisms. Using an F2-intercross (n = 283 offspring) between behaviorally divergent nine-spined stickleback (Pungitius pungitius) populations, we tested for and explored the genetic basis of different behavioral traits with the aid of quantitative trait locus (QTL) analyses based on 226 microsatellite markers. The behaviors were analyzed both separately (viz. feeding activity, risk-taking and exploration) and combined in order to map composite behavioral type. Two significant QTL-explaining on average 6 % of the phenotypic variance-were detected for composite behavioral type on the experiment-wide level, located on linkage groups 3 and 8. In addition, several suggestive QTL located on six other linkage groups were detected on the chromosome-wide level. Apart from providing evidence for the genetic basis of behavioral variation, the results provide a good starting point for finer-scale analyses of genetic factors influencing behavioral variation in the nine-spined stickleback.

A perspective on micro-evo-devo: progress and potential

The term "micro-evo-devo" refers to the combined study of the genetic and developmental bases of natural variation in populations and the evolutionary forces that have shaped this variation. It thus represents a synthesis of the fields of evolutionary developmental biology and population genetics. As has been pointed out by several others, this synthesis can provide insights into the evolution of organismal form and function that have not been possible within these individual disciplines separately. Despite a number of important successes in micro-evo-devo, however, it appears that evo devo and population genetics remain largely separate spheres of research, limiting their ability to address evolutionary questions. This also risks pushing contemporary evo devo to the fringes of evolutionary biology because it does not describe the causative molecular changes underlying evolution or the evolutionary forces involved. Here we reemphasize the theoretical and practical importance of micro-evo-devo as a strategy for understanding phenotypic evolution, review the key recent insights that it has provided, and present a perspective on both the potential and the remaining challenges of this exciting interdisciplinary field.

Genetic architecture of parallel pelvic reduction in ninespine sticklebacks

Teleost fish genomes are known to be evolving faster than those of other vertebrate taxa. Thus, fish are suited to address the extent to which the same vs. different genes are responsible for similar phenotypic changes in rapidly evolving genomes of evolutionary independent lineages. To gain insights into the genetic basis and evolutionary processes behind parallel phenotypic changes within and between species, we identified the genomic regions involved in pelvic reduction in Northern European ninespine sticklebacks (Pungitius pungitius) and compared them to those of North American ninespine and threespine sticklebacks (Gasterosteus aculeatus). To this end, we conducted quantitative trait locus (QTL) mapping using 283 F₂ progeny from an interpopulation cross. Phenotypic analyses indicated that pelvic reduction is a recessive trait and is inherited in a simple Mendelian fashion. Significant QTL for pelvic spine and girdle lengths were identified in the region of the Pituitary homeobox transcription factor 1 (Pitx1) gene, also responsible for pelvic reduction in threespine sticklebacks. The fact that no QTL was observed in the region identified in the mapping study of North American ninespine sticklebacks suggests that an alternative QTL for pelvic reduction has emerged in this species within the past 1.6 million years after the split between Northern European and North American populations. In general, our study provides empirical support for the view that alternative genetic mechanisms that lead to similar phenotypes can evolve over short evolutionary time scales.

Ecological causes of morphological evolution in the three-spined stickleback

The central assumption of evolutionary theory is that natural selection drives the adaptation of populations to local environmental conditions, resulting in the evolution of adaptive phenotypes. The three-spined stickleback (Gasterosteus aculeatus) displays remarkable phenotypic variation, offering an unusually tractable model for understanding the ecological mechanisms underpinning adaptive evolutionary change. Using populations on North Uist, Scotland we investigated the role of predation pressure and calcium limitation on the adaptive evolution of stickleback morphology and behavior. Dissolved calcium was a significant predictor of plate and spine morph, while predator abundance was not. Stickleback latency to emerge from a refuge varied with morph, with populations with highly reduced plates and spines and high predation risk less bold. Our findings support strong directional selection in three-spined stickleback evolution, driven by multiple selective agents.

Molecular evolutionary and population genomic analysis of the nine-spined stickleback using a modified restriction-site-associated DNA tag approach

In recent years, the explosion of affordable next generation sequencing technology has provided an unprecedented opportunity to conduct genome-wide studies of adaptive evolution in organisms previously lacking extensive genomic resources. Here, we characterize genome-wide patterns of variability and differentiation using pooled DNA from eight populations of the nine-spined stickleback (Pungitius pungitius L.) from marine, lake and pond environments. We developed a novel genome complexity reduction protocol, defined as paired-end double restriction-site-associated DNA (PE dRAD), to maximize read coverage at sequenced locations. This allowed us to identify over 114,000 short consensus sequences and 15,000 SNPs throughout the genome. A total of 6834 SNPs mapped to a single position on the related three-spined stickleback genome, allowing the detection of genomic regions affected by divergent and balancing selection, both between species and between freshwater and marine populations of the nine-spined stickleback. Gene ontology analysis revealed 15 genomic regions with elevated diversity, enriched for genes involved in functions including immunity, chemical stimulus response, lipid metabolism and signalling pathways. Comparisons of marine and freshwater populations identified nine regions with elevated differentiation related to kidney development, immunity and MAP kinase pathways. In addition, our analysis revealed that a large proportion of the identified SNPs mapping to LG XII is likely to represent alternative alleles from divergent X and Y chromosomes, rather than true autosomal markers following Mendelian segregation. Our work demonstrates how population-wide sequencing and combining inter- and intra-specific RAD analysis can uncover genome-wide patterns of differentiation and adaptations in a non-model species.

Parallel reduction in expression of the eye development gene hedgehog in separately derived cave populations of the amphipod Gammarus minus

Caves provide excellent settings to examine evolutionary questions. Subterranean environments are characterized by similar and consistent conditions. Cave-adapted species often share characteristics such as diminished pigmentation, elongated limbs and reduced or absent eyes. Relatively little is known about the evolution and development of troglomorphic traits in invertebrates. In this study, we compare expression of the eye development genes hedgehog, pax6, sine oculis and dachshund in individuals from multiple independently derived cave populations of the amphipod Gammarus minus. hedgehog expression was significantly reduced in cave populations, compared to genetically related surface populations. Interestingly, no differences were found in pax6, sine oculis or dachshund expression. Because hedgehog-related genes are also involved in eye reduced in Astyanax mexicanus, these genes may be consistent targets of evolution during cave adaptation. These results provide support for the hypothesis of genomic 'hotspots' of evolution and allow comparison of adaptive mechanisms among diverse animals in subterranean environments.

Ecological transcriptomics of lake-type and riverine sockeye salmon (Oncorhynchus nerka)

BACKGROUND

There are a growing number of genomes sequenced with tentative functions assigned to a large proportion of the individual genes. Model organisms in laboratory settings form the basis for the assignment of gene function, and the ecological context of gene function is lacking. This work addresses this shortcoming by investigating expressed genes of sockeye salmon (Oncorhynchus nerka) muscle tissue. We compared morphology and gene expression in natural juvenile sockeye populations related to river and lake habitats. Based on previously documented divergent morphology, feeding strategy, and predation in association with these distinct environments, we expect that burst swimming is favored in riverine population and continuous swimming is favored in lake-type population. In turn we predict that morphology and expressed genes promote burst swimming in riverine sockeye and continuous swimming in lake-type sockeye.

RESULTS

We found the riverine sockeye population had deep, robust bodies and lake-type had shallow, streamlined bodies. Gene expression patterns were measured using a 16 k microarray, discovering 141 genes with significant differential expression. Overall, the identity and function of these genes was consistent with our hypothesis. In addition, Gene Ontology (GO) enrichment analyses with a larger set of differentially expressed genes found the "biosynthesis" category enriched for the riverine population and the "metabolism" category enriched for the lake-type population.

CONCLUSIONS

This study provides a framework for understanding sockeye life history from a transcriptomic perspective and a starting point for more extensive, targeted studies determining the ecological context of genes.

The role of gene expression in ecological speciation

Ecological speciation is the process by which barriers to gene flow between populations evolve due to adaptive divergence via natural selection. A relatively unexplored area in ecological speciation is the role of gene expression. Gene expression may be associated with ecologically important phenotypes not evident from morphology and play a role during colonization of new environments. Here we review two potential roles of gene expression in ecological speciation: (1) its indirect role in facilitating population persistence and (2) its direct role in contributing to genetically based reproductive isolation. We find indirect evidence that gene expression facilitates population persistence, but direct tests are lacking. We also find clear examples of gene expression having effects on phenotypic traits and adaptive genetic divergence, but links to the evolution of reproductive isolation itself remain indirect. Gene expression during adaptive divergence seems to often involve complex genetic architectures controlled by gene networks, regulatory regions, and “eQTL hotspots.” Nonetheless, we review how approaches for isolating the functional mutations contributing to adaptive divergence are proving to be successful. The study of gene expression has promise for increasing our understanding ecological speciation, particularly when integrative approaches are applied.

Pelvic fins in teleosts: structure, function and evolution

The pelvic fins of teleosts are paired appendages that are considered to be homologous to the hind limbs of tetrapods. Because they are less important for swimming, their morphology and function can be flexibly modified, and such modifications have probably facilitated the adaptations of teleosts to various environments. Recently, among these modifications, pelvic-fin loss has gained attention in evolutionary developmental biology. Pelvic-fin loss, however, has only been investigated in a few model species, and various biological aspects of pelvic fins in teleosts in general remain poorly understood. This review summarizes the current state of knowledge regarding pelvic fins, such as their structure, function and evolution, to elucidate their contribution to the considerable diversity of teleosts. This information could be invaluable for future investigations into various aspects of pelvic fins, which will provide clues to understanding the evolution, diversity and adaptations of teleosts.

Perspectives on the genetic architecture of divergence in body shape in sticklebacks

The body shape of fishes encompasses a number of morphological traits that are intrinsically linked to functional systems and affect various measures of performance, including swimming, feeding, and avoiding predators. Changes in shape can allow a species to exploit a new ecological niche and can lead to ecological speciation. Body shape results from the integration of morphological, behavioral and physiological traits. It has been well established that functional interdependency among traits plays a large role in constraining the evolution of shape, affecting both the speed and the repeated evolution of particular body shapes. However, it is less clear what role genetic or developmental constraints might play in biasing the rate or direction of the evolution of body shape. Here, we suggest that the threespine stickleback (Gasterosteus aculeatus) is a powerful model system in which to address the extent to which genetic or developmental constraints play a role in the evolution of body shape in fishes. We review the existing data that begins to address these issues in sticklebacks and provide suggestions for future areas of research that will be particularly fruitful for illuminating the mechanisms that contribute to the evolution of body shape in fishes.

Evolutionary tinkering with conserved components of a transcriptional regulatory network

A surprising level of evolutionary plasticity is revealed by analysis of differences between related yeasts in the mechanisms regulating the essential cellular process of ribosomal gene expression.

Gene expression variation between species is a major contributor to phenotypic diversity, yet the underlying flexibility of transcriptional regulatory networks remains largely unexplored. Transcription of the ribosomal regulon is a critical task for all cells; in S. cerevisiae the transcription factors Rap1, Fhl1, Ifh1, and Hmo1 form a multi-subunit complex that controls ribosomal gene expression, while in C. albicans this regulation is under the control of Tbf1 and Cbf1. Here, we analyzed, using full-genome transcription factor mapping, the roles, in both S. cerevisiae and C. albicans, of each orthologous component of this complete set of regulators. We observe dramatic changes in the binding profiles of the generalist regulators Cbf1, Hmo1, Rap1, and Tbf1, while the Fhl1-Ifh1 dimer is the only component involved in ribosomal regulation in both fungi: it activates ribosomal protein genes and rDNA expression in a Tbf1-dependent manner in C. albicans and a Rap1-dependent manner in S. cerevisiae. We show that the transcriptional regulatory network governing the ribosomal expression program of two related yeast species has been massively reshaped in cis and trans. Changes occurred in transcription factor wiring with cellular functions, movements in transcription factor hierarchies, DNA-binding specificity, and regulatory complexes assembly to promote global changes in the architecture of the fungal transcriptional regulatory network.

Conserved metabolic machineries direct energy production and investment in most life forms. However, variation in the transcriptional regulation of the genes that encode this machinery has been observed and shown to contribute to phenotypic differences between species. Here, we show that the regulatory circuits governing the expression of central metabolic components (in this case the ribosomes) in different yeast species have an unexpected level of evolutionary plasticity. Most transcription factors involved in the regulation of expression of ribosomal genes have in fact been reused in new ways during the evolutionary time separating S. cerevisiae and C. albicans to generate global changes in transcriptional network structures and new ribosomal regulatory complexes.

Phenotypic plasticity facilitates recurrent rapid adaptation to introduced predators

A central role for phenotypic plasticity in adaptive evolution is often posited yet lacks empirical support. Selection for the stable production of an induced phenotype is hypothesized to modify the regulation of preexisting developmental pathways, producing rapid adaptive change. We examined the role of plasticity in rapid adaptation of the zooplankton Daphnia melanica to novel fish predators. Here we show that plastic up-regulation of the arthropod melanin gene dopa decarboxylase (Ddc) in the absence of UV radiation is associated with reduced pigmentation in D. melanica. Daphnia populations coexisting with recently introduced fish exhibit environmentally invariant up-regulation of Ddc, accompanied by constitutive up-regulation of the interacting arthropod melanin gene ebony. Both changes in regulation are associated with adaptive reduction in the plasticity and mean expression of melanin. Our results provide evidence that the developmental mechanism underlying ancestral plasticity in response to an environmental factor has been repeatedly co-opted to facilitate rapid adaptation to an introduced predator.

Adaptive evolution of pelvic reduction in sticklebacks by recurrent deletion of a Pitx1 enhancer

The molecular mechanisms underlying major phenotypic changes that have evolved repeatedly in nature are generally unknown. Pelvic loss in different natural populations of threespine stickleback fish has occurred through regulatory mutations deleting a tissue-specific enhancer of the Pituitary homeobox transcription factor 1 (Pitx1) gene. The high prevalence of deletion mutations at Pitx1 may be influenced by inherent structural features of the locus. Although Pitx1 null mutations are lethal in laboratory animals, Pitx1 regulatory mutations show molecular signatures of positive selection in pelvic-reduced populations. These studies illustrate how major expression and morphological changes can arise from single mutational leaps in natural populations, producing new adaptive alleles via recurrent regulatory alterations in a key developmental control gene.

Implications of a fossil stickleback assemblage for Darwinian gradualism

Darwin postulated that a complete fossil record would contain numerous gradual transitions between ancestral and descendant species, but 150 years after publication of The Origin of Species, few such transitions have materialized. The fossil stickleback Gasterosteus doryssus and the deposit in which it occurs provide excellent conditions to detect such transitions. Abundant, well-preserved fossils occur in a stratigraphic setting with fine temporal resolution. The paleoecology of G. doryssus resembles the ecology of modern lakes that harbour the phenotypically similar three-spined stickleback Gasterosteus aculeatus. Gasterosteus aculeatus are primitively highly armoured, but G. doryssus comprised two contemporaneous biological species with relatively weak armour, including a near-shore, benthic feeder (benthic) and an offshore planktivore (limnetic). The benthic species expanded its range into the limnetic zone of the lake, where it apparently switched to planktivory and evolved reduced armour within c. 5000 years in response to directional selection. Although gradual evolution of mean phenotypes occurred, a single major gene caused much of evolutionary change of the pelvic skeleton. Thus, Darwin's expectation that transitions between species in the fossil record would be gradual was met at a fine time scale, but for pelvic structure, a well-studied trait, his expectation that gradual change would depend entirely on numerous, small, heritable differences among individuals was incorrect.

The Darwinian stickleback Gasterosteus aculeatus: a history of evolutionary studies

The history of studies on the taxonomy and evolutionary biology of the three-spined stickleback Gasterosteus aculeatus from the 18th century to the present is reviewed. After the publication of Darwin's Origin of the Species, four important dates, 1925, 1947, 1967 and 2001, are identified as marking major gains in the understanding of the evolution of the diversity in morphological, life-history, physiological and behavioural traits that characterizes G. aculeatus. The period 1925-1970 led to the identification of the main themes of research: status and adaptive significance of lateral-plate morphs; inter and intrapopulation trait variation in freshwater resident G. aculeatus and the adaptive significance of the variation. Between 1970 and 2001, these themes were investigated using variation observed particularly along the Pacific coast of the U.S.A. and Canada, notably in the Cook Inlet region of Alaska and the Haida Gwaii Archipelago. Studies on adaptive radiation and reproductive isolation in lacustrine, ecomorph pairs (limnetics and benthics) discovered in the Strait of Georgia region have been particularly productive. From 2001, the application of genomic studies to these problems began to open up the study of the relationships between genotype, phenotype and selective advantage to causal analysis.

The genetic architecture of skeletal convergence and sex determination in ninespine sticklebacks

The history of life offers plentiful examples of convergent evolution, the independent derivation of similar phenotypes in distinct lineages. The emergence of convergent phenotypes among closely related lineages (frequently termed "parallel" evolution) is often assumed to result from changes in similar genes or developmental pathways, but the genetic origins of convergence remains poorly understood. Ninespine (Pungitius pungitius) and threespine (Gasterosteus aculeatus) stickleback fish provide many examples of convergent evolution of adaptive phenotypes, both within and between genera. The genetic architecture of several important traits is now known for threespine sticklebacks; thus, ninespine sticklebacks provide a unique opportunity to critically test whether similar or different chromosome regions control similar phenotypes in these lineages. We have generated the first genome-wide linkage map for ninespine sticklebacks and used quantitative trait locus mapping to identify chromosome regions controlling several skeletal traits and sex determination. In ninespine sticklebacks, these traits mapped to chromosome regions not previously known to control the corresponding traits in threespine sticklebacks. Therefore, convergent morphological evolution in these related, but independent, vertebrate lineages might have different genetic origins. Comparative genetics in sticklebacks provides an exciting opportunity to study the mechanisms controlling similar phenotypic changes in different animal groups.

Asymmetric lower-limb malformations in individuals with homeobox PITX1 gene mutation

Clubfoot is one of the most common severe musculoskeletal birth defects, with a worldwide incidence of 1 in 1000 live births. In the present study, we describe a five-generation family with asymmetric right-sided predominant idiopathic clubfoot segregating as an autosomal-dominant condition with incomplete penetrance. Other lower-limb malformations, including patellar hypoplasia, oblique talus, tibial hemimelia, developmental hip dysplasia, and preaxial polydactyly, were also present in some family members. Genome-wide linkage analysis with Affymetrix GeneChip Mapping 10K mapping data from 13 members of this family revealed a multipoint LOD(max) of 3.31 on chromosome 5q31. A single missense mutation (c.388G-->A) was identified in PITX1, a bicoid-related homeodomain transcription factor critical for hindlimb development, and segregated with disease in this family. The PITX1 E130K mutation is located in the highly conserved homeodomain and reduces the ability of PITX1 to transactivate a luciferase reporter. The PITX1 E130K mutation also suppresses wild-type PITX1 activity in a dose-dependent manner, suggesting dominant-negative effects on transcription. The propensity for right-sided involvement in tibial hemimelia and clubfoot suggests that PITX1, or pathways involving PITX1, may be involved in their etiology. Implication of a gene involved in early limb development in clubfoot pathogenesis also suggests additional pathways for future investigations of idiopathic clubfoot etiology in humans.

The genetic basis of phenotypic convergence in beach mice: similar pigment patterns but different genes

Convergent evolution is a widespread phenomenon seen in diverse organisms inhabiting similar selective environments. However, it is unclear if similar phenotypes are produced by the same or different genes and mutations. Here we analyze the molecular mechanisms underlying convergent pigment pattern among subspecies of the beach mouse (Peromyscus polionotus) inhabiting the Gulf and Atlantic coasts of Florida. In these two geographic regions, separated by more than 300 km, "beach mice" have lighter colored coats than do their mainland counterparts, produced by natural selection for camouflage against the pale coastal sand dunes. We measured color pattern in eight beach mouse subspecies and showed that three of the Gulf Coast subspecies are more phenotypically similar to an Atlantic coast subspecies than to their Gulf Coast neighbors. However, light-colored beach mice do not form a monophyletic group. Previous results implicated a single derived amino acid change in the melanocortin-1 receptor (Mc1r) as a major contributor to pigment pattern in the Gulf Coast beach mice; despite phenotypic similarities, the derived Mc1r allele was not found in the Atlantic coast beach mouse populations. Here we show that Atlantic coast beach mice have high levels of Mc1r polymorphism but they lack unique alleles. Functional assays revealed that single amino acid mutations segregating in Atlantic coast beach mice do not cause any change in Mc1r activity compared with the activity of Mc1r from dark-colored mice. These joint results show that convergent pigment patterns in recently diverged beach mouse subspecies--whose developmental constraints are presumably similar--have evolved through a diversity of genetic mechanisms.