Pleiotropic effects of a single gene on skeletal development and sensory system patterning in sticklebacks

Freshwater Colonization, Adaptation, and Genomic Divergence in Threespine Stickleback

The Threespine Stickleback is ancestrally a marine fish, but many marine populations breed in fresh water (i.e., are anadromous), facilitating their colonization of isolated freshwater habitats a few years after they form. Repeated adaptation to fresh water during at least 10 My and continuing today has led to Threespine Stickleback becoming a premier system to study rapid adaptation. Anadromous and freshwater stickleback breed in sympatry and may hybridize, resulting in introgression of freshwater-adaptive alleles into anadromous populations, where they are maintained at low frequencies as ancient standing genetic variation. Anadromous stickleback have accumulated hundreds of freshwater-adaptive alleles that are disbursed as few loci per marine individual and provide the basis for adaptation when they colonize fresh water. Recent whole-lake experiments in lakes around Cook Inlet, Alaska have revealed how astonishingly rapid and repeatable this process is, with the frequency of 40% of the identified freshwater-adaptive alleles increasing from negligible (∼1%) in the marine founder to ≥50% within ten generations in fresh water, and freshwater phenotypes evolving accordingly. These high rates of genomic and phenotypic evolution imply very intense directional selection on phenotypes of heterozygotes. Sexual recombination rapidly assembles freshwater-adaptive alleles that originated in different founders into multilocus freshwater haplotypes, and regions important for adaptation to freshwater have suppressed recombination that keeps advantageous alleles linked within large haploblocks. These large haploblocks are also older and appear to have accumulated linked advantageous mutations. The contemporary evolution of Threespine Stickleback has provided broadly applicable insights into the mechanisms that facilitate rapid adaptation.

Leg length and bristle density, both necessary for water surface locomotion, are genetically correlated in water striders

Access to hitherto unexploited ecological opportunities is associated with phenotypic evolution and often results in significant lineage diversification. Yet our understanding of the mechanisms underlying such adaptive traits remains limited. Water striders have been able to exploit the water-air interface, primarily facilitated by changes in the density of hydrophobic bristles and a significant increase in leg length. These two traits are functionally correlated and are both necessary for generating efficient locomotion on the water surface. Whether bristle density and leg length have any cellular or developmental genetic mechanisms in common is unknown. Here, we combine comparative genomics and transcriptomics with functional RNA interference assays to examine the developmental genetic and cellular mechanisms underlying the patterning of the bristles and the legs in Gerris buenoi and Mesovelia mulsanti, two species of water striders. We found that two duplication events in the genes beadex and taxi led to a functional expansion of the paralogs, which affected bristle density and leg length. We also identified genes for which no function in bristle development has been previously described in other insects. Interestingly, most of these genes play a dual role in regulating bristle development and leg length. In addition, these genes play a role in regulating cell division. This result suggests that cell division may be a common mechanism through which these genes can simultaneously regulate leg length and bristle density. We propose that pleiotropy, through which gene function affects the development of multiple traits, may play a prominent role in facilitating access to unexploited ecological opportunities and species diversification.

Pleiotropy or linkage? Their relative contributions to the genetic correlation of quantitative traits and detection by multitrait GWA studies

Genetic correlations between traits may cause correlated responses to selection. Previous models described the conditions under which genetic correlations are expected to be maintained. Selection, mutation, and migration are all proposed to affect genetic correlations, regardless of whether the underlying genetic architecture consists of pleiotropic or tightly linked loci affecting the traits. Here, we investigate the conditions under which pleiotropy and linkage have different effects on the genetic correlations between traits by explicitly modeling multiple genetic architectures to look at the effects of selection strength, degree of correlational selection, mutation rate, mutational variance, recombination rate, and migration rate. We show that at mutation-selection(-migration) balance, mutation rates differentially affect the equilibrium levels of genetic correlation when architectures are composed of pairs of physically linked loci compared to architectures of pleiotropic loci. Even when there is perfect linkage (no recombination within pairs of linked loci), a lower genetic correlation is maintained than with pleiotropy, with a lower mutation rate leading to a larger decrease. These results imply that the detection of causal loci in multitrait association studies will be affected by the type of underlying architectures, whereby pleiotropic variants are more likely to be underlying multiple detected associations. We also confirm that tighter linkage between nonpleiotropic causal loci maintains higher genetic correlations at the traits and leads to a greater proportion of false positives in association analyses.

