Genetic mapping of natural variation in schooling tendency in the threespine stickleback

Extreme intraspecific divergence in mitochondrial haplotypes makes the threespine stickleback fish an emerging evolutionary mutant model for mito-nuclear interactions

Mitochondrial DNA is primarily maternally inherited in most animals and evolves about 10 times faster than biparentally inherited nuclear DNA. Mitochondrial dysfunction (mt-dys) arises when interactions between the co-evolving mitochondrial and nuclear genomes are perturbed in essential processes like oxidative phosphorylation (OXPHOS). Over time mt-dys can lead to mitochondrial diseases (mt-diseases), which are surprisingly prevalent and include common diseases such as Alzheimer’s, Parkinson’s, and diabetes. Unfortunately, the strong impact that intraspecific mitochondrial and nuclear genetic variation has on mt-disease complicates its study and the development of effective treatments. Animal models have advanced our understanding of mt-disease but their relevance to human conditions is often limited by their relatively low nuclear genetic diversity. Many traditional laboratory models also typically have a single mitochondrial haplotype (mitotype), in stark contrast to over 5,000 mitotypes in humans worldwide. The threespine stickleback fish has an evolutionary history that has made it a favorable evolutionary mutant model (EMM) for studying mito-nuclear interactions and possibly mt-diseases. EMMs are species with naturally evolved states that mimic maladaptive human diseases. In threespine stickleback, a period of isolation followed by introgression of the mitochondrial genome from a sister species resulted in the maintenance of two distinct mitochondrial haplotypes which continue to segregate within many populations of wild stickleback. The existence of two mitogenomes segregating in numerous genetically diverse populations provides a unique system for exploring complex mito-nuclear dynamics. Here we provide the first complete coding region analysis of the two threespine stickleback mitotypes, whose mitogenomic divergence exceeds that of other mammalian models for mitochondrial disease and even that between ancient and modern humans. We find that divergence is not uniform across the mitogenome, but primarily impacts protein coding genes, and significantly impacts proteins in Complex I of OXPHOS. The full characterization of these highly divergent intraspecific mitotypes provides a foundation for the development of threespine stickleback as an EMM for mito-nuclear interactions.

Kinematic analysis of social interactions deconstructs the evolved loss of schooling behavior in cavefish

Fish display a remarkable diversity of social behaviors, both within and between species. While social behaviors are likely critical for survival, surprisingly little is known about how they evolve in response to changing environmental pressures. With its highly social surface form and multiple populations of a largely asocial, blind, cave-dwelling form, the Mexican tetra, Astyanax mexicanus, provides a powerful model to study the evolution of social behavior. Here we use motion tracking and analysis of swimming kinematics to quantify social swimming in four Astyanax mexicanus populations. In the light, surface fish school, maintaining both close proximity and alignment with each other. In the dark, surface fish no longer form coherent schools, however, they still show evidence of an attempt to align and maintain proximity when they find themselves near another fish. In contrast, cavefish from three independently-evolved populations (Pachón, Molino, Tinaja) show little preference for proximity or alignment, instead exhibiting behaviors that suggest active avoidance of each other. Two of the three cave populations we studied also slow down when more fish are present in the tank, a behavior which is not observed in surface fish in light or the dark, suggesting divergent responses to conspecifics. Using data-driven computer simulations, we show that the observed reduction in swimming speed is sufficient to alter the way fish explore their environment: it can increase time spent exploring away from the walls. Thus, the absence of schooling in cavefish is not merely a consequence of their inability to see, but may rather be a genuine behavioral adaptation that impacts the way they explore their environment.

The Collective Behavior of Ant Groups Depends on Group Genotypic Composition