Threespine Stickleback: A Model System For Evolutionary Genomics

The repeated adaptation of oceanic threespine sticklebacks to fresh water has made it a premier organism to study parallel evolution. These small fish have multiple distinct ecotypes that display a wide range of diverse phenotypic traits. Ecotypes are easily crossed in the laboratory, and families are large and develop quickly enough for quantitative trait locus analyses, positioning the threespine stickleback as a versatile model organism to address a wide range of biological questions. Extensive genomic resources, including linkage maps, a high-quality reference genome, and developmental genetics tools have led to insights into the genomic basis of adaptation and the identification of genomic changes controlling traits in vertebrates. Recently, threespine sticklebacks have been used as a model system to identify the genomic basis of highly complex traits, such as behavior and host-microbiome and host-parasite interactions. We review the latest findings and new avenues of research that have led the threespine stickleback to be considered a supermodel of evolutionary genomics.

Diversity and sexual dimorphism in the head lateral line system in North Sea populations of threespine sticklebacks, Gasterosteus aculeatus (Teleostei: Gasterosteidae)

The mechanosensory lateral line of fishes is a flow sensing system and supports a number of behaviors, e.g. prey detection, schooling or position holding in water currents. Differences in the neuromast pattern of this sensory system reflect adaptation to divergent ecological constraints. The threespine stickleback, Gasterosteus aculeatus, is known for its ecological plasticity resulting in three major ecotypes, a marine type, a migrating anadromous type and a resident freshwater type. We provide the first comparative study of the pattern of the head lateral line system of North Sea populations representing these three ecotypes including a brackish spawning population. We found no distinct difference in the pattern of the head lateral line system between the three ecotypes but significant differences in neuromast numbers. The anadromous and the brackish populations had distinctly less neuromasts than their freshwater and marine conspecifics. This difference in neuromast number between marine and anadromous threespine stickleback points to differences in swimming behavior. We also found sexual dimorphism in neuromast number with males having more neuromasts than females in the anadromous, brackish and the freshwater populations. But no such dimorphism occurred in the marine population. Our results suggest that the head lateral line of the three ecotypes is under divergent hydrodynamic constraints. Additionally, sexual dimorphism points to divergent niche partitioning of males and females in the anadromous and freshwater but not in the marine populations. Our findings imply careful sampling as an important prerequisite to discern especially between anadromous and marine threespine sticklebacks.

Allometric relationships in morphological traits associated with foraging, swimming ability, and predator defense reveal adaptations toward brackish and freshwater environments in the threespine stickleback

Freshwater colonization by threespine stickleback has led to divergence in morphology between ancestral marine and derived freshwater populations, making them ideal for studying natural selection on phenotypes. In an open brackish-freshwater system, we previously discovered two genetically distinct stickleback populations that also differ in geometric shape: one mainly found in the brackish water lagoon and one throughout the freshwater system. As shape and size are not perfectly correlated, the aim of this study was to identify the morphological trait(s) that separated the populations in geometric shape. We measured 23 phenotypes likely to be important for foraging, swimming capacity, and defense against predation. The lateral plate morphs in freshwater displayed few significant changes in trait sizes, but the low plated expressed feeding traits more associated with benthic habitats. When comparing the completely plated genetically assigned populations, the freshwater, the hybrids, the migrants and the lagoon fish, many of the linear traits had different slopes and intercepts in trait-size regressions, precluding our ability to directly compare all traits simultaneously, which most likely results from low variation in body length for the lagoon and migrant population. We found the lagoon stickleback population to be more specialized toward the littoral zone, displaying benthic traits such as large, deep bodies with smaller eyes compared to the freshwater completely plated morph. Further, the lagoon and migrant fish had an overall higher body coverage of lateral plates compared to freshwater fish, and the dorsal and pelvic spines were longer. Evolutionary constraints due to allometric scaling relationships could explain the observed, overall restricted, differences in morphology between the sticklebacks in this study, as most traits have diversified in common allometric trajectories. The observed differences in foraging and antipredation traits between the fish with a lagoon and freshwater genetic signature are likely a result of genetic or plastic adaptations toward brackish and freshwater environments.

Phosphorus limitation does not drive loss of bony lateral plates in freshwater stickleback (Gasterosteus aculeatus)