Recently, researchers have documented variation between groups in collective behavior. However, how genetic variation within and between groups contributes to population-level variation for collective behavior remains unclear. Understanding how genetic variation of group members relates to group-level phenotypes is evolutionarily important because there is increasing evidence that group-level behavioral variation influences fitness and that the genetic architecture of group-level traits can affect the evolutionary dynamics of traits. Social insects are ideal for studying the complex relationship between individual and group-level variation because they exhibit behavioral variation at multiple scales of organization. To explore how the genetic composition of groups affects collective behavior, we constructed groups of pharaoh ants (Monomorium pharaonis) from 33 genetically distinct colonies of known pedigree. The groups consisted of either all workers from the same single colony or workers from two genetically different colonies, and we assayed the exploration and aggression of the groups. We found that collective exploration, but not aggression, depended on the specific genotypic combination of group members, i.e., we found evidence for genotype-by-genotype epistasis for exploration. Group collective behavior did not depend on the pedigree relatedness between genotypes within groups. Overall, this study highlights that specific combinations of genotypes influence group-level phenotypes, emphasizing the importance of considering nonadditive effects of genotypic interactions between group members.

Expanding evolutionary neuroscience: insights from comparing variation in behavior

Neuroscientists have long studied species with convenient biological features to discover how behavior emerges from conserved molecular, neural, and circuit level processes. With the advent of new tools, from viral vectors and gene editing to automated behavioral analyses, there has been a recent wave of interest in developing new, "nontraditional" model species. Here, we advocate for a complementary approach to model species development, that is, model clade development, as a way to integrate an evolutionary comparative approach with neurobiological and behavioral experiments. Capitalizing on natural behavioral variation in and investing in experimental tools for model clades will be a valuable strategy for the next generation of neuroscience discovery.

QTL Mapping of Intestinal Neutrophil Variation in Threespine Stickleback Reveals Possible Gene Targets Connecting Intestinal Inflammation and Systemic Health

Selection, via host immunity, is often required to foster beneficial microbial symbionts and suppress deleterious pathogens. In animals, the host immune system is at the center of this relationship. Failed host immune system-microbial interactions can result in a persistent inflammatory response in which the immune system indiscriminately attacks resident microbes, and at times the host cells themselves, leading to diseases such as Ulcerative Colitis, Crohn’s Disease, and Psoriasis. Host genetic variation has been linked to both microbiome diversity and to severity of such inflammatory disease states in humans. However, the microbiome and inflammatory states manifest as quantitative traits, which encompass many genes interacting with one another and the environment. The mechanistic relationships among all of these interacting components are still not clear. Developing natural genetic models of host-microbe interactions is therefore fundamental to understanding the complex genetics of these and other diseases. Threespine stickleback (Gasterosteus aculeatus) fish are a tractable model for attacking this problem because of abundant population-level genetic and phenotypic variation in the gut inflammatory response. Previous work in our laboratory identified genetically divergent stickleback populations exhibiting differences in intestinal neutrophil activity. We took advantage of this diversity to genetically map variation in an emblematic element of gut inflammation - intestinal neutrophil recruitment - using an F2-intercross mapping framework. We identified two regions of the genome associated with increased intestinal inflammation containing several promising candidate genes. Within these regions we found candidates in the Coagulation/Complement System, NFkB and MAPK pathways along with several genes associated with intestinal diseases and neurological diseases commonly accompanying intestinal inflammation as a secondary symptom. These findings highlight the utility of using naturally genetically diverse ‘evolutionary mutant models’ such as threespine stickleback to better understand interactions among host genetic diversity and microbiome variation in health and disease states.

Quantitative Trait Loci (QTL) Mapping

Quantitative trait loci (QTL) are genetic regions that influence phenotypic variation of a complex trait, often through genetic interactions with each other and the environment. These are commonly identified through a statistical genetic analysis known as QTL mapping. Here, I present a step-by-step, practical approach to QTL mapping along with a sample data file. I focus on methods commonly used and discoveries that have been made in fishes, and utilize a multiple QTL mapping (MQM) approach in the free software package R/qtl.

Social environment determines the effect of boldness and activity on survival