Connecting the selective forces that drive the evolution of phenotypes to their underlying genotypes is key to understanding adaptation, but such connections are rarely tested experimentally. Threespine stickleback (Gasterosteus aculeatus) are a powerful model for such tests because genotypes that underlie putatively adaptive traits have been identified. For example, a regulatory mutation in the Ectodysplasin (Eda) gene causes a reduction in the number of bony armor plates, which occurs rapidly and repeatedly when marine sticklebacks invade freshwater. However, the source of selection on plate loss in freshwater is unknown. Here, we tested whether dietary reduction of phosphorus can account for selection on plate loss due to a growth advantage of low-plated fish in freshwater. We crossed marine fish heterozygous for the 16 kilobase freshwater Eda haplotype and compared the growth of offspring with different genotypes under contrasting levels of dietary phosphorus in both saltwater and freshwater. Eda genotype was not associated with growth differences in any treatment, or with mechanisms that could mitigate the impacts of phosphorus limitation, such as differential phosphorus deposition, phosphorus excretion, or intestine length. This study highlights the importance of experimentally testing the putative selective forces acting on phenotypes and their underlying genotypes in the wild.

Adaptation via pleiotropy and linkage: Association mapping reveals a complex genetic architecture within the stickleback Eda locus

Genomic mapping of the loci associated with phenotypic evolution has revealed genomic "hotspots," or regions of the genome that control multiple phenotypic traits. This clustering of loci has important implications for the speed and maintenance of adaptation and could be due to pleiotropic effects of a single mutation or tight genetic linkage of multiple causative mutations affecting different traits. The threespine stickleback (Gasterosteus aculeatus) is a powerful model for the study of adaptive evolution because the marine ecotype has repeatedly adapted to freshwater environments across the northern hemisphere in the last 12,000 years. Freshwater ecotypes have repeatedly fixed a 16 kilobase haplotype on chromosome IV that contains Ectodysplasin (Eda), a gene known to affect multiple traits, including defensive armor plates, lateral line sensory hair cells, and schooling behavior. Many additional traits have previously been mapped to a larger region of chromosome IV that encompasses the Eda freshwater haplotype. To identify which of these traits specifically map to this adaptive haplotype, we made crosses of rare marine fish heterozygous for the freshwater haplotype in an otherwise marine genetic background. Further, we performed fine-scale association mapping in a fully interbreeding, polymorphic population of freshwater stickleback to disentangle the effects of pleiotropy and linkage on the phenotypes affected by this haplotype. Although we find evidence that linked mutations have small effects on a few phenotypes, a small 1.4-kb region within the first intron of Eda has large effects on three phenotypic traits: lateral plate count, and both the number and patterning of the posterior lateral line neuromasts. Thus, the Eda haplotype is a hotspot of adaptation in stickleback due to both a small, pleiotropic region affecting multiple traits as well as multiple linked mutations affecting additional traits.

Spatial, sexual, and rapid temporal differentiation in neuromast expression on lateral plates of Haida Gwaii threespine stickleback (Gasterosteus aculeatus)

Lateral lines, a major sensory modality in fishes, are diverse among taxa, but their intraspecific variation has received limited attention. We examined numbers of superficial neuromasts on the buttressing lateral plates (LP) of 1910 threespine stickleback (Gasterosteus aculeatus Linnaeus, 1758) from 26 ecologically and morphologically diverse populations on the Haida Gwaii archipelago, western Canada. Extending from previous studies, we predicted that (i) highly stained dystrophic localities would have threespine stickleback with elevated numbers of neuromasts per plate due to a greater reliance on non-visual sensory modalities and (ii) that LP count and neuromast numbers per plate would functionally covary with predatory assemblage. We found that there were no differences in neuromast count across major habitats (marine, lake, stream), but clear-water populations and those with predatory fish had significantly more neuromasts per plate than most populations in highly stained dystrophic lakes, the effects being accentuated on the first buttressing plate (LP4). We also report the first evidence that neuromast counts per plate are sexually dimorphic, with males having a greater density of neuromasts in most populations. Two transplant experiments between ecologically opposite habitats indicate that within 12 generations, neuromast counts per plate can rapidly shift in response to a change in habitat.

Investigating the evolution and development of biological complexity under the framework of epigenetics

Biological complexity is a key component of evolvability, yet its study has been hampered by a focus on evolutionary trends of complexification and inconsistent definitions. Here, we demonstrate the utility of bringing complexity into the framework of epigenetics to better investigate its utility as a concept in evolutionary biology. We first analyze the existing metrics of complexity and explore the link between complexity and adaptation. Although recently developed metrics allow for a unified framework, they omit developmental mechanisms. We argue that a better approach to the empirical study of complexity and its evolution includes developmental mechanisms. We then consider epigenetic mechanisms and their role in shaping developmental and evolutionary trajectories, as well as the development and organization of complexity. We argue that epigenetics itself could have emerged from complexity because of a need to self‐regulate. Finally, we explore hybridization complexes and hybrid organisms as potential models for studying the association between epigenetics and complexity. Our goal is not to explain trends in biological complexity but to help develop and elucidate novel questions in the investigation of biological complexity and its evolution.

This manuscript argues that biological complexity is better understood under the framework of epigenetics and that the epigenetic interactions emerge from the self‐regulation of complex systems. Hybrids are offered as models to study these properties.

Adaptive Diversification of the Lateral Line System during Cichlid Fish Radiation

The mechanosensory lateral line system is used by fishes to sense hydrodynamic stimuli in their environment. It provides information about flow regimes, proximity to substrate, and the presence and identity of prey and predators and represents a means of receiving communication signals from other fish. Thus we may expect lateral line system structures to be under strong divergent selection during adaptive radiation. Here, we used X-ray micro-computed tomography scans to quantify variation in cranial lateral line canal morphology within the adaptive radiation of Lake Malawi cichlids. We report that cranial lateral line canal morphology is strongly correlated with diet and other aspects of craniofacial morphology, including the shape of oral jaws. These results indicate an adaptive role for the lateral line system in prey detection and suggest that diversification of this system has taken an important role in the spectacular evolution of Lake Malawi's cichlid fish diversity.

Efficient CRISPR-Cas9 editing of major evolutionary loci in sticklebacks

BACKGROUND

Stickleback fish are widely used to study the genetic and ecological basis of phenotypic evolution. Although several major loci have now been identified that contribute to evolutionary differences between wild populations, further study of the phenotypes associated with particular genes and mutations has been limited by the difficulty of generating targeted mutations at precise locations in the stickleback genome.

APPROACH AND AIMS

We compared different methods of expressing single-guide RNAs (sgRNAs) and Cas9 activity in fertilized stickleback eggs. We used an easily scored pigmentation gene (SLC24A5) to screen for molecular lesions, phenotypic effects, and possible germline transmission of newly induced alleles. We then used the optimized CRISPR methods to target two major evolutionary loci in sticklebacks, KITLG and EDA. We hypothesized that coding region mutations in the KITLG gene would alter body pigmentation and possibly sex determination, and that mutations in the EDA gene would disrupt the formation of most armor plates, fin rays, spines, teeth, and gill rakers.

RESULTS

Targeted deletions were successfully induced at each target locus by co-injecting one-cell stage stickleback embryos with either Cas9 mRNA or Cas9 protein, together with sgRNAs designed to protein-coding exons. Founder animals were typically mosaic for multiple mutations, which they transmitted through the germline at overall rates of 21 to 100%. We found that the copy of KITLG on the X chromosome (KITLGX) has diverged from the KITLG on the Y chromosome (KITLGY). Predicted loss-of-function mutations in the KITLGX gene dramatically altered pigmentation in both external skin and internal organ, but the same was not true for KITLGY mutations. Predicted loss-of-function mutations in either the KITLGX or KITLGY genes did not lead to sex reversal or prevent fertility. Homozygous loss-of-function mutations in the EDA gene led to complete loss of armor plates, severe reduction or loss of most soft rays in the dorsal, anal, and caudal fins, and severe reductions in tooth and gill raker number. In contrast, long dorsal and pelvic spines remained intact in EDA mutant animals, suggesting that common co-segregation of plate loss and spine reduction in wild populations is unlikely to be due to pleiotropic effects of EDA mutations.

CONCLUSION

CRISPR-Cas9 approaches can be used to induce germline mutations in key evolutionary loci in sticklebacks. Targeted coding region mutations confirm an important role for KITLG and EDA in skin pigmentation and armor plate reduction, respectively. They also provide new information about the functions of these genes in other body structures.

Investigating the association between armour coverage and parasite infection in an estuarine population of stickleback

BACKGROUND