Populations of animals are composed of individuals that differ in ecologically relevant behaviors. Building evidence also suggests that individuals occupy different social niches. Here, in a mark-recapture experiment, we show evidence of an interacting effect of behavior and social niche on survival in the wild: bold individuals had higher survival if they were initially captured in groups while shy, inactive individuals had higher survival if they were initially captured when alone. These findings provide support for the hypothesis that behavioral type-environment correlations can be favored by natural selection.

The yellow gene influences Drosophila male mating success through sex comb melanization

Drosophila melanogaster males perform a series of courtship behaviors that, when successful, result in copulation with a female. For over a century, mutations in the yellow gene, named for its effects on pigmentation, have been known to reduce male mating success. Prior work has suggested that yellow influences mating behavior through effects on wing extension, song, and/or courtship vigor. Here, we rule out these explanations, as well as effects on the nervous system more generally, and find instead that the effects of yellow on male mating success are mediated by its effects on pigmentation of male-specific leg structures called sex combs. Loss of yellow expression in these modified bristles reduces their melanization, which changes their structure and causes difficulty grasping females prior to copulation. These data illustrate why the mechanical properties of anatomy, not just neural circuitry, must be considered to fully understand the development and evolution of behavior.

More than 100 years ago, Nobel-prize winning geneticist Thomas Hunt Morgan and his colleagues discovered that some fruit flies inherited genetic mutations that caused their body color to change. The yellow flies had a mutation in one specific gene and these mutants did not only look different from normal flies, they behaved differently too. Specifically, yellow males were far less successful at mating than normal males, demonstrating for the first time that some behaviors had a genetic basis.

Since then it has remained a mystery how the genetic mutations that cause yellow coloration in fruit flies lead to unsuccessful mating attempts. Geneticists have long suggested that mutations in insect pigment genes cause changes in the fly’s brain because these pigments are made from dopamine, a chemical messenger that acts in the brain. They proposed that yellow flies must have altered levels of dopamine in their brains which was causing them to fail at mating.

To solve this mystery, Massey et al. used a series of genetic experiments and high speed-videos to assess how mutations in male yellow fruit flies affected their mating behavior. The experiments showed that yellow fruit flies mated poorly not because of changes in their brain but because of changes in specialized structures on their legs called sex combs. The yellow males lack melanin pigments in their sex combs, which changes their structure. As a result, the yellow males would court female flies but were then unable to grab and mount them. This explains why yellow flies often fail to mate and why fruit flies have sex combs in the first place.

The study reveals the importance of scientists considering that genes that affect behavior may do so by changing anatomy rather than by altering the brain. The results also may benefit those working to control insect pests. For example, they could help insect pest managers to develop strategies that prevent reproduction in other insects that spread disease or destroy crops.

Evolutionary diversity as a catalyst for biological discovery

The tremendous diversity of animal behaviors has inspired generations of scientists from an array of biological disciplines. To complement investigations of ecological and evolutionary factors contributing to behavioral evolution, modern sequencing, gene editing, computational and neuroscience tools now provide a means to discover the proximate mechanisms upon which natural selection acts to generate behavioral diversity. Social behaviors are motivated behaviors that can differ tremendously between closely related species, suggesting phylogenetic plasticity in their underlying biological mechanisms. In addition, convergent evolution has repeatedly given rise to similar forms of social behavior and mating systems in distantly related species. Social behavioral divergence and convergence provides an entry point for understanding the neurogenetic mechanisms contributing to behavioral diversity. We argue that the greatest strides in discovering mechanisms contributing to social behavioral diversity will be achieved through integration of interdisciplinary comparative approaches with modern tools in diverse species systems. We review recent advances and future potential for discovering mechanisms underlying social behavioral variation; highlighting patterns of social behavioral evolution, oxytocin and vasopressin neuropeptide systems, genetic/transcriptional "toolkits," modern experimental tools, and alternative species systems, with particular emphasis on Microtine rodents and Lake Malawi cichlid fishes.

Genetic analysis reveals candidate genes for activity QTL in the blind Mexican tetra, Astyanax mexicanus