When threespine stickleback colonized fresh water, they repeatedly evolved reduced armour plating via changes in Eda allele frequency. This evolution is typically attributed to differential predation pressure between marine and freshwater environments. However, the chromosomal region containing Eda is associated with many other phenotypes, including schooling, antipredator behaviour, and immunity. Consequently, the evolution of armour plating may be driven by multiple selective pressures acting on Eda or linked genes.

QUESTION

Is parasite infection associated with armour phenotype?

HYPOTHESIS

Parasite load differs between stickleback armour plate morphs.

ORGANISMS

An armour-polymorphic population of threespine stickleback (Gasterosteus aculeatus), and their parasites.

FIELD SITE

In June 2009 and 2012, we sampled stickleback from a single human-made salt-marsh pool in the Campbell River Estuary on Vancouver Island.

METHODS

We counted macroparasites on approximately 100 fish per year and counted lateral armour plates. We used generalized linear models to test for correlations between armour morph and parasite load.

RESULTS

Most parasite species were not associated with armour. The gill parasite Thersitina was more abundant on more fully armoured fish in both years. The nematode Eustrongylides also exhibited a marginally significant positive trend. If parasitic infections reduce stickleback fitness, this positive covariance between armour and infection would accelerate the loss of armour plating in stickleback colonizing fresh water.

Embryonic origin and genetic basis of cave associated phenotypes in the isopod crustacean Asellus aquaticus

Characteristics common to animals living in subterranean environments include the reduction or absence of eyes, lessened pigmentation and enhanced sensory systems. How these characteristics have evolved is poorly understood for the majority of cave dwelling species. In order to understand the evolution of these changes, this study uses an invertebrate model system, the freshwater isopod crustacean, Asellus aquaticus, to examine whether adult differences between cave and surface dwelling individuals first appear during embryonic development. We hypothesized that antennal elaboration, as well as eye reduction and pigment loss, would be apparent during embryonic development. We found that differences in pigmentation, eye formation, and number of segments of antenna II were all present by the end of embryonic development. In addition, we found that cave and surface hatchlings do not significantly differ in the relative size of antenna II and the duration of embryonic development. To investigate whether the regions responsible for eye and pigment differences could be genetically linked to differences in article number, we genotyped F2 hybrids for the four previously mapped genomic regions associated with eye and pigment differences and phenotyped these F2 hybrids for antenna II article number. We found that the region previously known to be responsible for both presence versus absence of pigment and eye size also was significantly associated with article number. Future experiments will address whether pleiotropy and/or genetic linkage play a role in the evolution of cave characteristics in Asellus aquaticus.

Contemporary ancestor? Adaptive divergence from standing genetic variation in Pacific marine threespine stickleback

Background

Populations that have repeatedly colonized novel environments are useful for studying the role of ecology in adaptive divergence – particularly if some individuals persist in the ancestral habitat. Such “contemporary ancestors” can be used to demonstrate the effects of selection by comparing phenotypic and genetic divergence between the derived population and their extant ancestors. However, evolution and demography in these “contemporary ancestors” can complicate inferences about the source (standing genetic variation, de novo mutation) and pace of adaptive divergence. Marine threespine stickleback (Gasterosteus aculeatus) have colonized freshwater environments along the Pacific coast of North America, but have also persisted in the marine environment. To what extent are marine stickleback good proxies of the ancestral condition?

Results

We sequenced > 5800 variant loci in over 250 marine stickleback from eight locations extending from Alaska to California, and phenotyped them for platedness and body shape. Pairwise F_ST varied from 0.02 to 0.18. Stickleback were divided into five genetic clusters, with a single cluster comprising stickleback from Washington to Alaska. Plate number, Eda, body shape, and candidate loci showed evidence of being under selection in the marine environment. Comparisons to a freshwater population demonstrated that candidate loci for freshwater adaptation varied depending on the choice of marine populations.

Conclusions

Marine stickleback are structured into phenotypically and genetically distinct populations that have been evolving as freshwater stickleback evolved. This variation complicates their usefulness as proxies of the ancestors of freshwater populations. Lessons from stickleback may be applied to other “contemporary ancestor”-derived population studies.

Electronic supplementary material

The online version of this article (10.1186/s12862-018-1228-8) contains supplementary material, which is available to authorized users.

Evolution and Genetics of Precocious Burrowing Behavior in Peromyscus Mice

A central challenge in biology is to understand how innate behaviors evolve between closely related species. One way to elucidate how differences arise is to compare the development of behavior in species with distinct adult traits [1]. Here, we report that Peromyscus polionotus is strikingly precocious with regard to burrowing behavior, but not other behaviors, compared to its sister species P. maniculatus. In P. polionotus, burrows were excavated as early as 17 days of age, whereas P. maniculatus did not build burrows until 10 days later. Moreover, the well-known differences in burrow architecture between adults of these species-P. polionotus adults excavate long burrows with an escape tunnel, whereas P. maniculatus dig short, single-tunnel burrows [2-4]-were intact in juvenile burrowers. To test whether this juvenile behavior is influenced by early-life environment, we reciprocally cross-fostered pups of both species. Fostering did not alter the characteristic burrowing behavior of either species, suggesting that these differences are genetic. In backcross hybrids, we show that precocious burrowing and adult tunnel length are genetically correlated and that a P. polionotus allele linked to tunnel length variation in adults is also associated with precocious onset of burrowing in juveniles, suggesting that the same genetic region-either a single gene with pleiotropic effects or linked genes-influences distinct aspects of the same behavior at these two life stages. These results raise the possibility that genetic variants affect behavioral drive (i.e., motivation) to burrow and thereby affect both the developmental timing and adult expression of burrowing behavior.

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'.

Eda haplotypes in three-spined stickleback are associated with variation in immune gene expression

Haplotypes underlying local adaptation and speciation are predicted to have numerous phenotypic effects, but few genes involved have been identified, with much work to date concentrating on visible, morphological, phenotypes. The link between genes controlling these adaptive morphological phenotypes and the immune system has seldom been investigated, even though changes in the immune system could have profound adaptive consequences. The Eda gene in three-spined stickleback is one of the best studied major adaptation genes; it directly controls bony plate architecture and has been associated with additional aspects of adaptation to freshwater. Here, we exposed F2 hybrids, used to separate Eda genotype from genetic background, to contrasting conditions in semi-natural enclosures. We demonstrate an association between the Eda haplotype block and the expression pattern of key immune system genes. Furthermore, low plated fish grew less and experienced higher burdens of a common ectoparasite with fitness consequences. Little is currently known about the role of the immune system in facilitating adaptation to novel environments, but this study provides an indication of its potential importance.

Trait Correlations in the Genomics Era

Thinking about the evolutionary causes and consequences of trait correlations has been dominated by quantitative genetics theory that is focused on hypothetical loci. Since this theory was initially developed, technology has enabled the identification of specific genetic variants that contribute to trait correlations. Here, we review studies of the genetic basis of trait correlations to ask: What has this new information taught us? We find that causal variants can be pleiotropic and/or linked in different ways, indicating that pleiotropy and linkage are not alternative genetic mechanisms. Further, many trait correlations have a polygenic basis, suggesting that both pleiotropy and linkage likely contribute. We discuss implications of these findings for the evolutionary causes and consequences of trait correlations.

The Genetic Basis of Upland/Lowland Ecotype Divergence in Switchgrass (Panicum virgatum)

The evolution of locally adapted ecotypes is a common phenomenon that generates diversity within plant species. However, we know surprisingly little about the genetic mechanisms underlying the locally adapted traits involved in ecotype formation. The genetic architecture underlying locally adapted traits dictates how an organism will respond to environmental selection pressures, and has major implications for evolutionary ecology, conservation, and crop breeding. To understand the genetic architecture underlying the divergence of switchgrass (Panicum virgatum) ecotypes, we constructed a genetic mapping population through a four-way outbred cross between two northern upland and two southern lowland accessions. Trait segregation in this mapping population was largely consistent with multiple independent loci controlling the suite of traits that characterizes ecotype divergence. We assembled a joint linkage map using ddRADseq, and mapped quantitative trait loci (QTL) for traits that are divergent between ecotypes, including flowering time, plant size, physiological processes, and disease resistance. Overall, we found that most QTL had small to intermediate effects. While we identified colocalizing QTL for multiple traits, we did not find any large-effect QTL that clearly controlled multiple traits through pleiotropy or tight physical linkage. These results indicate that ecologically important traits in switchgrass have a complex genetic basis, and that similar loci may underlie divergence across the geographic range of the ecotypes.

Variation in Lateral Plate Quality in Threespine Stickleback from Fresh, Brackish and Marine Water: A Micro-Computed Tomography Study