Animal models provide useful tools for exploring the genetic basis of morphological, physiological and behavioral phenotypes. Cave-adapted species are particularly powerful models for a broad array of phenotypic changes with evolutionary, developmental and clinical relevance. Here, we explored the genetic underpinnings of previously characterized differences in locomotor activity patterns between the surface-dwelling and Pachón cave-dwelling populations of Astyanax mexicanus. We identified multiple novel QTL underlying patterns in overall levels of activity (velocity), as well as spatial tank use (time spent near the top or bottom of the tank). Further, we demonstrated that different regions of the genome mediate distinct patterns in velocity and tank usage. We interrogated eight genomic intervals underlying these activity QTL distributed across six linkage groups. In addition, we employed transcriptomic data and draft genomic resources to generate and evaluate a list of 36 potential candidate genes. Interestingly, our data support the candidacy of a number of genes, but do not suggest that differences in the patterns of behavior observed here are the result of alterations to certain candidate genes described in other species (e.g., teleost multiple tissue opsins, melanopsins or members of the core circadian clockwork). This study expands our knowledge of the genetic architecture underlying activity differences in surface and cavefish. Future studies will help define the role of specific genes in shaping complex behavioral phenotypes in Astyanax and other vertebrate taxa.

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.

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.

A standardized method for quantifying consistent individual differences in schooling behaviour

A method for quantifying consistent individual differences in schooling behaviour is presented. This method, which utilizes a school of models, improves on previous methods by removing the unwanted variation that is introduced by live stimulus fish while still providing the physiological experience of schooling to the focal fish. Three-spined stickleback Gasterosteus aculeatus observed in the model school assay exhibited consistent individual differences in schooling behaviour.

Understanding how animal groups achieve coordinated movement

Moving animal groups display remarkable feats of coordination. This coordination is largely achieved when individuals adjust their movement in response to their neighbours' movements and positions. Recent advancements in automated tracking technologies, including computer vision and GPS, now allow researchers to gather large amounts of data on the movements and positions of individuals in groups. Furthermore, analytical techniques from fields such as statistical physics now allow us to identify the precise interaction rules used by animals on the move. These interaction rules differ not only between species, but also between individuals in the same group. These differences have wide-ranging implications, affecting how groups make collective decisions and driving the evolution of collective motion. Here, I describe how trajectory data can be used to infer how animals interact in moving groups. I give examples of the similarities and differences in the spatial and directional organisations of animal groups between species, and discuss the rules that animals use to achieve this organisation. I then explore how groups of the same species can exhibit different structures, and ask whether this results from individuals adapting their interaction rules. I then examine how the interaction rules between individuals in the same groups can also differ, and discuss how this can affect ecological and evolutionary processes. Finally, I suggest areas of future research.

Quantitative trait locus analysis of body shape divergence in nine-spined sticklebacks based on high-density SNP-panel

Heritable phenotypic differences between populations, caused by the selective effects of distinct environmental conditions, are of commonplace occurrence in nature. However, the actual genomic targets of this kind of selection are still poorly understood. We conducted a quantitative trait locus (QTL) mapping study to identify genomic regions responsible for morphometric differentiation between genetically and phenotypically divergent marine and freshwater nine-spined stickleback (Pungitius pungitius) populations. Using a dense panel of SNP-markers obtained by restriction site associated DNA sequencing of an F2 recombinant cross, we found 22 QTL that explained 3.5-12.9% of phenotypic variance in the traits under investigation. We detected one fairly large-effect (PVE = 9.6%) QTL for caudal peduncle length-a trait with a well-established adaptive function showing clear differentiation among marine and freshwater populations. We also identified two large-effect QTL for lateral plate numbers, which are different from the lateral plate QTL reported in earlier studies of this and related species. Hence, apart from identifying several large-effect QTL in shape traits showing adaptive differentiation in response to different environmental conditions, the results suggest intra- and interspecific heterogeneity in the genomic basis of lateral plate number variation.

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.

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.