Introduction

It is important to understand the drivers leading to adaptive phenotypic diversity within and among species. The threespine stickleback (Gasterosteus aculeatus) has become a model system for investigating the genetic and phenotypic responses during repeated colonization of fresh waters from the original marine habitat. During the freshwater colonization process there has been a recurrent and parallel reduction in the number of lateral bone plates, making it a suitable system for studying adaptability and parallel evolution.

Objective

The aim of this study was to investigate an alternative evolutionary path of lateral plate reduction, where lateral plates are reduced in size rather than number.

Materials and Methods

A total of 72 threespine stickleback individuals from freshwater (n = 54), brackish water (n = 27) and marine water (n = 9) were analysed using microcomputed tomography (μCT) to determine variation in size, thickness and structure of the lateral plates. Furthermore, whole-body bone volume, and bone volume, bone surface and porosity of lateral plate number 4 were quantified in all specimens from each environment.

Results

The results showed a significant difference in plate size (area and volume) among populations, where threespine stickleback from polymorphic freshwater and brackish water populations displayed lateral plates reduced in size (area and volume) compared to marine stickleback

Conclusions

Reduction of lateral plates in threespine stickleback in fresh and brackish water occurs by both plate loss and reduction in plate size (area and volume).

The genetic architecture of local adaptation and reproductive isolation in sympatry within the Mimulus guttatus species complex

The genetic architecture of local adaptation has been of central interest to evolutionary biologists since the modern synthesis. In addition to classic theory on the effect size of adaptive mutations by Fisher, Kimura and Orr, recent theory addresses the genetic architecture of local adaptation in the face of ongoing gene flow. This theory predicts that with substantial gene flow between populations local adaptation should proceed primarily through mutations of large effect or tightly linked clusters of smaller effect loci. In this study, we investigate the genetic architecture of divergence in flowering time, mating system-related traits, and leaf shape between Mimulus laciniatus and a sympatric population of its close relative M. guttatus. These three traits are probably involved in M. laciniatus' adaptation to a dry, exposed granite outcrop environment. Flowering time and mating system differences are also reproductive isolating barriers making them 'magic traits'. Phenotypic hybrids in this population provide evidence of recent gene flow. Using next-generation sequencing, we generate dense SNP markers across the genome and map quantitative trait loci (QTLs) involved in flowering time, flower size and leaf shape. We find that interspecific divergence in all three traits is due to few QTL of large effect including a highly pleiotropic QTL on chromosome 8. This QTL region contains the pleiotropic candidate gene TCP4 and is involved in ecologically important phenotypes in other Mimulus species. Our results are consistent with theory, indicating that local adaptation and reproductive isolation with gene flow should be due to few loci with large and pleiotropic effects.

Ectodysplasin signalling genes and phenotypic evolution in sculpins (Cottus)

Despite their deeply conserved function among vertebrates, ectodysplasin (Eda) signalling genes are involved in microevolutionary change in humans and sticklebacks. If such a dual role is common, Eda signalling genes constitute hotspots for morphological evolution. Variation in sculpin (Cottus) skin prickling and body shape resembles patterns caused by variation in Eda signalling in sticklebacks. We mapped Eda signalling genes and performed quantitative trait locus mapping in crosses between Cottus rhenanus and Cottus perifretum. A genomic region containing the Eda receptor (Edar) was strongly associated with prickling and contributed to shape. The expression of Edar in developing prickles and skeletal elements in Cottus was confirmed by in situ hybridization. Coding sequence changes between Edar alleles in C. rhenanus and C. perifretum exceeded sequence differentiation in other vertebrates. However, it is likely that additional genetic elements besides coding changes affect the phenotypic variation. Although the phenotype in a natural hybrid lineage between C. rhenanus and C. perifretum resembles C. perifretum, the respective coding Edar alleles are not fully fixed (88.6%). Hence, our results support an involvement of Eda signalling in microevolutionary changes, but imply that the Edar gene is affected by multiple evolutionary processes that vary among freshwater sculpins.

Factors influencing the effect size distribution of adaptive substitutions

The distribution of effect sizes of adaptive substitutions has been central to evolutionary biology since the modern synthesis. Early theory proposed that because large-effect mutations have negative pleiotropic consequences, only small-effect mutations contribute to adaptation. More recent theory suggested instead that large-effect mutations could be favoured when populations are far from their adaptive peak. Here we suggest that the distributions of effect sizes are expected to differ among study systems, reflecting the wide variation in evolutionary forces and ecological conditions experienced in nature. These include selection, mutation, genetic drift, gene flow, and other factors such as the degree of pleiotropy, the distance to the phenotypic optimum, whether the optimum is stable or moving, and whether new mutation or standing genetic variation provides the source of adaptive alleles. Our goal is to review how these factors might affect the distribution of effect sizes and to identify new research directions. Until more theory and empirical work is available, we feel that it is premature to make broad generalizations about the effect size distribution of adaptive substitutions important in nature.

Evolution of Schooling Behavior in Threespine Sticklebacks Is Shaped by the Eda Gene

Despite longstanding interest in the genetic mechanisms that underlie behavioral evolution, very few genes that underlie naturally occurring variation in behavior between individuals or species are known, particularly in vertebrates. Here, we build on our previous forward genetic mapping experiments and use transgenic approaches to identify Ectodysplasin as a gene that causes differences in schooling behavior between wild populations of threespine stickleback (Gasterosteus aculeatus) fish. This work provides rare insight into the proximate mechanisms that have shaped the evolution of vertebrate behavior.

How the effects of aging and stresses of life are integrated in mortality rates: insights for genetic studies of human health and longevity

Increasing proportions of elderly individuals in developed countries combined with substantial increases in related medical expenditures make the improvement of the health of the elderly a high priority today. If the process of aging by individuals is a major cause of age related health declines then postponing aging could be an efficient strategy for improving the health of the elderly. Implementing this strategy requires a better understanding of genetic and non-genetic connections among aging, health, and longevity. We review progress and problems in research areas whose development may contribute to analyses of such connections. These include genetic studies of human aging and longevity, the heterogeneity of populations with respect to their susceptibility to disease and death, forces that shape age patterns of human mortality, secular trends in mortality decline, and integrative mortality modeling using longitudinal data. The dynamic involvement of genetic factors in (i) morbidity/mortality risks, (ii) responses to stresses of life, (iii) multi-morbidities of many elderly individuals, (iv) trade-offs for diseases, (v) genetic heterogeneity, and (vi) other relevant aging-related health declines, underscores the need for a comprehensive, integrated approach to analyze the genetic connections for all of the above aspects of aging-related changes. The dynamic relationships among aging, health, and longevity traits would be better understood if one linked several research fields within one conceptual framework that allowed for efficient analyses of available longitudinal data using the wealth of available knowledge about aging, health, and longevity already accumulated in the research field.

Commentary: When does understanding phenotypic evolution require identification of the underlying genes?

Adaptive evolution is fundamentally a genetic process. Over the past three decades, characterizing the genes underlying adaptive phenotypic change has revealed many important aspects of evolutionary change. At the same time, natural selection is often fundamentally an ecological process that can often be studied without identifying the genes underlying the variation on which it acts. This duality has given rise to disagreement about whether, and under what circumstances, it is necessary to identify specific genes associated with phenotypic change. This issue is of practical concern, especially for researchers who study nonmodel organisms, because of the often enormous cost and labor required to "go for the genes." We here consider a number of situations and questions commonly addressed by researchers. Our conclusion is that although gene identification can be crucial for answering some questions, there are others for which definitive answers can be obtained without finding underlying genes. It should thus not be assumed that considerations of "empirical completeness" dictate that gene identification is always desirable.

Discriminating selection on lateral plate phenotype and its underlying gene, Ectodysplasin, in threespine stickleback

When a signature of natural selection is discovered on a gene that is pleiotropic or in tight linkage with other genes, it is challenging to determine which of the affected phenotypes is under selection. One way to make progress is to employ methods for analyzing natural selection on correlated traits, including both genotype and phenotype. We used this approach in threespine stickleback to estimate selection on a rapidly evolving trait, lateral armor plates, while controlling for variation at its major underlying gene, Ectodysplasin (Eda), and vice versa. This allowed for independent estimates of selection on lateral plates and on Eda via other traits. Previously, we demonstrated allele frequency changes at Eda in a pond experiment. Here we show that this resulted from selection on both plates and on Eda, implying additional selection on other phenotypic traits affected by the same gene. This represents the first evidence for direct selection on lateral plates independent of selection on the Eda locus and highlights the value of measuring selection on both phenotypes and genotypes in studies of adaptation.

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.